Next Article in Journal
An Innovative Grey Approach for Group Multi-Criteria Decision Analysis Based on the Median of Ratings by Using Python
Next Article in Special Issue
Critical Indices and Self-Similar Power Transform
Previous Article in Journal
New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method
Previous Article in Special Issue
Classical Partition Function for Non-Relativistic Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach

1
Department of Mathematics, CHRIST, Bengaluru 560029, India
2
Faculty of Education, Harran University, Sanliurfa 63050, Turkey
3
School of Information Science and Technology, Yunnan Normal University, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(2), 123; https://doi.org/10.3390/axioms10020123
Submission received: 4 May 2021 / Revised: 6 June 2021 / Accepted: 13 June 2021 / Published: 16 June 2021
(This article belongs to the Special Issue Modern Problems of Mathematical Physics and Their Applications)

Abstract

:
The main objective of the present study is to analyze the nature and capture the corresponding consequences of the solution obtained for the Gardner–Ostrovsky equation with the help of the q -homotopy analysis transform technique ( q -HATT). In the rotating ocean, the considered equations exemplify strong interacting internal waves. The fractional operator employed in the present study is used in order to illustrate its importance in generalizing the models associated with kernel singular. The fixed-point theorem and the Banach space are considered to present the existence and uniqueness within the frame of the Caputo–Fabrizio (CF) fractional operator. Furthermore, for different fractional orders, the nature has been captured in plots. The realized consequences confirm that the considered procedure is reliable and highly methodical for investigating the consequences related to the nonlinear models of both integer and fractional order.

1. Introduction

The concept of fractional calculus (FC) has received much consideration in recent years due to its wide applicability and ability to capture more consequences of real-world problems. Even though the concept of classical calculus can exemplify most real-life problems and aid us in predicting the nature of complex phenomena, it has become a narrow subset of FC. Even though pioneers propose many new notions, many things need to be derived in order to ensure all classes of phenomena, which will be achieved by overcoming the limitations raised by mathematicians and scientists [1,2,3,4,5,6,7]. This is particularly true when researchers try to study, analyze and predict behaviors related to history, long-range memory, heritage, chaos, epidemiology, and other such subjects. There is always a door open for innovation, novelty, improvisation, and modifications in the research when it comes to the investigation of consequences that help us to solve real-world problems (the present pandemic, for example) and in this connection, many researchers have derived stimulating results with the aid of FC, and by using efficient techniques with the aid of fundamental results of FC [7,8,9,10].
The investigation of physical models and analysis of the associated properties are always a hot topic in applied mathematics with appropriate tools. In particular, this is the case for research on earthquakes, electrochemistry, signal processing, viscoelastic models, fluid dynamics, plasma physics, and many others. In the present investigation, we consider the nonlinear model describing the strong interacting internal waves in the ocean. In 1995, Oregon Bay experimented Coastal Ocean Probe Experiment (COPE) [11], found that the internal shallow waves are very strong, and also, it is stated that the shallow water wave phenomena are not completely nurtured with one nonlinear term. Furthermore, there is a need for dual nonlinear terms to examine the physical behaviors of the model. These experiments on shallow-water waves aid in formulating the model with dual-power law nonlinearity.
The equation that exemplifies the above phenomena is familiarly known as the Gardner equation. Generally, the Korteweg de-Vries (KdV) equations are widely used to exemplify the behaviors of long waves and their physical significance. Here, we consider the modified KdV equation, called the Gardner–Ostrovsky (GO) equation, due to the large amplitude for long internal waves and extended rotational effects [12]. Now, we consider the following two equations respectively exemplifying Ostrovsky and Gardner-Ostrovsky (GO) equations with the isopycnic surface ( u ( x ,   t ) ) :
d d x ( d u ( x , t ) d t + κ 1 u d u d x + β d 3 u d x 3 ) = c u ( x , t ) ,
d u ( x , t ) d t + κ 1 u d u d x + κ 2 u 2 d u d x + β d 3 u d x 3 = 0 ,
where κ 1 and κ 2 are, respectively, quadratic and cubic nonlinearities coefficients, β is the coefficient of small-scale dispersion, and c is the velocity of dispersion-less linear waves. Here, κ 1 and κ 2 are two proportional terms involving in the above equations, and u d u d x is due to the nonlinear hydrodynamic system, and it comes traditionally. Moreover, the term u 2 d u d x appears when u d u d x is a bird that is arbitrarily small. The GO equation is employed, and formally known to exemplify by large-amplitude waves with strong nonlinearity effects [13,14,15,16,17,18].
Recently, many authors illustrated the fractional operator’s essence in capturing and understanding the more relative consequences of physical phenomena with higher nonlinearity and complexity with long-range memory and history-based properties [19]. In this paper, we consider the fractional-order version of Equations (1) and (2) with Caputo–Fabrizio fractional operator with order α , fractional Ostrovsky (FO) and fractional Gardner–Ostrovsky (FGO) equations given by:
d d x ( D 0 C F t α u ( x , t ) + κ 1 u d u d x + β d 3 u d x 3 ) = c u ( x , t ) ,
D 0 C F t α u ( x , t ) + κ 1 u d u d x + κ 2 u 2 d u d x + β d 3 u d x 3 = 0 .
The vital and challenging job is finding the solution for nonlinear problems. Even though we have distinct schemes for finding the solution for these problems, they may possess a conversion of nonlinear to linear, partial to ordinary, perturbation, assumption of additional parameters, or discretization. Here, we consider the method proposed based on a topological concept called the homotopy analysis technique, which is nurtured to solve complicated nonlinear problems, and for 20 years, numerous scholarly works have been carried out with the aid of this [20,21]. Mankind is always searching for innovation, modernization and improvisation to increase accuracy and reduce the complexity associated with the method. In association with this, the authors in [22] came with the improvisation associated with Laplace transform (LT).
The projected scheme offers more freedom to choose the initial conditions, and the novelty is that it offers a simple solution procedure. This method is perceptible and includes all merits achieved by HAM, and also it attracted many researchers to analyze the diverse class of models and systems. The novelty of the projected scheme is the ability to solve the nonlinear problems without any assumption, perturbation, or conversion from nonlinear to linear, or partial to ordinary differential equations. Moreover, it is associated with parameters that are very helpful in converging the results attained towards a beneficial solution. Furthermore, it is related to well-posed transformation, which aims to reduce the complexity and increase the applicability and reliability of the method.
The considered technique ( q -HATT) is widely used for examining diverse classes of nonlinear problems and investigating the distinct mathematical models associated with real life [23,24,25,26,27,28,29,30]. For instance, the system of nonlinear differential equations exemplifying the outbreak of COVID-19 in India is examined within the frame of a nonsingular derivative in [24] using the projected scheme; HIV infection of lymphocyte cells is numerically analyzed by researchers in [25]; some interesting results have been derived with the help of the considered algorithm by authors in [29], in which the system exemplified the wind-influenced projectile motion; and the numerical simulation is presented by authors in [28] for the coupled Korteweg-de Vries system and this derived some interesting consequences in terms of numerical simulation.

