You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
This commit was created on GitHub.com and signed with GitHub’s verified signature.
The key has expired.
v2023.0.0 - 10th Mar 2023
Physics
Dust nucleation: chemical network allows for self-consistent (carbon-rich) dust formation in simulations once the temperature drops below the condensation temperature, e.g. in AGB star winds (Siess et al. 2022). Just set DUST_NUCLEATION=yes
Implementation of implicit flux limited diffusion solver for radiation based on Whitehouse & Bate 2004, Whitehouse, Bate & Monaghan (2005) and Bate & Keto (2015). The implementation is a direct port of the solver used in Matthew Bate's sphNG code, generously contributed by Matthew and ported to phantom by Mike Lau and Daniel Price.
Compute and store a mass accretion rate for each sink particle based on the accreted mass over one dtmax interval. This can be used in MCFOST to provide accretion luminosity from each sink particle (Borchert et al. 2022a, 2022b)
Allow for use of PdV work and shock heating as source terms in MCFOST when performing live-coupled simulations with phantom+MCFOST (thanks to Elli Borchert and Sahl Rowther, #344; see appendix of Borchert et al. 2022b)
New equation of state that implements gas plus radiation pressure plus recombination of Hydrogen and Helium (Lau et al. 2022; #267)
New stratified disc locally isothermal equation of state added (thanks to Caitlyn Hardiman; #268, #278, #308)
GR simulations now work with the gas+radiation equation of state by implementing an option to evolve "s" as the entropy instead of P/rho^gamma (thanks to Fitz Hu; #250, #324)
New spherical raytracer for sink particle radiation / dust acceleration, particularly in the context of stellar winds (thanks to Mats Esseldeurs, Lionel Siess, Ward Homan)
option to use reconstruction on velocities to reduce dissipation away from shocks in the GR code (where the standard Cullen & Dehnen shock switch cannot be used)
added conservation check which stops the MHD code if hdivB/B is too large AND hdivbbmax_max == 1 (thanks to James Wurster)
added an override to force sink creation if density is too high; also an addition step in the barotropic equation of state at very high densities (thanks to James Wurster)
option for heating from sink particles, which adds heats all particles within the softening radius of the sink according to the sink particle luminosity (#216)
mass-weighted interpolations of dust-gas quantities are now default in dust growth (#377; thanks to Stephane Michoulier)
option for Aeolian erosion in dust growth module (#365, #372; thanks to Stephane Michoulier)
Added possibility of varying beta cooling with radius as a power law (#361; thanks to Cristiano Longarini)
option for damping boundary conditions for comparing phantom to simulations of planet-disc interaction with grid codes (#351)
Setup
can now set up hierarchical quadruple star systems in setup_disc (thanks to Amena Faruqi; #355)
setup and simple cooling function added for cooling shock problem from Creasey et al. (2011)
can now give --maxp=1e7 flag to phantomsetup instead of compiling with MAXP= which is now obsolete
major reorganisation of setup_star into neat subroutines with high level functions (#297, also #203)
dumps written during relax_star now usable as starting file
various improvements to setting up stars with variable composition (thanks to Mike Lau)