arxiv: v1 [math.ap] 20 Aug 2014

Size: px
Start display at page:

Download "arxiv: v1 [math.ap] 20 Aug 2014"

Transcription

1 Initial Layer an Relaxation Limit of Non-Isentroic Comressible Euler Equations with Daming arxiv: v1 [math.ap] 2 Aug 214 Fuzhou Wu Mathematical Sciences Center, Tsinghua University Beijing 184, China Abstract In this aer, we stuy the relaxation limit of the relaxing Cauchy roblem for non-isentroic comressible Euler equations with aming in multi-imensions. We rove that the velocity of the relaxing equations converges weakly to that of the relaxe equations, while other variables of the relaxing equations converges strongly to the corresoning variables of the relaxe equations. We show that as relaxation time aroaches, there exists an initial layer for the illreare ata, the convergence of the velocity is strong outsie the layer; while there is no initial layer for the well-reare ata, the convergence of the velocity is strong near t =. Keywors: non-isentroic Euler equation with aming, relaxation limit, ill-reare ata, initial layer Contents 1 Introuction 2 2 Preliminaries 7 3 A Priori Estimates for the Relaxe Equations 1 4 Existence in [, T] of Classical Solutions to the Relaxe Equations 23 5 Uniform A Priori Estimates for the Relaxing Equations 25 6 Initial Layer an Relaxation Limit of the Relaxing Equations 32 7 Extensions 37 References 4 michael8723@gmail.com; fuzhou.wu@yahoo.com 1

2 1 Introuction In this aer, we use,u,s, to enote the ressure, velocity, entroy an ensity of ieal gases resectively with the equation of state = (,S) := 1 γ A 1 γ ex{ S γ }, where A >, γ = C C V > 1 are constants. Assume the numbers,, S satisfy >, >, = A γe S. Then we stuy the relaxation limit of the relaxing Cauchy roblem for 3D non-isentroic comressible Euler equations with aming: t +u +γ u =, τ 2 u t +τ 2 u u+ 1 +u =, (1.1) S t +u S =, (,u,s)(x,) = ( (x,τ), u(x,τ) τ,s (x,τ)), where ( (x,τ),u (x,τ),s (x,τ)) are small erturbations of (,, S). In (1.1), (,u,s, ) (,, S, ) as x +. τ is a small ositive arameter reresenting the relaxation time, let τ (,1]. The ensity satisfies the equation t + u + u = by (1.1), an (x,τ) = ( (x,τ),s (x,τ)) is small erturbation of = (, S). The equations (1.1) are erive from ˆ t +û ˆ+γˆ û =, û t +û û+ 1ˆ ˆ+ τû 1 =, Ŝ t +û Ŝ =, with the rescaling of variables: (ˆ,û,Ŝ)(x,) = ( (x,τ),u (x,τ),s (x,τ)), (1.2) t = τt,u(x,t) = û τ (x,t ),(x,t) = ˆ(x,t ), (x,t) = ˆ (x,t ),S(x,t) = Ŝ(x,t ), (1.3) then (,u,s, ) satisfy the equations (1.1). Let τ in the relaxing equations (1.1), we formally obtain the following relaxe equations t +u +γ u =, 1 +u =, (1.4) S t +u S =, where = (,S). (,S)(x,) = (lim τ (x,τ), lim τ S (x,τ)), Due to its funamental imortance in both alication an nonlinear PDE theory, the relaxation limit roblems have been attracting much attention. We survey there some results closely relate to this aer. 2

