The Cauchy-Lagrangian method for numerical analysis of Euler flow

Abstract

A novel semi-Lagrangian method is introduced to solve numerically the Euler equation for ideal incompressible flow in arbitrary space dimension. It exploits the time-analyticity of fluid particle trajectories and requires, in principle, only limited spatial smoothness of the initial data. Efficient generation of high-order time-Taylor coefficients is made possible by a recurrence relation that follows from the Cauchy invariants formulation of the Euler equation (Zheligovsky & Frisch, J. Fluid Mech. 2014, 749, 404-430). Truncated time-Taylor series of very high order allow the use of time steps vastly exceeding the Courant-Friedrichs-Lewy limit, without compromising the accuracy of the solution. Tests performed on the two-dimensional Euler equation indicate that the Cauchy-Lagrangian method is more - and occasionally much more - efficient and less prone to instability than Eulerian Runge-Kutta methods, and less prone to rapid growth of rounding errors than the high-order Eulerian time-Taylor algorithm. We also develop tools of analysis adapted to the Cauchy-Lagrangian method, such as the monitoring of the radius of convergence of the time-Taylor series. Certain other fluid equations can be handled similarly.

30 pp., 13 figures, 45 references. Minor revision. In press in Journal of Scientific Computing

References in corpus (3)

• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• Cauchy's almost forgotten Lagrangian formulation of the Euler equation for 3D incompressible flow
• How smooth are particle trajectories in a $\Lambda$CDM Universe?

Cited by in corpus (17)

• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Geometric formulation of the Cauchy invariants for incompressible Euler flow in flat and curved spaces
• Shell-crossing in quasi-one-dimensional flow
• Shell-crossing in a $\Lambda$CDM Universe
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Origin of Enhanced Skin Friction at the Onset of Boundary-Layer Transition
• A Characteristic Mapping Method for the two-dimensional incompressible Euler equations
• Lagrangian theory for cosmic structure formation with vorticity: Newtonian and post-Friedmann approximations
• A Characteristic Mapping Method for the three-dimensional incompressible Euler equations
• Insights from a pseudospectral study of a potentially singular solution of the three-dimensional axisymmetric incompressible Euler equation
• Quasi-spherical collapse of matter in $\Lambda$CDM
• Eye of the Tyger: early-time resonances and singularities in the inviscid Burgers equation
• A constructive approach to regularity of Lagrangian trajectories for incompressible Euler flow in a bounded domain
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• Computation of kinematic and magnetic $\alpha$-effect and eddy diffusivity tensors by Padé approximation
• Fast and accurate collapse-time predictions for collisionless matter
• Singularity formation of vortex sheets in 2D Euler equations using the characteristic mapping method