Startseite New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation
Artikel Open Access

New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation

  • Ahmed Hussein Msmali , Mohammad Tamsir , Neeraj Dhiman EMAIL logo und Mohammed A. Aiyashi
Veröffentlicht/Copyright: 6. Dezember 2021

Abstract

The objective of this work is to propose a collocation technique based on new cubic trigonometric B-spline (NCTB-spline) functions to approximate the regularized long-wave (RLW) equation. This equation is used for modelling numerous problems occurring in applied sciences. The NCTB-spline collocation method is used to integrate the spatial derivatives. We use the Rubin–Graves linearization technique to linearize the non-linear term. The accuracy and efficiency of the technique are examined by employing it on three important numerical examples which have three invariants of motion viz. mass, momentum, and energy. It is observed that the error norms of the present method are less than the error norms of the methods available in the literature. The numerical values of these invariants have also been approximated, which remain conserved during the program run which shows that the propagation of the solitary wave is represented perfectly. The propagation of one and two solitary waves and undulations of waves are depicted graphically. The stability analysis shows that the method is unconditionally stable.

1 Introduction

The regularized long-wave (RLW) equation is used as a model to describe long-wave behavior. Peregrine [1] introduced this equation to demonstrate the development of undular bore. Benjamin [2] discussed the properties of this equation. The aforesaid equation has been used to model various phenomena in science such as hydromagnetic waves, longitudinal dispersive waves, shallow water waves, anharmonic lattice, and rotating flow down a tube. Olver [3] proved that the aforesaid equation has only three conservation laws. The fixed conservation laws show that this equation is non-integrable. In fact, there are various initial and boundary conditions for that the general solution has not been obtained. Hence, there are only a few analytical solutions with a limited set of these conditions. So, there is a need for numerical methods for the RLW equation with different sets of initial and boundary conditions.

Various approximations have been found in the literature for the RLW equation with several sets of initial and boundary conditions. For instance, Gardner et al. [4] solved it numerically by the least square method based on linear space-time finite elements, whereas Dogan [5] used linear finite element-based Galerkin method. The authors of Ref. [6] used the splitting method in combination with the CB-spline method for solving non-linear RLW equation. The authors of Refs. [7,8] proposed linearized implicit and three-level finite difference methods, respectively, while a new finite difference method based on the quintic spline and splitting method has been proposed by Raslan [9]. Petrov-Galerkin finite element method, least-square quadratic finite element method and least-square cubic B-spline finite element method have been proposed by authors of Refs. [10,11,12], respectively, for solving the aforesaid equation. Saka et al. [13] and Zaki [14] used splitting methods in the combination of cubic and quadratic B-spline finite elements, respectively. The collocation method based on both quadratic and quintic B-splines has been proposed by Dağ et al. [15]. Mei and Chen [16] presented explicit multistep Galerkin method to solve this equation, whereas Dag and Dereli [17] investigated it numerically using radial basis-based meshless method. Eilbeck and McGuire [18] studied the solitary wave solutions numerically. Additionally, Cimpoiasu [19] investigated some traveling wave solutions for the long–short wave resonance model.

In Ref. [20], the authors presented a differential quadrature method based on cubic B-spline basis functions, whereas Saka and Dag [21] presented a collocation method based on quartic B-spline basis functions for the numerical solution of the aforesaid equation. The authors of Refs. [22,23,24,25] used Galerkin’s method using linear finite elements, quadratic B-spline finite element-based lumped Galerkin method, B-spline finite element method, and quadratic B-spline Galerkin finite element method in combination with space-splitting technique, respectively, to investigate the motion of a single solitary wave, development of two solitary wave interaction and an undular bore, numerically. The authors of Ref. [26] proposed a fourth-order collocation method based on cubic B-splines to solve RLW and modified RLW equations. Guo et al. [27] applied multiple integral FVM based on Lagrange interpolation for the Rosenau-RLW equation. Recently, Dhiman and Tamsir [28] applied a trigonometric cubic B-spline collocation technique for Fisher s reaction-diffusion equation.

In this article, we consider the RLW equation as follows:

(1) ( u μ u x x ) t + ( ε u + α ) u x = 0 , x [ a , b ] , t ( 0 , T ] ,

subject to

(2) u ( x , 0 ) = φ ( x ) , x [ a , b ] ,

and the boundary conditions

(3) u ( a , t ) = β 1 , u ( b , t ) = β 2 , u x ( a , t ) = 0 , u x ( b , t ) = 0 , t ( 0 , T ] ,

where μ , ε , α are known positive parameters and the wave amplitude u = u ( x , t ) is to be obtained. The terms β 1 and β 2 are given constants. Bona and Byrant [29] already proved that there exists a unique solution for the problem (1)–(3). The problem (1)–(3) conserve mass, momentum, and energy. These terms are physical quantities of motion, which are given by

(4) I 1 = a b u d x , I 2 = a b ( u 2 + μ u x 2 ) d x , I 3 = a b ( u 3 + 3 u x 3 ) d x .

The rest of the article is organized as follows. The procedure of the new trigonometric B-spline collocation technique is given in Section 2. In Section 3, a brief discretization of the problems using the proposed technique is given. The von Neumann stability is discussed in Section 4. In Section 5, numerical results are given. Finally, Section 6 focuses on the conclusions.

2 New cubic trigonometric B-spline (NCTB-spline) collocation technique

In this section, we introduce the collocation method based on NCTB-spline functions. First, we partition the problem domain x [ a , b ] into mesh of equal length h = x i + 1 x i for i = 0 , 1 , , M . For discrete form, we use u ( x i , t j ) = u i j for i = 0 , 1 , , M and j = 0 , 1 , , N . The NCTB-spline functions NB i ( x ) for i = 1 , 0 , , M , M + 1 are given by Ref. [30]

(5) NB i ( x ) = d i 2 ( g 3 ( z i 2 ) ) , x [ x i 2 , x i 1 ) , k = 0 3 c i 1 , k g k ( z i 1 ) , x [ x i 1 , x i ) , k = 0 3 b i , k g k ( z i 1 ) , x [ x i , x i + 1 ) , a i + 1 g 0 ( z i 1 ) , x [ x i + 1 , x i + 2 ) , 0 , else ,

where

z i ( x ) = π 2 x x i x i + 1 x i , i = 1 , 2 , , a i = b i , 3 = c i , 0 = d i = v l i 2 μ i ψ i , b i , 0 = v ψ i ( α i l i + β i l i 1 ) μ i ψ i 2 , b i , 1 = v β i ψ i , b i , 2 = v l i + 1 ψ i , c i , 1 = v l i 1 ψ i , c i , 2 = v α i ψ i , c i , 3 = v ψ i ( α i l i + 1 + β i l i ) μ i ψ i 2 ,

