arxiv: v1 [math.ap] 4 Jan 2018

Transcription

1 A SCALING LIMIT FROM THE WAVE MAP TO THE HEAT FLOW INTO S 2 NING JIANG, YI-LONG LUO, SHAOJUN TANG, AND ARGHIR ZARNESCU arxiv: v math.ap] 4 Jan 208 Abstract. In this paper we study a limit connecting a scaled wave map with the heat flow into the unit sphere S 2. We show quantitatively how that the two equations are connected by means of an initial layer correction. This limit is motivated as a first step into understanding the limit of zero inertia for the hyperbolic-parabolic Ericksen-Leslie s liquid crystal model.. Introduction.. Wave map and heat flow. We consider a hyperbolic system for functions d : R + R 3 S 2 : t d = d+( d 2 t d 2 )d, (.) subject to initial data: for any x R 3, d t=0 = d 0 (x) S 2, t d t=0 = d 0 (x) R 3,d 0 (x) d 0 (x) = 0, (.2) where = tt is the standard wave operator, and the compatibility condition d 0 d 0 = 0 on the initial data is due to the fact that d 0. The system (.) is a wave map from R 3 to the unit sphere S 2, with a damping term t d. One way of interpreting this system is as follows: setting the righthand side of (.) equal to 0, we obtain d = ( d 2 t d 2 )d. This is the well-known wave map, which can be characterized variationally as a critical point of the functional A(d) = ( d 2 t d 2 )dxdt, (.3) 2 among maps d satisfying the target constraint, d : R + R 3 S 2. Thus the full system (.) can be viewed as a gradient flow of the functional (.3). Another gradient flow can be obtained by formally dropping some terms out of the previous system, and obtaining the heat flow t d = d+ d 2 d. (.4) Similarly as before setting the right-hand side equal to zero we obtain the equations for the harmonic map from R 3 to the unit sphere S 2 namely d+ d 2 d = 0, (.5) which is a critical point of the energy functional E(d) = ˆ d 2 dx. (.6) 2 There exist deep relations between the two systems, (.) and (.4) and one way to see this is by considering the following parabolic scaling: d (t,x) := d( t, x ), (.7) January 8, 208.

2 2 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU Then, d satisfies the following scaled wave map: t d = ( tt )d +( d 2 t d 2 )d, (.8) on R + R 3. For this scaled system we take the initial values independent of, namely: d t=0 = d in (x) S 2, t d t=0 = d in (x) R 3. (.9) The finite-time behaviour of the limit 0 for the system (.8) with initial data (.9) is the focus of this paper. It is easy to see that letting = 0 in (.8) will formally give the heat flow (.4). However, a refined analysis and the introduction of an initial layer is needed in order to overcome the singular character of this limit and understand the relationship between the system (.8) and its formal limit, as it will be seen in the Theorem. below. Notations and conventions: Throughout this paper, we use the following standard notations: ˆ d p L = d p dx, d p H k = γ d L 2, d Ḣk = γ d L 2. R 3 γ k γ k Additionally, for the Hilbert space L 2 L 2 (dx,r 3 ), we use the following notation to denote the standard inner product: ˆ f,g = fgdx. R 3 Furthermore, if there is a generic constant C > 0 such that the inequality f(t) Cg(t) holds for all t 0, we denote this inequality by f g..2. Ericksen-Leslie s hyperbolic liquid crystal model. Our motivation for considering the previously mentioned limit comes from the hydrodynamic theory of nematic liquid crystals. The most widely accepted equations of nematics were proposed by Ericksen 5, 6, 7] and Leslie 0, ] in the 960 s (see for more details Section 5. of 2] ). The general hyperbolicparabolic Ericksen-Leslie system consists of an equation for velocity u of the centers of mass of the rod-like molecules, coupled with an equation for the direction d of these molecules. More specifically we have the following equations (in non-dimensional form): t u+u u 2 µ 4 u+ p = div( d d)+div σ, divu = 0, (.0) Dud 2 = d+γd+λ (D u d Bd)+λ 2 Ad on R + R 3 with constraint d =, where A = 2 ( u+ u ) and B = 2 ( u u ), D u f = t f+u f is the material derivative of f respect to the vector u, D 2 ud = t D u d+u D u d. The Lagrangian γ that enforces the unit-length constraint d = is given by γ γ(u,d,d u d) = D u d 2 + d 2 λ 2 d Ad. The stress tensor appearing in the equation for u is given by: σ ij ( σ(u,d,d u d) ) ij =µ d k d p A kp d i d j +µ 2 d j ((D u d) i +B ki d k ) +µ 3 d i ((D u d) j +B kj d k )+µ 5 d j d k A ki +µ 6 d i d k A kj. The constant > 0 measures the inertial effects. The constants µ i ( i 6) are known as Leslie coefficients and one has µ 4 > 0. Furthermore, we have: λ = µ 2 µ 3,λ 2 = µ 5 µ 6,µ 2 +µ 3 = µ 6 µ 5, where the last relation is called Parodi relation. For the more background and derivation of (.0), see 0] and 9]. For any fixed > 0, in 9] the first two named authors of the current paper proved the local well-posedness of the system (.0) under assumptions on the Leslie coefficients which

3 WAVE MAP LIMIT 3 ensure the dissipativity of the basic energy law, and global well-posedness with small initial data under further damping effect, i.e. λ < 0. As noted inthe Conclusion section of9], theinertial constant > 0is, physically, inmost common non-dimensionalisations and materials, very small. Formally, letting = 0 will give the parabolic Ericksen-Leslie system which is basically a coupling of Navier-Stokes equations and an extension of the heat flow to the unit sphere. However it is a very challenging task to obtaining estimates uniform in for the full system (.0), in order to understand the limit 0. In the current paper, the problem we consider what appears as a simple instance of this general problem, namely the case where the bulk velocity u = 0 and the coefficient λ = in (.0). For this case, the system (.0) is reduced to the scaled wave map (.8), i.e. the wave map (.8) with a damping can be regarded as an Ericksen-Leslie s liquid crystal flow unaffected by the fluid velocity..3. Initial layer and the main result. As mentioned in the previous two subsections the formal limit of the equation (.8), obtained by setting = 0 is provided by the heat flow for functions with values into S 2 : t d 0 = d 0 + d 0 2 d 0,d 0 S 2, (.) The limit we consider is a singular limit, as the character of the equations changes, from a hyperbolic-typesystem for > 0 to a parabolic system for = 0. An immediate manifestation of the difference between the two types of equations is related to the initial conditions, which for the limit equation take the form: d 0 t=0 = d in (x) S 2. (.2) Thus, we note that the wave map is a system of hyperbolic equations with two initial conditions, while the heat flow is a parabolic system with only one initial condition. Usually the solution of the heat flow does not satisfy the second initial condition in(.9). This disparity between the initial conditions of the wave map (.8) and of the heat flow (.) indicates that in one should expect an initial layer in time, appearing in the limiting process 0. A formal derivation (postponed for later, in Section 2) indicates that this should be of the form: where D(x) is defined as d I 0( t,x) = ( d in (x) t d 0 (0,x))exp( t ) = ( d in (x) d in (x) d in (x) 2 d in (x))exp( t ) (.3) D(x)exp( t ), D(x) d in (x) d in (x) d in (x) 2 d in (x). Our study of the limit from the wave map (.8) to the heat flow (.5) is inspired by the classical approach of Caflisch on the compressible Euler limit of the Boltzmann equation ]. This approach is based on the Hilbert expansion in which the leading term is given by solutions of the limit equation. The Caflisch s approach assumes that a solution of the limiting equation (which in our case is the heat flow (.5)) is known beforehand. Then the solution to the original equation (which in our case is the wave map (.8)) can be constructed around the limiting equation with perturbations as expansions in powers of. Based on the arguments above and the formal analysis in Section 2, in the expansions, besides the heat flow, the leading term should also include an initial layer. More specifically, we take the following ansatz of the solution d to the system (.8): d (t,x) = d 0 (t,x)+d I 0( t,x) + d R(t,x), (.4)