2. Preliminaries

The basic definitions are presented in this segment for the CF and Laplace transform. Specifically, we recall the notions related to Caputo–Fabrizio fractional operator [19,31].
Definition 1. 
The CF fractional derivative for f H 1 ( a ,   b ) is presented as [19]:
D t α ( f ( t ) ) = ( α ) 1 α a t f ( t ) exp [ α t ϑ 1 α ] d ϑ ,   b > a ,
where ( α ) is a normalization function and admits ( 0 ) = ( 1 ) = 1 . Further, if f H 1 ( a ,   b ) , then we have:
D t α ( f ( t ) ) = α ( α ) 1 α t ( f ( t ) f ( ϑ ) ) exp [ α t ϑ 1 α ] d ϑ .
Definition 2. 
The CF fractional integralfor f H 1 ( a ,   b ) is presented as [31]:
I   C F α ( f ( t ) ) = 2 ( 1 α ) ( 2 α ) ( α ) f ( t ) + 2 α ( 2 α ) ( α ) 0 t f ( ϑ ) d ϑ , 0 < α < 1 ,   t 0 .
Note
According to [31], the following must hold:
2 ( 1 α ) ( 2 α ) ( α ) + 2 α ( 2 α ) ( α ) = 1 ,   0 < α < 1 ,
which gives ( α ) = 2 2 α . With the help of the above equation, researchers in [31] proposed a novel Caputo derivative as follows
D t α ( f ( t ) ) = 1 1 α 0 t f ( t ) exp [ α t ϑ 1 α ] d ϑ ,   0 < α < 1 .
Definition 3. 
The LT for a CF derivative D 0   C F t α f ( t ) is presented as [19] below:
[ D 0             C F C t ( α + n ) f ( t ) ] = s n + 1 [ f ( t ) ] s n f ( 0 ) s n 1 f ( 0 ) f ( n ) ( 0 ) s + ( 1 s ) α .

3. Fundamental Procedure of the Considered Method