where

l i = x i + 1 x i , μ i = l i 1 + l i , α i = l i 1 + v l i , β i = v l i 1 + l i , v = 1 + 2 η , ψ i = v l i 2 + v 2 l i 1 + v l i .

The new trigonometric functions are defined as:

(6) g 0 ( z ) = ( 1 sin z ) 2 ( 1 + ( 1 η ) sin z ) , g 1 ( z ) = sin z ( 1 sin z ) ( 2 ( 1 η ) ( 1 sin z ) ) , g 2 ( z ) = cos z ( 1 cos z ) ( 2 ( 1 η ) ( 1 cos z ) ) , g 3 ( z ) = ( 1 cos z ) 2 ( 1 + ( 1 η ) cos z ) ,

where η 1 2 , 2 is the shape parameter. The set of NCTB-spline basis functions { NB i ( x ) : i = 1 , 0 , , M , M + 1 } forms a basis over the problem domain. Table 1 gives the values of NB i ( x ) and it is two derivatives at the knots.

Table 1

Values of NB i ( x ) , NB i ( x ) , and NB i ( x ) at the knots

x i 2 x i 1 x i x i + 1 x i + 2
NB i ( x ) 0 ξ 1 ξ 2 ξ 3 0
N B i ( x ) 0 ξ 4 ξ 5 ξ 6 0
NB i ( x ) 0 ξ 7 ξ 8 ξ 9 0

We assume that the approximation u M at ( x i , t j ) is written as a linear combination of the NCTB-spline functions and unknown time-dependent coefficients as:

(7) u M ( x , t ) = m = 1 M + 1 C m ( t ) NB m ( x ) .

The NCTB-spline and its four principle derivatives vanish outside of the region [ x i 2 , x i + 2 ] . Since each NCTB-spline covers four elements, each element is covered by four NCTB-splines. Hence, principle four NCTB-splines cover the interval [ x i , x i + 1 ]. So, the variation of the U M ( x , t ) , over the element, can be expressed as

(8) u M ( x , t ) = m = i 1 i + 1 C m ( t ) NB m ( x ) .

From Eq. (8) and Table 1, we get the approximated u ( x , t ) , u x ( x , t ) , and u x x ( x , t ) as follows:

(9) u i j = ξ 1 C i 1 j + ξ 2 C i j + ξ 3 C i + 1 j ,

(10) ( u x ) i j = ξ 4 C i 1 j + ξ 5 C i j + ξ 6 C i + 1 j ,

(11) ( u x x ) i j = ξ 7 C i 1 j + ξ 8 C i j + ξ 9 C i + 1 j ,

where C i j = C i ( t j ) .

3 Discretization of the RLW equation

This section considers the problem (1)–(3) when μ , ε , α are given. First, we discretize Eq. (1) as follows:

(12) ( u μ ( u x x ) ) i j + 1 ( u μ ( u x x ) ) i j Δ t + ε 2 ( ( u u x ) i j + 1 + ( u u x ) i j ) + α 2 ( ( u x ) i j + 1 + ( u x ) i j ) = 0 .

We linearize the non-linear term by Rubin–Graves technique [31] as follows:

(13) u j + 1 ( u x ) j + 1 = u j ( u x ) j + 1 + u j + 1 ( u x ) j u j ( u x ) j .

Using Eq. (13) in Eq. (12), we get

(14) ( u μ ( u x x ) ) i j + 1 + Δ t ε 2 ( u i j ( u x ) i j + 1 + u i j + 1 ( u x ) i j ) + Δ t α 2 ( u x ) i j + 1 = ( u μ ( u x x ) ) i j Δ t α 2 ( u x ) i j .

Simplifying the above equation, we get

(15) P i j u i j + 1 + Q i j ( u x ) i j + 1 μ ( u x x ) i j + 1 = u i j Δ t α 2 ( u x ) i j μ ( u x x ) i j ,

where

P i j = 1 + Δ t ε 2 ( u x ) i j , Q i j = Δ t ε 2 u i j + Δ t α 2 .

Now, using equations (9)–(11), we get

(16) ( P i j ξ 1 + Q i j ξ 4 μ ξ 7 ) C i 1 j + 1 + ( P i j ξ 2 + Q i j ξ 5 μ ξ 8 ) C i j + 1 + ( P i j ξ 3 + Q i j ξ 6 μ ξ 9 ) C i + 1 j + 1 = ξ 1 Δ t α 2 ξ 4 μ ξ 7 C i 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 C i j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 C i + 1 j , i = 1 , 2 , , M 1 , j = 0 , 1 , , N .

Now, we discretize the boundary conditions as follows:

(17) ξ 1 C 1 j + ξ 2 C 0 j + ξ 3 C 1 j = β 1 , j = 1 , 2 , , N ,

(18) ξ 1 C M 1 j + 1 + ξ 2 C M j + 1 + ξ 3 C M + 1 j + 1 = β 2 , j = 1 , 2 , , N .

From equations (17) and (18), we get

(19) C 1 j = β 1 ξ 1 ξ 2 ξ 1 C 0 j ξ 3 ξ 1 C 1 j , j = 1 , 2 , , N ,

(20) C M + 1 j = β 2 ξ 3 ξ 1 ξ 3 C M 1 j ξ 2 ξ 3 C M j , j = 1 , 2 , , N .

For i = 0 , using (19) in (16), we get

(21) P 0 j ξ 2 + Q 0 j ξ 5 μ ξ 8 ξ 2 ξ 1 ( P 0 j ξ 1 + Q 0 j ξ 4 μ ξ 7 ) C 0 j + 1 + ( P 0 j ξ 3 + Q 0 j ξ 6 μ ξ 9 ξ 3 ξ 1 ( P 0 j ξ 1 + Q 0 j ξ 4 μ ξ 7 ) ) C 1 j + 1 = ξ 2 Δ t α 2 ξ 5 μ ξ 8 ξ 2 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 C 0 j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 ξ 3 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 C 1 j + α 1 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 α 1 ξ 1 ( P 0 j ξ 1 + Q 0 j ξ 4 μ ξ 7 ) ) , j = 0 , 1 , , N .

For i = M , using (20) in (16), we get

(22) ( P M j ξ 1 + Q M j ξ 4 μ ξ 7 ) ξ 1 ξ 3 ( P M j ξ 3 + Q M j ξ 6 μ ξ 9 ) C M 1 j + 1 + ( ( P M j ξ 2 + Q M j ξ 5 μ ξ 8 ) ξ 2 ξ 3 ( P M j ξ 3 + Q M j ξ 6 μ ξ 9 ) ) C M j + 1 = ξ 1 Δ t α 2 ξ 4 μ ξ 7 ξ 1 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 C M 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 ξ 2 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 C M j + α 2 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 ( P M j ξ 3 + Q M j ξ 6 μ ξ 9 ) , j = 0 , 1 , , N .