3 For the relaxing isothermal comressible Euler equations with aming: { t + ( u) =, u t + u u+ σ τ u =, (1.5) where σ 2 = Rθ is constant. In the Sobolev sace H s ((R) ),s Z,s > 1+ 2 (see [3]), u τ σ2 in D (R R + ), converges strongly to the solution of the heat equation in C([,T],H s (B r )), where < s < s, B r is a ball with raius r. The results of [3] was extene in [24] to more general Sobolev sace of fractional orer. In [8], (1.5) was stuie in one imension with BV large ata away from vacuum, an it was rove in [8] that converges strongly to the solution of the heat equation in L 2 (R [,T]) (global in sace) by using the stream function. For the relaxing isentroic comressible Euler equations with aming { t + ( u) =, u t + u u+ + 1 τ u =, (1.6) where (ρ) = Aρ γ. In the Sobolev sace H s ((R) ),s Z,s > (see [13]), u τ in D (R R + ), converges strongly to the solution of the orous meia equation in C([,T],H s (B r )), where < s < s. In the Besov sace B2,1 σ (R ),σ = 1+ 2 (see [26]) an in the Chemin-Lerner sace (see [25]), the ensity of (1.6) converges strongly to the solution of the orous meia equation. Relaxation limit roblem also aears in Euler-Poisson equations, see [15, 9, 1, 7] for weak solutions an [2, 23] for smooth solutions. It has been rove that the current ensity, which is the rouct of the electron ensity an electron velocity, converges weakly to that of the rift-iffusion moel. If the initial ata are well-reare, the current ensity converges strongly to that of the rift-iffusion moel (see [28]). If the initial ata are ill-reare, the authors (see [16]) rove the ifference between the current ensity of 1D hyroynamic moel an that of the rift-iffusion moel ecays exonentially fast in the large time interval [, 1 β log( 1 τ )] with λ (,1),β >. The key of the roof in [16] is λ that the solutions of the relaxing an relaxe equations converge to the corresoning stationary solutions exonentially fast while both stationary solutions are close to each other. As to the relaxation limit of weak solutions (see [4]) an classical solutions (see [1, 27, 11]) to non-isentroic Euler-Poisson equations, the current ensity converges weakly to that of the energy-transortation moel or rift-iffusion moel. However, there have been no rigorous analysis of the initial layer an strong convergence of the velocity for the ill-reare ata in the above mentione aers. A main istinction of results in this aer is that we give results on the initial layer an strong convergence of the velocity. Our main concern is the non-isentroic flow (1.1), but our results are vali for the isentroic flow (1.6) an isothermal flow (1.5) (assuming no vacuum). We show that for the ill-reare initial ata, there exists an initial layer [,t ] for the velocity, where < t < Cτ 2 δ,δ >. Outsie the initial layer, the velocity of the relaxing equations converge strongly to that of the relaxe equations. Only for the wellreare initial ata, there is no initial layer, the convergence of the velocity 3

4 is strong near t =. The key of our analysis in this aer is uniform a riori estimates with resect to τ an ointwise ecay of the quantity u+ 1. Also, the methos in this aer can be alie to the relaxation limit roblems for Euler-Poisson equations an Euler-Maxwell equations. The first result in this aer is the following convergence result: Theorem 1.1. Suose that the initial ata for the relaxing Cauchy roblem (1.1) satisfy ( (x,τ),u (x,τ),s (x,τ)) H 4 ( ), inf lim (x,τ) >, x τ inf x (x,τ) >, ( (x,τ),u (x,τ),s (x,τ) S) H4 ( ) is sufficiently small for some constants > an S. Then for any finite T >, the roblem (1.1) amits a unique solution (,u,s, ) in [,T] satisfying l t( ), τ l tu, l t(s S), l t( ) L ([,T],H 4 l ( )), l 2,, u, S S, l 2 Hl ([,T],H 4 l ( )), such that as τ, in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, (1.7) S S in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, u ũ in l 2 Hl ([,T],H 4 l ( )), (1.8) for some function (, S,,ũ) which is a weak solution to the relaxe equations (1.4) with the ata (lim (x,τ), lim S (x,τ)), where K enotes any comact τ τ subset of. (, S,,ũ) is the classical solution to (1.4), if (lim (x,τ) τ, lim S (x,τ) S) H 5 ( ) H 4 ( ) is satisfie. τ The main results concerne with the initial layer an strong convergence of the velocity are state in the following Theorem: Theorem 1.2. Let (,u,s, ) an (,ũ, S, ) be the solutions obtaine in Theorem 1.1 in [,T] an K enotes any comact subset of. Then it hols that, 1 for the ill-reare ata, i.e., lim τ τ u (x,τ)+ 1 (x,τ) (x,τ), there exists an initial layer [,t ] with < t < Cτ 2 δ for the velocity u, where C >, δ >, such that as τ, u ũ in C((,T], C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6; (1.9) 1 for the well-reare ata, i.e., lim τ τ u (x,τ)+ 1 (x,τ) (x,τ) =, H 2 ( ) as τ, u ũ in C([,T],C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6. (1.1) Remark 1.3. (i) Comare (1.9) with (1.1). In (1.9), (,T] can not be relace by [,T]. That is, the convergence is not uniform near t =, since the initial layer evelos for the ill-reare ata. 4

5 (ii) For any fixe t (,T), u ũ in C([t,T], C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, (1.11) while (1.11) is equivalent to (1.9) ue to the arbitrariness of t. (iii) If one relaces the initial ata ( (x,τ),u (x,τ),s (x,τ), (x,τ)) with ( (x),u (x),s (x), (x)) which are ineenent of τ, then for the well-reare ata, (x) const an u (x) (equilibrium states); while for the ill-reare ata, (x) const or u (x) (non-equilibrium states). In the following, we give more comments on Theorem 1.1 an Theorem 1.2. If τ > is fixe, the global existence of classical solutions to the equations (1.2) is rove in [21, 22, 29]. However, for the relaxation limit roblem in this aer, τ > is variant an aroaches, so we nee the uniform existence of the solutions an uniform regularities (1.7) which are ifferent from [21, 22, 29]. The uniform a riori estimates with resect to τ rouce the convergence results. To make sure that the solutions of the relaxe equations remain classical an the estimates relate to the initial layer can be roceee, ( (x,τ),u (x,τ),s (x,τ) are require to be in H 4 ( ). The uniform a riori estimates for the equations (1.1) imly the uniform regularities (1.7). Passing to the limit, we have the convergence results (1.8). Let us give some comments an remarks on the initial layer as follows. The asymtotic exansions of the solutions to (1.1) give us some inication of the initial layer. We illustrate this as follows: if u (x,τ) = O(τ), we assume that the initial ata have asymtotic exansion ( (x,τ), u (x,τ),s (x,τ), (x,τ)) = τ m τ 2m ( m,um,sm, m ), an solutions of the equations (1.1) have the asymtotic exansion ((x,t),u(x,t),s(x,t), (x,t)) = m τ 2m ( m (x,t),u m (x,t),s m (x,t), m(x,t)), then the leaing orer rofiles satisfy the equations t +u +γ u =, u + 1 =, t S +u S =, (,S )(x,) = (,S), (1.12) if the initial velocity is well-reare, i.e., u = 1, where = (,S ). Note that in this case, u + 1 = matches u + 1 =. However, if the initial velocity is ill-reare, i. e., u 1, we may assume that the velocity has the asymtotic exansion u(x,t) = (u m (x,t)+û m (x,z)), z = t τ 2, m 5

6 then the leaing orer rofile of the initial layer correction û satisfies the equation { z û +û =, û (x,) = u + 1. (1.13) Then û(x,z) = (u + 1 )e z. Thus, the ifference between u an 1 ecays exonentially within the initial layer. The above arguments of (1.12) an (1.13) inicate the relationshi between the existence of initial layer an a class of initial ata on a formal level. We are not concerne with the asymtotic exansions in this aer (as to asymtotic exansion analysis in relaxation limit roblem for Euler-Poisson equations, see [6, 2, 12]; for Euler-Maxwell equations, see [17, 5]) an only focus on rigorous analysis of the initial layer an relaxation limit of the relaxing equations (1.1). The relaxation limit is a singular limit, since (, 1,S, ) converge to (,ũ, S, ), instea of u ũ. In orer to measure the ifference between u an 1, we introuce a quantity: where k 1 > is a constant efine in Section 2. η = 1 k 1 (u+ 1 ), (1.14) The ointwise ecay of η outsie the initial layer is the key to the strong convergence of the velocity. η satisfies the following transortation equations with aming an forcing η t +u η + 1 τ2η = forcing terms, (1.15) where τ [forcing terms] is boune uniform with resect to τ, the aming effect becomes stronger as τ ecreases. Also, η satisfies another equation: η = τ2 k 1 (u t +u u). (1.16) Then the equations (1.15) an (1.16) rouce the following estimates resectively: η 2 C η 2 H 2 ( ) C η t= 2 H 2 ( ) ex{ t τ }+Cτ 2, 2 T η t 2 H 1 ( ) x Cτ2. (1.17) Therefore, for the well-reare ata, lim τ η(t) = for any t [,T], there is no iscreancy between η(t) t> an η t=, an thus there is no initial layer, u ũ strongly in C([,T],C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6. While for the ill-reare ata, lim η t= while lim η(t) = for any τ τ t [t,t] where t = Cτ 2 δ. The iscreancy between η t= an η(t) t=t make the initial layer exist. Within the layer [,t ], η ecreases raily, 1 is uniform boune, then u changes ramatically. Outsie the layer, u converges 6

7 strongly to ũ. Namely, as τ, u ũ strongly in C((,T],C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6. Finally, the main results of this aer can be extene to erioic omains. Due to the convenience of erioic bounary conitions, a riori estimates for T 3 are similar to those for, the results for T 3 are state as follows: Suose that the initial ata for (1.1) satisfy ( (x,τ),u (x,τ),s (x,τ)) H 4 (T 3 ) an other similar conitions given in Theorem 1.1. Then for any finite T >, the roblem (1.1) amits a unique solution (,u,s, ) in [,T] such that as τ, (,S, ) (, S, ) in C([,T],C 2+µ1 (T 3 ) W 3,µ2 (T 3 )),µ 1 (, 1 2 ),2 µ 2 < 6, u ũ in l 2 Hl ([,T],H 4 l (T 3 )), for some function (, S,,ũ) which is a weak solution to (1.4) with the ata (lim (x,τ), lim S (x,τ) S). (, S,,ũ) is the classical solution to (1.4), τ τ if the ata is in H 4 (T 3 ) H 3 (T 3 ). For the ill-reare ata, there exists an initial layer for the velocity u, such that as τ, u ũ in C((,T], C +µ1 (T 3 ) W 1,µ2 (T 3 )),µ 1 (, 1 2 ),2 µ 2 < 6; for the well-reare ata, as τ, u ũ in C([,T],C +µ1 (T 3 ) W 1,µ2 (T 3 )),µ 1 (, 1 2 ),2 µ 2 < 6. The rest of this aer is organize as follows: In Section 2, we reformulate the equations into aroriate forms an erive the equations of η. In Section 3, we rove a riori estimates in [,T] for the relaxe equations. In Section 4, we rove the existence in [,T] of classical solutions to the relaxe equations. In Section 5, we rove uniform a riori estimates for the relaxing equations. In Section 6, we iscuss the relaxation limit an initial layer of the relaxing equations. In Section 7, we exten the main results of this aer to erioic omains. 2 Preliminaries In this section, we will reformulate the equations (1.1) into aroriate forms, efine the energy quantities an erive the equations of the quantity η. Fortherelaxingequations(1.1)togetherwiththeirinitialata(,u,S, ) an constants, S,, we introuce the constants: k 1 = 1 γ, k 2 = γ, efine the variables: ξ =, φ = S S, ζ =, v = 1 k 1 u. 7

8 In orer to have the solutions of the relaxing equations (1.1) in any finite time interval [, T] an obtain uniform a riori estimates, we reformulate the equations (1.1) into the following form: ξ t +k 2 v = γk 1 ξ v k 1 v ξ, τ 2 v t +k 2 ξ +v = k 1 τ 2 v v + 1 k 1 ( 1 1 ) ξ, (2.1) φ t = k 1 v φ, (ξ,v,φ)(x,) = ( (x,τ), 1 k 1τ u (x,τ),s (x,τ) S), where = ζ + = (ξ +,φ+ S) an ζ satisfies the equation ζ t +k 1 v ζ + k 1 v = by (2.1). (ξ,v,φ,ζ) (,,,) as x +. Letτ = inthe relaxingequations(2.1), weformallyobtainethefollowing relaxe equations, which are equivalent to the relaxe equations (1.4). ξ t +k 2 v = γk 1 ξ v k 1 v ξ, k 2 ξ +v = 1 k 1 ( 1 1 ) ξ, (2.2) φ t = k 1 v φ, (ξ,φ)(x,) = (lim τ (x,τ), lim τ S (x,τ) S), where = ζ + = (ξ +,φ+ S). The global existence of classical solutions to IBVP (2.2) with bounary conition v n Ω = has been rove in [21]. However, in this aer, we rove the existence in any finite time interval [,T] of classical solutions to Cauchy roblem (2.2) with small ata in H 5 ( ) H 4 ( ). Note that the regularity inex for is one orer higher than that for erioic or boune omains. In orer to rove the existence in any finite time interval [,T] of classical solutions to (2.2) via energy metho, we efine the following energy quantities: Definition 2.1. Define F[ξ](t) := D α ξ(t) 2 L 2 ( ) + α 4 F[ξ](t) := α 5 F X [ξ](t) := F 1 [ξ](t) = 1 α 4 α 4 F[v](t) = + α 4 F[φ](t) := α 2 α 4 F[ζ](t) := α 4 α 2 D α ξ(t) 2 L 2 ( ) + α 3 D α ξ(t) 2 L 2 ( ) + D α ξ(t) 2 L 2 ( ) α 2 α =4 D α t ξ(t) 2 L 2 ( ) D α v(t) 2 L 2 ( ) + α 2 D α φ(t) 2 L 2 ( ) + α 2 D α ζ(t) 2 L 2 ( ) + α 2 D α t ξ(t) 2 L 2 ( ), D α t ξ(t) 2 L 2 ( ), D α t ξ(t) 2 L 2 ( ), ξ (Dα ξ) 2 x α =2 ξ (Dα t ξ) 2 x, D α t v(t) 2 L 2 ( ), D α t φ(t) 2 L 2 ( ), D α t ζ(t) 2 L 2 ( ), F[ξ,v,φ,ζ](t) := F[ξ](t)+F[v](t)+F[φ](t)+F[ζ](t). (2.3) 8

9 Moreover, we can erive the evolution equations of v from (2.2), which is useful for roving a riori estimate for the L boun of F[v](t). Aly i to (2.2) 1 an then substitute i ξ for k 1 v i, we have v t = k 1 (1 γ)v( v) k 1 v v k1 2 ( v 2 )+ γ ( v). (2.4) To rove the uniform regularities (1.7) is equivalent to rove the following uniform regularities: l t ξ, τ l t v, l t φ, l t ζ L ([,T],H 4 l ( )), l 2, ξ, v, φ, ζ l 2 Hl ([,T],H 4 l ( )). (2.5) In orer to use the energy metho to erive uniform a riori estimates with resect to τ, we efine the following energy quantities: Definition 2.2. Define E[ξ](t) := E X [ξ](t) = l 2, l+ α 4 l 2,<l+ α 4 E 1 [ξ](t) = E[ξ](t) E[v](t) := l 2, l+ α 4 E 1 [v](t) = E[v](t)+ E[φ](t) := E[ζ](t) := l 2, l+ α 4 l 2, l+ α 4 l t Dα ξ(t) 2 L 2 ( ), l t Dα ξ 2 L 2 ( ), l 2,l+ α =4 ξ ( l t Dα ξ) 2 x, l td α v(t) 2 L 2 ( ), l 2,l+ α =4 E[ξ,v](t) := E[ξ](t)+E[v](t), ( 1) l t Dα v 2 x, l t Dα φ(t) 2 L 2 ( ), l t Dα ζ(t) 2 L 2 ( ), E[ξ,v,φ,ζ](t) := E[ξ,v](t)+E[φ](t)+E[ζ](t). (2.6) To rove the convergence results (1.8) is equivalent to rove the following convergence results: as τ, ξ ξ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, φ φ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, ζ ζ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, v ṽ in l 2 Hl ([,T],H 4 l ( )), (2.7) where K enotes any comact subset of. For the ill-reare ata, lim η t=, while lim η(t) t t = for any τ τ small t >. This iscreancy between η t= an η(t) t=t makes the initial layer for the velocity generate. 9

10 In orer to rove the strong convergence outsie the initial layer, we nee the ointwise ecay of the quantity: 1 η(x,t) = v(x,t)+ k 1 (x,t) ξ(x,t). (2.8) Differentiate (2.8) wit resect to t, we have η t = v t + 1 k 1 ξ t ζt k 1 2 ξ, then η satisfies the following transortation equation with aming an forcing: η t +k 1 v η + 1 τ η = v ( 1 ) ξ 1 ( v) ξ γ ξ v 2 γ ( v)+ 1 2(v (2.9) ζ + v) ξ, where τ [forcing terms] is boune uniform with resect to τ, ( v) is a matrix, the aming effect for η becomes stronger as τ ecreases. The equation (2.9) rouces the following estimate: η 2 C η 2 H 2 ( ) C η t= 2 H 2 ( ) ex{ t τ 2 }+Cτ 2. (2.1) Also, η satisfies another equation: which rouces the following estimate: η = τ 2 (v t +k 1 v v), (2.11) T η t 2 H 1 ( ) x Cτ2. (2.12) Base on the estimates (2.1) an (2.12), we have the strong convergence of v outsie the initial layer. In, we have Gagliao-Nirenberg tye inequalities, which are useful in estimating the nonlinear terms. { L 6 ( ) H 1 ( ), U 2 L 4 ( ) U (2.13) 2 U L2 ( ), (see [18]), where U is a vector or scalar function, = ( t, ). In the sequent sections, we will use the following notations: X Y enotes the estimate X CY for some imlie constant C > which may ifferent line by line. ( ) k enotes a vector in, for instance, ω k = δ ijk i v j, where δ ijk is totally anti-symmetric tensor such that δ 123 = δ 231 = δ 312 = 1,δ 213 = δ 321 = δ 132 = 1, others are. R.H.S. is the abbreviation for right han sie. 3 A Priori Estimates for the Relaxe Equations Inthissection, weerivearioriestimatesin[,t]fortherelaxeequations (2.2), where T > is finite an ineenent of τ. The following lemma concerns the estimates relate to the ensity, the roof is omitte for its simlicity. 1

11 Lemma 3.1. For any given T (,+ ), if where < ε 1, then su F[ξ,v,φ,ζ](t) ε, t T su ζ t ε, su ζ ε, su ζ ε. (3.1) t T t T t T The following lemma states that F[ξ](t) an F 1 [ξ](t) are equivalent. Lemma 3.2. For any given T (,+ ), there exists ε 1 > which is ineenent of (ξ,v,φ ), such that if su F[ξ,v,φ,ζ](t) ε 1, then ξ 3, ζ t T 2 an there exist c 1 >,c 2 > such that c 1 F[ξ](t) F 1 [ξ](t) c 2 F[ξ](t). (3.2) To make calculations simler, we calculate t F 1[ξ](t) an tf[v](t) searately. The following lemma concerns the estimate for L boun of F 1 [ξ](t). Lemma 3.3. For any given T (,+ ), if where < ε 1, then for t [,T], su F[ξ,v,φ,ζ](t) ε, t T t F 1[ξ](t)+2F[v](t) C ε(f X [ξ](t)+f[v](t)). (3.3) Proof. Let (2.2) (ξ,v), we get ( ξ 2 ) t +2k 2 ξ v+2k 2 v ξ +2 v 2 = 2γk 1 ξ 2 v 2k 1 ξv ξ + 2 k 1 ( 1 1 ) ξ v. (3.4) Integrate (3.4) in an note that (ξv)x =, we get t ξ 2 x+2 v 2 x = 2γk 1 v (ξ 2 ) 2k 1 ξv ξ + 2 k 1 ( 1 1 ) ξ vx ε ξ L 2 ( ) v L 2 ( ) (3.5) ε(f X [ξ](t)+f[v](t)). Aly D α to (2.1), where 1 α 4, we get { (D α ξ) t +k 2 (D α v) = γk 1 D α (ξ v) k 1 D α (v ξ), k 2 (D α ξ)+d α v = 1 k 1 D α [( 1 1 ) ξ]. (3.6) 11

12 Let (3.6) (D α ξ,d α v), we get ( D α ξ 2 ) t +2k 2 D α ξ (D α v)+2k 2 D α v (D α ξ)+2 D α v 2 = 2γk 1 (D α ξ)d α (ξ v) 2k 1 (D α ξ)d α (v ξ) + 2 k 1 (D α v) D α [( 1 1 ) ξ]. (3.7) Integrate (3.7) in an note that (D α ξd α v)x =, we get t D α ξ 2 x+2 D α v 2 x = 2γk 1 (D α ξ)d α (ξ v) 2k 1 (D α ξ)d α (v ξ) (3.8) + 2 k 1 (D α v) D α [( 1 1 ) ξ]x := I 1. When 1 α 3, it is easy to check that I 1 ε(f X [ξ](t)+f[v](t)). When α = 4, we estimate the quantity I 1 ξ t (Dα ξ) 2 x, then I 1 ξ t (Dα ξ) 2 x = 2γk 1 (D α ξ)ξ (D α v)x 2k 1 (D α ξ)v (D α ξ)x + 2 k 1 ( 1 1 )(Dα v) (D α ξ)x 2 ξ (Dα ξ)(d α ξ t )x t ( ξ )(Dα ξ) 2 x ε(f X [ξ](t)+f[v](t))+k 1 D α ξ 2 vx 2γk 1 ξ(d α ξ) (D α v)x+ 2 k 1 ( 1 1 )(Dα v) (D α ξ)x 2 ξ (Dα ξ)(d α ξ t )x 2 ξ (Dα ξ)[d α ξ t +k 1 γ (D α v)]x + 2 k 1 ( 1 1 )(Dα v) (D α ξ)x+ ε(f X [ξ](t)+f[v](t)). Aly D α to (2.1) 1, where α = 4, we get D α ξ t +k 1 γ (D α v) = k 1 D α (v ξ) k 1 γ Plug (3.1) into the following integral, we get ξ (Dα ξ)[d α ξ t +k 1 γ (D α v)]x α 1 > = ξ (Dα ξ)[r.h.s. of (3.1)]x ε(f X [ξ](t)+f[v](t))+ k1 2 D α ξ 2 ( ξ v)x ε(f X [ξ](t)+f[v](t)). D α1 ξ (D α2 v). (3.9) (3.1) (3.11) 12

13 Aly D α to k 1 v + ξ =, where α = 4, we get (D α ξ) = k 1 D α v l 1+ α 1 > k 1 l1 t D α1 ζ l2 t D α2 v. (3.12) Plug (3.12) into the following integral, we get ( 1 1 )(Dα v) (D α ξ)x = ( 1 1 )(Dα v) [R.H.S. of (3.12)]x εf[v](t). (3.13) t Plug (3.11) an (3.13) into (3.9), we get I 1 t ξ (Dα ξ) 2 x ε(f X [ξ](t)+f[v](t)). (3.14) After summing α, where α 4, we have ( D α ξ 2 x ) ξ (Dα ξ) 2 x +2 < α 4 α =4 ε(f X [ξ](t)+f[v](t)). < α 4 D α v 2 x (3.15) Next, we estimate the art of F 1 [ξ](t) which contains the time erivatives. Aly D α t to (2.1), where α 2, we get { (D α ξ t ) t +k 2 (D α v t ) = γk 1 D α t (ξ v) k 1 D α t (v ξ), k 2 (D α ξ t )+D α v t = 1 k 1 D α t [( 1 1 ) ξ]. (3.16) Let (3.16) (D α ξ t,d α v t ), we get ( D α ξ t 2 ) t +2k 2 D α ξ t (D α v t )+2k 2 D α v t (D α ξ t )+2 D α v t 2 = 2γk 1 (D α ξ t )D α t (ξ v) 2k 1 (D α ξ t )D α t (v ξ) + 2 k 1 (D α v t ) D α t [( 1 1 ) ξ]. (3.17) Integrate (3.17) in an note that (D α ξ t D α v t )x =, we get t D α ξ t 2 x+2 D α v t 2 x = 2γk 1 (D α ξ t )D α t (ξ v) 2k 1 (D α ξ t )D α t (v ξ) (3.18) + 2 k 1 (D α v t ) D α t [( 1 1 ) ξ]x := I 2. When α 1, it is easy to check that I 2 ε(f X [ξ](t)+f[v](t)). 13

14 When α = 2, we estimate the quantity I 2 t ξ (Dα ξ t ) 2 x, then I 2 ξ t (Dα ξ t ) 2 x = 2γk 1 (D α ξ t )ξ (D α v t )x 2k 1 (D α ξ t )v (D α ξ t )x + 2 k 1 ( 1 1 )(Dα v t ) (D α ξ t )x 2 ξ (Dα ξ t )(D α ξ tt )x t ( ξ )(Dα ξ t ) 2 x ε(f X [ξ](t)+f[v](t))+k 1 D α ξ t 2 vx 2γk 1 ξ(d α ξ t ) (D α v t )x+ 2 k 1 ( 1 1 )(Dα v t ) (D α ξ t )x 2 ξ (Dα ξ t )(D α ξ tt )x 2 ξ (Dα ξ t )[D α ξ tt +k 1 γ (D α v t )]x + 2 k 1 ( 1 1 )(Dα v t ) (D α ξ t )x+ ε(f X [ξ](t)+f[v](t)). Aly D α t to (2.1) 1, where α = 2, we get (3.19) D α ξ tt +k 1 γ (D α v t ) = k 1 D α (v t ξ) k 1 D α (v ξ t ) k 1 γ D α1 ξ t (D α2 v) k 1 γ α 1 α 1 > Plug (3.2) into the following integral, we get ξ (Dα ξ t )[D α ξ tt +k 1 γ (D α v t )]x D α1 ξ (D α2 v t ). = ξ (Dα ξ t )[R.H.S. of (3.2)]x ε(f X [ξ](t)+f[v](t))+ k1 2 D α ξ t 2 ( ξ v)x ε(f X [ξ](t)+f[v](t)). Aly D α t to k 1 v + ξ =, where α = 2, we get (D α ξ t ) = k 1 D α v t k 1 D α1 ζd α2 v t k 1 D α1 ζ t D α2 v. α 1 > α 1 Plug (3.22) into the following integral, we get ( 1 1 )(Dα v t ) (D α ξ t )x = ( 1 1 )(Dα v t ) [R.H.S. of (3.22)]x εf[v](t). (3.2) (3.21) (3.22) (3.23) Plug (3.21) an (3.23) into (3.19), we get I 2 ξ t (Dα ξ t ) 2 x ε(f X [ξ](t)+f[v](t)). (3.24) 14

15 t After summing α, where α 2, we have ( D α ξ t 2 x ) ξ (Dα ξ t ) 2 x +2 α 2 α =2 ε(f X [ξ](t)+f[v](t)). By (3.15)+(3.25), we obtain Thus, Lemma 3.3 is rove. α 2 D α v t 2 x (3.25) t F 1[ξ](t)+2F[v](t) ε(f X [ξ](t)+f[v](t)). (3.26) The following lemma concerns the estimate for L boun of F[v](t). Lemma 3.4. For any given T (,+ ), if su F[ξ,v,φ,ζ](t) ε, t T where < ε 1, then for t [,T], t F[v](t)+ γ a D α v 2 x+ γ a C εf[v](t). α 4 α 2 D α v t 2 x (3.27) Proof. Firstly, we estimate the L boun of v H4 ( ): Let D α v D α (2.4), where α 4, we get t D α v 2 = 2D α v D α [k 1 (1 γ)v( v) k 1 v v k1 2 ( v 2 )] [D α1 ( )Dα2 ( v)] D α v + 2γ a α 1 > + 2γ a Dα ( v) D α v. Integrate (3.28) in, we get t D α v 2 x = 2 D α v D α [k 1 (1 γ)v( v) k 1 v v k1 2 ( v 2 )]x + 2γ a [D α1 ( )Dα2 ( v)] D α vx α 1 > + 2γ a Dα ( v) D α vx := I 3. When α 3, 2γ a Dα v ( D α v)x = 2γ a Dα v 2 x 2γ a ( D α v)d α v ( )x C εf[v](t) 2γ a Dα v 2 x. (3.28) (3.29) (3.3) 15

16 Then t D α v 2 x+ 2γ a Dα v 2 x εf[v](t). (3.31) When α = 4, we estimate each term of I 3 searately. Assume the ositive 4 constants Λ,α,1,Λ,α,2,Λ,α,3,Λ,α,4 are so small that Λ,α,j γ a. The first term of I 3 : 2k 1 (1 γ) D α v D α [v( v)]x C εf[v](t)+λ,α,1 Dα v 2 x+c (v Dα v) 2 x C εf[v](t)+λ,α,1 Dα v 2 x. j=1 (3.32) The secon term of I 3 : 2k 1 D α v D α (v v)x C εf[v](t)+k 1 v D α v 2 x C εf[v](t). The thir term of I 3 : k 1 D α v D α ( v 2 )x = k 1 D α vd α ( v 2 )x Λ,α,2 Dα v 2 x+c Dα ( v 2 ) 2 x C εf[v](t)+λ,α,2 Dα v 2 x. The fourth term of I 3 : [D α1 ( )Dα2 ( v)] D α vx 2γ a 2γ a α 1 > + 2γ a + 2γ a α 1 =1 α 1 =2 α 1 =3 D α1 ( ) D α2 ( v) L 2 ( ) D α v L 2 ( ) D α1 ( ) D α2 ( v) L2 ( ) D α v L2 ( ) D α1 ( ) L 4 ( ) D α2 ( v) L 4 ( ) D α v L 2 ( ) (3.33) (3.34) (3.35) + 2γ a ( v) D α ( ) L 2 ( ) D α v L2 ( ) C εf[v](t)+λ,α,3 Dα v 2 x. The fifth term of I 3 : 2γ a Dα v ( D α v)x = 2γ a Dα v 2 x 2γ a ( D α v)d α v ( )x C εf[v]+( 2γ a +Λ,α,4) Dα v 2 x+c [D α v ( )]2 x C εf[v]+( 2γ a +Λ,α,4) Dα v 2 x. (3.36) 16

17 By (3.32)+(3.33)+(3.34)+(3.35)+(3.36), we have that when α = 4, t D α v 2 x+ γ a Dα v 2 x εf[v](t). (3.37) By (3.31) an (3.37), we have that when α 4, t α 4 D α v 2 x+ γ a α 4 Dα v 2 x εf[v](t). (3.38) Next, we estimate the L boun of v t H 2 ( ): Let D α v t D α t (2.4), where α 2, we get t D α v t 2 = 2D α v t D α t [k 1 (1 γ)v( v) k 1 v v k1 2 ( v 2 )] [D α1 ( )Dα2 ( v t )] D α v t + 2γ a + 2γ a α 1 > α 1 [D α1 t ( )Dα2 ( v)] D α v t + 2γ a Dα ( v t ) D α v t. Integrate (3.39) in, we get t D α v t 2 x = 2 D α v t D α t [k 1 (1 γ)v( v) k 1 v v k1 2 ( v 2 )]x + 2γ a [D α1 ( )Dα2 ( v t )] D α v t x α 1 > [D α1 t ( )Dα2 ( v)] D α v t x + 2γ a α 1 + 2γ a Dα ( v t ) D α v t x := I 4. When α 1, 2γ a Dα v t ( D α v t )x = 2γ a Dα v t 2 x 2γ a ( D α v t )D α v t ( )x C εf[v](t) 2γ a Dα v t 2 x. (3.39) (3.4) (3.41) Then t D α v t 2 x+ 2γ a Dα v t 2 x εf[v](t). (3.42) When α = 2, we estimate each term of I 4 searately. Assume the ositive 4 constants Λ 1,α,1,Λ 1,α,2,Λ 1,α,3,Λ 1,α,4 are so small that Λ 1,α,j γ a. j=1 17

18 The first term of I 4 : 2k 1 (1 γ) D α v t D α t [v( v)]x ( = 2k 1 (1 γ) D α v t D α [v t ( v)]x+ ) D α v t D α [v( v t )]x C εf[v](t)+λ 1,α,1 Dα v t 2 x+c (v Dα v t ) 2 x C εf[v](t)+λ 1,α,1 Dα v t 2 x. The secon term of I 4 : 2k 1 D α v t D α t (v v)x 2k 1 ( D α v t D α (v t v)x+ D α v t D α (v v t )x C εf[v](t)+k 1 C εf[v](t). v D α v t 2 x ) (3.43) (3.44) The thir term of I 4 : k 1 D α v t D α t ( v 2 )x = 2k 1 D α v t D α (v v t )x Λ 1,α,2 Dα v t 2 x+c Dα (v v t ) 2 x C εf[v](t)+λ 1,α,2 Dα v 2 x. The fourth term of I 4 : [D α1 ( )Dα2 ( v t )] D α v t x 2γ a 2γ a α 1 > α 1 =1 D α1 ( ) D α2 ( v t ) L2 ( ) D α v L2 ( ) + 2γ a Dα ( ) ( v t ) L 2 ( ) D α v L 2 ( ) C εf[v](t)+λ 1,α,3 Dα v t 2 x. (3.45) (3.46) The fifth term of I 4 : [D α1 t ( )Dα2 ( v)] D α v t x 2γ a α 1 2γ a t( ) D α ( v) L 2 ( ) D α v L 2 ( ) D α1 t ( ) L 4 ( ) D α2 ( v) L 4 ( ) D α v L 2 ( ) + 2γ a α 1 =1 (3.47) + 2γ a ( v) D α t ( ) L 2 ( ) D α v L2 ( ) C εf[v](t). 18

19 The sixth term of I 4 : 2γ a Dα ( v t ) D α v t x = 2γ a Dα v t 2 x 2γ a ( D α v t )D α v t ( )x C εf[v]+( 2γ a +Λ 1,α,4) Dα v t 2 x+c [D α v t ( )]2 x C εf[v]+( 2γ a +Λ 1,α,4) Dα v t 2 x. (3.48) By (3.43) + (3.44) + (3.45) + (3.46) + (3.47) + (3.48), we have that when α = 2, t D α v t 2 x+ γ a Dα v t 2 x εf[v](t). (3.49) t By (3.42) an (3.49), we have that when α 2, α 2 D α v t 2 x+ γ a α 2 Dα v t 2 x εf[v](t). (3.5) Sum (3.38)+(3.5), we obtain (3.27). Thus, Lemma 3.4 is rove. The structure of the equations (2.1) imlies F X [ξ](t) can be estimates by F[v](t), as the following lemma state: Lemma 3.5. For any given T (,+ ), if su F[ξ,v,φ,ζ](t) ε, t T there exists ε 2 >, which is ineenent of (ξ,v,φ ), such that if < ε min{1,ε 2 }, then for t [,T], for some c 3 >. F X [ξ](t) c 3 F[v](t), (3.51) Proof. Aly D α to ξ = k 1 v, where α 3, we get D α ξ = k 1 D α ( v), D α ξ 2 L 2 ( ) Dα ( v) 2 L 2 ( ) F[v](t)+F[ζ]F[v](t). (3.52) Aly D α to ξ t = k 1 v ξ k 1 γ v, where α 2, we get D α ξ t = k 1 D α (v ξ) k 1 γd α ( v), D α ξ t 2 L 2 ( ) Dα (v ξ) 2 L 2 ( ) + Dα ( v) 2 L 2 ( ) F[v](t)+F X [ξ]f[v](t). (3.53) By (3.52)+(3.53), we have F X [ξ](t) C 1 F[v](t)+C 1 F[ζ]F[v](t)+C 1 F X [ξ]f[v](t) C 1 F X [ξ]f[v](t)+c 1 (1+ε)F[v](t). (3.54) 19

20 for some C 1 >. Let ε 2 1 2C 1, c 3 = 2C 1 (1+ε 2 ) 2C 1 +1, then F X [ξ](t) c 3 F[v](t). Thus, Lemma 3.5 is rove. Remark 3.6. The roof of Lemma 3.5 inicates that uner the small ata assumtion su F[ξ,v,φ,ζ](t) ǫ, t T for some sufficiently small ε, we have for t [,T]. ξ 2 H 4 ( ) + ξ t 2 H 3 ( ) F[v](t)+F[ζ](t), F[v](t) F[ξ](t)+F[ζ](t), The following lemma gives not only the L bouns of F 1 [ξ](t),f[v](t), but also the boun of T F[ v](s)s. Lemma 3.7. For any given T (,+ ), there exists ε 3 > which is ineenent of (ξ,v,φ ), such that if su F[ξ,v,φ,ζ](t) ε, t T where ε min{1,ε 1,ε 2,ε 3 }, then for t [,T], F[ξ](t)+F[v](t) β 1 ( ξ H 5 ( ) + φ H 4 ( )), T T F[v](s)s+ F[ v](s)s β 2 ( ξ H 5 ( ) + φ H 4 ( )), (3.55) for some β 1 >,β 2 >. Proof. In view of Lemmas 3.3, 3.4 an 3.5, we have a riori estimates as follows: t F 1[ξ](t)+2F 1 [v](t) C 2 ǫ(fx [ξ](t)+f[v](t)), t F[v](t)+ γ a + γ a F X [ξ](t) c 3 F[v](t). α 4 α 2 By (3.56) 1 +(3.56) 2, we get D α v 2 x D α v t 2 x C 2 εf[v](t), (3.56) t (F 1[ξ](t)+F[v](t))+2F[v](t) + γ a D α v 2 x+ γ a α 4 C 2 ǫ(fx [ξ](t)+f[v](t)). α 2 D α v t 2 x (3.57) 2

21 lug (3.56) 3 into (3.57), we get [R.H.S. of (3.57)] C 2 (1+c 3 ) ǫf[v](t). (3.58) Take ε 3 = 1 C 2 2 (1+c3)2, when ε ε 3, we have t (F 1[ξ](t)+F[v](t))+F[v](t) + γ a D α v 2 x+ γ a α 4 α 2 D α v t 2 x. (3.59) Integrate (3.59) in, when t [,T], t F 1 [ξ](t)+f[v](t)+ F[v](s)s t + γ a α 4 D α v 2 xs+ γ a t α 2 D α v t 2 xs F 1 [ξ]()+f[v]() C 3 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ), (3.6) for some C 3 >. Take β 1 = C 3 min{c 1,1},β 2 = C 3 min{1, 4γ Thus, Lemma 3.7 is rove. 3a }. The following lemma concerns the L boun of F[φ](t) in [,T]. Lemma 3.8. For any given T (,+ ), if su F[ξ,v,φ,ζ](t) ε, t T where < ε min{1,ε 1,ε 2,ε 3 }, then for t [,T], an F[φ](t) has L boun in [,T]: t F[φ](t) β 3F[v](t) 1 2 F[φ](t), (3.61) F[φ](t) C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 3 (R 4 ))}, (3.62) for some β 3 >,C 6 >. Proof. Let D α φ D α (2.2) 3, where α 4, we get ( D α φ 2 ) t = 2 D α φd α1 v (D α2 φ) v D α φ 2. α 1 > Integrate (3.63) in, we have t D α φ 2 x = 2 α 1 > D α φd α1 v (D α2 φ)x+ D α φ 2 vx := I 5. (3.63) (3.64) 21

22 When α 4, it is easy to check I 5 F[v](t) 1 2F[φ](t). Sum α, we have t D α φ 2 x C 4 F[v](t) 1 2F[φ](t). (3.65) α 4 Let D α φ t D α t (2.2) 3, where α 2, we get ( D α φ t 2 ) t = 2 D α φ t D α1 v t (D α2 φ) α 1 2 α 1 > D α φ t D α1 v (D α2 φ t ) v D α φ t 2. (3.66) Integrate (3.66) in, we have t td l α φ 2 x = 2 D α φ t D α1 v t (D α2 φ)x 2 α 1 + D α φ t 2 vx := I 6. α 1 > D α φ t D α1 v (D α2 φ t )x (3.67) When α 2, it is easy to check I 6 F[v](t) 1 2F[φ](t). Sum α, we have t t ldα φ 2 x C 5 F[v](t) 1 2F[φ](t). (3.68) α 2 Let β 3 = C 4 +C 5, sum (3.65) an (3.68), we have t F[φ](t) β 3F[v](t) 1 2F[φ](t). (3.69) Since T (,+ ) is finite, integrate (3.68) from to t, when t [,T], we obtain the a riori estimate for F[φ](t). F[φ](t) C 6 φ 2 t H 4 ( ) ex{ β 3 F[v](s) 1 2 s} T C 6 φ 2 H 4 ( ) ex{β 3T F[v](s)s} C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )}. Thus, Lemma 3.8 is rove. (3.7) Due to ζ = (ξ +,φ + S), we can estimate F[ζ](t) in the following lemma: Lemma 3.9. For any given T (,+ ), if su F[ξ,v,φ,ζ](t) ε, t T where < ε min{1,ε 1,ε 2,ε 3 }, then for t [,T], F[ζ](t) C 7 C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )} +C 7 β 1 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ), (3.71) for some C 7 >. 22

23 Proof. Since ζ = γ 1 A ( +ξ) 1 γ ex{ S+φ γ } γ 1 A 1 γ ex{ S γ }, we have F[ζ](t) F[ξ](t)+F[φ](t). (3.72) The following lemma concerns L boun of F[ξ](t). Lemma 3.1. For any given T (,+ ), if su F[ξ,v,φ,ζ](t) ε, t T where < ε min{1,ε 1,ε 2,ε 3 }, then for t [,T], F[ξ](t) C 8 C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )} +(2C 8 +1)β 1 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ). (3.73) Proof. We have the L boun of F 1 [ξ](t), then for t [,T], ξ 2 L 2 ( ) β 1( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ). As Remark 3.6 sate, ξ 2 H 4 ( ) + ξ t 2 H 3 ( ) F[v](t) + F[ζ](t), while the L bouns of F[v](t) an F[ζ](t) in [,T] have been rove, then ξ 2 H 4 ( ) + ξ t 2 H 3 ( ) C 8 C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )} +2C 8 β 1 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ), for some C 8 >. Thus, Lemma 3.1 is rove. 4 Existence in [,T] of Classical Solutions to the Relaxe Equations In this section, we rove the existence in [,T] of classical solutions to the relaxe equations (2.2) uner small ata assumtion. After eliminating v from (2.2), we have the following arabolic-hyerbolic equations: ξ t = γ a ξ + a ξ φ, φ t = a ξ 1 φ, (ξ,φ)(x,) = (lim (x,τ), lim S (x,τ) S), τ τ where = ζ + = (ξ +,φ+ S). (4.1) The roof of the local existence of classical solutions to (4.1) is stanar (using the linearization-iteration-convergence scheme), so we give a lemma on the local existence without roof here. 23