Here, we hired the differential equation to present the basic algorithm of the projected method with the initial condition:
D 0 C F t α v ( x , t ) +   v ( x , t ) + N   v ( x , t ) = f ( x , t ) ,         0 < α 1 ,
and
v ( x , 0 ) = ( x ) .
We obtained this by applying LT on Equation (11):
[ v ( x , t ) ] ( x ) s + s + ( 1 s ) α s { [ v ( x , t ) ] + [ N   v ( x , t ) ] [ f ( x , t ) ] } = 0 .
For φ ( x , t ; q ) , N is contracted as follows:
N [ φ ( x , t ; q ) ] = [ φ ( x , t ; q ) ] ( x ) s + s + ( 1 s ) α s { [   φ ( x , t ; q ) ] + L [ N φ ( x , t ; q ) ] L [ f ( x , t ) ] }
where q [ 0 , 1 n ] . Then, the homotopy is defined by results in:
( 1 n q ) [ φ ( x ,   t ; q ) v 0 ( x , t ) ] = q N [ φ ( x ,   t ; q ) ] .
For q = 0 and q = 1 n , we have:
φ ( x , t ; 0 ) = v 0 ( x , t ) ,     φ ( x , t ; 1 n ) = v ( x , t ) .
By using the Taylor theorem, we consider:
φ ( x , t ; q ) = v 0 ( x , t ) + m = 1 v m ( x , t ) q m ,
where
v m ( x , t ) = 1 m ! m φ ( x ,   t ; q ) q m | q = 0 .
For the appropriate chaise of v 0 ( x , t ) ,   n   and the series (15) converges at q = 1 n . Then:
v ( x , t ) = v 0 ( x , t ) + m = 1 v m ( x , t ) ( 1 n ) m .
After differentiating Equation (15) m -times with q and multiplying by 1 m ! and substituting q = 0 , one can obtain:
[ v m ( x , t ) k m v m 1 ( x , t ) ] = m ( v m 1 ) ,
where
v m = { v 0 ( x , t ) , v 1 ( x , t ) , , v m ( x , t ) } .
Equation (20) reduces after employing inverse L T to:
v m ( x , t ) = k m v m 1 ( x , t ) + 1 [ m ( v m 1 ) ] ,
where
m ( v m 1 ) = L [ v m 1 ( x , t ) ] ( 1 k m n ) ( ( x ) s + s + ( 1 s ) α s L [ f ( x , t ) ] ) + s + ( 1 s ) α s L [ R v m 1 + m 1 ] ,
and
k m = { 0 ,       m 1 , n ,       m > 1 .
Here, m is homotopy polynomial and presented as:
m = 1 m ! [ m φ ( x ,   t ; q ) q m ] q = 0 and   φ ( x ,   t ; q ) = φ 0 + q φ 1 + q 2 φ 2 +   .
By the help of Equations (22) and (23), we found:
v m ( x , t ) = ( k m + ) v m 1 ( x , t ) ( 1 k m n ) 1 ( ( x ) s + s + ( 1 s ) α s L [ f ( x , t ) ] ) + 1 { s + ( 1 s ) α s L [ R v m 1 + m 1 ] } .
With the help of q -HATM, the series solution is:
v ( x , t ) = v 0 ( x , t ) + m = 1 v m ( x , t ) ( 1 n ) m .

4. Solution for FO and FGO Equations

The above-illustrated scheme is employed for the two fractional-order equations (namely, FO and FGO equations).
Case I
For FO equation
Consider the equation defined in Equation (3):
D 0 C F t α u ( x , t ) + κ 1 u d u d x + β d 3 u d x 3 c   u ( x , t ) d x = 0 ,
with
u ( x ,   0 ) = ρ sech 2 ( x ) .
Taking L T on Equation (28) and with the help of Equation (29), we obtain:
[ u ( x , t ) ] = 1 s ( ρ sech 2 ( x ) ) s + ( 1 s ) α s { κ 1 u d u d x + β d 3 u d x 3 c   u ( x , t ) d x } .
Now, N is defined as:
N [ φ ( x , t ; q ) ] = [ φ ( x , t ; q ) ] 1 s ( ρ sech 2 ( x ) ) + s + ( 1 s ) α s { κ 1 φ ( x , t ; q ) φ ( x , t ; q ) x + β d 3 φ ( x , t ; q ) d 3 x c   φ ( x , t ; q ) d x } .
At ( x , t ) = 1 , the deformation equation presented as:
[ u m ( x , t ) k m u m 1 ( x , t ) ] = m [ u m 1 ] ,
where
m [ u m 1 ] = [ u m 1 ( x , t ) ] ( 1 k m n ) { 1 s ( ρ sech 2 ( x ) ) } + s + ( 1 s ) α s { κ 1 i = 0 m 1 u i u m 1 i x + β d 3 u m 1 d x 3 c   u m 1 ( x , t ) d x } .
On employing inverse L T on Equation (32), it simplifies to:
u m ( x ,   t ) = k m u m 1 ( x ,   t ) + L 1 { m [ u m 1 ] } .
Now, by simplifying the above system we can evaluate the terms of the series solution:
u ( x ,   t ) = u 0 ( x ,   t ) + m = 1 u m ( x ,   t ) ( 1 n ) m .
Now, by the help of the initial condition, we can derive the terms of Equation (29) as:
u 1 ( x ,   t ) = ( 1 α + α t ) ( ρ tanh ( x ) ( c 16 β   sech 4 ( x ) + 2 ρ sech 4 ( x ) κ 1 + 8 β sech 2 ( x ) tanh 2 ( x ) ) ) ,

4.1. Existence of Solution for Ostrovsky Equation

Here, the existence and uniqueness are illustrated with CF operator for the considered Equation (28) as:
D 0 C F t α [ u ( x ,   t ) ] = Q 1 ( x ,   t ,   u ) .
Now, using Equation (36) and results derived in [31], we obtained:
u ( x ,   t ) u ( x ,   0 ) = I 0 C F t α { κ 1 u d u d x β d 3 u d x 3 + c   u ( x , t ) d x } .
Then we have from [31] the follows:
u ( x ,   t ) u ( x ,   0 ) = 2 ( 1 α ) ( α ) Q ( x , t , u ) + 2 α ( 2 α ) ( α ) 0 t Q 1 ( x ,   ζ , u ) d ζ .
Theorem 1. 
The kernel Q admits the Lipschitz condition and contraction if 0   ( κ 1 2 Λ ( a 1 + a 2 ) + β Λ 3 c ξ ) < 1 satisfies.
Proof. 
Consider the two functions A and A 1 to prove the theorem, then
Q 1 ( x ,   t ,   u ) Q 1 ( x ,   t ,   u 1 ) = κ 1 [ u ( x , t ) u ( x , t ) x u ( x , t 1 ) u ( x , t 1 ) x ] + β 3 x 3 [ u ( x , t ) u ( x , t 1 ) ] c b ( x ) = κ 1 2 [ x [ u 2 ( x , t ) u 2 ( x , t 1 ) ] ] + β 3 x 3 [ u ( x , t ) u ( x , t 1 ) ] c b ( x ) κ 1 2 Λ [ u ( x , t ) u ( x , t 1 ) ] + β Λ 3 c b ( x ) u ( x , t ) u ( x , t 1 ) ( κ 1 2 Λ ( a 1 + a 2 ) + β Λ 3 c ξ ) u ( x , t ) u ( x , t 1 ) ,
where a 1 = u and a 2 = u 1 are the bounded function and b ( x ) = ξ is also a bounded function. Set Ψ 1 = κ 1 2 Λ ( a 1 + a 2 ) + β Λ 3 c ξ in Equation (39), then:
Q 1 ( x ,   t , u ) Q 1 ( x ,   t , u 1 ) Ψ 1 u ( x ,   t ) u ( x ,   t 1 ) .
Equation (40) provides the Lipschitz condition for Q 1 . Similarly, we can see that if 0   κ 1 2 Λ ( a 1 + a 2 ) + β Λ 3 c ξ < 1 , then it implies the contraction. With the help of the above equations, Equation (38) simplifies to:
u ( x ,   t ) = u ( x ,   0 ) + 2 ( 1 α ) ( 2 α ) ( α ) Q 1 ( x ,   t , u ) + 2 α ( 2 α ) ( α ) 0 t Q 1 ( x ,   ζ , u ) d ζ .
u n ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) Q 1 ( x ,   t , u n 1 ) + 2 α ( 2 α ) ( α ) 0 t Q 1 ( x ,   ζ , u n 1 ) d ζ ,
Now, between the terms the successive difference is defined as:
ϕ 1 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x , t , u n 1 ) Q 1 ( x , t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x , t , u n 1 ) Q 1 ( x , t , u n 2 ) ) d ζ .
Notice that:
u n ( x ,   t ) = i = 1 n ϕ 1 i ( x ,   t ) .
Then we have:
ϕ 1 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x ,   t , u n 1 ) Q 1 ( x ,   t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x ,   t , u n 1 ) Q 1 ( x ,   t , u n 2 ) ) d ζ .
Application of the triangular inequality, Equation (45) reduces to:
ϕ 1 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x ,   t , u n 1 ) Q 1 ( x ,   t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x ,   t , u n 1 ) Q 1 ( x ,   t , u n 2 ) ) d ζ .
The Lipschitz condition satisfied by the kernel t 1 , then:
ϕ 1 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 ϕ ( n 1 ) ( x ,   t ) + 2 α ( 2 α ) ( α ) Ψ 1 0 t ϕ ( n 1 ) ( x ,   t ) d ζ .
By the aid of the above result, we state the following result:
Theorem 2. 
If we have specific t 0 , then the solution for Equation (28) will exist and be unique. Further, we have:
2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 + 2 α ( 2 α ) ( α ) Ψ 1 t 0 < 1 .
Proof. 
Let u ( x , t ) be the bounded functions admitting the Lipschitz condition. Then, we obtain, using Equations (46) and (47):
ϕ 1 i ( x ,   t ) u n ( x ,   0 ) [ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 + 2 α ( 2 α ) ( α ) Ψ 1 t ] n .
Therefore, for the obtained solution, continuity and existence are verified. Now, to prove the Equation (48) is a solution for Equation (28), we consider:
u ( x , t ) u ( x , 0 ) = u n ( x ,   t ) K 1 n ( x ,   t ) .
Let us consider:
K 1 n ( x , t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x ,   t , u ) Q 1 ( x ,   t ,   u n 1 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x ,   ζ , u ) Q 1 ( x ,   ζ , u n 1 ) ) d ζ   2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x ,   t , u ) Q 1 ( x ,   t ,   u n 1 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x ,   ζ , u ) Q 1 ( x ,   ζ , u n 1 ) ) d ζ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 u u n 1 + 2 α ( 2 α ) ( α ) Ψ 1 u u n 1 t .
This process gives:
K 1 n ( x , t ) ( 2 ( 1 α ) ( 2 α ) ( α ) + 2 α ( 2 α ) ( α ) t ) n + 1 Ψ 1 n + 1 M
Similarly, at t 0 we can obtain:
K 1 n ( x ,   t ) ( 2 ( 1 α ) ( 2 α ) ( α ) + 2 α ( 2 α ) ( α ) t 0 ) n + 1 Ψ 1 n + 1 M .
Next, for the solution of the projected model, we prove its uniqueness. Suppose u * ( x ,   t ) is another solution, then:
u ( x , t ) u * ( x , t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x , t , u ) Q 1 ( x , t ,   u * ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x , ζ , u ) Q 1 ( x , ζ , u * ) ) d ζ .
Now, employing the norm on the above equation, we obtain:
u ( x , t ) u * ( x , t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 1 ( x , t , u ) Q 1 ( x ,   t ,   u * ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 1 ( x , ζ , u ) Q 1 ( x , ζ , u * ) ) d ζ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 u ( x , t ) u * ( x , t ) + 2 α ( 2 α ) ( α ) Ψ 1 t u ( x , t ) u * ( x , t ) .
On simplification, this becomes:
u ( x , t ) u * ( x , t ) ( 1 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 2 α ( 2 α ) ( α ) Ψ 1 t ) 0 .
From the above condition, it is clear that u ( x , t ) = u * ( x , t ) , if:
( 1 2 ( 1 α ) ( 2 α ) ( α ) Ψ 1 2 α ( 2 α ) ( α ) Ψ 1 t ) 0 .
Hence, Equation (55) proves our required result. □
Case II
For the FGO equation
Consider the equation defined in Equation (4):
D 0 C F t α u ( x , t ) + κ 1 u d u d x + κ 2 u 2 d u d x + β d 3 u d x 3 = 0 ,
with condition defined in Equation (29). Now, taking L T on Equation (56) and by the assist of Equation (29), we obtain:
[ u ( x , t ) ] = 1 s ( ρ sech 2 ( x ) ) s + ( 1 s ) α s { κ 1 u d u d x + κ 2 u 2 d u d x + β d 3 u d x 3 } .
Now, N is defined as:
N [ φ ( x , t ; q ) ] = [ φ ( x , t ; q ) ] 1 s ( ρ sech 2 ( x ) ) + s + ( 1 s ) α s { κ 1 φ ( x , t ; q ) φ ( x , t ; q ) x + κ 2 φ 2 ( x , t ; q ) φ ( x , t ; q ) x + β d 3 φ ( x , t ; q ) d 3 x } .
Here:
m [ u m 1 ] = [ u m 1 ( x , t ) ] ( 1 k m n ) { 1 s ( ρ sech 2 ( x ) ) } + s + ( 1 s ) α s { κ 1 i = 0 m 1 u i u m 1 i x + κ 1 i = 0 m 1 j = 0 i u i j u i u m 1 i x + β d 3 u m 1 d x 3 } .

4.2. Existence of Solution for Fractional Gardner’s Ostrovsky Equation

Here, the existence and uniqueness are illustrated with CF operator for the considered Equation (56) as:
D 0 C F t α [ u ( x ,   t ) ] = Q 2 ( x ,   t ,   u ) .
Now, using Equation (56) and results derived in [31], we obtained:
u ( x ,   t ) u ( x ,   0 ) = I 0 C F t α { κ 1 u d u d x κ 2 u 2 d u d x β d 3 u d x 3 } .
Then, we have from [31] the following:
u ( x ,   t ) u ( x ,   0 ) = 2 ( 1 α ) ( α ) Q 2 ( x , t , u ) + 2 α ( 2 α ) ( α ) 0 t Q 2 ( x ,   ζ , u ) d ζ .

4.3. Existence of Solution for the Fractional Gardner–Ostrovsky Equation

Here, we will state the following results in a similar manner carried out for the above case.
Theorem 3. 
The kernel Q admits the Lipschitz condition and contraction if 0   ( κ 1 2 Λ ( a 1 + a 2 ) + κ 2 Λ ( a 1 + a 2 + a 1 a 2 ) + β Λ 3 ) < 1 satisfies for Equation (56).
Proof. 
Consider the two functions A and A 1 to prove the theorem, then:
Q 2 ( x ,   t ,   u ) Q 2 ( x ,   t ,   u 1 ) = κ 1 [ u ( x , t ) u ( x , t ) x u ( x , t 1 ) u ( x , t 1 ) x ] + κ 2 [ u 2 ( x , t ) u ( x , t ) x u 2 ( x , t 1 ) u ( x , t 1 ) x ] + β 3 x 3 [ u ( x , t ) u ( x , t 1 ) ] = κ 1 2 [ x [ u 2 ( x , t ) u 2 ( x , t 1 ) ] ] + κ 2 [ u 2 ( x , t ) u ( x , t ) x u 2 ( x , t 1 ) u ( x , t 1 ) x ] + β 3 x 3 [ u ( x , t ) u ( x , t 1 ) ] ( κ 1 2 Λ ( a 1 + a 2 ) + κ 2 Λ ( a 1 + a 2 + a 1 a 2 ) + β Λ 3 ) u ( x , t ) u ( x , t 1 ) ,
where a 1 = u and a 2 = u 1 is the bounded function. Set Ψ 2 = κ 1 2 Λ ( a 1 + a 2 ) + β Λ 3 c ξ in Equation (63), then:
Q 2 ( x ,   t , u ) Q 2 ( x ,   t , u 1 ) Ψ 2 u ( x ,   t ) u ( x ,   t 1 ) .
Equation (64) provides the Lipschitz condition for Q 2 . Similarly, we can see that if 0   κ 1 2 Λ ( a 1 + a 2 ) + κ 2 Λ ( a 1 + a 2 + a 1 a 2 ) + β Λ 3 < 1 , then it implies the contraction. By the assist of the above equation, Equation (62) simplifies to
u ( x ,   t ) = u ( x ,   0 ) + 2 ( 1 α ) ( 2 α ) ( α ) Q 2 ( x ,   t , u ) + 2 α ( 2 α ) ( α ) 0 t Q 2 ( x ,   ζ , u ) d ζ .
Then we obtain the recursive form as follows:
u n ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) Q 2 ( x ,   t , u n 1 ) + 2 α ( 2 α ) ( α ) 0 t Q 2 ( x ,   ζ , u n 1 ) d ζ ,
Now, between the terms, the successive difference is defined as:
ϕ 2 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x , t , u n 1 ) Q 2 ( x , t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x , t , u n 1 ) Q 2 ( x , t , u n 2 ) ) d ζ .
Notice that:
u n ( x ,   t ) = i = 1 n ϕ 2 i ( x ,   t ) .
Then, we have:
ϕ 2 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x ,   t , u n 1 ) Q 2 ( x ,   t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x ,   t , u n 1 ) Q 2 ( x ,   t , u n 2 ) ) d ζ
Application of the triangular inequality, Equation (69) reduces to:
ϕ 2 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x ,   t , u n 1 ) Q 2 ( x ,   t , u n 2 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x ,   t , u n 1 ) Q 2 ( x ,   t , u n 2 ) ) d ζ .
The Lipschitz condition satisfied by the kernel t 1 , then:
ϕ 2 n ( x ,   t ) = u n ( x ,   t ) u n 1 ( x ,   t ) 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 ϕ ( n 1 ) ( x ,   t ) + 2 α ( 2 α ) ( α ) Ψ 1 0 t ϕ ( n 1 ) ( x ,   t ) d ζ .
By the aid of the above result, we state the following result:
Theorem 4. 
If we have specific t 0 , then the solution for Equation (56) will exist and be unique. Further, we have:
Proof. 
Let u ( x , t ) be the bounded functions admitting the Lipschitz condition. Then, from Equations (70) and (71), we obtain:
ϕ 2 i ( x ,   t ) u n ( x ,   0 ) [ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 + 2 α ( 2 α ) ( α ) Ψ 2 t ] n
Therefore, for the obtained solution, continuity and existence are verified. Now, to prove that Equation (72) is a solution for Equation (56), we consider:
u ( x , t ) u ( x , 0 ) = u n ( x ,   t ) K 1 n ( x ,   t ) .
Let us consider:
K 1 n ( x ,   t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x ,   t , u ) Q 2 ( x ,   t ,   u n 1 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x ,   ζ , u ) Q 2 ( x ,   ζ , u n 1 ) ) d ζ   2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x ,   t , u ) Q 2 ( x ,   t ,   u n 1 ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x ,   ζ , u ) Q 2 ( x ,   ζ , u n 1 ) ) d ζ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 u u n 1 + 2 α ( 2 α ) ( α ) Ψ 2 u u n 1 t .
This process gives:
K 2 n ( x ,   t ) ( 2 ( 1 α ) ( 2 α ) ( α ) + 2 α ( 2 α ) ( α ) t ) n + 1 Ψ 2 n + 1 M
Similarly, at t 0 we can obtain:
K 2 n ( x ,   t ) ( 2 ( 1 α ) ( 2 α ) ( α ) + 2 α ( 2 α ) ( α ) t 0 ) n + 1 Ψ 2 n + 1 M .
Next, for the solution of the projected model, we prove the uniqueness. Suppose u * ( x ,   t ) is another solution, then:
u ( x , t ) u * ( x , t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x , t , u ) Q 2 ( x , t ,   u * ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x , ζ , u ) Q 2 ( x , ζ , u * ) ) d ζ .
Now, employing the norm on the above equation we obtain:
u ( x , t ) u * ( x , t ) = 2 ( 1 α ) ( 2 α ) ( α ) ( Q 2 ( x , t , u ) Q 2 ( x ,   t ,   u * ) ) + 2 α ( 2 α ) ( α ) 0 t ( Q 2 ( x , ζ , u ) Q 2 ( x , ζ , u * ) ) d ζ 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 u ( x , t ) u * ( x , t ) + 2 α ( 2 α ) ( α ) Ψ 2 t u ( x , t ) u * ( x , t ) .
On simplification:
u ( x , t ) u * ( x , t ) ( 1 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 2 α ( 2 α ) ( α ) Ψ 2 t ) 0 .
From the above condition, it is clear that u ( x , t ) = u * ( x , t ) , if
( 1 2 ( 1 α ) ( 2 α ) ( α ) Ψ 2 2 α ( 2 α ) ( α ) Ψ 2 t ) 0 .
Hence, Equation (79) proves our required result. □