At t j + 1 , j = 0 , 1 , , N , Eqs. (21), (16), and (22) can be written into the form of a linear system of order ( M + 1 ) × ( M + 1 ) as follows:

q ˆ 0 j ξ 2 ξ 1 p ˆ 0 j r ˆ 0 j ξ 3 ξ 1 p ˆ 0 j p ˆ 1 j q ˆ 1 j r ˆ 1 j p ˆ 2 j q ˆ 2 j r ˆ 2 j p ˆ M 2 j q ˆ M 2 j r ˆ M 2 j p ˆ M j ξ 1 ξ 3 r ˆ M j q ˆ M j ξ 2 ξ 3 r ˆ M j C 0 j + 1 C 1 j + 1 C 2 j + 1 C M 2 j + 1 C M 1 j + 1 C M j + 1 = R 0 j R 1 j R 2 j R M 2 j R M 1 j R M j ,

where

p ˆ i j = ( P i j ξ 1 + Q i j ξ 4 μ ξ 7 ) , q ˆ i j = ( P i j ξ 2 + Q i j ξ 5 μ ξ 8 ) , r ˆ i j = ( P i j ξ 3 + Q i j ξ 6 μ ξ 9 ) , R 0 j = ξ 2 Δ t α 2 ξ 5 μ ξ 8 ξ 2 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 C 0 j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 ξ 3 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 C 1 j + α 1 ξ 1 ξ 1 Δ t α 2 ξ 4 μ ξ 7 ( P 0 j ξ 1 + Q 0 j ξ 4 μ ξ 7 ) , j = 0 , 2 , , N , R i j ξ 1 Δ t α 2 ξ 4 μ ξ 7 C i 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 C i j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 C i + 1 j , i = 1 , 2 , , M 1 , j = 0 , 1 , , N , R M j = ξ 1 Δ t α 2 ξ 4 μ ξ 7 ξ 1 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 C M 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 ξ 2 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 C M j + α 2 ξ 3 ξ 3 Δ t α 2 ξ 6 μ ξ 9 ( P M j ξ 3 + Q M j ξ 6 μ ξ 9 ) , j = 0 , 1 , , N .

Now, we find the vector ( C 1 0 , C 0 0 , , C M 0 , C M + 1 0 ) from the initial condition which gives M + 1 equation in N + 3 unknowns. To eliminate the unknowns C 1 0 and C M + 1 0 , we use u x ( a , t ) = 0 and u x ( b , t ) = 0 , which gives

(23) C 1 0 = ξ 5 ξ 4 C 0 0 ξ 6 ξ 4 C 1 0 ,

(24) C M + 1 0 = ξ 4 ξ 6 C M 1 0 ξ 5 ξ 6 C M 0 .

Removing C 1 0 and C M + 1 0 from the system, we have a tridiagonal system of order ( M + 1 ) × ( M + 1 ) as follows:

ξ 2 ξ 1 ξ 5 ξ 4 ξ 3 ξ 1 ξ 6 ξ 4 ξ 1 ξ 2 ξ 3 ξ 1 ξ 2 ξ 3 ξ 1 ξ 2 ξ 3 ξ 1 ξ 1 ξ 4 ξ 6 ξ 2 ξ 1 ξ 5 ξ 6 C 0 0 C 1 0 C 2 0 C M 1 0 C M 0 = φ ( x 0 ) ξ 1 ξ 4 φ ( x 0 ) φ ( x 1 ) φ ( x 2 ) φ ( x M 1 ) φ ( x M ) ξ 1 ξ 6 φ ( x M ) .

4 Stability analysis

Now, the von-Neumann stability analysis [32,33,34,35,36] is carried out for the RLW equation. We consider the discretized system (12) as

(25) ( u μ ( u x x ) ) i j + 1 + Δ t ε 2 ( u i j ( u x ) i j + 1 + u i j + 1 ( u x ) i j ) + Δ t α 2 ( u x ) i j + 1 = ( u μ ( u x x ) ) i j Δ t α 2 ( u x ) i j .

Using equations (9)–(11), we get

(26) 1 + Δ t ε 2 ( u x ) i j ξ 1 + Δ t ε 2 u i j + Δ t α 2 ξ 4 μ ξ 7 C i 1 j + 1 + 1 + Δ t ε 2 ( u x ) i j ξ 2 + Δ t ε 2 u i j + Δ t α 2 ξ 5 μ ξ 8 C i j + 1 + 1 + Δ t ε 2 ( u x ) i j ξ 3 + Δ t ε 2 u i j + Δ t α 2 ξ 6 μ ξ 9 C i + 1 j + 1 = ξ 1 Δ t α 2 ξ 4 μ ξ 7 C i 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 C i j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 C i + 1 j .

Now, we assume u = k ¯ 1 and u x = k ¯ 2 as local constants for known level to get

(27) ( A ˆ ξ 1 + B ˆ ξ 4 μ ξ 7 ) C i 1 j + 1 + ( A ˆ ξ 2 + B ˆ ξ 5 μ ξ 8 ) C i j + 1 + ( A ˆ ξ 3 + B ˆ ξ 6 μ ξ 9 ) C i + 1 j + 1 = ξ 1 Δ t α 2 ξ 4 μ ξ 7 C i 1 j + ξ 2 Δ t α 2 ξ 5 μ ξ 8 C i j + ξ 3 Δ t α 2 ξ 6 μ ξ 9 C i + 1 j ,

where

A ˆ = 1 + Δ t ε 2 k 2 , B ˆ = Δ t ε 2 k 1 + Δ t α 2 .

Next, we consider C i j = δ j e k i ϕ as trial solutions at a given point x i , where ϕ = θ h and k = 1 . The θ is the mode number and h is the element size. Substituting the trial solution in the above equation, we get

(28) ( A ˆ 1 cos ( k ϕ ) + A ˆ 2 + A ˆ 3 cos ( k ϕ ) i ( A ˆ 1 + A ˆ 3 ) sin ( k ϕ ) ) δ B ˆ = B ˆ 1 cos ( k ϕ ) + B ˆ 2 + B ˆ 3 cos ( k ϕ ) i ( B ˆ 1 + B ˆ 3 ) sin ( k ϕ ) ,

where

A ˆ 1 = A ˆ ξ 1 + B ˆ ξ 4 μ ξ 7 , A ˆ 2 = A ˆ ξ 2 + B ˆ ξ 5 μ ξ 8 , A ˆ 3 = A ˆ ξ 3 + B ˆ ξ 6 μ ξ 9 , B ˆ 1 = ξ 1 Δ t α 2 ξ 4 μ ξ 7 , B ˆ 2 = ξ 2 Δ t α 2 ξ 5 μ ξ 8 , B ˆ 3 = ξ 3 Δ t α 2 ξ 6 μ ξ 9 .

