arxiv: v1 [math.ap] 28 May 2009

Transcription

1 Wave breaking in the short-pulse equation arxiv: v [math.ap] 8 May 9 Yue Liu, Dmitry Pelinovsky, and Anton Sakovich Department of Mathematics, University of Teas at Arlington, Arlington, TX, 769, USA Department of Mathematics, McMaster University, Hamilton, ON, L8S 4K, Canada Abstract Sufficient conditions for wave breaking are found for the short-pulse equation describing wave packets of few cycles on the ultra-short pulse scale. The analysis relies on the method of characteristics and conserved quantities of the short-pulse equation and holds both on an infinite line and in a periodic domain. Numerical illustrations of the finite-time wave breaking are given in a periodic domain. Introduction The short-pulse equation,.) u t = u+ 6 u3 ), t >,, is a useful and simple approimation of nonlinear wave packets in dispersive media in the limit of few cycles on the ultra-short pulse scale [, ]. This equation is a dispersive generalization of the following advection equation.) u t = u u, t >,. According to the method of characteristics, the advection equation.) ehibits wave breaking in a finite time for any initial data u,) = u ) on an infinite line if u ) is continuously differentiable and there is a point such that u )u ) >. This follows from the implicit solution uξ,t) = u ξ), ξ,t) = ξ u ξ)t, t >, ξ, for any given u ) C ). We say that the finite-time wave breaking occurs if there eists a finite time T, ) such that.3) lim sup t T u,t)u,t) =, while lim For the simple advection equation.), { T = inf ξ sup t T u ξ)u ξ) : u ξ)u ξ) > u,t) <. }.

2 In view of this result, we address the question if the dispersion term u in the short-pulse equation.) can stabilize global dynamics of the advection equation.) at least for small initial data. Local well-posedness of the short-pulse equation on an infinite line was proven in []. Theorem Schäfer & Wayne, 4). Let u H ). There eists a T > such that the short-pulse equation.) admits a unique solution ut) C[,T),H )) C [,T),H )) satisfying u) = u. Furthermore, the solution ut) depends continuously on u. To etend local solutions into a global solution, Pelinovsky & Sakovich [9] used the following conserved quantities of the short-pulse equation:.4) E := u d, u E := +u ).5) d = + d, +u [ )] u E := +u u.6) d = +u +u )5/d. If ut) is a local solution in Theorem, then E, E, and E are bounded and constant in time for all t [,T). Global well-posedness of the short-pulse equation on an infinite line was proven in [9] for small initial data in H satisfying.7) E +E u L + u L <. The condition.7) can be sharpen using the scaling transformation for the short-pulse equation.). Let α + be an arbitrary parameter. If u,t) is a solution of.), then ũ, t) is also a solution of.) with.8) = α, t = α t, ũ, t) = αu,t). The conserved quantities transform as follows: ) Ẽ = +ũ d = α Ẽ = +u ) d = αe, ũ u = α +ũ )5/d +u )5/d = α E. Finding the minimum of Ẽ+Ẽ = αe +α E in α gives a sharper sufficient condition on global well-posedness [9]. Theorem Pelinovsky & Sakovich, 9). Let u H ) and E E <. Then the short-pulse equation.) admits a unique global solution ut) C +,H )) satisfying u) = u.

3 Theorem does not eclude wave breaking in a finite time for large initial data and this paper gives a proof that the wave breaking may occur in the short-pulse equation.). Negating the sufficient condition for global well-posedness in Theorem, a necessary condition for the wave breaking follows: the wave breaking may occur in the short-pulse equation.) with the initial data u H only if E E. We shall find a sufficient condition for the wave breaking in the short-pulse equation.). Unlike the previous work in [9], we will not be using conserved quantity E but will rely on the conservation of E, E, and the energy [ E := u ) ].9) u4 d. Here u is defined from a local solution u by u := u,t)d = u,t)d = ) u,t)d, thanks to the zero-mass constraint u,t)d = for all t,t). Note that the initial data u does not have generally to satisfy the zero-mass constraint u )d =.) Thanks to the Sobolev inequality, u L 4 C u H for some C >, the quantity E is bounded if u H ) Ḣ ), where Ḣ is defined by its norm u Ḣ := u L. Note that if u H ) Ḣ ), then u)d = is satisfied.) Our main result on the wave breaking on an infinite line is formulated as follows. Theorem 3. Let u H ) Ḣ ) and T be the maimal eistence time of Theorem. Let F := E + 8E E +E 4 ) / ) /, F := E +E + ) / E F, and assume that there eists such that u )u ) > and F either u ) > 4F ) /3, u ) u ) > F + F u ) 3 / ) F, or F u ) 4F ) /3, u ) u ) > F. Then T <, so that the solution ut) C[,T),H ) Ḣ )) of the short-pulse equation.) blows up in the sense of limsup t T u,t)u,t) =, while lim t T u,t) L F. 3