5. Results and Discussion

Here, we capture the physical nature of two cases with different fractional order, small-scale dispersion, and homotopy parameters. It is essential to illustrate the effect of the fractional operator incorporated in the considered model. For Case I, we present the surfaces of fractional values of the order in Figure 1, and then it is illustrated in Figure 2 with a particular value of time ( t = 0.5 ) . In the Ostrovsky equation, the dispersion term plays a pivotal role and its effect is presented in Figure 3 with attained results. The main advantage of the considered scheme is the association of homotopy parameters (n and ). They play a critical role in the archiving solution for the modulation of convergence. For Case I we presented -curves in Figure 4.
Similarly, we present Case II (i.e., fractional Gardner–Ostrovsky equation). For the change in fractional values of the order, we capture the response in Figure 5 and Figure 6. The effect of the small dispersion term in the FGO equation is cited in Figure 7. In Figure 2, we have a diverse change in between the range of 1 and 1 with x . With the change in β within the range of 0.4 and 1 for this case, we can observe peak variations with the same values that we observe for x with fractional order in Figure 3, and similarly for Case II in Figure 6 and Figure 7. For Case I and II, we respectively presented -curves in Figure 8.These types of studies can influence researchers to investigate physical phenomena.

6. Conclusions

In this study, we investigated the equations exemplifying the rotating ocean with strong interacting internal waves with the fractional operator using q -HATT. With the help of the fixed-point hypothesis, the existence of the solution of the two cases is illustrated associated with the Caputo–Fabrizio operator. Furthermore, the effect of the coefficients of small-scale dispersion is illustrated with change in space x for both the cases, and we found noticeable variations in the archived results. The considered scheme is free from many limitations, including conversion from PDE to ODE, nonlinear to linear, or discretization. Finally, the present study demonstrates the effect of fractional order, parameters associated with models, and methods with their corresponding consequences.

