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.

ON THE GLOBAL EXISTENCE OF SMOOTH SOLUTIONS TO THE MULTI-DIMENSIONAL COMPRESSIBLE EULER EQUATIONS WITH TIME-DEPENDING DAMPING IN HALF SPACE

This paper is a continue work of [4, 5]. In the previous two papers, we studied the Cauchy problem of the multi-dimensional compressible Euler equations with time-depending damping term --u/（1＋t）λpu, where λ≥ 0 and μ 〉 0 are constants. We have showed that, for all λ ≥ 0 andμ 〉 0 the smooth solution to the Cauchy problem exists globally or blows up in finite time. In the present paper, instead of the Cauchy problem we consider the initial- boundary value problem in the half space R＋^d with space dimension d = 2, 3. With the help of the special structure of the equations and the fluid vorticity, we overcome the difficulty arisen from the boundary effect. We prove that there exists a global smooth solution for 0 ≤λ 〈 1 when the initial data is close to its equilibrium state. In addition, exponential decay of the fluid vorticity will also be established.

The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate

In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.

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.

Variational iteration method for solving compressible Euler equations

This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.

TWO-DIMENSIONAL RIEMANN PROBLEMS:FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS

In this paper we survey the authors＇ and related work on two-dimensional Riemann problems for hyperbolic conservation laws, mainly those related to the compressible Euler equations in gas dynamics. It contains four sections： 1. Historical review. 2. Scalar conservation laws. 3. Euler equations. 4. Simplified models.

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 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.

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.

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.

The Lifespan for 3-D Spherically Symmetric Compressible Euler Equations

In this paper, we study the lifespan of solutions for three dimensional compressible Euler equations with spherically symmetric initial data that is a small perturbation of amplitude ε from a constant state. From our result, the classical solutions have to blow up in finite time in spite of any small ε.

GLOBAL ENTROPY SOLUTIONS TO AN INHOMOGENEOUS ISENTROPIC COMPRESSIBLE EULER SYSTEM

In this article, we develop a new technique to prove the global existene of entropy solutions to an inhomogeneous isentropic compressible Euler equations through the compensated compactness and vanishing viscosity method. In particular, the entropy solutions are uniformly bounded independent of time.

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.

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.

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.

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.

ENTROPIES AND FLUX-SPLITTINGS FOR THE ISENTROPIC EULER EQUATIONS

The authors establish the existence of a large class of mathematical entropies (the so-called weak entropies) associated with the Euler equations for an isentropic, compressible fluid governed by a general pressure law. A mild assumption on the behavior of the pressure law near the vacuum is solely required. The analysis is based on an asymptotic expansion of the fundamental solution (called here the entropy kernel) of a highly singular Euler- Poisson-Darboux equation. The entropy kernel is only Holder continuous and its regularity is carefully inversti- gated. Relying on a notion introduced earlier by the authors, it is also proven that, for the Euler equations, the set of eatropy flux-splittings coincides with the set of entropies-entropy fluxes. There results imply the existence of a flux-splitting consistent with all of the entropy inequalities.