stammler / dustpy

Python Package for Modelling Dust Evolution in Protoplanetary Disks.

DustPy

Dust Coagulation and Evolution in Protoplanetary Disks

DustPy is a Python package to simulate the evolution of dust in protoplanetary disks.

DustPy simulates the radial evolution of gas and dust in a protoplanetary disk, involving viscous evolution the the gas disk, advection and diffusion of the dust disk, as well as dust growth by solving the Smoluchowski equation.

Please read the documentation for a detailed description. By using any version of DustPy you agree to these terms of usage.

Installation

DustPy can be installed via the Python Package Index

pip install dustpy

Requirements

DustPy needs a Python3 distribution and a Fortran compiler installed on your system.

Documentation

https://stammler.github.io/dustpy/

• 1. Basics
• 2. Simple Customization
• 3. Advanced Customization
• 4. The Standard Model
• 5. Dust Coagulation
• 6. Dust Evolution
• 7. Gas Evolution
   
• Library: dustpylib
   
• Test: Analytical Coagulation Kernels
• Test: Gas Evolution
   
• Example: Ice Lines
• Example: Planetary Gaps
• Example: Planetesimal Formation
   
   
• A. Citation
• B. List of Publications
• C. Contributing/Bugs/Features
• D. DustPy Discussions
• E. Changelog

Module Reference

Framework

DustPy is using the Simframe framework for scientific simulations (Stammler & Birnstiel 2022)

Acknowledgements

DustPy has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 714769.

DustPy was developed at the University Observatory of the Ludwig Maximilian University of Munich.

dustpy's Issues

Tracking changes in isotopic content of grid

I notice that Stammler has already studied "How fragmentation impacts dust filtering by planets", which is related to NN and CN dichotomy. However, I find it knotty to directly track how isotopic changes in the grid along with the evolution of mass bin.

Since the calculation of the changes in the mass bin has already been encapsulated and changes of isotopic correspond to mass bin, I'm not sure if simply adding a subgroup to Grid can achieve this purpose.

I'm not sure if there's a better alternative yet.

Another approach to achieve planetary gap

Example in documentation has demonstrated how to construct planetary gap by changing the alpha viscosity.

Here's another way to achieve planetary gap, which is directly adding a torch in evolution equation(Lin & Papaloizou(1986)):
$\frac{\partial \Sigma}{\partial t} = \frac{1}{R}\frac{\partial}{\partial R}[3R^{1/2} \frac{\partial }{\partial R}(\nu \Sigma R^{1/2}) - \frac{2 A\Sigma R^{3/2}}{(G M_\star )^{1/2}}]$
where A is torque density distribution.

The core issue is that Dustpy does not seem to allow changing the form of the evolution equation (torque cannot be treated as an external source as well).

Is it possible to implement a torque driven planetary gap without changing the fundamental code logic?

Installation not working on Python 3.12

The installation is currently not working on Python 3.12 or higher due to the removal of numpy.distutils, which was used to compile the Fortran files.