Now, simplifying the terms in equation (28), we get

(29) δ 2 = ( B ˆ 1 B ˆ 2 + B ˆ 3 ) 2 + 4 ( B ˆ 1 + B ˆ 3 ) B ˆ 2 cos 2 ( k ϕ ) ( A ˆ 1 A ˆ 2 + A ˆ 3 ) 2 + 4 ( A ˆ 1 + A 3 ) A ˆ 2 cos 2 ( k ϕ ) ,

where

A ˆ 1 A ˆ 2 + A ˆ 3 = ( ξ 1 ξ 2 + ξ 3 ) Δ t α 2 ( ξ 4 ξ 5 + ξ 6 ) μ ( ξ 7 ξ 8 + ξ 9 ) + Δ t α 2 ( k 2 ( ξ 1 ξ 2 + ξ 3 ) + α ( k 1 + 2 ) ( ξ 4 ξ 5 + ξ 6 ) ) , B ˆ 1 B ˆ 2 + B ˆ 3 = ( ξ 1 ξ 2 + ξ 3 ) Δ t α 2 ( ξ 4 ξ 5 + ξ 6 ) μ ( ξ 7 ξ 8 + ξ 9 ) .

Since ξ 1 ξ 2 + ξ 3 > 0 and ξ 4 ξ 5 + ξ 6 > 0 , we get ( A ˆ 1 A ˆ 2 + A ˆ 3 ) 2 ( B ˆ 1 B ˆ 2 + B ˆ 3 ) 2 . Also, A ˆ 1 B ˆ 1 , A ˆ 2 B ˆ 2 , and A ˆ 3 B ˆ 3 . Therefore, we get δ 1 , and hence the method is unconditionally stable for the discretized system of the RLW equation.

5 Results and discussion

In this section, we consider three numerical examples in order to check the accuracy and efficiency of the proposed method. The L 2 and L errors are given by

(30) L 2 = h i = 0 M u i u i e 2 1 / 2 , L = max 0 i M u i u i e ,

where u i and u e are numerical and analytical solutions at x i .

5.1 Example 1

First, we consider the RLW equation (1) in x [ 40 , 60 ] with analytical solution [8,26,37]

(31) u ( x , t ) = 3 c sech 2 ( k x k v t k x 0 ) ,

and the boundary conditions (3), where β 1 = β 2 = 0 . This exact solution corresponds to the single solitary wave centrally located at x 0 initially with amplitude 3 c and with k , where k = 1 2 c ε μ v 1 2 . The term v = 1 + c ε stands for the wave velocity. The exact values of mass, momentum, and energy are given by

(32) I 1 = 6 c k , I 2 = 12 c 2 k + 48 k c 2 μ 5 , I 3 = 36 c 2 k + 144 c 3 5 k .

We fix ε = 1 , μ = 1 , α = 1 , and x 0 = 0 to approximate Example 1. First, we approximate the error norms for Example 1 with parameters M = 200 , Δ t = 0.001 for x = c = 0.1 and 0.3 at different t , which are presented in Table 2. It is clear from this table that the obtained results are better than the results presented by Mittal and Rohila [26]. Next, we computed the L error norms with M = 200 , Δ t = 0.1 , t = 20 , which are presented in Tables 3 and 4 for c = 0.1 and 0.3, respectively, together with the exact values of mass, energy, and momentum. We found from these tables that the present results are better than the results obtained by Dag and Ozer [12], Jain et al. [6], Kutluay and Esen [8], Mittal and Rohila [26], Raslan [9], two methods of Saka and Dag [21] viz. QBCM1 and QBCM2, and Zaki [14]. Moreover, during the program run, we observe that the approximated values of mass, energy, and momentum remain almost constants, i.e., propagation of the solitary wave represented accurately. The comparison between numerical and analytical solutions together with absolute error norms for M = 200 , Δ t = 0.1 for c = 0.1 at t = 4 , 8, and 16 are depicted in Figure 1. One can see from these figures that there is an excellent agreement between approximate and analytical solutions. Figure 2 shows the profiles of single solitary wave at various time with M = 200 , Δ t = 0.1 for c = 0.1 and 0.3, respectively, at t = 1 , 5 , 10 , 15 , and 20.

Table 2

Comparison of the present solutions with the solutions obtained by Mittal and Rohila [26], in terms of L , for Example 1

c = 0.1 c = 0.3
t Mittal and Rohila [26] Present Mittal and Rohila [26] Present
1 4.57 × 1 0 7 2.01 × 1 0 7 1.00 × 1 0 5 5.40 × 1 0 6
2 3.79 × 1 0 5 1.83 × 1 0 7 1.84 × 1 0 5 5.00 × 1 0 6
3 3.80 × 1 0 5 1.01 × 1 0 7 6.09 × 1 0 5 5.50 × 1 0 6
4 3.80 × 1 0 5 1.06 × 1 0 7 6.08 × 1 0 5 5.25 × 1 0 6
5 2.23 × 1 0 6 1.90 × 1 0 7 6.08 × 1 0 5 4.57 × 1 0 6
6 2.75 × 1 0 6 1.07 × 1 0 7 6.08 × 1 0 5 3.10 × 1 0 6
7 3.14 × 1 0 6 1.05 × 1 0 7 6.09 × 1 0 5 2.25 × 1 0 6
8 3.00 × 1 0 6 1.54 × 1 0 7 6.09 × 1 0 5 2.12 × 1 0 6
9 3.82 × 1 0 6 1.50 × 1 0 7 6.08 × 1 0 5 1.80 × 1 0 6
10 3.78 × 1 0 6 1.46 × 1 0 7 6.09 × 1 0 5 1.35 × 1 0 6
Table 3

Comparison of the present method and existing method results with c = 0.03 , for Example 1

Methods L I 1 I 2 I 3
Analytical 3.97992 0.81046 2.57900
Present method 1.12 × 1 0 6 3.97993 0.81046 2.57901
Dag and Ozer [12] 1.50 × 1 0 3 3.96466 0.80461 2.56971
Jain et al. [6] 6.70 × 1 0 2 4.41219 0.89734 2.85361
Kutluay and Esen [8] 2.10 × 1 0 4 3.97997 0.81045 2.57901
Mittal and Rohila [26] 7.81 × 1 0 5 3.97993 0.81046 2.57901
Raslan [9] 2.27 × 1 0 4 3.97995 0.80972 2.57607
QBCM1 [21] 8.30 × 1 0 5 3.97995 0.81046 2.57901
QBCM2 [21] 1.29 × 1 0 4 3.97995 0.81046 2.57901
Zaki [14] 1.81 × 1 0 4 3.97781 0.80963 2.5762
Table 4