24 Lemma 4.1. (Local Existence) If (lim ξ (x,τ), lim φ (x,τ)) H 5 ( ) H 4 ( ), inf lim τ τ (x,τ) >, then x τ there exists a finite time T >, such that Cauchy roblem (4.1) amits a unique local classical solution (ξ, φ) satisfying { (ξ,φ) l 1 Cl ([,T ),H 5 l ( ) H 4 l ( )), (4.2) ξ C( [,T )). The above lemma imlies the local existence of classical solutions to Cauchy roblem (2.2) for (ξ,φ) l 1 Cl ([,T ),C 2 l ( ) C 1 l ( )). Base on the a riori estimates for (ξ,v,φ,ζ), the solution of (2.2) can be extene from [,T ) to any finite time interval [,T]. Theorem 4.2. (Existence in [, T]) Assume (lim ξ (x,τ), lim φ (x,τ)) H 5 ( ) H 4 ( ), inf lim τ τ (x,τ) >. x τ There exists a sufficientlysmall number δ 1 >, suchthat if lim ξ (x,τ) 2 τ H 5 ( ) + lim φ (x,τ) 2 τ H 4 ( ) δ 1, then Cauchy roblem for the relaxe equations (2.2) amits a unique classical solution (ξ, v, φ, ζ) satisfying ξ l 1 Cl ([,T],H 5 l ( )), (v,φ,ζ) l 1 Cl ([,T],H 4 l ( ) 3 ), v, ξ C( [,T]). (4.3) Proof. In view of Lemmas 3.7,3.8,3.9,3.1, we have the following a riori estimates: for any given T (,+ ), if where < ε min{1,ε 1,ε 2,ε 3 }, then su F[ξ,v,φ,ζ](t) ε, (4.4) t T F[ξ](t)+F[v](t) β 1 ( ξ H5 ( ) + φ H4 ( )), E[φ](t) C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 3 (R 4 ) )}, E[ζ](t) C 7 C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )} +C 7 β 1 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ), (4.5) F[ξ](t) C 8 C 6 φ 2 H 4 ( ) ex{β 2β 3 T( ξ 2 H 5 ( ) + φ 2 H 4 ( ) )} +(2C 8 +1)β 1 ( ξ 2 H 5 ( ) + φ 2 H 4 ( ) ). Since ε is ineenent of ξ (x,τ) an φ (x,τ), there exists δ 1 > such that if lim ξ (x,τ) 2 τ H 5 ( ) + lim φ (x,τ) 2 τ H 4 ( ) δ 1, then lim ξ (x,τ) 2 τ H 5 ( ) + lim φ (x,τ) 2 τ H 4 ( ) min{ ε ε 2β 1, 4C 7β 1 }, lim τ φ (x,τ) 2 H 4 ( ) min{ ε 4C 6 ex{ β2β3tε 2β 1 }, (1 2C7)ε 4C 7 ex{ β2β3tε 2β 1 }}. (4.6) 24