4 4 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU where d 0 (t,x) obeys the heat flow (.) and the initial layer d I t 0(,x) is defined in (.3). Plugging (.4) into the system (.8), the remainder term d R (t,x) must satisfy the system with the initial conditions where the singular term S(d R ) is tt d R + td R d R = S(d R )+R(d R ) (.5) d R(0,x) = D(x), t d R(0,x) = 0, (.6) S(d R ) = tt d 0 D(x)exp( t )+ d 0 2 d R + d R 2 d 0 t d 0 +D(x)exp( t ) 2 d 0 d 0 2 D(x)exp( t ) + 2 ( d 0 d R)d R 2 ( d 0 D(x))exp( t )d ( d 0 d R)d 0, and the regular term R(d R ) is R(d R ) = t d 0 +D(x)exp( t )+ t d 2 ] R D(x)exp( t )+d R 2( t d 0 +D(x)exp( t )) td R + t d R 2] d 0 2( D(x) d R )exp( t )d 0 + D(x) 2 exp( 2t )d 0 + d R 2 d R 2( d 0 D(x))exp( t )d R 2( d 0 d R)D(x)exp( t ) 2 ( D(x) d R)exp( t )d R d R 2 D(x)exp( t )+ D(x) 2 exp( 2t )d R 2 D(x) 2 3 exp( 3t )D(x) +2 ( d 0 D(x))exp( 2t )D(x)+2( D(x) d R)exp( 2t )D(x). According to Eells-Sampson s classical result in 4], for the heat flow (.) on the unit sphere S 2, one can have the following results of local well-posedness: Proposition.. For any given d in S 2 satisfying d in Ḣk (R 3 ) for any integer k > 2, there exists a time T = T( d in Ḣk) > 0 such that (.) admits a unique classical solution d 0 L (0,T;Ḣk ) L 2 (0,T;Ḣk+ ). Moreover, there is a constant C = C ( d in Ḣk,T) > 0 such that the solution d 0 satisfies d 0 2 L (0,T;Ḣk ) + d 0 2 L 2 (0,T;H k ) C. The proof can be found in Chapter 5 in the book 3]. Now we state the main result of this paper: Theorem.. We consider vector fields d in : R 3 S 2 and d in : R 3 R 3, satifying the compatibility condition d in d in 0. Assume that d in H 6, d in H 5, and let T > 0 be the time interval of existence of the solution of the heat flow (.) with initial condition d in, determined in Proposition.. Then, there exists an 0 0 ( d in H 6, d in H 5,T) (0, 2 ) such that for all (0, 0) we have that on the interval 0,T] the wave map equation (.8) with the initial conditions (.9) admits a unique solution with the form (.4), i.e., d (t,x) = d 0 (t,x)+d I 0( t,x) + d R(t,x), where d 0 is the solution of the heat flow (.) with initial condition d in and d I 0 (t,x) is the initial layer (.3). Moreover, there exists a positive constant C 0 = C 0 (d in, d in,t) > 0, such that the remainder term d R satisfies the bound for all (0, 0 ). t d R 2 L (0,T;H 2 ) + d R 2 L (0,T;H 3 ) C 0 (.7)

5 WAVE MAP LIMIT 5 Remark.. The rate of convergence we obtain is optimal. Indeed, in order to see this, it suffices to note that the limit we study contains as a particular case the linear scalar case of the singular limit of the damped wave equation to the heat equation. Indeed, let us consider a solution of the scalar damped wave equation: tt θ + t θ = θ (.8) for θ : R R 3 R with initial datas: θ (0,x) = θ 0 (x), t θ (0,x) = θ (x) (.9) all smooth functions. Also consider the solution of the heat equation: for θ 0 : R R 3 R with initial data: t θ 0 = θ 0 (.20) θ 0 (0,x) = θ 0 (x), (.2) all smooth functions. Denoting n 0 (t,x) := (cosθ 0 (t,x),sinθ 0 (t,x),0) and n (t,x) := (cosθ (t,x),sinθ (t,x),0) we have that n 0 is a solution of the heat-flow: with initial data t n 0 = n 0 +n 0 n 0 2 (.22) n 0 (0,x) = n 0 (x) := (cosθ 0,sinθ 0,0) (.23) while n is a solution of the wave-map flow: with initial datas: ( tt n +n t n 2 )+ t n = n +n n 2 (.24) n (0,x) = n 0 (x) := (cosθ 0 (x),sinθ 0 (x),0), (.25) t n (0,x) = n (x) = ( sinθ 0 (x),cosθ 0 (x),0)θ (x) (.26) Taking θ = θ 0 the claimed optimality of the rate of converge is shown in 2]. A rigorous justification of the formal expansion (.4) in the context of classical solutions is provided in this paper. For the original wave map (.8) with a damping the energy bounds of d uniform in small > 0 do not seem available. By taking the expansion (.4) of the solutions d to the system (.8) with the initial conditions (.9), one can yield a remainder system (.5)-(.6) of d R. Although the remainder system (.5) with the initial data conditions (.6) is still nonlinear and singular (with singular terms of the type ), it has weaker nonlinearities than the original system (.8). More precisely, after using the expansion (.4), the nonlinear term ( t d 2 + d 2 )d is replaced by either linear terms (the unknown d is superseded by the known d 0 ) or a nonlinear term with the same form but with some higher order powers of in front. So, by standard energy estimates, we can get uniform bounds in small > 0 of the remainder system (.5)-(.6). The organization of this paper is as follows: in next section, we give the formal analysis for the asymptotic behavior of the wave map (.8) with a damping and initial conditions (.9) as the inertia density 0 by constructing the initial layer d I 0 (t,x) to deal with the compatibility of the original initial conditions(.9) and the initial condition of the limit system (.). In Section 3, we estimate the uniform energy bounds on small > 0 of the remainder

