In this work we are interested in developing approximate models for water waves equation. We present the derivation of the new equations uses approximation of the phase velocity that arises in the linear water wave theory. We treat the (KdV) equation and similarly the C-H equation. Both of them describe unidirectional shallow water waves equation.

At the same time, together with the (BBM) equation we propose, we provide the best approximation of the phase velocity for small wave numbers that can be obtained with second and third-order equations. We can extend the results of [3, 4]. A comparison between the methods is mentioned in this article.

In many cases the numerical methods gives us a good approximation related to approximation of solutions. The finite-element method is a special case of the Galerkin method in which the base functions are chosen such that each base function becomes 1 at the corresponding nodes but otherwise 0 at other nodes. The Galerkin method permitted finite-element techniques to be extended into areas such as fluid mechanics, especially in water waves theory. Milder, Miles and Broer deal with the water waves on the Hamiltonian theory of surface waves (see

The KdV equation was originally derived by Korteweg and de Vries

We are concerned with a two-dimensional, irrotational flow of an incompressible ideal fluid with a free surface under the gravitational field. The domain occupied by the fluid is bounded from below by a solid bottom. The upper surface is a free boundary and we take the influence of the gravitational field into account on the free surface. Our main interest is motion of the free surface, which is called a gravity wave.

The KdV equation was obtained by Benjamin

The nonlinear surface water waves in perturbed problem are discussed by

The main aim in this work is to effectively derive the KdV equations and employ the modified form of Laplace decomposition method introduced by Khuri

The Galerkin method permitted finite-element techniques to be extended into areas such as fluid mechanics and series solution of some problems in elastic equilibrium of rods and plates (see

The equation for the surface of water is y=η(x,t) where y=0 represents the equilibrium situation. We assume that the motion is irrotational within the wave. Therefore, we can describe the wave inside the water by a velocity potential Φ(x,y,t) whose gradient is the velocity field

(∇Φ = ∂Φ/∂x,∂Φ/∂y ) ( 2.1)

The divergence-free condition on the velocity field implies that the velocity potential Φ satisfies the Laplace's equation:

(∂^2 Φ)/〖∂x〗^2 +(∂^2 Φ)/〖∂y〗^2 =0,for-h<y<η(x,t) (2.2)

We consider the unsteady two-dimensional flow of inviscid, incompressible fluid in a constant gravitational field. The space coordinates are (𝑥,y) and the gravitational acceleration 𝑔 is in the negative 𝑦 direction. Let ℎ be the undisturbed depth of the fluid. The bottom of the fluid is assumed to have no topography at 𝑦 = − ℎ. This problem describes the interface dynamics, between air and water waves, under the gravity 𝑔.

Since the majority of shallow water ﬂow of practical interest is two-dimensional, the governing differential equation can be written as:

(∂^2 (Uu))/(∂x^2 )+(∂^2 (Uv))/∂x∂y+(∂^2 (Uτ))/(∂y^2 )=(p_z (x,y))/S+(∂(Uh))/(∂x^2 ) (2.3)

where, is the lateral pressure that is being applied, is the flexural deformation of the wave.

Here, and denote the velocity components in and directions, respectively, and denotes the depth of the free surface.

On a solid fixed boundary, the normal velocity of the fluid must vanish aty= -h

∂Φ/∂y =0 at y=-h, (2.3)

which dictates that there is no flow perpendicular to the bottom. The boundary conditions at the free surface y=η(x,t) are given by

∂η/∂t+∂η/∂x ∂Φ/∂x= ∂Φ/∂y (2.4)

The Camassa-Holm (CH) equation:

〖 (u-u_xx)〗_t+3uu_(x )=ε (〖2 u_(x ) u〗_(xx )+〖uu〗_(xxx )) , xϵ Ω,t>0 (2.5)

or u_(t )-u_txx+3uu_(x )- 〖2 u_(x ) u〗_(xx )-〖uu〗_(xxx )=0 (2.6)

In this section, we show that it is indeed possible to derive KdV equations

ϵ (∂^2 Φ)/〖∂x〗^2 +(∂^2 Φ)/〖∂y〗^2 =0 , for -1<y<ϵ η (x,t) (3.1.1)

∂Φ/∂y =0 at y=-1, (3..1.2)

ϵ ∂η/∂t+ϵ^2 ∂η/∂x ∂Φ/∂x= ∂Φ/∂y at y=ϵ η (3.1.3)

To get a suitable expansion for the velocity potential Φ, we note that from the Laplace equation (3.1) and the boundary condition (3.2), we may derive the following representation for Φ(𝑥,𝑦,𝑡) as a series in 𝑦+1. Since 𝜕Φ/𝜕𝑦 = 0 at 𝑦 = −1,

Ф(x,y,t)=∑_(n=0)^∞▒〖(y+1)^n Ф_n (x,t) (3.1.4)

Recall that the conditions on a simply supported edge parallel to the y-axis at

U|_(X=a)=0 (3.2.1)

〖 M〗_x |_(X=a)=-S((∂^2 w)/(∂x^2 )+ϑ (∂^2 w)/(∂y^2 ) )_(X=a)=0 (3.2.2)

On a simply supported edge parallel to the x-axis at

U|_(Y=a)=0 (3.2.3)

〖 M〗_y |_(Y=a)=-S((∂^2 w)/(∂y^2 )+ϑ (∂^2 w)/(∂x^2 ) )_(Y=a)=0 (3.2.4)

The following part presents the numerical results achieved for the problems of approximation. All of the computation was carried out with the help of a symbolic algebra software system

In this section, a brief outline of LDM is explained. For this, we consider the general nonlinear partial differential equation of first order (without loss of generality) in the following form:

Lu(x,t)+Ru(x,t)+N(u(x,t))=h(x,t) (4.1)

with the following initial condition:

u(x,0)=f(x) (4.2)

where 𝐿 is the first-order differential operator, 𝐿=𝜕/𝜕𝑡, 𝑅 is linear differential operator, 𝑁(𝑢) presents the nonlinear term, and ℎ(𝑥,𝑡) is the source term. The methodology consists of applying Laplace transform first on both sides of (4.1)

£[Lu(x,t)]+£[Ru(x,t)]+£[N(u(x,t))]=£[h(x,t)] (4.3)

Using the differentiation property of Laplace transform, we get

s£[u(x,t)]-f(x)+£[Ru(x,t)]+£[N(u(x,t))]=£[h(x,t)] (4.4)

Now, we will define the solution (𝑥,) by the series in the following form:

u(x,t)=∑_(n=0)^∞ u_n (x,t) (4.5)

and the nonlinear operator 𝑁(𝑢) represented by an infinite series of the so-called Adomian's polynomials:

N(u)=∑_(n=0)^∞▒〖A_n〗 (4.6)

where u_n (x,t), n≥0 are the components of 𝑢(𝑥,𝑡) that will be elegantly determined and A_n are called Adomian's polynomials and defined by

〖 A〗_n=1/n! [d^n/(dλ^n ) N(∑_(i=0)^∞▒〖λ^i u_i 〗)]_(λ=0), n≥0 (4.7)

Using (4.5) and (4.6) in (4.4), we get

£⌊∑_(n=0)^∞▒〖u_n (x,t)〗⌋=1/s f(x) 1/s £[h(x,t)]-1/s £⌊R∑_(n=0)^∞▒〖u_n (x,t)〗⌋+1/s £⌊∑_(n=0)^∞▒A_n ⌋ (4.8)

For numerical comparisons purpose, we construct the solution (𝑥,) such that

〖 lim┬(n→∞)〗〖U_n (x,t)=u(x,t)〗, where 〖 U〗_n (x,t)=∑_(i=0)^(n-1)▒〖u_i (x,t), n≥0〗(4.9)

The method has been used to solve problems in mechanical engineering such as structural mechanics, dynamics, fluid flow, heat and mass transfer, acoustics and other related fields. We consider the idea to approximate the solution with a polynomial involving a set of parameters. The polynomial is made to satisfy both the differential equation and the associated boundary conditions. Using orthogonality of polynomials, we can approximate the solution to the differential equation on shallow water wave equation.

Let´s we have the second differential equation for a horizontal surface for water waves as follow

(∂^2 U)/(∂x^2 )+(∂^2 U)/∂x∂y+(∂^2 U)/(∂y^2 )=(p_z (x,y))/S+∂U/(∂x^2 ) (5.1)

Consider a differential operator L, defined as

L=[∂^2/(∂x^2 )+∂^2/∂x∂y+∂^2/(∂y^2 )] (5.2)

So, eq. (2.1) becomes, LU=P_z/S (5.3)

We can express the solution to the above equation in terms of the eigenfunction and eigenvalues, which are defined as:

Le_nm (x,y)=λ_nm Le_nm (x,y) (5.4)

where e_nm (x,y) are the eigenfunctions and λ_nm are the corresponding eigenvalues. Once the eigenfunctions and eigenvalues are known, it is possible to express U(x,y) as,

U(x,y)=∑_(n=1)^∞▒∑_(m=1)^∞▒〖u_nm e^(λ_nm ) e_nm (x,y)〗 (5.5)

Polynomial approximating functions will be used to represent the lateral deformation of the horizontal zone water. In this discussion, the trial function (

∅_i (x,y)=∑_(i=1)^N▒∑_(j=1)^N▒〖a[i,j] u_j (x,y)〗 (5.6)

where

u_j (x,y)=x^(L_j ) 〖∙y〗^(M_j ) (5.7)

and, L_j and M_j are positive integers and ∅_i (x,y) are coefficients to be determined.

We see that the method mentioned in this article is one of the methods which were used to obtain the KdV equations. It is relatively short in mathematical calculation and more effective.

The Galerkin method can be used to approximate the solution to ordinary differential equations, partial differential equations and integral equations. The Galerkin method is used to determine the coefficients of these polynomials. The finite element and Galerkin methods are currently the standard numerical technique in use to solve various nonlinear problems.

In this study, we present model equations for surface water waves by using Laplace decomposition technique. It is used to estimate the KdV equations for the nonlinear theory, describing the behaviour of the perturbed system.

While the Hamiltonian expansions are complex in mathematical calculation. For the same, the Dirichlet-Neumann operator expansions is relatively a long method when compared with the Galerkin method. Special attention is given to derive the solutions for the KdV equations, which describe the water waves propagation by using the MLDM, and then analyzed and discussed theoretically and computationally.

We observe that the elevation of the water waves is in form of traveling solitary waves, which increases in amplitude as the wave number increases.

Finally, the horizontal and vertical velocities have a nonlinear characters, which describe the physical situation of the system for free surface between air and water.

1. LDM does not require small parameters which are needed in perturbation method;

2. the method is very promising for solving wide application in nonlinear differential equations.

The extension of the general finite element method of structural analysis to water waves studies, since almost any particular situation can now be analyzed using a general purpose computer program.

In the last years, with the advent of the digital computer, other methods such as finite differences and finite elements have become practical. In 1956 Clough, Martin and Topp introduced the finite element method.