4 emark. The quantities F and F in Theorem 3 can be defined by F := sup u,t) L, F := sup u,t) L. t [,T) t [,T) Note that the scaling transformation.8) gives ũ, t) L = α u,t) L, ũ, t) L = α u,t) L, so that the sufficient condition of Theorem 3 with new definitions of F and F is invariant in α. We note, however, that, while the bound on F in Theorem 3 scales correctly as F = αf, the bound on F is not correctly scaled in α because Ẽ = α 3 E, Ẽ = α 5 E. This is an artefact of using Sobolev embedding in Lemma 3 below. To prove Theorem 3, we shall adopt the method of characteristics and proceed with apriori differential estimates. Our results remain valid in a periodic domain, where Theorem 3 is replaced by Theorem 4 below. The technique of characteristics and apriori differential estimates has been applied for wave breaking in other nonlinear wave equations, see [, 3, 4, 5, 6, 7, ] for an incomplete list of references. We emphasize that, unlike the previous work in [9], we avoid using a transformation between the short-pulse equation.) and the integrable sine Gordon equation in characteristic coordinates. Our proof of the wave breaking for the short-pulse equation.) does not suggest that there eists a similar wave breaking for the sine Gordon equation, it is rather the breaking of the coordinate transformation between the two equations. On a similar note, we do not use the integrability properties of the short-pulse equation.) such as the La pair, the inverse scattering transform method, the bi-hamiltonian formulation, and the eistence of eact soliton solutions. The article is constructed as follows. The proof of Theorem 3 is given in Section. Section 3 reports etension of Theorem 3 to a periodic domain. Section 4 contains numerical evidences of the finite-time wave breaking in a periodic domain. Wave breaking on an infinite line Let us rewrite the Cauchy problem for the short-pulse equation on an infinite line in the form.) { ut = u u + u,, t > u,) = u ),, where u := u,t)d. In what follows, we use both notations ut) and u,t) for the same solution of the Cauchy problem.). Local eistence of solutions with the conservation of E and E is described by the following result. 4

5 Lemma. Let u H s ) Ḣ ), s. There eist a maimal time T = Tu ) > and a unique solution u,t) to the Cauchy problem.) such that ut) C[,T),H s ) Ḣ )) C [,T),H s )) satisfying u) = u. Moreover, the solution ut) depends continuously on the initial data u and the values of E, E, and E in.4),.5), and.9) are constant on [,T). Proof. If u H s ) Ḣ ), s, then u H s+ ), so that u )d =. By the theorem of Schäfer & Wayne [], there eists a solution of the short-pulse equation.), so that ut) C[,T),H s )) C [,T),H s )) ut) := u t u u C,T),H s )). Therefore, ut) C[,T),H s ) Ḣ )) in view of continuity of u Ḣ at t =. Because f H s ), s implies lim f) =, the zero-mass constraint holds in the form.) u,t)d =, t [,T). Let us define ut) := u ) t 6 u3. By the zero-mass constraint.) and uniqueness of the solution ut) for any t [,T), we obtain lim u,t) =, t [,T). Using balance equations for the densities of E, E, and E, we write [ u) ] [ u4 = u) ] +u ) t u ) t = t = 36 u6 u) + 4 u4 ) u ) +u Integrating the balance equation in for any t [,T), we complete the proof that E, E, and E are bounded and constant on [,T). emark. The maimal eistence time T > in Lemma is independent of s in the following sense. If u H s ) H s ) Ḣ ) for s,s and s s, then and ut) C[,T),H s ) Ḣ )) C [,T),H s )) ut) C[,T ),H s ) Ḣ )) C [,T ),H s )) with the same T = T. See Yin [] for standard arguments..,, 5

