A Primer on Eulerian Computational Fluid Dynamics for Astrophysics

Abstract

We present a pedagogical review of some of the methods employed in Eulerian computational fluid dynamics (CFD). Fluid mechanics is governed by the Euler equations, which are conservation laws for mass, momentum, and energy. The standard approach to Eulerian CFD is to divide space into finite volumes or cells and store the cell-averaged values of conserved hydro quantities. The integral Euler equations are then solved by computing the flux of the mass, momentum, and energy across cell boundaries. We review both first-order and second-order flux assignment schemes. All linear schemes are either dispersive or diffusive. The nonlinear, second-order accurate total variation diminishing (TVD) approach provides high resolution capturing of shocks and prevents unphysical oscillations. We review the relaxing TVD scheme, a simple and robust method to solve systems of conservation laws like the Euler equations. A 3-D relaxing TVD code is applied to the Sedov-Taylor blast wave test. The propagation of the blast wave is accurately captured and the shock front is sharply resolved. We apply a 3-D self-gravitating hydro code to simulating the formation of blue straggler stars through stellar mergers and present some numerical results. A sample 3-D relaxing TVD code is provided in the appendix.

23 pages, 12 figures; includes sample 3-D hydro code and new section on stellar mergers; accepted by PASP

Cited by in corpus (29)

• Nonthermal Radiation of Young Supernova Remnants
• Axisymmetric simulations of magneto--rotational core collapse: dynamics and gravitational wave signal
• GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics
• GAMER-2: a GPU-accelerated adaptive mesh refinement code -- accuracy, performance, and scalability
• On Aspects Pertaining to the Perpendicular Diffusion of Solar Energetic Particles
• Precision era of the kinetic Sunyaev-Zeldovich effect: simulations, analytical models and observations and the power to constrain reionization
• Numerical simulations of diffusive shock acceleration in SNRs
• A Free, Fast, Simple and Efficient TVD MHD Code
• Shocks in the Early Universe
• The intergalactic medium temperature and Compton y parameter
• Cosmic ray acceleration in magnetic circumstellar bubbles
• Numerical simulations of buoyancy instabilities in galaxy cluster plasmas with cosmic rays and anisotropic thermal conduction
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (I)
• Out-of-Core Hydrodynamic Simulations for Cosmological Applications
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (II)
• A comparison of hydrodynamics techniques for modelling collisions between main sequence stars
• The angular momentum structure of cosmic ray driven galactic outflows triggered by stream accretion
• Magnetohydrodynamics with GAMER
• The phase structure of cosmic ray driven outflows in stream fed disc galaxies
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Numerical simulations of composite supernova remnants for small $\sigma$ pulsar wind nebulae
• Hydrodynamical Simulations of Colliding Jets
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (III)
• RAPID: A fast, high resolution, flux-conservative algorithm designed for planet-disk interactions
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (IV)
• Differentiable Cosmological Hydrodynamics for Field-Level Inference and High Dimensional Parameter Constraints
• GAMER with out-of-core computation
• Boltzmann-Poisson-like approach to simulating the galactic halo response to satellite accretion Dependence on the halo density profile
• Kiloparsec-scale turbulence driven by reionization may grow intergalactic magnetic fields