Lagrangian methods are widely used in many fields for multi-material compressible flow simulations such as in astrophysics and inertial confinement fusion (ICF), due to their distinguished advantage in capturing material interfaces automatically. In some of these applications, multiple internal energy equations such as those for electron, ion and radiation are involved. In the past decades, several staggered-grid based Lagrangian schemes have been developed which are designed to solve the internal energy equation directly. These schemes can be easily extended to solve problems with multiple internal energy equations. However such schemes are typically not conservative for the total energy. Recently, significant progress has been made in developing cell-centered Lagrangian schemes which have several good properties such as conservation for all the conserved variables and easiness for remapping. However, these schemes are commonly designed to solve the Euler equations in the form of the total energy, therefore they cannot be directly applied to the solution of either the single internal energy equation or the multiple internal energy equations without significant modifications. Such modifications, if not designed carefully, may lead to the loss of some of the nice properties of the original schemes such as conservation of the total energy. In this paper, we establish an equivalency relationship between the cell-centered discretizations of the Euler equations in the forms of the total energy and of the internal energy. By a carefully designed modification in the implementation, the cell-centered Lagrangian scheme can be used to solve the compressible fluid flow with one or multiple internal energy equations and meanwhile it does not lose its total energy conservation property. An advantage of this approach is that it can be easily applied to many existing large application codes which are based on the framework of solving multiple internal energy equations. Several two dimensional numerical examples for both Euler equations and three-temperature hydrodynamic equations in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry preserving, accuracy and non-oscillatory performance.

A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.

In this paper, we propose and analyze a nodal-continuous and $H^1$-conforming finite element method for the numerical computation of Maxwell's equations, with singular solution in a fractional order Sobolev space $H^r(\Omega)$, where $r$ may take any value in the most interesting interval $(0, 1)$. The key feature of the method is that mass-lumping linear finite element $L^2$ projections act on the curl and divergence partial differential operators so that the singular solution can be sought in a setting of $L^2(\Omega)$ space. We shall use the nodal-continuous linear finite elements, enriched with one element bubble in each element, to approximate the singular and non-$H^1$ solution. Discontinuous and nonhomogeneous media are allowed in the method. Some error estimates are given and a number of numerical experiments for source problems as well as eigenvalue problems are presented to illustrate the superior performance of the proposed method.

The logarithmic nonlinearity has been used in many partial differential equations (PDEs) for modeling problems in various applications. Due to the singularity of the logarithmic function, it introduces tremendous difficulties in establishing mathematical theories, as well as in designing and analyzing numerical methods for PDEs with such nonlinearity. Here, we take the logarithmic Schrödinger equation (LogSE) as a prototype model. Instead of regularizing f(ρ)=lnρ in the LogSE directly and globally as being done in the literature, we propose a local energy regularization (LER) for the LogSE by first regularizing F(ρ)=ρlnρ−ρ locally near ρ=0+ with a polynomial approximation in the energy functional of the LogSE and then obtaining an energy regularized logarithmic Schrödinger equation (ERLogSE) via energy variation. Linear convergence is established between the solutions of ERLogSE and LogSE in terms of a small regularization parameter 0<ε≪1. Moreover, the conserved energy of the ERLogSE converges to that of LogSE quadratically, which significantly improves the linear convergence rate of the regularization method in the literature. Error estimates are also presented for solving the ERLogSE by using Lie–Trotter splitting integrators. Numerical results are reported to confirm our error estimates of the LER and of the time-splitting integrators for the ERLogSE. Finally, our results suggest that the LER performs better than regularizing the logarithmic nonlinearity in the LogSE directly.

No abstract uploaded!