THE RIEMANN PROBLEM FOR ISENTROPIC COMPRESSIBLE EULER EQUATIONS WITH DISCONTINUOUS FLUX

We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.

Radon Measure Solutions to Riemann Problems for Isentropic Compressible Euler Equations of Polytropic Gases

We solve the Riemann problems for isentropic compressible Euler equations of polytropic gases in the class of Radon measures,and the solutions admit the concentration of mass.It is found that under the requirement of satisfying the over-compressing entropy condition:(i)there is a unique delta shock solution,corresponding to the case that has two strong classical Lax shocks;(ii)for the initial data that the classical Riemann solution contains a shock wave and a rarefaction wave,or two shocks with one being weak,there are infinitely many solutions,each consists of a delta shock and a rarefaction wave;(iii)there are no delta shocks for the case that the classical entropy weak solutions consist only of rarefaction waves.These solutions are self-similar.Furthermore,for the generalized Riemann problem with mass concentrated initially at the discontinuous point of initial data,there always exists a unique delta shock for at least a short time.It could be prolonged to a global solution.Not all the solutions are self-similar due to the initial velocity of the concentrated point-mass(particle).Whether the delta shock solutions constructed satisfy the over-compressing entropy condition is clarified.This is the first result on the construction of singular measure solutions to the compressible Euler system of polytropic gases,that is strictly hyperbolic,and whose characteristics are both genuinely nonlinear.We also discuss possible physical interpretations and applications of these new solutions.

Blowup of Solutions to the Non-Isentropic Compressible Euler Equations with Time-Dependent Damping and Vacuum

This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data is assumed to be radially symmetric and the initial density contains vacuum, we obtain that classical solution, especially the density, will blow up on finite time. The results also reveal that damping can really delay the singularity formation.

Steady Compressible Euler Equations of Concentration Layers for Hypersonic-limit Flows Passing Three-dimensional Bodies and Generalized Newton-Busemann Pressure Law

For stationary hypersonic-limit Euler flows passing a solid body in three-dimensional space,the shock-front coincides with the upwind surface of the body,hence there is an infinite-thin layer of concentrated mass,in which all particles hitting the body move along its upwind surface.By proposing a concept of Radon measure solutions of boundary value problems of the multi-dimensional compressible Euler equations,which incorporates the large-scale of three-dimensional distributions of upcoming hypersonic flows and the small-scale of particles moving on two-dimensional surfaces,the authors derive the compressible Euler equations for flows in concentration layers,which is a stationary pressureless compressible Euler system with source terms and independent variables on curved surface.As a by-product,they obtain a formula for pressure distribution on surfaces of general obstacles in hypersonic flows,which is a generalization of the classical Newton-Busemann law for drag/lift in hypersonic aerodynamics.

A Third Order Conservative Lagrangian Type Scheme on Curvilinear Meshes for the Compressible Euler Equations

Based on the high order essentially non-oscillatory(ENO)Lagrangian type scheme on quadrilateral meshes presented in our earlier work[3],in this paper we develop a third order conservative Lagrangian type scheme on curvilinear meshes for solving the Euler equations of compressible gas dynamics.The main purpose of this work is to demonstrate our claim in[3]that the accuracy degeneracy phenomenon observed for the high order Lagrangian type scheme is due to the error from the quadrilateral mesh with straight-line edges,which restricts the accuracy of the resulting scheme to at most second order.The accuracy test given in this paper shows that the third order Lagrangian type scheme can actually obtain uniformly third order accuracy even on distorted meshes by using curvilinear meshes.Numerical examples are also presented to verify the performance of the third order scheme on curvilinear meshes in terms of resolution for discontinuities and non-oscillatory properties.

A Runge Kutta Discontinuous Galerkin Method for Lagrangian Compressible Euler Equations in Two-Dimensions

This paper presents a new Lagrangian type scheme for solving the Euler equations of compressible gas dynamics.In this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method,and the mesh moves with the fluid flow.The scheme is conservative for the mass,momentum and total energy and maintains second-order accuracy.The scheme avoids solving the geometrical part and has free parameters.Results of some numerical tests are presented to demonstrate the accuracy and the non-oscillatory property of the scheme.

Analytical Blowup Solutions to the Compressible Euler Equations with Time-depending Damping

In this paper,the analytical blowup solutions of the N-dimensional radial symmetric compressible Euler equations are constructed.Some previous results of the blowup solutions for the compressible Euler equations with constant damping are generalized to the time-depending damping case.The generalization is untrivial because that the damp coefficient is a nonlinear function of time t.

L^(2)-CONVERGENCE TO NONLINEAR DIFFUSION WAVES FOR EULER EQUATIONS WITH TIME-DEPENDENT DAMPING

In this paper,we are concerned with the asymptotic behavior of L^(∞) weak-entropy solutions to the compressible Euler equations with a vacuum and time-dependent damping-m/(1+t)^(λ).As λ∈(0,l/7],we prove tht the L^(∞) weak-entropy solution converges to the nonlinear diffusion wave of the generalized porous media equation(GPME)in L^(2)(R).As λ∈(1/7,1),we prove that the L^(∞) weak-entropy solution converges to an expansion around the nonlinear diffusion wave in L^(2)(R),which is the best asymptotic profile.The proof is based on intensive entropy analysis and an energy method.

EXISTENCE AND UNIQUENESS OF THE GLOBAL L^(1) SOLUTION OF THE EULER EQUATIONS FOR CHAPLYGIN GAS

