Access the full text.
Sign up today, get DeepDyve free for 14 days.
Jean-Paul Chehab LAMFA, Université de Picardie Jules Verne, France Denys Dutykh CNRS–LAMA, Université Savoie Mont Blanc, France On time relaxed schemes and formulations for dispersive wave equations arXiv.org / hal arXiv:1903.02212v1 [physics.class-ph] 6 Mar 2019 Last modified: March 7, 2019 On time relaxed schemes and formulations for dispersive wave equations Jean-Paul Chehab and Denys Dutykh Abstract. The numerical simulation of nonlinear dispersive waves is a central research topic of many investigations in the nonlinear wave community. Simple and robust solvers are needed for numerical studies of water waves as well. The main diﬃculties arise in the numerical approximation of high order derivatives and in severe stability restrictions on the time step, when explicit schemes are used. In this study we propose new relaxed system formulations which approximate the initial dispersive wave equation. However, the resulting relaxed system involves ﬁrst order derivatives only and it is written in the form of an evolution problem. Thus, many standard methods can be applied to solve the relaxed problem numerically. In this article we illustrate the application of the new relaxed scheme on the classical Korteweg–de Vries equation as a prototype of stiﬀ dispersive PDEs. Key words and phrases: dispersive wave equations; shallow water ﬂows; relaxation; quasi-compressibility. MSC: [2010] 74J15 (primary), 74S10, 74J30 (secondary) PACS: [2010] 47.35.Bb (primary), 47.35.Fg, 47.85.Dh (secondary) Key words and phrases. dispersive wave equations; shallow water ﬂows; relaxation; quasi- compressibility. Corresponding author. Contents 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2 Mathematical models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1 Korteweg–de Vries equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 3 Properties . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 3.1 Dispersion relation analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 4 Numerical methods and results . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 4.1 Linear stability analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 4.2 Enhanced stability analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 4.3 Simulation of a soliton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 4.4 Second order relaxed times schemes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 5 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 5.1 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 5.2 Perspectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 A Other types of dispersive equations . . . . . . . . . . . . . . . . . . . . . . . 22 A.1 Benjamin–Bona–Mahony equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 A.2 Peregrine system . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 J.-P. Chehab & D. Dutykh 4 / 26 1. Introduction Nonlinear dispersive waves arise in various ﬁelds of science including the solid [29] and ﬂuid mechanics [20]. Below we focus on models stemming from the free surface hydrody- namics mainly due to scientiﬁc interests of the authors. In this ﬁeld, analytical solutions (such as solitary or cnoidal waves [8, 10]) are available only for some simpliﬁed models. That is why the numerical simulation remains one of the main tools in nonlinear dispersive wave studies. In this work the authors have an ambition to propose a novel strategy to tackle numerically these problems. Let us review ﬁrst the modern numerical approaches to this problem. Historically, the ﬁnite diﬀerences were applied to dispersive wave problems [19]. Today, the trend is to apply Galerkin–type methods (e.g. FEM) in smooth situations [1, 26] and ﬁnite volumes [14, 15] or coupled ﬁnite volumes (for the hyperbolic part)/ﬁnite diﬀerences (for the dispersive part) via the operator splitting [6]. Notice that there is an eﬀort to develop discontinuous Galerkin–type methods for dispersive wave equations as well [16, 36], which combine the robustness of ﬁnite volumes with the accuracy of ﬁnite elements. However, the ‘price to pay’ (i.e. the CPU time) remains excessively high turning these methods into a luxury limousine rather than a working horse of numerical analysis. Finally, for periodic situations one can apply highly eﬃcient Fourier-type pseudo-spectral methods [13]. On the borderline between ﬁnite diﬀerences and pseudo-spectral methods there exist compact ﬁnite diﬀerence schemes proposed by Lele (1992) with spectral-like resolution. To give an example, these schemes were applied to the Serre equations in [9]. In the present study we propose to change the numerical strategy in contrast to studies described above. The main idea consists in modifying the governing equations by slightly perturbing them with some ad-hoc terms. In this way, we gain the structure suitable for numerical simulations and we hope that solutions of the perturbed and un-perturbed prob- lems will remain close provided that perturbation parameter is small. This idea is inspired by pseudo-compressibility methods proposed for the numerical simulation of incompress- ible Navier–Stokes equations (see e.g. [21]). The main philosophy of this approach is as follows. When we solve numerically a diﬀerential equation, some errors are introduced due to the underlying discretization process. So, perhaps, for the sake of convenience, one can perturb also the governing equations within the same discretization error which is already introduced. And if everything is done judiciously, the end user will not even see the diﬀerence between the relaxed and the original problems solutions. It makes a lot of ‘ifs’, but nevertheless we undertake this programme below for some well-known dispersive wave equations. We would like to mention here also another relaxation technique based on the application of local time-averaging operators [2]. A relaxation scheme for the KdV nd equation was proposed in [4]. However, their formulation contains 2 order derivatives. Below, we will propose an alternative formulation which has an advantage to involve only Relaxed formulations for dispersive wave equations 5 / 26 st ∗ the 1 order derivatives. The idea we employ can be traced back at least to various local formulations used in continuous [33, 34] and discontinuous [24, 36] Galerkin methods. However, we push it one step further towards the so-called relaxed local formulations. The present study is organized as follows. In Section 2 we propose the relaxed formu- lations for several well-known dispersive wave models. Some mathematical properties of a relaxed formulation for the KdV equation are discussed in Section 3. Our numerical approach based on this relaxed formulation of the KdV equation is outlined in Section 4. The ﬁrst numerical results are presented in the same Section. Finally, in Section 5 we outline the main conclusions and perspectives of our study. The relaxation of some other dispersive wave equations is discussed in Appendix A. 2. Mathematical models Instead of working out the most general situation, in the present study we follow the so- called Gelfand principle which states that a theory should be illustrated on the simplest non-trivial example. Below we consider one such example and two others are treated in Appendix A. All cases which steem from the water wave theory [32]. For each model we develop the corresponding relaxed formulation. 2.1. Korteweg–de Vries equation The celebrated Korteweg–de Vries (KdV) equation was derived independently by th J. Boussinesq [7] and D. Korteweg & G. de Vries [22] at the end of the XIX century and in the scaled form it reads: u + u u + u = 0 . (2.1) t x xxx The variable u (x, t) may be interpreted physically as the free surface elevation or the horizontal depth-integrated ﬂuid velocity. There exist transformations in order to obtain the canonical form (4.1) [20]. We shall rewrite the last equation in the conservative form in order to obtain the ﬂux of the evolution variable u : 1 2 u + u + u = 0 . (2.2) t xx 2 x In the case of continuous Galerkin methods this formulation is also sometimes called the modiﬁed one. Here we mean I. M. Gelfand to avoid any confusion. J.-P. Chehab & D. Dutykh 6 / 26 Now, the order of derivatives has to be lowered. It can be done similarly to the ODE theory by introducing additional variables: 1 2 u + u + w = 0 , 2 x u = v , v = w . At this step we are compatible with local formulations used in the framework of discontinu- ous Galerkin methods [24, 36]. It is obvious that the last system is completely equivalent to equation (2.2) (at least for the smooth solutions). However, the last system is not of the evolution type. It is here that the relaxation comes into the play in order to correct this shortcoming: 1 2 u + u + w = 0 , (2.3) 2 x δ v + v − u = 0 , (2.4) t x δ w + w − v = 0 , (2.5) t x where δ ≪ 1 is a small parameter. It can be seen that in the limit δ → 0 we recover the original formulation (4.1) (at least formally). Equations (2.4), (2.5) allow also to assess the physical meaning of the parameter δ . Indeed, the terms v and δ v (and w with δ w ) need to have the same units. Thus, the t t multiplication by δ has to compensate the derivative with respect to time. Henceforth, from dimensional arguments, δ has the unit of time. 3. Properties Above we presented relaxed formulations for three well-known equations. However, the properties of new formulations have to be studied. From now on we focus on the KdV equation only (2.3) – (2.5) in order to understand the relaxation eﬀects exhaustively. Other equations will be tackled in following studies. 3.1. Dispersion relation analysis In order to study the dispersive properties, ﬁrst we have to linearize equations (2.3) – (2.5): u + w = 0 , (3.1) t x δ v + v − u = 0 , (3.2) t x δ w + w − v = 0 . (3.3) t x Now, we look for plane-wave solutions of the form i (k x − ω t) i (k x − ω t) i (k x − ω t) u (x, t) = u · e , v (x, t) = v · e , w (x, t) = w · e , 0 0 0 Relaxed formulations for dispersive wave equations 7 / 26 where k is the wavenumber, ω is the frequency and u , v , w ∈ R are some real 0 0 0 amplitudes. By substituting the plane-wave ansatzä into a linear system (3.1) – (3.3), we obtain a system of linear algebraic equations, where unknowns are u , v , w : 0 0 0 −i ω u + i k w = 0 , 0 0 −i δ ω v + v − i k u = 0 , 0 0 0 −i δ ω w + w − i k v = 0 . 0 0 0 The necessary condition to have non-trivial solutions to the above linear system of equations is that its determinant is equal to zero: − i ω 0 i k det = 0 . − i k 1 − i δ ω 0 0 − i k 1 − i δ ω After some simpliﬁcations we obtain the desired dispersion relation between ω and k : 3 2 2 3 ω + k − 2 i δ ω − δ ω = 0 . (3.4) So, one can see that by taking the limit δ → 0 in the last equation, we recover the classical KdV dispersion relation ω = − k . This polynomial equation (3.4) in wave frequency ω admits in general three complex roots. Their general expression is rather cumbersome. However, we can construct their asymptotic expansions: 3 6 9 2 3 ω (k) = − k + 2 i δ k + 7 k δ + O (δ ) , (3.5) 3/2 √ √ i k 1 3 5 9/2 6 3/2 ω (k) = − ± i − i + k ± − i k δ − i k δ + O (δ ) . (3.6) 2, 3 2 8 These expansions (3.5), (3.6) might be used to interpret the behaviour of numerical solu- tions. Since we deal with singular perturbations in equation (3.4), the passage to the limit δ → 0 in solutions ω (k) is far from being trivial. However, from the expansion (3.5) 2, 3 we learn that parameter δ plays the both rôles of dissipation and dispersion since it enters into real and imaginary parts of the regular branch ω (k) . We would like to highlight the fact that the relaxation parameter δ aﬀects the dispersive properties of the regular branch (i.e. Re ω (k)) only to the second order in δ . Thus, we can conclude that the proposed relaxation approach is very careful with respect to the dispersion relation of the original KdV equation. Remark 1. Taking the limit in the polynomial equation (3.4) turns out to be an easy task, as it often happens, than taking the limit δ → 0 in the solutions ω (k) . On the other 2, 3 hand, the branch ω (k) seems to be completely regular from this point of view. Two other branches ω (k) are artiﬁcial modes introduced by relaxation. 2, 3 This operation can be seen as the Fourier transform in both space and time. J.-P. Chehab & D. Dutykh 8 / 26 4. Numerical methods and results In this Section we are going to study numerically the relaxed formulations. For this purpose we propose three diﬀerent discretisations of the KdV equation. One of them is based on the relaxed formulation proposed above and two others are popular ﬁnite diﬀer- ence schemes (Crank–Nicolson and Sanz–Serna) based on the original KdV equation. The numerical advantages of the relaxed formulation are illustrated below. Before considering the time schemes, we focus on the discretization in space which is realized with ﬁnite diﬀerences. The use of compact schemes as presented in [23] allows to reach a high order of accuracy, with a spectral-like resolution, while preserving the capabilities of the ﬁnite diﬀerences approach, including the implementation on a general Cartesian grids and for a variety of boundary conditions. When dealing with surface wave equations, the high level of accuracy is an important property to capture correctly the active frequencies of the solution without working on very ﬁne grids. Hence, the computation can be done with a reasonable computational time. At ﬁrst we recall brieﬂy the principle of the Compact Schemes (CS) and we give the discretization matrices that will be used for the simulation of the models presented in Section 2. In two words, the compact schemes consist in approaching a linear operator (diﬀerenti- ation as well as interpolation) by a rational (instead of polynomial-like) ﬁnite diﬀerences scheme. These schemes are then implicit, this allows to increase the accuracy and mimic the spectral global dependence, see [23] for more details. The approximations at grid points of the operators applied to a regular function u is realized as follows. def Let U := (U , . . . , U ) be a vector whose the components are the approximations u 1 N at (regularly spaced) grid points x = i h , i = 1, . . . , N , here h = is the spatial step-size. We compute approximations of V = L (u) (x ) as solution of a system i i P · V = Q U , −1 the approximation matrix is then formally B = P · Q . For simplicity we present here fourth order accurate CS. When considering periodic boundary conditions, the matrices P and Q are for the ﬁrst order derivative in space: 3 3 0 − 1 1 2 2 4 4 3 3 − 0 2 2 def . def 1 . . . 4 . . . P := , Q := , . . . . . . . 2 h . . 3 3 − 0 2 2 1 1 3 3 4 4 − 0 2 2 Relaxed formulations for dispersive wave equations 9 / 26 For the second order derivative in space we have the following pair of matrices: 12 6 6 − − 5 5 5 1 1 6 12 6 1 − − 10 10 5 5 5 1 1 6 12 6 1 − − 10 10 5 5 5 . . . . . . . . . . . . P = , Q = . . . . . . . 2 1 1 6 12 6 1 − − 10 10 5 5 5 1 1 6 12 6 1 − − 10 10 5 5 5 6 6 12 − − 5 5 5 Finally, for the third order derivative in space we have 1 1 1 1 1 0 −1 − 1 2 2 2 2 1 1 1 1 1 1 0 −1 − 2 2 2 2 . . . . . . . . . . . . P = , Q = . . . . . . . 3 1 1 1 1 1 − 1 0 −1 2 2 2 2 1 1 1 1 1 −1 − 1 0 2 2 2 2 We now present the semi-implicit relaxed time scheme we will use for the simulation, stressing out its enhanced stability as compared to classical semi-implicit schemes (Crank– Nicolson’s for the linear terms and forward Euler’s for the nonlinear ones) while giving comparable results to the ones obtained by fully nonlinear Sanz–Serna’s which is second order accurate in time and unconditionally stable. Let D , D and D be the discretization matrices of the ﬁrst, second and third order x xx xxx derivative on a N regularly spaced grid points, with periodic boundary conditions; these matrices will be obtained in practice with compact schemes as presented above; I will denote the N × N identity matrix. The relaxed scheme produces the iterations (k+1) (k) U − U 1 (k+1) (k) (k) 2 + D (W + W ) + (U ) = 0 , Δt 2 (k+1) (k) V − V 1 (k+1) (k) (k+1) (k) δ + (V + V ) − D (U + U ) = 0 , Δt 2 (k+1) (k) W − W 1 (k+1) (k) (k+1) (k) δ + (W + W ) − D (V + V ) = 0 . Δt 2 We rewrite this coupled system by blocs as Δt (k+1 I 0 U ) Δt Δt (k+1 · = − D (δ + ) I 0 V ) x N 2 2 Δt Δt (k+1 0 − D (δ + ) I W ) x N 2 2 Δt (k) Δt (k) 2 I 0 − U − D (U ) N x 2 2 Δt Δt (k) · + . D (δ − ) I 0 V 0 x N 2 2 Δt Δt (k) 0 D (δ − ) I W 0 x N 2 2 J.-P. Chehab & D. Dutykh 10 / 26 We can now introduce the iteration matrix def i −1 e M := (M ) · M = r r −1 Δt Δt I 0 I 0 − N N 2 2 Δt Δt Δt Δt · . − D (δ + ) I 0 D (δ − ) I 0 x N x N 2 2 2 2 Δt Δt Δt Δt 0 − D (δ + ) I 0 D (δ − ) I x N x N 2 2 2 2 We can resume the relaxed scheme as follows Algorithm 1 Relaxed Scheme for the KdV equation. (0) 1: U is given (0) (0) (0) (0) 2: Set V = D U , W = D V x x (0) (0) (0) (0) 3: Set Z = U , V , W 4: for k = 0, 1, . . . do Δt (k) (k) 2 5: Set F = − D (U ) , 0, 0 i (k + 1) e (k) (k) 6: Solve M · Z = M · Z + F r r (k+1) (k + 1) (k + 1) (k + 1) 7: Set U = Z (1 : N) , V = Z (N + 1 : 2 N) , (k + 1) (k + 1) W = Z (2 N + 1 : 3 N) 8: end for We now present the two references schemes to which we will compare the relaxed one. We set def Δt def Δt i e M := I + D , M := I − D N x x x N x x x CN CN 2 2 and we deﬁne the corresponding iteration matrix M as CN def i −1 e M := (M ) · M . CN CN CN Algorithm 2 Fully nonlinear Sanz–Serna’s scheme. (0) 1: U is given 2: 3: for k = 0, 1, . . . do (k + 1) (k) i (k + 1) e (k) Δt U + U 4: Solve M U = M U − D CN CN 2 2 5: end for We will use also 4.1. Linear stability analysis Before comparing the stability properties of the classical CN and of the relaxed schemes, we give hereafter a simple but instructive result illustrating the advantage of the compact schemes in a context in which the spectral properties of the operators must be restituted Relaxed formulations for dispersive wave equations 11 / 26 Algorithm 3 Semi-implicit Crank–Nicolson scheme. (0) 1: U is given 2: 3: for k = 0, 1, . . . do i (k + 1) e (k) Δt (k) 4: Solve M U = M U − D U CN CN 5: end for at the discrete level. Writing the CN scheme for the Airy equation, we have the following th induction relation among Fourier coeﬃcients attached to the m frequency (k + 1) (k) uˆ − uˆ i m m 3 (k + 1) (k) − m uˆ + uˆ = 0 , m m Δt 2 hence Δt 1 + i m (k + 1) 2 (k) uˆ = uˆ . m m Δt 1 − i m Δt 1 + i m The Fourier symbol of the CN operator is and its values are all displayed Δt 1 − i m on the unit circle, for any Δt > 0 . It is then an important property to be captured by the spectrum of the Crank–Nicolson matrix, when dealing with ﬁnite diﬀerences. We nd give hereafter in Figure 1 the spectra of the CN matrix when the 2 order FD and the th 4 order compact FD are used. We observe that the compact scheme allows to capture numerically the spectral properties of the linear operator as in the analytical Fourier nd analysis, while the 2 order ﬁnite diﬀerences succeed to capture correctly only a small th part of the spectrum. Comparable results are obtained in Figure 2 when considering 6 order compact FD: they illustrate the capabilities of the compact schemes to mimmic the spectral properties of the linear propagation operator. 4.2. Enhanced stability analysis A ﬁrst way to understand the eﬀect of the relaxation in the stability of the time marching th schemes is to consider the linear system (3.1). We have used the 6 order compact schemes for the discretization of the spatial derivatives, as presented in the beginning of Section 4. We compare here in the complex plane the eigenvalues of the two matrices used in the numerical schemes: M and M . We recall that these matrices are not of the same r CN size so there is no reason that their relative spectra coincide. The eigenvalues of M are CN mapped on the unit circle (see Figure 1), we present hereafter comparison with those of M for diﬀerent values of N , Δt and δ (the relaxation parameter). At ﬁrst, we take a very −17 small value of δ = 10 , however the results we obtain are identical all for larger but J.-P. Chehab & D. Dutykh 12 / 26 Eigenvalues Iter. Mat. Fourier Symbol of the CN Operator 1 1 CN with 2d order FD Fourier symbol of CN CN with 4th order FD Unit circle 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.5 0 0.5 1 1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −3 −3 (a) Δt = 10 CN Matrices spectra (b) Δt = 10 Fourier symbol Eigenvalues Iter. Mat. Fourier Symbol of the CN Operator 1 1 CN with 2d order FD Fourier symbol of CN CN with 4th order FD Unit circle 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.5 0 0.5 1 1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −2 −2 (c) Δt = 10 CN Matrices spectra (d) Δt = 10 Fourier symbol nd th Figure 1. Spectra of the M matrices for the 2 (red markers) and 4 order CN −3 (black markers) ﬁnite diﬀerences schemes: ℓ = 20 , N = 100 , Δt = 10 and −2 Δt = 10 . −5 still small values of δ , up to 10 , as illustrated in Figure 3. For larger values of δ , say −4 −3 −1 δ = 10 , δ = 10 , δ = 10 , eigenvalues of M are placed outside the unit disk making the relaxed scheme unstable, see Figures 4, 5 and 6. Below, in Table 1, we give the maximum and the minimum of the modulus of the eigenvalues of the M and M matrices for diﬀerent values of the time step Δt and r CN relaxation parameter δ . We can observe the inﬂuence of δ on the spectral radius of M : the relaxed scheme becomes unstable, i.e. ρ (M ) > 1 , for not suﬃciently small values of −4 δ , say, e.g. δ > 10 ; for small values of δ , we have ρ (M ) = 1 , the eigenvalues of M r r are perfectly matched on the unit circle and the scheme is then unconditionally stable. Relaxed formulations for dispersive wave equations 13 / 26 Eigenvalues Iter. Mat. Fourier Symbol of the CN Operator 1 1 Fourier symbol of CN CN with 2d order FD Unit circle 0.8 0.8 CN with 6th order FD 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −3 −3 (a) Δt = 10 CN Matrices spectra (b) Δt = 10 Fourier symbol Eigenvalues Iter. Mat. Fourier Symbol of the CN Operator 1 1 CN with 2d order FD Fourier symbol of CN CN with 6th order FD Unit circle 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −2 −2 (c) Δt = 10 CN Matrices spectra (d) Δt = 10 Fourier symbol nd th Figure 2. Spectra of the M matrices for the 2 (red markers) and 6 order CN −3 (black markers) ﬁnite diﬀerences schemes: ℓ = 100 , N = 400 , Δt = 10 −2 and Δt = 10 . 4.3. Simulation of a soliton To illustrate the eﬀect of the time relaxation, we consider the simulation of a soliton on the domain I = 0, ℓ , for diﬀerent values of δ and of Δt , namely u (x, t) = a sech κ (x − x − c t) , (4.1) with c = and κ = c . In our experiments presented below we use a = 0.8 , ℓ = 100 and x = . The interval I is discretized with N = 400 equidistant points and the space th discretization is realized by using 6 compact schemes. Of course, due to the nonlinear orbital stability properties of the solitons [35], the approximation to the solitary wave by J.-P. Chehab & D. Dutykh 14 / 26 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. 1 1 relaxed scheme relaxed scheme 0.8 0.8 CN scheme CN scheme 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −3 −2 (a) Δt = 10 (b) Δt = 10 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. 1 1 relaxed scheme 0.8 0.8 relaxed scheme CN scheme CN scheme 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 −1 −1 (c) Δt = 10 (d) Δt = 5 × 10 Figure 3. Spectra of M and M matrices. For the spatial discretisation we r CN −17 use ℓ = 100, N = 400 points and the relaxation parameter δ = 10 . Various values of the time step are taken. the numerical solution of the relaxed system cannot be considered in a long time interval, making the validation delicate. Indeed, the classical error measures such as L norms are irrelevant for orbits, meanwhile a numerical method approximates a classical solution and not an orbit. We observe that the relaxation allows to approach the exact solution with an expected − 17 level of accuracy. For instance, when δ = 10 , the solutions computed by the diﬀerent − 4 schemes coincide on a fairly long time interval. With δ = 10 and Δt = 0.1 , the solution coincide on a smaller time interval as it is expected. More generally, we notice in our numerical experiments that the pairs (δ , Δt) that make the relaxed scheme stable in the linear case are operant for the nonlinear KdV equation as well . The contrary is also However, we do not have a mathematical proof for this statement. Relaxed formulations for dispersive wave equations 15 / 26 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. relaxed scheme CN scheme relaxed scheme 0.8 CN scheme 0.6 0.5 0.4 0.2 0 0 -0.2 -0.4 -0.5 -0.6 -0.8 -1 -1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.5 0 0.5 1 −3 −2 (a) Δt = 10 (b) Δt = 10 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. relaxed scheme relaxed scheme 1 CN scheme 1 CN scheme 0.5 0.5 0 0 -0.5 -0.5 -1 -1 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 −1 −1 (c) Δt = 10 (d) Δt = 5 × 10 Figure 4. Spectra of M and M matrices. For the spatial discretisation we r CN −4 use ℓ = 100, N = 400 points and the relaxation parameter δ = 10 . Various values of the time step are taken. observed, a numerical instability holds, for example, for Δt = 0.005 and δ = 0.001 , see Table 1. The numerical result for the soliton propagation are shown in Figure 7 for two − 17 − 4 values of the relaxation parameter δ = 10 and δ = 10 . On the right panels we show the evolution of the L error computed thanks to the exact Solution (4.1). 4.4. Second order relaxed times schemes As illustrates above, relaxed schemes allow a correct numerical integration using only approximation to ∂ which is generally easily available in a number of situations, while the approximation to ∂ can be very tricky to build. However, the relaxation scheme is x x x only ﬁrst order accurate in time and as the semi-implicit Euler scheme, it reveals not to J.-P. Chehab & D. Dutykh 16 / 26 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. 1.5 relaxed scheme relaxed scheme CN scheme CN scheme 0.5 0.5 0 0 -0.5 -0.5 -1 -1 -1.5 -1 -0.5 0 0.5 1 -1.5 -1 -0.5 0 0.5 1 1.5 −3 −2 (a) Δt = 10 (b) Δt = 10 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. relaxed scheme relaxed scheme 1 CN scheme 1 CN scheme 0.5 0.5 0 0 -0.5 -0.5 -1 -1 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 −1 −1 (c) Δt = 10 (d) Δt = 5 × 10 Figure 5. Spectra of M and M matrices. For the spatial discretisation we r CN −3 use ℓ = 100, N = 400 points and the relaxation parameter δ = 10 . Various values of the time step are taken. be suited for long time interval simulations. To overcome this drawback we propose here to reach a second order of accuracy by using a classical Richardson extrapolation. In two words, the numerical time integration of ODE du = F (u) , dt by the forward Euler scheme deﬁnes the iterations k + 1 k k k u = u + Δt F (u ) = G (u ) , Δt Relaxed formulations for dispersive wave equations 17 / 26 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. 1.5 relaxed scheme relaxed scheme CN scheme CN scheme 0.5 0.5 0 0 -0.5 -0.5 -1 -1 -1.5 -1 -0.5 0 0.5 1 -1.5 -1 -0.5 0 0.5 1 1.5 −3 −2 (a) Δt = 10 (b) Δt = 10 Eigenvalues Iter. Mat. Eigenvalues Iter. Mat. 3 2 relaxed scheme relaxed scheme CN scheme CN scheme 1.5 0.5 0 0 -0.5 -1 -1 -2 -1.5 -3 -2 -3 -2 -1 0 1 2 3 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 −1 −1 (c) Δt = 10 (d) Δt = 5 × 10 Figure 6. Spectra of M and M matrices. For the spatial discretisation we r CN −1 use ℓ = 100, N = 400 points and the relaxation parameter δ = 10 . Various values of the time step are taken. which are ﬁrst order accurate approximations to u (k Δt) . The Richardson extrapolated sequence is deﬁned by v = G (u ) , 1 Δt v = G (u ) , 2 , 0 Δt/2 v = G (v ) , 2 , 1 Δt/2 2 , 0 k + 1 u = 2 v − v , 2 , 1 1 and is second order accurate in time. We will start here from a simple IMEX method, says backward Euler for the linear terms and forward for the nonlinear ones: in other words if F (u) writes as F (u) = − A u + H (u) , then the propagator is formally deﬁned as k −1 k k G (u ) = (I + Δt A) · u + Δt H (u ) . Δt J.-P. Chehab & D. Dutykh 18 / 26 Table 1. Minimal and maximal values of eigenvalues of iteration matrices for th N = 400 and ℓ = 100 . Here, 6 order Compact Schemes are used. N Δt δ M : σ , ρ M : σ , ρ CN CN CN r r r −3 −17 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1 CN CN r r −2 −17 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1 CN CN r r −1 −17 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1 CN CN r r −1 −17 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1 CN CN r r −3 −10 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.000000047 CN CN r r −2 −10 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.00000007 CN CN r r −1 −10 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.000000008 CN CN r r −1 −10 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.0000000015 CN CN r r −3 −8 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.000004775 CN CN r r −2 −8 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.00000691 CN CN r r −1 −8 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.000000798 CN CN r r −1 −8 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.0000001599 CN CN r r −3 −4 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.047217 CN CN r r −2 −4 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.070529 CN CN r r −1 −4 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.00796 CN CN r r −1 −4 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.001598 CN CN r r −3 −3 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.178956 CN CN r r −2 −3 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.56593 CN CN r r −1 −3 400 10 10 σ = 1, ρ = 1 σ = 1 , ρ = 1.07633 CN CN r r −1 −3 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.015844 CN CN r r −3 −1 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.02591 CN CN r r −2 −1 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.290614 CN CN r r −1 −1 400 10 10 σ = 1, ρ = 1 σ = 1, ρ = 2.90107786 CN CN r r −1 −1 400 5 × 10 10 σ = 1, ρ = 1 σ = 1, ρ = 1.909388 CN CN r r The IMEX relaxed scheme consists then in solving the following linear system at each step: (k + 1 (k) Δt (k) 2 I 0 −ΔtD U ) U − D (U ) N x x (k + 1 (k) · = + . − ΔtD (δ + Δt) I 0 V ) V 0 x N (k + 1 (k) 0 − ΔtD (δ + Δt) I W ) W 0 x N Relaxed formulations for dispersive wave equations 19 / 26 Comparison of the computed solutions at t=81.5 Error vs time on time interval [0,81.5] -2 0.9 Relaxed 0.8 Sanz-Serna Semi-explicit CN 0.7 Exact Sol -4 Err. Relaxed Err. Semi-explicit CN 0.6 Err. Sanz-Serna 0.5 -6 0.4 0.3 0.2 -8 0.1 -10 -0.1 0 20 40 60 80 100 0 10 20 30 40 50 60 70 80 90 time −3 −17 −3 −17 (a) Δt = 5 × 10 , δ = 10 , t = 139 (b) Δt = 5 × 10 , δ = 10 , L error Comparison of the computed solutions at t=80 Error vs time on time interval [0,80] -1 0.9 0.8 -2 0.7 Relaxed Sanz-Serna Err. Relaxed 0.6 -3 Err. Semi-explicit CN Semi-explicit CN 10 Err. Sanz-Serna Exact Sol 0.5 -4 0.4 0.3 -5 0.2 0.1 -6 -7 -0.1 0 20 40 60 80 100 0 10 20 30 40 50 60 70 80 90 time −1 −4 −1 −4 (c) Δt = 10 , δ = 10 , t = 120 (d) Δt = 10 , δ = 10 , L error Figure 7. Propagation of a solitary wave with relaxed schemes. (IMEX) The matrix of the system is noted by M . The Extrapolated Relaxed schemes writes Δt , δ as Algorithm 4 : Extrapolated Relaxed Scheme (0) 1: u given 2: 3: for k = 0, 1, ··· until convergence do (IMEX) Δt (k 4: Solve M v = − F (u ) , Δt/2, δ (IMEX) Δt 5: Solve M v = − F (u ) , 2 1 Δt/2, δ (IMEX) (k) 6: Solve M v = − Δt F (u ) , Δt, δ (k + 1) 7: Set u = 2 u − u . 2 3 8: end for Error in L norm Error in L norm J.-P. Chehab & D. Dutykh 20 / 26 Comparison of the computed solutions at t=140 Error vs time on time interval [0,140] -2 0.8 Err. Relaxed Relaxed 0.7 Err. Semi-explicit CN Sanz-Serna Err. Sanz-Serna Semi-explicit CN -3 0.6 Exact Sol 0.5 -4 0.4 0.3 -5 0.2 0.1 -6 -7 -0.1 0 20 40 60 80 100 0 50 100 150 time −1 −4 −1 −4 (a) Δt = 10 , δ = 10 , t = 140 (b) Δt = 10 , δ = 10 , t = 140 Comparison of the computed solutions at t=140 Error vs time on time interval [0,140] -2 0.8 Relaxed Err. Relaxed 0.7 Sanz-Serna Err. Semi-explicit CN -3 Semi-explicit CN Err. Sanz-Serna Exact Sol 0.6 -4 0.5 0.4 -5 0.3 -6 0.2 0.1 -7 -8 -0.1 0 20 40 60 80 100 0 50 100 150 time −2 −4 −2 −4 (c) Δt = 10 , δ = 10 , t = 140 (d) Δt = 10 , δ = 10 , t = 140 Figure 8. Comparison of the exact and numerical solutions to the KdV using Sanz–Serna and extrapolated Relaxed and IMEX schemes. Hereafter, in Figure 8 – 9, we present the comparison of the evolution of the soliton and its numerical approximations (Sanz–Serna and extrapolated Relaxed and IMEX schemes). We observe that the time extrapolation allows to approach the exact solution with a good accuracy on longer time intervals, when considering diﬀerent values of δ and −4 −1 Δt ; this is notable particularly, e.g. when taking δ = 10 and Δt = 10 . 5. Discussion Above we presented some rationale behind time-relaxed formulations for several well- known dispersive wave equations. The main conclusions and perspectives of our study are outlined below. Error in L norm Error in L norm Relaxed formulations for dispersive wave equations 21 / 26 Comparison of the computed solutions at t=140 Error vs time on time interval [0,140] -2 0.8 Relaxed Err. Relaxed 0.7 Sanz-Serna Err. Semi-explicit CN Semi-explicit CN Err. Sanz-Serna -3 Exact Sol 0.6 0.5 -4 0.4 0.3 -5 0.2 0.1 -6 -7 -0.1 0 20 40 60 80 100 0 50 100 150 time −1 −4 −1 −4 (a) Δt = 10 , δ = 10 , t = 140 (b) Δt = 10 , δ = 10 , t = 140 Comparison of the computed solutions at t=140 Error vs time on time interval [0,140] -2 0.8 Relaxed Err. Relaxed Sanz-Serna 0.7 Err. Semi-explicit CN Semi-explicit CN Err. Sanz-Serna -3 Exact Sol 0.6 0.5 -4 0.4 0.3 -5 0.2 0.1 -6 -7 -0.1 0 20 40 60 80 100 0 50 100 150 time −2 −4 −2 −4 (c) Δt = 10 , δ = 10 , t = 140 (d) Δt = 10 , δ = 10 , t = 140 Figure 9. Comparison of the exact and numerical solutions to the KdV using Sanz-Serna and extrapolated Relaxed and IMEX schemes. 5.1. Conclusions In this study we presented an approximate reformulation for several dispersive wave equations. This formulation was inspired somehow by quasi– (or pseudo–) compressibility methods to solve incompressible Navier–Stokes equations [21]. So, by following this ‘philosophy’ we proposed relaxed formulations for celebrated Korteweg–de Vries (KdV), Benjamin–Bona–Mahony (BBM) and Peregrine’s system of equations. However, it is obvious that the same technique can be extended to many other scalar and vectorial dispersive wave equations. This formulation has a simple advantage to involve ﬁrst order derivatives only while being in the form of a coupled evolution problem. This is the main diﬀerence with various local (or modiﬁed) reformulations used in continuous and discontinuous Galerkin methods [24, 33, 34, 36]. Many standard numerical methods can Error in L norm Error in L norm J.-P. Chehab & D. Dutykh 22 / 26 be applied to solve the proposed relaxed formulation numerically. In the present study we applied compact ﬁnite diﬀerences (with spectral-like resolution) to discretize the problem in space along with a simple time stepping. The presented numerical tests and validations for the KdV equation show that the relaxed discrete formulation possesses a larger CFL–type stability limit comparing to similar compact discretizations of the standard KdV equation (without relaxation). This preliminary conclusion indicates that relaxed schemes might be good candidates for the numerical simulation of notoriously stiﬀ systems of equations such as the KdV–KdV Boussinesq-type system extensively studied numerically in e.g. [5]. 5.2. Perspectives Above we presented some numerical illustrations for the classical KdV equation only. As it was mentioned in the previous Section, this formulation was extended to other weakly- nonlinear models as well. The main goal consists in extending these time-relaxation numer- ical methods to Boussinesq-type systems of weakly nonlinear weakly dispersive equations such as the celebrated classical Peregrine [28], modiﬁed Peregrine [12] or even fully nonlinear Serre–Green–Naghdi (SGN) equations [31] (some more conventional numer- ical strategies for these equations were outlined in e.g. [13, 25, 26]). On a more theoretical side, we would like to obtain the estimations of the diﬀerence between the perturbed and unperturbed problems solutions. In other words, we would like to have theoretical arguments to state that the proposed relaxation process generates solutions which remain close to those of the original equation. This is the main point on our ‘theoretical’ agenda. Acknowledgments D. Dutykh would like to acknowledge the hospitality of the Laboratory LAMFA and of the University of Picardie Jules Verne during his visit in November 2016. Reciprocally, J.-P. Chehab acknowledges the hospitality of the Laboratory of Mathematics (LAMA UMR #5127) and of the University Savoie Mont Blanc during his visit in December A. Other types of dispersive equations In this Appendix we provide two other examples of relaxation for some widely used dispersive wave equations (scalar and system case). Relaxed formulations for dispersive wave equations 23 / 26 A.1. Benjamin–Bona–Mahony equation The celebrated Benjamin–Bona–Mahony (BBM) equation was derived ﬁrst in [3, 27] and can be recast in the following dimensionless form: u + u u − u = 0 . t x xxt Similarly to the KdV case the variable u (x, t) can be the free surface elevation or a horizontal ﬂuid velocity either. It is possible to propose a similar relaxed formulation for the BBM equation as well. Indeed, let us recast it ﬁrst in a conservative form: 1 2 u − u + u = 0 . xx t x nd The 2 derivative can be lowered by introducing additional variables: 1 2 u − w + u = 0 , t 2 x u = v , v = w . def It is convenient to introduce a new variable β (x, t) := u − w (x, t) to simplify the ﬁrst equation: 1 2 β + ( β + w) = 0 , t 2 (β + w) = v , v = w . The last step consists in adding relaxation terms to obtain an evolutionary system in all variables: β + ( β + w) = 0 , t 2 δ v + v − ( β + w) = 0 , t x δ w + w − v = 0 , t x with δ ≪ 1 being again a small parameter. A.2. Peregrine system The Peregrine system was proposed by D. H. Peregrine (1967) in [28]. In a par- ticular case of a ﬂuid layer of constant depth (i.e. the even bottom case) this system reads η + (d + η ) u = 0 , (A.1) 1 2 u + u u + g η − d u = 0 , (A.2) t x xxt x 3 where d > 0 is the constant ﬂuid depth and g > 0 is the gravity acceleration. The vari- ables η (x, t) and u (x, t) are the free surface elevation and the depth-averaged horizontal velocity correspondingly. The sketch of the ﬂuid domain is shown in Figure 10. J.-P. Chehab & D. Dutykh 24 / 26 Figure 10. Sketch of the ﬂuid domain with a free surface. As the ﬁrst step, the Peregrine system (A.1), (A.2) is rewritten in the following conservative form: η + (d + η ) u = 0 , 1 2 1 2 u − d u + u + g η = 0 . xx 3 2 t x The last conservative Peregrine system can be lowered and relaxed similar to the BBM case (see the Section A.1) by introducing new variables: 1 2 η + (d + η ) (ρ + d w) = 0 , 3 x 1 1 2 2 ρ + (ρ + d w) + g η = 0 , 2 3 1 2 δ v + v − ρ + d w = 0 , δ w + w − v = 0 , t x def 1 2 where δ ≪ 1 and ρ (x, t) := u (x, t) − d w (x, t) is an auxiliary variable. If one takes the limit δ → 0 , we shall recover a formulation similar to what is used in local (or modiﬁed) continuous Galerkin methods [33, 34]. Remark 2. The celebrated KdV and BBM models can be written for the variable u (x, t) being the free surface elevation η (x, t) or the horizontal velocity (after a few changes of variables). In this way these equations appear in the theory of water waves. However, the KdV and BBM-type equations appear in many other physical settings as well (see e.g. [11, 17, 18, 20, 30]). References [1] D. C. Antonopoulos, V. A. Dougalis, and D. E. Mitsotakis. Galerkin approximations of the periodic solutions of Boussinesq systems. Bulletin of Greek Math. Soc., 57:13–30, 2010. 4 Relaxed formulations for dispersive wave equations 25 / 26 [2] M. Antuono, V. Y. Liapidevskii, and M. Brocchini. Dispersive Nonlinear Shallow-Water Equations. Stud. Appl. Math., 122(1):1–28, 2009. 4 [3] T. B. Benjamin, J. L. Bona, and J. J. Mahony. Model equations for long waves in nonlinear dispersive systems. Philos. Trans. Royal Soc. London Ser. A, 272:47–78, 1972. 23 [4] F. Benkhaldoun and M. Seaïd. New ﬁnite-volume relaxation methods for the third-order diﬀerential equations. Commun. Comput. Phys., 4:820–837, 2008. 4 [5] J. L. Bona, V. A. Dougalis, and D. E. Mitsotakis. Numerical solution of KdV-KdV systems of Boussinesq equations: I. The numerical scheme and generalized solitary waves. Mat. Comp. Simul., 74:214–228, 2007. 22 [6] P. Bonneton, F. Chazel, D. Lannes, F. Marche, and M. Tissier. A splitting approach for the fully nonlinear and weakly dispersive Green-Naghdi model. J. Comput. Phys., 230:1479–1498, 2011. 4 [7] J. V. Boussinesq. Essai sur la théorie des eaux courantes. Mémoires présentés par divers savants à l’Acad. des Sci. Inst. Nat. France, XXIII:1–680, 1877. 5 [8] H. Chen, M. Chen, and N. Nguyen. Cnoidal Wave Solutions to Boussinesq Systems. Non- linearity, 20:1443–1461, 2007. 4 [9] R. Cienfuegos, E. Barthélémy, and P. Bonneton. A fourth-order compact ﬁnite volume scheme for fully nonlinear and weakly dispersive Boussinesq-type equations. Part I: Model development and analysis. Int. J. Numer. Meth. Fluids, 51:1217–1253, 2006. 4 [10] D. Clamond. Cnoidal-type surface waves in deep water. J. Fluid Mech, 489:101–120, jul 2003. 4 [11] A. Duran, D. Dutykh, and D. Mitsotakis. On the Galilean Invariance of Some Nonlinear Dispersive Wave Equations. Stud. Appl. Math., 131(4):359–388, nov 2013. 24 [12] A. Durán, D. Dutykh, and D. Mitsotakis. Peregrine’s System Revisited. In N. Abcha, E. N. Pelinovsky, and I. Mutabazi, editors, Nonlinear Waves and Pattern Dynamics, pages 3–43. Springer International Publishing, Cham, 2018. 22 [13] D. Dutykh, D. Clamond, P. Milewski, and D. Mitsotakis. Finite volume and pseudo-spectral schemes for the fully nonlinear 1D Serre equations. Eur. J. Appl. Math., 24(05):761–787, 2013. 4, 22 [14] D. Dutykh, T. Katsaounis, and D. Mitsotakis. Finite volume schemes for dispersive wave propagation and runup. J. Comput. Phys., 230(8):3035–3061, apr 2011. 4 [15] D. Dutykh, T. Katsaounis, and D. Mitsotakis. Finite volume methods for unidirectional dispersive wave models. Int. J. Num. Meth. Fluids, 71:717–736, 2013. 4 [16] C. Eskilsson and S. J. Sherwin. Discontinuous Galerkin Spectral/hp Element Modelling of Dispersive Shallow Water Systems. J. Sci. Comput., 22:269–288, 2005. 4 [17] F. Fedele and D. Dutykh. Vortexons in axisymmetric Poiseuille pipe ﬂows. EPL, 101(3):34003, feb 2013. 24 [18] R. Grimshaw. Internal Solitary Waves. In R. Grimshaw, editor, Environmental Stratiﬁed Flows, pages 1–27. Springer US, 2002. 24 [19] M. S. Ismail. A ﬁnite diﬀerence method for Korteweg-de Vries like equation with nonlinear dispersion. International Journal of Computer Mathematics, 74(2):185–193, 2000. 4 [20] R. S. Johnson. A modern introduction to the mathematical theory of water waves. Cambridge University Press, Cambridge, 1997. 4, 5, 24 [21] M. Kameyama, A. Kageyama, and T. Sato. Multigrid iterative algorithm using pseudo- compressibility for three-dimensional mantle convection with strongly variable viscosity. J. J.-P. Chehab & D. Dutykh 26 / 26 Comput. Phys, 206:162–181, 2005. 4, 21 [22] D. J. Korteweg and G. de Vries. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Phil. Mag., 39(5):422–443, 1895. 5 [23] S. K. Lele. Compact ﬁnite diﬀerence schemes with spectral-like resolution. J. Comp. Phys., 103(1):16–42, nov 1992. 8 [24] D. Levy, C.-W. Shu, and J. Yan. Local discontinuous Galerkin methods for nonlinear disper- sive equations. J. Comput. Phys., 196(2):751–772, 2004. 5, 6, 21 [25] D. Mitsotakis, D. Dutykh, and J. Carter. On the nonlinear dynamics of the traveling-wave solutions of the Serre system. Wave Motion, 70:166–182, apr 2017. 22 [26] D. Mitsotakis, B. Ilan, and D. Dutykh. On the Galerkin/Finite-Element Method for the Serre Equations. J. Sci. Comput., 61(1):166–195, feb 2014. 4, 22 [27] D. H. Peregrine. Calculations of the development of an undular bore. J. Fluid Mech., 25(02):321–330, mar 1966. 23 [28] D. H. Peregrine. Long waves on a beach. J. Fluid Mech., 27:815–827, 1967. 22, 23 [29] A. V. Porubov and G. A. Maugin. Propagation of localized longitudinal strain waves in a plate in the presence of cubic nonlinearity. Phys. Rev. E, 74:46617, 2006. 4 [30] H. Schamel. A modiﬁed Korteweg-de Vries equation for ion acoustic wavess due to resonant electrons. J. Plasma Phys., 9(03):377–387, mar 1973. 24 [31] F. Serre. Contribution à l’étude des écoulements permanents et variables dans les canaux. La Houille blanche, 8:374–388, 1953. 22 [32] J. J. Stoker. Water Waves: The mathematical theory with applications. Interscience, New York, 1957. 5 [33] M. Walkley and M. Berzins. A ﬁnite element method for the one-dimensional extended Boussinesq equations. Int. J. Num. Meth. Fluids, 29(2):143–157, jan 1999. 5, 21, 24 [34] M. Walkley and M. Berzins. A ﬁnite element method for the two-dimensional extended Boussinesq equations. Int. J. Num. Meth. Fluids, 39(10):865–885, aug 2002. 5, 21, 24 [35] M. I. Weinstein. Existence and dynamic stability of solitary wave solutions of equations arising in long wave propagation. Comm. Partial Diﬀ. Eqns., 12(10):1133–1173, jan 1987. [36] J. Yan and C.-W. Shu. A local discontinuous Galerkin method for KdV type equations. SIAM J. Num. Anal., 40:769–791, 2002. 4, 5, 6, 21 J.-P. Chehab: Université de Picardie Jules Verne, LAMFA CNRS UMR 7352, 33, rue Saint-Leu, 80039 Amiens, France E-mail address: Jean-Paul.Chehab@u-picardie.fr URL: http://www.lamfa.u-picardie.fr/chehab/ D. Dutykh: Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LAMA, 73000 Chambéry, France and LAMA, UMR 5127 CNRS, Université Savoie Mont Blanc, Campus Scientifique, 73376 Le Bourget-du-Lac Cedex, France E-mail address: Denys.Dutykh@univ-savoie.fr URL: http://www.denys-dutykh.com/
Nonlinear Sciences – arXiv (Cornell University)
Published: Mar 6, 2019
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.