L 2 and L error norms with c = 0.1 , Δ t = 0.1 , and Δ x = 0.125 , for Example 1

t QBSCM L 2 × 1 0 4 QBSCM L × 1 0 4 Present method L 2 × 1 0 4 Present method L × 1 0 4
0 0 0 0 0
4 0.1757 0.0693 0.1258 0.0548
8 0.2249 0.0887 0.1565 0.0652
12 0.3355 0.1072 0.2547 0.0923
16 0.4075 0.1224 0.3582 0.1567
20 0.4315 0.1321 0.3649 0.1825
Figure 1 
                  Numerical and analytical solutions together with absolute error norms at: (a) 
                        
                           
                           
                              t
                              =
                              4
                           
                           t=4
                        
                     , (b) 
                        
                           
                           
                              t
                              =
                              8
                           
                           t=8
                        
                     , and (c) 
                        
                           
                           
                              t
                              =
                              16
                           
                           t=16
                        
                     , for Example 1.
Figure 1

Numerical and analytical solutions together with absolute error norms at: (a) t = 4 , (b) t = 8 , and (c) t = 16 , for Example 1.

Figure 2 
                  Single solitary wave profiles with 
                        
                           
                           
                              M
                              =
                              200
                           
                           M=200
                        
                      for: (a) 
                        
                           
                           
                              c
                              =
                              0.1
                           
                           c=0.1
                        
                      and (b) 
                        
                           
                           
                              c
                              =
                              0.3
                           
                           c=0.3
                        
                     , for Example 1.
Figure 2

Single solitary wave profiles with M = 200 for: (a) c = 0.1 and (b) c = 0.3 , for Example 1.

5.2 Example 2

Now, we consider the RLW (1) in interval x [ 0 , 120 ] with the initial condition [8,9,26,37]

(33) u ( x , 0 ) = j = 1 2 3 C j sech 2 ( k j ( x x j ) ) ,

where C j = 4 k j 2 1 4 k j 2 , j = 1 , 2 . The RLW Eq. (1) with this initial condition represents two positive solitary wave equations. For the numerical computation of this example, we choose k 1 = 0.4 , k 2 = 0.3 , x 1 = 15 , x 2 = 35 , ε = 1 , α = 1 , and μ = 1 . Figure 3 depicts the propagation of two positive solitary waves with N = 300 for t = 0 , 5 , 10 , 15 , 20 , and 25, respectively. We also check the efficiency of the proposed method by computing the invariant mass, momentum, and energy. Tables 5 and 6 show the comparison of the computed invariants with the invariants presented in Refs. [8,26]. It is clear from these figures and tables that an excellent agreement is found between present and existing results. Also, it is noticed that the exact values of mass, energy, and momentum remain conserved during the program run which shows that the propagation of the solitary wave represented perfectly.

Figure 3 
                  Plot of the numerical solutions for 
                        
                           
                           
                              M
                              =
                              300
                           
                           M=300
                        
                      at different time levels: (a) 
                        
                           
                           
                              t
                              =
                              0
                           
                           t=0
                        
                     , (b) 
                        
                           
                           
                              t
                              =
                              5
                           
                           t=5
                        
                     , (c) 
                        
                           
                           
                              t
                              =
                              10
                           
                           t=10
                        
                     , (d) 
                        
                           
                           
                              t
                              =
                              15
                           
                           t=15
                        
                     , (e) 
                        
                           
                           
                              t
                              =
                              20
                           
                           t=20
                        
                     , and (f) 
                        
                           
                           
                              t
                              =
                              25
                           
                           t=25
                        
                     , for Example 2.
Figure 3

Plot of the numerical solutions for M = 300 at different time levels: (a) t = 0 , (b) t = 5 , (c) t = 10 , (d) t = 15 , (e) t = 20 , and (f) t = 25 , for Example 2.

Table 5

Numerical values of mass and momentum at different time levels, for Example 2

t I 1 (Mass) Present I 2 (Momentum) Present I 1 (Mass) [26] I 2 (Momentum) [26] I 1 (Mass) [8] I 2 (Momentum) [8]
0 37.9165 120.4126 37.9165 120.5220 37.9165 120.3515
2 37.9168 120.4115 37.9169 120.5230 37.9168 120.3571
4 37.9170 120.4356 37.9170 120.5240 37.9170 120.3584
6 37.9171 120.4360 37.9171 120.5320 37.9170 120.3586
8 37.9173 120.4275 37.9173 120.5560 37.9172 120.3583
10 37.9174 120.4596 37.9174 120.5950 37.9173 120.3570
12 37.9174 120.4675 37.9174 120.5340 37.9173 120.3915
14 37.9174 120.3570 37.9174 120.2570 37.9174 120.4156
16 37.9175 120.3045 37.9175 120.1940 37.9174 120.3886
18 37.9175 120.3495 37.9175 120.3570 37.9174 120.3653
20 37.9174 120.4169 37.9175 120.4620 37.9174 120.3599
22 37.9174 120.4168 37.9175 120.4620 37.9174 120.3599
24 37.9175 120.4215 37.9175 120.5160 37.9175 120.3595
Table 6

Numerical values of energy, for Example 2

t I 3 (Energy) Present I 3 (Energy) [26] I 3 (Energy) [8]
0 744.0811 744.0810 744.0814
2 744.0760 744.0770 744.0387
4 744.0560 744.0660 744.0110
6 744.0201 744.0320 744.9796
8 743.8815 743.8910 743.8679
10 743.4185 743.4130 743.4202
12 742.3056 742.2550 742.3387
14 741.5525 741.5420 741.5781
16 742.5314 742.5810 742.4890
18 7435232 743.5710 743.4752
20 743.9081 743.9290 743.8638
22 743.9081 743.9290 743.8638
24 743.8676 743.9175 744.0037

5.3 Example 3

Finally, we consider the RLW Eq. (1) with initial condition [8,26,37,38,39]

(34) u ( x , 0 ) = 0.5 β 1 1 tanh x x c d ,

to study the development of the undular bore together with the boundary conditions

(35) u ( a , t ) = β 1 , u ( b , t ) = 0 .

The term d denotes the slope between the shallow and deeper water, whereas u ( x , 0 ) = 0 represents the elevation of the water above the equilibrium surface at initial time. The change in water level of magnitude u ( x , 0 ) is centered at x = x c . We choose the parameters β 1 = 0.1 , ε = 1.5 , μ = 1 , α = 1 , x c = 0 , a = 36 , b = 300 , Δ t = 0.1 , Δ x = 0.24 , d = 2 , and 5 and we run the program until t = 250 . The initial functions at t = 0 are depicted in Figure 4 for d = 5 , 2 while the undulation profiles at t = 250 are shown in Figure 5 for d = 5 and 2, respectively. Graphs in Figure 5 represent the undular bore profiles when the gentle slope d = 5 and steep slope d = 2 , respectively, are used. A few more undulations can be seen when a steep slope, i.e., d = 0 is used as the initial condition. It is noticed that an increase in number of undulations is because of the form of the initial function. Also, we see that the magnitudes of leading undulations for the d = 5 and 2 become closer and closer for the longer run. An excellent agreement between the present and existing results [1,4,15] is found.