6 6 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU system (.5)-(.6). Finally, based on the uniform energy estimates in the previous section, Theorem. of this paper is proved in Section Formal Analysis In this section we present the formal analysis of the limit 0 for the damped wave map (.8) with the initial conditions (.9). Out of the equation (.8) we note that the formal limit, obtaining by setting = 0, is the heat flow system (.) for functions taking values into S 2. We can then naturally take the ansatz d (t,x) = d 0 (t,x)+ ˆd R (t,x), (2.) where d 0 (t,x) is a solution of the heat flow system (.) and ˆd R (t,x) satisfies a hyperbolic system, formally similar to (.8) but without the geometric constraint of taking values into S 2. If the ansatz (2.) were reasonable then ˆd R (t,x) = O(α ) in some sense for some α > 0 as > 0 is small enough. However, by the second initial condition in (.9) and the heat flow system (.), we know that tˆd R (0,x) = t d (0,x) t d 0 (0,x) = d in (x) t d 0 (0,x), which will not go to 0 as 0 for arbitrarily given vectors d in (x) and d in (x). As a consequence, ˆd R (t,x) O(α ) uniformly in time for any α > 0, and then the ansatz (2.) is not satisfactory. Therefore, in order to compensate the effect of the initial data, we need to introduce a correction term d I t 0(,x ) for some β > 0 to be determined, called initial layer, such that β the second initial condition in (.9) is satisfied by d 0 (t,x)+d I t 0(,x ) without disturbing too β much the first initial condition in (.9), namely d I 0( 0 β,x ) = O( α ) for some α > 0 as 0. Thus we take the alternative ansatz d (t,x) = d 0 (t,x)+d I t 0(,x ) + d β R (t,x), (2.2) where the power in front of the remainder term is motivated by the scaling we chose. This measures the rate of convergence, as it will be shown in the proof of Theorem.. Recalling that d 0 is a solution of the heat flow (.), we plug (2.2) into the system (.8) and obtain: tt (d 0 + d R)+ t d R + ( tt d I 0 + t d I 0) (d I 0 + d R) = t (d 0 +d I 0 + d R) 2 (d 0 +d I 0 + d R) + (d 0 +d I 0 + d R ) 2 (d I 0 + d R )+ 2 d 0 (d I 0 + d R )+ (di 0 + d R ) 2] d 0. Then we construct the initial layer in order to cancel certain time-dependent terms in the previous equations and to accommodate the discrepancy in the initial data, namely we take d I 0 satisfying the x-dependent ODE and the initial-data condition: { tt d I t 0(,x ) + β t d I ( t 0,x ) = 0, β t d I 0 0(,x ) = d in (x) β t d 0 (0,x). Furthermore, since d I 0 is an initial layer, the following condition at infinity is required: d I 0 (,x) = lim y di 0 (y,x) = 0. (2.3)

7 WAVE MAP LIMIT 7 By solving the ODE system with the given boundary conditions we have d I t 0(,x ) = ( din (x) β t d 0 (0,x) ) exp ( t ). (2.4) We remark that the initial layer d I t 0(,x ) in (2.4) is, in fact, independent of β > 0 and β d I 0 0(,x ) = ( din (x) β t d 0 (0,x) ) 0 for any given d in (x) and d in (x) as 0. Consequently, the ansatz (2.2) is reasonable. Without loss of generality, we take β = in the ansatz (2.2). Thus, by substituting (.) and (.4) into the system (.8), we derive the equation satisfied by the remainder d R (t,x) as follows: 3 2 tt d R + tt d 0 + t d R d R + D(x)exp ( t ) ( = t d0 D(x)exp ( ) t ) + d 2( R d0 D(x)exp ( t ) ) + d R + ( d 0 D(x)exp ( t ) ) + d 2( ( ) ) R D(x)exp t + d R + 2 d 0 ( D(x)exp( t )+ d R) + ( D(x)exp( t )+ d R) 2] d 0, which, after multiplication by 3 2 is the equation (.5) we used before. 3. Uniform Energy Estimates In this section, we will provide, by energy methods, bounds that are uniform with respect to small inertia constant > 0, for the remainder system (.5)-(.6). By Proposition., the d 0, which obeys the heat flow (.) into the unit sphere S 2 is regarded as a known quantity in the remainder system (.5)-(.6). To conveniently state our results, we need to introduce the following energy functionals: E (t) = t d R 2 H 2 + ( ) d R 2 H d R 2 H 2 + t d R +d R 2 H 2, F (t) = ( 2) t d R 2 H d R 2 H 2. The following lemma provides the claimed uniform energy estimates : Lemma 3.. Let d in Ḣ7 (R 3 ;S 2 ) and 0,T] be the interval of existence of the the solution of the heat flow with initial data d in, as provided in Proposition.. For d in H 5 assume that there exists a 0 0 ( d in H 6, d in H 5,T) (0, 2 ) such that for all (0, 0 ) we that d R L (0,T);H 3 ) and t d L (0,T);H 2 ) is a solution to the remainder system (.5)-(.6). Then there exists a positive constant C = C( d in H 6, d in H 5,T) > 0 such that the inequality holds for all (0, 2 ) and t 0,T). d dt E (t)+3f (t) C +E (t) ] +E (t) ] (3.) We remark that the condition 0 < < 2 guarantees the relation 2 < <, which makes the energy functionals E (t) and F (t) non-negative. Since our goal is to rigorously analyze the asymptotic behavior as 0 for the wave map (.8)-(.9), the condition 0 < < 2 is sufficient.

8 8 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU Proof. For the convenience of notations, we rewrite the singular terms of the remainder system (.5) as S(d R ) = ( tt d 0 + D(x)exp( t )+ td 0 +D(x)exp( t ) 2 d 0 ) + d 0 2 D(x)exp( t )+2 d 0 D(x)exp( t )d S 0 + ( 2( d0 d ) } R )d R + d R 2 d 0 (3.2) + d0 2 d R +2( d 0 d R )d 0] } S3 S +S 2 +S 3 and the regular terms as R(d R ) = t d 0 +D(x)exp( t ) 2 D(x)exp( t )+ D(x) 2 exp( 2t )d 0 2 D(x) 3 2 exp( 3t )D(x)+2 ( d 0 D(x))exp( 2t )D(x) 2 ( t d 0 +D(x)exp( t )) t d R d 0 t d 0 +D(x)exp( t ) 2 d R +2 ( t d 0 +D(x)exp( t )) t d R D(x)exp( t )+ D 2 exp( 2t )d R 2( d 0 D(x))exp( t )d R 2( d 0 d R )D(x)exp( t ) 2( D(x) d R )exp( t )d 0 +2( D(x) d R )exp( 2t )D(x) + 3 2D(x)exp( t ) td R 2 2 ( t d 0 +D(x)exp( t )) t d R d R t d R 2 d 0 2 ( D(x) d R )exp( t )d R d R 2 D(x)exp( t ) t d R 2 d R + d R 2 d R S 2 } } R R 2 R 3 } R4 R +R 2 +R 3 +R 4. (3.3) Step. L 2 -estimates. Multiplying the remainder equation (.5) by t d R, integrating over R 3 and by parts, we obtain the following equation: ( d 2 dt t d R 2 L + ) 2 d R 2 L + 2 td R 2 L = S(d 2 R ), td R + R(d R ), td R. (3.4) (I) Estimates for the singular terms S(d R ), td R : For estimating S, t d R, we use the Hölder inequality, the Sobolev embedding theorems, the facts that exp( t ) and d 0 = to obtain: tt d 0 + D(x)exp( t ), td R ( tt d 0 L t L 2 + D(x) x L 2) td R L 2 and d 0 2 D(x)exp( t )+2 d 0 D(x)exp( t )d 0, t d R ( d 0 2 L D(x) L 2 + d 0 L D(x) L 2) td R L 2 Similarly: ( d 0 2 L + d 0 t H2 x L t Hx) 2 D(x) H td R L 2. t d 0 +D(x)exp( t ) 2 d 0, t d R t d 0 +D(x)exp( t ) 2 L td 4 R L 2 ( t d 0 2 ) L t L4 + D(x) 2 x L td 4 R L 2 ( t d 0 2 ) L t H + D(x) 2 x H td R L 2. Summarizing, we estimate S, t d R as follows: S, t d R tt d 0 L t L 2 + x td 0 2 L t H + D(x) 2 x H ] +( d 0 2 L +) D(x) td t H2 H x 2 R L 2. x (3.5)