6 By using the local well-posedness result in Lemma and energy estimates, we obtain a precise blow-up scenario of the solutions to the Cauchy problem.). Lemma. Let u H ) Ḣ ) and ut) be a solution of the Cauchy problem.) in Lemma. The solution blows up in a finite time T, ) in the sense of lim t T u,t) H = if and only if limsup t T u,t)u,t) = +. Proof. We only need to prove the necessary condition, since the singularity in u,t)u,t) as t T implies the singularity in u,t) H as t T. Assume a finite maimal eistence time T, ) and suppose, by the contradiction, that there is M > such that.3) supu,t)u,t) M <, t [,T). Applying density arguments, we approimate the initial value u H ) by functions u n H3 ), n, so that lim n u n = u. Furthermore, write u n t) for the solution of the Cauchy problem.) with initial data u n. Using the regularity result proved in Lemma for s = 3, it follows from Sobolev s embedding that, if u n t) C[,T),H 3 ) Ḣ )), then u n,t) is a twice continuously differentiable function of on for any t [,T). It is then deduced from the short-pulse equation.) that d u n dt ) d = u n u n ) 3 d M u n ) d and d u n dt ) d = 5 u n u n un ) d 5M u n ) d. The Gronwall inequality yields for all t [,T), u n,t) L u n ) L e M t, u n,t) L u n ) L e 5 Mt. Since u n H converges to u H as n, we infer from the continuous dependence of the local solution ut) on initial data u that u,t) H remains bounded on [,T) for the solution ut) in Lemma. Therefore, the contradiction is obtained and either T is not a maimal eistence time or the bound.3) is not valid on [,T). We also show that the blow-up of Lemma is the wave breaking in the sense of condition.3). In other words, both u,t) L and u,t) L are uniformly bounded for all t [,T). Lemma 3. Let u H ) Ḣ ) and T > be the maimal eistence time of the solution u,t) in Lemma. Then,.4) u,t) L F, u,t) L F, t [,T), 6

7 where.5).6) F := E + 8E E +E 4 ) / ) /, F := E +E + ) / E F. Proof. For all t [,T) and the solution u,t), we have u,t) = uu d uu d u u + + +u d +u E / u + +u )d E / E +E u,t) ) /. L As a result, we obtain u,t) 4 L E E +E u,t) L, so that bound.5) is foundfrom the quadratic equation on u,t) L. On the other hand, we have u,t) H ) / = u,t) L + u,t) L = E +E + u,t) 4 L 4 E +E + E u,t) L. Using the Sobolev inequality u L u H and bound.5), we obtain bound.6). Let us introduce a continuous family of characteristics for solutions of the Cauchy problem.). Let ξ, t [,T), where T is the maimal eistence time in Lemma, and denote.7) = Xξ,t), u,t) = Uξ,t), u,t) = Gξ,t), so that.8) { Ẋt) = U, X) = ξ, { Ut) = G, U) = u ξ), where dots denote derivatives with respect to time t on a particular characteristics = Xξ,t) for a fied ξ. Applying classical results in the theory of ordinary differential equations, we obtain two useful results on the solutions of the initial-value problem.8). Conserved quantities E and E of the Cauchy problem.) are used to control values of U and G on the family of characteristics. 7