Figure 4 
                  Initial profile at 
                        
                           
                           
                              t
                              =
                              0
                           
                           t=0
                        
                      for: (a) 
                        
                           
                           
                              d
                              =
                              5
                           
                           d=5
                        
                      and (b) 
                        
                           
                           
                              d
                              =
                              2
                           
                           d=2
                        
                     , for Example 3.
Figure 4

Initial profile at t = 0 for: (a) d = 5 and (b) d = 2 , for Example 3.

Figure 5 
                  Undulation profiles at time 
                        
                           
                           
                              t
                              =
                              250
                           
                           t=250
                        
                      for: (a) 
                        
                           
                           
                              d
                              =
                              5
                           
                           d=5
                        
                      and (b) 
                        
                           
                           
                              d
                              =
                              2
                           
                           d=2
                        
                     , for Example 3.
Figure 5

Undulation profiles at time t = 250 for: (a) d = 5 and (b) d = 2 , for Example 3.

6 Conclusion

An NCTB-spline collocation technique has been proposed for the numerical computation of the RLW equation. The proposed collocation technique is used to integrate the spatial derivatives, whereas the discretization of the time derivative is done by the usual finite difference method. The numerical results have been presented in tabular form as well as graphically at different time levels. The present solutions have been compared with the solutions obtained by methods of Dag and Ozer [12], Jain et al. [6], Kutluay and Esen [8], Mittal and Rohila [26], Raslan [9], two methods of Saka and Dag [21] viz. QBCM1 and QBCM2, and Zaki [14] as well as with the exact solutions. It has been observed that the present error norms are less than the error norms presented in Refs. [6,8,9,12,14,21,26]. We noticed that the physical quantities of motion remain conserved during the program run, which shows that the propagation of the solitary wave is represented perfectly. The stability analysis of the discretized system of the RLW equation shows that the method is unconditionally stable.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

References

[1] Peregrine DH. Calculations of the development of an undular bore. J Fluid Mech. 1966;25:321–30. 10.1017/S0022112066001678Suche in Google Scholar

[2] Benjamin TB, Bona JL, Mahony JJ. Model equations for long waves in non-linear dispersive systems. Philos Trans R Soc. 1972;272:47–78. 10.1098/rsta.1972.0032Suche in Google Scholar

[3] Olver PJ. Euler operators and conservation laws of the BBM equation. Math Proc Cambridge Philos Soc. 1979;85:143–59. 10.1017/S0305004100055572Suche in Google Scholar

[4] Gardner LRT, Gardner GA, Dogan A. A least squares finite element scheme for the RLW equation. Commun Numer Methods Eng. 1996;12:795–804. 10.1002/(SICI)1099-0887(199611)12:11<795::AID-CNM22>3.0.CO;2-OSuche in Google Scholar

[5] Dogan A. Numerical solution of the RLW equation using linear finite elements within Galerkins method. Appl Math Model. 2002;26:771–83. 10.1016/S0307-904X(01)00084-1Suche in Google Scholar

[6] Jain PC, Shankar R, Singh TV. Numerical solutions of RLW equation. Commun Numer Methods Eng. 1993;9:579–86. 10.1002/cnm.1640090705Suche in Google Scholar

[7] Eilbeck JC, McGuire GR. Numerical study of the regularized long wave equation I: Numerical methods. J Comput Phys. 1975;19:43–57. 10.1016/0021-9991(75)90115-1Suche in Google Scholar

[8] Kutluay S, Esen A. A finite difference solution of the regularized long-wave equation. Math Probl Eng. 2006;2006:1–14. 10.1155/MPE/2006/85743Suche in Google Scholar

[9] Raslan KR. A computational method for the regularized long wave (RLW) equation. Appl Math Comput. 2005;167:1101–18. 10.1016/j.amc.2004.06.130Suche in Google Scholar

[10] Avilez-Valente P, Seabra-Santos FJ. A Petrov-Galerkin finite element scheme for the regularized long wave equation. Comput Mech. 2004;34:256–70. 10.1007/s00466-004-0570-4Suche in Google Scholar

[11] Dag I. Least squares quadratic B-spline finite element method for the regularized long wave equation. Comp Methods Appl Mech Eng. 2000;182:205–15. 10.1016/S0045-7825(99)00106-1Suche in Google Scholar

[12] Dag I, Ozer MN. Approximation of the RLW equation by the least square cubic B-spline finite element method. Appl Math Model. 2001;25:221–31. 10.1016/S0307-904X(00)00030-5Suche in Google Scholar

[13] Saka B, Dag I, Dogan I. Galerkin method for the numerical solution of the RLW equation using quadratic B-splines. Int J Comput Math. 2004;81:727–39. 10.1080/00207160310001650043Suche in Google Scholar

[14] Zaki SI. Solitary waves of the split RLW equation. Comp Phys Commun. 2001;138:80–91. 10.1016/S0010-4655(01)00200-4Suche in Google Scholar

[15] Dağ İ, DoGan A, Saka B. B-Spline Collocation Methods For Numerical Solutions Of The Rlw Equation. Int J Comp Math. 2003;80:743–57. 10.1080/0020716021000038965Suche in Google Scholar

[16] Mei L, Chen Y. Explicit multistep method for the numerical solution of RLW equation. Appl Math Comput. 2012;218:9547–54. 10.1016/j.amc.2012.03.050Suche in Google Scholar

[17] Dag I, Dereli Y. Numerical solution of RLW equation using radial basis functions. Int J Comp Math. 2010;87(1):63–76. 10.1080/00207160801965255Suche in Google Scholar

[18] Eilbeck JC, McGuire GR. Numerical study of the regularized long-wave equation II: Interaction of solitary waves. J Comput Phys. 1977;23:63–73. 10.1016/0021-9991(77)90088-2Suche in Google Scholar

[19] Cimpoiasu R. Travelling wave solutions for the Long-Short wave resonance model through an improved (G’/G)-expansion method. Rom J Phys. 2018;63:111. Suche in Google Scholar

[20] Korkmaz A, Dag I. Numerical simulations of boundary-forced RLW equation with cubic B-spline-based differential quadrature methods. Arab J Sci Eng. 2013;38:1151–60. 10.1007/s13369-012-0353-8Suche in Google Scholar