In this paper, we establish the global existence and uniqueness of the solution of the Cauchy problem of a one-dimensional compressible isentropic Euler system for a Chaplygin gas with large initial data in the space Lloc1. The hypotheses on the initial data may be the least requirement to ensure the existence of a weak solution in the Lebesgue measurable sense. The novelty and also the essence of the difficulty of the problem lie in the fact that we have neither the requirement on the local boundedness of the density nor that which is bounded away from vacuum. We develop the previous results on this degenerate system.The method used is Lagrangian representation, the essence of which is characteristic analysis.The key point is to prove the existence of the Lagrangian representation and the absolute continuity of the potentials constructed with respect to the space and the time variables.We achieve this by finding a property of the fundamental theorem of calculus for Lebesgue integration, which is a sufficient and necessary condition for judging whether a monotone continuous function is absolutely continuous. The assumptions on the initial data in this paper are believed to also be necessary for ruling out the formation of Dirac singularity of density. The ideas and techniques developed here may be useful for other nonlinear problems involving similar difficulties.

Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics

We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.

Preventing Pressure Oscillations Does Not Fix Local Linear Stability Issues of Entropy-Based Split-Form High-Order Schemes

Recently,it was discovered that the entropy-conserving/dissipative high-order split-form discontinuous Galerkin discretizations have robustness issues when trying to solve the sim-ple density wave propagation example for the compressible Euler equations.The issue is related to missing local linear stability,i.e.,the stability of the discretization towards per-turbations added to a stable base flow.This is strongly related to an anti-diffusion mech-anism,that is inherent in entropy-conserving two-point fluxes,which are a key ingredi-ent for the high-order discontinuous Galerkin extension.In this paper,we investigate if pressure equilibrium preservation is a remedy to these recently found local linear stability issues of entropy-conservative/dissipative high-order split-form discontinuous Galerkin methods for the compressible Euler equations.Pressure equilibrium preservation describes the property of a discretization to keep pressure and velocity constant for pure density wave propagation.We present the full theoretical derivation,analysis,and show corresponding numerical results to underline our findings.In addition,we characterize numerical fluxes for the Euler equations that are entropy-conservative,kinetic-energy-preserving,pressure-equilibrium-preserving,and have a density flux that does not depend on the pressure.The source code to reproduce all numerical experiments presented in this article is available online(https://doi.org/10.5281/zenodo.4054366).

A Reconstructed Discontinuous Galerkin Method for the Euler Equations on Arbitrary Grids

A reconstruction-based discontinuous Galerkin(RDG(P1P2))method,a variant of P1P2 method,is presented for the solution of the compressible Euler equations on arbitrary grids.In this method,an in-cell reconstruction,designed to enhance the accuracy of the discontinuous Galerkin method,is used to obtain a quadratic polynomial solution(P2)from the underlying linear polynomial(P1)discontinuous Galerkin solution using a least-squares method.The stencils used in the reconstruction involve only the von Neumann neighborhood(face-neighboring cells)and are compact and consistent with the underlying DG method.The developed RDG method is used to compute a variety of flow problems on arbitrary meshes to demonstrate its accuracy,efficiency,robustness,and versatility.The numerical results indicate that this RDG(P1P2)method is third-order accurate,and outperforms the third-order DG method(DG(P2))in terms of both computing costs and storage requirements.

Global Stability to Steady Supersonic Rayleigh Flows in One-Dimensional Duct

Heat exchange plays an important role in hydrodynamical systems,which is an interesting topic in theory and application.In this paper,the authors consider the global stability of steady supersonic Rayleigh flows for the one-dimensional compressible Euler equations with heat exchange,under the small perturbations of initial and boundary conditions in a finite rectilinear duct.

Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Quasi-neutral limit of the full bipolar Euler-Poisson system

The quasi-neutral limit of the multi-dimensional non-isentropic bipolar Euler-Poisson system is considered in the present paper. It is shown that for well-prepared initial data the smooth solution of the non-isentropic bipolar Euler-Poisson system converges strongly to the compressible non-isentropic Euler equations as the Debye length goes to zero.

Continuous Finite Element Subgrid Basis Functions for Discontinuous Galerkin Schemes on Unstructured Polygonal Voronoi Meshes

We propose a new high order accurate nodal discontinuous Galerkin(DG)method for the solution of nonlinear hyperbolic systems of partial differential equations(PDE)on unstructured polygonal Voronoi meshes.Rather than using classical polynomials of degree N inside each element,in our new approach the discrete solution is represented by piecewise continuous polynomials of degree N within each Voronoi element,using a continuous finite element basis defined on a subgrid inside each polygon.We call the resulting subgrid basis an agglomerated finite element(AFE)basis for the DG method on general polygons,since it is obtained by the agglomeration of the finite element basis functions associated with the subgrid triangles.The basis functions on each sub-triangle are defined,as usual,on a universal reference element,hence allowing to compute universal mass,flux and stiffness matrices for the subgrid triangles once and for all in a pre-processing stage for the reference element only.Consequently,the construction of an efficient quadrature-free algorithm is possible,despite the unstructured nature of the computational grid.High order of accuracy in time is achieved thanks to the ADER approach,making use of an element-local space-time Galerkin finite element predictor.The novel schemes are carefully validated against a set of typical benchmark problems for the compressible Euler and Navier-Stokes equations.The numerical results have been checked with reference solutions available in literature and also systematically compared,in terms of computational efficiency and accuracy,with those obtained by the corresponding modal DG version of the scheme.