8 Lemma 4. Let u H ) Ḣ ) and T > be the maimal eistence time of the solution ut) in Lemma. Then there eists a unique solution Xξ,t) C [,T)) to the initial-value problem.8). Moreover, the map X, t) : is an increasing diffeomorphism for any t [,T) with ξ Xξ,t) = ep t ) uxξ,s),s)u Xξ,s),s)ds >, t [,T), ξ. Proof. EistenceanduniquenessofXξ,t) C,[,T)) followsfromtheintegral equation Xξ,t) = ξ t U ξ,s)ds, t [,T), ξ, since Uξ,t) C,[,T)) for the solution ut) in Lemma. Using the chain rule, we obtain t ) ξ Ẋξ,t) = Wξ,t) ξ Xξ,t) ξ Xξ,t) = ep Wξ,s)ds, where Wξ,t) = uxξ,t),t)u Xξ,t),t) C,[,T)). Therefore, ξ Xξ,t) > for all t [,T) and ξ. Let Vξ,t) = u Xξ,t),t), Wξ,t) = uxξ,t),t)u Xξ,t),t) Uξ,t)Vξ,t) and compute their rate of changes along the family of characteristics { V = VW +U,.9) Ẇ = W +VG+U. Let F,F > be fied in terms of conserved quantities E, E, and E as in Lemma 3 and assume that there eists ξ such that Wξ,) > and ) F /3 either Vξ,) >, Vξ,) Wξ,) > F + F Vξ,) 3 / 4F ) F, or ) F /3 Vξ,), Vξ,) Wξ,) > F. 4F Because of the invariance of system.9) with respect to G G, U U, V V, W W, itissufficienttoconsiderthecasewithvξ,) >. Wewillprovethat, undertheconditions above, Vξ,t) and Wξ,t) remain positive and monotonically increasing functions for all 8

9 t >, forwhichtheyarebounded,sothatvξ,t)andwξ,t)satisfytheaprioridifferential estimates { V VW F,.) Ẇ W VF. In what follows, we use Vt) and Wt) instead of Vξ,t) and Wξ,t) for a particular ξ. The following lemma establishes sufficient conditions on the initial point V), W)) that ensure that a lower solution satisfying { V = VW F,.) Ẇ = W VF, goes to infinity in a finite time. Lemma 5. Assume that the initial data for system.) satisfy.) either F ) /3 V) >, V)W) > F + F V 3 ) / 4F ) F,.3) or ) F /3 < V), V)W) > F. 4F Then the trajectory of system.) blows up in a finite time t, ) such that Vt) and Wt) are positive and monotonically increasing for all t [,t ) and there is C > such that.4) lim t t Vt) =, lim t t Wt) =, and lim t t t t)vt) = C. Moreover, t is bounded by V).5) t { / }. min V), V ) F V 3 )+ ) F Proof. Let us first consider the homogeneous version of system.) for F =, that is.6) { V = VW, Ẇ = W VF. Of course, F is never zero, otherwise E =. This case is used merely for illustration, since eplicit solutions can be obtained for F =, whereas qualitative analysis has to be developed for F. System.6) is integrable since W = V V d dt V V ) = F V = V C F t), where C = W)/V). Integrating the last equation for Vt), we obtain the eplicit solution of the truncated system, Vt) = V) CV)t+ F V)t, Wt) = C F t)vt). 9

10 The solution reaches infinity in a finite time t, ) if V) > and C V) > F. Note that these conditions coincide with condition.) for F =. Also note that t is the first positive root of CV)t+ F V)t = so that Ut ) = C F t > and t = W) W ) F V) F V) = V) V)+ V ) F V 3 ) V) V ) F V 3 ). Consider now the full system.). Let and rewrite the system in the form V =, W = y.7) { ẋ = F y, ẏ = F y F. Epressing y from the system, we can rewrite it in the scalar form.8) y = F ẋ ẍ = f)+3f ẋ, where f) = F F 3. The only critical point of equation.8) is,ẋ) =,) where = F /F )/3 is the root of f). We need to show that there is a domain D in the phase plane,ẋ), so that all initial data in D generate trajectories in D that cross the vertical line = in a finite time. To do so, we construct a Lyapunov function for equation.8) in the form E,ẋ) = ẋ F + 4 F 4. The function E,ẋ) has a global minimum at,). For any solution t), we have d dt E,ẋ) = 3F ẋ > for >. The zero level of the Lyapunov function E,ẋ) = passes through the points,) and,), where = 4 /3 > see Figure ). It is clear that E,ẋ) > in the domain D = {,ẋ) : >, ẋ < σ)}, where σ) = { F F 4) /,, ),, [, ). We note that the condition,ẋ) D is equivalent to the condition { < <, y > F + F F 4) /, >, y > F,