Author Contributions

Conceptualization, methodology, software, writing: P.V.; validation, formal analysis: H.M.B.; investigation, resources: W.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This projected work was partially (not financial) supported by the Scientific Research Project Fund of Harran University with the project number HUBAP:21132.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liouville, J. Memoire surquelques questions de geometrieet de mecanique, et sur un nouveau genre de calcul pour resoudreces questions. J. Ec. Polytech. 1832, 13, 1–69. [Google Scholar]
  2. Riemann, G.F.B. Versuch Einer Allgemeinen Auffassung der Integration und Differentiation, Gesammelte Mathematische Werke; Leipzig, Germany, 1896. [Google Scholar]
  3. Caputo, M. Elasticita e Dissipazione; Zanichelli: Bologna, Italy, 1969. [Google Scholar]
  4. Miller, K.S.; Ross, B. An Introduction to Fractional Calculus and Fractional Differential Equations; A Wiley: New York, NY, USA, 1993. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  6. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  7. Baleanu, D.; Guvenc, Z.B.; Tenreiro Machado, J.A. New Trends in Nanotechnology and Fractional Calculus Applications; Springer: Dordrecht, The Netherlands; Heidelberg, Germany; London, UK; New York, NY, USA, 2010. [Google Scholar]
  8. Baishya, C.; Achar, S.J.; Veeresha, P.; Prakasha, D.G. Dynamics of a fractional epidemiological model with disease infection in both the populations. Chaos 2021, 31. [Google Scholar] [CrossRef]
  9. Akinyemi, L.; Veeresha, P.; Senol, M. Numerical solutions for coupled nonlinear schrodinger-Korteweg-de Vries and Maccari’s systems of equations. Mod. Phys. Lett. B 2021. [Google Scholar] [CrossRef]
  10. Baishya, C.; Jaipala. Numerical solution of fractional predator-prey model by trapezoidal based homotopy perturbation method. Int. J. Math. Arch. 2018, 9, 252–259. [Google Scholar]
  11. Antonova, M.; Biswas, A. Adiabatic parameter dynamics of perturbed solitary waves. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 734–748. [Google Scholar] [CrossRef]
  12. Holloway, P.; Pelinovsky, E.; Talipova, T. A generalised Korteweg-de Vries model of internal tide transformation in the coastal zone. J. Geophys. 1999, 104, 333–350. [Google Scholar] [CrossRef] [Green Version]
  13. Wazwaz, A.M. The variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients. Cent. Eur. J. Eng. 2014, 4, 64–71. [Google Scholar] [CrossRef]
  14. Apel, J.; Ostrovsky, L.A.; Stepanyants, Y.A.; Lynch, J.F. Internal solitons in the ocean and their effect on underwater sound. J. Acoust. Soc. Am. 2007, 121, 695–722. [Google Scholar] [CrossRef]
  15. Ostrovsky, L.A.; Pelinovsky, E.N.; Shrira, V.I.; Stepanyants, Y.A. Beyond the KDV: Post-explosion development. Chaos 2015, 25, 097620. [Google Scholar] [CrossRef]
  16. Biswas, A.; Krishnan, E.V. Exact solutions for Ostrovsky equation. Indian J. Phys. 2011, 85, 1513–1521. [Google Scholar] [CrossRef]
  17. Grimshaw, R.; Stepanyants, Y.; Alias, A. Formation of wave packets in the Ostrovsky equation for both normal and anomalous dispersion. Proc. R. Soc. A 2015, 472. [Google Scholar] [CrossRef]
  18. Stepanyants, Y.A. Nonlinear Waves in a rotating ocean the Ostrovsky equation and its generalizations and Applications. Atmos. Ocean. Phys. 2020, 56, 16–32. [Google Scholar] [CrossRef]
  19. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  20. Liao, S.J. Homotopy analysis method and its applications in mathematics. J. Basic Sci. Eng. 1997, 5, 111–125. [Google Scholar]
  21. Liao, S.J. Homotopy analysis method: A new analytic method for nonlinear problems. Appl. Math. Mech. 1998, 19, 957–962. [Google Scholar]
  22. Singh, J.; Kumar, D.; Swroop, R. Numerical solution of time- and space-fractional coupled Burgers’ equations via homotopy algorithm. Alex. Eng. J. 2016, 55, 1753–1763. [Google Scholar] [CrossRef] [Green Version]
  23. Srivastava, H.M.; Kumar, D.; Singh, J. An efficient analytical technique for fractional model of vibration equation. Appl. Math. Model. 2017, 45, 192–204. [Google Scholar] [CrossRef]
  24. Safare, K.M.; Betageri, V.S.; Prakasha, D.G.; Veeresha, P.; Kumar, S. A mathematical analysis of ongoing outbreak COVID-19 in India through nonsingular derivative. Numer. Methods Partial Differ. Equ. 2021, 37, 1282–1298. [Google Scholar] [CrossRef]
  25. Bulut, H.; Kumar, D.; Singh, J.; Swroop, R.; Baskonus, H.M. Analytic study for a fractional model of HIV infection of CD4+T lymphocyte cells. Math. Nat. Sci. 2018, 2, 33–43. [Google Scholar] [CrossRef]
  26. Prakasha, D.G.; Malagi, N.S.; Veeresha, P.; Prasannakumara, B.C. An efficient computational technique for time-fractional Kaup-Kupershmidt equation. Numer. Methods Partial Differ. Equ. 2021, 37, 1299–1316. [Google Scholar] [CrossRef]
  27. Kumar, D.; Agarwal, R.P.; Singh, J. A modified numerical scheme and convergence analysis for fractional model of Lienard’s equation. J. Comput. Appl. Math. 2018, 399, 405–413. [Google Scholar] [CrossRef]
  28. Akinyemi, L.; Huseen, S.N. A powerful approach to study the new modified coupled Korteweg-de Vries system. Math. Comput. Simul. 2021, 177, 556–567. [Google Scholar] [CrossRef]
  29. Veeresha, P.; Ilhan, E.; Baskonus, H.M. Fractional approach for analysis of the model describing wind-influenced projectile motion. Phys. Scr. 2021, 96, 075209. [Google Scholar] [CrossRef]
  30. Akinyemi, L.; Şenol, M.; Huseen, S.N. Modified homotopy methods for generalized fractional perturbed Zakharov–Kuznetsov equation in dusty plasma. Adv. Differ. Equ. 2021, 45. [Google Scholar] [CrossRef]
  31. Losada, J.; Nieto, J.J. Properties of the new fractional derivative without singular Kernel. Prog. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
