Appendix E: Changelog
This is a list of changes made to DustPy since version v1.0.0 including discussions of their influence on the simulations.
v1.0.6
Release date: 30th June 2024
Simframe compatibility
DustPy compatible to the latest Simframe v1.0.5.Simframe has to be updated to at least version 1.0.5.Due to the changes in Simframe, the snapshots stored in sim.t.snapshots now need to additionally contain the time, at which the initial conditions are written. By default this is for t=0. Furthermore, 0-dimensional Fields exist now, such as sim.grid.Nr or sim.grid.Nm.
Bugfix to initial fluxes and source terms
A bug was writing false gas and dust fluxes and source terms into the output files. The results of previous simulations were not affected by this issue.
Bugfix to implicit fluxes and source terms
When using implicit integration, the gas and dust fluxes and source terms are computed after the integration step from the new surface density. A bug in referencing to the old values of the surface densities caused fluxes and source terms to be computed incorrectly. The simulations are not affected, since these quantities are only computed for convenience and post-simulation analysis, but values written to the output files may have been incorrect.
Modification to initial dust distribution
If sim.ini.dust.allowDriftingParticle is False (default) particles above a critical Stokes number are removed from the initial dust size distribution. There was an error in the computation of the Stokes number, which has been fixed in this version. However, the critical Stokes number has been adjusted, such that the initial conditions are almost identical to previous versions.
Modification to total dust flux function
The default function to compute the total dust flux has been modified to be downstream compatible to other software packages using DustPy. Simulations that are not customizing this function are not affected by this change.
v1.0.5
Release date: 3rd December 2023
Using Meson as build system
Due to the deprecation of numpy.distutils, DustPy is now using Meson as build system.
Bugfix to velocity distribution
For zero relative collision velocities a not-a-number error could have been triggered.
Bugfix to plotting script
A bug in the plotting script has been fixed causing lines not to be updated. May not be fixed for all Matplotlib version. Please use the most recent one.
Preparation for the addition of multiple gas species
In order to add multiple gas species in future versions, the Jacobian of the gas surface density has been modified. All previous models that have not specifically customized the gas Jacobian should be compatible with this version.
v1.0.4
Release date: 17th June 2023
Bugfix to boundary conditions
A bug prevented the boundaries to be set correctly for the boundary conditions "val", "grad", and "pow". The boundary conditions "const_val", "const_grad", and "const_pow" were not affected by this bug.
v1.0.3
Release date: 7th February 2023
Correction to inital particle size distribution
If the parameter Simulation.ini.dust.allowDriftingParticles is set to False before Simulation.initialize(), DustPy is removing particles in the outer disk that are initially close to the drift barrier. An error has been corrected in the conversion of Stokes number to particles size, in which the particle bulk density Simulation.dust.rhos has been in the numerator instead of the denominator. The effect of this correction on simulations is negligible, since
Simulation.dust.rhos ~ 1.
Removal of non-ASCII characters
Due to potential incompabilities all non-ASCII character have been removed from the docstrings in the Fortran files.
v1.0.2
Release date: 18th November 2022
Change in default temperature profile
A factor of \(\frac{1}{2}\) was added to the luminosisty in the standard update function of the gas temperature Simulation.gas.T, when calculating the midplane equilibrium temperature in accordance with e.g. Chiang & Goldreich (1997) or Dullemond et al. (2001):
\(T\left( r \right) = \sqrt[4]{\frac{1}{2}\frac{0.05\,L_*}{4\,\pi\,r^2\,\sigma_\mathrm{SB}}}\)
Compared to previous versions the midplane temperature is therefore reduced by a factor of \(\sqrt[4]{\frac{1}{2}}\approx0.84\). For a comparison with radiative transfer calculation with RADMC-3D please have a look at this example in the dustpylib documentation.
v1.0.1
Change to Collision Kernel
A factor was added to the standard update function of Simulation.dust.kernel, that reduces the collision rates of equal size collisions by a factor of \(\frac{1}{2}\):
\(K_{ij} = \frac{1}{1+\delta_{ij}} \frac{\pi\left( a_i + a_j \right)^2}{\sqrt{2\pi\left( H_i^2 + H_j^2 \right)}} v_\mathrm{rel}\)
Collisions between particles of two separate reservoirs \(i\) and \(j\) occur at a rate \(\propto N_i N_j\), since every particle in one reservoir can collide with every particle in the other reservoir. Collisions of particles in the same reservoir occur at reduces rates: the first particle can collide with \(N_i-1\) other particles, the second particle with \(N_i-2\) particles, and so on, to avoid double counting of collisions. In that case the collision rates are \(\propto \frac{1}{2}N_i\left(N_i-1\right)\), which is \(\frac{1}{2}N_i^2\) in the limit of large \(N_i\).
This change has no noticable effect on simulations performed with previous versions, since it is only reducing the collision rate for a single type of collisions.
v1.0.0
Release date: 3rd July 2022