Abstract:The radial one-dimensional Euler equations are often rewritten in what is known as the geometric source form. The differential operator is identical to the Cartesian case, but source terms result. Since the theory and numerical methods for the Cartesian case are well-developed, they are often applied without modification to cylindrical and spherical geometries. However, numerical conservation is lost. In this article, AUSM+\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$^+$$\end{document}-up is applied to a numerically conservative (discrete) form of the Euler equations labeled the geometric form, a nearly conservative variation termed the geometric flux form, and the geometric source form. The resulting numerical methods are compared analytically and numerically through three types of test problems: subsonic, smooth, steady-state solutions, Sedov’s similarity solution for point or line-source explosions, and shock tube problems. Numerical conservation is analyzed for all three forms in both spherical and cylindrical coordinates. All three forms result in constant enthalpy for steady flows. The spatial truncation errors have essentially the same order of convergence, but the rate constants are superior for the geometric and geometric flux forms for the steady-state solutions. Only the geometric form produces the correct shock location for Sedov’s solution, and a direct connection between the errors in the shock locations and energy conservation is found. The shock tube problems are evaluated with respect to feature location using an approximation with a very fine discretization as the benchmark. Extensions to second order appropriate for cylindrical and spherical coordinates are also presented and analyzed numerically. Conclusions are drawn, and recommendations are made. A derivation of the steady-state solution is given in the Appendix.

A Proper Discretization of Hydrodynamic Equations in Cylindrical Coordinates for Astrophysical Simulations

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

A New Coupled Computational Method in Conjunction with Three-Dimensional Finite Volume Schemes for Nonlinear Coupled Constitutive Relations

A Computational Strategy for Nonlinear Coupled Constitutive Relations of Rarefied Nonequilibrium Flows on Unstructured Grids

Cartesian grid method with a consistent direct discretization approach for numerical simulation of non-Newtonian fluid flow

Dynamics of Cylindrical Particles in the Contraction Flow of a Second-Order Fluid

A Primer on Eulerian Computational Fluid Dynamics for Astrophysics

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

Direct Numerical Simulation of Turbulent Pipe Flow.

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Multiscale Radiative Transfer in Cylindrical Coordinates

Analytical solutions for mantle flow in cylindrical and spherical shells

Discretization form of the continuity condition at the polar axis, with application to the gyrokinetic simulation in a magnetic fusion torus

A Compact Scheme for the Streamfunction-Velocity Formulation of the 2D Steady Incompressible Navier-Stokes Equations in Polar Coordinaes

Numerical Simulation of Viscous Flow Around Unrestrained Cylinders

The impact of the form of the Euler equations for radial flow in cylindrical and spherical coordinates on numerical conservation and accuracy

A new transient higher order compact scheme for computation of flow and heat transfer in nonuniform polar grids

The Athena Astrophysical MHD Code in Cylindrical Geometry

Cartesian Grid Method for Gas Kinetic Scheme

A high‐order compact scheme for solving the 2D steady incompressible Navier‐Stokes equations in general curvilinear coordinates