STUDIES IN RECURRENCE IN NONLINEAR DISPERSIVE WAVE EQUATIONS
|
|
- Bertram Austin
- 5 years ago
- Views:
Transcription
1 STUDIES IN RECURRENCE IN NONLINEAR DISPERSIVE WAVE EQUATIONS T. Arbogast, J. L. Bona, and J.-M. Yuan December 25, 29 Abstract This paper is concerned with the study of recurrence phenomena in nonlinear dispersive wave equations using numerical simulation. Two classes of equations are investigated, namely, the generalized KdV-equations, u t + u p u x + ɛu xxx =, and the generalized RLW- or BBM-equations, u t + u x + u p u x ɛu xxt =. The numerical evidence presented here shows that even non-integrable equations exhibit recurrence phenomena. We also determine an approximate relationship between recurrence time and amplitude of the initial data. 1 Introduction In the last few decades, a new paradigm has arisen in mathematics. In this conception, numerical simulation and experimentation have been added to the arsenal of techniques available to aid our understanding of mathematical issues. Such a paradigm is also common in physics, and more recently the other sciences and even some of the social sciences. It is distinct from the use of computation in engineering and also in science and social science, where the goal is to make direct, numerically-valued predictions which are used in forecasting, design, control and the like. Early examples of the use of computers to uncover mathematical information include the famous Birch and Swinnerton-Dyer conjectures in number theory and the pioneering work on what came to be known as solitons by Fermi, Pasta, and Ulam. Despite the relative lack of computational power in the 195 s, these works and others like them have proved to be of lasting importance. In the present study, we continue in the line pioneered by Fermi, Pasta and Ulam [16] and later by Zabusky and Kruskal [38]. In [16], numerical simulations indicated that a mass and spring model with non-hooksian springs (modelled as an an-harmonic lattice) features near recurrence of an initial perturbation. This model, which had small cubic and quartic additions to the usual harmonic potential governing the springs, was initially conceived as a possible model for heat conduction. With a sinusoidal initial disturbance, one expected a kind of thermalization, or equipartition of energy into all the modes of oscillation. Thus 1
2 the near recurrence of the initial state was unexpected. Upon reflection, this outcome might not seem surprising in the light of Poincaré s famous observations of chaotic patterns in the dynamics of planetary interactions and his associated recurrence theorem. However, the recurrence observed by Fermi, Pasta and Ulam was different from Poincaré recurrence. For one thing, the recurrence they observed happens on a much shorter time scale. Zabusky and Kruskal [38] revisited the Fermi-Pasta-Ulam system (henceforth the FPU system) in They took a continuum limit of the FPU lattice and came to the Korteweg-de Vries equation (KdV-equation) u t + uu x + ɛu xxx =. (1) Following Fermi et al., they posed this equation as an initial-value problem with periodic boundary conditions and simulated the problem numerically with ɛ =.22 and u(x, ) = cos(πx), x 2. They also found recurrence of initial states for (1). Later, Zabusky [37] pursued in more detail the recurrence phenomenon for the KdV-equation. Fixing the initial data as u(x, ) = cos(2πx), he denoted by t B, the time until shock formulation for (1) with ɛ =, and by t R = t R ( ɛ) the recurrence time for (1) with the given value of ɛ. On the basis of numerical studies of Tappert [33] using a split-step Fourier method, Zabusky suggested the relationships and t R t B =.7 ɛ (2) the number of emergent solitons =.2 ɛ. (3) Indeed, this was another finding of Zabusky and Kruskal s numerical study [38], that general classes of initial data resolved into the special travelling waves called solitary waves by Scott Russell [3], first written analytically by Boussinesq [13], and later termed solitons by Zabusky and Kruskal. Theoretical support for these relationships was later put forward by Lax [21], Novikov [27], Toda [34], McKean & Trubowitz [22], McKean [23], and Bourgain [12]. The recurrence phenomenon for the periodic initial-value problem for the KdV-equation (1) has been linked to the existence of quasi-periodic motion of completely integrable Hamiltonian systems. It is natural to inquire as to whether or not other infinite-dimensional Hamiltonian systems might have the same property. This is the principal perspective of the present study, which is made in the context of the generalized KdVequation (GKdV equation henceforth), u t + u p u x + u xxx =, (4) 2
3 and its counterpart, the generalized BBM-equation (GBBM-equation) u t + u x + u p u x u xxt =. (5) These equations are all Hamiltonian. However, except for (4) with p = 1 or 2, they do not appear to be integrable. In particular, their solitary-wave solutions do not interact exactly as they do for (4) with p = 1 or 2 (see e.g. [8]). However, the GKdV equation for p = 3 and the GBBM-equations for all p seem to have the property of resolution into solitary waves. More precisely, and following Lax [2], for suitably restricted initial values u(, ) = ψ, the solution u emanating therefrom has the property that u(x, t) = k φ ci (x c i t + θ i (t)) + r(x, t) i=1 where φ ci is a solitary-wave solution of speed c i, θ i (t) θ i ( ) as t for 1 i k, and lim t sup x R r(x, t) =. (Thus, while the L 2 -norm of r may remain positive as t, its L -norm tends algebraically to zero, as is characteristic of dispersive phenomena.) It is thought that the property of near recurrence is closely related to the property of resolution into solitary waves. In a different, completely integrable context, Camassa and Lee [14] have numerical simulations indicative of this point. It will follow from this that in cases where the remainder r has significant energy (i.e. its L 2 -norm is a significant fraction of the L 2 -norm of the initial data), the recurrence is less noticeable. The paper is organized as follows. In Section 2, the numerical methods are introduced and convergence and accuracy ascertained. Section 3 is devoted to the recurrence for the generalized KdV equations. In Section 4 the recurrence for the generalized BBM-equation is presented. The overall conclusions are summarized in Section 5. 2 A Numerical Method It is our purpose to investigate numerically recurrence phenomenon for the nonlinear dispersive wave equations (4) and (5), with p a positive integer. The numerical schemes used are based on a q-stage Gauss-Legendre method for the temporal discretization and a finite-element method for the spatial discretization (see [7] and its precursor [6]). The following convergence theorem for such a method applied to (4) was proved by Bona et al. in [7]. Theorem 2.1. Suppose that, as the spatial discretization parameter h +, the temporal discretization k obeys the hypotheses 1. if p = 1, k = O(h 3/(2(q+2)) ) for q 2 or k = O(h 3/4 ) for q = 1, 2. if p = 2, kh 1/2 is sufficiently small for q 2 or k = O(h 3/4 ) for q = 1, 3
4 3. if p 3, kh 1 is sufficiently small for all q 1. Then for h sufficiently small, there exists a unique solution U n of the fully discrete numerical scheme, which is a q-stage Gauss-Legendre time stepping coupled with smooth splines of degree r 1 finite element approximation in space, such that max U n u n c(k 2 + h r ) for q = 1, (6) n J max U n u n c(k q+2 + h r ) for q 2, (7) n J where J is fixed, r 2, U n represents the numerical solution at time t = nk and u n represents the analytical solution at the same moment. A theorem of this nature is not only of academic interest. It tells us what to expect from a properly written and coded numerical scheme. If we do not see these rates of convergence, we can be pretty sure something is awry. In any event, as we are asking the reader to believe the results of our simulations, it behooves us to provide evidence. In the earlier work [7], such tests were performed on solitary-wave solutions. That was appropriate, as perturbations of solitary waves figured prominently. Here the simulations feature interaction of solitary-wave structures and so we deemed it more appropriate to test the code on multi-soliton solutions. Consider first the KdV equation rescaled in the form and pose it with the initial datum u t + 1uu x +.25u xxx = (8) u(x, ) = cosh(4 1(x.49)) + cosh(8 1(x.13)) (3 cosh(2 1(x +.23)) + cosh(6 1(x.25))) 2. Using Inverse Scattering theory [17, 18], the exact solution on R of (8) with the above initial datum is seen to be u(x, t) = cosh(4 1(x 4t.49)) + cosh(8 1(x 16t.13)) (3 cosh(2 1(x 28t +.23)) + cosh(6 1(x 12t.25))). 2 (9) As t, and along the line x 16t = constant, the asymptotic behavior of u(x, t) is 4.8 sech 2 (4 1(x 16t x 2 ). If instead, x 4t is fixed, the asymptotic behavior of u(x, t) will look like 1.2 sech 2 (2 1(x 4t x 1 )). Since the error term is exponentially small, it is straightforward to deduce that as t, u(x, t) 4.8 sech 2 (4 1(x 16t x 2 )) sech 2 (2 1(x 4t x 1 )) (1) The amplitudes of the two emerging solitons are thus 4.8 and 1.2, respectively. As is well-known, phase shifts appear during the interaction. The taller soliton has moved forward and the shorter one backward, compared to the positions 4
5 they would have reached if there had been no interaction. The phase shifts for the solution in (9) from the formula in [35] are forward shift of the larger soliton after interaction;.88685, backward shift of the smaller soliton after interaction; As a promising preliminary indication of accuracy, the phase shifts after the interaction in our numerical tests turn out to be forward shift of the larger soliton after interaction;.9, backward shift of the smaller soliton after interaction;.17. We compare the exact 2-soliton solution (9) of the KdV-equation on the real line with the numerical solution on the periodic domain [,1] for accuracy (see [4]). The outcome of this study is reported in Table 1 with h =.1 and k =.25. It shows good accuracy as expected from Theorem 2.1. Note that the taller soliton catches up and interacts with the shorter one around t =.31. Normalized L 2 error. 2.6E E E E E E E-6 Table 1: Error between the exact 2-soliton solution and the numerical solution 3 Recurrence for the generalized KdV equations Equation (4) with p = 2 is usually called the modified-kdv-equation, and like the KdV-equation (3), it is in a particular sense completely integrable. For p 3, as far as one can tell, there is no analog of the complete integrability that holds when p = 1 or 2 (see e.g.[9]). Results of McLeod and Olver [28] indicate that (5) with p = 1 is not integrable, and numerical simulations show convincingly that the solitary-wave solutions of (5) for various values of p are not solitons, which is to say they do not interact exactly. Indeed, they are not all stable for p > 4 (see [8], [11], and [25]). On the other hand, the KAM-theory of Kolmogoroff, Arnold, Moser and Siegel (see [2]) states that a weakly perturbed, completely integrable system retains the property of possessing quasi-periodic solutions. Hence at least in certain parameter regimes, the recurrence of initial states might be expected to persist for suitable Hamiltonian perturbations of (4) or (5) with p = 2, say. A Hamiltonian for the generalized KdV-equation (4) is 5
6 H(u, u x ) = 1 [ ] ɛu2 x (x, t) p + 2 up+2 (x, t) dx. (11) It is far from clear that the case p = 3, say, corresponds to a small perturbation of an integrable case p = 1 or 2! We begin with some numerical simulations involving the exact multi-soliton solutions obtained from Inverse Scattering theory. The relevant formulas may be found, for example, in Whitham s text [35]. The exact solutions are ideal for testing the implementation of the numerical scheme in a context that is quite relevant to recurrence phenomena. The following natural quantity will be used throughout in our discussion of recurrence. It is the L 2 (, 1)-norm of the difference between the solution at time t and the initial value, modulo translation and periodicity. Definition 3.1. Let u(x, t) be a solution to one of the evolution equations studied here, emanating from nonzero initial data u(x, ). The shape difference is f(t) u(x,) 2, where f(t) = min θ(t) 1 1 ( u(x + θ(t), ) u(x, t)) 2dx. (12) We take θ() = so that f() =. A recurrence time T r > is a time at which the shape difference is again small. Of course, T r will be at least a local minimum of f. While it is not obvious how to make this notion more precise, it will become clear that in practice these recurrence times do exist. The recurrence time t r > closest to is the first recurrence time. It will transpire that later recurrence times are often more or less integer multiples of t r. 3.1 The KdV equation (p = 1) In the following experiments, we report recurrence results for these test cases: 1. two-solition, 2. four-soliton with dispersive tail. Test 1. (two-soliton) By rescaling, (1) is written in the form u t + 1uu x u xxx = (13) and posed with initial datum u(x, ) = 3.6sech 2 (1 1(x.5)). 6
7 As in (1), the exact solution of (13) on R is u(x, t) = cosh(2 1(x 4t.5)) + cosh(4 1(x 16t.5)) (3 cosh(1 1(x 28t.5)) + cosh(3 1(x 12t.5))) 2. As in (11), as t, (14) u(x, t) 4.8sech 2 (2 1(x 16t x 2 )) + 1.2sech 2 (1 1(x 4t x 1 )). (15) The recurrence results are shown in Table 2 when the problem is posed as a periodic initial-value problem with the initial datum effectively having compact support in the period domain. T r shape difference Table 2: Two-soliton initial data. T r is the computed value of the recurrence time with h =.2 and k = The positions of 2 solitons Figure 1: Two-soliton initial data. The dark line in the figure at the right represents the peak of the smaller soliton. The lighter line tracks the peak of the taller soliton. Remark 3.1. From Table 2 and the graph of the shape difference (see Figure 1), one sees a definite value.79 of the recurrence time and a small linear growth in the shape difference at multiples of the recurrence time. The number of solitons is now increased to ascertain how well the program functions in a more complex environment. Test 2. (four-soliton plus dispersive wave) We solve u u + 1uu x u xxx = (16) 7
8 with u(x, t) = 9.sech 2 (1 1(x.5)). Using the Inverse Scattering theory, the solution of (16) is seen to have the form u(x, t) sech 2 ( (x 36t x 1 )) sech 2 ( (x 16t x 2 )) sech 2 ( (x 4t x 3 )) sech 2 ( (x 4t x 4 )) + dispersive tail (17) as t. The recurrence result for (16) is in Table 3 and Figure 2. T r shape difference Table 3: Test 2. T r is the computed value of the recurrence time with h =.1 and k = Figure 2: Test 2. Recurrence results for initial data with four-soliton plus a dispersive tail in the KdV-equation 3.2 The modified KdV equation (p = 2) In this subsection, recurrence is studied in the context of the modified KdV equation. Directing attention back to the generalized KdV-equation, the case p = 2 is known to be completely integrable. In each numerical experiment, we report numerically generated approximations just as in the previous section. We approximated solutions of the mkdv equation u t + u 2 u x u xxx = (18) 8
9 with the initial condition u(x, ) = a cos(2π(x.5)) 2, x 1, (19) taking the six cases: a=.6,.75,.875, 1., 1.1, and The recurrence times are given in Tables 4-6. Snapshots of the solutions at recurrence times are available in [36]. a=.6 a= Table 4: a=.6 and.75. T r is the computed value of the recurrence time with h=.2 and k= shape difference Figure 3: The left figure represents the recurrence result for u(x, ) =.6 cos 2 (2π(x.5)). The right figure represents the recurrence result for u(x, ) =.75 cos 2 (2π(x.5)). a=.875 a= Table 5: a=.875 and 1.. T r is the computed value of the recurrence time with h=.2 and k=.1. 9
10 Figure 4: The left figure represents the recurrence result for u(x, ) =.875 cos 2 (2π(x.5)). The right figure represents the recurrence result for u(x, ) = 1. cos 2 (2π(x.5)). a=1.1 a= Table 6: a=1.1 and T r is the computed value of the recurrence time with h=.2 and k= Figure 5: The left figure represents the recurrence result for u(x, ) = 1.1 cos 2 (2π(x.5)). The right figure represents the recurrence result for u(x, ) = 1.25 cos 2 (2π(x.5)). 1
11 amplitude of initial data t r Table 7: The amplitudes of the initial data a cos 2 (2π(x.5)) and their corresponding first recurrence times. A search is made for a relationship between the first recurrence time and the amplitude a. Such a relation is obtained for the fixed value.121 of ɛ and variations of a as indicated. The data is summarized in Table 7. Using least squares fitting of t r = A a B + C, the following equation emerges: t r = (amplitude) (2) It appears that B = 4 3 and C = 22. This result is shown in Figure y = *x**( ) y = *x**( 4/3) The fisrt recurrence time The fisrt recurrence time Amplitude Amplitude Figure 6: The relationship between the first recurrence time of the mkdv equation and various amplitudes of cos 2 (2π(x.5)), according to the data in Table 7 with the solid line being the least squares fit. The left plot shows the full fit to t r = A a B + C; the right plot fits only t r = A a
12 3.3 The generalized KdV equation with p = 3 In this subsection, recurrence is studied in the context of the generalized KdV equation with p = 3. The case p = 2 is known to be integrable, as mentioned before. However, the evidence suggests that for p 3, the evolution equation (4) is not integrable. We use the two-stage diagonally implicit Runge-Kutta (DIRK) method in time and the finite-element method with 1-periodic smooth splines of order 4 (degree 3) in space. We approximated solutions of the equation u t + u 3 u x u xxx = (21) numerically via the scheme just mentioned with the initial conditions u(x, ) = a cos 2 (2π(x.5)), x 1. (22) The amplitude a was given the values: a=.8,.825,.85,.875,.9,.925, and.95. The results are shown in Tables 8-11 and Figures 7-1. a=.8 a= Table 8: a=.8 and.825. T r is the computed value of the recurrence time with h =.2 and k = Figure 7: The left figure is the recurrence result for a=.8. The right figure is the recurrence result for a=
13 a=.85 a= Table 9: a=.85 and.875. T r is the computed value of the recurrence time with h =.2 and k = Figure 8: The left figure is the recurrence result for a=.85. The right figure is the recurrence result for a=.875. a=.9 a= Table 1: a=.9 and.925. T r is the computed value of the recurrence time with h =.2 and k = Figure 9: The left figure is the recurrence result for a=.9. The right figure is the recurrence result for a=
14 T r shape difference Table 11: a=.95. T r is the computed value of the recurrence time with h =.1 and k = amplitude t=2.331 Figure 1: The left figure is the recurrence result for a=.925. The right figure is the solution at the first recurrence time. amplitude of initial data t r Table 12: The amplitudes of the cos 2 initial data and the corresponding first recurrence time of the solution. 14
15 There is clear evidence of recurrence for (22)-(23). Using least squares fitting of t r = A a B +C on the data summarized in Table 12, there appears the equation for which B 11 1 t r = (amplitude) , (23) and C 3 (see Fig. 11). 3.5 y = *x**( ) y = 5.5*x**( 11/1) 3 The fisrt recurrence time The fisrt recurrence time Amplitude Amplitude Figure 11: The relationship between the first recurrence of GKdV and various amplitudes of cos 2 (2π(x.5)) according to the data in Table 12 with least squares fitting of t r = A a B + C on the left and t r = A a on the right. Remark 3.2. From the information acquired in Sections 3.2 and 3.3 (Figs. 6 and 11), one is tempted to hazard the guess t r c p a p 2 +2p 4 p 2 +2p 5 + dp, where c p and d p are constants related to p. 4 Recurrence for the generalized BBM-equation In this section, recurrence is studied in the context of the BBM-equation. In 1972, Benjamin et al. [3] found that, under the same approximations and assumptions that originally led Boussinesq and Korteweg and de Vries to their equation, the partial differential equation u t + u x + uu x ɛu xxt = (24) could equally well be justified as a model. The above equation is called the BBMequation or regularized-long-wave (RLW) equation. While the KdV equation 15
16 is an integrable equation and has infinitely many conservation laws [26], Olver [28] proved that the BBM-equation has only three conservation laws. Bona, Pritchard, and Scott [1] presented a quantitative comparison between the solutions to the initial-value problem for KdV and BBM. In [9], the same authors studied the interaction between two solitary waves of the BBM-equation. They showed that there were dispersive tails generated after the interaction. Thus it is provisionally inferred that the BBM-equation is not integrable. BBM-type equations are not temporally stiff as are explicit discretizations of the KdV-type equations. J. Pasciak [29], B.-Y. Guo and V. Manoranjan [19], and D. Sloan [31] have developed spectral methods for BBM-equation. Pasciak [29] used a fourth order Runge-Kutta method for the temporal discretization. Here we use a finite element method in space and a 2-stage Gauss-Legendre method in time. The three conservation laws for the GBBM-equation u t + u x + u p u x ɛu xxt = (25) are (see [32]) 1 u(x, t)dx = c 1, (26) 1 1 u 2 (x, t)dx + ɛ u 2 (x, t) + 1 u 2 x (x, t)dx = c 2, (27) 2 (p + 1)(p + 2) up+2 (x, t)dx = c 3. (28) As the second conservation law (27) is equivalent to the H 1 -norm, it seems natural to use that functional to measure recurrence in the equations (25). Definition 4.1. Let u(x, t) be a solution emanating from initial data u(x, ). g(t) The shape difference for the BBM-type equation is u(x,), where H 1 and g(t) = +ɛ { 1 min θ(t) 1 1 u(x, ) H 1 = ( u(x + θ(t), ) u(x, t) ) 2dx ( ux (x + θ(t), ) u x (x, t) ) 2 dx } The BBM-equation (p=1) First, solutions of the BBM-equation (29) (u 2 (x, ) + ɛu 2 x (x, ) ) dx. (3) u t + u x + uu x 1 4 u xxt = (31) 16
17 were approximated numerically with the initial condition u(x, ) = a cos 2 (2π(x.5)), x 1. (32) Experiments were run with a =.75,.875, 1, 1.25, 1.5, and 2. Results appear in Tables and Figures a=.75 a= Table 13: a =.75 and.875. T r is the computed value of the recurrence time with h =.4 and k = Figure 12: The left figure is the recurrence result for a=.75. The right figure is the recurrence result for a=.875. a=1. a= Table 14: a=1. and T r is the computed value of the recurrence time with h =.2 and k =.1. 17
18 Figure 13: The left figure is the recurrence result for a=1.. The right figure is the recurrence result for a=1.25. a=1.5 a= Table 15: a=1.5 and 2.. T r is the computed value of the recurrence time with h =.2 and k = Figure 14: The left figure is the recurrence result for a=1.5. The right figure is the recurrence result for a=2.. 18
19 amplitude of the initial data t r Table 16: The amplitudes of the initial data and their corresponding first recurrence times. A relationship may be deduced between the first recurrence time and the amplitude of the cos 2 initial data for the BBM-equation with fixed ɛ =.1. From the data, summarized in Table 16, least squares fitting gives the equation (see Figure 15) of t r = A a B + C, t r = (amplitude) (33) 2 y = *x**( ) The fisrt recurrence time Amplitude Figure 15: The relationship between the first recurrence of BBM and the amplitude of cos 2 (2π(x.5)). 4.2 The modified BBM-equation (p=2) In this subsection, recurrence is studied in the context of the modified BBMequation. Solutions of the mbbm-equation u t + u x + u 2 u x ɛu xxt =, (34) were approximated numerically, where ɛ = 1 4, with the initial condition u(x, ) = a cos 2 (2π(x.5)), x 1. (35) 19
20 Six simulations were run with a =.75,.85,.9, 1, 1.15, and Results are shown in Tables and Figures a=.75 a= Table 17: a=.75 and.85. T r is the computed value of the recurrence time with h =.4 and k = Figure 16: The left figure is the recurrence result for a=.75. The right figure is the recurrence result for a=.85. a=.9 a= Table 18: a=.9. T r is the computed value of the recurrence time with h =.2 and k =.1. As before, a relationship is inferred between the first recurrence time and the amplitude of the cos 2 initial data for the modified BBM-equation with fixed ɛ =.1. From the data in Table 2, least squares fitting yields which is plotted in Figure 19. t r = (amplitude) , (36) 2
21 Figure 17: The left figure is the recurrence result for a=.9. The right figure is the recurrence result for a=1.. a=1.15 a= Table 19: a=1.15 and T r is the computed value of the recurrence time with h =.2 and k = Figure 18: The left figure is the recurrence result for a=1.15. The right figure is the recurrence result for a=1.25. amplitude of the initial data t r Table 2: The amplitudes of the initial data and their corresponding first recurrence times. 21
22 8.5 y = *x**( ) The fisrt recurrence time Amplitude Figure 19: The relationship between the first recurrence of mbbm and the amplitude of cos 2 (2π(x.5)). The solid line is the least-squares approximation of the data. 4.3 The generalized BBM-equation with p = 3 In this subsection, the generalized BBM-equation u t + u x + u 3 u x ɛu xxt = (37) is considered with u(x, ) = cos 2 (2π(x.5)), different values of ɛ and periodic boundary conditions on [,1]. In the cubic case where the nonlinearity is u 3 u x, this term is much smaller that the other terms when u < 1, and so the equation is essentially linear. Recurrence on relatively short time scales requires nonlinearity, however. In this section we have elected to vary the parameter ɛ instead of the amplitude of the initial data. Simple rescaling shows that varying ɛ is equivalent to varying the period of the initial data. The 9 cases ɛ =.5,.6,.75,.825,.95,.1,.175,.35, and.5 are integrated numerically. The outcomes are displayed in Tables and Figures
23 ɛ=.5 ɛ= Table 21: ɛ=.5 and ɛ=.6. T r is the computed value of the recurrence time with h=.2 and k= Figure 2: The left figure is the recurrence result for ɛ=.5. The right figure is the recurrence result for ɛ=.6. ɛ=.75 ɛ= Table 22: ɛ=.75 and ɛ=.825. T r is the computed value of the recurrence time with h =.2 and k = Figure 21: The left figure is the recurrence result for ɛ=.75. The right figure is the recurrence result for ɛ=
24 ɛ=.95 ɛ= Table 23: ɛ=.95 and ɛ=.1. T r is the computed value of the recurrence time with h=.2 and k= Figure 22: The left figure is the recurrence result for ɛ=.95. The right figure is the recurrence result for ɛ=.1. ɛ=.175 ɛ= Table 24: ɛ=.175 and ɛ=.35. T r is the computed value of the recurrence time with h=.2 and k=.1. The data is summarized in Table 26. Using least squares fitting, there obtains the equation: t r = 4.84(1ɛ) (38) This is plotted with the data from Table 26 in Figure
25 Figure 23: The left figure is the recurrence result for ɛ=.175. The right figure is the recurrence result for ɛ=.35. T r shape difference Table 25: ɛ=.5. T r is the computed value of the recurrence time with h =.2 and k = amplitude t=.96 Figure 24: The left figure is the recurrence result for ɛ=.5. The right figure is the solution at the first recurrence time. 25
26 1ɛ t r Table 26: The relationship between the first recurrence time of GBBM and various values of ɛ. 8 y = 4.84*x**(.841) The fisrt recurrence time ε Figure 25: The various ɛ and their corresponding first recurrence times. 26
27 5 Conclusions For the KdV- and mkdv-equations, recurrence is connected to the property of resolution into solitary waves together with the exact interaction of solitary waves. Studies not reported here show that we do have resolution into solitary waves for the GKdV-equation with p = 3 and for the GBBM-equations. However, it is also known that the solitary-wave solutions of these equations do not interact exactly. Rather, a very small amount of energy is lost to dispersion each time they interact. However, the numerical results presented herein support the notion that these solitary-wave solutions obtain near recurrence of initial states. Because of dissipative energy loss, we expect the recurrence to eventually weaken as less energy is concentrated in solitary waves. Because of the very small relative loss of energy, this decay of recurrence is expected to occur on much longer time scales than those appearing here. Indeed, as far as modelling real wave phenomena is concerned, the time scale for decay of recurrence appears considerably longer than the time scale of formal validity of the model (see e.g. [5]). References [1] M. J. Ablowitz and H. Segur, Solitons and the Inverse Scattering Transform, SIAM Studies in Appl. Math.: Philadelphia, [2] V. I. Arnold, Mathematical methods of classical mechanics, Springer- Verlag: Berlin, [3] T. B. Benjamin, J. L. Bona and J. J. Mahony, Model equations for long waves in nonlinear dispersive systems, Phil. Trans. Royal Soc. London. Ser. A 272 (1972) [4] J. L. Bona, Convergence of periodic wavetrains in the limit of large wavelength, Appl. Sci. Res. 37 (1981) [5] J. L. Bona, T. Colin and D. Lannes, Long wave approximations for water waves, preprint (24). [6] J. L. Bona, V. A. Dougalis and O. A. Karakashian, Fully discrete Galerkin methods for the Korteweg-de Vries equation, J. Comp. Math. Applic. 12A (1986) [7] J. L. Bona, V. A. Dougalis, O. A. Karakashian and W. R. McKinney, Conservative, high-order numerical schemes for the generalized Kortewegde Vries equation, Philos. Trans. Royal Soc. London. Ser. A 351 (1995) [8] J. L. Bona, W. R. McKinney and J. M. Restrepo, Stable and unstable solitary-wave solutions of the generalized regularized long-wave equation, J. Nonlinear Science 1 (2)
28 [9] J. L. Bona, W. G. Pritchard and L. R. Scott, Solitary-wave interaction, Phys. Fluids 23 (198) [1] J. L. Bona, W. G. Pritchard and L. R. Scott, A comparison of solutions of two model equations for long waves, Lectures in Applied Mathematics vol. 2 (ed. N. Lebovitz) (1983) [11] J. L. Bona, P. E. Souganidis and W. A. Strauss, Stability and instability of solitary waves of KdV type, Proc. Royal Soc. London Ser. A 411 (1987) [12] J. Bourgain, Fourier transform restriction phenomena for certain lattice subsets and applications to nonlinear equations, part II; the KdV equation, Geom. Func. Anal. 3 (1993) [13] J. Boussinesq, Théorie de l intumescence liquid appelée onde solitaire ou de translation, se propageant dans un canal rectangulaire, Comptes Rendus Acad. Sci. 72 (1871) [14] R. Camassa and L. Lee, Complete integrable particle methods and the recurrence of initial states for a nonlinear shallow-water wave equation, J. Comp. Phys. 227 (28) [15] P. G. Drazin and R. S. Johnson, Solitons: an introduction, Cambridge Texts in Appl. Math. Cambridge University Press, [16] E. Fermi, J. Pasta, and S. M. Ulam, Studies in nonlinear problems, Los Alamos Report LA194 (1955), also in Nonlinear Wave Motion, Lectures in Applied Math vol. 15 (1974) pp (ed. A.C. Newell) American Math. Soc.: Providence. [17] C. S. Gardner, J. M. Greene, M. D. Kruskal and R. M. Miura, Method for solving the Korteweg-de Vries equation, Phys. Rev. Lett. 19 (1967) [18] C. S. Gardner, J. M. Greene, M. D. Kruskal and R. M. Miura, The Korteweg-de Vries equation and generalizations (VI): Methods for exact solutions, Comm. Pure Appl. Math. 27 (1974) [19] B.-Y. Guo and V. Manoranjan, A spectral method for solving the RLW equation, IMA J. Num. Anal. 5 (1985) [2] P. D. Lax, Integrals of nonlinear equations of evolution and solitary waves, Comm. Pure Appl. Math. 21 (1968) [21] P. D. Lax, Periodic solutions of the KdV equation, Comm. Pure Appl. Math. 28 (1975) [22] H. P. McKean and E. Trubowitz, Hill s operator and hyperelliptic function theory in the presence of infinitely many branch points, Comm. Pure Appl. Math. 29 (1976)
29 [23] H. P. McKean and P. van Moerbeke, The spectrum of Hill s equation, Inventiones Math 3 (1975) [24] J. B. McLeod and P. J. Olver, The connection between partial differential equations solvable by inverse scattering and ordinary differential equations of Painlevé type, SIAM J. Math. Anal. 14 (1983) [25] Y. Martel and F. Merle, Instability of solitons for the critical generalized Korteweg-de Vries equation, Geom. Funct. Anal. 11 (21) [26] R. M. Miura, C. S. Gardner and M. D. Kruskal, Korteweg-de Vries equation and generalizations II. Existence of conservation laws and constants of motion, J. Math. Phys. 9 (1968) [27] S. P. Novikov, The periodic problem for the Korteweg-de Vries equation, Funct. Anal. Appl. 8 (1974) [28] P. J. Olver, Euler operators and conservation laws of the BBM equation, Math. Proc. Cambridge Philos. Soc. 85 (1979) [29] J. E. Pasciak, Spectral methods for a nonlinear initial value problem involving pseudo differential operators, SIAM J. Numer. Anal. 19 (1982) [3] J. S. Russell, Report on waves, Rep. 14th Meet. Brit. Assoc. Adv. Sci. (1844) [31] D. Sloan, Fourier pseudospectral solution of the regularized long wave equation, Comp. Appl. Math. 36 (1996) [32] P. E. Souganidis and W. A. Strauss, Instability of a class of dispersive solitary waves, Proc. Royal Soc. Edinburgh 114A (199) [33] F.D. Tappert, Numerical solutions of the Korteweg-de Vries equation and its generalizations by the split-step Fourier method, in Nonlinear Wave Motion, Lecture Notes in Applied Math. vol. 15 (1974) , American. Math. Soc. : Providence R. I. [34] M. Toda, Studies of a nonlinear lattice, Phys. Rep. 18 (1975) [35] G. B. Whitham, Linear and Nonlinear Waves, Wiley: New York, [36] J.-M. Yuan, Studies in recurrence and singularity formation in nonlinear dispersive wave equations, Ph.D. Thesis, University of Texas at Austin, 21. [37] N. J. Zabusky, Nonlinear lattice dynamics and energy sharing, J. Phys. Soc. Japan 26 (1968) [38] N. J. Zabusky and M. D. Kruskal, Interactions of solitons in a collisionless plasma and the recurrence of initial states, Phys. Rev. Lett. 15 (1965)
NUMERICAL SOLITARY WAVE INTERACTION: THE ORDER OF THE INELASTIC EFFECT
ANZIAM J. 44(2002), 95 102 NUMERICAL SOLITARY WAVE INTERACTION: THE ORDER OF THE INELASTIC EFFECT T. R. MARCHANT 1 (Received 4 April, 2000) Abstract Solitary wave interaction is examined using an extended
More informationB-splines Collocation Algorithms for Solving Numerically the MRLW Equation
ISSN 1749-889 (print), 1749-897 (online) International Journal of Nonlinear Science Vol.8(2009) No.2,pp.11-140 B-splines Collocation Algorithms for Solving Numerically the MRLW Equation Saleh M. Hassan,
More informationDerivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations
Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations H. A. Erbay Department of Natural and Mathematical Sciences, Faculty of Engineering, Ozyegin University, Cekmekoy 34794,
More informationWhat Is a Soliton? by Peter S. Lomdahl. Solitons in Biology
What Is a Soliton? by Peter S. Lomdahl A bout thirty years ago a remarkable discovery was made here in Los Alamos. Enrico Fermi, John Pasta, and Stan Ulam were calculating the flow of energy in a onedimensional
More informationNumerical Methods for Generalized KdV equations
Numerical Methods for Generalized KdV equations Mauricio Sepúlveda, Departamento de Ingeniería Matemática, Universidad de Concepción. Chile. E-mail: mauricio@ing-mat.udec.cl. Octavio Vera Villagrán Departamento
More informationFORCED OSCILLATIONS OF A CLASS OF NONLINEAR DISPERSIVE WAVE EQUATIONS AND THEIR STABILITY
Jrl Syst Sci & Complexity (2007) 20: 284 292 FORCED OSCILLATIONS OF A CLASS OF NONLINEAR DISPERSIVE WAVE EQUATIONS AND THEIR STABILITY Muhammad USMAN Bingyu ZHANG Received: 14 January 2007 Abstract It
More informationExact and approximate nonlinear waves generated by the periodic superposition of solitons
Journal of Applied Mathematics and Physics (ZAMP) 0044-2275/89/060940-05 $ 1.50 + 0.20 Vol. 40, November 1989 9 1989 Birkhguser Verlag, Basel Exact and approximate nonlinear waves generated by the periodic
More informationNonlinear Fourier Analysis
Nonlinear Fourier Analysis The Direct & Inverse Scattering Transforms for the Korteweg de Vries Equation Ivan Christov Code 78, Naval Research Laboratory, Stennis Space Center, MS 99, USA Supported by
More informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationRecurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University
Nonlinear Physics: Theory and Experiment IV Gallipoli, Lecce, Italy 27 June, 27 Recurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University david.trubatch@usma.edu
More informationExact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially Integrable Equations
Thai Journal of Mathematics Volume 5(2007) Number 2 : 273 279 www.math.science.cmu.ac.th/thaijournal Exact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationChebyshev Collocation Spectral Method for Solving the RLW Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.7(2009) No.2,pp.131-142 Chebyshev Collocation Spectral Method for Solving the RLW Equation A. H. A. Ali Mathematics
More informationarxiv: v1 [physics.flu-dyn] 14 Jun 2014
Observation of the Inverse Energy Cascade in the modified Korteweg de Vries Equation D. Dutykh and E. Tobisch LAMA, UMR 5127 CNRS, Université de Savoie, Campus Scientifique, 73376 Le Bourget-du-Lac Cedex,
More informationDynamics of Solitary Waves Induced by Shock Impulses in a Linear Atomic Chain*
Dynamics of Solitary Waves Induced by Shock Impulses in a Linear Atomic Chain* PHUOC X. TRAN, DONALD W. BRENNER, and C. T. WHITE Naval Research Laboratory, Washington, DC 20375-5000 Abstract The propagation
More informationFreedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation
Freedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation Alex Veksler 1 and Yair Zarmi 1, Ben-Gurion University of the Negev, Israel 1 Department
More informationA Note On Solitary Wave Solutions of the Compound Burgers-Korteweg-de Vries Equation
A Note On Solitary Wave Solutions of the Compound Burgers-Korteweg-de Vries Equation arxiv:math/6768v1 [math.ap] 6 Jul 6 Claire David, Rasika Fernando, and Zhaosheng Feng Université Pierre et Marie Curie-Paris
More informationWhat we do understand by the notion of soliton or rather solitary wave is a wave prole which is very stable in the following sense
Introduction to Waves and Solitons What we do understand by the notion of soliton or rather solitary wave is a wave prole which is very stable in the following sense i) It is localized, which means it
More informationAPPROXIMATE MODEL EQUATIONS FOR WATER WAVES
COMM. MATH. SCI. Vol. 3, No. 2, pp. 159 170 c 2005 International Press APPROXIMATE MODEL EQUATIONS FOR WATER WAVES RAZVAN FETECAU AND DORON LEVY Abstract. We present two new model equations for the unidirectional
More information2. The generalized Benjamin- Bona-Mahony (BBM) equation with variable coefficients [30]
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.12(2011) No.1,pp.95-99 The Modified Sine-Cosine Method and Its Applications to the Generalized K(n,n) and BBM Equations
More informationSoliton Solutions of a General Rosenau-Kawahara-RLW Equation
Soliton Solutions of a General Rosenau-Kawahara-RLW Equation Jin-ming Zuo 1 1 School of Science, Shandong University of Technology, Zibo 255049, PR China Journal of Mathematics Research; Vol. 7, No. 2;
More informationarxiv: v1 [nlin.ps] 18 Sep 2008
Asymptotic two-soliton solutions solutions in the Fermi-Pasta-Ulam model arxiv:0809.3231v1 [nlin.ps] 18 Sep 2008 Aaron Hoffman and C.E. Wayne Boston University Department of Mathematics and Statistics
More informationSymmetry reductions and travelling wave solutions for a new integrable equation
Symmetry reductions and travelling wave solutions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX 0, 50 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es
More informationDispersive and soliton perturbations of finite-genus solutions of the KdV equation: computational results
Dispersive and soliton perturbations of finite-genus solutions of the KdV equation: computational results Thomas Trogdon and Bernard Deconinck Department of Applied Mathematics University of Washington
More informationarxiv:nlin/ v1 [nlin.si] 25 Sep 2006
Remarks on the conserved densities of the Camassa-Holm equation Amitava Choudhuri 1, B. Talukdar 1a and S. Ghosh 1 Department of Physics, Visva-Bharati University, Santiniketan 73135, India Patha Bhavana,
More informationContinuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China
Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China the 3th GCOE International Symposium, Tohoku University, 17-19
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationHOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (21), 89 98 HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS Hossein Jafari and M. A. Firoozjaee Abstract.
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationFourier spectral methods for solving some nonlinear partial differential equations
Int. J. Open Problems Compt. Math., Vol., No., June ISSN 99-; Copyright ICSRS Publication, www.i-csrs.org Fourier spectral methods for solving some nonlinear partial differential equations Hany N. HASSAN
More informationBasic Concepts and the Discovery of Solitons p. 1 A look at linear and nonlinear signatures p. 1 Discovery of the solitary wave p. 3 Discovery of the
Basic Concepts and the Discovery of Solitons p. 1 A look at linear and nonlinear signatures p. 1 Discovery of the solitary wave p. 3 Discovery of the soliton p. 7 The soliton concept in physics p. 11 Linear
More informationDISSIPATIVE INITIAL BOUNDARY VALUE PROBLEM FOR THE BBM-EQUATION
Electronic Journal of Differential Equations, Vol. 8(8), No. 149, pp. 1 1. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) DISSIPATIVE
More informationLecture 10: Whitham Modulation Theory
Lecture 10: Whitham Modulation Theory Lecturer: Roger Grimshaw. Write-up: Andong He June 19, 2009 1 Introduction The Whitham modulation theory provides an asymptotic method for studying slowly varying
More informationA MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE
Pacific Journal of Applied Mathematics Volume 1, Number 2, pp. 69 75 ISSN PJAM c 2008 Nova Science Publishers, Inc. A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Wen-Xiu Ma Department
More informationGeneration of undular bores and solitary wave trains in fully nonlinear shallow water theory
Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK
More informationBreaking soliton equations and negative-order breaking soliton equations of typical and higher orders
Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID
More informationINVERSE SCATTERING TRANSFORM, KdV, AND SOLITONS
INVERSE SCATTERING TRANSFORM, KdV, AND SOLITONS Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762, USA Abstract: In this review paper, the
More informationApplication of linear combination between cubic B-spline collocation methods with different basis for solving the KdV equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 016, pp. 191-04 Application of linear combination between cubic B-spline collocation methods with different basis
More informationSolitons : An Introduction
Solitons : An Introduction p. 1/2 Solitons : An Introduction Amit Goyal Department of Physics Panjab University Chandigarh Solitons : An Introduction p. 2/2 Contents Introduction History Applications Solitons
More informationRelation between Periodic Soliton Resonance and Instability
Proceedings of Institute of Mathematics of NAS of Ukraine 004 Vol. 50 Part 1 486 49 Relation between Periodic Soliton Resonance and Instability Masayoshi TAJIRI Graduate School of Engineering Osaka Prefecture
More informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationThe Traveling Wave Solutions for Nonlinear Partial Differential Equations Using the ( G. )-expansion Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.8(009) No.4,pp.435-447 The Traveling Wave Solutions for Nonlinear Partial Differential Equations Using the ( )-expansion
More informationExact solutions through symmetry reductions for a new integrable equation
Exact solutions through symmetry reductions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX, 1151 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es
More informationPrimary Results on Simulation of a New Kind Traveling Wave. of Permanent Form
Primary Results on Simulation of a New Kind Traveling Wave of Permanent Form Min CHEN y Abstract This talk concerns the nonlinear dispersive waves in a water channel. Results from the study of a new version
More informationAn integrable shallow water equation with peaked solitons
An integrable shallow water equation with peaked solitons arxiv:patt-sol/9305002v1 13 May 1993 Roberto Camassa and Darryl D. Holm Theoretical Division and Center for Nonlinear Studies Los Alamos National
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationExact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Published by Canadian Center of Science and Education Exact Solutions for the Nonlinear +-Dimensional Davey-Stewartson Equation
More informationMATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS
MATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS Poisson Systems and complete integrability with applications from Fluid Dynamics E. van Groesen Dept. of Applied Mathematics University oftwente
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
More informationSome Soliton Solutions of Non Linear Partial Differential Equations by Tan-Cot Method
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578,p-ISSN: 319-765X, 6, Issue 6 (May. - Jun. 013), PP 3-8 Some Soliton Solutions of Non Linear Partial Differential Equations by Tan-Cot Method Raj Kumar
More informationShock Waves & Solitons
1 / 1 Shock Waves & Solitons PDE Waves; Oft-Left-Out; CFD to Follow Rubin H Landau Sally Haerer, Producer-Director Based on A Survey of Computational Physics by Landau, Páez, & Bordeianu with Support from
More information12 Lecture 12. Lax pair and Toda lattice
128 12 Lecture 12. Lax pair and Toda lattice The story in this lecture goes as follows chronologically. Scott-Russel observed (1834) the first solitary wave propagation along the Union Canal, Scotland
More informationGeneration of Secondary Solitary Waves in the Variable-Coefficient Korteweg-de Vries Equation
Loughborough University Institutional Repository Generation of Secondary Solitary Waves in the Variable-Coefficient Korteweg-de Vries Equation This item was submitted to Loughborough University's Institutional
More informationThe superposition of algebraic solitons for the modified Korteweg-de Vries equation
Title The superposition of algebraic solitons for the modified Korteweg-de Vries equation Author(s) Chow, KW; Wu, CF Citation Communications in Nonlinear Science and Numerical Simulation, 14, v. 19 n.
More informationThe General Form of Linearized Exact Solution for the KdV Equation by the Simplest Equation Method
Applied and Computational Mathematics 015; 4(5): 335-341 Published online August 16 015 (http://www.sciencepublishinggroup.com/j/acm) doi: 10.11648/j.acm.0150405.11 ISSN: 38-5605 (Print); ISSN: 38-5613
More informationTHE LAX PAIR FOR THE MKDV HIERARCHY. Peter A. Clarkson, Nalini Joshi & Marta Mazzocco
Séminaires & Congrès 14, 006, p. 53 64 THE LAX PAIR FOR THE MKDV HIERARCHY by Peter A. Clarkson, Nalini Joshi & Marta Mazzocco Abstract. In this paper we give an algorithmic method of deriving the Lax
More informationHamiltonian partial differential equations and Painlevé transcendents
Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Cauchy problem for evolutionary PDEs with
More informationAvailable online at J. Math. Comput. Sci. 2 (2012), No. 1, ISSN:
Available online at http://scik.org J. Math. Comput. Sci. 2 (2012), No. 1, 15-22 ISSN: 1927-5307 BRIGHT AND DARK SOLITON SOLUTIONS TO THE OSTROVSKY-BENJAMIN-BONA-MAHONY (OS-BBM) EQUATION MARWAN ALQURAN
More informationCONVERGENCE OF SOLITARY-WAVE SOLUTIONS IN A PERTURBED BI-HAMILTONIAN DYNAMICAL SYSTEM. I. COMPACTONS AND PEAKONS.
CONVERGENCE OF SOLITARY-WAVE SOLUTIONS IN A PERTURBED BI-HAMILTONIAN DYNAMICAL SYSTEM. I. COMPACTONS AND PEAKONS. Y. A. Li 1 and P. J. Olver 1, Abstract. We investigate how the non-analytic solitary wave
More informationKink, singular soliton and periodic solutions to class of nonlinear equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic
More informationCHAPTER 1 INTRODUCTION. 1.1 KORTEWEG - de VRIES EQUATION. It is well known that the celebrated Korteweg-de Vries. (KdV) equation.
CHAPTER 1 INTRODUCTION 1.1 KORTEWEG - de VRIES EQUATION It is well known that the celebrated Korteweg-de Vries (KdV) equation U + 6UU t x xxx is the simplest model equation describing a nonlinear dispersive
More informationKeywords: Exp-function method; solitary wave solutions; modified Camassa-Holm
International Journal of Modern Mathematical Sciences, 2012, 4(3): 146-155 International Journal of Modern Mathematical Sciences Journal homepage:www.modernscientificpress.com/journals/ijmms.aspx ISSN:
More informationModified Simple Equation Method and its Applications for some Nonlinear Evolution Equations in Mathematical Physics
Modified Simple Equation Method and its Applications for some Nonlinear Evolution Equations in Mathematical Physics Elsayed M. E. Zayed Mathematics department, Faculty of Science Zagazig University, Zagazig,
More informationNumerical Study of Oscillatory Regimes in the KP equation
Numerical Study of Oscillatory Regimes in the KP equation C. Klein, MPI for Mathematics in the Sciences, Leipzig, with C. Sparber, P. Markowich, Vienna, math-ph/"#"$"%& C. Sparber (generalized KP), personal-homepages.mis.mpg.de/klein/
More informationAre Solitary Waves Color Blind to Noise?
Are Solitary Waves Color Blind to Noise? Dr. Russell Herman Department of Mathematics & Statistics, UNCW March 29, 2008 Outline of Talk 1 Solitary Waves and Solitons 2 White Noise and Colored Noise? 3
More informationQuasipatterns in surface wave experiments
Quasipatterns in surface wave experiments Alastair Rucklidge Department of Applied Mathematics University of Leeds, Leeds LS2 9JT, UK With support from EPSRC A.M. Rucklidge and W.J. Rucklidge, Convergence
More informationNew Exact Travelling Wave Solutions for Regularized Long-wave, Phi-Four and Drinfeld-Sokolov Equations
ISSN 1749-3889 print), 1749-3897 online) International Journal of Nonlinear Science Vol.008) No.1,pp.4-5 New Exact Travelling Wave Solutions for Regularized Long-wave, Phi-Four and Drinfeld-Sokolov Euations
More informationNotes on the Inverse Scattering Transform and Solitons. November 28, 2005 (check for updates/corrections!)
Notes on the Inverse Scattering Transform and Solitons Math 418 November 28, 2005 (check for updates/corrections!) Among the nonlinear wave equations are very special ones called integrable equations.
More informationDeterminant Expressions for Discrete Integrable Maps
Typeset with ps2.cls Full Paper Determinant Expressions for Discrete Integrable Maps Kiyoshi Sogo Department of Physics, School of Science, Kitasato University, Kanagawa 228-8555, Japan Explicit
More informationNUMERICAL SIMULATION OF NONLINEAR DISPERSIVE QUANTIZATION
NUMERICAL SIMULATION OF NONLINEAR DISPERSIVE QUANTIZATION GONG CHEN AND PETER J. OLVER Abstract. When posed on a periodic domain in one space variable, linear dispersive evolution equations with integral
More informationA new integrable system: The interacting soliton of the BO
Phys. Lett., A 204, p.336-342, 1995 A new integrable system: The interacting soliton of the BO Benno Fuchssteiner and Thorsten Schulze Automath Institute University of Paderborn Paderborn & Germany Abstract
More informationCollisionless Shocks and the Earth s Bow Shock
Collisionless Shocks and the Earth s Bow Shock Jean-Luc Thiffeault AST 381 Gas Dynamics 17 November 1994 1 Introduction The mean free path for particle collisions in the solar wind at the boundary of the
More informationMultisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system
arxiv:407.7743v3 [math-ph] 3 Jan 205 Multisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system L. Cortés Vega*, A. Restuccia**, A. Sotomayor* January 5,
More informationMultiple-Soliton Solutions for Extended Shallow Water Wave Equations
Studies in Mathematical Sciences Vol. 1, No. 1, 2010, pp. 21-29 www.cscanada.org ISSN 1923-8444 [Print] ISSN 1923-8452 [Online] www.cscanada.net Multiple-Soliton Solutions for Extended Shallow Water Wave
More informationA Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion
Journal of Computational Physics 155, 456 467 (1999) Article ID jcph.1999.6351, available online at http://www.idealibrary.com on A Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion
More informationTravelling Wave Solutions for the Gilson-Pickering Equation by Using the Simplified G /G-expansion Method
ISSN 1749-3889 (print, 1749-3897 (online International Journal of Nonlinear Science Vol8(009 No3,pp368-373 Travelling Wave Solutions for the ilson-pickering Equation by Using the Simplified /-expansion
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationComputing Time-Periodic Solutions of Nonlinear Systems of Partial Differential Equations
1 Computing ime-periodic Solutions of Nonlinear Systems of Partial Differential Equations David M. Ambrose Department of Mathematics, Drexel University, Philadelphia, PA 191 USA E-mail: ambrose@math.drexel.edu
More informationOn the Whitham Equation
On the Whitham Equation Henrik Kalisch Department of Mathematics University of Bergen, Norway Joint work with: Handan Borluk, Denys Dutykh, Mats Ehrnström, Daulet Moldabayev, David Nicholls Research partially
More informationOn the N-tuple Wave Solutions of the Korteweg-de Vnes Equation
Publ. RIMS, Kyoto Univ. 8 (1972/73), 419-427 On the N-tuple Wave Solutions of the Korteweg-de Vnes Equation By Shunichi TANAKA* 1. Introduction In this paper, we discuss properties of the N-tuple wave
More informationarxiv:math/ v1 [math.ap] 1 Jan 1992
APPEARED IN BULLETIN OF THE AMERICAN MATHEMATICAL SOCIETY Volume 26, Number 1, Jan 1992, Pages 119-124 A STEEPEST DESCENT METHOD FOR OSCILLATORY RIEMANN-HILBERT PROBLEMS arxiv:math/9201261v1 [math.ap]
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationINTERACTIONS OF SOLITARY WAVES IN HIERARCHICAL KdV-TYPE SYSTEM Research Report Mech 291/08
INTERACTIONS OF SOLITARY WAVES IN HIERARCHICAL KdV-TYPE SYSTEM Research Report Mech 291/8 Lauri Ilison, Andrus Salupere Centre for Nonlinear Studies, Institute of Cybernetics at Tallinn University of Technology,
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationLecture Notes on the Fermi-Pasta-Ulam-Tsingou Problem
Lecture Notes on the Fermi-Pasta-Ulam-Tsingou Problem Andrew Larkoski November 9, 016 Because this is a short week, we re going to just discuss some fun aspects of computational physics. Also, because
More informationA Novel Nonlinear Evolution Equation Integrable by the Inverse Scattering Method
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part, 384 39 A Novel Nonlinear Evolution Equation Integrable by the Inverse Scattering Method Vyacheslav VAKHNENKO and John PARKES
More informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationGevrey regularity in time for generalized KdV type equations
Gevrey regularity in time for generalized KdV type equations Heather Hannah, A. Alexandrou Himonas and Gerson Petronilho Abstract Given s 1 we present initial data that belong to the Gevrey space G s for
More informationNumerical solution of General Rosenau-RLW Equation using Quintic B-splines Collocation Method
Available online at www.ispacs.com/cna Volume 2012, Year 2012 Article ID cna-00129, 16 pages doi:10.5899/2012/cna-00129 Research Article Numerical solution of General Rosenau-RLW Equation using Quintic
More informationZero dispersion and viscosity limits of invariant manifolds for focusing nonlinear Schrödinger. equations
J. Math. Anal. Appl. 315 (2006) 642 655 www.elsevier.com/locate/jmaa Zero dispersion and viscosity limits of invariant manifolds for focusing nonlinear Schrödinger equations Y. Charles Li Department of
More informationLecture17: Generalized Solitary Waves
Lecture17: Generalized Solitary Waves Lecturer: Roger Grimshaw. Write-up: Andrew Stewart and Yiping Ma June 24, 2009 We have seen that solitary waves, either with a pulse -like profile or as the envelope
More informationBuilding Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients
Journal of Physics: Conference Series OPEN ACCESS Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients To cite this article: M Russo and S R Choudhury 2014
More informationKdV equation obtained by Lie groups and Sturm-Liouville problems
KdV equation obtained by Lie groups and Sturm-Liouville problems M. Bektas Abstract In this study, we solve the Sturm-Liouville differential equation, which is obtained by using solutions of the KdV equation.
More informationSoliton solutions of some nonlinear evolution equations with time-dependent coefficients
PRAMANA c Indian Academy of Sciences Vol. 80, No. 2 journal of February 2013 physics pp. 361 367 Soliton solutions of some nonlinear evolution equations with time-dependent coefficients HITENDER KUMAR,
More informationExtended Jacobi Elliptic Function Expansion Method for Nonlinear Benjamin-Bona-Mahony Equations
International Mathematical Forum, Vol. 7, 2, no. 53, 239-249 Extended Jacobi Elliptic Function Expansion Method for Nonlinear Benjamin-Bona-Mahony Equations A. S. Alofi Department of Mathematics, Faculty
More informationTHE TODA LATTICE. T (y) = α(e βy 1)
THE TODA LATTICE Background The Toda lattice is named after Morikazu Toda, who discovered in the 1960 s that the differential equations for a lattice with internal force T (y) = α(e βy 1) (with α, β constant)
More informationA small dispersion limit to the Camassa Holm equation: A numerical study
Available online at www.sciencedirect.com Mathematics and Computers in Simulation 80 (2009) 120 130 A small dispersion limit to the Camassa Holm equation: A numerical study Jennifer Gorsky a,, David P.
More informationEnergy transfer model and large periodic boundary value problem for the quintic NLS
Energy transfer model and large periodic boundary value problem for the quintic NS Hideo Takaoka Department of Mathematics, Kobe University 1 ntroduction This note is based on a talk given at the conference
More information