11 which is nothing but the set of conditions.) and.3) at t =. By continuity, if,ẋ) D at t =, then,ẋ) remains in D for some time t >. No critical points of system.8) are located in D and t) is decreasing function for any t > as long as the trajectory stays in D. ecall that E,ẋ) > and d dt E,ẋ) > for any,ẋ) D. A trajectory in D can not cross ẋ = σ) because E,σ)) = for, ) and ẍ = f)+3f ẋ < for > and ẋ <. Therefore, the trajectory either reaches = in a finite time t, ) or escapes to ẋ = for >. To eliminate the last possibility, we note that d ẋ 3F ) dt = f), so that ẋt) ẋ) 3F )+tf)) > for any finite time interval. Moreover, ẋ is bounded from zero in D by the level curve E,ẋ) = E),ẋ)), which is a conve curve in D. Therefore, ẋt) ma{ẋ),ρ}, t >, as long as >, where ρ < is uniquely found from E, ρ) = E), ẋ)), that is from the point of intersection of the level curve of E, ẋ) = E), ẋ)) with the negative ẋ-ais. Therefore, t) )+tma{ẋ),ρ}, so that t) reaches = in a finite time t, ) for any trajectory in D. Moreover, finding ρ eplicitly gives the bound on the blow-up time < <, t ) ẋ ) F )+ F 4 )) /, >, t ) ẋ), which becomes bound.5) after the return back to variable Vt). Since V =, W = y, and y = F ẋ >, we have lim t t Vt) = and lim t t Wt) = for any trajectories in D. Since ẋ < for the trajectory in D, we also have t) t t) as t t so that there eists C > such that lim t t t t)vt) = C. It remains to show that Vt) and Wt) are monotonically increasing functions on [,t ). To do so, we write where V = VW F = ẋ, Ẇ = W F V = g,ẋ), g,ẋ) = ẋ F ẋ f).

12 ẋ E, ẋ) = σ) D Figure : Domain D in the phase plane,ẋ) of equation.8). For any trajectory in D, ẋt) < so that Vt) >. Furthermore, since g,ẋ) is zero at a curve outside the domain D, because g,σ)) = F + F 4 +F σ) >,, ], then Ẇt) > for any trajectory in D. ecalling that W = UV in system.9) and using the bound.4), we obtain the upper bound for any solution at the family of characteristics.9) V = UV +U F V +F. We can now show that any upper solution with V) > goes to infinity in a finite time. Lemma 6. Consider.) Vt) = F V +F with V) >. There eists t, ) such that Vt) is positive, monotonically increasing for all t [,t ) and there eists C > such that Moreover, t /W). limt t)vt) = C. t t Proof. Since V > for any V, Vt) is monotonically increasing function. To show that Vt) reaches in a finite time, one can integrate the separable equation.) eplicitly and obtain Vt) = tanarctanv)+f t), so that since sup + cot ). t = π/ arctanv) F F V) W),

13 Applying results of Lemmas 5 and 6, we conclude the proof of Theorem 3. Proof of Theorem 3. Let V, W) satisfy system.9) corresponding to the characteristics with ξ. Let V,W) be the lower solution of system.) in Lemma 5 with V) = V) and W) = W). Let V be the upper solution of equation.) in Lemma 6 with V) = V). Let t be the blow-up time of the lower solution and t be the blow-up time of the upper solution. The upper bound for the solution of system.9) follows from the comparison principle for the differential equations since VW +U = V +) U V +)F which implies that Vt) Vt) for all t [,t ), for which Vt) remains bounded. To obtain the lower bound, we note that V V, W W Let V = [V,W] T, V = [V,W] T, and { V VW F VW F = V, Ẇ +F V W W = Ẇ +F V. L = [ F ] be a nilpotent matri of order one, so that e tl = I +tl. Thus, we write V V, W W d e tl V ) d e tl V ). dt dt Integrating this equation in t >, we infer that e tl Vt) e tl Vt), and since e tl is invertible for any t, we conclude that Vt) Vt), Wt) Wt), for all t [,T) [,t ), for which V,W) remain finite. Therefore, V,W) become infinite as t T and T [t,t ]. emark 3. The bounds on t and t in Lemmas 5 and 6 are inconclusive to compare T with the eact time of blow-up T := W) along a particular characteristic of the dispersionless advection equation.). 3 Wave breaking in a periodic domain Consider now the Cauchy problem for the short-pulse equation.) in a periodic domain 3.) { ut = u u + u, S, t >, u,) = u ), S, 3

