arxiv: v1 [math.ap] 27 Jun 2017

Transcription

1 ON WHITHAM AND RELATED EQUATIONS CHRISTIAN KLEIN, FELIPE LINARES, DIDIER PILOD, AND JEAN-CLAUDE SAUT arxiv:76.87v [math.ap] 7 Jn 7 Abstract. The aim of this paper is to stdy, via theoretical analysis and nmerical simlations, the dynamics of Whitham and related eqations. In particlar we establish rigoros bonds between soltions of the Whitham and KdV eqations and provide some insights into the dynamics of the Whitham eqation in different regimes, some of them being otside the range of validity of the Whitham eqation as a water waves model.. introdction The Whitham eqation was introdced formally in [6] as an alternative to the KdV eqation, by keeping the eact dispersion of the water waves system. No rigoros derivation from the water waves system is known, and there is no consistent scaling allowing to derive it. Actally the idea in [6] was to propose a model describing also the occrrence of waves of greatest height with the Stokes degrees angle at the crest which is clearly impossible with the KdV eqation. As was conjectred in [6] and later proved in [4], this is actally the case (with a different angle thogh) bt this wave breaking phenomena seems to be irrelevant in the (KdV) regime where the Whitham eqation can be really viewed as a consistent water wave model, as will be discssed in the present paper. In short, and as was noticed by Whitham, the desired qalitative effect is in bt this wave breaking effect cannot be qantitatively linked with the water waves system and ths to real waves. From the modeling point of view, the advantage of the Whitham eqation with respect to KdV seems to be the enlargement of the range of validity of this asymptotic model in terms of freqencies thogh this is not so easy to qantify. On the other hand, de to the very different behavior of the dispersion at low and large freqencies, the Whitham eqation is a fascinating mathematical object since, as will be discssed later, it has three interesting different asymptotic regimes. One shold ths distingish between the seflness of the Whitham eqation as a relevant water waves model which seems to be poor (see [5, 4] for some preliminary comparisons with eperiments) and its mathematical properties that are rich, in particlar as a sefl and relevant toy model to Date: Jne 5, 7.

2 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT provide some insights into the effect of adding a weak dispersion term to a nonlinear hyperbolic eqation. We will restrict in this paper to the Whitham eqation posed on the real line. Many interesting aspects of the periodic problem, in particlar concerning periodic solitary waves, can be fond in [, 4, 8, 53, 55]. However or reslts concerning the Cachy problem and its link with the KdV one are straightforwardly valid in the periodic case. In order to make the link with the KdV eqation (and ths with the fll water wave system) we will write the Whitham eqation as in [4] : Whit (.) t + L ɛ + ɛ =, where the non local operator L ɛ is related to the dispersion relation of the (linearized) water waves system and is defined by L ɛ = l( ( ) / tanh ɛ D ɛd) := and D = i = i ɛ D. The (small) parameter ɛ measres the comparable effects of nonlinearity and dispersion (see below). The Whitham eqation is spposed to be an approimation of the fll water waves system on time scales of order ɛ, see [4] and Section below. Taking the formal limit ɛ ξ in L ɛ, (.) redces to the KdV eqation KdV (.) t + + ɛ + ɛ 6 =. We do not know of a complete rigoros jstification of the Whitham eqation from the water waves system. More precisely no correct scaling seems to eist allowing to connect directly and rigorosly the Whitham eqation and the water wave system. See however [4] pages 3-4 and Section below where a comparison with the KdV eqation is displayed, jstifying ths the Whitham eqation via the KdV approimation of weakly nonlinear long srface water waves. On the other hand the Whitham eqation can be viewed as the onedimensional restriction of the fll dispersion KP eqation introdced in [4] to overcome the bad behavior of the dispersion relation of the sal KP eqations at low freqencies in (see also the analysis in [4]). We refer to [44] for a frther stdy of the Cachy problem and to [] for the eistence of localized solitary waves, close to the sal KP I ones in the case of strong srface tension): FDbis (.3) t + c W W ( ( ) / ɛ D ɛ ) + ɛ D D + ɛ 3 =,

3 THE WHITHAM EQUATION 3 with c W W ( ɛk) = ( + βɛk ) ( tanh ɛk ɛk ) /, where β is a dimensionless coefficient measring the srface tension effects and D ɛ = D + ɛd, D = i, D = i y. The link with the fll water wave system is via the choice of parameters (see [4]). Denoting by h a typical depth of the flid layer, a a typical amplitde of the wave, λ and λ y typical wave lengths in and y respectively, the relevant regime here is when ɛ a ( ) ( ) h λ h. λ y λ For prely gravity waves, β =, (.3) becomes FDter (.4) t + c W W ( ( ) / ɛ D ɛ ) + ɛ D D + ɛ 3 =, with c W W ( ɛk) = ( tanh ɛk ɛk ) /, which redces to the Whitham eqation (.) when does not depend on y. In presence of srface tension (β > ) (.3) redces when does not depend on y to the capillary Whitham eqation (see [53, 6]) WhitST (.5) t + L ɛ + ɛ =, where ( ) / tanh ɛ D dispst (.6) Lɛ = ( + βɛ D ) /, β >. ɛ D Taking the formal limit ɛ ξ in L ɛ, (.5) redces to the KdV eqation KdVST (.7) t + + ɛ ( ) 3 β + ɛ =. Note also that the Whitham eqation with srface tension looks, for high freqencies, like the following fractional KdV (fkdv) eqation with α = / fkdvst (.8) t + + D α =. We refer to [37, 4, 43] for some properties of the fkdv eqations viewed as toy models to stdy the inflence of a weak dispersive pertrbation on the dynamics of the Brgers eqation.

