Self-adjusting entropy-stable scheme for compressible Euler equations

In this work,a self-adjusting entropy-stable scheme is proposed for solving compressible Euler equations.The entropy-stable scheme is constructed by combining the entropy conservative flux with a suitable diffusion operator.The entropy has to be preserved in smooth solutions and be dissipated at shocks.To achieve this,a switch function,which is based on entropy variables,is employed to make the numerical diffusion term be automatically added around discontinuities.The resulting scheme is still entropy-stable.A number of numerical experiments illustrating the robustness and accuracy of the scheme are presented.From these numerical results,we observe a remarkable gain in accuracy.     

Trajectory attractors for dissipative 2D Euler and Navier-Stokes equations

A trajectory attractor is constructed for the 2D Euler system containing an additional dissipation term −ru, r > 0, with periodic boundary conditions. The corresponding dissipative 2D Navier-Stokes system with the same term −ru and with viscosity v > 0 also has a trajectory attractor, . Such systems model large-scale geophysical processes in atmosphere and ocean (see [1]). We prove that → as v → 0+ in the corresponding metric space. Moreover, we establish the existence of the minimal limit of the trajectory attractors as v → 0+. We prove that is a connected invariant subset of . The connectedness problem for the trajectory attractor by itself remains open. Dedicated to the memory of Leonid Volevich Partially supported by the Russian Foundation for Basic Research (projects no 08-01-00784 and 07-01-00500). The first author has been partially supported by a research grant from the Caprio Foundation, Landau Network-Cento Volta.     

Stabilized non-dissipative approximations of Euler equations in generalized curvilinear coordinates

We discuss stabilization strategies for finite-difference approximations of the compressible Euler equations in generalized curvilinear coordinates that do not rely on explicit upwinding or filtering of the physical variables. Our approach rather relies on a skew-symmetric-like splitting of the convective derivatives, that guarantees preservation of kinetic energy in the semi-discrete, low-Mach-number limit. A locally conservative formulation allows efficient implementation and easy incorporation into existing compressible flow solvers. The validity of the approach is tested for benchmark flow cases, including the propagation of a cylindrical vortex, and the head-on collision of two vortex dipoles. The tests support high accuracy and superior stability over conventional central discretization of the convective derivatives. The potential use for DNS/LES of turbulent compressible flows in complex geometries is discussed.     

A compact finite difference scheme for the fractional sub-diffusion equations

In this paper, a compact finite difference scheme for the fractional sub-diffusion equations is derived. After a transformation of the original problem, the L1 discretization is applied for the time-fractional part and fourth-order accuracy compact approximation for the second-order space derivative. The unique solvability of the difference solution is discussed. The stability and convergence of the finite difference scheme in maximum norm are proved using the energy method, where a new inner product is introduced for the theoretical analysis. The technique is quite novel and different from previous analytical methods. Finally, a numerical example is provided to show the effectiveness and accuracy of the method.     

An Asymptotic Preserving scheme for the Euler equations in a strong magnetic field

This paper is concerned with the numerical approximation of the isothermal Euler equations for charged particles subject to the Lorentz force (the ‘Euler–Lorentz’ system). When the magnetic field is large, or equivalently, when the parameter ε$\epsilon$ representing the non-dimensional ion cyclotron frequency tends to zero, the so-called drift-fluid (or gyro-fluid) approximation is obtained. In this limit, the parallel motion relative to the magnetic field direction splits from perpendicular motion and is given implicitly by the constraint of zero total force along the magnetic field lines. In this paper, we provide a well-posed elliptic equation for the parallel velocity which in turn allows us to construct an Asymptotic-Preserving (AP) scheme for the Euler–Lorentz system. This scheme gives rise to both a consistent approximation of the Euler–Lorentz model when ε$\epsilon$ is finite and a consistent approximation of the drift limit when ε→0$\epsilon \to 0$. Above all, it does not require any constraint on the space and time-steps related to the small value of ε$\epsilon$. Numerical results are presented, which confirm the AP character of the scheme and its Asymptotic Stability.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

One and two soliton solutions for seventh-order Caudrey-Dodd-Gibbon and Caudrey-Dodd-Gibbon-KP equations

In this work, we explore more applications of the simplified form of the bilinear method to the seventhorder Caudrey-Dodd-Gibbon (CDG) and the Caudrey-Dodd-Gibbon-KP (CDG-KP) equation. We formally derive one and two soliton solutions for each equation. We also show that the two equations do not show resonance.     

A third-order finite-volume residual-based scheme for the 2D Euler equations on unstructured grids

A residual-based (RB) scheme relies on the vanishing of residual at the steady-state to design a transient first-order dissipation, which becomes high-order at steady-state. Initially designed within a finite-difference framework for computations of compressible flows on structured grids, the RB schemes displayed good convergence, accuracy and shock-capturing properties which motivated their extension to unstructured grids using a finite volume (FV) method. A second-order formulation of the FV–RB scheme for compressible flows on general unstructured grids was presented in a previous paper. The present paper describes the derivation of a third-order FV–RB scheme and its application to hyperbolic model problems as well as subsonic, transonic and supersonic internal and external inviscid flows.