9 WAVE MAP LIMIT 9 It is easy to derive the estimates of S 2, t d R and S 3, t d R as follows: S 2, t d R = 2( d 0 d R )d R + d R 2 d 0, t d R ( d 0 L d R L d R L 2 + d R 2 L ) td 4 R L 2 (3.6) ( d 0 L t Hx 2 +) d R 2 H td 2 R L 2 and S 3, t d R = d 0 2 d R +2( d 0 d R )d 0, t d R ( d 0 2 L dr L 2 + d 0 L d R L 2) td R L 2 ( d 0 L t H 2 x + d 0 2 L t H2 x)( d R L 2 + d R L 2) td R L 2. (3.7) Hence we have the estimate of S(d R ), td R by combining the inequalities (3.5), (3.6) and (3.7): where the constant S(d R), t d R C ( + d R 2 H 2 + dr L 2 + d R L 2) td R L 2, (3.8) C = C tt d 0 L t L 2 x + td 0 2 L t H x + D(x) 2 H +( d 0 2 L t H2 x +)( D(x) H 2 +)] > 0 for some computable positive constant C > 0. (II) Estimates for the regular terms R(d R ), td R : We have divided theregular terms R(d R ) into fourparts, which we will estimate separately. The estimate of R, t d R. By thehölderinequality andthesobolevembeddingtheorems we have t d 0 +D(x)exp( t ) 2 D(x)exp( t ), td R t d 0 +D(x)exp( t ) 2 L D(x) L 2 td R L 2 ( t d 0 2 L t H2 x + D(x) 2 H 2 ) D(x) L 2 t d R L 2 and ( d0 D(x))exp( 2t )D(x), td R d 0 L D(x) L D(x) L 2 t d R L 2 d 0 L t H 2 x D(x) 2 H 2 t d R L 2, where we have used the fact that exp( t ). The other two terms in R, t d R are similarly estimated, as follows: D(x) 2 exp( 2t )d 0, t d R D(x) 2 L 4 t d R L 2 D(x) 2 H 2 t d R L 2, 3 2 D(x) 2 exp( 3t )D(x), td R 3 2 D(x) L D(x) 2 L 4 t d R L D(x) 3 H 2 t d R L 2. For any small enough such that (0, 2 ], we have So from the above inequalities, we obtain the following estimate: R, t d R ( t d 0 2 L t H2 x + d 0 2 L t H2 x + D(x) 2 H 2) D(x) H 2 + D(x) 2 H 2 ] t d R L 2. (3.9)

10 0 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU The estimate of R 2, t d R. We have the following estimates for the first three terms in R 2, t d R : t d 0 +D(x)exp( t ) 2 d R, td R td 0 +D(x)exp( t ) 2 L d R L 2 td R L 2 ( t d 0 2 ) L t H2 + D(x) 2 x H d 2 R L 2 t d R L 2, ( t d 0 +D(x)exp( t )) t d Rd 0, t d R ( t d 0 L + D(x) L ) t d R 2 L 2 ( t d 0 L t Hx 2 + D(x) ) H 2 t d R 2 L 2 and ( t d 0 +D(x)exp( t )) D(x)exp( t ) td R, t d R t d 0 +D(x)exp( t ) L D(x) L td R 2 L 2 ( t d 0 L t H 2 x + D(x) H 2 ) D(x) H 2 t d R 2 L 2, where we have used the Hölder inequality, the Sobolev embedding theorems and the fact that d 0 =. As for the following three terms, one can easily obtain: D(x) 2 exp( 2t )d R, td R D(x) 2 L d R L 2 td R L 2 D(x) 2 H 3 d R L 2 td R L 2, ( d 0 D(x))exp( t )d R, t d R d 0 L D(x) L d R L 2 t d R L 2 d 0 L t H 2 x D(x) H 3 d R L 2 t d R L 2, ( d 0 d R )D(x)exp( t ), td R d 0 L D(x) L d R L 2 td R L 2 d 0 L t H 2 x D(x) H 2 d R L 2 td R L 2, where we have used the Hölder inequality, the Sobolev embedding theorems and the bound exp( t ). Similarly as before, we estimate the last two terms, as follows: ( D(x) d R)exp( t )d 0, t d R D(x) L d R L 2 t d R L 2 D(x) H 3 d R L 2 t d R L 2, ( D(x) d R )exp( 2t )D(x), td R D(x) L D(x) L d R L 2 td R L 2 D(x) 2 H 3 d R L 2 t d R L 2. Combining the above estimates and using that (0, 2 ], we get R 2, t d R ( + t d 0 L t H 2 x + d 0 L t H 2 x + D(x) H 3 )( t d 0 L t H 2 x + D(x) H 3 ) ( d R L 2 + d R L 2 + td R L 2 ) t d R L 2. (3.0) The estimate of R 3, t d R.One can easily derive the following estimates ( t d 0 +D(x)exp( t )) td R d R, td R ( t d 0 L t H2 x + D(x) H 2 ) d R H 2 t d R 2 L 2, ( D(x) d R )exp( t )d R, td R D(x) H 3 d R H 2 d R L 2 td R L 2,