4 4 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT There is a Bossinesq like, system version of (.) for waves propagating in both directions. As for the Whitham eqation it cannot been derived directly from the water wave system by a consistent asymptotic analysis bt by a rather heristic and formal argment. Actally, one obtains a fll dispersion system when in the Bossinesq regime one keeps (formally) the original dispersion of the water waves system (see [4], [8], and [] where interesting nmerical simlations of the propagation of solitary waves are performed). ( ) / Setting again L ɛ = tanh ɛ D ɛ D, we get with D = i or i : FDd (.9) FDd (.) when d = and when d =. { ηt + L ɛ + ɛ(η) = t + η + ɛ =, { ηt + L ɛ + ɛ (η) = t + η + ɛ =, Taking the limit ɛ ξ in L ɛ, (.9) redces formally to { ηt + Kap (.) + ɛ 3 + ɛ(η) = t + η + ɛ =, while in the two-dimensional case, (.) redces in the same limit to { ηt + + ɛ FDdbis (.) 3 + ɛ (η) = t + η + ɛ =, that is to the (linearly ill-posed) system one gets by epanding to first order the Dirichlet to Nemann operator with respect to ɛ in the fll water wave system (see [4]). On the other hand, (.9) and (.) are linearly well-posed. As will be seen below the Cachy problem for the nonlinear system is locally well-posed nder a non physical condition on the initial data thogh. System (.) is also known in the Inverse Scattering commnity as the Kap-Kpperschmidt system (see [3, 39]). It is completely integrable thogh linearly ill-posed since the eigenvales of the dispersion matri are ±iξ( ɛ 3 ξ ) / (we refer to [4] for an analysis of the ill-posedness of the nonlinear Kap and related systems). It has eplicit solitary waves (see [3]). The Bossinesq system (.9) can therefore be seen as a (well-posed) reglarization of the Kap-Kpperschmidt system. It is not known to be completely integrable (this is nlikely, see in particlar the simlations in Section 7). The fll dispersion Bossinesq system has the following Hamiltonian strctre ( ) η t + Jgrad H ɛ (η, ) = As noticed in [], the se of nonlocal models for shallow water waves is also sggested in [6].

5 THE WHITHAM EQUATION 5 where when d = and H ɛ (U) = J = y y, ( Lɛ + η + ɛη ) ddy, R ( η U =, ) ( ) η t + Jgrad H ɛ (η, ) = where ( ) J = and H ɛ (η, ) = ( L ɛ + η + ɛ η)d, R when d =. We will see in the net section that the fll dispersion Bossinesq system has mathematical properties that make it dobtfl as a relevant water wave model. When srface tension is taken ( into accont, one shold replace the operator L ɛ by P ɛ = (I + βɛ D ) ) tanh( ɛ D ) ɛ D where again the parameter β > measres srface tension (see [4]), yielding a more dispersive fll dispersion Bossinesq system. When β > 3, this fll dispersion Bossinesq system yields in the limit ɛ ξ in P ɛ, Bossinesq systems of the class a <, b = c = d = (see [7]) for which long time (that is on time scales of order /ɛ) well-posedness is established in [58], Theorem 4.5. When β < 3, the fll dispersion Bossinesq system redces in the formal limit ɛ ξ in P ɛ, to an ill-posed system, analogos to the Kap system in dimension. Remark.. Another Whitham-Bossinesq system is introdced in [9]. It writes HPFDd (.3) HPFDd (.4) when d = and { ηt + + ɛ(η) = t + P ɛ η + ɛ =, { ηt + + ɛ (η) = t + P ɛ η + ɛ =, This system seems a bit artificial contrary to (.9), (.) which are the eact conter parts of the Whitham eqation with respect to the original Bossinesq system

6 6 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT when d =. In the limit ɛ ξ in L ɛ, those systems redce formally to limhpd (.5) { ηt + + ɛ(η) = t + η + ɛ( 3 β)η + ɛ =, limhpd (.6) in dimension one and to { ηt + + ɛ (η) = t + η + ɛ( 3 β) η + ɛ =, in the two-dimensional case. Both systems are ill-posed when β < 3 while when β > 3 they belong to the class of a = b = d =, c < of classical (a,b,c,d) Bossinesq systems for which eistence on time scales of order /ɛ is established in [58], Theorems 4.6 and 4.7. The local well-posedness for (.3) is established in [9]. According to previos theoretical reslts and nmerical simlations, one epects at least three different regimes for the Whitham eqation (.) (withot srface tension):. Scattering for small initial data (see the simlations in [37] and in the last Section 7 of the present paper).. Finite time blow-p (csplike), see [37] for varios simlations displaying the strctre of a shock like blow-p (blow-p of the gradient with bonded sp-norm of the soltion). The occrrence of sch phenomena is rigorosly proven in [3] for a class of fractional KdV eqations and in [6] for the Whitham eqation itself (see also the nmerical simlations in [37]) bt these phenomena have probably nothing to do with the breaking of real water waves. In fact, when one keeps the small parameter ɛ in the eqation (that is not done in [6, 3]), the blow-p shold occr on time scales mch larger than /ɛ, the time scale on which the Whitham eqation is spposed to approimate the fll water wave system via the KdV eqation. We refer again to the simlations in Section A KdV, long wave regime. In fact it is shown in [] that (.) possesses specific solitary waves, close to those of KdV and formally stable. In this regime one can epect a KdV like behavior, namely the soliton resoltion, at least on sfficiently long time scales. Those solitary waves and their pertrbations are investigated nmerically in Section 7 of the present paper. One aim of the present paper is to give frther evidence of the relevance of those conjectres. The dynamics of the Whitham eqation with srface tension (.5) shold be different becase of the different behavior of the dispersion at high freqencies. This makes the eqation more dispersive and the epected dynamics is that of L critical KdV type eqations. In particlar the (epected)

7 THE WHITHAM EQUATION 7 finite time blow-p shold be similar to that of the L critical generalized KdV eqation cgkdv (.7) t =, or of the modified Benjamin-Ono eqation (also L critical) mbo (.8) t + H =, where H denotes the Hilbert transform. The rigoros analysis of blow-p for those eqations can be fond respectively in [46] and [48]. As aforementioned one aim of the present paper is to give some evidence for the above conjectres via mathematical analysis and mainly by carefl nmerical simlations. We will also give some hints on the qalitative behavior of the fll-dispersion Bossinesq systems (.9), (.). The paper is organized as follows. In a first section we give the epected error estimates on the correct time scales, between the soltions of the Whitham and KdV eqation for smooth initial data. Together with the classical reslts on the KdV approimation of srface water waves (see [4]) this implies a rigoros jstification of the Whitham eqation in the Bossinesq-KdV regime. The net section concerns the Cachy problem for the Whitham eqation with srface tension. Contrary to (.) dispersive estimates can be sed here to enlarge the space of resoltion to the Cachy problem. We then comment on the local well-posedness of the Cachy problem for the Bossinesq fll dispersion systems (.9) and (.) and of their capillary waves versions. We will see that the well-posedness of the system for pre gravity waves is obtained nder a very restrictive, non physical condition (positivity of the wave elevation). When this condition is not satisfied the system is (Hadamard) ill posed and those facts invalidate it as a relevant model for water waves. The presence of srface tension on the other hand prevents the appearance of Hadamard nstable modes when the initial elevation is not positive, the possible nstable modes being bonded then. The two net sections review (and comment on) known reslts concerning finite time blow-p and solitary wave soltions. Finally the last two sections display many accrate nmerical simlations aiming to illstrate and to detail varios properties of the Whitham eqations and systems, allowing to propose convincing conjectres on their dynamics. Section 7 is devoted to the Whitham eqation. We first constrct nmerically the solitary waves to the Whitham eqation withot srface tension and simlate their pertrbations for varios vales of the small parameter ɛ.

8 maintheo 8 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT Then we solve nmerically the Cachy problem for Gassian initial data λ ep( ). Depending on the size of ɛ and λ a finite time blow-p may occr, at a time otside the physically relevant time scales O(/ɛ) thogh. The sitation is qite different when srface tension is inclded. Actally one shows a finite time blow-p very similar to the one of the L critical KdV eqation t =. In Section 8 we consider the one-dimensional Whitham-Bossinesq systems, with and withot srface tension. We constrct nmerically solitary waves and stdy their stability. Then we simlate soltions of the Cachy problem with initial data satisfying or not the well-posedness condition.. Comparison between the Whitham and KdV eqations We will compare here the soltions v and of respectively the Whitham eqation (.) and the KdV eqation (.) with the same initial data v = = φ H (R). Here is the main reslt of this section. Theorem.. Let φ H (R). Then, for all j N, j, there eists M j = M j ( φ H j+8) > sch that maintheo. (.) ( v)(t) H j M j ɛ t, for all t ɛ. Remark.. The implicit constant in the notation t ɛ depends on φ H for j = and and on φ H j+ for j. Remark.. A similar resl holds mtatis mtandi for the periodic problem since we se energy type methods. Recall that sch a theorem was proved by Bona, Pritchard and Scott [8] for the comparison between BBM and KdV. We also refer to Albert and Bona [] for other comparison reslts in the long wave regime... The linear case. In this section, we compare the linear versions of (.) and (.), i.e. the Airy eqation Airy (.) t + + ɛ 6 3 =, and the linear Whitham eqation linwhitham (.3) t v + l( ɛd) v =, where l( ɛd) denotes the Forier mltiplier of symbol l( ɛξ) defined by l( ( ) / tanh ɛξ ɛξ) =. ɛξ Then, we have the following reslt

9 lintheo THE WHITHAM EQUATION 9 Theorem.. Let φ H (R). Then, for all j N, j, there eists N j = N j ( φ H j+7) > sch that the soltions of (.) and v of (.3) associated to the same initial datm φ satisfy lintheo. (.4) ( j v)(t) L N j ɛ ( + t), symbcomp symbcomp. (.5) for all t R +. Remark.3. Note that we cold also obtain bonds for the differences of and v in L arging as in the proof of Theorem.. The following technical reslt which compares the symbol of l( ɛd) and + ɛ 6 3 will be needed below. Lemma.3. Assme that ɛ ξ. Then, l( ɛξ)ξ ( ξ ɛ 6 ξ3) ɛ ξ 5. Proof. The proof follows directly from the epansion ( ) tanh() = 6 + O( 4 ), for <. Proof of Theorem.. The soltions of (.) and (.3) are respectively given by the nitary grops Airygrop (.6) (, t) = e t( + ɛ 6 3 ) φ() = e it(ξ ɛ 6 ξ3) e iξ φ(ξ) dξ, and linwhithamgrop (.7) v(, t) = e tl( ɛd) φ() = R R e itl( ɛξ)ξ e iξ φ(ξ) dξ. Then, it follows from the mean vale ineqality and Lemma.3 that ( j ) (, t) v(, t) e ɛ ξ < it(ξ ɛ6 ξ3) e itl( ɛξ)ξ ξ j φ(ξ) dξ + ξ j φ(ξ) ɛ ξ > tɛ ξ j+5 φ(ξ) dξ + ɛ ξ j+4 φ(ξ). ɛ ξ < ɛ ξ > Then, we dedce from the Cachy-Schwarz ineqality that j ( v)(, t) L ɛ ( + t) φ H j+7, which concldes the proof of the theorem. Remark.4. Note that fndamental soltions of both the linear Whitham and KdV eqations have a qite different behavior. The KdV one is given by the Airy fnction G KdV (, t) = C ( Ai ɛt) /3 ( t ( ɛt) /3 ).

10 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT Set 3 K(, t) = ( ) / e iξ e itξ tanh ξ ξ dξ. π R The fndamental soltion of the Whitham eqation is ths K ɛ (, t) = ( ) K ɛ ɛ, t. ɛ One can easily establish that contrary to the Airy fnction, K is an nbonded fnction. Here we follow the analysis in [9] in a different setting. Actally, we write ( tanh ξ ξ ξ ) / = (sign ξ) ξ / ( + e ξ ) / = (sign ξ) ξ / + r(ξ) where r is a continos fnction eponentially decaying to zero at ±. Using the elementary identity one obtains the decomposition K(, t) = e it(signξ) i ξ / e iξ dξ+ where R e ia e ib = ( + i sin a eia/ )e ib/, R f t (ξ) = i sin tr(ξ) f t (ξ)e it(signξ) ξ / e iξ dξ = I t ()+I t (), r(ξ) eit. Since f t decays eponentially to zero when ξ, Riemann-Lebesge lemma implies that for every t >, I t is a continos fnction of decaying to zero at infinity. On the other hand, following the analysis in Section 3 of [9] one can prove that It decays algebraically to when for t fied while for instance I () 3/ ep( iξ ) when for some non zero ξ. Althogh we will not se it, we recall for the sake of completeness a dispersive estimate derived in [5], Theorem.5, on the free Whitham grop (see also [49], Lemma.4). Note the difference with the classical L L estimates on the Airy grop. We denote S ɛ (t) = e itl( ɛd). Meso Theorem.4. There eists C > independent of ɛ sch that for any φ S(R ) the following estimates hold : ( ). S ɛ (t)φ) C ɛ /4 ( + t/ ɛ) /8 + ( + t/ ɛ) / ( φ H + φ ). 3 Note that this fnction is different from the one sed in [4] to stdy the properties of the solitary waves of the Whitham eqation.

11 apriorikdv THE WHITHAM EQUATION ( ). S ɛ (t)φ) C ɛ 3/4 ( + t/ ɛ) /3 φ L + ( + t/ ɛ) / ( φ H + φ ). ( ) 3. S ɛ (t)φ) C ɛ 3/4 ( + t/ ɛ) /3 φ L + ( + t/ ɛ) / ( φ H + φ )... A priori estimates on and v. It is well-known that the KdV eqation is well-posed in H (R). Moreover, by sing the complete integrability of KdV and in particlar the fact that KdV possesses an infinite nmber of conserved qantities, one can get global bonds at the H j -level for any j. We refer for eample to Sat [57], Bona and Smith [] and Bona, Pritchard and Scott [8] and the references therein. Proposition.5. Let φ H (R). Then there eists a niqe soltion C([, + ) : H (R)) to (.) sch that (, ) = φ. Moreover, the flow map data-soltion φ is continos from H (R) into C([, + ) : H (R)). Frthermore, the following bonds hold tre. For every j, there eists C j = C j ( φ H j) (note that the C j can be chosen to be non-increasing fnctions of their argments) sch that bondkdv (.8) (t) H j C j ( φ H j), t. aprioriwhitham Being a skew-adjoint pertrbation of the Brgers eqation, the Whitham eqation is also trivially well-posed on H (R) (and also in H s (R) for s > 3 )4 bt on a time of interval of length /ɛ. Proposition.6. Let φ H (R). Then there eist a positive time T ɛ, a niqe soltion v C([, T ] : H (R)) to (.) sch that v(, ) = φ. Moreover, the flow map data-soltion φ v is continos from H (R) into C([, T ] : H (R)). Frthermore, the following bonds hold tre. For every j N, j bondwhitham (.9) (t) H j φ H j, t ɛ. Note that the implicit constant in (.9) depends on φ H j as in (.). Proof. We only eplain how to prove (.9). Let v C([, T ] : H (R)) be a soltion of (.) with initial datm v(, ) = φ. Let J s denote the Bessel potential of order s, i.e. (J s f) (ξ) = ( + ξ ) s f(ξ). Then, it follows from the Kato-Ponce commtator estimates [33] and integrations by parts that d dt J s v L cε v L J s v L, 4 Similarly to the Brgers eqation, (see eg [4]) the Cachy problem for the Whitham eqation is epected to be ill-posed in H 3 (R) (see [7] for weaker reslts on related fractional KdV eqations).

12 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT for any s >. If s > 3, we dedce from the Sobolev embedding Hs (R) L (R) that d dt J s v L cε J s v 3 L. Hence, we dedce from a classical ODE argment that bondwhitham. (.) v(t) H s φ H s, if t (cɛ φ H s). aprioripi This finishes the proof of Proposition.6. Finally, we also need a bond on π ɛ (D) where is the soltion of (.) and π ɛ (D) is the Forier mltiplier of symbol π ɛ (ξ) defined by ( pi (.) π ɛ (ξ) = i (ξ ɛ 6 ξ3 ) l( ) ɛξ)ξ. Proposition.7. Let φ H (R) and let be the soltion of (.) evolving from φ obtained in Proposition.5. Then, for all j N, j, there eists A j = A j ( φ H j+8) > sch that aprioripi. (.) π j ɛ (D)(t) L A j ɛ, aprioripi. (.3) aprioripi.3 (.4) for all t ɛ. Proof. We apply the operator π j ɛ (D) to (.), mltiply the eqation by π j ɛ (D) integrate in the space variable over R and integrate by parts and se the Cachy-Schwarz ineqality to dedce that d dt j π ɛ (D) L = ɛ R π ɛ (D) j+ ( ) π ɛ (D) j d ɛ j+ π ɛ (D)( ) L j π ɛ (D)() L. Now, we introdce a ct-off fnction η C (R) sch that spp (η) [, ] and η = over [ /, /] and we define the Forier mltiplier P </ ɛ of symbol η( ɛ ) and P / ɛ = P </ ɛ. Then, the triangle ineqality yields j+ π ɛ (D)( ) L P </ ɛ j+ π ɛ (D)( ) L + P / ɛ j+ π ɛ (D)( ) L, and we need to estimate both terms on the right-hand side of (.4). To control the first one, we se Plancherel s identity, Lemma.3 and the fact that H s (R) is a Banach algebra for s > to dedce that P </ ɛ j+ π ɛ (D)( ) L ɛ H j+6 ɛ H, j+6 which together with Proposition.5 gives that aprioripi.4 (.5) P </ ɛ j+ π ɛ (D)( ) L ɛ C j+6 ( φ H j+6). To control the second one, we se again Plancherel s identity to obtain that P / ɛ j+ π ɛ (D)( ) L ( + ξ ) 3 ξ (j+) ( ) (ξ) dξ. ξ / ɛ

13 THE WHITHAM EQUATION 3 Now, observe that ɛ 4 ξ 8 on the spport of the integral, so that aprioripi.5 (.6) P / ɛ j+ π ɛ (D)( ) L ɛ H j+8 ɛ H ɛ C j+8 j+8 ( φ H j+8), since H 8 (R) is a Banach algebra and where we sed Proposition.5 on the last ineqality. Then, we dedce gathering (.3) (.6) that d aprioripi.6 (.7) dt j π ɛ (D) L ɛ 3 C j+8 ( φ H j+8) π j ɛ (D) L. Therefore, we dedce integrating between and t that aprioripi.7 (.8) π j ɛ (D)(t) L π j ɛ (D)φ L + C j+8 ( φ H j+8) ɛ 3 t. aprioripi.8 (.9) Finally, arging as above, we get from Lemma.3 and Plancherel s identity that j π ɛ (D)φ L P </ ɛ j π ɛ (D)φ L + P / ɛ j π ɛ (D)φ L ɛ φ H j+5 + ɛ φ H j+7. Then, we dedce combining (.8) and (.9) that aprioripi.9 (.) π j ɛ (D)(t) L A j ( φ H j+8)ɛ ( + ɛt) A j ( φ H j+8)ɛ, as soon as t t/ɛ..3. Proof of Theorem.. Let z = v and j N, j. It is dedced from eqations (.) and (.) that z solves the initial vale problem { t z + l( ɛd) z + π ɛ (D) + ɛv z ɛz = diff (.) w(, ) =, where π ɛ (D) is the Forier mltiplier of symbol π ɛ (ξ) defined in (.). Note that the eqation of z is well defined for < t /ɛ. Differentiate j times, mltiply the eqation (.) by z, j integrate in the space variable over R and integrate by parts to dedce that d ( j dt z) d = π j ɛ (D) z j d ɛ (v j z) z j d R R R maintheo.3 (.) ɛ (z j ) z j d. Ths, the Cachy-Schwarz ineqality, the Leibniz rle and integration by parts yield maintheo.4 (.3) d dt j z L π j ɛ (D) L z j L + ɛ ( L + v L ) z j L in the cases j = and j = and maintheo.4b (.4) d j dt j z L π j ɛ (D) L z j ( L + ɛ k L + v k L ) z H j R k=

14 maintheo.5 (.5) 4 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT in the cases where j. Now, on the one hand, we get by sing the Sobolev embedding and Propositions.5 and.6 that j ( k (, t) L + v(, k t) L ) k= (, t) H κ + v(, t) H κ K ( φ H κ) =: K, for all t ɛ (where the implicit constant depends on the H κ norm of φ as eplained in the proof of proposition.6). Here, we sed the notation κ = for j = or and κ = j + for j. On the other other hand, by applying Proposition.7, it follows that maintheo.6 (.6) π j ɛ (D) L ɛ A j+8 ( φ H j+8) =: ɛ K, for all t ɛ. Therefore, we dedce gathering (.3), (.5), (.6) and sing Gronwall s ineqality that e Kɛt maintheo.7 (.7) z(t) H j K ɛ K e K ɛ t, K ɛ whenever t ɛ, which concldes the proof of Theorem.. Remark.5. A similar comparison reslt can be established between the Whitham eqation with srface tension (.5) and the KdV eqation (.7) by sing the epansion ( + β ) / ( tanh ) / = 6 ( 3β) + O(4 ). 3. The Cachy problem for the capillary Whitham eqation As was already noticed the Cachy problem for the Whitham eqations (.) and (.5) is trivially well-posed in H s (R), s > 3 on time scales of order ɛ. This reslt can be improved (by enlarging the space of resoltion) for the capillary Whitham eqation (.5) by sing its dispersive properties. Actally one gets, see [43] where general fractional KdV (fkdv) eqations are considered: WST Theorem 3.. Assme that s > 6. Then, for every H s (R), there eist a positive time T ɛ = T ɛ ( H s) = (/ ɛ) (which can be chosen as a non-increasing fnction of its argment), and a niqe soltion to (.5) satisfying (, ) = sch that maintheo (3.) C([, T ɛ ] : H s (R)) and L ([, T ɛ ] : L (R)).

15 THE WHITHAM EQUATION 5 Moreover, for any < T < T ɛ, there eists a neighborhood U of in H s (R) sch that the flow map data-soltion maintheo (3.) ST s : U C([, T ]; H s (R)),, is continos. Remark 3... Since the vale α = is L critical for the fkdv eqation we conjectre that the Cachy problem for (.5) is globally well posed for initial data in the energy space H 4 (R) having a sfficiently small L norm.. This reslt was recently improved in [5] where the local well-posedness is obtained for s > For large enogh initial data, one epects a finite time blow-p silmilar to that of the L -critical generalized KdV eqation or to the cbic- Benjamin-Ono eqation proven respectively in [46] and [48], as displayed in the simlations of Section 7. Remark 3.. Using the same energy methods one can prove a reslt similar to that of Theorem. between soltions of (.5) and (.7). 4. The Cachy problem for the Bossinesq-Whitham systems In both cases (withot or with srface tension) the Bossinesq-Whitham systems are linearly well-posed since the linearized systems write, say in dimension one ( ) ( η η t + A = ) where the Forier transform of the matri operator A has real eigenvales ( ) / tanh ɛ ξ λ ± (ξ) = ± ɛ ξ in absence of srface tension and in presence of srface tension. ( ) / tanh ɛ ξ λ ± (ξ) = ±( + ɛβ ξ ) / ɛ ξ The theory of the Cachy problem for the nonlinear system withot srface tension is relatively straightforward and we focs on the case of space dimension one. Since the operator L ɛ has order zero, (.9) is an order zero pertrbation of a first order system, more precisely it writes hyp (4.) t U + A ɛ (U) U + H ɛ U = ( ( ) ( ) η ɛ ɛη L where U =, A ) ɛ (U) =, H ɛ ɛ U = ɛ. The matri A ɛ (U) has eigenvales λ satisfying (ɛ λ) = ɛη

16 6 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT so that the system t U + A ɛ (U) U = is hyperbolic in the regions where hypcond (4.) η >, and Hadamard ill-posed if this condition is violated. Actally, nder the condition hypcond (4.3) η C >, the ( system ) can be symmetrized via the positive definite symmetrizer S(U) =. ɛη The same process can be applied to (.9) yielding a symmetric hyperbolic system pertrbed by the order zero operator H ɛ namely ( ) η whitsym (4.4) S(U)U t + ɛ U η ɛη + H ɛ U =. Setting ζ = η + N where N > is a fied constraint, the standard theory of symmetrizable hyperbolic systems (see eg [45]) imply the local well-posedness of the Cachy problem for (ζ, ) H s (R) H s (R), s > 3 for initial data ζ sch that η is sfficiently small. Similar argments in the two-dimensional case yield well-posedness in H s (R ) H s (R ), s >. Note that the condition (4.) implies that the wave is always of elevation which seems to be the case of solitary wave soltions (see Section 7). On the other hand, condition (4.3) implies that the wave cannot tend to zero at infinity and ths the pertrbations of the solitary wave soltions are eclded from the range of well-posed initial data. Those facts are not physically realistic, invaliding the Bossinesq-Whitham system as a relevant water waves model. Remark 4.. We do not know of a local well-posedness reslt for (.9) or (.) nder the assmption (4.). To illstrate the (Hadamard) instabilities occrring when condition (4.3) is violated, let consider the linearization arond (η, ) = ( c, ) where c is a positive constant. The linearized system has eigenvales λ ± where ( / tanh ɛξ λ ± = ±iξ cɛ). ɛξ All modes are nstable when ɛ > /c while when ɛ < /c all modes corresponding to ξ c,ɛ are nstable where c,ɛ is the niqe positive ɛ soltion of tanh = cɛ. Note that c,ɛ as ɛ. Ths the set of stable modes get larger and larger when for a fied c, the small parameter ɛ tends to.

17 THE WHITHAM EQUATION 7 When the nonlocal term is removed, all modes are nstable, for any c >. The effect of dispersion is ths here to create a range of linearly stable modes when c <, getting larger and larger when ɛ. Things are different for Bossinesq-Whitham systems in presence of srface tension since dispersive effects play a significant role here. In fact, srface tension prevents the appearance of (Hadamard) nstable modes for the linearized system at (η, ) = ( c, ). Actally the linearized system eigenvales are now ( λ ± = ±iξ ( + βɛξ ) tanh / ɛξ cɛ). ɛξ Whatever the vales of c and ɛ the possible nstable modes are bonded and there are no more Hadamard instabilities. The effect of srface tension is ths to sppress the Hadamard instabilities. A convenient way to get an idea of the natre of the fll system is to derive an eqivalent system by diagonalizing the linear part of the system. More precisely, we define ( ( + βɛξ Â(ξ) = iξ ) tanh ) ɛ ξ ɛ ξ the Forier transform of the dispersion matri with eigenvales ±iξ( + βɛξ ) / ( tanh ɛ ξ ɛ ξ ) /. Setting ( η U = ) and ( ( ζ W = = P v) U, P = Lɛ L ɛ the linear part of the system is diagonalized as where W t + DW =, ( ) Lɛ D = L, ɛ that is two dispersive eqations of order. ), P = In the W variable the complete system writes now newfdd (4.5) t W + DW + ɛn ɛ (W ) =, ( ) Lɛ Lɛ

18 8 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT where ( ) L N ɛ (W ) = ɛ [(ζ v)l ɛ (ζ + v)) + L ɛ (ζ + v) (ζ v)] + (ζ v) (ζ v) L. ɛ [(ζ v)l ɛ (ζ + v)) + L ɛ (ζ + v) (ζ v)] (ζ v) (ζ v) While the linear part is of order, the nonlinear terms are formally of order. However the strctre of the nonlocal nonlinear terms might prevent the implementation of the methods sed in [43] for fractional type KdV eqations to obtain the local well-posedness of the Cachy problem on time scales of order / ɛ. On the other hand one might think of proving the local well-posedness by sing an elementary energy method, reminiscent of the one sed in [58] for the (,,, ) Bossinesq system (see [58] Theorem 3.) or a symmetrization techniqe as in [58] Theorem 4.5. However those methods do not seem to work since a control of the H s H s+/, s > norm on (η, ) is not enogh to control the nonlinear terms (in case of the (,,, ) Bossinesq system one gets a control on the H s H s+ norm). 5. Finite time blow-p for the Whitham eqation It has been established in [6] that soltions of the Whitham eqation (.) may blow-p in finite time, the blow-p being "shock-like", that is blow-p of the gradient with bonded L norm. The nmerical simlations in [37] display the csp natre of the singlarity depending on the sign of the initial data. On the other hand those works consider the Whitham eqation (.) with ɛ =, ths they do not address the long-wave, KdV limit. In particlar the crcial dependence of the blow-p time with respect to ɛ is not obtained. Note that this blow-p time shold be larger that /ɛ, the time scale on which the KdV and Whitham eqations are close and asymptotic models for the propagation of long, weakly nonlinear srface water waves for which no singlarities are epected. Note also that a qite different finite time blow-p is epected for the Whitham eqation with srface tension (.5) that cannot be anymore viewed as a "weak" dispersive pertrbation of the Brgers eqation. In fact, as was already noticed, it is reminiscent for large freqencies of the L critical fractional KdV eqation cfkdv (5.) t + + D = for which a finite time blow-p à la Martel-Merle [46, 48] is epected (bt not yet proven), see the nmerical simlations in [35, 36, 37]. 6. Solitary waves 6.. Whitham withot srface tension. We will focs here on non periodic solitary wave soltions of the Whitham eqation (.), that is soltions of the form ( ct).

19 THE WHITHAM EQUATION 9 We refer to [,, 3, 8, 55] for interesting (theoretical and nmerical) stdies on periodic traveling waves. In particlar the eistence of a global bifrcation branch of π periodic smooth traveling wave soltions is established in [3]. The eistence of solitary wave soltions to the Whitham eqation (.) decaying to zero to infinity and close to the KdV soliton has been proven in [] by variational methods. Symmetry and decay properties of Whitham solitary waves are established in [3]. The analysis in [] is made on (.) with ɛ = nder the scaling () = ɛ α w(ɛ β ) where α β = so that R = ɛ. An interesting isse is that of the transverse stability of the Whitham solitary wave in the framework of the fll dispersion KP eqation (.3). Considering the similar problem for the sal KP I/II eqations (see [56] and the references therein) one can conjectre that the Whitham solitary wave is transversally stable when the srface tension parameter β is less than 3 and nstable otherwise. We intend to go back to this isse in a sbseqent paper. We do not know of any rigoros reslt on the eistence of solitary wave soltions to the fll dispersion Bossinesq system. The following comptations give some evidence to the eistence of solitary wave soltions "close" to the KdV soliton. Recall that the system reads η t + T ɛ + ɛ(η) =, sys (6.) t + η + ɛ ( ) =, where ˆT ɛ = (tanh(k ɛ)/(k ɛ)) /. In the formal limit ɛ, the system redces to the wave eqation. We are here interested in localized travelling wave soltions to the system (6.). Ptting (, t) = U( ct) and η(, t) = N( ct), we find after integration cn + T ɛ U + ɛnu =, travelsys (6.) N = cu ɛ U. Eliminating N from the first eqation of (6.) via the second, we find travel (6.3) (T ɛ c )U + 3ɛc U ɛ U 3 =.

20 travele (6.4) C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT Writing c = + αɛ with α >, we get, by performing a formal epansion in ɛ and neglecting terms of order ɛ or higher in (6.3) which gives after integration 3 U αu + 3 U =, traveleint (6.5) (U ) = 6αU 3U 3 which has the KdV soliton ( ) 3α sol (6.6) U = αsech ( ct) wtravel (6.7) as a soltion. Note that for a soltion (N, U) with a given velocity c of the above eqations, ( N, U) provides a soltion for the same system with velocity c. Similarly, a solitary wave (, t) = U( ct) of the Whitham eqation (.) satisfies the eqation ɛ U + (T ɛ c)u =. Remark 6.. Following the method sed in [3] one can prove that soltions of (6.7) that tend to as decay to eponentially in the sense that for some ν >, e ν φ L (R) L (R). Writing c = + δɛ where δ > is a constant independent of ɛ, we get, neglecting terms of order ɛ or higher in the formal epansion in ɛ, wtravel (6.8) U 6 U + δu =, which gives after integration wtravel3 (6.9) (U ) = U 3 + 6δU. Ths we find again the KdV soliton ( ) 3δ wtravel4 (6.) U = 3δsech ( ct). Remark 6.. The above formal considerations sggest that one cold obtain the eistence of solitary wave soltions to both the Whitham eqation and the Whitham system by pertrbations argments starting from the KdV soliton. In particlar in the case of the Whitham eqation this wold give an alternative proof to the one in [].

21 SWST (6.) THE WHITHAM EQUATION 6.. Whitham with srface tension. The sitation is qite different for the Whitham eqation with srface tension (.5). A solitary wave soltion (, t) = U( ct) of (.5) shold satisfy the eqation ɛ U + ( T ɛ c)u = where T ɛ is defined by (.6). The eistence of localized non trivial soltions of (6.) reslts from Theorem. in [6]. Since (.5) (as the fkdv eqation with α = ) is L critical one epects the instability of those solitary waves by blow-p, similarly to the generalized KdV eqation with p = 4, see [46] or to the modified Benjamin- Ono eqation ([48]). On the other hand, a solitary wave soltion (η( ct), ( ct)) of the fll-dispersion Bossinesq system shold satisfy SWFDd (6.) { cη + P + ɛη = c + η ɛ =, ( where P ɛ = (I + βɛ D ) ) tanh( ɛ D ) ɛ D. Eliminating η from the second eqation in (6.) yields the single eqation SWFD (6.3) c + P ɛ + 3 c ɛ 3 =. One cannot apply directly the reslt of [6] to prove the eistence of non trivial soltions to (6.3) since the nonlinear term is not homogeneos.we plan to come back to this isse in a net paper. 7. Nmerical simlations for the Whitham eqation The simlations in this section will illstrate several aspects of the dynamics of the Whitham eqation allded to in the previos sections and will also give evidence for some facts not yet rigorosly proven. 7.. Nmerical constrction of solitary waves for the Whitham eqation. In this sbsection we constrct nmerically solitary wave soltions to the Whitham eqations, i.e., localized soltions of (6.7). To this end we se the same techniqe as in [37]: since the wanted soltions are epected to be rapidly decreasing, it is convenient to formlate the problem as an essentially periodic one with a sfficiently large period sch that the soltions and the nmerically interesting derivatives vanish with nmerical accracy at the domain bondaries ±Lπ (we typically choose L = 5 in the following). The Forier transforms are approimated in standard way via

22 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT discrete Forier transforms conveniently implemented via the Fast Forier transform (FFT). With this approach eqation (6.7) is approimated via a finite dimensional (we se N Forier modes) nonlinear eqation system formally written as NK (7.) F(Û) =, where Û denotes the discrete Forier transform of U. This nonlinear system is solved with a Newton-Krylov iteration. This means the action of the inverse of the Jacobian of F on F in the standard Newton iteration is determined iteratively via GMRES [54]. As the initial iterate we choose the KdV ( soliton written as U = 3δsech 3δ ). ( ct) Note that the reality of U has to be enforced dring the iteration. As an eample we stdy the case ɛ =. and vary the constant δ. This is mainly done in order to have for small δ an initial iterate close to the Whitham soliton. For larger δ, we se as an initial iterate the soltion for a slightly smaller vale of δ constrcted before. Since we vary δ, the choice of ɛ is, however, no restriction for the general case since the only important qantity is the velocity c. It is to be epected that there will be solitons for c becase of the analogy to KdV, bt that for larger vales of c the dispersion of the Whitham eqation becomes in contrast to the KdV eqation too weak to compensate the nonlinearity. Ths we epect that in contrast to KdV there might be an pper limit of c for solitons to the Whitham eqation. In fact we obtain the sitation shown in Fig.. For c =.,.,.5,.,. we se N = 4 Forier modes, for c =. N = QU 4 QU Figre. Solitary waves for the Whitham eqation (6.7): on the left the Whitham solitons for c =.,.,.5,.,. (in the order of growing maima), on the right the Whitham soliton for c =. in ble and the corresponding KdV soliton in red. figwhitsol It can be seen that maima of the soltion grow as epected, bt that the solitons become more peaked and more compressed compared to the corresponding KdV case. In fact we did not scceed to constrct solitons for

23 THE WHITHAM EQUATION 3 c mch greater than.. The failre of an iteration to converge obviosly does not imply that there will be no longer a soliton, it jst means that the nmerical approach no longer can be sed. However the fact that the iteration did not converge even for larger resoltions, for initial iterates being the nmerical soltion for a slightly smaller δ and for an iteration with relaation is a strong indication that there might not be a Whitham solitary waves for mch larger velocities. There are additional differences between Whitham and KdV solitons. The mass (the L norm) of the KdV soliton (6.) is given by M = 4 6δ 3/, the mass of the Whitham solitons can be seen on the left of Fig.. As epected the mass is identical for small δ to the one of the KdV soliton. Bt for larger speeds c, the mass grows more slowly. The energy of the KdV soliton (6.) is proportional to δ 5/, the energy of the Whitham soliton can be seen on the right of Fig.. The energy crve appears to flatten near the reachable maimal velocities M 6 E c c Figre. Mass of the solitary waves for the Whitham eqation (6.7) in dependence of c (the continos line gives the mass of the KdV solitons (6.)) on the left, and the corresponding energy on the right. figwhitsolmass 7.. Nmerical stdy of pertrbed Whitham solitons. In this sbsection, we stdy the time evoltion of pertrbations of the Whitham solitons constrcted above. To this end we write the Whitham eqation for the Forier transform in a commoving frame, i.e., whithamforier (7.) û t + ɛ ikû + ( ˆT ɛ c)ikû =, and approimate the Forier transform via a discrete Forier transform as before. For the time integration we se as in [37] an implicit Rnge-Ktta method of forth order with a fied point iteration, see [37] for details. The accracy of the soltion is controlled via the Forier coefficients which shold decrease eponentially to machine precision (we work here in doble precision which allows a maimal accracy of 6 ) for smooth fnctions

24 4 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT and via the conservation of the energy [ whithame (7.3) E = R ( T ɛ c) ɛ ] 6 3 d. De to navoidable nmerical errors, this energy does nmerically evolve with time, bt ths provides an estimate of the nmerical error. As shown in [34], conserved qantities typically overestimate the L error (the maimm of the difference between nmerical and eact soltion) by - orders of magnitde. As a first test of the code, we take the nmerically constrcted soliton with ɛ =. and c =. of the previos sbsection as initial data. In the sed commoving frame, the solitons shold correspond to a stationary soltion. We se N t = 4 time steps for t [, ] and find that the difference between the evolved soltion and the initial data is of the order of 3, the order of the error with which the eqation (6.7) had been solved in the previos section. The nmerically compted relative energy E(t)/E() is of the order of 5.5. We first consider as initial data the soliton with ɛ =. c =. pls a Gassian pertrbation which is a pertrbation of the order of % of the maimm of the soliton. We se N t = 4 time steps for t [, ]. The Forier coefficients in this case decrease to the order of 4 dring the whole comptation, and the relative energy is conserved to the order better than 4. In this case there appears a slightly faster soliton travelling to the right in the commoving frame, as can be seen in Fig. 3. Bt in addition to radiation travelling to the left, a slower soliton travelling to the left in the frame commoving with the npertrbed soliton seems to emanate from the initial data. This is even more visible in the plot on the right of Fig. 3 where a close p of the smaller soliton is shown for t =. The soliton appears to be stable, the pertrbation leads to a slightly larger soliton and possibly a smaller soliton pls radiation. The sitation changes if the same soliton as in Fig. 3 is considered, bt this time with a pertrbation 3 ep( ) corresponding to roghly 4% of the maimm of the soliton. We se N t = 4 time steps for t [, 8]. In this case the soltion appears to develop a csp. We show in Fig. 4 the soltion and the modls of its Forier coefficients at t = 7.5. It can be seen that the soltion clearly rns ot of resoltion in Forier space, already at t 7, the Forier coefficients decrease only to. In Fig. 5 it can be seen that the L norm of the soltion grows clearly beyond the initial data, bt it is not clear whether it diverges. The same behavior can be seen for the L norm of on the right of the same figre. A fit of the Forier coefficients as in [59, 37] reveals that a singlarity appears to approach the real ais in the comple -plane. The idea of this approach is that an essential singlarity of the form ( s ) µ, µ / Z, in the comple plane

25 THE WHITHAM EQUATION Figre 3. Soltion to the Whitham eqation with =. for initial data the soliton with c =. pls a Gassian pertrbation; on the left the soltion in dependence of time, on the right a close p of the soltion at the final time t =. figwhitsolcgass log k Figre 4. Soltion to the Whitham eqation with =. for initial data the soliton with c =. pls a pertrbation 3 ep( ); on the left the soltion for t = 7.5, on the right the modls of the Forier coefficients. for a fnction leads to an asymptotic behavior of the Forier transform fit (7.4) (k) k µ+ e δk, k, where δ = Ims. For t = 7.547, this singlarity seems to hit the real ais (δ ) indicating the formation of a csp (the fitted coefficient µ of the essential singlarity is positive as in [37]). Ths there seems to be a hyperbolic blow-p in this case. The soliton for c =. appears therefore to be nstable against blow-p in the form of a csp for sfficient size of the pertrbation. Since the solitons of the Whitham eqation are rapidly decreasing, it is possible to stdy soliton interactions. In Fig. 6, we show the solitons with c =.5 and c =., the latter being centered at = 4 by simply mltiplying its Forier coefficients by e4ik. The sm of the soltions gives two soliton initial data since they both vanish with nmerical precision where figwhitsolc3gass

26 6 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT t t Figre 5. Norms of the soltion to the Whitham eqation with ɛ =. for initial data the soliton with c =. pls a pertrbation 3 ep( ); on the left the L norm, on the right the L norm of normalized to for t =. figwhitsolc3gassnorm the other takes vales above the nmerical error. On the left of Fig. 6 we show the soltion to the Whitham eqation for t [, ]. The modls of the Forier coefficients of the soltion decreases to dring the whole comptation. It can be seen that the soliton interaction resembles the KdV two-soliton: the solitons have almost the same shape after the collision, there is jst a phase shift (note that the comptation is performed in a frame commoving with the smaller soliton at c =.5). The close-p of the soltion at the final time on the right of Fig. 6 reveals, however, that this is not an eact two-soliton since there is dispersive radiation propagating to the left. This indicates that the Whitham eqation is as epected not integrable Figre 6. Soltion to the Whitham eqation for ɛ =. for initial data being the sm of the solitons with c =.5 centered at = and for c =. centered at = 4 in dependence of time on the left and a close-p of the soltion at t = on the right. figsatwhithamsol 7.3. Nmerical stdy of Gassian initial data for different vales of ɛ. In this section we stdy nmerically the time evoltion of Gassian

27 THE WHITHAM EQUATION initial data (, ) = ep( ) for different vales of the parameter. Since we are interested in stdying the soltions on time scales of order /, we introdce the rescaled time τ = t in which the Whitham eqation (.) reads (7.5) τ + + (T ) =, where we have sed a frame commoving with velocity. This eqation is solved with the same nmerical approach as in the previos sbsection. For = we se N = 4 Forier modes and Nt = 4 time steps for τ [,.]. The soltion appears to develop a csp for τ >.75 as can be seen in Fig. 7 on the left. The Forier coefficients on the right of the same figre indicate a loss of resoltion. A fitting of the Forier coefficients to (7.4) indicates that for τc =.75, the parameter δ indicates the distance between a singlarity in the comple -plane vanishes, and the corresponding parameter µ.3635 indicating the formation of a csp proportional to ( s )/ log 6 whithamta k Figre 7. Soltion to the Whitham eqation with = for initial data (, ) = ep( ); on the left the soltion for τ =.75, on the right the modls of the corresponding Forier coefficients. Note that the critical time τc above does not indicate the eact blow-p time, this jst indicates that the singlarity in the comple plane leading to the asymptotic behavior (7.4) of the Forier coefficients is too close to the ais to be nmerically distingished from. If one considers the same initial data as in Fig. 7 for smaller vales of, one obtains the critical times and eponents shown in the table in Fig.. It can be seen that the critical times τc are always greater than the τc for =. In fact the critical times appear to grow with decreasing. This implies, however, that the soltions do not have a blow-p on time scales of order / (recall that τ = t ). For vales of smaller than.4, there does not appear to be a blow-p at all. In Fig. 8 we show the soltion to the same initial data as in Fig. 7 for =.. We se N = Forier modes for L[ π, π] with L = and Nt = 4 time steps for τ [, ]. It can be seen that two stable solitons figwhitgasse

28 8 C. KLEIN, F. LINARES, D. PILOD, AND J.-C. SAUT figwhitgasstable ɛ τ c µ Table. Critical times τ c = ɛt c and eponents according to (7.4) for the soltions to the Whitham eqation for initial data (, ) = ep( ) for different vales of ɛ. appear which gives spport to the soliton resoltion conjectre in the contet of the Whitham eqation: stable solitons seem to appear in the long time behavior of the soltion. In addition there is the sal dispersive radiation also known from KdV soltions. The L norm of the soltion on the right of Fig. 8 also seems to correspond to a soliton with speed greater than (we are in a commoving frame with c = ). Neither the Forier coefficients nor the L norm of indicate the formation of a singlarity. The latter again indicates the appearance of a soliton t Figre 8. Soltion to the Whitham eqation with ɛ =. for initial data (, ) = ep( ); on the left the soltion for τ =, on the right the L norm of the soltion in dependence of time. figwhitgasse The soltion in Fig.8 for small ɛ becomes obviosly closer to the soltion to the KdV eqation τ + + /6, which can be obtained from the Whitham eqation in the formal limit ɛ for the same initial data. This KdV soltion is shown in Fig Nmerical stdy of the Whitham eqation with srface tension. In this sbsection we will stdy the Whitham eqation with srface