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.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

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.

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.

Klein-Gordon-Schrodinger方程的几种差分格式及比较

探究在特定的初值和边界条件下一维Klein-Gordon-Schrodinger方程的几种差分格式并进行比较。利用经典的向前差分算子、中心差分算子、Crank-Nicolson方法和紧差分算子分别为Klein-Gordon-Schrodinger方程构造向前Euler式、Crank-Nicolson格式及紧差分格式。结果表明:Crank-Nicolson格式及紧差分格式能够精确地保持离散电荷和能量守恒。数值实验验证了理论结果的正确性。

A Fourth-Order Unstructured NURBS-Enhanced Finite Volume WENO Scheme for Steady Euler Equations in Curved Geometries

In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.

Conical Sonic-Supersonic Solutions for the 3-D Steady Full Euler Equations

This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.

MULTIPLE STATIONARY SOLUTIONS OF EULER-POISSON EQUATIONS FOR NON-ISENTROPIC GASEOUS STARS

The motion of the self-gravitational gaseous stars can be described by the Euler-Poisson equations. The main purpose of this paper is concerned with the existence of stationary solutions of Euler-Poisson equations for some velocity fields and entropy functions that solve the conservation of mass and energy. Under different restriction to the strength of velocity field, we get the existence and multiplicity of the stationary solutions of Euler-Poisson system.

CONSTRAINT VIOLATION STABILIZATION OF EULER-LAGRANGE EQUATIONS WITH NON-HOLONOMIC CONSTRAINTS

Two constraint violation stabilization methods are presented to solve the Euler Lagrange equations of motion of a multibody system with nonholonomic constraints. Compared to the previous works, the newly devised methods can deal with more complicated problems such as those with nonholonomic constraints or redundant constraints, and save the computation time. Finally a numerical simulation of a multibody system is conducted by using the methods given in this paper.

EULER SCHEME AND MEASURABLE FLOWS FOR STOCHASTIC DIFFERENTIAL EQUATIONS WITH NON-LIPSCHITZ COEFFICIENTS

For a stochastic differential equation with non-Lipschitz coefficients, we construct, by Euler scheme, a measurable flow of the solution, and we prove the solution is a Markov process.

SERIES PERTURBATIONS APPROXIMATE SOLUTIONS TO N-S EQUATIONS AND MODIFICATION TO ASYMPTOTIC EXPANSION MATCHED METHOD

A method that series perturbations approximate solutions to N-S equations with boundary conditions was discussed and adopted. Then the method was proved in which the asymptotic solutions of viscous fluid flow past a sphere were deducted. By the ameliorative asymptotic expansion matched method, the matched functions, are determined easily and the ameliorative curve of drag coefficient is coincident well with measured data in the case that Reynolds number is less than or equal to 40 000.

Single-Step Arbitrary Lagrangian-Eulerian Discontinuous Galerkin Method for 1-D Euler Equations

We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity modified by some smoothing,which is found to con-siderably reduce the numerical dissipation introduced by Riemann solvers.The scheme preserves constant states for any mesh motion and we also study its positivity preservation property.Local grid refinement and coarsening are performed to maintain the mesh qual-ity and avoid the appearance of very small or large cells.Second,higher order methods are developed and several test cases are provided to demonstrate the accuracy of the proposed scheme.

带副翼的翼身组合体绕流的Euler和N-S方程解

将对接分区网格与分区求解算法结合 ,有效地求解了带副翼偏转的翼身组合体绕流的N S方程 .数值方法中选用VanLeer分裂格式离散无粘通量项 ,采用中心差分格式来离散粘性通量项 .分区交界面采用了一种满足通量守恒的内边界耦合条件 .数值算例表明该方法是求解带操纵面偏转的翼身组合体绕流的有效方法 .

非结构网格多段翼型绕流Euler及N-S方程数值模拟

描述二维非结构无粘网格和粘性网格的生成方法 ,对二维翼型绕流进行了 Euler及 N- S方程数值模拟。运用推进面方法 ,生成无粘网格 ,进一步结合推进层方法 ,生成粘性网格。应用中心有限体积法 ,采用 B- L代数湍流模型 ,完成了对多段翼型流场的 Euler方程和 N- S方程的数值模拟计算。计算与实验结果对比表明 ,所采用的网格生成技术及流场计算方法是正确。

基于Euler/N-S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jam eson的中心差分有限体积法求解Eu ler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合。

高阶熵条件格式下Euler方程与N-S方程的混合算法

针对在工程问题中使用N S方程计算量大的问题 ,提出只在粘性作用显著的局部区域求解N S方程 ,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题 ,并提出了一种基于高阶熵条件格式的算子分裂算法 .应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算 .计算结果与试验数据的对比表明 ,所提出的方法是切实可行的 .

GLOBAL EXISTENCE AND ASYMPTOTIC BEHAVIOR FOR THE 3D COMPRESSIBLE NON–ISENTROPIC EULER EQUATIONS WITH DAMPING

We investigate the global existence and asymptotic behavior of classical solutions for the 3D compressible non-isentropic damped Euler equations on a periodic domain. The global existence and uniqueness of classical solutions are obtained when the initial data is near an equilibrium. Furthermore, the exponential convergence rates of the pressure and velocity are also proved by delicate energy methods.

UNCONDITIONAL STABLE SOLUTIONS OF THE EULER EQUATIONS FOR TWO－AND THREE－D WINGS IN ARBITRARY MOTION

The work presented here shows the unsteady inviscid results obtained for the twoand three-dimensional wings which are in rigid and flexible osciliations.The results are generated by a finite volume Euler method. It is based on theRunge- Kutta time stepping scheme developed by Jameson et al.. To increase the timestep which is limited by the stability of Runge-Kutta scheme, the implicit residualsmoothing which is modified by using variable coefficients io prerent the loss of flowphysics for the unsteady flows is engaged in the calculations. With this unconditionalstable solver the unsteady flws about the wings in arbitrary motion can be receivedefficiently.The two- and three-dimensional rectangular wings which are in rigid andflexible pitching oscillations in the transonic flow are invesigated here, some of thecomputational results are compared with the experimental data. The influence of thereduced frequency for the two kinds of the wings are researched. All the results givenin this work are reasonable.

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.

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.