Capturing Turbulence with Numerical Dissipation: a Simple Dynamical Model for Unresolved Turbulence in Hydrodynamic Simulations

Vadim A. Semenov
2024-10-31
Abstract:Modeling unresolved turbulence in astrophysical gasdynamic simulations can improve the modeling of other subgrid processes dependent on the turbulent structure of gas: from flame propagation in the interiors of combusting white dwarfs to star formation and chemical reaction rates in the interstellar medium, and non-thermal pressure support of circum- and intergalactic gas. We present a simple method for modeling unresolved turbulence in hydrodynamic simulations via tracking its sourcing by local numerical dissipation and modeling its decay into heat. This method is physically justified by the generic property of turbulent flows that they dissipate kinetic energy at a rate set by the energy cascade rate from large scales, which is independent of fluid viscosity regardless of its nature, be it physical or numerical. We calibrate and test our model against decaying supersonic turbulence simulations. Despite its simplicity, the model quantitatively reproduces multiple non-trivial features of the high-resolution turbulence run: the temporal evolution of the average small-scale turbulence, its dependence on spatial scale, and the slope and scatter of the local correlation between subgrid turbulent velocities, gas densities, and local compression rates. As an example of practical applications, we use our model in isolated galactic disk simulations to model locally variable star formation efficiency at the subresolution scale. In the supersonic, star-forming gas, the new model performs comparably to a more sophisticated model where the turbulent cascade is described by explicit subgrid terms. Our new model is straightforward to implement in many hydrodynamic codes used in galaxy simulations as it utilizes already existing infrastructure to implicitly track the numerical dissipation in such codes.
Astrophysics of Galaxies,Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?