[21] Saka B, Dag I. Quartic B-spline collocation algorithms for numerical solution of the RLW equation. Numer Methods Partial Differ Equ. 2007;23:731–51. 10.1002/num.20201Suche in Google Scholar

[22] Dogan A. Application of Galerkin’s method to equal width wave equation. Appl Math Comput. 2005;160(1):65–76. 10.1016/j.amc.2003.08.105Suche in Google Scholar

[23] Esen A. A numerical solution of the equal width wave equation by a lumped Galerkin method. Appl Math Comput. 2005;168(1):270–82. 10.1016/j.amc.2004.08.013Suche in Google Scholar

[24] Raslan KR. A computational method for the equal width equation. Int J Comp Math. 2004;81(1):63–72. 10.1080/00207160310001614963Suche in Google Scholar

[25] Saka B. A finite element method for equal width equation. Appl Math Comput. 2006;175(1):730–47. 10.1016/j.amc.2005.07.034Suche in Google Scholar

[26] Mittal RC, Rohila R. A fourth order cubic B-spline collocation method for the numerical study of the RLW and MRLW equations. Wave Motion. 2018;80:47–68. 10.1016/j.wavemoti.2018.04.001Suche in Google Scholar

[27] Guo C, Li F, Zhang W, Luo Y. A conservative numerical scheme for Rosenau-RLW equation based on multiple integral finite volume method. Boundary Value Problems. 2019;2019:168.10.1186/s13661-019-1273-2Suche in Google Scholar

[28] Dhiman N, Tamsir M. A collocation technique based on modified form of trigonometric cubic B-spline basis functions for Fisher’s reaction-diffusion equation. Multidis Model Mater Struct. 2018;14:923–39. 10.1108/MMMS-12-2017-0150Suche in Google Scholar

[29] Bona JL, Bryant PJ. A mathematical model for long waves generated by wavemakers in non-linear dispersive systems. Math Proc Cambridge Philosoph Soc. 1973;73(2):391. 10.1017/S0305004100076945Suche in Google Scholar

[30] Majeed A, Qayyum F. New rational cubic trigonometric B-spline curves with two shape parameters. Comput Appl Math. 2020;39:198. 10.1007/s40314-020-01197-5Suche in Google Scholar

[31] Rubin SG, Graves RA. A cubic spline approximation for problems in fluid mechanics. Washington, DC: Nasa TR R-436; 1975. Suche in Google Scholar

[32] Caglar HN, Caglar SH, Twizell EH. The numerical solution of fifthorder boundary value problems with sixth-degree B-spline functions. Appl Math Lett 199912:25–30. 10.1016/S0893-9659(99)00052-XSuche in Google Scholar

[33] Chattaraj P, RaoKoneru S, Deb B. Stability analysis of finite difference schemes for quantum mechanical equations of motion. J Comput Phys. 1987;72:504–12. 10.1016/0021-9991(87)90098-2Suche in Google Scholar

[34] Dhiman N, Tamsir M. Re-modified quintic B-spline collocation method for the solution of Kuramoto-Sivashinsky type equations. Multidis Model Mater Struct. 2018. 10.1108/MMMS-06-2018-0111. Suche in Google Scholar

[35] O’Brien GG, Hyman MA, Kaplan S. A study of the numerical solution of partial differential equations. J Math Phys. 1950;29:223–51. 10.1002/sapm1950291223Suche in Google Scholar

[36] Vichnevetsky R. Stability charts in the numerical approximation of partial differential equations: a review. Math Comp Simulat. 1979;21:170–7. 10.1016/0378-4754(79)90130-7Suche in Google Scholar

[37] Dağ İ, Saka B. A Cubic B-Spline Collocation Method for the EW Equation. Math Comput Appl. 2004;9(3):381–92. 10.3390/mca9030381Suche in Google Scholar

[38] Esen A, Kutluay S. A linearized implicit finite-difference method for solving the equal width wave equation. Int J Comp Math. 2006;83(3):319–30. 10.1080/00207160600740958Suche in Google Scholar

[39] Gardner LRT, Gardner GA, Ayoub FA, Amein NK. Simulation of the EW undular bore. Commun Numer Methods Eng. 1997;13:583–92. 10.1002/(SICI)1099-0887(199707)13:7<583::AID-CNM90>3.0.CO;2-ESuche in Google Scholar

Received: 2021-09-02
Revised: 2021-11-06
Accepted: 2021-11-08
Published Online: 2021-12-06