14 wheresisaunitcircleequippedwithperiodicboundaryconditionsand is the mean-zero anti-derivative in the form u := u,t)d u,t)d d. S Local well-posedness and useful conserved quantities for the Cauchy problem 3.) in a periodic domain are obtained in the following lemma. Lemma 7. Assume that u H s S), s and S u )d =. Then there eist a maimal time T > such that the Cauchy problem 3.) admits a unique solution ut) C[,T),H s S)) C [,T),H s S)) satisfying u,) = u ) and Su,t)d = for all t [,T). Moreover, the solution ut) depends continuously on the initial data u and the quantities E = u d, E = +u d are constant for all t [,T). S Proof. Eistence of the solution u,t) and continuous dependence on u is proved on S similarly to what is done in Lemma on. To prove the zero-mass constraint, we note u t t) C,T),H s S)), u u t) C[,T),H s S)), so that for all t,t), we have u,t)d = t d+ S Su u u ) d =. S Initial values of E and E are bounded if u H s S), s. Conservation of E and E on [,T) follows from the balance equations u ) = t u) + ) 4 u4, ) ) +u = t u +u thankstothecontinuity andtheperiodicboundaryconditionsfor ut) C,T),Hs+ S)), ut) C,T),H s S)), and u t) C,T),H s S)) in on S if s. emark 4. The assumption S u )d = on the initial data u in the periodic domain S is necessary as it follows from the following apriori estimate u,t)d u )d ut) u L S), t,t). S S Note that S u,t)d =, for all t,t) and ut) C[,T),Hs S)), s. Hence the above estimate implies that S u )d =. Note that no zero-mass constraint is necessary on an infinite line in Theorem. 4 S,

15 The blow-up scenario for the solutions to the Cauchy problem 3.) coincides with the one in Lemma after the change of by S. The main result of this section is the proof of the finite-time wave breaking in a periodic domain, according to the following theorem. Theorem 4. Let u H S) and S u )d =. Assume that there eists such that u )u ) > and either u ) > or u ) E 4E / E 4E / ) /3, u ) u ) > E + E / u ) 3 / ) E, ) /3, u ) u ) > E. Then there eists a finite time T, ) such that the solution ut) C[,T),H S)) of the Cauchy problem 3.) blows up with the property limsup t T S u,t)u,t) = +, while lim t T u,t) L E. Proof. Let T > be the maimal time of eistence of the solution ut) C[,T),H S)) to the Cauchy problem 3.) constructed in Lemma 7. Since Su,t)d =, for each t [,T) there is a ξ t [,] such that uξ t,t) =. Then for S and t [,T), we have u,t) = u,t)d u,t) d E. ξ t S Since ut) C[,T),H 3 S)) is the mean-zero periodic function of for each t [,T), there eists another ξ t [,] such that u ξ t,t) =. Then for S and t [,T), we have u,t) = u,t)d u,t) d E. ξ t Therefore, bounds.4) are rewritten with F := E, F := E The rest of the proof follows the proof of Theorem 3. S 4 Numerical evidence of wave breaking The goal of this section is to complement the analytic results by several eamples and numerical computations. More specifically, we first show that the sufficient condition for wave breaking in Theorem 3 is not satisfied for the eact modulated pulse solution to the short-pulse equation which is known to remain globally bounded in space and time. Then we consider the interplay between global well-posedness and wave breaking of Theorems and 3 for a class of decaying data on an infinite line. Finally, we perform numerical simulations 5