11 WAVE MAP LIMIT by using the Hölder inequality and the Sobolev embedding theorems. Recalling that d 0 = and using the Sobolev embeddings H L 4 and H 2 L we get 3 2 D(x)exp( t ) td R 2, t d R 3 2 D(x) L t d R 2 L 4 t d R L D(x) H 2 t d R 2 H t d R L 2, t d R 2 d 0, t d R t d R 2 L 4 t d R L 2 t d R 2 H t d R L 2, d R 2 D(x)exp( t ), td R D(x) L d R 2 L 4 t d R L 2 D(x) H 2 d R 2 H t d R L 2, The above estimates immediately give the bound on R 3, t d R : R 3, t d R ( + t d 0 L t H 2 x + D(x) H 3 )( t d R 2 H + d R 2 H + d R H 2 d R L 2 + d R H 2 td R L 2 ) t d R L 2 ( + t d 0 L t H 2 x + D(x) H 3 )( t d R 2 H + d R 2 H 2 ) t d R L 2. (3.) The estimate of R 4, t d R. The first term in R 4, t d R can be bounded as t d R 2 d R, td R d R L td R 2 L 4 t d R L 2 d R H 2 t d R 2 H t d R L 2 by using the Hölder inequality and the Sobolev embedding theorems. The other term can be bounded in a similar way: d R 2 d R, td R d R L d R 2 L 4 t d R L 2 d R H 2 d R 2 H t d R L 2. Hence we obtain the estimate of R 4, t d R as follows: R 4, t d R ( t d R 2 H + d R 2 H ) d R H 2 t d R L 2. (3.2) Summing up the inequalities (3.9), (3.0), (3.) and (3.2), we get where R(d R ), td R C 2 + d R H + t d R L 2 + t d R 2 H + d R 2 H 2 + d R H 2( t d R 2 H + d R 2 H ) ] t d R L 2, C 2 = C ( + t d 0 2 L t H2 x + d 0 2 L t H2 x + D(x) 2 H 3 )( + t d 0 L t H 2 x + D(x) H 3 ) > 0, and C is a positive computable constant. Therefore, plugging the estimates (3.8) and (3.3) into the equality (3.4), we have ( d 2 dt t d R 2 L + ) 2 d R 2 L + 2 td R 2 L 2 { (+ d C R 2 H + d 2 R L 2 + d R L 2) td R + d L 2 + R H + t d R L 2 + t d R 2 H + d R 2 H + d R 2 H 2( t d R 2 H + d R 2 ) ] } H t d R L 2, where the constant C = C(+ D(x) 3 H 3 + tt d 0 L (0,T;L 2 ) + t d 0 3 L (0,T;H 2 ) + d 0 4 L (0,T;H 2 ) ) > 0 for some computable positive constant C. (3.3) (3.4) (III) Estimates of the norm d R L 2: Observing that the norm d R L2 appearing on the right hand side of (3.4) is not yet controlled, we need additional work to estimate d R L2. In order to do this it is natural to

12 2 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU multiply the equation of the reminder term (.5) by d R, integrate over R3 and by parts, and use the identity: tt d R,d R = d dt td R,d R td R 2 L 2 = d 2 dt( t d R +d R 2 L 2 t d R 2 L d ) 2 R 2 L 2 t d R 2 L, 2 in order to get d 2 dt t d R +d R 2 L + ( 2 ) d R 2 L 2 t d R 2 L ] 2 t d R 2 L + 2 d R 2 L 2 (3.5) = S(d R ),d R + R(d R ),d R. Using the estimates (3.8) and (3.3) previously derived for bounding the terms S(d R ), td R and R(d R ), td R, we can analogously estimate the terms S(d R ),d R and R(d R ),d R as follows: and S(d R ),d R C ( + d R 2 H + dr 2 L 2 + d R d L 2) R L 2, (3.6) R(d R),d R C 2 + d R H + t d R L 2 + t d R 2 H + d R 2 H 2 + d R H 2( t d R 2 H + d R 2 H ) ] d R L 2. So plugging (3.6) and (3.7) into (3.5) we obtain: d 2 dt t d R +d R 2 L + ( 2 ) d R 2 L 2 t d R 2 L ] 2 t d R 2 L + 2 d R 2 L 2 { (+ d C R 2 H + d 2 R L 2 + d R d L 2) R + d L 2 + R H + t d R L 2 + t d R 2 H + d R 2 H + d R 2 H 2( t d R 2 H + d R 2 ) ] } H d R L 2. (3.7) (3.8) Multiplying the inequality (3.8) by 2 and then adding it to the inequality (3.4), we get the L 2 -energy estimate: d 4 dt t d R 2 L + 2 t d R +d R 2 L + ( 2 ) d R 2 L d R 2 L 2 ]+( 2 ) td R 2 L d R 2 L 2 C { (+ d R 2 H 2 + dr L 2 + d R L 2)( td R L 2 + dr L 2) (3.9) where the constant + + d R H + t d R L 2 + t d R 2 H + d R 2 H 2 + d R H 2 ( t d R 2 H + d R 2 H ) ] ( t d R L 2 + d R L 2 ) }, C = C(+ D(x) 3 H 3 + tt d 0 L (0,T;L 2 ) + t d 0 3 L (0,T;H 2 ) + d 0 4 L (0,T;H 2 ) ) > 0 for some computable positive constant C. Step 2. Higher order estimates. In order to use the inequality (3.9) we also need a higher order estimate. To obtain this we take k (k =,2) in the equation (.5), we multiply it by k t d R, integrate over R3 and by parts, thus obtaining the following equality ( d 2 dt t k d R 2 L + 2 k+ d ) R 2 L + 2 t k d R 2 L 2 = k S(d R), t k d R + k R(d R), t k d (3.20) R (I) Estimates of the singular terms k S(d R ), t k d R :

13 WAVE MAP LIMIT 3 The singular terms can be divided into three parts: k S i, t k d R (i =,2,3) which we estimate separately. For the term k S, t k d R, by using the Hölder inequality and the Sobolev embedding theorems, we obtain: k( tt d 0 + D(x)exp( t )), t k d R ( tt d 0 L t Hx 2 + D(x) H 4) td R H 2 and k ( t d 0 +D(x)exp( t ) 2 d 0 ), t k d R i ( t d 0 +D(x)) j ( t d 0 +D(x)) e d 0, t k d R e + i+j=k i ( t d 0 +D(x)) j ( t d 0 +D(x))d 0, t k d R (+ d 0 L t Hx 2)( td 0 2 L t H2 + D(x) 2 x H ) td 2 R H 2. Similarly as for estimating k ( t d 0 +D(x)exp( t ) 2 d 0 ), t k d R, we can easily get the following estimates: k( d 0 2 D(x)exp( t )), t k d R i+ d 0 j+ d 0 e D(x), t k d R d 0 2 L D(x) t H2 H x 3 td R H 2 and k( d 0 D(x)exp( t )d 0), t k d R i+ d 0 j+ D(x) e d 0, t k d R e + i+j=k i+ d 0 j+ D(x)d 0, t k d R (+ d 0 L t Hx 2) d 0 L t Hx 3 D(x) H 3 td R H 2. Thus we have the estimate of k S, t k d R as follows: k S, t k d R {(+ d 0 L t H 2 x )( td 0 2 L t H3 x + d 0 2 L t H3 x + D(x) 2 H 3 ) } + tt d 0 L t Hx 2 + D(x) td H 4 R H 2. (3.2) For the term k S 2, t k d R, we can also use the Hölder inequality and Soblev embedding theorems to get k( ( d 0 d R R) )d, t k d R i+ d 0 j+ d R e d R, t k d R e + i+j=k i+ d 0 j+ d R d R, t k d R d 0 L t Hx 4 d R H 2 d R H 2 td R H 2,

