Computation of fluid flows in non-inertial contracting, expanding, and rotating reference frames

We present the method for computation of fluid flows that are characterized by the large degree of expansion/contraction and in which the fluid velocity is dominated by the bulk component associated with the expansion/contraction and/or rotation of the flow. We consider the formulation of Euler equations of fluid dynamics in a homologously expanding/contracting and/or rotating reference frame. The frame motion is adjusted to minimize local fluid velocities. Such approach allows to accommodate very efficiently large degrees of change in the flow extent. Moreover, by excluding the contribution of the bulk flow to the total energy the method eliminates the high Mach number problem in the flows of interest. An important practical advantage of the method is that it can be easily implemented with virtually any Eulerian hydrodynamic scheme and adaptive mesh refinement (AMR) strategy. We also consider in detail equation invariance and existence of conservative formulation of equations for special classes of expanding/contracting reference frames. Special emphasis is placed on extensive numerical testing of the method for a variety of reference frame motions, which are representative of the realistic applications of the method. We study accuracy, conservativity, and convergence properties of the method both in problems which are not its optimal applications as well as in systems in which the use of this method is maximally beneficial. Such detailed investigation of the numerical solution behavior is used to define the requirements that need to be considered in devising problem-specific fluid motion feedback mechanisms.

[1]  Symmetries of fluid dynamics with polytropic exponent , 2000, hep-th/0009092.

[2]  F. S. Sherman,et al.  Unsteady motion of continuous media , 1960 .

[3]  Michael L. Norman,et al.  The Formation of the First Star in the Universe , 2001, Science.

[4]  R. LeVeque,et al.  Adaptive Mesh Refinement Using Wave-Propagation Algorithms for Hyperbolic Systems , 1998 .

[5]  Alexei M. Khokhlov,et al.  Fully Threaded Tree Algorithms for Adaptive Refinement Fluid Dynamics Simulations , 1997, astro-ph/9701194.

[6]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

[7]  H. Trac,et al.  A moving frame algorithm for high Mach number hydrodynamics , 2003, astro-ph/0309599.

[8]  W. Hillebrandt,et al.  Full-star type Ia supernova explosion models , 2005 .

[9]  A. I. Saichev,et al.  The large-scale structure of the Universe in the frame of the model equation of non-linear diffusion , 1989 .

[10]  R. LeVeque.  Wave Propagation Algorithms for Multidimensional Hyperbolic Systems , 1997 .

[12]  A. Frank,et al.  AstroBEAR: AMR for Astrophysical Applications - I: Methods , 2005 .

[13]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[14]  R. Salmon,et al.  Geophysical Fluid Dynamics , 2019, Classical Mechanics in Geophysical Fluid Dynamics.

[15]  William H. Press,et al.  Numerical recipes in Fortran 77 : the art of scientificcomputing. , 1992 .

[16]  Riccardo Fazio.  Numerical Applications of the Scaling Concept , 1999 .

[17]  M. L. Norman,et al.  Simulating Radiating and Magnetized Flows in Multiple Dimensions with ZEUS-MP , 2005, astro-ph/0511545.

[18]  Michael L. Norman,et al.  Adaptive-mesh radiation hydrodynamics. I - The radiation transport equation in a completely adaptive coordinate system. II - The radiation and fluid equations in relativistic flows , 1984 .

[19]  Following multi-dimensional type Ia supernova explosion models to homologous expansion , 2004, astro-ph/0408296.

[20]  The Maximal Kinematical Invariance Group of Fluid Dynamics and Explosion–Implosion Duality , 2000, hep-th/0007199.

[21]  O. Jahn,et al.  Symmetries of discontinuous flows and the dual Rankine–Hugoniot conditions in fluid dynamics , 2004, math-ph/0407061.

[22]  R. Anderson,et al.  An arbitrary Lagrangian-Eulerian method with adaptive mesh refinement for the solution of the Euler equations , 2004 .

[23]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[24]  R.,et al.  Moving-Mesh Methods for One-Dimensional Hyperbolic Problems Using CLAWPACK , 2003 .

[25]  P. Rentrop,et al.  Generalized Runge-Kutta methods of order four with stepsize control for stiff ordinary differential equations , 1979 .

[26]  J. T. Mendonça,et al.  Explosion implosion duality and the laboratory simulation of astrophysical systems , 2000, astro-ph/0003385.

[27]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[28]  The maximal invariance group of Newton’s equations for a free point particle , 2001, math-ph/0102011.

[29]  S. Penner.  Physics of shock waves and high-temperature hydrodynamic phenomena - Ya.B. Zeldovich and Yu.P. Raizer (translated from the Russian and then edited by Wallace D. Hayes and Ronald F. Probstein); Dover Publications, New York, 2002, 944 pp., $34. , 2003 .