The Numerical Solution of the MRLW Equation Using the Multigrid Method

In this paper, we obtained the numerical solutions of the modified regularized long-wave (MRLW) equation t x x xxt u u u u μu 2 0 + + − = α , by using the multigrid method and finite difference method. The solitary wave motion, interaction of two and three solitary waves, and development of the Maxwellian initial condition into solitary waves are studied using the proposed method. The numerical solutions are compared with the known analytical solutions. Using L L 2 , ∞ error norms and conservative properties of mass, momentum and energy, accuracy and efficiency of the mentioned method will be established through comparison with other techniques.


Introduction
The numerical solution of partial differential equations requires some discretization of the domain into a collection of points.A large system of equations comes out from discretization of the same partial differential equations and the optimal method for solving these problems is multigrid method, see [1]- [9].
Consider the following one-dimensional modified regularized long-wave (MRLW) equation: equation: 2 0, where t is the time, x is the space coordinate, , α µ are positive constants and u is the wave amplitude with the physical boundary conditions 0 u → as x → ±∞ .This equation was first introduced to describe the development of an undular bore by Peregrine [10] and later by Benjamin et al. [11].Equation (1) has various applications as in physics media since it describes the phenomena with weak nonlinearity and dispersion waves, including nonlinear transverse waves in shallow water, ion-acoustic and magneto hydrodynamic waves in plasma, and phonon packets in nonlinear crystals [11].
Although the analytical solutions of the MRLW equation, with a limited set of boundary and initial conditions, have been existed, many authors are recently interested in the numerical solutions of this equation.Gardner et al. [12] introduced a collocation solution to the MRLW equation using quintic B-spline finite elements.Khalifa et al. [13] [14] applied the finite difference and cubic B-spline collocation finite element method to obtain the numerical solutions of the MRLW equation.Solutions based on collocation method with quadratic B-spline finite elements and the central finite difference method for time are investigated by Raslan [15].Raslan and Hassan [16] solved the MRLW equation by a collocation finite element method using quadratic, cubic, quartic, and quintic B-spline to obtain the numerical solutions of the single solitary wave.Ali [17] has formulated a classical radial basis function collocation method for solving the MRLW equation.Haq et al. [18] have developed a numerical scheme based on quartic B-spline collocation method for the numerical solution of MRLW equation.Karakoc and Geyikli [19] solved the MRLW equation by using the Petrov-Galerkin finite element method.
An outline of this paper is as follows: we begin in Section 2 by reviewing the analytical solution of the MRLW equation.In Section 3, we derive a new numerical method based on the multigrid technique and finite difference method for obtaining the numerical solution of MRLW equation.Finally, in Section 4, we introduce the numerical results for solving the MRLW equation through some well known standard problems.

The Analytical Solution
The exact solution of Equation ( 1) can be written in the form [12] [14]: which represents the motion of a single solitary wave with amplitude c , where ( ) 0 The conservation properties of the MRLW equation related to mass, momentum and energy are determined by following three invariants on the region a x b ≤ ≤ : ( )

Numerical Method
The basic idea of multigrid techniques is illustrated by Brandt [1].In this section we apply this method for initial boundary value problem, except that, the upper boundary conditions change with time, in which the initial condition is ( ) ( ) for 0 t T < < .Dividing the interval of time to K parts, we obtain the solutions of the partial differential equation at time t 1 and use these solutions as initial values for the next level ( ) ( ) , and for the other, we obtain the solutions at time T .The numbers of points in a coarse grid for this domain are two points.We apply the full multigrid algorithm for the MRLW equation.Assuming the initial condition ( ) ( ) ,0 u x f x = and the solution ( ) has the usual partition with a space step size x ∆ and a time step size t ∆ ( ) We start handling the non-linear term 2 x u u by ex- pressing in the form Step 2: Starting from 1 k = in the coarse grid, we can calculate the approximate value , i n u at two points using Equation ( 5) leading to: The right hand side for the last equation can be computed using the initial and boundary conditions.
Step 3: Interpolating the grid functions from the coarse grid to fine grid using linear interpolation that can be written explicitly as: Step 4: Doing relaxation sweep on Step 5: Computing the residuals ( ) Step 6: Computing an approximate solution of error e k .
Step 7: Interpolating the solution of error e k onto Step 8: 1 K K = + , go to step 2 (lead to the solution at higher time level as needed).