14 4 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU and Summarizing, we obtain k( d ) R 2 d 0, t k d R i+ d R j+ d R e d 0, t k d R e + i+j=k i+ d R j+ d Rd 0, t k d R (+ d 0 L t Hx 3) d R 2 H td 2 R H 2. k S 2, t k d R (+ d 0 L t Hx 4)( d R H 2 + d R H 2) d R H 2 td R H 2. (3.22) For the estimate of k S 3, t k d R, we get the estimate of the first term by using again the Hölder inequality and the Sobolev embedding theorems: k ( d 0 2 d R), t k d R i+ d 0 j+ d 0 e d R, t k d R = d 0 2 L t H4 x dr H 2 td R H 2. Recalling that d 0 =, one can easily estimate the second term k( ( d 0 d ) R)d 0, t k d R i+ d 0 j+ d R e d 0, t k d R e + i+j=k Thus by the above two estimates we have i+ d 0 j+ d Rd 0, t k d R (+ d 0 L t Hx 3) d 0 L t Hx 4 d R H 2 td R H 2. k S 3, t k d R (+ d 0 L t H 4 x ) d 0 L t H 4 x( d R H 2 + d R H 2) td R H 2. (3.23) Then the inequalities (3.2), (3.22) and (3.23) give the following estimate k S(d R), t k d R C k ( + d R 2 H 2 + d R 2 H 2 + dr H 2 + d R H 2) td R H 2, (3.24) where the positive constant C k is C k = C{ (+ d0 L t H4 x)( + t d 0 2 L t H3 x + d 0 2 L t H4 x + D(x) 2 H 3 ) + tt d 0 L t H 2 x + D(x) H 4 } for some computable positive constant C. (II) Estimates of the regular terms k R(d R ), t k d R : Finally, we turn to estimating the regular terms k R(d R ), t k d R, which are divided into four parts: k R i, t k d R (i =,2,3). For the terms k R, t k d R, by the Hölder inequality and the Sobolev embedding theorems we have: k ( t d 0 +D(x)exp( t ) 2 D(x)exp( t )), t k d R i ( t d 0 +D(x)) j ( t d 0 +D(x)) e D(x), t k d R D(x) H 3( D(x) 2 H 2 + t d 0 2 L t H2 x ) td R H 2,

15 and WAVE MAP LIMIT 5 k ( D(x) 2 exp( 2t )d 0), t k d R i+ D(x) j+ D(x) e d 0, t k d 0 e + i+j=k i+ D(x) j+ D(x)d 0, t k d R (+ d 0 L t H 2 x ) D(x) 2 H 3 t d R H 2. We can estimate the following two terms in a similar way, hence we get the following inequalities: 3 2 k ( D(x) 2 exp( 3t )D(x)), t k d R 3 2 i+ D(x) j+ D(x) e D(x), t k d R and 3 2 D(x) 3 H 3 t d R H 2 k ( d 0 D(x))exp( 2t )D(x)), t k d R i+ d 0 j+ D(x) e D(x), t k d R d 0 L t H 4 x D(x) 2 H 3 t d R H 2. So we have the estimate of k R, t k d R as follows: k R, t k d R (+ d 0 L t H 4 x + D(x) H 3) ( t d 0 2 L t H2 x + D(x) 2 H 3 ) t d R H 2. (3.25) For the terms k R 2, t k d R, by using yet again the Hölder inequality and the Sobolev embedding theorems, we have k( t d 0 +D(x)exp( t )) t d RD(x)exp( t ), t k d R ( i t d 0 + i D(x)exp( t )) j t d R e D(x)exp( t ), t k d R e + i+j=k ( i t d 0 + i D(x)exp( t )) j t d RD(x)exp( t ), t k d R ( t d 0 L t H 4 x + D(x) H 4) D(x) H 4 td R 2 H 2, Similarly as for estimating the term k( t d 0 +D(x)exp( t )) t d R D(x)exp( t ), t k d R we can easily obtain the following estimates: k( t d 0 +D(x)exp( t ) 2 d R), t k d R ( t d 0 2 L t H4 x + D(x) 2 H 4 ) d R H 2 t d R H 2, k( d 0 d R D(x)exp( t )), t k d R ( d 0 2 L t H4 x + D(x) 2 H 4 ) d R H 2 td R H 2, k( d 0 D(x)exp( t )d R), t k d R ( d 0 2 L t H4 x + D(x) 2 H 4 ) d R H 2 td R H 2. Observingthestructureoftheterms k( D(x) 2 exp( 2t R) )d, t k d R and k( D(x) d R exp( 2t )D(x)), t k d R one can similarly estimate the following terms: k( D(x) 2 exp( 2t )d R), t k d R D(x) 2 H 5 d R H 2 t d R H 2, k( D(x) d R exp( 2t )D(x)), t k d R D(x) 2 H 5 d R H 2 td R H 2.

16 6 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU Furthermore we get: k ( t d 0 +D(x)exp( t )) t d Rd 0 ], t k d R i ( t d 0 +D(x)) j t d R e d 0, t k d R e + i+j=k i ( t d 0 +D(x)) j t d Rd 0, t k d R (+ d 0 L t H 3 x )( td 0 L t H 4 x + D(x) H 4) td R 2 H 2 Similarly as for estimating k ( t d 0 +D(x)exp( t )) t d R d 0], t k d R, one can also get: k( D(x) d R exp( t )d 0), t k d R (+ d 0 L t H3 x ) D(x) H 5 d R H 2 td R H 2. Thus we have the following estimate of k R 2, t k d R : k R 2, t k d R ( td 0 L t H 4 x + D(x) H 5 + td 0 2 L t H4 x + D(x) 2 H 5 + d 0 2 L t H4 x ) ( d R H 2 + d R H 2 + td R H 2) td R H 2. (3.26) For the terms k R 3, t k d R, one can use an estimate similar to the one for the term k( t d 0 D(x)exp( t ) td R), t k d R to get: 3 2 k ( D(x)exp( t ) td R 2), t k d R 3 2 D(x) H 4 t d R 3 H, 2 k ( D(x) d ) R exp( t )d R, t k d R D(x) H 5( d R 2 H + d 2 R 2 H ) 2 t d R H 2, k ( d R 2 D(x)exp( t )), t k d R D(x) H 4 d R 2 H 2 t d R H 2. Reasoning analogously as in estimating k ( t d 0 +D(x)exp( t )) t d R d R ], t k d R we have k ( t d 0 +D(x)exp( t )) t d R d R ], t k d R e + i+j=k i ( t d 0 +D(x)) j t d R e d R, t k d R i ( t d 0 +D(x)) j t d R d R, t k d R ( t d 0 L t H 4 x + D(x) H 4)( d R H 2 + d R H 2) td R 2 H 2 and furthermore, using that d 0 =, it is easy to obtain k ( t d R 2 d 0 ), t k d R (+ d 0 L t H3 x ) td R 3 H 2. Summarizing, we get the estimate of k R 3, t k d R as follows: k R 3, t k d R (+ t d 0 L t H 4 x + d 0 L t H 3 x + D(x) H 5) (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 ) t d R H 2. (3.27) For the terms k R 4, t k d R, we get, by the Hölder inequality and Sobolev embedding theorems: k ( t d R 2 d R ), t k d R = i t d R j t d R e d R, t k d R e + i t d R j t d R d R, t k d R i+j=k