Figure 1. Surfaces of q -HATT solution for (a) α = 0.50 , (b) α = 0.75 and (c) α = 1 at n = 1 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Figure 1. Surfaces of q -HATT solution for (a) α = 0.50 , (b) α = 0.75 and (c) α = 1 at n = 1 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Axioms 10 00123 g001
Figure 2. Response of the obtained solution with distinct α and time at n = 1 ,   t = 0.5 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Figure 2. Response of the obtained solution with distinct α and time at n = 1 ,   t = 0.5 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Axioms 10 00123 g002
Figure 3. Response of the obtained solution with distinct β and time at n = 1 ,   t = 0.5 ,   = 1 ,   α = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Figure 3. Response of the obtained solution with distinct β and time at n = 1 ,   t = 0.5 ,   = 1 ,   α = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case I.
Axioms 10 00123 g003
Figure 4. -curves for Case I at x = 1 ,   t = 0.01 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for (a) n = 1 and (b) n = 2 .
Figure 4. -curves for Case I at x = 1 ,   t = 0.01 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for (a) n = 1 and (b) n = 2 .
Axioms 10 00123 g004
Figure 5. Surfaces of q -HATT solution for (a) α = 0.50 , (b) α = 0.75 and (c) α = 1 at n = 1 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case II.
Figure 5. Surfaces of q -HATT solution for (a) α = 0.50 , (b) α = 0.75 and (c) α = 1 at n = 1 ,   = 1 ,   β = 1 ,   c = 1 ,   κ 1 = 1   and ρ = 0.001 for Case II.
Axioms 10 00123 g005
Figure 6. Response of the obtained solution with distinct α and time at n = 1 ,   t = 0.5 ,   = 1 ,   β = 1 ,   κ 1 = κ 2 = 1 and ρ = 0.001 for Case II.
Figure 6. Response of the obtained solution with distinct α and time at n = 1 ,   t = 0.5 ,   = 1 ,   β = 1 ,   κ 1 = κ 2 = 1 and ρ = 0.001 for Case II.
Axioms 10 00123 g006
Figure 7. Response of the obtained solution with distinct β and time at n = 1 ,   t = 0.5 ,   = 1 ,   α = 1 , κ 1 = κ 2 = 1 and ρ = 0.001 for Case II.
Figure 7. Response of the obtained solution with distinct β and time at n = 1 ,   t = 0.5 ,   = 1 ,   α = 1 , κ 1 = κ 2 = 1 and ρ = 0.001 for Case II.
Axioms 10 00123 g007
Figure 8. Curves for Case II at x = 1 ,   t = 0.01 ,   = 1 ,   β = 1 ,   κ 1 = κ 2 = 1 ,   n = 1 and ρ = 0.001 .
Figure 8. Curves for Case II at x = 1 ,   t = 0.01 ,   = 1 ,   β = 1 ,   κ 1 = κ 2 = 1 ,   n = 1 and ρ = 0.001 .
Axioms 10 00123 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Veeresha, P.; Baskonus, H.M.; Gao, W. Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach. Axioms 2021, 10, 123. https://doi.org/10.3390/axioms10020123

AMA Style

Veeresha P, Baskonus HM, Gao W. Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach. Axioms. 2021; 10(2):123. https://doi.org/10.3390/axioms10020123

Chicago/Turabian Style

Veeresha, Pundikala, Haci Mehmet Baskonus, and Wei Gao. 2021. "Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach" Axioms 10, no. 2: 123. https://doi.org/10.3390/axioms10020123

APA Style

Veeresha, P., Baskonus, H. M., & Gao, W. (2021). Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach. Axioms, 10(2), 123. https://doi.org/10.3390/axioms10020123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop