• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    THE RIEMANN PROBLEM FOR ISENTROPIC COMPRESSIBLE EULER EQUATIONS WITH DISCONTINUOUS FLUX*

    2024-03-23 08:06:06孫印正屈愛芳
    關(guān)鍵詞:海棠

    (孫印正) (屈愛芳)

    Department of Mathematics, Shanghai Normal University, Shanghai 200234, China E-mail: sunyz1742@163.com; afqu@shnu.edu.cn

    Hairong YUAN (袁海棠)

    School of Mathematical Sciences and Shanghai Key Laboratory of Pure Mathematics and Mathematical Practice, East China Normal University, Shanghai 200241, China E-mail: hryuan@math.ecnu.edu.cn

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

    Key words compressible Euler equations; Riemann problem; Radon measure solution;delta shock; discontinuous flux; wave interactions

    1 Introduction

    For many applications in mathematical analysis and physical modeling,considerable attention has been paid to the Euler equations.The one dimensional isentropic compressible Euler equations read as

    for (x,t)∈R×[0,∞), andρstands for the mass density,ufor the velocity, andpfor the pressure.For a two-phase flow consisting of pressureless flow and polytropic flow separated by a free boundaryx=x(t), the state equation can be written as

    where 1<γ <3 for most gases is the adiabatic exponent.If there is concentration of mass onx=x(t), that is called a delta shock.Since the flux is discontinuous along the curvex=x(t), system (1.1), together with (1.2), comprises a typical system of conservation laws with a discontinuous flux.

    So far,much effort has been devoted to studying conservation laws with discontinuous flux.This is especially so for the scalar case, on account of its applications to a wide range of fields including traffic flow, the sedimentation of suspensions, elasticity, multi-phase flow in porous media, polymer flooding, etc.A feature of these works is that the discontinuity of the flux is fixeda priori.More precisely,researchers have focused mainly on the Cauchy problem of scalar equations like

    whereF(x,u) =H(x)f(u)+(1-H(x))g(u), and whereH(x) is the Heaviside function.The existence of weak solutions has been obtained by various methods, including, for example,vanishing viscosity,front tracking,and numerical schemes.These weak solutions are,in general,Lebesgue integrable functions; see [1-3, 5, 6, 13] and the references cited therein.Recently Aektaet al.[3]established a solution of the scalar conservation law(1.3)with this concentration supported on a fixed discontinuity of flux.Fewer results for the system of conservation laws with discontinuous flux have been obtained.Fazio and LeVeque[9], and LeVeque[17, section 14.14],discussed a multifluid problem for non-isentropic compressible Euler equations by numerical simulation, where the interface between two distinct polytropic gases is a contact discontinuity.Readers may catch a glimpse of the system of conservation laws with discontinuous flux in[12, 16, 27].As far as we know, there is no result on the discontinuous flux problem for isentropic flow with concentration.

    Usually,the failure of strict hyperbolicity of the system in some regions may lead to the loss of the spatial total variation bound for some conserved quantities, which makes the problem sophisticated.For example, the adhesion property of the pressureless flow may lead to mass concentration.To characterize the concentration phenomenon of system (1.1), one needs an appropriate definition of a measure solution rather than a Lebesgue integrable weak solution.On the other hand, choosing an appropriate entropy condition to deal with the non-uniqueness of weak solutions is also a thorny issue.Several criteria have been proposed for this, depending on the flux function; see,for example,[4,7,11,29].Among these criteria,the over-compressing condition is usually used in the construction of solutions with concentration for other conservation laws [22].A delta shock of (1.1), (1.2) is said to satisfy the over-compressing condition,if

    whereul(t)=u(x(t)-,t),ur(t)=u(x(t)+,t), and wherex=x(t) is the delta shock front.

    Recently,Jin,Qu and Yuan[15]studied the singular Riemann problem for the Euler system of polytropic gases.They admitted the Radon measure solution proposed by Qu, Yuan and Zhao [25] for a steady Euler system with a general state equation, which generalizes previous definitions of measure solutions of the Euler system for the pressureless case, and they constructed non-unique piecewise constant solutions of the classical Riemann problem under the over-compressing condition, as well as certain existence results for singular Riemann problems.The discontinuity of the flux that might support concentration is determined together with the gas [18]; see [10, 14, 23-26] for more physical background.This means that one needs to propose new criteria on the discontinuity of the flux in order to ensure a unique solution.

    In this paper, we consider the singular Riemann problem for the conservation laws with discontinuous flux(1.1), (1.2).The initial data for the singular Riemann problem is formulated as

    whereρ0is a nonnegative number.Notice that the singular Riemann problem(1.1),(1.2),(1.6)is consistent with the classical Riemann problem in the caseρ0=u0= 0, while it is totally different in the caseρ0>0.To study this singular Riemann problem, we admit the concept of a Radon measure solution of the Euler system.Besides the global existence of infinite solutions of the classical Riemann problem satisfying the over-compressing condition (1.5) for any initial left and right states, we also establish the global existence of infinite solutions forρ0>0 with any initial left and right states under the over-compressing condition (1.5).To get the unique solution, we introduce some new criteria on the discontinuity of the flux based on the natural connection of the singular Riemann problem and the piston problem (cf.[10]); that is, we have the following:

    Sticking and Slipping Condition (abbreviated as SS condition): when the pressureless gas flows into the delta shock, it sticks to the delta shock and becomes a part of it, while the polytropic gas on the right of the delta shock has the same speed as it.

    We obtain the unique local piecewise smooth solution as well as its clear geometrical structure under the SS condition for a large range of initial data.The global existence of such solutions is considered in a future work.For initial data not in the relevant range, further exploration is required, for various technical reasons.

    As an application of the singular Riemann problem, we consider a free piston problem.Initially, there is a free piston with massρ0and velocityu0, and it separates the uniform pressureless gas on the left from the uniform polytropic gas on the right.The piston is accelerated by the difference of the forces acting on each side.We prove that the existence of this piston problem can be reduced to the existence of a singular Riemann problem, and therefore we get the local existence of a piecewise smooth solution of the piston problem.

    For the free piston problem with more general initial data, Liu [20] obtained a global solution with a small bounded variation by studying two boundary value problems of the one dimensional full Euler equations in Lagrangian coordinates, coupled with an ordinary differential equation satisfied by the piston.Takeno [28] considered a free piston problem for onedimensional isentropic gas dynamics, with the state equation on each side of the piston being a different polytropic gas, by the method of compensated compactness.More results on piston problems can be found in the review article of Ding and Li [8].

    The main contributions of this paper include, first, showing that the discontinuity of the flux is free, which is determined together with the solution, rather than being fixeda priorias was done in previous works.Second, the existence of multiple solutions satisfying the over compressing condition for the singular Riemann problem is demonstrated, and solutions with more abundant wave structures are found, in particular, the structure composed of a delta shock and a rarefaction wave.Third, the interaction of a rarefaction wave and a delta shock is analyzed.Fourth, we propose an additional physical boundary condition on the discontinuity of the flux function to obtain the unique solution; this has applications for studying classical piston problems.Finally, we establish a direct relation between the singular Riemann problem and the classical fluid-structure coupling problem, and get the local existence as well as the structure of a classical piston problem.

    The organization of the rest of this paper is as follows: in Section 2,we propose a definition of Radon measure solutions for the Riemann problem of the conservation laws with discontinuous flux (1.1), (1.2), (1.6).The generalized Rankine-Hugoniot conditions of the delta shock are also derived.In Section 3, we construct the Radon measure solutions of the Riemann problem(1.1), (1.2), (1.6) withρ0=0; these consist of a delta shock and classical waves in a polytropic gas.The structure of solutions depends on the relations of the initial states and whether the over-compressing condition (1.5) or the SS condition is under consideration.Solutions consisting of three waves are also considered.In Section 4, we focus onρ0>0.Because of the initial mass concentration, we have to tackle the problem of the interactions of delta shocks and classical waves, especially rarefaction waves, and thus consider the existence of compound waves.Whether these solutions satisfy the over-compressing condition is also verified.In Section 5,under the SS condition, by regarding the delta shock as a free boundary, the local existence of a unique piecewise smooth solution is established.As an application of the singular Riemann problem, we consider, in the end, some related free piston problems, which leads to the equivalence of delta shocks and free pistons.We briefly summarize the results and propose some unsolved problems in a short Section 6.

    2 Definition of Radon Measure Solutions and Delta Shock

    We now present the definition of Radon measure solutions for the Riemann problem (1.1),(1.2), (1.6).

    2.1 Definition of Radon Measure Solutions

    LetL1be the Lebesgue measure, andM(R) the space of Radon measures on the real line R.Consider a mapping

    is continuous fort.We denote the collection of all these mappings byC([0,∞);M(R)).In particular, we say thatm(t) is a Dirac measure with weightωL(t) supported on a Lipschitz curveL.={x=x(t)|t ∈[0,T)}, and denoted specifically asωL(t)δx(t), ifωL(t) is continuous,and if

    In the sequel,C1candCcdenote, as usual, the set of continuously differentiable and continuous functions with compact supports, respectively.The derivative of a functionφ(x,t)with respect totis written asφt, and similarly,φx.=?xφ.We also use the standard notationm ??to mean that the measuremis absolutely continuous with respect to the measure?.

    Definition 2.1 Let? ∈C([0,∞);M(R)), and letu(·,t) be a?(t)-measurable function on R for allt ≥0.The pair (?,u) is called a Radon measure solution of the Riemann problem(1.1), (1.2) and (1.6) if there existm,n,?∈C([0,∞);M(R)) such that

    1) for any test functionφ ∈C1c(R2), it holds that

    3) the classical Lax entropy condition is valid for all discontinuities in the flow field (ρ,u)in the space-time region, where?(t)?L1.

    Obviously, if(?(t),u(·,t))is a Radon measure solution to the Riemann problem(1.1), (1.2)and (1.6), then the initial data (1.6) should be written rigourously as

    whereL1is the Lebesgue measure restricted on anL1-measurable setA, and IAis the indicator function of the setA.

    2.2 Delta Shock

    We call a Radon measure solutiona delta shockif it is of the form

    with Ωl.={x ∈R|x <x(t)}, Ωr.={x ∈R|x >x(t)},x(0) = 0, and where (ul,ρl) (resp.(ur,ρr)) is a smooth solution of (1.1) in Ωl(resp.Ωr) with the state functionp= 0 (resp.p(ρ) =ργ).We now consider the Radon measure solutions such that the initial left state and right state are connected by only a single delta shock.Whether these solutions satisfy the over-compressing condition (1.5) or the SS condition is to be verified.

    As usual,we start with the generalized Rankine-Hugoniot conditions along the discontinuity of the flux and get the delta shock curve.By substituting(2.9)into(2.4),we obtain the following jump conditions along the discontinuityx=x(t):

    Here [ρ].=[ρ](x(t),t)=ρ(x(t)+,t)-ρ(x(t)-,t).Hence, we have that [ρ]=ρr-ρl.These four equations in (2.10), together with the Radon-Nikodym derivatives’ relations (2.5),

    Then, the delta shock is determined and given by (2.12) and (2.13), once the front of delta shockx=x(t) is solved.

    We go on to determine the curvex=x(t).Equations (2.11) and (2.13) give that

    One may refer to [15, Lemma 2] for more detailed calculations, though keep in mind that here the coefficient of[ρ]inaispr,rather than[p],as it is in that paper.This difference creates a big difference in the analysis.

    Now, we have the following conclusion regarding the existence ofx=x(t):

    Lemma 2.2 Under the condition thatwρ(t), given by (2.13)1, is nonnegative, we have the following solvability conclusion regarding the equation (2.15) with respect tox(t):

    It follows from this lemma that, for the Riemann problem (1.1), (1.2), (1.6), ifρ0=0, the delta shock front is a straight line.For [ρ]=0, the difference of pressure only affects the speed of the wave, and the difference of speed affects both the weight and the speed of the wave.For[ρ]/=0, both the speed and the weight of the wave depend on all of the initial data.

    The proofs of these two lemmas can be obtained by some tedious and straightforward calculations similar to those presented in [15, Lemma 2, Propositions 1, 2, 3] (with some minor modifications), so we omit them here.Special attention should be given to the case [ρ] = 0,[u]=0 andρ0=0, which makes the equation (2.14) unsolvable; that is, the solution consisting only of a single delta shock does not exist due to the pressure of the polytropic gas.This is another difference from [15].The main results of Lemmas 2.3 and 2.4 are illustrated in Table 1.

    Remark 2.5 Fix(ul,ρl)withρl ≥0.In view ofa(ul,ρl;u,ρ)=ρlρ(u-ul)2-ργ(ρ-ρl),it is easy to check that

    Table 1 Delta shock solution

    For the classical Riemann problem (ρ0= 0), it follows from (2.18) that, forρl= 0, the curvea= 0 reduces to theu-axis.Forρl >0, the monotonicity and convexity of these two curves can be analyzed by standard differentiation.

    As shown in Lemmas 2.3 and 2.4,for a given initial left stateUlin the phase(u,ρ)-plane,the domain to whichUrbelongs such that the solution of the classical Riemann problem consists of a single delta shock is the region{(u,ρ)|a ≥0}{ρ=ρl,u >ul}.Meanwhile, for the singular Riemann problem, whenUrbelongs to a portion of the region{(u,ρ)|a ≥0}, the delta shock exists globally.Otherwise, the delta shock exists locally.This feature may lead,significantly, to the possibility of multi-solutions in the construction.Additional conditions,according to certain physical mechanisms, should be proposed in order to ensure a unique solution.In the sequel, we will see that to construct a global solution to the singular Riemann problem without any additional condition on the delta shock is straightforward.Once extra conditions are proposed on the delta shock, however, the problem becomes quite complicated.We will discuss the problem under respectively the over-compressing condition (1.5) and the SS condition.

    2.2.1 Discontinuity of Flux That Satisfies the Over-compressing Condition (1.5)

    In what follows, we check the over-compressing condition (1.5) for the delta shocks constructed in Lemmas 2.3 and 2.4.Note thatbcan be rewritten as

    Thus, a necessary condition for the over-compressing (1.5) is thatb ≥0.That is to say, whenb <0, the over-compressing condition (1.5) initially fails.Also,asatisfies that

    Then,a ≤ρ2r[u]2when [ρ]≥0, anda ≥ρ2r[u]2when [ρ]≤0.By the relations amongx′(t),ulandur, we then have

    Lemma 2.6 Considering the Radon measure solutions of the singular Riemann problem(1.1), (1.2) and (1.6) obtained by Lemmas 2.3 and 2.4, we have the following results:

    1) forρ0=0, the solutions obtained by Lemma 2.3 satisfy the over-compressing condition(1.5) if one of the following conditions is satisfied:

    (i) [ρ]=0, [u]≤0,pr ≤ρl[u]2,

    (ii) [ρ]>0, [u]≤0,a ≥ρ2l[u]2,

    (iii) [ρ]<0, [u]≤0, 0≤a ≤ρ2l[u]2;

    2) forρ0>0, the solutions obtained by Lemma 2.4 satisfy the over-compressing condition(1.5) globally if one of the following conditions is satisfied:

    (i) [ρ]=0,ur ≤u0≤ul,pr ≤ρl[u]2,

    (ii) [ρ]>0,ur ≤u0≤ul,a ≥ρ2l[u]2,

    (iii) [ρ]<0,ur ≤u0≤ul, 0≤a ≤ρ2l[u]2;

    3) forρ0>0, the solutions obtained by Lemma 2.4 satisfy the over-compressing condition(1.5) locally if one of the following conditions is satisfied:

    (i) [ρ]=0,ur ≤u0≤ul,pr >ρl[u]2,

    (ii) [ρ]>0,ur ≤u0≤ul,a <ρ2l[u]2,

    (iii) [ρ]<0,ur ≤u0≤ul,a >ρ2l[u]2.

    Lemma 2.6 implies that the curvea=ρ2l[u]2in the phase (u,ρ)-plane also plays an important role in constructing solutions; this will be specified in Subsection 3.2.

    2.2.2 Discontinuity of Flux That Satisfies SS Condition

    We go on to check the SS condition for the delta shocks constructed in Lemmas 2.3 and 2.4.We have the following conclusions:

    Lemma 2.7 (1)Forρ0=0,the problem(1.1),(1.2)and(1.6)admits a global delta shock with the discontinuityx=x(t) satisfying thatx′(t)=urif and only if

    Proof (1) Assume that the initial left and right states are connected by a single delta shockx=x(t).For the classical Riemann problem, that is,ρ0=0, when [ρ]=0, [u]≤0, and by (2.16), the relationx′(t)=urgives that

    withur ≤ul.By replacing the expressionadefined in (2.18), we have that

    3 Radon Measure Solutions of the Classical Riemann Problem (ρ0 =0)

    In this section, we investigate the Radon measure solutions of the classical Riemann problem (1.1), (1.2) and (1.6), that is,ρ0= 0, for general initial left and right constant data(ρl,r,ul,r).As is shown in Section 2, the initial left and right constant states cannot always be connected by just a single delta shock globally, or even locally.To construct the solution, more elementary waves should be taken into account.We search now for solutions with elementary waves connected by constant states.Notice that the flow on the left of the discontinuity of the flux is pressureless, while it is a polytropic gas on the right.Recall that for the pressureless gas, the elementary waves are a contact discontinuity and a delta shock; for the polytropic gas,the elementary waves are the classical shock wave, the rarefaction wave and the delta shock,as indicated in [15].Since this still lacks a general entropy condition to ensure the unique solution, we search for the solutions with the least number of waves in the construction of the Riemann problem first for all initial data.Also, we will give some examples which show that,for some initial data, there exist multiple solutions consisting of different numbers of waves.As mentioned above, the over-compressing condition (1.5) and the SS condition are addressed to ensure uniqueness in the construction of solutions consisting of delta shocks.We will further pick up those solutions under either of these two conditions.

    3.1 Multi-solution of the Classical Riemann Problem

    Lemma 2.3 shows that, when the initial dataUl,rsatisfies one of (i) [ρ] = 0, [u]<0; (ii)[ρ]<0; (iii) [ρ]>0,a ≥0, it can be connected by the discontinuity of the flux with the concentration of mass and momentum.What left is to consider the cases (iv) [ρ]>0,a <0;(v) [ρ]=0, [u]≥0.Recall that (see the appendix in [15]), for the polytropic gas, the set of left states that can be connected to the right stateUrby a forward shock or a forward rarefaction wave in state space are, respectively,

    We consider first the case for the initial data satisfying that [ρ]>0,a <0.We seek solutions that consist of a delta shock followed by an elementary wave in polytropic flow, and state the following:

    Lemma 3.1 Consider the classical Riemann problem (1.1), (1.2), (1.6) withρ0= 0.If the initial states satisfy that [ρ]>0,a <0, then

    (i) there exist infinitely many Radon measure solutions consisting of a delta shock followed by a forward rarefaction waveR2;

    (ii) constructing Radon measure solutions consisting of a delta shock followed by a forward shockS2is inappropriate;

    (iii) neither constructing solutions consisting of a delta shock followed a backward rarefaction waveR1, nor a delta shock followed by a backward shockS1, is appropriate.

    which is consistent with (3.7).

    For [ρ]<0, we can get an inequality (3.7) similar to that for [ρ]>0.By standard calculations, we solve (3.7) and obtain that

    which is unsolvable forρ >0, and therefore illustrates our claim.

    that is, the curvesC2andD2are tangential atP1.

    LetQ1be the domain surrounded by the curvesC1,C2,D1andD2.Then, forUr ∈{(u,ρ)|a(ul,ρl;u,ρ)<0}, anyUm ∈R-12(Ur)∩(Q1{ρ=ρl,u ≥ul}) can be chosen as intermediate states such thatUlcan be connected toUmby a delta shock andUmcan be connected toUrby a forward rarefaction waveR2.In addition,we may check,without any difficulty,thatx′(t)≤λ2(Um) whenρm=0,um ≥ul; see Figure 1.

    Figure 1 When Ur ∈{(u,ρ)|a <0}, there exist Radon measure solutions consisting of a delta shock followed by a forward rarefaction wave R2, with the intermediate states being in the region surrounded by C1, C2, D1 and D2, as well as the u-axis.While constructing a solution consisting of a delta shock and a forward shock S2, is inappropriate.

    At the end of the proof, we show the inappropriateness of constructing solutions consisting of a delta shock followed by a backward shockS1forUr ∈{(u,ρ)|a(ul,ρl;u,ρ)<0}.Assume,to the contrary, thatUlcan be connected to an intermediate stateUmby a delta shock and thatUmcan be connected toUrby a backward shockS1.Then

    which contradicts (3.24).Then, we claim that the curvesD′1andD′2must be in the interior of the region{(u,ρ)|a(ul,ρl;u,ρ)<0}, which implies that solutions consisting of a delta shock and a backward shockS1do not exist because of the non-existence of admissible intermediate states.This completes the proof.□

    Consider now the case for the initial data satisfying that [ρ] = 0, [u]≥0.We go on to seek singular Riemann solutions that consist of a delta shock followed by elementary waves in a polytropic flow.We have the following conclusions:

    Lemma 3.2 Consider the classical Riemann problem (1.1), (1.2), (1.6) withρ0= 0.If the initial states satisfy that [ρ]=0, [u]≥0, then

    (i) there exist Radon measure solutions consisting of a delta shock followed by a forward rarefaction waveR2;

    (ii) there exist Radon measure solutions consisting of a delta shock followed by a backward rarefaction waveR1;

    (iii)neither constructing solutions consisting of a delta shock followed by a backward shockS1nor a delta shock followed by a forward shockS2is appropriate;

    (iv) there exist Radon measure solutions with three waves.

    Proof (i) For givenUl,rsatisfying that [ρ]=0, [u]≥0, these two states can be joined by a delta shock and a forward rarefaction waveR2.In fact, for anyUm ∈R-12(Ur)∩(Q1∪{ρ=0,u ≥ul}),Ulcan be connected toUmas above by a delta shock, andUmcan be connected toUrby a forward rarefaction waveR2; see Figure 2.More precisely, assume that the forward rarefaction wave curveR2(ul,0) passes through (ul,0) on the left, and intersects the curve{ρ=ρl,u ≥ul}at ~U=(~u,ρl), where

    Then, ifur ≥~u, the left stateUlis in fact connected uniquely to the intermediate stateUm=(um,0) by a contact discontinuity, which can be regarded as a delta shock with a weight of zero.Ifur <~u, the intermediate statesUmare neither vacuums, nor unique.

    Figure 2 When the initial states satisfy[ρ]=0,[u]≥0,there exist Radon measure solutions consisting of a delta shock followed by a forward rarefaction wave R2,with the intermediate states being in the region surrounded by C1, C2, D1 and D2, as well as the u-axis.

    (ii) To construct solutions consisting of a delta shock followed by a backward rarefaction waveR1, similarly to Lemma 3.1, we determine the region to which the intermediate stateUmbelongs by proving thatx′(t)≤λ1(Um), whereλ1=u-c.

    Figure 3 When the initial states satisfy[ρ]=0,[u]≥0,there exist Radon measure solutions consisting of a delta shock followed by a backward rarefaction wave R1, with the intermediate states being in the region surrounded by C1, C2, D3 and D4.

    (iii) Suppose, to the contrary, thatUmis an intermediate state such thatUlis connected toUmby a delta shock andUmis connected toUrby a forward shockS2.Then

    which contradicts (3.34).Therefore, there is no solution consisting of a delta shock and a forward shockS2because of the non-existence of an admissible intermediate state.

    (iv) It is easy to check that when the initial states satisfy that [ρ]=0 and [u]≥0, we can choose a statem ∈R-12(Ur) such thatUlandare connected by a delta shock followed by a backward rarefaction waveR1, andmis connected toUrby a forward rarefaction waveR2;that is, this solution contains three waves; see Figure 4.This completes the proof.□

    Figure 4 Solutions with three waves can be constructed by choosing the intermediate state Um as belonging to the region Q2, which is surrounded by curves C1, C2, D3 and D4.

    Note that, so long as the intermediate stateUmthat we have chosen belongs to the regionQ2, solutions containing three waves can be constructed for general initial data rather than satisfying that [ρ] = 0, [u]≥0.This fact also implies that the statesUlandcan be connected by a delta shock followed by a backward rarefaction waveR1as well, in addition to the single delta shock constructed in Lemma 2.3.

    Remark 3.3 Ifρl= 0, the curve{(u,ρ)|a(ul,ρl;u,ρ) = 0}degenerates to theu-axis.Lemma 3.1 implies that there exists a unique Radon measure solution consisting only of a forward rarefaction waveR2for the classical Riemann problem (1.1), (1.2), (1.6) withρ0= 0 andρr >0.

    Combining Lemmas 2.3, 3.1, 3.2, we arrive at

    Theorem 3.4 The classical Riemann problem(1.1),(1.2),(1.6)withρ0=0 admits global solutions forul,ur ∈R andρl,r ≥0.Moreover, we have that

    (i)when the initial dataUl,rsatisfies one of(a)[ρ]=0,[u]<0;(b)[ρ]<0;(c)[ρ]>0,a ≥0,then there exists a solution consisting only of the discontinuity of the flux with concentration of mass and momentum(i.e., delta shock).There also exist infinitely many solutions consisting of a delta shock followed by a backward rarefaction waveR1for (a) and (b);

    (ii) when the initial states satisfy that [ρ] = 0, [u]≥0, then there exist infinitely many Radon measure solutions consisting of a delta shock followed by a forward rarefaction waveR2,or of a delta shock followed by a backward rarefaction waveR1.Solutions consisting of three waves can be constructed as well;

    (iii) when the initial states satisfy that [ρ]>0,a <0, then there exist infinitely many Radon measure solutions consisting of a delta shock followed by a forward rarefaction waveR2.

    3.2 Multi-solution Under Over-compressing Condition

    As indicated by Theorem 3.4,the uniqueness of a solution to the classical Riemann problem is lost, due to the appearance of delta shocks.We go on in this subsection to single out those solutions constructed in Theorem 3.4 that satisfy the over-compressing condition (1.5).

    For given left state(ul,ρl),we divide the upper(u,ρ)-plane into several parts by the curves defined in Section 2, according to the results of Lemma 2.6.In addition, define that

    By direct calculations, we obtain that, for a given left stateUlwithρl >0, the curvesC1,E1do not intersect whenρ >0, and so neither do the curvesC2,E2; see Figure 5.Forρl=0, the curvesC1andE1, as well asC2andE2, coincide, which is consistent with Remark 3.3.

    Figure 5 When Ur ∈I, there exists a Radon measure solution consisting of a single delta shock.When Ur ∈II, there exist Radon measure solutions consisting of a delta shock followed by a forward rarefaction wave R2.When Ur ∈III, there exists a unique Radon measure solution consisting of a delta shock followed by a forward rarefaction wave R2.

    We introduce some notations for the partitions as follows:

    Theorem 3.5 The classical Riemann problem(1.1),(1.2),(1.6),for whichρ0=0,admits global solutions that satisfy the over-compressing condition (1.5) forul,r ∈R andρl,r ≥0.Moreover, we have that

    (i) ifUr ∈I, there exists a Radon measure solution consisting of a single delta shock;

    (ii) ifUr ∈II, there exist infinitely many Radon measure solutions consisting of a delta shock followed by a forward rarefaction waveR2, with intermediate states lying in I1∪E1;

    (iii) ifUr ∈III, there exists a unique Radon measure solution consisting of a contact discontinuity followed by a forward rarefaction waveR2, with the intermediate state being a vacuum.

    Proof (i) For the case whenUr ∈I, the results have be proven in Lemma 2.6.

    (ii) IfUr ∈II, we now show the conclusion by constructing solutions consisting of a delta shock followed by a forward rarefaction waveR2.For the curveE1, we have that

    By continuity, the curvesE1andR-12(u,ρr) must intersect at a unique point ~U.It is easy to check thatx′(t)≤~u+~c.This implies the existence result in region II.

    We now determine the region to which the intermediate stateUmbelongs,under the further requirement thatx′(t)≤um+cm.As was discussed in Lemma 3.1, this region is precisely that surrounded by the curvesD1andD2in region I.Solving the equations

    so that these solutions satisfy the over-compressing condition (1.5) as well.□

    Remark 3.6 The proof of the above theorem also implies that under the over-compressing condition (1.5), the structure of the classical Riemann solution is unique, in the sense that the solutions consist of the same types of waves,but these may differ from each other in terms of the strength of the corresponding waves, so that the intermediate state of the solution may be nonunique.In addition, solutions with three waves do not exist here, since the over-compressing condition(1.5)rules out the possibility of wave structures that consist of a delta shock followed by a backward wave in the polytropic flow.

    3.3 Unique Solution Under SS Condition

    In this subsection,we continue to consider the classical Riemann problem(1.1),(1.2),(1.6)withρ0=0, and present the following conclusion:

    Theorem 3.7 Under the SS condition, the classical Riemann problem (1.1), (1.2), (1.6)withρ0= 0 admits a unique global Radon measure solution consisting of waves that connect piecewise constant states.More precisely,

    (i) ifUr ∈E1, there exists a unique Radon measure solution consisting of a single delta shock;

    (ii) ifUr ∈I, there exists a unique Radon measure solution consisting of a delta shock followed by a forward shockS2;

    (iii) ifUr ∈II∪III, there exists a unique Radon measure solution consisting of a delta shock followed by a forward rarefaction waveR2.

    Proof For any given left stateUl, denote that

    which is precisely the curveE1defined in (3.38).Then the curveF(Ul) andR2(ul,0) divide the upper-half state plane into the same three parts as those specified in (3.39).

    IfUr ∈F(Ul), it follows from Lemma 2.7 thatUrcan be connected to the left stateUlby a delta shockx=urt.

    IfUr ∈I, the forward shock curveS-12(Ur) viaUron the right will intersect withF(Ul) at a unique point (um,ρm).Thus,Ulcan be connected toUmby a delta shockx=umt, whileUmcan be connected toUrby a forward shockS2.

    IfUr ∈II, the forward rarefaction wave curveR-12(Ur) viaUron the right will intersect withF(Ul) at a unique point (um,ρm).Thus,Ulcan be connected toUmby a delta shockx=umt, whileUmcan be connected toUrby a forward rarefaction waveR2.

    IfUr ∈III, the curveR-12(Ur) viaUrwill intersect with theu-axis at a unique pointUm= (0).Thus,Ulcan be connected to the vacuum stateUmby a contact discontinuityx=ult, whileUmcan be connected toUrby a forward rarefaction waveR2.□

    4 Multi-solution of the Singular Riemann Problem under Over-compressing Condition

    In this section, we consider the singular Riemann problem (1.1), (1.2), (1.6) withρ0>0.We first restrict ourselves to the solution with waves connected by constant states.Under the over-compressing condition(1.5)on a delta shock,we find that there is no solution consisting of waves connected by constant states for some initial data.Solutions with more abundant wave structures should be taken into consideration.

    4.1 Global Solutions with Waves Connecting Piecewise Constant States

    Note that we can always construct a local delta shock for the singular Riemann problem withρ0>0 which can be prolonged; we can do this by solving a classical Riemann problem withρ0= 0 at the time when the delta shock disappears.We have the following result about the existence of Radon measure solutions:

    Theorem 4.1 The singular Riemann problem (1.1), (1.2), (1.6) admits infinitely many global Radon measure solutions, the structure of which near the origin consists only of a delta shock.

    Figure 6 Global delta shock

    Figure 7 [ρ]=0, [u]>0

    Figure 8 [ρ]>0, a <0

    Figure 9 [ρ]/=0, a >0, b <-

    By Theorem 3.4 and Theorem 4.1, we know that the singular Riemann problem admits infinitely many solutions.Next, as was done in Section 3 forρ0=0, we would like to construct solutions satisfying the over-compressing condition (1.5).Then, interactions between the delta shock and the classical waves are inevitable.

    4.2 Multi-solution Under Over-compressing Condition

    In this subsection, we mainly consider the singular Riemann problem forρ0>0 under the over-compressing condition (1.5).Depending on the initial states on both sides of the concentration point that satisfies (1.5) or not, we divide the discussion aboutu0≤ur+crinto four cases.Case 1: the initial states satisfy the over-compressing condition (1.5); namely,ur ≤u0≤ul.Case 2:u0≤min{ul,ur}.Case 3:u0≥max{ul,ur}.Case 4:ul ≤u0≤ur.After that, we discuss the caseu0>ur+cr.

    4.2.1 The Initial States Satisfying the Over-compressing Conditionur ≤u0≤ul

    Recall that,for any non-vacuum left initial stateUl,the upper state space(u,ρ)-plane is divided into three parts, I, II, III, given in (3.39) by the curvesE1via (ul,0) andR2(ul,0).Moreover,I is divided into three parts too; see Figure 5.Denote by II1the region II∩{u ≤ul}, for simplicity.Asur ≤ul, the initial right stateUrbelongs either to II1or to I.

    If the initial data satisfies thatur ≤u0≤min{ul,ur+cr}andUr ∈I,then,by Lemma 2.6,the singular Riemann problem (1.1), (1.2), (1.6) admits a solution consisting only of a single delta shock.It also admits infinitely many solutions consisting of a delta shock and a forward rarefaction waveR2.

    Lemma 4.2 Under the over-compressing condition (1.5), with a given left stateUl, ifUr ∈II1andur ≤u0≤min{ul,ur+cr}, then the singular Riemann problem (1.1), (1.2), (1.6)admits infinitely many solutions consisting of a delta shock and a forward rarefaction waveR2,with the constant intermediate stateUm ∈I1satisfying thatum+cm ≥u0; see Figure 10.

    Proof ForUr ∈II1,one may find,by Lemma 2.6,an intermediate stateUm ∈R-12(Ur)∩I to get the local solution.The assumption thatur <ulimplies thatR-12(Ur)∩I/=?.Thus,anyUm ∈R-12(Ur)∩I satisfying thatu0≤um+cmcan be chosen as an intermediate state to connect the left stateUlby a delta shock and the right stateUrby a forward rarefaction waveR2at first.Since the delta shockx=x(t) is no longer straight, we need to check whether it will interact with the rarefaction waveR2later.It is enough to calculate the interaction time ofx=x(t) and the characteristicx=(um+cm)t.Supposing that they interact at some timet, then we have that, if [ρ]=ρm-ρl=0, the interaction timetshould satisfy that

    which, combined with (2.19), gives thatt=0, or that,

    which is consistent with (4.3) for [ρ]=0.

    Sinceu0≤um+cm, the numerator of (4.3) and (4.5) is non-positive.The denominator of(4.3) and (4.5) is non-negative if and only if

    It follows from (4.6) that, forUm ∈R-12(Ur)∩I1andu0≤um+cm, the value oftgiven by (4.3) and (4.5) is negative; that is, the solutions consist of a delta shock connectingUlandUm, followed by a forward rarefaction waveR2connectingUmandUr, without any interaction of waves (see Figure 10).This completes the proof.□

    Figure 10 Um ∈I1, u0 ≤um+cm

    Figure 11 Um ∈I1, u0 >um+cm

    Under the over-compressing condition(1.5),we consider global solutions in which the delta shock is coupled with a rarefaction wave.

    Lemma 4.3 Under the same assumptions as those used in Lemma 4.2, the singular Riemann problem (1.1), (1.2), (1.6) admits infinitely many global solutions.These solutions consist of a compound wave at first, which is the coupling of a delta shock and a forward rarefaction waveR2.In some solutions, these two waves will be separated from each other by an intermediate constant stateUm ∈I1satisfying thatum+cm ≤u0.In some other solutions,this kind of couple will never terminate in finite time, withx′(∞) =u∞+c∞satisfying thatU∞=(u∞,ρ∞)∈D1∪D2and thatu∞+c∞<u0; see Figures 11 and 12.

    Proof For any pointUm ∈R-12(Ur)∩I1, withu0>um+cm, although there is no local solution with the same structure as that in Lemma 4.2, we can construct a solution by first finding a compound wave of a delta shock and a forward rarefaction wave.

    Assume that the front of the delta shock in the compound wave isx=x(t).By the generalized Rankine-Hugoniot conditions (2.10), we have that

    whereu+(t),ρ+(t) are defined as in (4.1).Notice that the right side of the delta shock is assumed to be a centred rarefaction wave; it is a system of equations with unknown functionsωρ(t),ωm(t),x(t),u+(t) andρ+(t).We drop all of the subscripts “+”, without confusion.The first two equations of (4.7) yield that

    which determines the existence ofρ+(t)in(4.7).The other unknown functions are consequently determined by (4.7).

    We now determine the time when the delta shock separates from the rarefaction waveR2.Lettingρ=ρmin (4.17), we obtain that

    Therefore, whenUm ∈I1, the timet6above is positive.Whent >t6, it becomes a case ofu0<um+cm, as was discussed in Lemma 4.2.Thus, the delta shock and rarefaction wave are separated by a constant intermediate stateUm; see Figure 11.

    As the selection ofUm ∈R-12(Ur)∩I1tends to the boundaryD1andD2, the timet6defined by (4.18)tends to infinity,which implies that the coupling will never terminate in finite time; see Figure 12.□

    Figure 12 U∞∈D1 ∪D2,u0 >u∞+c∞

    Figure 13 Um ∈I2 ∪I3, u0 ≤um+cm

    We next search for solutions whereby a delta shock interacts with the forward rarefaction waveR2.

    Lemma 4.4 Under the same assumptions as those used in Lemma 4.2, the singular Riemann problem(1.1),(1.2),(1.6)admits infinitely many solutions consisting of a delta shock and a forward rarefaction waveR2.They are separated from each other at first by an intermediate constant stateUm ∈I2∪I3withum+cm ≥u0.After a period of time, they intersect and then interact with each other.This interaction will never terminate in finite time.The solution is illustrated in Figure 13.

    Proof For anyUm ∈R-12(Ur)∩(I2∪I3), as in Lemma 4.2, we first to show that there exists an intermediate constant state that connects the pressureless flow on the left by a delta shock satisfying the over-compressing condition and connects the uniform polytropic gas on the right by a forward rarefaction waveR2.We justify that these two waves intersect at a finite time.Similarly to the analysis in Lemma 4.2, the assumption thatur <ulimplies that(Ur)∩(I2∪I3)/=?.Ifum+cm ≥u0, then both (4.3) and (4.5) are positive, which means that the solutions consisting of a delta shock connectingUlandUm, followed by a forward rarefaction waveR2connectingUmandUr, exist in finite time, and these two waves are away from each other first, but they will intersect at the timet7given by (4.3) or (4.5).When these two waves meet, one must determine the interaction of a delta shock and a forward rarefaction wave.We now turn to solving this problem.

    Whent >t7, the delta shock and forward rarefaction waveR2interact with each other.By the generalized Rankine-Hugoniot conditions (2.10), the delta shock frontx=x(t) is given by

    whereτ=t-t7, and where (u1,ρ1) are the speed and weight of the delta shock att=t7.Similarly to Lemma 4.3, we drop all of the subscripts “+”, without confusion.Then, by (4.7)and (4.8), the first two equations of (4.19) give that

    Therefore, the last equation of (4.19) becomes

    Note that

    which determines the existence ofρ+(t)in(4.19).Consequently, the system of equations(4.19)is solved.

    We next check whether this interaction will terminate in finite time.It is enough to calculate the interaction time ofx=x(t) and the characteristic thatx=(um+cm)tfort >t7.Lettingρ=ρmandτ=t-t7in the (4.25), above, we have that

    which implies that the delta shock does not intersect the characteristic linex=(um+cm)tfort >t7, so this interaction will never terminate.The solution is illustrated by Figure 13.□

    Lemma 4.5 Under the over-compressing condition (1.5), with a given left stateUl, ifUr ∈II1and the initial velocities of the singular Riemann problem (1.1), (1.2), (1.6) satisfy thatul ≥u0>ur+cr, then this admits infinitely many solutions consisting of a local delta shock which then can be prolonged by a solution with the singular Riemann initial dataur=u0≤min{ul,ur+cr}.

    Proof WhenUr ∈II1,ul ≥u0>ur+cr,by Lemma 2.6,we construct a local delta shock int ∈[0,t8] such thatur=x′(t8)<ul, which becomes the case discussed above fort >t8.□

    Remark 4.6 The Lemmas 4.2-4.5 above show that by choosing a different intermediate stateUm,we can construct solutions with different wave structures for the same initial data.All of these solutions satisfy the over-compressing condition (1.5) on the delta shock.In addition,the proof of Lemma 4.3 implies that, so long as the intermediate stateUmsatisfies thatu0>um+cm,we can construct such compound wave solutions without the assumption thatu0≥ur.

    4.2.2 The Initial States not Satisfying Over-compressing Condition

    As discussed in Remark 4.6, whenρl=0, the region I degenerates into theu-axis, and the vacuum stateUm=(,0) must be taken into account, so that whenu0>, we can construct the compound wave solution as in Lemma 4.3.This fact will be used in the sequel.

    Lemma 4.7 Under the over-compressing condition (1.5), for the singular Riemann problem (1.1), (1.2), (1.6), with given left stateUl, if the initial velocity satisfies thatur ≤ul <u0,then it admits global solutions consisting of a delta shock,a contact discontinuity and a forward rarefaction waveR2.

    Proof For a given right stateUr, assume that the rarefaction wave curveR-12(Ur) in the phase (u,ρ)-plane intersects with theu-axis at (0) defined as above.Thenul ≥ur ≥and the left stateUlandUm= (,0) are connected by a contact discontinuityx=ultin a zero pressure flow, which is denoted byJ.Ifur ≤u0≤ur+cr, the vacuum stateUmand the right stateUrare connected by a delta shockx=x(t) coupled with a forward rarefaction waveR2.Meanwhile, this delta shock will intersect with the contact discontinuityx=ultin finite timet9, at whichu(x(t9)+,t9)≤x′(t9)≤ul.Thus this becomes the case discussed in Section 4.2.1 fort >t9;see Figure 14.Note that the solutions of the new singular Riemann problem att=t9constructed as in Section 4.2.1 are actually unaffected by the rarefaction waves connecting on the right of the state (u(x(t9),t9),ρ(x(t9),t9)).

    Figure 14 ur ≤ul <u0 ≤ur +cr

    Figure 15 ur ≤ul <u0

    Ifur ≤u0, we construct a solution withUmandUrbeing connected by a local delta shockx=x(t) fort ∈[0,t10], where

    which implies that the delta shock intersects with the contact discontinuityx=ultbefore the over-compressing condition (1.5) fails.Thus this also becomes the case discussed in Section 4.2.1 after the interaction timet10, sinceur ≤x′(t10)≤ul; see Figure 15.□

    Lemma 4.8 Under the over-compressing condition (1.5), for the singular Riemann problem(1.1),(1.2),(1.6),with a given left stateUl,if the initial velocity satisfies thatu0<ur <ul,then it admits global solutions consisting of a delta shock and a forward rarefaction waveR2.

    Proof Choosing the intermediate stateUm=(0),if≥u0,sinceul >u0,there exists,at least in finite time, a solution consisting of a delta shock connecting the statesUlandUmand a forward rarefaction waveR2connecting the statesUmandUr.The delta shock front is

    tends toulast →∞.Thus, as was discussed in Section 4.2.1, the delta shock and the rarefaction wave will intersect with each other in finite time, and this interaction will never terminate, although they do not interact initially; see Figure 13.

    Lemma 4.9 Under the over-compressing condition (1.5), for the singular Riemann problem (1.1), (1.2), (1.6), with a given left stateUl, if the initial velocity satisfiesu0≤ul ≤ur,then it admits global solutions consisting of a delta shock and a forward rarefaction waveR2.

    Lemma 4.10 Under the over-compressing condition(1.5),for the singular Riemann problem (1.1), (1.2), (1.6), with given left stateUl, if the initial velocity satisfies thatul <u0≤ur,then it admits global solutions consisting of a delta shock,a contact discontinuity and a forward rarefaction waveR2.

    Figure 16 u0 <

    Figure 17 ul ≤ ≤u0

    Lemma 4.11 Under the over-compressing condition(1.5),for the singular Riemann problem(1.1),(1.2),(1.6),with a given left stateUl,if the initial velocity satisfies thatul <ur <u0,then it admits global solutions consisting of a delta shock,a contact discontinuity and a forward rarefaction waveR2.

    By the discussion above, combining Lemmas 2.6 and 4.2-4.11, we arrive at the following conclusion:

    Theorem 4.12 Under the over-compressing condition (1.5), the singular Riemann problem (1.1), (1.2), (1.6) admits global Radon measure solutions for any initial data.These solutions are comprised of a delta shock followed by or coupled with a forward rarefaction waveR2if the initial data satisfies thatu0≤ul,while they comprise a delta shock followed by or coupled with a forward rarefaction waveR2and a contact discontinuity if the initial data satisfies thatu0>ul.

    Based on the analysis above,we give a corollary on a special case where the initial left state is a vacuum, which admits simpler wave structures.

    Corollary 4.13 Under the over-compressing condition(1.5), the singular Riemann problem(1.1),(1.2),(1.6)withρ0>0,ρl=0,ρr >0 admits global Radon measure solutions.More precisely,

    (i) ifu0≤the unique Radon measure solution consists of a global delta shock followed by a forward rarefaction waveR2, with the intermediate state being a vacuum;

    (iii)ifu0>ur+cr,the Radon measure solutions consist of a local delta shock fort ∈[0,t5]and a delta shock coupled with a rarefaction waveR2fort >t5.

    5 Unique Solution under the SS Condition and an Application to the Piston Problem

    As shown above, the over-compressing condition is not enough to single out a unique solution foru0>.Thus, we would like to seek solutions under the SS condition.As an application, we will illustrate the equivalence of the singular Riemann problem and a free piston problem by showing that the front of the delta shock and the trajectory of the piston are the same.Thus one could, for simplicity, use the former to solve the latter.

    5.1 Piecewise Smooth Solutions to the Singular Riemann Problem Under the SS Condition

    Lemma 2.7 implies that, generally, there is no solution with constant states connected by a delta shock, due to the arbitrariness ofu0.One needs to seek piecewise smooth solutions for this singular Riemann problem.

    Consider the caseu0≤with (,0) being the intersection ofR-12(Ur) and theu-axis.Iful <u0, there exists a global solution consisting of a contact discontinuityx=ultin a pressureless flow, a delta shock with the states on both sides of it being vacuums, and a forward rarefaction waveR2connecting the vacuum state on the left and polytropic gasUron the right.This solution is illustrated in Figure 16.

    Iful ≥u0, there exists, at first, a solution consisting of a delta shockx=x(t) connecting the pressureless flowUland the vacuum state and a forward rarefaction waveR2connecting the vacuum state and the polytropic gasUr.The velocity of this delta shock satisfies that

    Therefore, whenul ≤this solution is global, as is illustrated in Figure 10, withUmbeing a vacuum.Whenul >, this delta shock overtakes the forward rarefaction wave at

    After timet11, it becomes the case thatu0>, and the initial state is non-constant; this is complicated, and will be discussed in the future.We have therefore proved the following:

    Lemma 5.1 Under the SS condition, the singular Riemann problem (1.1), (1.2), (1.6)admits a unique global piecewise smooth solution for max{ul,u0} ≤, and a local piecewise smooth solution foru0≤<ul.

    and thus is uniquely determined.

    Define an angular domain

    withg(t)=(u1+c1)tandδbeing a sufficiently small positive constant.

    Now,the matter of determining a local classical solution of the singular Riemann problem is reduced to determining the classical solution of a free boundary problem of the Euler equations in ?1(δ) as follows:

    On the free boundaryx=x(t), we have the generalized Rankine-Hugoniot conditions (2.10),which are now reduced t o

    while on the other boundaryx=g(t), we have that

    We draw the following conclusion about this problem:

    Lemma 5.2 If the initial data satisfies thatρ0,1>0 andu0=u1,then the free boundary problem(5.5),(5.6),(5.7)admits a unique classical solution in ?1(δ)for sufficiently smallδ >0.

    Proof The Euler equations (5.5) can be rewritten in the characteristic form [19, p.6]

    Based on this lemma, we immediately have the following conclusion about the singular Riemann problem (1.1), (1.2), (1.6):

    Lemma 5.3 Under the SS condition, the singular Riemann problem (1.1), (1.2), (1.6)admits a unique local piecewise smooth solution foru0≤ur.

    Consider the caseur <u0<~ur, with ~urdetermined by (5.12) and (5.19).We seek a smooth solution in an angular domain surrounded by a delta shock and a forward shock wave.In this case, the polytropic gas is compressed by the delta shock on the left, and there may be a forward shock locally in front of the delta shock.We proceed to find this kind of solution.

    First of all,by the Rankine-Hugoniot condition,it is easy to get the shock solutionx=s′(0)twith left state being (~ρ1,u0), ~ρ1>ρrand right state being (ρr,ur) at the original point, by solving

    Define an angular domain

    where bothx=x(t) andx=s(t) are unknown curves to be determined.Then we formulate a free boundary problem in ?2(δ) with equations (5.5) in it, and the boundary conditions (5.6)onx=x(t), while we have the following boundary condition onx=s(t):

    Here [u] =ur-u(s(t)-,t), and so on.Also, the Lax entropy conditionρ(s(t)-,t)>ρris required on the shockx=s(t).

    In summary,by searching for the forward shock,we have assumed the following conditions:

    Here ~ρ1is determined uniquely by the first equation of (5.12).Under these assumptions, the second condition of (5.14) can be rewritten as

    which is attainable by the intermediate value theorem.When the equality (5.19) holds, the minimal characterizing number of matrixH2satisfies that|H2|min= 1.Thus, ifu0<~ur, we have that

    Therefore,it follows by[19,Theorem 3.1.1 and Theorem 3.5.1]that we have the existence of a uniqueC1solution to the free boundary problem (5.5), (5.6), (5.14) on ?2(δ) for sufficiently smallδ >0.Thus, we have

    Lemma 5.4 If the initial data satisfiesur <u0<~ur, with ~urdetermined by (5.12) and(5.19), then the free boundary problem (5.5), (5.6), (5.14) admits a unique classical solution in?2(δ) for sufficiently smallδ >0.

    Based on this lemma, we immediately have the following conclusion about the singular Riemann problem:

    Lemma 5.5 Under the SS condition, the singular Riemann problem (1.1), (1.2), (1.6)admits a unique local piecewise smooth solution forur <u0<~ur.

    Combining Lemmas 2.7, 5.1, 5.3 and 5.5, we arrive at

    Theorem 5.6 Under the SS condition, the singular Riemann problem (1.1), (1.2), (1.6)admits a unique local piecewise smooth solution foru0<~ur.More precisely, this solution consists of a delta shock followed by a forward rarefaction wave ifu0≤ur, and by a forward shock wave if ~ur >u0>ur.

    In the above discussions, the SS condition requires thatul ≥u0, i.e.the pressureless gas flows into the delta shock.Otherwise, whenul <u0, we add a contact discontinuityx=ulton the left of the delta shock, with the state between these two waves being a vacuum.Then, we obtain the unique local piecewise smooth solution for any initial data by the following corollary:

    Corollary 5.7 Under the SS condition, the singular Riemann problem (1.1), (1.2), (1.6)withρl=0 admits a unique local piecewise smooth solution foru0<~ur.Moreover,the solution consists of a delta shock followed by a forward rarefaction wave ifu0≤ur, and by a forward shock wave if ~ur >u0>ur.

    5.2 Application to the Piston Problem

    In this subsection, we consider a free piston problem.More precisely, letting thex-axis be a tube with constant cross sections,there is a piston with massm0>0 moving in the tube that separates the tube into two parts.There is a pressureless gas in the left part and a polytropic gas in the right.Without loss of generality, we suggest that the thickness of the piston is zero,so it reduces to a point mass.Att= 0, assume that the piston lies atx= 0, and that the states of gas on both sides are uniform, namely, (ρl,ul) and (ρr,ur), respectively, withρl,r ≥0 as well asul,rbeing constants.

    For the pressureless flow, we assume that the particles on the left side will stick to the piston once they collide.Meanwhile, the polytropic gas satisfies the slip boundary condition on the right surface of the piston.A special interesting case is where there is no gas on the left side of the piston; things are then reduced to a free piston problem with the free piston taking a positive mass and velocity moving between the polytropic gas and the vacuum.

    Both of the equations characterizing the pressureless flow on the left hand side of the piston and the polytropic gas on the right hand side of it are compressible isentropic Euler systems given by (1.1), where the state equation of the former isp=0 and that of the latter isp=ργ.Assume that the trajectory of the piston isx=y(t) withy(0) = 0.Then we have a typical problem of fluid-solid coupling.For the pressureless flow on the left of the piston, we have the free boundary problem

    Note that the term on the right hand side of the second equation in (5.21) represents the force acting on the particles sticking on the piston.This forcewp(t) is due to the reacting forces resulting from the change in the momentum of those particles hitting and then sticking to the piston, according to Newton’s Third Law.Obviously, this only acts on the particles that stick to the surface of the piston.Meanwhile,ωρ(t) represents the total mass of particles that stick to the piston from 0 to timet.

    For the polytropic flow on the right hand side of the piston, the free boundary problem is formulated as

    The fourth line is an impermeability condition established on the free boundary.

    Moreover, the motion of the piston satisfies Newton’s Second Law,

    where the subscript “+” denotes the right limit state near the piston, and “-” denotes the left one in the sequel.

    Determining the solution of the piston problem is reduced to solving the coupled problems(5.21)-(5.23).By the theorem of momentum, we can get the force enacted by the piston from(5.21) as follows:

    Comparing (5.25) and (5.27), the same equation for the trajectory of a piston and a delta shock front means that solving the piston problem (5.21)-(5.23) is equivalent to solving the singular Riemann problem, just by regarding the delta shock as a free piston.

    We conclude the main results of this section as follows:

    Theorem 5.8 The singular Riemann problem of one-dimensional isentropic compressible Euler equations with discontinuous flux(1.1),(1.2),(1.6)is equivalent to the free piston problem characterized by(5.21),(5.22)and(5.23)in the sense that the weight of the delta shockωρ(t)is exactly the mass of the piston at timet, and the function of the delta shockx=x(t)is nothing but the trajectory of the piston, provided that the same boundary conditions are proposed on the delta shock and the piston.

    Then it follows from Theorem 5.6 that we have the following conclusion about the existence of the piston problem:

    Theorem 5.9 The free piston problem (5.21)-(5.23) admits a unique local piecewise smooth solution ifu0<~ur.More precisely, this solution consists of a forward rarefaction wave in front of the piston whenu0≤ur, and a shock wave when ~ur >u0>ur.

    Furthermore,by the analysis regarding the singular Riemann problem,we also can conclude some properties of the movement of the piston for different initial data.For example, iful ≤u0≤ur,there is a vacuum state on the left side of the piston.If,in addition,u0>the piston moves with a reduced acceleration until the density of the gas on its right hand side decreases to zero.

    Remark 5.10 For the singular solutions constructed in Section 4, similar results can be derived so long as we consider the properties of the piston, such as whether it can release or absorb gases on either side.More precisely, we have the following two cases:

    Case 1 The left surface of the piston can release, if the gas leaves it, or absorb, if the gas runs to it, pressureless gas that keeps the flow near it uniform, while the right surface can release, if the gas leaves it, or absorb, if the gas runs to it, a polytropic gas that keeps the flow near it uniform.

    Case 2 The piston absorbs all of the particles that run to it, but it never releases any gas on either side.This corresponds to the over-compressing condition on the piston.

    For any given timet >0, during the time [0,t], we assume that the mass concentrated on the left and right surface of the piston are, respectively,

    The positiveml,r(t) means the absorption of the piston from the left or right side, and the negative one means the release.By the theorem of momentum,

    where we have chosen thex-axis direction as the positive direction of velocity.This, together with (5.28), gives that

    which is equivalent to (2.14) as well.

    6 Discussion

    So far, we have finished the construction of solutions to the singular Riemann problem for isentropic compressible Euler equations with discontinuous flux (1.1), (1.2), (1.6), which consists of a delta shock, a rarefaction wave, a shock wave, a contact discontinuity and a vacuum.Note that, whenρ0= 0, we refer to (1.6) as the classical Riemann initial data, but there are still differences here.The classical Riemann problem corresponds to the shock tube experiment, where the interface between two initial states will disappear.In this paper, this interface will always exist, and it has some special properties, such as releasing or absorbing gases on either side.Under the over-compressing condition (1.5), the solution may still be nonunique, due to the multiplicity of intermediate states between waves.It is necessary to propose a suitable condition on the discontinuity of the flux to get the unique solution.Fortunately,the combination of a slip boundary condition for a polytropic gas and an adhesion condition for a pressureless gas works well, and we can get a unique local piecewise smooth solution of a free piston problem under this condition by showing the equivalence of the singular Riemann problem and the piston problem.This indicates that by proposing different suitable conditions on the discontinuity of the flux, one may get solutions corresponding to different physical problems.

    As was shown in Section 5 above, we obtained the local existence of piecewise smooth solutions with clear geometrical structure to a free piston problem, under the sufficient but not necessary condition thatu0<~ur.Whenu0≥~ur,some other theories or methods are necessary to obtain corresponding solutions that might contain more abundant wave structures.We will,in future,conduct work regarding the global existence of the solutions we constructed in Section 5.

    Conflict of InterestThe authors declare no conflict of interest.

    猜你喜歡
    海棠
    海棠春
    讀《明天要遠(yuǎn)足》,寫“我要去遠(yuǎn)足”
    《大嘴狗》
    找春天①
    我喜歡的季節(jié)
    我有一個(gè)特別的月餅
    海棠林
    藝術(shù)百家:海棠
    戊戌海棠雅集
    中華詩詞(2018年7期)2018-11-08 12:47:28
    知音少,便做一朵未開海棠
    亚洲精品久久午夜乱码| 久久久精品94久久精品| 欧美老熟妇乱子伦牲交| 美女xxoo啪啪120秒动态图| 亚洲伊人久久精品综合| 大香蕉久久成人网| 波野结衣二区三区在线| 高清午夜精品一区二区三区| 日本与韩国留学比较| 亚洲综合色惰| 精品久久久久久电影网| 一区在线观看完整版| 亚洲欧美色中文字幕在线| 亚洲国产色片| 国产亚洲午夜精品一区二区久久| 汤姆久久久久久久影院中文字幕| av又黄又爽大尺度在线免费看| 日韩成人av中文字幕在线观看| 亚洲av免费高清在线观看| 老司机亚洲免费影院| 国产日韩一区二区三区精品不卡 | 交换朋友夫妻互换小说| 尾随美女入室| 午夜福利,免费看| 亚洲三级黄色毛片| 精品国产露脸久久av麻豆| 国产极品粉嫩免费观看在线 | 国产淫语在线视频| 久久热精品热| 久久人人爽人人爽人人片va| 亚洲精品久久午夜乱码| 在线观看免费高清a一片| 色5月婷婷丁香| 高清毛片免费看| 久久久久久久国产电影| 免费黄网站久久成人精品| 免费少妇av软件| 国产精品久久久久久久电影| 亚洲欧美成人综合另类久久久| 18禁裸乳无遮挡动漫免费视频| 亚洲一级一片aⅴ在线观看| 777米奇影视久久| 日产精品乱码卡一卡2卡三| 国产成人精品在线电影| 不卡视频在线观看欧美| 成年人免费黄色播放视频| 女性生殖器流出的白浆| 亚洲在久久综合| 两个人的视频大全免费| 久久久久精品久久久久真实原创| 亚洲av.av天堂| 又粗又硬又长又爽又黄的视频| 国产精品偷伦视频观看了| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区三区四区免费观看| 波野结衣二区三区在线| 美女脱内裤让男人舔精品视频| 日本-黄色视频高清免费观看| 七月丁香在线播放| av女优亚洲男人天堂| 国产色婷婷99| 久久99一区二区三区| 亚洲美女视频黄频| 免费av中文字幕在线| 国产日韩一区二区三区精品不卡 | 18在线观看网站| 亚洲欧美一区二区三区黑人 | 欧美日韩视频高清一区二区三区二| 精品一区二区三区视频在线| 在线播放无遮挡| 免费人妻精品一区二区三区视频| 人人澡人人妻人| 日本黄色片子视频| 午夜免费观看性视频| 国产白丝娇喘喷水9色精品| 全区人妻精品视频| 久久亚洲国产成人精品v| 男人添女人高潮全过程视频| 男人添女人高潮全过程视频| av视频免费观看在线观看| 免费大片黄手机在线观看| 欧美xxⅹ黑人| 国国产精品蜜臀av免费| 99热6这里只有精品| 久久毛片免费看一区二区三区| av在线观看视频网站免费| 国产精品一区二区在线不卡| 天堂俺去俺来也www色官网| 欧美+日韩+精品| 亚洲天堂av无毛| 久久影院123| 欧美另类一区| 男女边摸边吃奶| 婷婷色综合大香蕉| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久av不卡| 国产精品一二三区在线看| 成年女人在线观看亚洲视频| 少妇高潮的动态图| 人成视频在线观看免费观看| 久久免费观看电影| av国产精品久久久久影院| 国产深夜福利视频在线观看| 亚洲精品aⅴ在线观看| 国产精品免费大片| 天堂俺去俺来也www色官网| 国产在线视频一区二区| 黄片无遮挡物在线观看| 国产精品国产三级专区第一集| 免费少妇av软件| 日韩视频在线欧美| 国产一区有黄有色的免费视频| 国产淫语在线视频| 国产又色又爽无遮挡免| 亚洲经典国产精华液单| 欧美xxⅹ黑人| 国产极品天堂在线| 日韩制服骚丝袜av| 国产高清有码在线观看视频| 另类精品久久| 日韩一区二区视频免费看| 日本黄色日本黄色录像| 久久精品熟女亚洲av麻豆精品| 亚洲av成人精品一区久久| 大又大粗又爽又黄少妇毛片口| 精品人妻熟女毛片av久久网站| 热99久久久久精品小说推荐| 精品酒店卫生间| 桃花免费在线播放| 99久国产av精品国产电影| 狠狠精品人妻久久久久久综合| 夜夜爽夜夜爽视频| 美女中出高潮动态图| 国产 精品1| 国产精品嫩草影院av在线观看| 久久人人爽人人片av| 美女cb高潮喷水在线观看| 亚洲国产精品专区欧美| 欧美另类一区| 国产乱来视频区| 蜜桃在线观看..| 国产免费一级a男人的天堂| 亚洲美女黄色视频免费看| 亚洲色图 男人天堂 中文字幕 | 国产日韩欧美视频二区| 亚洲美女搞黄在线观看| 午夜视频国产福利| 亚洲性久久影院| 高清午夜精品一区二区三区| 国产欧美日韩综合在线一区二区| 制服人妻中文乱码| 国产高清有码在线观看视频| 一区二区三区乱码不卡18| 一本一本综合久久| 精品一区二区三区视频在线| 精品99又大又爽又粗少妇毛片| 夜夜爽夜夜爽视频| 国产男女超爽视频在线观看| 中文欧美无线码| 夫妻午夜视频| 热re99久久国产66热| 亚洲精品av麻豆狂野| 美女福利国产在线| 国产精品久久久久成人av| 天堂俺去俺来也www色官网| 十八禁高潮呻吟视频| 一个人看视频在线观看www免费| 只有这里有精品99| 亚洲综合色网址| 亚洲国产最新在线播放| 亚洲精品成人av观看孕妇| 特大巨黑吊av在线直播| 少妇精品久久久久久久| 久热这里只有精品99| 亚洲成色77777| 日韩av不卡免费在线播放| 菩萨蛮人人尽说江南好唐韦庄| av卡一久久| 久久精品国产a三级三级三级| 国产成人精品福利久久| freevideosex欧美| 男人操女人黄网站| 精品国产一区二区久久| 日韩大片免费观看网站| 国产视频首页在线观看| 午夜福利影视在线免费观看| 久久久精品区二区三区| 久久精品国产亚洲av涩爱| 国产精品秋霞免费鲁丝片| 色婷婷久久久亚洲欧美| 久久女婷五月综合色啪小说| 少妇人妻久久综合中文| 国产成人a∨麻豆精品| 欧美精品亚洲一区二区| 狠狠精品人妻久久久久久综合| 亚洲国产毛片av蜜桃av| 欧美人与善性xxx| a级毛片在线看网站| 人人妻人人澡人人看| 中国国产av一级| 欧美日韩av久久| 日韩电影二区| 欧美日韩一区二区视频在线观看视频在线| 肉色欧美久久久久久久蜜桃| 欧美日韩成人在线一区二区| 日本色播在线视频| 男女高潮啪啪啪动态图| 亚洲综合精品二区| av线在线观看网站| 亚洲av二区三区四区| av免费在线看不卡| 欧美日韩国产mv在线观看视频| 如何舔出高潮| 免费黄频网站在线观看国产| 久久婷婷青草| 国产精品不卡视频一区二区| 中文字幕人妻丝袜制服| 免费播放大片免费观看视频在线观看| 观看美女的网站| 亚洲av综合色区一区| 国产男女超爽视频在线观看| 男女国产视频网站| 久久国产精品男人的天堂亚洲 | 国产熟女欧美一区二区| 国产午夜精品久久久久久一区二区三区| 国产精品一国产av| 草草在线视频免费看| 精品久久久久久久久av| 一本久久精品| 黑人巨大精品欧美一区二区蜜桃 | 伦精品一区二区三区| 中文欧美无线码| 成人18禁高潮啪啪吃奶动态图 | av女优亚洲男人天堂| 亚洲av综合色区一区| 一区二区日韩欧美中文字幕 | 亚洲精品一区蜜桃| 国产av码专区亚洲av| 久久久久久久久大av| 日韩欧美精品免费久久| 国产69精品久久久久777片| 亚洲av在线观看美女高潮| 考比视频在线观看| 伊人亚洲综合成人网| 免费不卡的大黄色大毛片视频在线观看| 成人无遮挡网站| 如何舔出高潮| 99国产精品免费福利视频| 久久精品久久久久久噜噜老黄| 久久久久精品性色| 大又大粗又爽又黄少妇毛片口| 亚洲精品,欧美精品| 亚洲欧美日韩卡通动漫| 特大巨黑吊av在线直播| 伊人久久精品亚洲午夜| 久久久午夜欧美精品| 99热网站在线观看| 国国产精品蜜臀av免费| 亚洲婷婷狠狠爱综合网| 久久久久久久国产电影| 亚洲久久久国产精品| 汤姆久久久久久久影院中文字幕| a级毛片免费高清观看在线播放| 热99国产精品久久久久久7| 久久精品国产亚洲av天美| 亚洲第一av免费看| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 精品人妻在线不人妻| 国产成人免费观看mmmm| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 夜夜骑夜夜射夜夜干| 男人添女人高潮全过程视频| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区三区| 欧美 亚洲 国产 日韩一| 如何舔出高潮| 日韩在线高清观看一区二区三区| 人体艺术视频欧美日本| 女人精品久久久久毛片| av国产久精品久网站免费入址| 免费黄网站久久成人精品| 天堂中文最新版在线下载| 久久ye,这里只有精品| 久热久热在线精品观看| 国产在线免费精品| 欧美bdsm另类| 中文字幕最新亚洲高清| 成年人午夜在线观看视频| 久久久久久久亚洲中文字幕| 亚洲第一区二区三区不卡| 亚洲怡红院男人天堂| 性色avwww在线观看| 夜夜爽夜夜爽视频| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说| a级片在线免费高清观看视频| 日韩av在线免费看完整版不卡| 亚洲精品乱久久久久久| 亚洲精品久久午夜乱码| 国产黄色视频一区二区在线观看| 一个人免费看片子| 老司机影院毛片| 国产成人免费观看mmmm| 日韩亚洲欧美综合| 少妇被粗大的猛进出69影院 | 亚洲欧美清纯卡通| 亚洲少妇的诱惑av| 黑丝袜美女国产一区| 久久久久视频综合| 婷婷色综合www| 精品一区二区三卡| 国产黄频视频在线观看| 国产精品99久久99久久久不卡 | av不卡在线播放| 久久久国产精品麻豆| 大码成人一级视频| 国产 精品1| 日韩在线高清观看一区二区三区| 亚洲国产精品一区二区三区在线| 久久久久久久久大av| 国产在视频线精品| 新久久久久国产一级毛片| 欧美日韩一区二区视频在线观看视频在线| 日本猛色少妇xxxxx猛交久久| a级毛片黄视频| 精品久久久噜噜| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| 大码成人一级视频| 高清在线视频一区二区三区| 99热国产这里只有精品6| a级毛色黄片| 人成视频在线观看免费观看| 亚洲国产精品成人久久小说| 久久热精品热| 国产精品99久久99久久久不卡 | 亚洲精品视频女| 有码 亚洲区| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式| 国产精品.久久久| 精品一区二区免费观看| 色婷婷久久久亚洲欧美| h视频一区二区三区| 国产女主播在线喷水免费视频网站| 制服人妻中文乱码| 久久毛片免费看一区二区三区| 一级a做视频免费观看| 国产一区有黄有色的免费视频| 九九久久精品国产亚洲av麻豆| 亚洲国产成人一精品久久久| 汤姆久久久久久久影院中文字幕| 在线观看免费高清a一片| 日本91视频免费播放| 免费高清在线观看日韩| 国产黄频视频在线观看| 国产亚洲最大av| 曰老女人黄片| 久久精品熟女亚洲av麻豆精品| 在线观看免费日韩欧美大片 | 午夜免费观看性视频| 亚洲精品色激情综合| 少妇精品久久久久久久| 在线观看三级黄色| 久久狼人影院| 国产亚洲午夜精品一区二区久久| 国产又色又爽无遮挡免| 丰满迷人的少妇在线观看| 久久狼人影院| 色婷婷av一区二区三区视频| 九色亚洲精品在线播放| 亚洲中文av在线| 中国三级夫妇交换| 99久久精品国产国产毛片| 18禁动态无遮挡网站| 美女脱内裤让男人舔精品视频| 777米奇影视久久| 国产精品国产三级国产专区5o| 国产伦理片在线播放av一区| 99热网站在线观看| 成人二区视频| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 高清黄色对白视频在线免费看| 高清欧美精品videossex| 国产有黄有色有爽视频| 美女cb高潮喷水在线观看| 天堂中文最新版在线下载| 美女福利国产在线| 日本欧美视频一区| 中文字幕人妻丝袜制服| 女人精品久久久久毛片| 亚洲av免费高清在线观看| 成年女人在线观看亚洲视频| 中文字幕人妻丝袜制服| 91午夜精品亚洲一区二区三区| 国产亚洲欧美精品永久| 性色avwww在线观看| 如何舔出高潮| 九九在线视频观看精品| 男女高潮啪啪啪动态图| 免费黄网站久久成人精品| 久久久久久久久久成人| 亚洲av男天堂| 丰满饥渴人妻一区二区三| 汤姆久久久久久久影院中文字幕| 午夜激情福利司机影院| 亚洲美女黄色视频免费看| 天堂俺去俺来也www色官网| 日本与韩国留学比较| 人妻系列 视频| 国产黄频视频在线观看| 亚洲av不卡在线观看| 99热国产这里只有精品6| 久久久久久人妻| 九九久久精品国产亚洲av麻豆| 91成人精品电影| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 国产精品国产三级国产av玫瑰| 成人综合一区亚洲| 男女国产视频网站| 91精品三级在线观看| 国产成人av激情在线播放 | 日本vs欧美在线观看视频| 久久久国产欧美日韩av| 男女边吃奶边做爰视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 22中文网久久字幕| 久久久精品区二区三区| 如日韩欧美国产精品一区二区三区 | 亚洲少妇的诱惑av| 爱豆传媒免费全集在线观看| 人妻 亚洲 视频| 国产精品欧美亚洲77777| 日韩在线高清观看一区二区三区| 国产av国产精品国产| .国产精品久久| 欧美 亚洲 国产 日韩一| 国产女主播在线喷水免费视频网站| 少妇 在线观看| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 少妇人妻久久综合中文| 男人添女人高潮全过程视频| 久久久久人妻精品一区果冻| 久久这里有精品视频免费| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 十八禁高潮呻吟视频| av在线观看视频网站免费| 国产精品嫩草影院av在线观看| 人人妻人人澡人人爽人人夜夜| 中国国产av一级| 日韩电影二区| 国产一区二区在线观看av| 国产亚洲一区二区精品| 午夜福利网站1000一区二区三区| 免费黄色在线免费观看| 欧美bdsm另类| av在线老鸭窝| 激情五月婷婷亚洲| 国产精品一区二区三区四区免费观看| 春色校园在线视频观看| 午夜日本视频在线| 99热这里只有精品一区| xxx大片免费视频| av视频免费观看在线观看| 亚洲精品久久午夜乱码| 99精国产麻豆久久婷婷| a 毛片基地| 亚洲精品日韩在线中文字幕| 久热久热在线精品观看| 一本色道久久久久久精品综合| 中文字幕久久专区| 中国国产av一级| 性高湖久久久久久久久免费观看| 丝袜美足系列| 久久久久久久久久久免费av| 日本av手机在线免费观看| 成人午夜精彩视频在线观看| 久久精品久久久久久噜噜老黄| 狂野欧美白嫩少妇大欣赏| 亚洲人成77777在线视频| 国产综合精华液| 久久这里有精品视频免费| 久久久久精品久久久久真实原创| 免费黄色在线免费观看| 国产高清不卡午夜福利| 日本爱情动作片www.在线观看| 另类亚洲欧美激情| 亚洲综合色网址| 一级毛片aaaaaa免费看小| 青青草视频在线视频观看| 亚洲精华国产精华液的使用体验| 国产精品欧美亚洲77777| 亚洲av成人精品一区久久| 69精品国产乱码久久久| 国内精品宾馆在线| 亚洲,欧美,日韩| 国产亚洲午夜精品一区二区久久| 欧美日韩亚洲高清精品| 91午夜精品亚洲一区二区三区| 在现免费观看毛片| 免费大片黄手机在线观看| 国国产精品蜜臀av免费| 国产免费又黄又爽又色| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久人妻| 精品一区二区免费观看| 少妇被粗大的猛进出69影院 | 久久99精品国语久久久| 亚洲精品第二区| 满18在线观看网站| 亚洲,一卡二卡三卡| 午夜免费观看性视频| 桃花免费在线播放| 国产爽快片一区二区三区| 日韩av不卡免费在线播放| 热99国产精品久久久久久7| 狠狠婷婷综合久久久久久88av| 亚洲成人手机| videosex国产| 国产成人免费无遮挡视频| 黄色一级大片看看| 国产片特级美女逼逼视频| 国产 一区精品| 国产成人精品一,二区| 亚洲av电影在线观看一区二区三区| 久久久久久久久久人人人人人人| 下体分泌物呈黄色| 国产一区二区在线观看日韩| 日日啪夜夜爽| 中国国产av一级| 99精国产麻豆久久婷婷| 青春草国产在线视频| av免费在线看不卡| 色5月婷婷丁香| 免费观看的影片在线观看| 成人手机av| 精品亚洲成a人片在线观看| 成人18禁高潮啪啪吃奶动态图 | 欧美精品一区二区免费开放| 欧美丝袜亚洲另类| 欧美日韩成人在线一区二区| 麻豆乱淫一区二区| 久久热精品热| 国产精品熟女久久久久浪| 欧美三级亚洲精品| 午夜福利影视在线免费观看| 女性生殖器流出的白浆| 免费观看性生交大片5| 欧美+日韩+精品| 天美传媒精品一区二区| 91精品三级在线观看| 亚洲情色 制服丝袜| 麻豆成人av视频| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放| 国产成人精品一,二区| 中国三级夫妇交换| 日日啪夜夜爽| 十分钟在线观看高清视频www| 国产在线一区二区三区精| av福利片在线| 18禁观看日本| 少妇被粗大猛烈的视频| 国产精品蜜桃在线观看| 国产 一区精品| 91精品国产九色| 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 水蜜桃什么品种好| 亚洲精品亚洲一区二区| 国产视频内射| 日本色播在线视频| 免费黄网站久久成人精品| 精品一区二区免费观看| 亚洲少妇的诱惑av| 久久久国产一区二区| 女性生殖器流出的白浆| 99热6这里只有精品| 国产一区二区三区av在线| 久久久久久久久久久丰满| www.色视频.com| 丝袜喷水一区| 亚洲综合色网址| 男女无遮挡免费网站观看| 亚洲内射少妇av| h视频一区二区三区| 久久久精品94久久精品| h视频一区二区三区| 午夜久久久在线观看| 久热这里只有精品99| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 精品少妇内射三级| 亚洲成人av在线免费| 高清欧美精品videossex| 美女大奶头黄色视频| 韩国高清视频一区二区三区| 99热国产这里只有精品6| 国产午夜精品一二区理论片| 一本大道久久a久久精品| 美女主播在线视频| 午夜免费男女啪啪视频观看| 最近手机中文字幕大全| 少妇被粗大猛烈的视频| 热99久久久久精品小说推荐| 亚洲av免费高清在线观看| 视频区图区小说| 亚洲性久久影院|