Fourier spectral methods for solving some nonlinear partial differential equations
|
|
- Jewel Potter
- 6 years ago
- Views:
Transcription
1 Int. J. Open Problems Compt. Math., Vol., No., June ISSN 99-; Copyright ICSRS Publication, Fourier spectral methods for solving some nonlinear partial differential equations Hany N. HASSAN Hassan K. SALEH Department of Basic Science, Benha Faculty of Engineering, Benha University, Benha, 5, Egypt Department of Engineering Mathematics and Physics, Faculty of Engineering, Cairo University, Giza, Egypt s, : h_nasr77@yahoo.com, hassan.kamel@yahoo.com Abstract The spectral collocation or pseudospectral (PS) methods (Fourier transform methods) combined with temporal discretization techniques to numerically compute solutions of some partial differential equations (PDEs). In this paper, we solve the Korteweg-de Vries (KdV) equation using a Fourier spectral collocation method to discretize the space variable, leap frog and classical fourth-order Runge-Kutta scheme (RK) for time dependence. Also, Boussinesq equation is solving by a Fourier spectral collocation method to discretize the space variable, finite difference and classical fourth-order Runge- Kutta scheme (RK) for time dependence. Our implementation employs the Fast Fourier Transform (FFT) algorithm. Keywords: Fourier spectral method, Fast Fourier transform, Boussinesq equation, Korteweg-de Vries equation; leap frog, finite difference. Introduction Let us consider the equation u u u u, () t The equation () is known as the Korteweg-de Vries (KdV) equation (with non linearity uu, and dispersion u ) with real constants,. The KdV equation is one of the popular soliton equations. More than a century later the KdV equation
2 Hany N. HASSAN et al. 5 has been found in many other applications such as magnetohydro dynamic waves in a cold plasma, longitudinal vibrations of an enharmonic discrete-mass string, ion-acoustic waves in a cold plasma, pressure waves in liquids-gas bubble miture, rotating flow down a tube, and longitudinal dispersive waves in elastic rods. The eact solutions of () is c c u, t sec h ( ct ), () where is an arbitrary integration constant. We investigate the numerical solution of the KdV equation using the Fourier-Leap Frog methods, due to Fornberg and Whitham [], and the Fourier based fourthorder Runge-Kutta (RK) method for better stability of the solution. Consider the equation utt u ( u ) u, () The equation () is known as the Boussinesq equation, where subscripts and t denote differentiation, was introduced to describe the propagation of long waves in shallow water. The Boussinesq equation also arises in several other physical applications including one-dimensional nonlinear lattice waves, vibrations in a nonlinear string, and ion sound waves in a plasma. It is well known that The Boussinesq equation () has a bidirectional solitary wave solution u(, t) b sec h ( b( ct)). () representing a solitary wave, where c b is the propagation speed and b, arbitrary constants determining the height and the position of the maimum height of the wave, respectively. From the form of c it is apparent that the solution can propagate in either direction (left or right). In the present work, we have applied to () two numerical methods to study soliton solutions and investigate their interactions upon collision: a combination of finite differences and a Fourier pseudospectral method Analysis the method for the Korteweg-de Vries (KdV) equation Now, consider the standard KdV equation in the form ut u u u, [ L, L ], (5) where u = u(, t), subscripts and t denote differentiation, where L is a given number, usually large. We have changed the solution interval from [ L, L ] to [,π], with the change of variable. L Thus, the equation (5) becomes
3 Hany N. HASSAN et al. ( ) ut uu u L L, [, ], () u(, t) is transformed into Fourier space with respect to, and derivatives (or other operators) with respect to. Applying the inverse Fourier transform n u n F {( ik) F( u)}, n n,,. (7) using this with n = and n =, u { ( )} F ikf u, u F { ik F( u)}. The equation () becomes ut uf { ikf( u)} ( ) F { ik F( u)}, L L () In practice, we need to discretize the equation (). For any integer N >, we consider j j j, j,, N, N. Let u (, t) be the solution of the KdV equation (5). Then, we transform it into the discrete Fourier space as N ik N N j uˆ( k, t) F( u) u( j, t) e, k. N (9) J From this, using the inversion formula, we get ( N, ) ( ˆ) ˆ j (, ) ik j,. kn u t F u u k t e j N () where we denote the discrete Fourier transform and the inverse Fourier transform by F and F, respectively. Now, replacing F and F in () by their discrete counterparts, and discretizing () gives du( j, t) u( j, t) F { ikf( u)} ( ) F { ik F( u)}, j N. () dt L L T T [ u(, t), u(, t),, u( N, t)]. Letting U = [ u(, t ), u(, t ),, u( N, t )]. The equation () can be written in the vector form U t = F (U) () where F defines the right hand side of ().. Fourier Leap-Frog Method for KdV Equation A time integration known as a Leap-Frog method (a two step scheme) is given by n n n U U t F where the superscripts denote the time level at which the term is evaluated. In this subsection, we use a Fourier collocation method and this scheme to solve the
4 Hany N. HASSAN et al. 7 system given in () numerically. Here, we follow the analysis of Fornberg and Whitham []. Using the Leap-Frog scheme to advance in time, we obtain U(t +Δt) = U(t -Δt)+ΔtF(U(t)) This is called the Fourier Leap-Frog (FLF) scheme for the KdV equation (5). FLF needs two levels of initial data. Usually, the first one is given by the initial condition u(, ) and the second level u(, Δt) can be obtained by using a higherorder one-step method, for eample, a fourth order Runge- Kutta method. Thus, the time discretization for () is given by u(, t t) u(, t t) t u( j, t) F { ikf( u)} t( ) F { ik F( u)}. () L L where we denote ( ) k. L In general, the Fourier-Leap Frog scheme defined in () is accurate for low enough wave numbers, but it loses accuracy rapidly for increasing wave numbers.. Fourier Based RK Method for KdV Equation Notice that the Fourier-Leap-Frog method is a second-order scheme and has some disadvantage in the way that the solution of the model problem is subject to a temporal oscillation with period Δt a commonly suggested alternative method is the Runge-Kutta methods. The classical fourth-order Runge-Kutta methods for the system () is given by K = F(U n,t n ) K = F(U n + ΔtK,t n + Δt) K = F(U n + ΔtK,t n + Δt) () K = F(U n + ΔtK,t n + Δt) Δt U n+ = U n + K + K + K + K where the superscripts denote the time level at which the term is evaluated.. Linear stability analysis For an analysis of stability we use the standard Fourier analysis to find the condition imposed on the time step Δt. For simplicity we let and consider the KdV equation in the form u u u, (5) t
5 Hany N. HASSAN et al... By using FLF scheme We approimate this equation by u(, t t) u(, t t) tf { i F( u)} tf { i F( u)}. () We look a solution to () of the form t t i u(, t) e Substitution in () gives ( t t ) t i ( t t ) t i t t i t t i e e i t e i t e, i.e. it i t if ( t,, ), where f ( t,, ) t t The scheme is conditionally stable if and only if f ( t,, ) is real and less than one in magnitude. Let us again assume for simplicity that the period is [,π] and that this interval is discretized with N equidistant mesh points, that is N The wave number takes the values,,,, N. We want to find the largest value of t such that f ( t,, ) < is true for all.the most severe restriction on t is imposed for the,which are largest in magnitude, i.e. for ma, ma N. Thus, we obtain f ( t, ma, ) t( ) t( ) If we assume that >, then f ( t,, ) < when < t ( ) So, the stability condition becomes t <.55.. By using RK method To do stability analysis of the RK scheme for the KdV equation, we could use the approach used in analyzing the stability of the Fourier Leap-Frog schemes as t t i in the previous subsection. To do this, we substitute u(, t ) e into ().
6 Hany N. HASSAN et al. 9 After a very tedious and long derivation, we are led to the stability condition or after we numerically eperiment with the scheme, we see that the appropriate time step is the one satisfying the condition t. The right-hand side of the system of ODES in time given in () F (U) is F(U j) u( j, t) F { ikf( u)} F { i k F( u)}, j N. L L. Numerical Results and Eamples In order to show how good the numerical solutions are in comparison with the eact ones, we will use the L and L error norms defined by N eact num eact num L u u u, i ui i L eact num eact num u u ma u u. (7) i To implement the performance of the method, three test problems will be considered: the motion of a single solitary wave, the interaction of two positive solitary waves, the interaction of three positive solitary waves and other solutions... The motion of a single solitary wave Consider the KdV equation (5) with =, = and L =. ut uu u, [ L, L ]. The simplest eact solution is u(, t) sec h ( t). and initial condition u(,) sec h ( ). Eample In this eample we compute the numerical solutions u(, ) using the Fourier Leap-Frog scheme. The numerical solution at t =,.5,,.5 and in Figure. with N =. It is clear from Figure that the proposed method performs the motion of propagation of a solitary wave satisfactorily, which moved to the right with the preserved amplitude. In Figure, we plot the eact solution and Numerical solution at t = with N =, Figure show the error at each collection point at t = with N =, Table displays the values of error norms obtained at N =, N = and N = 5. i i
7 u(,t) u(,t) Hany N. HASSAN et al. 5.5 t = t =.5 t = t =.5 t = Fig. : Fourier spectral solution of the KdV equation using FLF scheme with N = at t =,.5,,.5 and. N t L L Table : Error norms of FLF scheme of the KdV equation at t = with N =, and 5. Numerical Eact Fig. : Fourier spectral solution of the KdV equation using FLF scheme with N = at t =
8 Error Hany N. HASSAN et al Fig. : The error (error = eact numerical) distributions in FLF scheme for the KdV equation with N = at t =. Eample Now, we solve the same problem using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. Figure shows simulation of the solution computed using N =. In Figure 5, we plot the eact solution and Numerical solution at t = with N =, Figure show the error at each collection point at t = with N =. Table displays error norms obtained at N =, N = and N = 5. From these results we can see that by carefully choosing the time steps, RK is more accurate than the Fourier Leap-frog, but Fourier Leapfrog is easier than RK. Fig. : Numerical simulation of Fourier spectral solution of the KdV equation using RK scheme with N =.
9 Error u(,t) Hany N. HASSAN et al. 5 N Δt L L Table : Error norms of Fourier spectral solution of the KdV equation using RK scheme at t = with N =, and 5. Numerical Eact Fig. 5: Fourier spectral solution of the KdV equation using RK scheme at t = with N = Fig. : The error distributions in Fourier spectral solution of the KdV equation using RK scheme at t = with N =... The interaction of two positive solitary waves Consider the KdV equation (5) with =, =, and L =. The initial data is
10 u(,t) u(,t) u(,t) Hany N. HASSAN et al. 5 u(,) n( n )sec h ( ) () results in n solitons solution that propagate with different velocities. The initial for n = is u(,) sech ( ) (9) The eact solution can be epressed as [] cosh( t) cosh( t) u(, t) [cosh( ) cosh( )] t t Eample We solve equation (5) with initial solution (9) using the FLF scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. We plot the eact solution and the numerical solutions at t = -.,-.,,., and. at N = 5 and t. in Figure 7. 7 t = -. Numerical Eact 7 t = -. Numerical Eact (a) (b) 7 Numerical Eact t = (c)
11 u(,t) u(,t) Hany N. HASSAN et al. 5 7 Numerical Eact t =. 7 Numerical Eact t = (d) (e) Fig. 7: Fourier spectral solution for interaction of two waves of the KdV equation using FLF scheme with initial condition (9) and N =5. Eample We solve eample using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. Figure show simulation of the solution computed. We plot the eact solution and the numerical solutions at t = -.,-.,,., and. with N = 5 and t.7 in Figure 9. Fig. : Numerical simulation of Fourier spectral of interaction of two waves of the KdV equation using RK scheme with initial condition (9) and N =5.
12 u(,t) u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al t = -. Numerical Eact 7 t = -. Numerical Eact (a) (b) 5 Numerical Eact t = (c) 7 Numerical Eact t =. 7 Numerical Eact t = (d) (e) Fig. 9: Fourier spectral solution for interaction of two waves of the KdV equation using RK scheme with initial condition (9) and N =5.
13 u(,t) u(,t) Hany N. HASSAN et al. 5 In net method, we can study the interaction of n solitary waves by using the initial condition given by the linear sum of n separate solitary waves of various amplitudes n i i i. () u(,).5c sec h.5c i Interaction of two positive solitary waves can be studied by using the initial condition given by the linear sum of two separate solitary waves of various amplitudes u(,) u (,) u (,) u(,) sec h ( ).5sec h (.5( )) () Eample 5 The calculation is carried out with the time step t =.9 and N = 5 over the region -, we solve using the FLF scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. We plot the numerical solutions at t =,.5,.,. and.5 with N = 5, in Figure, respectively. The initial function was placed on the left side of the region with the larger wave to the left of the smaller one as seen in the Figure, both waves move to the right with velocities dependent upon their magnitudes. The shapes of the two solitary waves is graphed during the interaction at t =. and after the interaction at time t =.5, which is seen to have separated the larger wave from the smaller one. According to the figure, the larger wave catches up the smaller wave at about t =.5, the overlapping process continues until the time t = 5, then two solitary waves emerge from the interaction and resume their former shapes and amplitudes..5 t =.5 t = (a) (b)
14 u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al. 57 t =. t = (c) (d).5 t = 5.5 t = (e) (f) Fig. : Fourier spectral solution for interaction of two waves of the KdV equation using FLF scheme with initial condition () and N = 5. Eample We solve eample 5 using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. The calculation is carried out with the time step t =.9 and N = 5 over the region, Figure shows simulation of the solution computed. We plot the numerical solutions at t =,.5,.,. and.5 with N = 5, in Figure, respectively.
15 u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al. 5 Figure : Numerical simulation of Fourier spectral of interaction of two waves of the KdV equation using RK scheme with initial condition () and N =5..5 t =.5 t = (a) (b).5 t =..5 t = (c) (d)
16 u(,t) u(,t) Hany N. HASSAN et al t = 5.5 t = (e) (f) Fig. : Fourier spectral solution for interaction of two waves of the KdV equation using RK scheme with initial condition () and N =5... The interaction of three positive solitary waves Consider the KdV equation (5) with =, =, and L = ut uu u, [,], The initial for n = in () is u(,) sec h ( ) () Eample 7 We solve (5) with initial condition ().The calculation is carried out with the time step t =. and N = 5 over the region -, we solve using the FLF scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. We plot the numerical solutions at t = -., -., -.5,,.5 and. with N =5, in Figure, respectively. Eample We solve eample 7 using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. The calculation is carried out with the time step t =.7 and N = 5 over the region -, we plot the numerical solutions at t = -., -., -.5,,.5, and. with N =5, in Figure, respectively.
17 u(,t) u(,t) u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al. t = -. t = (a) (b) t = -.5 t = (c) (d) t =.5 t = (e) (f) Fig. : Fourier spectral solution for interaction of three waves of the KdV equation using FLF scheme with initial condition () and N =5.
18 u(,t) u(,t) u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al. t = -. t = (a) (b) t = -.5 t = (b) (c) (d) t =.5 t = (e) (f ) Fig. : Fourier spectral solution for interaction of three waves of the KdV equation using RK scheme with initial condition () and N = 5.
19 u(,t) u(,t) Hany N. HASSAN et al. In this interaction of three positive solitary waves is studied by using the initial condition given by the linear sum of three separate solitary waves of various amplitudes from () u(,) u (,) u (,) u (,) (,) sec ( 5) sec (.5() ( )).5sec (.5( 5)) u h h h () Eample 9 We solve (5) using initial condition ().The calculation is carried out with the time step t =.57 and N = 5 over the region -7 7, we solve using the FLF scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. Figure 5 shows the numerical solution at t =, 5,,,,.As can be seen the three pulses travel with time to the right. But the taller soliton moves faster, hence the three occasionally merge and then split apart again. Eample We solve eample 9 using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. The calculation is carried out with the time step t =. and N = 5 over the region 7 7, we solve using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. Figure and 7 shows simulation and plane view of the solution computed.. t =. t = (a) (b)
20 u(,t) u(,t) u(,t) u(,t) Hany N. HASSAN et al.. t =. t = (c) (d). t =. t = (e) (f) Fig. 5: Fourier spectral solution for interaction of three waves of the KdV equation using FLF scheme with initial condition () and N = 5. Fig. : Numerical simulation of Fourier spectral of interaction of three waves of the KdV equation using RK scheme with initial condition () and N =5.
21 Hany N. HASSAN et al. Fig. 7: Plan view of Fourier spectral of interaction of three waves of the KdV equation using RK scheme with initial condition () and N = 5... Other solutions Now suppose that the initial condition is such that is does not just produce one or more solitons. Let us choose as initial condition [] u(,) sec h ( ) () with =, =, L =,Δt =.9 and N = using the RK scheme to march the solution in time and the Fourier spectral method to take care of the spatial domain. Four is not of the form n (n + ).The peak moving to the right. In addition, there are waves moving to the left. These will disperse and lose their form with time.
22 u(,t) u(,t) u(,t) Hany N. HASSAN et al. 5 Fig. : Numerical simulation of Fourier spectral solution of the KdV equation using RK scheme with initial condition () and N =. 5 t = - - (a) 5 t =. 5 t = (b) (c) Fig. 9: Fourier spectral solution of the KdV equation using RK scheme with initial condition () and N =.
23 Hany N. HASSAN et al. Analysis the method for the Boussinesq equation. Fourier based RK method for Boussinesq equation Consider the Boussinesq equation () in the following form: utt u uu u u, [ a, z] (5) We solve this equation by combination of RK method with respect to t and a Fourier pseudospectral method with respect to. To prepare the equation for numerical solution we introduce the auiliary variable v ut. This reduces the second-order equation in time to the first order system ut v () vt u uu u u We need two initial conditions. The initial conditions we use to numerically solve equation (5) can thus be etracted from the above relation () for t = at u(, t) and ut (, t ) u(,) b sec h ( b( )) v u b c h b b (7) (,) t (,) sec ( ( )) tanh( ( )) We have changed the solution interval from [a, z] to [,π], with the change of variable. a L where L = z a, thus the equation () become u v t () vt u uu u u L L L L u(, t) is transformed into Fourier space with respect to, and derivatives (or other operators) with respect to. Applying the inverse Fourier transform using u F { ikf ( u )} The equation () becomes u F k F u u F { k F( u)} { ( )}
24 Hany N. HASSAN et al. 7 u t v vt F { k F( u)} uf { k F( u)} L L F { ikf( u)} F { k F( u)}, L L In practice, we need to discretize the equation (). For any integer N >, we consider j j, j,,, N. N Let u (, t), v(, t) be the solution equation (9). Then, we transform it into the Discrete Fourier space as N ik N N j uˆ( k, t) F( u) u( j, t) e, k N J N ik N N j vˆ( k, t) F( v) v( j, t) e, k. N J from this, using the inversion formula, we get ( N, ) ( ˆ) ˆ j (, ) ik j, kn u t F u u k t e j N ( N, ) ( ˆ) ˆ j (, ) ik j, kn v t F v v k t e j N replacing F and F in (9) by their discrete counterparts, and discretizing (9) gives du( j, t) v( j, t) dt dv( j, t) F { k F( u)} u( j, t) F { k F( u)} dt L L () F { ikf( u)} F { k F( u)} L L (9) Letting U = [ u(, t), u(, t),, u(, )] T N t, V = [ v(, t), v(, t),, v(, )] T N t, U Ut w=, w t = V. V t The system of equations () can be written in the vector form w t = F (U, V) () where F defines the right hand side of ().
25 Hany N. HASSAN et al.. Combination of finite differences and a Fourier pseudospectral method The numerical scheme used is based on a combination of finite differences and a Fourier pseudospectral method. After we have changed the solution interval from [a, z] to [,π] and u(, t) is transformed into Fourier space with respect to, and derivatives (or other operators) with respect to, equation (5) become utt F { k F( u)} uf { k F( u)} L L F { ikf( u)} F { k F( u)}. L L In practice, we need to discretize the equation (). For any integer N >, Let u (, t) be the solution equation (). Then, we transform it into the Discrete Fourier space using the inversion formula, we get d u( j, t) { ( )} (, ) { ( )} F k F u u j t F k F u dt L L () F { ikf( u)} F { k F( u)}. L L T [ u(, t), u(, t),, u( N, t)]. Letting U = The equation () can be written in the vector form U tt = F (U) () where F defines the right hand side of (). The time derivative in equation () is discretized using a finite difference approimation, in terms of central differences n n n u u u F (U) ( t) or u(, t t) u(, t) u(, t t) ( t) F { k F( u)} L L u( j, t) F { k F( u)} F { ikf( u)} F { k F( u)}. L L we need two level initial conditions u (,-Δt),u (,). The initial condition u (, ) given in (7), from the central difference of ut (, t ) () (5)
26 Hany N. HASSAN et al. 9 n n u u ut (, t) t where Δt is time step, we have the approimation n n u u tut (, t) we can get u (,-Δt) u(, t) u(, t) tu (,) t u(, t) t b c sec h ( b( )) tanh( b( )). then we substitute u (,-Δt) and u (,) in () to get u (, Δt) ( t) u(, t) u(,) t ut (,) F { k F( u)} L L () u( j,) F { k F( u)} F { ikf( u)} F { k F( u)}. L L So we substitute u (,) and u (, Δt) in () to get u (, Δt) and so on until we get u(, t) at time t. Various values of N ( to ) and time step Δt =. to... Numerical results and eamples In order to show how good the numerical solutions are in comparison with the eact ones, we will use L and L error norms. L u u u u N eact num eact num i i i eact num eact num L u u ma u u. To implement the performance of the method, test problems will be considered: the motion of a single solitary wave in right direction, the motion of a single solitary wave in left direction. Several tests have been made for the wave solution of the Boussinesq verifying that for various values of N ( to ) and time step t =. to. using Combination of finite differences and a Fourier pseudospectral method. As b increases, the stability of the wave propagating in time breaks down and for values of b close to.5 the wave blows up. For stability, we have found that the maimum value of b is.. We take the Boussinesq equation of form (5) with periodic boundary condition u( a, t) u( z, t) i i i
27 u(,t) Hany N. HASSAN et al. 7 Eample We solve the equation (5) with the initial conditions (7),with b =.,, c proposed method performs the motion of propagation of a solitary wave satisfactorily, which moved to the right with the preserved amplitude. b and Δ =, with Δt =..It is clear from Figure that the t L L Table : Error norms of Fourier spectral solution of the Boussinesq equation using RK scheme at t =,, and with N =. N L L Amplitude Table : Error norms of Fourier spectral solution of the Boussinesq equation using RK scheme with N =, 5 and t Fig. Numerical simulation of Fourier spectral solution of the Boussinesq equation using RK scheme with c b and N =.
28 Error u(,t) t Hany N. HASSAN et al Fig.. Plan view of Fourier spectral solution of the Boussinesq equation using RK scheme with c b and N =...7. Numerical Eact Fig.. Fourier spectral solution of the Boussinesq equation using RK scheme at t = with c b and N = Fig.. The error distributions in Fourier spectral solution of the Boussinesq equation using RK scheme at t = with c b and N =.
29 u(,t) u(,t) Hany N. HASSAN et al t = t = t = t = t = Fig.. Fourier spectral solution of the Boussinesq equation using a finite difference scheme at t =,,, and with c b and N =...7 Numerical Eact Fig.5. Fourier spectral solution of the Boussinesq equation using a finite difference scheme at t = with c b and N =.
30 Error Hany N. HASSAN et al Fig..The error distributions of Fourier spectral solution of the Boussinesq equation using a finite difference scheme at t = with c b and N =. t L L Table 5. Error norms of Fourier spectral solution of the Boussinesq equation using a finite difference scheme at t =,,, and with N =. N L L Amplitude Table. Error norms for the single soliton of Fourier spectral solution of the Boussinesq equation using a finite difference scheme with N =,5, 5 and. Now, we solve the same problem with b =., =, c b and Δ = with Δt =.. It is clear from Figure 7 that the proposed method performs the motion of propagation of a solitary wave satisfactorily, which moved to the left with the preserved amplitude.
31 t Hany N. HASSAN et al. 7 u(,t)..5 t Fig. 7. Numerical simulation of Fourier spectral solution of the Boussinesq equation using RK scheme with c b and N = Fig.. Plan view of Fourier spectral solution of the Boussinesq equation using RK scheme with c b and N =.
32 Error u(,t) Hany N. HASSAN et al Numerical Eact Fig.9. Fourier spectral solution of the Boussinesq equation using RK scheme at t = with c b and N = Fig.. The error distributions of Fourier spectral solution of the Boussinesq equation using RK scheme at t = with c b and N =.
33 u(,t) u(,t) Hany N. HASSAN et al t = t = t = t = t = Fig.. Fourier spectral solution of the Boussinesq equation using a finite difference scheme at c b, t =,,, and with N =...7 Numerical Eact Fig.. Fourier spectral solution of the Boussinesq equation using a finite difference scheme at c b, t = with N =.
34 Error Hany N. HASSAN et al Fig.. The error distributions in Fourier spectral solution of the Boussinesq equation using a finite difference scheme with c b at t = with N =. Conclusions In Our study, we applied Fourier spectral collocation methods to solve partial differential equations of the form ut F ( u) particularly the Korteweg-de Vries equation (KdV), and the form utt F ( u) particularly the Boussinesq equation. A finite differences or a fourth-order Runge-Kutta method is used to march our solution in time, with the help of the fast Fourier transform algorithm our methods cost only (number of operations) O N log N, provided that N is a power of. We applied the leap-frog scheme combined with the Fourier spectral collocation, called the Fourier Leap-Frog method, to find numerical solution of the KdV equation, with α =, and β =. Also we applied the fourth-order Runge-Kutta (RK) method combined with the Fourier spectral collocation to the same problem. We presented stability conditions for these schemes, and we found out that the Fourier Leap-Frog method is less stable and far less accurate compared to the classical RK scheme. Even though it is quite costly to use, RK is easy to implement and needs only one level of initial data, whereas two levels of the initial data are needed in the Fourier Leap-Frog methods. We then applied the finite difference scheme combined with the Fourier spectral collocation to find numerical solution of the Bossinesq equation. Also we applied the fourth-order Runge-Kutta (RK) method combined with the Fourier spectral collocation to the same problem. In order to show how good the numerical solutions are in comparison with the eact ones, we will use L and L error
35 Hany N. HASSAN et al. 7 norms. It is apparently seen that Fourier spectral collocation method is powerful and efficient technique in finding numerical solutions for wide classes of nonlinear partial differential equations. 5 Open Problem The work presented in this paper transform some types of partial differential equations like the Korteweg-de Vries equation (KdV), Boussinesq equation to simple ordinary differential equations using Fourier spectral method, can be solved by simple techniques ( Leap frog, finite difference, Runge Kutta, ), the question here what the results of using other spectral methods to solve these equations instead of Fourier spectral method. The other question, in this paper we applied the used method to solve solitons, what happened when we solve equations by non solitons spectral methods. References [] Fornberg B. A Practical Guide to Pseudospectral Methods. Cambridge University Press, New York. 99. [] Fornberg B and Whitham G B. A numerical and Theoretical Study of Certain Nonlinear Wave Phenomena. Philos. Trans. Soc. London, 9, May 97, pp. 7-. [] Hany N. Hassan and Hassan K. Saleh, The Solution of the Regularized Long Wave Equation Using the Fourier Leap-Frog Method, Z. Naturforsch. 5a, 7 () [] Boyd J P. Chebyshev and Fourier Spectral Methods. Second Edition Dover Publications Inc., New York,. [5] Whitham G B. Linear and Nonlinear Waves. John Wiley &Sons, 97. [] Weideman, J.A.C., and Reddy, S.C., A MATLAB Differentiation Matri Suite, ACM Transactions on Mathematical Software,, No., December, pp [7] Efstratios, E., Tzirtzilakis, Charalampos D., Skokos, and Tassos C. Bountis, Numerical Solution of the Bossinesq equation using Spectral Methods and Stability of Solitary Wave Propagation, first International Conference From Scientific Computing to Computational Engineering, Athens, - September,. [] Eric Gourgoulhon, Introduction to spectral methods, th EU Network Meeting, Palma de Mallorca, Sept.. [9] Tzirtzilakis, E., Marinakis, V., Apokis, C., and Bountis, T., Soliton-like solutions of higher order wave equations of the Korteweg de Vries type, Journal of Mathematical Physics, No., December, pp. 5 5.
36 Hany N. HASSAN et al. 79 [] Bengt Fornberg and Tobin A. Driscoll, A Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion, Journal of Computational Physics55, 999, pp.5 7. [] Randall, J. LeVeque, Finite Difference Methods for Differential Equations, A Math 55, Winter Quarter, University of Washington, Version of January,. [] Drazin, P.G. and Johnson, R.S., Solitons: an introduction, Cambridge University Press, Cambridge, 99. [] Idris Dag, M. Naci Ozer, Approimation of the RLW equation by the least square cubic B-spline finite element method Applied Mathematical Modelling 5 () - [] H. Jafari, H. Hosseinzadeh and S. Mohamadzadeh Numerical Solution of System of Linear Integral Equations by using Legendre Wavelets Int. J. Open Problems Compt. Math., Vol., No. 5, December [5] Hany N. Hassan and M. A. El-Tawil, A new technique of using homotopy analysis method for solving high-order nonlinear differential equations, Mathematical Methods in the Applied Sciences,, 7 7 ().
Galerkin 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 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 informationLong Time Dynamics of Forced Oscillations of the Korteweg-de Vries Equation Using Homotopy Perturbation Method
Studies in Nonlinear Sciences 1 (3): 57-6, 1 ISSN 1-391 IDOSI Publications, 1 Long Time Dynamics of Forced Oscillations of the Korteweg-de Vries Equation Using Homotopy Perturbation Method 1 Rahmat Ali
More informationBBM equation with non-constant coefficients
Turkish Journal of Mathematics http://journals.tubitak.gov.tr/math/ Research Article Turk J Math (3) 37: 65 664 c TÜBİTAK doi:.396/mat-3-35 BBM equation with non-constant coefficients Amutha SENTHILKUMAR
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 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 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 information(Received 05 August 2013, accepted 15 July 2014)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.18(2014) No.1,pp.71-77 Spectral Collocation Method for the Numerical Solution of the Gardner and Huxley Equations
More informationResearch Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal Polynomial
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 4673, pages doi:.54//4673 Research Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal
More informationNUMERICAL SOLUTIONS OF THE FORCED PERTURBED KORTEWEG-DE VRIES EQUATION WITH VARIABLE COEFFICIENTS
International Journal of Pure and Applied Mathematics Volume 112 No. 3 217, 557-569 ISSN: 1311-88 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 1.12732/ijpam.v112i3.8
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
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 informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
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 informationApplication of Laplace Adomian Decomposition Method for the soliton solutions of Boussinesq-Burger equations
Int. J. Adv. Appl. Math. and Mech. 3( (05 50 58 (ISSN: 347-59 IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Application of Laplace Adomian
More informationA robust uniform B-spline collocation method for solving the generalized PHI-four equation
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 11, Issue 1 (June 2016), pp. 364-376 Applications and Applied Mathematics: An International Journal (AAM) A robust uniform B-spline
More informationMaejo International Journal of Science and Technology
Full Paper Maejo International Journal of Science and Technology ISSN 905-7873 Available online at www.mijst.mju.ac.th New eact travelling wave solutions of generalised sinh- ordon and ( + )-dimensional
More informationSimulation for Different Order Solitons in Optical Fibers and the Behaviors of Kink and Antikink Solitons
Simulation for Different Order Solitons in Optical Fibers and the Behaviors of Kink and Antikink Solitons MOHAMMAD MEHDI KARKHANEHCHI and MOHSEN OLIAEE Department of Electronics, Faculty of Engineering
More informationExact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable Coefficients
Contemporary Engineering Sciences, Vol. 11, 2018, no. 16, 779-784 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.8262 Exact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable
More informationNew Approach of ( Ǵ/G ) Expansion Method. Applications to KdV Equation
Journal of Mathematics Research; Vol. 6, No. ; ISSN 96-9795 E-ISSN 96-989 Published by Canadian Center of Science and Education New Approach of Ǵ/G Expansion Method. Applications to KdV Equation Mohammad
More informationGroup analysis, nonlinear self-adjointness, conservation laws, and soliton solutions for the mkdv systems
ISSN 139-5113 Nonlinear Analysis: Modelling Control, 017, Vol., No. 3, 334 346 https://doi.org/10.15388/na.017.3.4 Group analysis, nonlinear self-adjointness, conservation laws, soliton solutions for the
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 informationNew approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations
Volume 28, N. 1, pp. 1 14, 2009 Copyright 2009 SBMAC ISSN 0101-8205 www.scielo.br/cam New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations HASSAN A. ZEDAN Mathematics
More informationRestrictive Taylor Approximation for Gardner and KdV Equations
Int. J. Adv. Appl. Math. and Mech. 1() (014) 1-10 ISSN: 47-59 Available online at www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Restrictive Taylor Approximation
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 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 informationDavydov Soliton Collisions
Davydov Soliton Collisions Benkui Tan Department of Geophysics Peking University Beijing 100871 People s Republic of China Telephone: 86-10-62755041 email:dqgchw@ibmstone.pku.edu.cn John P. Boyd Dept.
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 informationComputational Solutions for the Korteweg devries Equation in Warm Plasma
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 16(1, 13-18 (1 Computational Solutions for the Korteweg devries Equation in Warm Plasma E.K. El-Shewy*, H.G. Abdelwahed, H.M. Abd-El-Hamid. Theoretical Physics
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 informationSTUDIES IN RECURRENCE IN NONLINEAR DISPERSIVE WAVE EQUATIONS
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
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 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 informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationStrongly nonlinear long gravity waves in uniform shear flows
Strongly nonlinear long gravity waves in uniform shear flows Wooyoung Choi Department of Naval Architecture and Marine Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA Received 14 January
More informationA cubic B-spline Galerkin approach for the numerical simulation of the GEW equation
STATISTICS, OPTIMIZATION AND INFORMATION COMPUTING Stat., Optim. Inf. Comput., Vol. 4, March 216, pp 3 41. Published online in International Academic Press (www.iapress.org) A cubic B-spline Galerkin approach
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 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 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 informationSolving Poisson s Equations Using Buffered Fourier Spectral Method
Solving Poisson s Equations Using Buffered Fourier Spectral Method Yinlin Dong Hassan Abd Salman Al-Dujaly Chaoqun Liu Technical Report 2015-12 http://www.uta.edu/math/preprint/ Solving Poisson s Equations
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 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 informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationA Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs
Journal of Computational Physics 8, 57 67 (00) doi:0.006/jcph.00.77 A Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs Tobin A. Driscoll Department of Mathematical Sciences, University
More informationNumerical study of some nonlinear wave equations via Chebyshev collocation method
IJST (203) 37A4: 477-482 Iranian Journal of Science & Technology http://ijsts.shirazu.ac.ir Numerical study of some nonlinear wave equations via Chebyshev collocation method N. Y. A. Elazem and A. Ebaid*
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 informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationResearch Article A New Method for Riccati Differential Equations Based on Reproducing Kernel and Quasilinearization Methods
Abstract and Applied Analysis Volume 0, Article ID 603748, 8 pages doi:0.55/0/603748 Research Article A New Method for Riccati Differential Equations Based on Reproducing Kernel and Quasilinearization
More informationMaking Waves in Vector Calculus
Making Waves in Vector Calculus J. B. Thoo Yuba College 2014 MAA MathFest, Portland, OR This presentation was produced
More informationNew Solutions for Some Important Partial Differential Equations
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.2,pp.109-117 New Solutions for Some Important Partial Differential Equations Ahmed Hassan Ahmed Ali
More informationSolution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.3,pp.227-234 Solution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition
More informationA Good Numerical Method for the Solution of Generalized Regularized Long Wave Equation
Modern Applied Science; Vol. 11, No. 6; 2017 ISSN 1913-1844 E-ISSN 1913-1852 Published by Canadian Center of Science and Education A Good Numerical Method for the Solution of Generalized Regularized Long
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 informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationPeriodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type Equation
Contemporary Engineering Sciences Vol. 11 2018 no. 16 785-791 HIKARI Ltd www.m-hikari.com https://doi.org/10.12988/ces.2018.8267 Periodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationExact travelling wave solutions of a variety of Boussinesq-like equations
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 2 Exact travelling wave solutions of a variety
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More informationPeriodic and Solitary Wave Solutions of the Davey-Stewartson Equation
Applied Mathematics & Information Sciences 4(2) (2010), 253 260 An International Journal c 2010 Dixie W Publishing Corporation, U. S. A. Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
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 informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationNew explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation
Physics Letters A 07 (00) 107 11 www.elsevier.com/locate/pla New explicit solitary wave solutions for ( + 1)-dimensional Boussinesq equation and ( + 1)-dimensional KP equation Yong Chen, Zhenya Yan, Honging
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 informationChapter 3. Head-on collision of ion acoustic solitary waves in electron-positron-ion plasma with superthermal electrons and positrons.
Chapter 3 Head-on collision of ion acoustic solitary waves in electron-positron-ion plasma with superthermal electrons and positrons. 73 3.1 Introduction The study of linear and nonlinear wave propagation
More informationHere we consider soliton solutions of the Korteweg-de Vries (KdV) equation. This equation is given by
17 Solitons Lab Objective: We study traveling wave solutions of the Korteweg-de Vries (KdV) equation, using a pseudospectral discretization in space and a Runge-Kutta integration scheme in time. Here we
More information2007 Summer College on Plasma Physics
856-7 Summer College on Plasma Physics July - 4 August, 7 Numerical methods and simulations. Lecture : Simulation of partial differential equations. B. Eliasson Institut Fuer Theoretische Physik IV Ruhr-Universitaet
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 informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
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 informationPhase Plane Behavior of Solitary Waves in Nonlinear Layered Media
Phase Plane Behavior of Solitary Waves in Nonlinear Layered Media Randall J. LeVeque 1 and Darryl H. Yong 2 1 Department of Applied Mathematics, University of Washington, Box 35242, Seattle, WA 98195-242
More informationPseudo-Spectral Methods for Linear Advection and Dispersive Problems
Pseudo-Spectral Methods for Linear Advection and Dispersive Problems Riccardo Fazio and Salvatore Iacono Abstract The most important feature of numerical methods based on a spectral decomposition is the
More informationExact Solutions of the Generalized- Zakharov (GZ) Equation by the Infinite Series Method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 05, Issue (December 010), pp. 61 68 (Previously, Vol. 05, Issue 10, pp. 1718 175) Applications and Applied Mathematics: An International
More informationHere we consider soliton solutions of the Korteweg-de Vries (KdV) equation. This equation is given by u t + u u x + 3 u
Lab 3 Solitons Lab Objective: We study traveling wave solutions of the Korteweg-de Vries (KdV) equation, using a pseudospectral discretization in space and a Runge-Kutta integration scheme in time. Here
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 informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 7 (0) 3499 3507 Contents lists available at SciVerse ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Application of the
More informationSoliton and Periodic Solutions to the Generalized Hirota-Satsuma Coupled System Using Trigonometric and Hyperbolic Function Methods.
ISSN 1749-889 (print), 1749-897 (online) International Journal of Nonlinear Science Vol.14(01) No.,pp.150-159 Soliton and Periodic Solutions to the Generalized Hirota-Satsuma Coupled System Using Trigonometric
More informationCHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION
CHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION A version of this chapter has been published as Haar Wavelet Method for solving Fisher s equation, Appl. Math.Comput.,(ELSEVIER) 211 (2009)
More informationNUMERICAL INVESTIGATION OF THE DECAY RATE OF SOLUTIONS TO MODELS FOR WATER WAVES WITH NONLOCAL VISCOSITY
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 1, Number 2, Pages 333 349 c 213 Institute for Scientific Computing and Information NUMERICAL INVESTIGATION OF THE DECAY RATE OF SOLUTIONS
More informationBenha University Faculty of Science Department of Mathematics. (Curriculum Vitae)
Benha University Faculty of Science Department of Mathematics (Curriculum Vitae) (1) General *Name : Mohamed Meabed Bayuomi Khader *Date of Birth : 24 May 1973 *Marital Status: Married *Nationality : Egyptian
More informationNumerical Solutions of the Combined KdV-MKdV Equation by a Quintic B-spline Collocation Method
Appl. Math. Inf. Sci. Lett. 4 No. 1 19-24 (2016) 19 Applied Mathematics & Information Sciences Letters An International Journal http://d.doi.org/18576/amisl/040104 Numerical Solutions of the Combined KdV-MKdV
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationSoliton solutions of Hirota equation and Hirota-Maccari system
NTMSCI 4, No. 3, 231-238 (2016) 231 New Trends in Mathematical Sciences http://dx.doi.org/10.20852/ntmsci.2016115853 Soliton solutions of Hirota equation and Hirota-Maccari system M. M. El-Borai 1, H.
More informationThe (G'/G) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics
Vol.3, Issue., Jan-Feb. 3 pp-369-376 ISSN: 49-6645 The ('/) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics J.F.Alzaidy Mathematics Department, Faculty
More informationAnalysis of the equal width wave equation with the mesh-free reproducing kernel particle Ritz method
Analysis of the equal width wave equation with the mesh-free reproducing kernel particle Ritz method Cheng Rong-Jun( 程荣军 ) a) and Ge Hong-Xia( 葛红霞 ) b) a) Ningbo Institute of Technology, Zhejiang University,
More informationSolving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman
Solving Nonlinear Wave Equations and Lattices with Mathematica Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs home/whereman/
More informationA solitary-wave solution to a perturbed KdV equation
University of Warwick institutional repository: http://go.warwick.ac.uk/wrap This paper is made available online in accordance with publisher policies. Please scroll down to view the document itself. Please
More informationA numerical solution of a Kawahara equation by using Multiquadric radial basis function
Mathematics Scientific Journal Vol. 9, No. 1, (013), 115-15 A numerical solution of a Kawahara equation by using Multiquadric radial basis function M. Zarebnia a, M. Takhti a a Department of Mathematics,
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 informationParallel Numerical Methods for Solving Nonlinear Evolution Equations
Parallel Numerical Methods for Solving Nonlinear Evolution Equations Thiab R. Taha Department of Computer Science University of Georgia Athens, GA 362 email: thiab@cs.uga.edu June 1-14, 21 Abstract Nonlinear
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationAbsorbing Boundary Conditions for Nonlinear Wave Equations
Absorbing Boundary Conditions for Nonlinear Wave Equations R. B. Gibson December, 2010 1 1 Linear Wave Equations 1.1 Boundary Condition 1 To begin, we wish to solve where with the boundary condition u
More informationProgramming of the Generalized Nonlinear Paraxial Equation for the Formation of Solitons with Mathematica
American Journal of Applied Sciences (): -6, 4 ISSN 546-99 Science Publications, 4 Programming of the Generalized Nonlinear Paraxial Equation for the Formation of Solitons with Mathematica Frederick Osman
More informationPropagation of Solitons Under Colored Noise
Propagation of Solitons Under Colored Noise Dr. Russell Herman Departments of Mathematics & Statistics, Physics & Physical Oceanography UNC Wilmington, Wilmington, NC January 6, 2009 Outline of Talk 1
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationNEW PERIODIC WAVE SOLUTIONS OF (3+1)-DIMENSIONAL SOLITON EQUATION
Liu, J., et al.: New Periodic Wave Solutions of (+)-Dimensional Soliton Equation THERMAL SCIENCE: Year 7, Vol., Suppl., pp. S69-S76 S69 NEW PERIODIC WAVE SOLUTIONS OF (+)-DIMENSIONAL SOLITON EQUATION by
More informationOn weakly singular and fully nonlinear travelling shallow capillary gravity waves in the critical regime
arxiv:1611.115v3 [physics.class-ph] 8 Mar 17 Dimitrios Mitsotakis Victoria University of Wellington, New Zealand Denys Dutykh CNRS, Université Savoie Mont Blanc, France Aydar Assylbekuly Khoja Akhmet Yassawi
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 informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More information