17 WAVE MAP LIMIT 7 ( d R H 2 + d R H 2) td R 3 H 2 and then, similarly: k ( d R 2 d R), t k d R ( d R H 2 + d R H 2) d R 2 H 2 t d R H 2. So we obtain the estimate of k R 4, t k d R as follows: k R 4, t k d R ( d R H 2 + d R H 2)( d R 2 H 2 + t d R 2 H 2 ) t d R H 2. (3.28) Then the inequalities (3.25), (3.26), (3.27) and (3.28) give the estimate of the regular terms k R(d R ), t k d R as follows: k R(d R), t k d R C k2 ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 ) where the constant C k2 is + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 )+ t d R H 2 ] t d R H 2, (3.29) C k2 = C ( + t d 0 L t H 4 x + d 0 L t H 4 x + D(x) H 5 )( + t d 0 2 L t H4 x + D(x) 2 H 5 + d 0 L t H 4 x) > 0, and C is a computable positive constant. Therefore, by substituting the inequalities (3.24) and (3.29) into (3.20) one has d 2 dt ( t k d R 2 L 2 + k+ d R 2 L 2 ) + t k d R 2 L 2 C k { (+ d R 2 H 2 + d R 2 H 2 + dr H 2 + d R H 2) td R H 2 + ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 ) + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 )+ t d R H 2 ] t d R H 2 }, (3.30) where the positive constant C k is C k = C ( + tt d 0 L (0,T;H 2 ) + t d 0 3 L (0,T;H 4 ) + d 0 3 L (0,T;H 4 ) + D(x) 3 H 5 ) > 0 and C > 0 is a computable constant. (III) For the estimate of k d R L2 (k =,2): Applying k (k =,2) to the remainder equation (.5), multiplying by k d R, integrating over R 3 and by parts, we have d 2 dt k t d R + k d R 2 L + ( 2 ) k d R 2 L 2 t k d R 2 L ] 2 t k d R 2 L + 2 k+ d R 2 L 2 = k S(d R ), k d R + k R(d R ), k d R. (3.3) Similarly as in the estimates of the terms k S(d R ), t k d R and k R(d R ), t k d R in the inequalities(3.24)and(3.29),respectively, wecananalogouslyestimatetheterms k S(d R ), k d R and k R(d R ), k d R as follows: and k S(d R ), t k d R C k( + d R 2 H + d 2 R 2 H + dr 2 H 2 + d R d H 2) R H 2, (3.32) k R(d R ), t k d R C k2 ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 ) + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 )+ t d R H 2 ] d R H 2. (3.33)

18 8 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU By plugging the inequalities (3.32) and (3.33) into the equality (3.3), we get the following estimate: d 2 dt k t d R + k d R 2 L + ( 2 ) k d R 2 L 2 t k d R 2 L 2 ] t k d R 2 L 2 + k+ d R 2 L 2 { (+ d C k R 2 H + d 2 R 2 H + dr 2 H 2 + d R d H 2) R H 2 (3.34) + ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 ) + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 )+ t d R H 2] d R H 2 Multiplying the inequality (3.34) by 2 and adding it to the inequality (3.30), we obtain the higher order estimate: d 4 dt t k d R 2 L + ( ] 2 ) k d R 2 L k+ d R 2 L + k 2 t d R + k d R 2 L 2 + ( 2) t k d R 2 L k+ d R 2 L 2 3 (+ d k{ 2 C R 2 H + d R 2 2 H + dr 2 H 2 + d R d H 2)( R H 2 + td R H 2 + ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 )+ t d R H 2 + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 ) ]( d R H 2 + t d R H 2) }. }. ) (3.35) Therefore, combining the L 2 -estimate (3.9) and the kth-order estimate (3.35) for k =,2, we obtain: d 4 dt t d R 2 H + ( ] 2 ) d R 2 H d R 2 H + 2 t d R +d R 2 H 2 + ( 2) t d R 2 H d R 2 H 2 C { (+ d k R 2 H + d R 2 2 H + d 2 R H 2 + d R d H 2)( R H 2 + td R H 2 + ( d R H 2 + d R H 2)(+ d R 2 H 2 + t d R 2 H 2 )+ t d R H 2 + (+ t d R 2 H 2 + d R 2 H 2 + d R 2 H 2 ) ]( d R H 2 + t d R H 2) }, where the positive constant C k is C k = C ( + tt d 0 L (0,T;H 2 ) + t d 0 3 L (0,T;H 4 ) + d 0 3 L (0,T;H 4 ) + D(x) 3 H 5 ) > 0 ) (3.36) and C > 0 is a computable constant. Then, by the definition of the energy functionals E (t) and F (t), and the condition 0 < < 2, the H2 -estimate (3.36) implies that d dt E (t)+4f (t) C { } E 2 (t)+e (t)+ 2E 3 2 (t)+ 3 2E 2 (t)++e 2 (t)+e (t)]f 2 (t), where C > 0, which immediately implies the claimed inequality (3.) by using Young s inequality. Consequently, the proof of Lemma 3. is completed. 4. The Proof of Theorem. In this section we will provide the proof of Theorem., by using the uniform energy bounds (3.) in Section 3. Before doing this, we note that for any fixed inertia constant > 0 the well-posedness of the remainder system (.5)-(.6) can be stated as follows: Proposition 4.. Given d in : R 3 S 2 and d in : R 3 R 3 satisfying d in H 6, din H 4 with d in d in 0, we define D(x) d in (x) d in (x) d in (x) 2 d in (x) and denote M = D 2 H 2 +2 D 2 H 2 <.