© 2021 Ahmed Hussein Msmali et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Artikel in diesem Heft

  1. Regular Articles
  2. Circular Rydberg states of helium atoms or helium-like ions in a high-frequency laser field
  3. Closed-form solutions and conservation laws of a generalized Hirota–Satsuma coupled KdV system of fluid mechanics
  4. W-Chirped optical solitons and modulation instability analysis of Chen–Lee–Liu equation in optical monomode fibres
  5. The problem of a hydrogen atom in a cavity: Oscillator representation solution versus analytic solution
  6. An analytical model for the Maxwell radiation field in an axially symmetric galaxy
  7. Utilization of updated version of heat flux model for the radiative flow of a non-Newtonian material under Joule heating: OHAM application
  8. Verification of the accommodative responses in viewing an on-axis analog reflection hologram
  9. Irreversibility as thermodynamic time
  10. A self-adaptive prescription dose optimization algorithm for radiotherapy
  11. Algebraic computational methods for solving three nonlinear vital models fractional in mathematical physics
  12. The diffusion mechanism of the application of intelligent manufacturing in SMEs model based on cellular automata
  13. Numerical analysis of free convection from a spinning cone with variable wall temperature and pressure work effect using MD-BSQLM
  14. Numerical simulation of hydrodynamic oscillation of side-by-side double-floating-system with a narrow gap in waves
  15. Closed-form solutions for the Schrödinger wave equation with non-solvable potentials: A perturbation approach
  16. Study of dynamic pressure on the packer for deep-water perforation
  17. Ultrafast dephasing in hydrogen-bonded pyridine–water mixtures
  18. Crystallization law of karst water in tunnel drainage system based on DBL theory
  19. Position-dependent finite symmetric mass harmonic like oscillator: Classical and quantum mechanical study
  20. Application of Fibonacci heap to fast marching method
  21. An analytical investigation of the mixed convective Casson fluid flow past a yawed cylinder with heat transfer analysis
  22. Considering the effect of optical attenuation on photon-enhanced thermionic emission converter of the practical structure
  23. Fractal calculation method of friction parameters: Surface morphology and load of galvanized sheet
  24. Charge identification of fragments with the emulsion spectrometer of the FOOT experiment
  25. Quantization of fractional harmonic oscillator using creation and annihilation operators
  26. Scaling law for velocity of domino toppling motion in curved paths
  27. Frequency synchronization detection method based on adaptive frequency standard tracking
  28. Application of common reflection surface (CRS) to velocity variation with azimuth (VVAz) inversion of the relatively narrow azimuth 3D seismic land data
  29. Study on the adaptability of binary flooding in a certain oil field
  30. CompVision: An open-source five-compartmental software for biokinetic simulations
  31. An electrically switchable wideband metamaterial absorber based on graphene at P band
  32. Effect of annealing temperature on the interface state density of n-ZnO nanorod/p-Si heterojunction diodes
  33. A facile fabrication of superhydrophobic and superoleophilic adsorption material 5A zeolite for oil–water separation with potential use in floating oil
  34. Shannon entropy for Feinberg–Horodecki equation and thermal properties of improved Wei potential model
  35. Hopf bifurcation analysis for liquid-filled Gyrostat chaotic system and design of a novel technique to control slosh in spacecrafts
  36. Optical properties of two-dimensional two-electron quantum dot in parabolic confinement
  37. Optical solitons via the collective variable method for the classical and perturbed Chen–Lee–Liu equations
  38. Stratified heat transfer of magneto-tangent hyperbolic bio-nanofluid flow with gyrotactic microorganisms: Keller-Box solution technique
  39. Analysis of the structure and properties of triangular composite light-screen targets
  40. Magnetic charged particles of optical spherical antiferromagnetic model with fractional system
  41. Study on acoustic radiation response characteristics of sound barriers
  42. The tribological properties of single-layer hybrid PTFE/Nomex fabric/phenolic resin composites underwater
  43. Research on maintenance spare parts requirement prediction based on LSTM recurrent neural network
  44. Quantum computing simulation of the hydrogen molecular ground-state energies with limited resources
  45. A DFT study on the molecular properties of synthetic ester under the electric field
  46. Construction of abundant novel analytical solutions of the space–time fractional nonlinear generalized equal width model via Riemann–Liouville derivative with application of mathematical methods
  47. Some common and dynamic properties of logarithmic Pareto distribution with applications
  48. Soliton structures in optical fiber communications with Kundu–Mukherjee–Naskar model
  49. Fractional modeling of COVID-19 epidemic model with harmonic mean type incidence rate
  50. Liquid metal-based metamaterial with high-temperature sensitivity: Design and computational study
  51. Biosynthesis and characterization of Saudi propolis-mediated silver nanoparticles and their biological properties
  52. New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation
  53. Modal characteristics of harmonic gear transmission flexspline based on orthogonal design method
  54. Revisiting the Reynolds-averaged Navier–Stokes equations
  55. Time-periodic pulse electroosmotic flow of Jeffreys fluids through a microannulus
  56. Exact wave solutions of the nonlinear Rosenau equation using an analytical method
  57. Computational examination of Jeffrey nanofluid through a stretchable surface employing Tiwari and Das model
  58. Numerical analysis of a single-mode microring resonator on a YAG-on-insulator
  59. Review Articles
  60. Double-layer coating using MHD flow of third-grade fluid with Hall current and heat source/sink
  61. Analysis of aeromagnetic filtering techniques in locating the primary target in sedimentary terrain: A review
  62. Rapid Communications
  63. Nonlinear fitting of multi-compartmental data using Hooke and Jeeves direct search method
  64. Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger
  65. Knocking characteristics of a high pressure direct injection natural gas engine operating in stratified combustion mode
  66. What dominates heat transfer performance of a double-pipe heat exchanger
  67. Special Issue on Future challenges of advanced computational modeling on nonlinear physical phenomena - Part II
  68. Lump, lump-one stripe, multiwave and breather solutions for the Hunter–Saxton equation
  69. New quantum integral inequalities for some new classes of generalized ψ-convex functions and their scope in physical systems
  70. Computational fluid dynamic simulations and heat transfer characteristic comparisons of various arc-baffled channels
  71. Gaussian radial basis functions method for linear and nonlinear convection–diffusion models in physical phenomena
  72. Investigation of interactional phenomena and multi wave solutions of the quantum hydrodynamic Zakharov–Kuznetsov model
  73. On the optical solutions to nonlinear Schrödinger equation with second-order spatiotemporal dispersion
  74. Analysis of couple stress fluid flow with variable viscosity using two homotopy-based methods
  75. Quantum estimates in two variable forms for Simpson-type inequalities considering generalized Ψ-convex functions with applications
  76. Series solution to fractional contact problem using Caputo’s derivative
  77. Solitary wave solutions of the ionic currents along microtubule dynamical equations via analytical mathematical method
  78. Thermo-viscoelastic orthotropic constraint cylindrical cavity with variable thermal properties heated by laser pulse via the MGT thermoelasticity model
  79. Theoretical and experimental clues to a flux of Doppler transformation energies during processes with energy conservation
  80. On solitons: Propagation of shallow water waves for the fifth-order KdV hierarchy integrable equation
  81. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part II
  82. Numerical study on heat transfer and flow characteristics of nanofluids in a circular tube with trapezoid ribs
  83. Experimental and numerical study of heat transfer and flow characteristics with different placement of the multi-deck display cabinet in supermarket
  84. Thermal-hydraulic performance prediction of two new heat exchangers using RBF based on different DOE
  85. Diesel engine waste heat recovery system comprehensive optimization based on system and heat exchanger simulation
  86. Load forecasting of refrigerated display cabinet based on CEEMD–IPSO–LSTM combined model
  87. Investigation on subcooled flow boiling heat transfer characteristics in ICE-like conditions
  88. Research on materials of solar selective absorption coating based on the first principle
  89. Experimental study on enhancement characteristics of steam/nitrogen condensation inside horizontal multi-start helical channels
  90. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part I
  91. Numerical exploration of thin film flow of MHD pseudo-plastic fluid in fractional space: Utilization of fractional calculus approach
  92. A Haar wavelet-based scheme for finding the control parameter in nonlinear inverse heat conduction equation
  93. Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model
  94. Novel soliton solutions to the Atangana–Baleanu fractional system of equations for the ISALWs
  95. On the oscillation of nonlinear delay differential equations and their applications
  96. Abundant stable novel solutions of fractional-order epidemic model along with saturated treatment and disease transmission
  97. Fully Legendre spectral collocation technique for stochastic heat equations
  98. Special Issue on 5th International Conference on Mechanics, Mathematics and Applied Physics (2021)
  99. Residual service life of erbium-modified AM50 magnesium alloy under corrosion and stress environment
  100. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part I
  101. Diverse wave propagation in shallow water waves with the Kadomtsev–Petviashvili–Benjamin–Bona–Mahony and Benney–Luke integrable models
  102. Intensification of thermal stratification on dissipative chemically heating fluid with cross-diffusion and magnetic field over a wedge
Heruntergeladen am 16.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/phys-2021-0087/html
Button zum nach oben scrollen