The Cauchy-Lagrangian method for numerical analysis of Euler flow

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.

Amplitude equations for weakly nonlinear two-scale perturbations of free hydromagnetic convective regimes in a rotating layer

Weakly non-linear stability of regimes of free hydromagnetic thermal convection in a rotating horizontal layer with free electrically conducting boundaries is considered in the Boussinesq approximation. Perturbations are supposed to involve large spatial and temporal scales. Applying methods for homogenisation of parabolic equations, we derive the system of amplitude equations governing the evolution of perturbations under the assumption that the alpha-effect is insignificant in the leading order. The amplitude equations involve the operators of anisotropic combined eddy diffusivity correction and advection. The system is qualitatively different from the system of mean-field equations for large-scale perturbations of forced convective hydromagnetic regimes. It is mixed: equations for the mean magnetic perturbation are evolutionary, all the rest involve neither time derivatives, nor the molecular diffusivity operator.