16 in a periodic domain for a simple harmonic initial data and thus give illustrations to the sufficient condition for wave breaking in Theorem 4. Theorem 3 gives a sufficient condition for formation of shocks in the short-pulse equation.) on the infinite line. Let us show that this condition is not satisfied for eact modulated pulse solutions obtained in[8, ]. The simplest one-pulse solution is given in the parametric form u,t) = Uy,t), = Xy,t), where 4.) Uy,t) = 4mn msinψsinhφ+ncosψcoshφ m sin ψ +n cosh, φ msinψ nsinhφ Xy,t) = y +mn m sin ψ +n cosh φ, m,) is an arbitrary parameter, n = m, and φ = my +t), ψ = ny t). The pulse solution enjoys the periodicity property { Uy,t) = U y π m,t+ π m), Xy,t) = X y π m,t+ π m) + π m y,t). y,t), and an eponential decay in any direction transverse to the anti-diagonal on the y, t)-plane. Since X y = 8m n sin ψcosh φ m sin ψ +n cosh φ) = cos 4arctan msinψ ), ncoshφ the function = Xy,t) is invertible in y for all t if msinψ ncoshφ < tan π m 8 n tan π 8, that is for all m,m cr ), where m cr = sin π For these values of m, the pulse solution u,t) is analytic in variables,t), has the space-time periodicity u,t) = u π m,t+ π ),,t), m and the eponential decay in the transverse direction to the anti-diagonal in the, t)-plane. The graph of a nonsingular pulse solution for m =.3 is shown on Figure left). emark 5. Coordinate y in the eact solution 4.) is different from coordinate ξ in the method of characteristics because Xy,) y. Nevertheless, Xy,t) and Uy,t) satisfy the same set of equations X t = U y,t), U t = u =Xy,t), so that ξ and y are uniquely related by the representation ξ = Xy,). If y is found as a function of ξ, the initial data of the Cauchy problem.) is found from u ξ) = Uy,). 6

17 Since u,t) is analytic in for any fied t and decays to zero eponentially fast at infinity, it is clear that u,t) H ). Furthermore, since u = u t u u, it also follows that u,t) Ḣ ). We compute numerically bounds F and F using the eact solution 4.) and new definitions F := sup t u,t) L, F := sup u,t) L. t It follows from emark that bounds F and F defined above preserve the sufficient condition of Theorem 3 with respect to the scaling transformation.8). 4.) Let us define [ f := sup u ) u ) ] F, I [ f := sup I u ) u ) F F u ) 3 ) ] / F, where I = I = { F : u ) 4F { : u ) > F 4F ) /3, u )u ) >, } ) /3, u )u ) >, }., According to Theorem 3, wave breaking occurs if either f or f is positive. For the eact modulated pulse solution 4.) at t =, the numerical calculations show that the set I is empty and the quantity f is strictly negative for any m,m cr ), see Figure right). Therefore, the sufficient condition for the wave breaking in Theorem 3 is not satisfied, which corresponds to our understanding that the eact modulated pulse solutions 4.) remain bounded for all,t). We note, however, that the sufficient condition for the global well-posedness in Theorem is satisfied only for pulses with m, 3), since E E = 3m. This computation shows that the sufficient condition of Theorem is not sharp. Net we compare the sufficient conditions for the global well-posedness and wave breaking in Theorems and 3 for a class of initial data 4.3) u ) = a b )e b, a >, b >, where parameters a and b determine the amplitude and steepness of u. Note that the zero-mass constraint.) is satisfied by u and it is clear that u H ) Ḣ ). The conserved quantities E and E can be computed analytically, E = a π 56 5a b ) 48, E = 3a π b 3 8, b whereas the conserved quantities E and E are not epressed in a closed form. Using numerical approimations of the integrals, we determine the boundary of the well-posedness 7