25 Now, we can check a riori assumtion F[ξ,v,φ,ζ](t) ε is satisfie, then the valiity of the former a riori estimates is verifie. By Lemma 5.4, we have T E[ v](s)s β 2 ( ξ H 5 ( ) + φ H 4 ( )), (4.7) which imlies that for any given time T (,+ ), { v 2 L 2 ([,T],H 4 ( )) β 2( ξ H5 ( ) + φ H4 ( )), v t 2 L 2 ([,T],H 2 ( )) β 2( ξ H 5 ( ) + φ H 4 ( )). By Aubin-Lions Lemma, we obtain (4.8) v 2 C([,T],H 3 ( )) β 2( ξ H5 ( ) + φ H4 ( )), (4.9) which imlies that v C( [,T]) for any T (,+ ), then ξ = ak 1 v ak 1 v C( [,T]). (4.1) Due to the a riori estimates for (ξ,v,φ,ζ) an Lemma 4.1 on the local existence result, the solution of (2.2) can be extene from [,T ) to any finite time interval [,T]. Thus, Theorem 4.2 on the existence in [,T] of classical solutions to Cauchy roblem (2.2) is rove. 5 Uniform A Priori Estimates for the Relaxing Equations In this section, we erive unform a riori estimates for the relaxing Cauchy roblem (2.1). In orer to iscuss the initial layer an strong convergence of the velocity, we nee to estimate the higher orer time erivatives. Note that the notation of small quantity ǫ use in the following is ifferent from ε use before. The following lemma states that E[ξ](t) an E 1 [ξ](t) are equivalent, E[v](t) an E 1 [v](t) are equivalent. Lemma 5.1. For any given T (,+ ),τ (,1], there exists ǫ 1 > which is ineenent of (ξ,τv,φ ), such that if su E[ξ,τv,φ,ζ](t) ǫ 1, then ξ t T 3, ζ 2 an there exist c 4 >,c 5 > such that c 4 E[ξ](t) E 1 [ξ](t) c 5 E[ξ](t), c 4 E[v](t) E 1 [v](t) c 5 E[v](t). (5.1) ThefollowinglemmagivesuniformarioriestimateforE 1 [ξ](t)+τ 2 E 1 [v](t), which is equivalent to E[ξ](t)+τ 2 E[v](t). Lemma 5.2. For any given T (,+ ),τ (,1], if su E[ξ,τv,φ,ζ](t) ǫ, t T 25

26 where < ǫ 1, then for t [,T], t E 1[ξ](t)+τ 2 t E 1[v](t)+2E 1 [v](t) C ǫ(e X [ξ](t)+e[v](t)). (5.2) Proof. Let (2.1) (ξ,v), we get ( ξ 2 +τ 2 v 2 ) t +2k 2 ξ v+2k 2 v ξ +2 v 2 = 2γk 1 ξ 2 v 2k 1 ξv ξ 2k 1 τ 2 v v v + 2 k 1 ( 1 1 ) ξ v. (5.3) Integrate (5.3) in an note that (ξv)x =, we get t ξ 2 +τ 2 v 2 x+2 v 2 x = 2γk 1 v (ξ 2 ) 2k 1 ξv ξ 2k 1 τ 2 v v v+ 2 k 1 ( 1 1 ) ξ vx ǫ ξ L 2 ( ) v L 2 ( ) + ǫ v 2 L 2 ( ) ǫ(e X [ξ](t)+e[v](t)). (5.4) Aly td l α to (2.1), where l 2,1 l + α 4, then eliminate τ from both sies of the equations, we get ( td l α ξ) t +k 2 ( td l α v) = γk 1 td l α (ξ v) k 1 td l α (v ξ), τ 2 ( t l Dα v) t +k 2 ( t ldα ξ)+ t ldα v = k 1 τ 2 td l α (v v)+ 1 k 1 td l α [( 1 1 ) ξ]. (5.5) Let (5.5) ( l t Dα ξ, l t Dα v), we get ( l t Dα ξ 2 +τ 2 l t Dα v 2 ) t +2k 2 l t Dα ξ ( l t Dα v) +2k 2 l t Dα v ( l t Dα ξ)+2 l t Dα v 2 = 2γk 1 ( l td α ξ) l td α (ξ v) 2k 1 ( l td α ξ) l td α (v ξ) 2k 1 τ 2 ( td l α v) td l α (v v)+ 2 k 1 ( td l α v) td l α [( 1 1 ) ξ]. (5.6) t Integrate (5.6) in an note that ( l t Dα ξ l t Dα v)x =, we get td l α ξ 2 +τ 2 td l α v 2 x+2 td l α v 2 x = 2γk 1 ( t ldα ξ) t ldα (ξ v) 2k 1 ( t ldα ξ) t ldα (v ξ) 2k 1 τ 2 ( l td α v) l td α (v v)+ 2 k 1 ( l td α v) l td α [( 1 1 ) ξ]x := I 7. (5.7) When 1 l+ α 3, it is easy to check that I 7 ǫ(e X [ξ](t)+e[v](t)). When l + α = 4, we estimate the quantity I 7 ξ t ( l td α ξ) 2 x + 26

27 ( 1) l t Dα v 2 x, then I 7 ξ t ( l t Dα ξ) 2 x+τ 2 t ( 1) l t Dα v 2 x = 2γk 1 ( td l α ξ)ξ ( td l α v)x 2k 1 ( td l α ξ)v ( td l α ξ)x 2k 1 τ 2 v ( td l α v) ( td l α v)x+ 2 k 1 ( 1 1 )( l td α v) ( td l α ξ)x +2τ 2 ( 1)( l t Dα v) ( t ldα v t )x+τ 2 ζ t t ldα v 2 x τ 2 t 2 ξ ( l t Dα ξ)( t ldα ξ t )x t ( ξ )( l t Dα ξ) 2 x ǫ(e X [ξ](t)+e[v](t))+k 1 ( td l α ξ 2 +τ 2 td l α v 2 ) vx 2γk 1 ξ( td l α ξ) ( td l α v)x+ 2 k 1 ( 1 1 )( l td α v) ( td l α ξ)x 2 ξ ( l t Dα ξ)( t ldα ξ t )x+2τ 2 ( 1)( l t Dα v) ( t ldα v t )x 2 ξ ( l td α ξ)[ l td α ξ t +k 1 γ ( l td α v)]x + 2 k 1 ( 1 1 )( l t Dα v) [ ( l t Dα ξ)+k 1 τ 2 ( l t Dα v t )]x + ǫ(e X [ξ](t)+e[v](t)). Aly l t Dα to (2.1) 1, where l 2,l+ α = 4, we get td l α ξ t +k 1 γ ( td l α v) = k 1 td l α (v ξ) k 1 γ t l1 ξ ( Dα1 t l2 v). Dα2 l 1+ α 1 > (5.8) (5.9) Plug (5.9) into the following integral, we get ξ ( l t Dα ξ)( t ldα ξ t +k 1 γ ( t ldα v))x = ξ ( l t Dα ξ)[r.h.s. of (5.9)]x ǫ(e X [ξ](t)+e[v](t))+ k1 2 t ldα ξ 2 ( ξ v)x ǫ(e X [ξ](t)+e[v](t)). (5.1) Aly td l α to k 1 τ 2 v t + k1τ 2 2 v v + k 1 v + 1 ξ =, where l 2,l+ α = 4, we get 1 ( l td α ξ)+k 1 τ 2 ( l td α v t ) = k 1 l td α v k 2 1τ 2 l td α (v v) l 1+ α 1 > l1 t D α1 ( 1 ) l2 t D α2 ξ. (5.11) Plug (5.11) into the following integral, we get ( 1 1 )( l td α v) [ ( l td α ξ)+k 1 τ 2 l td α v t ]x = ( 1 1 )( l t Dα v) [R.H.S. of (5.11)]x k2 1 τ 2 l t Dα v 2 [( 1)τv]x+ ǫ(e X [ξ](t)+e[v](t)) ǫ(e X [ξ](t)+e[v](t)). (5.12) 27

28 Plug (5.1) an (5.12) into (5.8), we get I 7 t t ξ ( l td α ξ) 2 x+τ 2 t After Summing l an α, we have ( t ldα ξ 2 +τ 2 t ldα v 2 x l 2, l+ α 4 ξ ( l t Dα ξ) 2 x+τ 2 +2 l 2,l+ α =4 l 2, l+ α 4 ǫ(e X [ξ](t)+e[v](t)). t Then ( +τ 2 t +2 ( l 2, l+ α 4 ( l 2, l+ α 4 l 2, l+ α 4 ǫ(e X [ξ](t)+e[v](t)). l t Dα v 2 x l t Dα ξ 2 x ( 1) l td α v 2 x ǫ(e X [ξ](t)+e[v](t)). l t Dα v 2 x+ l td α v 2 x+ By (5.4)+(5.15), we obtain l 2,l+ α =4 l 2,l+ α =4 ( l 2,l+ α =4 (5.13) ) ( 1) l t Dα v 2 x ξ ( l t Dα ξ) 2 x ( l 2,l+ α =4 ) (5.14) 1) td α v 2 x 1) l td α v 2 x ) ) (5.15) t E 1[ξ](t)+τ 2 t E 1[v](t)+2E 1 [v](t) C ǫ(e X [ξ](t)+e[v](t)). (5.16) Thus, Lemma 5.2 is rove. The structure of the equations (2.1) imlies E X [ξ](t) can be estimates by E[v](t), as the following lemma state: Lemma 5.3. For any given T (,+ ),τ (,1], if su E[ξ,τv,φ,ζ](t) ǫ, t T there exists ǫ 2 >, which is ineenent of (ξ,v,φ ), such that if < ǫ min{1,ǫ 2 }, then for t [,T], E X [ξ](t) c 6 E[v](t), (5.17) Proof. Aly D α to ξ = k 1 τ 2 v t k 2 1 τ2 v v k 1 v, where α 3, we 28

29 get D α ξ = k 1 τ 2 D α ( v t ) k1 2τ2 D α ( v v) k 1 D α ( v), D α ξ 2 L 2 ( ) τ4 D α ( v t ) 2 L 2 ( ) +τ4 D α ( v v) 2 L 2 ( ) + Dα ( v) 2 L 2 ( ) D α v t 2 L 2 ( ) + Dα v 2 L 2 ( ) +E[ζ]E[v](t) +E[τv](t)E[v](t) +E[τ v](t)e[ζ]e[v](t). (5.18) Aly D α to ξ t = k 1 v ξ k 1 γ v, where α 3, we get D α ξ t = k 1 D α (v ξ) k 1 γd α ( v), D α ξ t 2 L 2 ( ) Dα (v ξ) 2 L 2 ( ) + Dα ( v) 2 L 2 ( ) D α v 2 L 2 ( ) +E X[ξ]E[v](t). (5.19) Aly D α to ξ tt = k 1 v t ξ k 1 v ξ t k 1 γ v t k 1 γξ t v, where α 2, we get D α ξ tt = k 1 D α (v t ξ) k 1 D α (v ξ t ) k 1 γd α ( v t ) k 1 γd α (ξ t v), D α ξ tt 2 L 2 ( ) Dα v t 2 L 2 ( ) +E X[ξ]E[v](t). By (5.18)+(5.19)+(5.2), we have (5.2) E X [ξ](t) C 8 E[v](t)+C 8 E[ζ]E[v](t)+C 8 E[τv](t)E[v](t) +C 8 E[τv](t)E[ζ]E[v](t) +C 8 E X [ξ]e[v](t) C 8 E X [ξ]e[v](t)+c 8 (1+2ǫ+ǫ 2 )E[v](t). (5.21) for some C 8 > Let ǫ 2 1 2C 8, c 6 = 2C 8 (1 + 2ǫ 2 + ǫ 2 2 ) 2C C 8. Then E X [ξ](t) c 6 E[v](t). Thus, Lemma 5.3 is rove. Base on the above a riori estimates, we rove not only the uniform L boun of E[ξ,τv](t), but also the uniform boun of E[v](s)s. Lemma 5.4. For any given T (,+ ),τ (,1], if su E[ξ,τv,φ,ζ](t) ǫ, t T there exists ǫ 3 >, which is ineenent of (ξ,v,φ ), such that if < ǫ min{1,ǫ 1,ǫ 2,ǫ 3 }, then for t [,T], E[ξ,τv](t) β 4, T E[v](s)s β 4, (5.22) for some β 4 > which is ineenent of τ. 29

30 Proof. In view of Lemmas 5.2 an 5.3, we have a riori estimates as follows: { t E 1[ξ](t)+τ 2 t E 1[v](t)+2E 1 [v](t) C 9 ǫ(ex [ξ](t)+e[v](t)), (5.23) E X [ξ](t) c 6 E[v](t). Plug (5.23) 2 into (5.23) 1, we get t E 1[ξ](t)+τ 2 t E 1[v](t)+2E 1 [v](t) C 9 (1+c 6 ) ǫe[v](t) ǫe1 [v](t). C9(1+c6) c 4 (5.24) Take ǫ 3 = c 2 4 C 2 9 (1+c6)2, when ǫ ǫ 3, t E 1[ξ](t)+τ 2 t E 1[v](t)+E 1 [v](t). (5.25) Integrate (5.25) from to t, we get E 1 [ξ](t)+τ 2 E 1 [v](t)+ t E 1 [v](s)s E 1 [ξ ](t)+e 1 [τv ](t) C 1 (ξ,u ) 2 H 4 ( ), (5.26) then t c 4 E[ξ](t)+c 4 E[τv](t) +c 4 E[v](s)s C 1 (ξ,u ) 2 H 4 ( ). (5.27) Thus, t E[ξ,τv](t)+ E[v](s)s C1 c 4 (ξ,u ) 2 H 4 ( ) β 4, (5.28) where β 4 > is ineenent of τ, since the initial ata are sufficiently small. Thus, Lemma 5.4 is rove. The following lemma concerns the uniform boun of E[φ](t). Here, the finiteness of T lays a key role in the roof. Lemma 5.5. For any given T (,+ ),τ (,1], if su E[ξ,τv,φ,ζ](t) ǫ, t T where < ǫ min{1,ǫ 1,ǫ 2,ǫ 3 }, then for t [,T], for some β 5 > which is ineenent of τ. E[φ](t) β 5, (5.29) Proof. Let l td α φ l td α (2.1) 3, where l 2, l+ α 4, we get ( td l α φ 2 ) t = 2 l 1+ α 1 > l td α φ l1 t D α1 v ( l2 t D α2 φ) v l td α φ 2. (5.3) 3

31 Integrate (5.3) in, we have t td l α φ 2 x = 2 td l α φ t l1 D α1 v ( t l2 D α2 φ)x l 1+ α 1 > + l t Dα φ 2 vx := I 8. (5.31) When l + α 4, it is easy to check I 8 E[v](t) 1 2E[φ](t). Sum l,α, we have t E[φ](t) C 11E[v](t) 1 2E[φ](t) (5.32) Integrate (5.32) in (,t), where t [,T], we obtain the uniform a riori estimate for E[φ](t) which is ineenent of τ. Here, τ > is variant, even τ aroaches +. E[φ](t) C 12 φ 2 t H 4 ( ) ex{ C 11 E[v](s) 1 2 s} T C 12 φ 2 H 4 ( ) ex{c 11T E[v](s)s} C 12 φ 2 H 4 ( ) ex{c 11Tβ 4 } β 5, (5.33) where β 5 > is ineenent of τ, since the initial ata are sufficiently small. Thus, Lemma 5.5 is rove. Due to ζ = (ξ +,φ + S), we can estimate E[ζ](t) in the following lemma: Lemma 5.6. For any given T (,+ ),τ (,1], if su E[ξ,τv,φ,ζ](t) ǫ, t T where < ǫ min{1,ǫ 1,ǫ 2,ǫ 3 }, then for t [,T], for some β 6 > which is ineenent of τ. Proof. E[ζ](t) CE[ξ](t)+CE[φ](t) E[ζ](t) β 6, (5.34) CC 12 φ 2 H 4 ( ) ex{c 11β 4 T (ξ,v ) 2 H 4 ( ) }+Cβ 4 (ξ,v ) 2 H 4 ( ) β 6, (5.35) where β 6 > is ineenent ofτ, since the initial ata aresufficiently small. By now, we have obtaine the L boun of E[ξ,τv,φ,ζ](t), namely, E[ξ](t) β 4, E[τv](t) β 4, E[φ](t) β 5, E[ζ](t) β 6. These bouns are ineenent of τ. Thus, as long as the initial ata are sufficiently small, the a riori assumtion su E[ξ,τv,φ,ζ](t) ǫ is vali. t T 31

32 6 Initial Layer an Relaxation Limit of the Relaxing Equations This section concerns the relaxation limit of the relaxing Cauchy roblem (2.1) an the initial layer for the velocity. In finite time interval [,T], the bouns of E[ξ](t),E[τv](t),E[φ](t),E[ζ](t) an T E[v](s)s are uniform with resect to τ, thus we can ass to the limit. The following theorem states the relaxation limit of the relaxing Cauchy roblem (2.1). Theorem 6.1. Suose that the initial ata for the relaxing Cauchy roblem 1 (2.1) satisfy ( (x,τ), k u 1τ (x,τ),s (x,τ) S) H 4 ( ), inf x (x,τ) >, inf lim 1 (x,τ) >, ( (x,τ), x τ k u 1τ (x,τ),s (x,τ) S) H4 ( ) is sufficiently small. Then for any finite T >, the roblem (2.1) amits a unique solution (ξ,v,φ,ζ) in [,T] satisfying such that as τ, l t ξ, τ l t v, l t φ, l t ζ L ([,T],H 4 l ( )), l 2, ξ, v, φ, ζ l 2 Hl ([,T],H 4 l ( )), (6.1) ξ ξ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, φ ξ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, ζ ζ in C([,T],C 2+µ1 (K) W 3,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6, v ṽ in l 2 Hl ([,T],H 4 l ( )), (6.2) for some function ( ξ,ṽ, φ, ζ) which is a weak solution to the relaxe equations (2.2), where K is arbitrary comact subset of. ( ξ,ṽ, φ, ζ) is the classical solution to (2.2), if (lim ξ (x,τ), lim φ (x,τ)) H 5 ( ) H 4 ( ) is satisfie.. τ τ Proof. In view of Lemmas 5.4, 5.5, 5.6, we have uniform regularities (6.1) uner small ata assumtion, then the solution to (1.2), i.e., (ˆ,û,Ŝ, ˆ )(x,t ) = (ξ +,k 1 τv,φ+ S,ζ + )(x,t) has uniform bouns with resect to τ: su E[ˆ,û,Ŝ S, ˆ ](t )+ 1 τ t T/τ T = su E[ξ,k 1 τv,φ,ζ](t) +k1 2 E[v](t)t t T (2+k 2 1 )β 4 +β 5 +β 6. T/τ E[û](t )t (6.3) By a riori estimate (6.3) an the local existence of (1.2) which is stanar (see [14]), we have the existence of classical solutions to (1.2) in the time interval [,T/τ], then the solutions of the relaxing equations (1.1) in [,T] can be constructe via the rescaling of variables (1.3). 32

33 Theexistenceoftherelaxingequations(1.1)in[,T]isuniformwithresect to τ, ue to the uniform regularities (1.7) an the fact that T is ineenent of τ. Thus, the uniform existence in [,T] an uniform regularities (1.7) for the relaxing Cauchy roblem (1.1) are rove. By Aubin s Lemma (see [19]): L ([,T],H 4 ( )) H 1 ([,T],H 3 ( )) C([,T],C 2+µ1 (K) W 3,µ2 (K)), where K is any comact subset of, then we have the convergence results as (6.2) state. For any ϕ 1 H 1 (K [,T)), it follows from (2.1) 1 that t t ϕ 1 ξ +γk 1 v (ξ ϕ 1 +ϕ 1 ξ)+γk 1 v ϕ 1 k 1 ϕ 1 v ξxs = ξ (x,τ)ϕ 1 (,)x. (6.4) In view of convergenceresults (6.2), let τ, we have ξ ξ, ξ ξ, v ṽ, then t t ϕ 1 ξ +γk1 ṽ ( ξ ϕ 1 +ϕ 1 ξ)+γk 1 ṽ ϕ 1 k 1 ϕ 1 ṽ ξxs = lim ξ (x,τ)ϕ 1 (,)x. τ (6.5) t For any ϕ 2 H 1 (K [,T))3, it follows from (2.1) 2 that v ϕ 2 1 ξ 1 k ϕ 2xs = τ 2 t v t ϕ 2 +k 1 v v ϕ 2 xs. (6.6) While as τ, [R.H.S. of (6.6)] τ 2 v t L 2 ( [,T)) ϕ 2 L 2 ( [,T)) +τ τv v L 2 ( [,T)) ϕ 2 L 2 ( [,T)). (6.7) Let τ in (6.6), we have t ṽ ϕ 2 1 k 1 ξ ϕ 2 xs =. (6.8) For any ϕ 3 H 1 (K [,T)), it follows from (2.1) 3 that t φ t ϕ 3 k 1 v φϕ 3 xs = φ (x,τ)ϕ 3 (,)x. (6.9) Let τ in (6.9), we have t φ t ϕ 3 k 1 ṽ φϕ 3 xs = lim τ φ (x,τ)ϕ 3 (,)x. (6.1) 33

34 While in arbitrary K, +ξ + ξ, A(ζ + ) γ ex{ S +φ} A( ζ + ) γ ex{ S + φ}, +ξ = A(ζ + ) γ ex{ S +φ}, imly = + ξ = A( ζ + ) γ ex{ S + φ} = A γe S, in K. (6.11) Since K an test functions ϕ 1, ϕ 2,ϕ 3 are arbitrary, (6.5), (6.8), (6.1) an (6.11) imly ( ξ,ṽ, φ, ) is the weak solution of the relaxe equations (2.2) with initial ata (lim ξ, lim φ ). τ τ As rove in Theorem 4.2, Cauchy roblem for the relaxe equations (2.2) with small ata (lim ξ, lim φ ) H 5 ( ) H 4 ( ) amits a unique classical τ τ solution. By the uniqueness of the weak solutions to (2.2), ( ξ,ṽ, φ, ) is classical solution of the relaxe equations (2.2) if the ata satisfy (lim ξ, lim φ ) τ τ H 5 ( ) H 4 ( ). Thus, Theorem 6.1 is rove. Besies the weak convergence result of the velocity, we have the strong convergence of the velocity outsie the initial layer for the ill-reare ata an strong convergence of the velocity near t = for the well-reare ata, as the following theorem state. Theorem 6.2. Let (ξ,v,φ,ζ) an ( ξ,ṽ, φ, ζ) be the solutions obtaine in Theorem 6.1, an K is any comact subset of. For the ill-reare ata, i.e., 1 lim v (x,τ)+ τ k ξ 1 (x,τ) (x,τ), there exists an initial layer [,t ] with < t < Cτ 2 δ for the velocity u, where C >, δ >, such that as τ, v ṽ in C((,T],C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6; v 1 for the well-reare ata, i.e., lim (x,τ)+ τ k ξ 1 (x,τ) (x,τ) as τ, H 2 ( ) =, v ṽ in C([,T],C +µ1 (K) W 1,µ2 (K)),µ 1 (, 1 2 ),2 µ 2 < 6. Proof. Aly D α to (2.9), where α 2, we have (D α η) t + 1 τ D α η = D α [v ( 1 ) ξ 1 ( v) ξ γ ξ v 2 γ ( v)+ 1 2(v ζ + v) ξ] k 1 D α 1 v D α2 η. (6.12) Let (6.12) D α η, we get ( D α η 2 ) t + 2 τ D α η 2 = 2D α η D α [v ( 1 ) ξ 1 ( v) ξ γ ξ v 2 γ ( v)+ 1 2(v ζ + v) ξ] 2k 1 D α 1 v D α2 η D α η. (6.13) 34

35 After integrating (6.13) in, we have t D α η 2 x+ 2 τ D α η 2 x 2 = 2 D α η D α [v ( 1 ) ξ 1 ( v) ξ γ ξ v γ ( v) + 1 2(v ζ + v) ξ]x k 1 v D α η 2 x 2k 1 D α1 v D α2 η D α ηx α 1> R 3 1 2τ D α η 2 x+2τ 2 D α [v ( 1 ) ξ 1 ( v) ξ γ ξ v 2 γ ( v)+ 1 2(v ζ + v) ξ] 2 L 2 ( ) +k 1 v D α η 2 x 2k 1 D α1 v D α2 η D α ηx α 1> R 3 1 2τ D α η 2 x+c D α ξ 2 2 L 2 ( ) +C Dα φ 2 L 2 ( ) +C D α ζ 2 L 2 ( ) [ τdα v 6 L 6 ( ) + Dα ( 1 ) 6 L 6 ( ) +(2+γ) ( 1 ) 6 L 6 ( ) + τdα ( v) 6 L 6 ( ) +(3+γ) Dα ξ 6 L 6 ( ) +(1+γ) τd α ( v) 6 L 6 ( ) + Dα ( τv 2) 6 L 6 ( ) + Dα ζ 6 L 6 ( ) ] +C τ ( v) ( ) 2 L 2 ( ) +C τd α ( v) 2 L 2 ( ) + 1 4τ D α η 2 x 4k 2 1 τ τ v 2k 1 4τ τd α1 1 v 4τ D α2 η 2 x 2 α 1> 1 τ D α η 2 x+c 2 13 E[ξ,τv,φ,ζ](t), (6.14) where < τ < min{1,τ } is small enough such that 4k 1 τ τ v +8k 1 τ τd α1 v CτE[τv](t) 1 2 Cτ ǫ 1. α 1> Thus, t t τ 2 D α η 2 x+ 1 D α η 2 x C 13 E[ξ,τv,φ,ζ](t), R( 3 ) ex{ t τ 2 } D α η 2 x C 13 ex{ t τ 2 }E[ξ,τv,φ,ζ](t). (6.15) After integrating from to t, we get D α η 2 x ex{ t τ } D α η 2 2 t x+c 13 ex{ t τ } D α η 2 2 x+c 13 ǫ t ex{ t τ } D α η 2 2 x+c 13 ǫτ 2. ex{ t s τ 2 }E[ξ,τv,φ,ζ](s)s ex{ t s τ 2 }s (6.16) For any fixe t (,T), when t t, we have D α η 2 x ex{ t τ } D α η 2 2 x+c 13 ǫτ 2. (6.17) 35

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Consistency and asymptotic normality

Consistency and asymptotic normality Consistency an asymtotic normality Class notes for Econ 842 Robert e Jong March 2006 1 Stochastic convergence The asymtotic theory of minimization estimators relies on various theorems from mathematical

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Colin Cameron: Brief Asymptotic Theory for 240A

Colin Cameron: Brief Asymptotic Theory for 240A Colin Cameron: Brief Asymtotic Theory for 240A For 240A we o not go in to great etail. Key OLS results are in Section an 4. The theorems cite in sections 2 an 3 are those from Aenix A of Cameron an Trivei

More information

Colin Cameron: Asymptotic Theory for OLS

Colin Cameron: Asymptotic Theory for OLS Colin Cameron: Asymtotic Theory for OLS. OLS Estimator Proerties an Samling Schemes.. A Roama Consier the OLS moel with just one regressor y i = βx i + u i. The OLS estimator b β = ³ P P i= x iy i canbewrittenas

More information

7. Introduction to Large Sample Theory

7. Introduction to Large Sample Theory 7. Introuction to Large Samle Theory Hayashi. 88-97/109-133 Avance Econometrics I, Autumn 2010, Large-Samle Theory 1 Introuction We looke at finite-samle roerties of the OLS estimator an its associate

More information

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL

More information

Abstract A nonlinear partial differential equation of the following form is considered:

Abstract A nonlinear partial differential equation of the following form is considered: M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent

More information

ASYMPTOTICS TOWARD THE PLANAR RAREFACTION WAVE FOR VISCOUS CONSERVATION LAW IN TWO SPACE DIMENSIONS

ASYMPTOTICS TOWARD THE PLANAR RAREFACTION WAVE FOR VISCOUS CONSERVATION LAW IN TWO SPACE DIMENSIONS TANSACTIONS OF THE AMEICAN MATHEMATICAL SOCIETY Volume 35, Number 3, Pages 13 115 S -9947(999-4 Article electronically publishe on September, 1999 ASYMPTOTICS TOWAD THE PLANA AEFACTION WAVE FO VISCOUS

More information

Consistency and asymptotic normality

Consistency and asymptotic normality Consistency an ymtotic normality Cls notes for Econ 842 Robert e Jong Aril 2007 1 Stochtic convergence The ymtotic theory of minimization estimators relies on various theorems from mathematical statistics.

More information

Convergence of random variables, and the Borel-Cantelli lemmas

Convergence of random variables, and the Borel-Cantelli lemmas Stat 205A Setember, 12, 2002 Convergence of ranom variables, an the Borel-Cantelli lemmas Lecturer: James W. Pitman Scribes: Jin Kim (jin@eecs) 1 Convergence of ranom variables Recall that, given a sequence

More information

Lecture 6 : Dimensionality Reduction

Lecture 6 : Dimensionality Reduction CPS290: Algorithmic Founations of Data Science February 3, 207 Lecture 6 : Dimensionality Reuction Lecturer: Kamesh Munagala Scribe: Kamesh Munagala In this lecture, we will consier the roblem of maing

More information

Global Solutions to the Coupled Chemotaxis-Fluid Equations

Global Solutions to the Coupled Chemotaxis-Fluid Equations Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ

MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ GENERALIZED UNIFORMLY CONTINUOUS SEMIGROUPS AND SEMILINEAR HYPERBOLIC SYSTEMS WITH REGULARIZED DERIVATIVES MARKO NEDELJKOV, DANIJELA RAJTER-ĆIRIĆ Abstract. We aopt the theory of uniformly continuous operator

More information

Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig

Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Max-Planck-Institut für Mathematik in en Naturwissenschaften Leipzig Simultaneous Blowup an Mass Separation During Collapse in an Interacting System of Chemotactic Species by Elio Euaro Espejo, Angela

More information

THE ZEROS OF A QUADRATIC FORM AT SQUARE-FREE POINTS

THE ZEROS OF A QUADRATIC FORM AT SQUARE-FREE POINTS THE ZEROS OF A QUADRATIC FORM AT SQUARE-FREE POINTS R. C. BAKER Abstract. Let F(x 1,..., x n be a nonsingular inefinite quaratic form, n = 3 or 4. Results are obtaine on the number of solutions of F(x

More information

Problem set 2: Solutions Math 207B, Winter 2016

Problem set 2: Solutions Math 207B, Winter 2016 Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

(1.1) (n 2 u 2 ) t + div(n 2 u 2 u 2 ) + q(n 2 ) = n 2 Ψ n 2 u 2,

(1.1) (n 2 u 2 ) t + div(n 2 u 2 u 2 ) + q(n 2 ) = n 2 Ψ n 2 u 2, LAGE TIME BEHAVIO OF SOLUTIONS TO N-DIMENSIONAL BIPOLA HYDODYNAMIC MODEL FO SEMICONDUCTOS FEIMIN HUANG, MING MEI, AND YONG WANG Deicate to Professor Akitaka Matsumura on his 6th birthay Abstract. In this

More information

ON THE AVERAGE NUMBER OF DIVISORS OF REDUCIBLE QUADRATIC POLYNOMIALS

ON THE AVERAGE NUMBER OF DIVISORS OF REDUCIBLE QUADRATIC POLYNOMIALS ON THE AVERAGE NUMBER OF DIVISORS OF REDUCIBLE QUADRATIC POLYNOMIALS KOSTADINKA LAPKOVA Abstract. We give an asymtotic formula for the ivisor sum c

More information

Econometrics I. September, Part I. Department of Economics Stanford University

Econometrics I. September, Part I. Department of Economics Stanford University Econometrics I Deartment of Economics Stanfor University Setember, 2008 Part I Samling an Data Poulation an Samle. ineenent an ientical samling. (i.i..) Samling with relacement. aroximates samling without

More information

Mod p 3 analogues of theorems of Gauss and Jacobi on binomial coefficients

Mod p 3 analogues of theorems of Gauss and Jacobi on binomial coefficients ACTA ARITHMETICA 2.2 (200 Mo 3 analogues of theorems of Gauss an Jacobi on binomial coefficients by John B. Cosgrave (Dublin an Karl Dilcher (Halifax. Introuction. One of the most remarkable congruences

More information

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

A transmission problem for the Timoshenko system

A transmission problem for the Timoshenko system Volume 6, N., pp. 5 34, 7 Copyright 7 SBMAC ISSN -85 www.scielo.br/cam A transmission problem for the Timoshenko system C.A. RAPOSO, W.D. BASTOS an M.L. SANTOS 3 Department of Mathematics, UFSJ, Praça

More information

He s Homotopy Perturbation Method for solving Linear and Non-Linear Fredholm Integro-Differential Equations

He s Homotopy Perturbation Method for solving Linear and Non-Linear Fredholm Integro-Differential Equations nternational Journal of Theoretical an Alie Mathematics 2017; 3(6): 174-181 htt://www.scienceublishinggrou.com/j/ijtam oi: 10.11648/j.ijtam.20170306.11 SSN: 2575-5072 (Print); SSN: 2575-5080 (Online) He

More information

Discrete Operators in Canonical Domains

Discrete Operators in Canonical Domains Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

NONLINEAR QUARTER-PLANE PROBLEM FOR THE KORTEWEG-DE VRIES EQUATION

NONLINEAR QUARTER-PLANE PROBLEM FOR THE KORTEWEG-DE VRIES EQUATION Electronic Journal of Differential Equations, Vol. 11 11), No. 113, pp. 1. ISSN: 17-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu NONLINEAR QUARTER-PLANE PROBLEM

More information

Stability of steady states in kinetic Fokker-Planck equations for Bosons and Fermions

Stability of steady states in kinetic Fokker-Planck equations for Bosons and Fermions www.oeaw.ac.at Stability of steay states in kinetic Fokker-Planck equations for Bosons an Fermions L. Neumann, C. Sarber RICAM-Reort 2006-34 www.ricam.oeaw.ac.at STABILITY OF STEADY STATES IN KINETIC FOKKER-PLANCK

More information

Multiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type

Multiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type Nonlinear Analysis 7 29 536 546 www.elsevier.com/locate/na Multilicity of weak solutions for a class of nonuniformly ellitic equations of -Lalacian tye Hoang Quoc Toan, Quô c-anh Ngô Deartment of Mathematics,

More information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

1 Riesz Potential and Enbeddings Theorems

1 Riesz Potential and Enbeddings Theorems Riesz Potential and Enbeddings Theorems Given 0 < < and a function u L loc R, the Riesz otential of u is defined by u y I u x := R x y dy, x R We begin by finding an exonent such that I u L R c u L R for

More information

On the Cauchy Problem for Von Neumann-Landau Wave Equation

On the Cauchy Problem for Von Neumann-Landau Wave Equation Journal of Applie Mathematics an Physics 4 4-3 Publishe Online December 4 in SciRes http://wwwscirporg/journal/jamp http://xoiorg/436/jamp4343 On the Cauchy Problem for Von Neumann-anau Wave Equation Chuangye

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED

More information

Problem Set 2 Solution

Problem Set 2 Solution Problem Set 2 Solution Aril 22nd, 29 by Yang. More Generalized Slutsky heorem [Simle and Abstract] su L n γ, β n Lγ, β = su L n γ, β n Lγ, β n + Lγ, β n Lγ, β su L n γ, β n Lγ, β n + su Lγ, β n Lγ, β su

More information

Discharge initiation and plasma column formation in aspect ratio A=2 tokamak.

Discharge initiation and plasma column formation in aspect ratio A=2 tokamak. Discharge initiation an lasma column formation in asect ratio A toama... Khayrutinov E.A. Azizov, A.D. Baralov, G.G.Glaush, I.L.Taibaeva, Ph.W.West 3 Troits, Moscow eg., ussia NNC, K 3 General Atomics,

More information

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions Journal of Multivariate Analysis 99 2008) 1383 1392 www.elsevier.com/locate/jmva Bivariate istributions characterize by one family of conitionals an conitional ercentile or moe functions Barry C. Arnol

More information

1 Extremum Estimators

1 Extremum Estimators FINC 9311-21 Financial Econometrics Handout Jialin Yu 1 Extremum Estimators Let θ 0 be a vector of k 1 unknown arameters. Extremum estimators: estimators obtained by maximizing or minimizing some objective

More information

arxiv: v1 [math.pr] 28 Jan 2018

arxiv: v1 [math.pr] 28 Jan 2018 Submitte to the Annals of Probability WAVELET ANALYSIS OF THE BESOV REGULARITY OF LÉVY WHITE NOISES arxiv:80.09245v [math.pr] 28 Jan 208 By Shayan Azizneja, Julien Fageot an Michael Unser Ecole olytechnique

More information

On a class of nonlinear viscoelastic Kirchhoff plates: well-posedness and general decay rates

On a class of nonlinear viscoelastic Kirchhoff plates: well-posedness and general decay rates On a class of nonlinear viscoelastic Kirchhoff plates: well-poseness an general ecay rates M. A. Jorge Silva Department of Mathematics, State University of Lonrina, 8657-97 Lonrina, PR, Brazil. J. E. Muñoz

More information

By completing this chapter, the reader will be able to:

By completing this chapter, the reader will be able to: hater 4. Mechanics of Particles Particle mechanics governs many rinciles of article measurement instrumentation an air cleaning technologies. Therefore, this chater rovies the funamentals of article mechanics.

More information

THE ZERO-ELECTRON-MASS LIMIT IN THE HYDRODYNAMIC MODEL FOR PLASMAS

THE ZERO-ELECTRON-MASS LIMIT IN THE HYDRODYNAMIC MODEL FOR PLASMAS THE ZERO-ELECTRON-MASS LIMIT IN THE HYDRODYNAMIC MODEL FOR PLASMAS GIUSEPPE ALÌ, LI CHEN, ANSGAR JÜNGEL, AND YUE-JUN PENG Abstract. The limit of vanishing ratio of the electron mass to the ion mass in

More information

Digitally delicate primes

Digitally delicate primes Digitally elicate rimes Jackson Hoer Paul Pollack Deartment of Mathematics University of Georgia Athens, Georgia 30602 Tao has shown that in any fixe base, a ositive roortion of rime numbers cannot have

More information

NONLOCAL p-laplace EQUATIONS DEPENDING ON THE L p NORM OF THE GRADIENT MICHEL CHIPOT AND TETIANA SAVITSKA

NONLOCAL p-laplace EQUATIONS DEPENDING ON THE L p NORM OF THE GRADIENT MICHEL CHIPOT AND TETIANA SAVITSKA NONLOCAL -LAPLACE EQUATIONS DEPENDING ON THE L NORM OF THE GRADIENT MICHEL CHIPOT AND TETIANA SAVITSKA Abstract. We are studying a class of nonlinear nonlocal diffusion roblems associated with a -Lalace-tye

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

Sums of independent random variables

Sums of independent random variables 3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where

More information

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Arch. Rational Mech. Anal. Digital Object Ientifier (DOI) https://oi.org/1.17/s25-18-1336-z Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Siming He & Eitan Tamor Communicate

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Generalized Tractability for Multivariate Problems

Generalized Tractability for Multivariate Problems Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION

THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION MARTIN BURGER, MARCO DI FRANCESCO, AND YASMIN DOLAK Abstract. The aim of this paper is to iscuss the

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION Electronic Journal of Differential Equations, Vol. 216 (216), No. 138, pp. 1 7. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN

More information

Gradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices

Gradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices Graient flow of the Chapman-Rubinstein-Schatzman moel for signe vortices Luigi Ambrosio, Eoaro Mainini an Sylvia Serfaty Deicate to the memory of Michelle Schatzman (1949-2010) Abstract We continue the

More information

Submitted to the Journal of Hydraulic Engineering, ASCE, January, 2006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS

Submitted to the Journal of Hydraulic Engineering, ASCE, January, 2006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS Submitte to the Journal of Hyraulic Engineering, ASCE, January, 006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS Gary Parker 1, Member, ASCE an Horacio Toniolo ABSTRACT This note is evote to the correction

More information

On a class of Rellich inequalities

On a class of Rellich inequalities On a class of Rellich inequalities G. Barbatis A. Tertikas Dedicated to Professor E.B. Davies on the occasion of his 60th birthday Abstract We rove Rellich and imroved Rellich inequalities that involve

More information

Well-posedness of hyperbolic Initial Boundary Value Problems

Well-posedness of hyperbolic Initial Boundary Value Problems Well-poseness of hyperbolic Initial Bounary Value Problems Jean-François Coulombel CNRS & Université Lille 1 Laboratoire e mathématiques Paul Painlevé Cité scientifique 59655 VILLENEUVE D ASCQ CEDEX, France

More information

Martin Luther Universität Halle Wittenberg Institut für Mathematik

Martin Luther Universität Halle Wittenberg Institut für Mathematik Martin Luther Universität alle Wittenberg Institut für Mathematik Weak solutions of abstract evolutionary integro-ifferential equations in ilbert spaces Rico Zacher Report No. 1 28 Eitors: Professors of

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

LOCAL SOLVABILITY AND BLOW-UP FOR BENJAMIN-BONA-MAHONY-BURGERS, ROSENAU-BURGERS AND KORTEWEG-DE VRIES-BENJAMIN-BONA-MAHONY EQUATIONS

LOCAL SOLVABILITY AND BLOW-UP FOR BENJAMIN-BONA-MAHONY-BURGERS, ROSENAU-BURGERS AND KORTEWEG-DE VRIES-BENJAMIN-BONA-MAHONY EQUATIONS Electronic Journal of Differential Equations, Vol. 14 (14), No. 69, pp. 1 16. ISSN: 17-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu LOCAL SOLVABILITY AND BLOW-UP

More information

Location of solutions for quasi-linear elliptic equations with general gradient dependence

Location of solutions for quasi-linear elliptic equations with general gradient dependence Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Boundary regularity for elliptic problems with continuous coefficients

Boundary regularity for elliptic problems with continuous coefficients Boundary regularity for ellitic roblems with continuous coefficients Lisa Beck Abstract: We consider weak solutions of second order nonlinear ellitic systems in divergence form or of quasi-convex variational

More information

On some parabolic systems arising from a nuclear reactor model

On some parabolic systems arising from a nuclear reactor model On some parabolic systems arising from a nuclear reactor moel Kosuke Kita Grauate School of Avance Science an Engineering, Wasea University Introuction NR We stuy the following initial-bounary value problem

More information

A Spectral Method for the Biharmonic Equation

A Spectral Method for the Biharmonic Equation A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic

More information

L p Theory for the Multidimensional Aggregation Equation

L p Theory for the Multidimensional Aggregation Equation L p Theory for the Multiimensional Aggregation Equation ANDREA L. BERTOZZI University of California - Los Angeles THOMAS LAURENT University of California - Los Angeles AND JESÚS ROSADO Universitat Autònoma

More information

3 The variational formulation of elliptic PDEs

3 The variational formulation of elliptic PDEs Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle

More information

BLOW-UP OF SOLUTIONS FOR WEAKLY COUPLED SYSTEMS OF COMPLEX GINZBURG-LANDAU EQUATIONS

BLOW-UP OF SOLUTIONS FOR WEAKLY COUPLED SYSTEMS OF COMPLEX GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Euations Vol. 27 27 No. 96. 8. ISSN: 72-669. URL: htt://eje.math.txstate.eu or htt://eje.math.unt.eu BLOW-UP OF SOLUTIONS FOR WEAKLY COUPLED SYSTEMS OF COMPLEX GINZBURG-LANDAU

More information

Nonlinear Stability for Multidimensional Fourth-Order Shock Fronts

Nonlinear Stability for Multidimensional Fourth-Order Shock Fronts Arch. Rational Mech. Anal. 26 Digital Object Identifier DOI.7/s25-5-9-y Nonlinear Stability for Multidimensional Fourth-Order Shock Fronts Peter Howard & Changbing Hu Communicated by C.M. Dafermos Abstract

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

arxiv: v1 [math.ap] 31 Dec 2018

arxiv: v1 [math.ap] 31 Dec 2018 OPTIMAL DECAY FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITHOUT ADDITIONAL SMALLNESS ASSUMPTIONS arxiv:18111714v1 [mathap] 31 Dec 018 ZHOUPING XIN AND JIANG XU Abstract We investigate the large time

More information

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne

More information

Centrum voor Wiskunde en Informatica

Centrum voor Wiskunde en Informatica Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng

More information

Noether s theorem applied to classical electrodynamics

Noether s theorem applied to classical electrodynamics Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance

More information

A LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM

A LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM

More information

ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL. Qiang Du, Manlin Li and Chun Liu

ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL. Qiang Du, Manlin Li and Chun Liu DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 8, Number 3, October 27 pp. 539 556 ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL Qiang

More information

Remarks on Singularities, Dimension and Energy Dissipation for Ideal Hydrodynamics and MHD

Remarks on Singularities, Dimension and Energy Dissipation for Ideal Hydrodynamics and MHD Commun. Math. Phys. 184, 443 455 (1997) Communications in Mathematical Physics c Springer-Verlag 1997 Remarks on Singularities, imension an Energy issipation for Ieal Hyroynamics an MH Russel E. Caflisch,

More information

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075

More information

A Simple Exchange Economy with Complex Dynamics

A Simple Exchange Economy with Complex Dynamics FH-Kiel Universitf Alie Sciences Prof. Dr. Anreas Thiemer, 00 e-mail: anreas.thiemer@fh-kiel.e A Simle Exchange Economy with Comlex Dynamics (Revision: Aril 00) Summary: Mukherji (999) shows that a stanar

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Dissipative numerical methods for the Hunter-Saxton equation

Dissipative numerical methods for the Hunter-Saxton equation Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a

More information

EXISTENCE OF SOLUTIONS TO WEAK PARABOLIC EQUATIONS FOR MEASURES

EXISTENCE OF SOLUTIONS TO WEAK PARABOLIC EQUATIONS FOR MEASURES EXISTENCE OF SOLUTIONS TO WEAK PARABOLIC EQUATIONS FOR MEASURES VLADIMIR I. BOGACHEV, GIUSEPPE DA PRATO, AND MICHAEL RÖCKNER Abstract. Let A = (a ij ) be a Borel mapping on [, 1 with values in the space

More information