Numerical Results
In this section, numerical solutions of MRLW equation are obtained for standard problems as: the motion of Y. M. Abo Essa et al.
single solitary wave, interaction of two and three solitary waves and development of Maxwellian initial condition into solitary waves. 2 L and L ∞ error norms are used to show how good the numerical results in compar- ison with the exact results.

The Motion of Single Solitary Wave
Consider equation ( 1) with boundary conditions x x u a t u b t u a t u b t t and the initial condition (4).
The analytical values of the invariants of this problem can be found as [12]: I at various times we use the numerical solutions by applying the multigrid method up to 10 t = .As reported in Table 1, the error norms 2 L , L ∞ are found to be small enough, and the computed val- ues of invariants are in good agreement with their analytical values 1

I =
Table 2 shows a comparison of the values of the invariants and error norms obtained by the present method with those obtained by other methods [12] [13] [17] [19].It is clearly seen from Table 2 that the error norms obtained by the present method are smaller than the other methods.

Interaction of Two Solitary Waves
Consider the interaction of two separated solitary waves having different amplitudes and travelling in the same direction as a second problem.For this problem, the initial condition is given by:   For the computational discussion, we use parameters I are listed in Table 3. Table 3 shows a comparison of the values of the invariants obtained by present method with those obtained in [19].It is seen that the numerical values of the invariants remain almost constant during the computer run.

Interaction of Three Solitary Waves
In this section, the behavior of the interaction of three solitary waves having different amplitudes and travelling in the same direction was studied.So, we consider Equation (1) with the initial condition given by the linear sum of three well-separated solitary waves of different amplitudes: For the computational work, we used parameters  4. Table 4 shows a comparison of the values of the invariants obtained by the present method with those obtained in [19].It is seen that the numerical values of the invariants remain almost constant during the computer run.

The Maxwellian Initial Condition
Finally, the development of the Maxwellian initial condition: into a train of solitary waves is discussed.It is known that the behavior of the solution with the Maxwellian condition (16) depends on the values of µ .So, we study each of two cases: 0.015 µ = and 0.004.

µ =
Table 5 contains the obtained numerical values of the invariants and a comparison of the values of the invariants obtained by present method with those obtained in [19].

Conclusion
In this work we extended the use of multigrid technique to initial boundary value problems, namely the MRLW problem.We tested our scheme through single solitary wave in which the analytic solution is known.Our scheme was extended to study the interaction of two and three solitary waves and Maxwellian initial condition where the analytic solutions are unknown during the interaction.The performance and accuracy of the method were shown by calculating the error norms 2 L , L ∞ and conservative properties of mass, momentum and energy.computed results showed that the present scheme is a successful numerical technique for solving the MRLW problem.
-time and centre-space difference for Equation (1) is is the approximate value of u on the fine grid with 2 k = .By taking this solution on coarse grid and repeating steps 3 -7, we obtain the approximate values of u on the is the solution at the time level 1 K + .
with earlier studies[12] [13] [17][19] we take the parameters 6the interval [0, 100].To find the error norms 2 L , L ∞ and the numerical inva- riants 1 2 , I I and 3 is run up to time 1 t = and numerical values of the in- variant quantities1 2   , I I and 3 I are displayed in Table

Table 1 .
Invariants and error norms for single solitary wave when

Table 2 .
Comparison of errors and invariants for single solitary wave when

Table 3 .
[19]arison of invariants for the interaction of two solitary waves with results from[19]when 1 0.03,

Table 4 .
[19]arison of invariants for the interaction of three solitary waves with results from[19]when 1 0.03, c =

Table 5 .
Invariants of MRLW equation using the Maxwellian condition.