18 u,t) 8 6. f.4 4 t m Figure : The eact modulated pulse solution 4.) of the short-pulse equation.) for m =.3 left). The quantity f is negative for any m,m cr ) right). regioninthea,b)-planebyfindingtheparametersaandbfromthecondition E E =. We also compute the boundaryof the wave breaking region in the a,b)-plane by computing f and f in 4.). Unlike the case of modulated pulses, we find that the set I is empty and f may change the sign along the curve on the a,b)-plane. The two boundaries are shown on Figure 3, where we can see that the two regions of global well-posedness and wave breaking are disjoint. Finally, we perform numerical simulations of the periodic Cauchy problem 3.) with the -periodic initial data 4.4) u ) = acosπ, a >. The two conserved quantities E and E in Lemma 7 are computed analytically as E = a, E = π Eπai), where E stands for a complete elliptic integral. Using the above conserved quantities we find out that the sufficient condition for the wave breaking in Theorem 4 is satisfied for a >.53. To illustrate the behaviour of a solution u,t) to the Cauchy problem 3.), we perform numerical simulations using a pseudospectral method. When the parameter a is sufficiently small, the value of Wt) := supu,t)u,t) S remains bounded as shown on the top panel of Figure 4 for a =.. On the other hand, when a becomes larger, the wave breaking occurs, as on the bottom panel of Figure 4 for a =.5. On the bottom right panel of Figure 4 we show using the linear regression that the curve W t) is fitted well with the straight line A + Bt for some coefficients A,B). 8

19 Figure 3: Boundaries of the global well-posedness and the wave breaking in the Cauchy problem.) with initial data 4.3): the global well-posedness occurs below the lower curve and the wave breaking occurs above the upper curve..8 u,t) supu u ). t t u,t) supu u ) t.5.5 T t Figure 4: Solution surface u,t) left) and the supremum norm Wt) right) for a =. top) and a =.5 bottom). The dashed curve on the bottom right picture shows the linear regression with C =.7, T =

20 T.3 C a a Figure 5: Time of wave breaking T versus a left). Constant C of the linear regression versus a right). Thus, we make a conclusion that Wt) C T t for < T t, where C = B and T = AB. Using the linear regression, we also obtain pairs T,C) for different values of a. The results are shown on Figure 5. Note that the constant C approach as a gets larger. This observation is consistent with the eact blow-up law Wt) = T t obtained for the dispersionless advection equation.) using the method of characteristics. eferences [] Y. Chung, C.K..T. Jones, T. Schäfer, and C.E. Wayne, Ultra-short pulses in linear and nonlinear media, Nonlinearity 8, ) [] A. Constantin, J. Escher, Well-posedness, global eistence, and blowup phenomenon for a periodic quasi-linear hyperbolic equation, Comm. Pure Appl. Math., 5 998), [3] A. Constantin and J. Escher, Wave breaking for nonlinear nonlocal shallow water equations, Acta Mathematica, 8, ). [4] D. Chae, A. Cordoba, D. Cordoba, and M.A. Fontelos, Finite time singularities in a D model of the quasi-geostrophic equation, Adv. Math., 94, 3 3 5). [5] J. Hunter, Numerical solutions of some nonlinear dispersive wave equations, Lectures in Appl. Math., 6 99), [6] Y. Liu and Z.Y. Yin, Global eistence and blow-up phenomena for the Degasperis Procesi equation, Comm. Math. Phys. 67, 8 8 6). [7] Y. Liu and Z.Y. Yin, On the blow-up phenomena for the Degasperis Procesi equation, Inter. Math. es. Not. 7, pages 7).

21 [8] Y. Matsuno, Periodic solutions of the short pulse model equation, J. Math. Phys. 49, 7358, 8 pp. 8) [9] D. Pelinovsky, A. Sakovich, Global well-posedness of the short-pulse and sine Gordon equations in energy space, perprint, 8), 4 pages. [] A. Sakovich and S. Sakovich, Solitary wave solutions of the short pulse equation, J. Phys. A: Math. Gen. 39, L36 L367 6). [] T. Schäfer and C. E. Wayne, Propagation of ultra-short optical pulses in cubic nonlinear media, Physica D, 96 4), 9 5. [] Z. Yin, On the Cauchy problem for an integrable equation with peakon solutions, Illinois J. Math. 47, ).