19 WAVE MAP LIMIT 9 Then, for any fixed (0, 2 ), there exists a time T = min{t, C ln( +M (+M)) } > 0, where T, C > 0 are provided in Proposition., Lemma 3., respectively, such that the remainder equation (.5) with the initial conditions (.6) admits a unique solution d R C(0,T );H 3 ) and t d R C(0,T );H 2 ). Moreover, the solution d R satisfies the inequality for all t 0,T ). t d R(t) 2 H 2 + d R(t) 2 H 3 2Me Ct +M (+M)e Ct (4.) Proof. We employ a mollifier argument to prove this proposition. For any fixed > 0 we can directly construct a system approximating (.5)-(.6) as follows: t wη = J ηwη +J η d R,η +J ηs(j η d R,η )+J ηr(j η d R,η ), d R,η where the mollifier operator J η is defined as t d R,η = w η, t=0 = J η D(x), wη t=0 = 0, J η f = F ( ξ F(f)(ξ) ), η where the symbol F denotes the standard Fourier transform operator and F is the inverse Fourier transform operator. By ODE theory in Hilbert spaces one can prove the existence and uniqueness of the approximate system (4.2) on the maximal time interval 0,T η). Then by the fact J 2 η = J η and the uniqueness of (4.2) we know that J η d R,η = d R,η and J ηw η = w η. Thus by the analogous energy estimate shown in Lemma 3. applied to the approximate system (4.2), one can obtain the following energy inequality for d R,η and w η d dt E,η(t)+3F,η (t) C +E,η (t) ] +E,η (t) ] (4.3) for all t 0,Tη ), where the positive constant C > 0 is independent of and η, and the energy functionals E,η (t), F,η (t) are of the same forms as E (t), F (t) defined in Section 3 (replacing d R by d R,η ), respectively. Since d R,η satisfies the initial conditions td R,η (0,x) = 0 and d R,η (0,x) = J η D(x), we know that for (0, 2 ) E,η (0) = t d R,η (0, ) 2 H 2 +( ) d R,η (0, ) 2 H 2 (4.2) + 2 d R,η (0, ) 2 H 2 + t d R,η (0, )+d R,η (0, ) 2 H 2 =( ) J η D 2 H J η D 2 H 2 + J η D 2 H 2 (4.4) ( ) D 2 H 2 +2 D 2 H 2 + D 2 H 2 = D 2 H 2 +2 D 2 H 2 = M <. Then, one can solve the ODE inequality (4.3) with the initial condition (4.4), obtaining that +E,η (t) +E,η (t) +E,η(0) +E,η (0) ec( )t +M +M ect holds for all t 0,Tη). Consequently, for all t 0,min{Tη,T, C ln( ) ) +M (+M) } E,η (t) we know that (+M)e Ct. (4.5) +M (+M)eCt Notice that the continuity of E,η (t) and the maximality of T η > 0 imply that T η C ln( +M (+M)) > 0. Hence the inequality (4.5) holds for all t 0,T ) uniformly in η > 0, where T = min{t, C ln( +M (+M)) } > 0.

20 20 N. JIANG, Y.-L. LUO, S. J. TANG, AND A. ZARNESCU Finally, we can finish the proof of this proposition by standard compactness methods and taking the limit as η 0. The uniqueness issue can be reduced to the uniqueness of the damped wave map system (.8) which can be obtained by methods analogous to those in the book of Shatah and Struwe 5]. For convenience, we omit the details of the proof. The Proof of Theorem.. Now, based on the energy estimate (4.) in Proposition 4., we verify Theorem.. We observe that the function f() := ( +M ) C ln (+M) is strictly decreasing in (0, 2 ) and limf() = +. Consequently, we can choose ց0 such that for any (0, 0 ) 0 = min { 2, (+M)e CT M } (0, 2 ) f() = ( +M ) C ln > T. (+M) As a result, for the number T determined in Proposition 4., we have that T T for all (0, 0 ). Then, the inequality (4.) in Proposition 4. implies that t d R 2 L (0,T;H 2 ) + d R 2 L (0,T;H 3 ) 2Me CT + 0 M 0 (+M)e CT := C 0 <, and the proof of Theorem. is completed. Acknowledgement The project of this paper was initialized when Ning Jiang visited Arghir Zarnescu at Basque Center for Applied Mathematics (BCAM) in March 207. They appreciate the hospitality of BCAM. The activity of Ning Jiang on this work was supported by Chinese NSF grant 478. The activity of Arghir Zarnescu on this work was partially supported by a grant of the Romanian National Authority for Scientific Research and Innovation, CNCS-UEFISCDI, project number PN-II-RU-TE ; by the Project of the Spanish Ministry of Economy and Competitiveness with reference MTM P; by the Basque Government through the BERC program; and by the Spanish Ministry of Economy and Competitiveness MINECO: BCAM Severo Ochoa accreditation SEV References ] R. Caflisch, The fluid dynamic limit of the nonlinear Boltzmann equation. Comm. Pure Appl. Math. 33 (980), no. 5, ] R. Chill, and A. Haraux, An optimal estimate for the time singular limit of an abstract wave equation. Funkcial. Ekvac. 47 (2004), no. 2, ] F. De Anna and A. Zarnescu, Global well-posedness and twist-wave solutions for the inertial Qian-Sheng model of liquid crystals. J. Differential Equations 264 (208), no. 2, ] J. Eells and J. H. Sampson, Harmonic mappings of Riemannian manifolds, American Journal of Mathematics, 86 (964), no., ] J. L. Ericksen, Conservation laws for liquid crystals. Trans. Soc. Rheology 5 96, ] J. L. Ericksen, Continuum theory of nematic liquid crystals. Res. Mechanica 2, (987), ] J. L. Ericksen, Liquid crystals with variable degree of orientation. Arch. Rational Mech. Anal. 3 (990), no. 2, ] E. Feireisl, E. Rocca, G. Schimperna and A. Zarnescu, On a hyperbolic system arising in liquid crystals modeling. arxiv: (206). 9] N. Jiang and Y. L. Luo, On well-posedness of Ericksen-Leslie s parabolic-hyperbolic liquid crystal model. arxiv: math.ap].

21 WAVE MAP LIMIT 2 0] F. M. Leslie, Some constitutive equations for liquid crystals. Arch. Rational Mech. Anal. 28 (968), no. 4, ] F. M. Leslie, Theory of flow phenomena in liquid crystals. The Theory of Liquid Crystals Vol. 4, pp. -8. Academic Press, London-New York, ] F-H. Lin and C. Liu, Static and dynamic theories of liquid crystal. J. Partial Differential Equations, 4 (200), no. 4, ] F.-H. Lin and C.-Y. Wang, The analysis of harmonic maps and their heat flows. World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, ] P-L. Lions, Mathematical Topics in Fluid Mechnics, Vol. : Incompressible Models. Clarendon Press. Oxford, ] J. Shatah and M. Struwe, Geometric Wave Equations. Courant Institute of Mathematical Sciences, New York University, New York and Amecican Mathematical Society, Providence, Rhode Island, (Ning Jiang) School of Mathematics and Statistics, Wuhan University, Wuhan, , P. R. China address: njiang@whu.edu.cn (Yi-Long Luo) School of Mathematics and Statistics, Wuhan University, Wuhan, , P. R. China address: yl-luo@whu.edu.cn (Shaojun Tang) School of Mathematics and Statistics, Wuhan University, Wuhan, , P. R. China address: shaojun.tang@whu.edu.cn (Arghir Zarnescu) BCAM, Basque Center for Applied Mathematics, Mazarredo 4, Bilbao, Spain address: azarnescu@bcamath.org