Next Article in Journal
Lyapunov Direct Method for Nonlinear Hadamard-Type Fractional Order Systems
Next Article in Special Issue
Deflection of Beams Modeled by Fractional Differential Equations
Previous Article in Journal
On the Approximate Solution of the Cauchy Problem in a Multidimensional Unbounded Domain
Previous Article in Special Issue
Second Derivative Block Hybrid Methods for the Numerical Integration of Differential Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Fractional Homotopy Perturbation Method for Solving Coupled Sine-Gordon Equations in Fractal Domain

1
School of Statistics and Mathematics, Inner Mongolia University of Finance and Economics, Hohhot 010070, China
2
School of Mathematical Sciences, Inner Mongolia University, Hohhot 010021, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(8), 404; https://doi.org/10.3390/fractalfract6080404
Submission received: 21 June 2022 / Revised: 14 July 2022 / Accepted: 19 July 2022 / Published: 22 July 2022
(This article belongs to the Special Issue Advances in Nonlinear Differential Equations with Applications)

Abstract

:
In this paper, the coupled local fractional sine-Gordon equations are studied in the range of local fractional derivative theory. The study of exact solutions of nonlinear coupled systems is of great significance for understanding complex physical phenomena in reality. The main method used in this paper is the local fractional homotopy perturbation method, which is used to analyze the exact traveling wave solutions of generalized nonlinear systems defined on the Cantor set in the fractal domain. The fractal wave with fractal dimension ε = l n 2 / l n 3 is numerically simulated. Through numerical simulation, we find that the obtained solutions are of great significance to explain some practical physical problems.

1. Introduction

In the field of science, most problems and phenomena are nonlinear, especially in mechanics [1], plasma physics [2], plasma waves, heat conduction [3] and chemical physics. Many scholars have made accurate analysis of various nonlinear problems and put forward some methods to solve them. Examples include Backlund transformation [4], bilinear method [5], sine-cosine method [6], Adomian decomposition method [7], variational iteration method [8], exp-function method [9] and homotopy perturbation method (HPM) [10,11]. We conclude that HPM is the most effective and simplest method to solve nonlinear problems.
HPM was first proposed by He [12]. He introduced the HPM into the solution of nonlinear wave equations. The advantage of this method is that it does not depend on the small parameters in the topology. There is a homotopy technique in topology that we can use to construct a homotopy with one parameter, and it is considered to be a small parameter [13]. HPM is widely used in nonlinear wave equations. It is used to solve water wave theory [14], heat conduction [15] and diffusion problems [16].
In recent years, Yang [17] redefined and generalized the local order, and constructed a complete framework of the local order defined on Cantor sets. Local fractions are usually used to describe various fractal problems in real life and science. Local fraction theory is a new efficient algorithm to obtain the exact solution of a local fraction problem. The solution of this theory is defined on the Cantor set. Yang used the local fractional homotopy perturbation (LFHP) method to solve the wave equations involving Cantor sets [18]. Zhang used the same method to solve the heat conduction equation in the fractal system [19]. The results show that the theoretical method is accurate and feasible, and has strong practical significance.
A coupled system of nonlinear equations is as follows:
ψ t t ψ x x = δ 2 sin ( ψ ω ) , ω t t c 2 ω x x = sin ( ψ ω ) ,
where δ ,   c are arbitrary constants, and δ > 0 ,   c > 0 .
The above coupled equations are an extension of the Frenkel–Kontorova model [20]. It simulates one-dimensional nonlinear wave processes in two-component media [21]. Coupled sine-Gordon equations with c = 1 were proposed to describe the open states in a DNA model. The coupled sine-Gordon equations can describe the propagation of optical pulses in fiber waveguide [22]. It can also describe long wave motion in shallow water. In recent years, many scholars have studied and analyzed the system (1). Salas gives its exact soliton solution and periodic solution [23]. Zhao et al. obtained the exact solution of the above equation by using the hyperbolic auxiliary function method [24]. Hosseini et al. solved the above coupling equations in nonlinear optics by using an improved Kudryashov method [22].
Our main purpose is to solve the traveling wave solutions of the following coupled local fractional sine-Gordon equations by using the LFHP method.
2 ε ψ t 2 ε 2 ε ψ x 2 ε = δ 2 sin ε ψ ω ε , 2 ε ω t 2 ε c 2 2 ε ω x 2 ε = sin ε ψ ω ε ,
where ε is the fractal dimension and 0 < ε 1 .
In this paper, the basic idea of the local fractional homotopy perturbation method is introduced and the traveling wave solutions of coupled local fractional sine-Gordon equations are obtained by using this method. The structure is as follows: The Section 2 introduces the basic theory and operation formula of local fractional calculus. The Section 3 introduces the basic idea and operation process of the LFHP method. In Section 4, we obtain the traveling wave solutions of coupled local fractional sinusoidal Gordon equations using the LFHP method, and obtain their wave images using MATLAB. Finally, we summarize the main conclusions in Section 5.

2. Local Fractional Calculus

2.1. Local Fractional Derivatives

Definition 1.
We define C ε ( μ , ν ) as a set of the nondifferentiable functions. Setting Ψ ( χ ) C ε ( μ , ν ) , local fractional derivative of Ψ ( χ ) of order ε ( 0 < ε 1 ) at χ = χ 0 is defined by [25]
D ( ε ) Ψ ( χ ) = d ε Ψ ( χ ) d χ ε χ = χ 0 = lim χ χ 0 Δ ε Ψ ( χ ) Ψ χ 0 χ χ 0 ε ,
where Δ ε Ψ ( χ ) Ψ χ 0 Γ ( 1 + ε ) Ψ ( χ ) Ψ χ 0 .
In addition, the local higher order fractional partial derivative is defined as
k ε χ k ε Ψ ( χ , γ ) = ε χ ε ε χ ε k times Ψ ( χ , γ ) .
The properties of the local fractional derivative are listed as follows [26]:
( R 1 ) D ( ε ) Ψ χ 0 ± χ 0 = D ( ε ) Ψ χ 0 ± D ( ε ) χ 0 , ( R 2 ) D ( ε ) Ψ χ 0 χ 0 = Ψ χ 0 D ( ε ) Ψ χ 0 + Ψ χ 0 D ( ε ) χ 0 , ( R 3 ) D ( ε ) Ψ χ 0 χ 0 = D ( ε ) Ψ χ 0 χ 0 Ψ χ 0 D ( ε ) χ 0 2 χ 0 ,
where ( χ 0 ) 0 , and Ψ ( χ 0 ) , ( χ 0 ) C ε ( μ , ν ) .

2.2. Local Fractional Integral

Definition 2.
The local fraction integral of Φ ( χ ) is defined as [27]
α I β ( ε ) Φ ( χ ) = 1 Γ ( 1 + ε ) μ ν Φ ( χ ) ( d χ ) ε = 1 Γ ( 1 + ε ) lim Δ χ σ 0 σ = 0 N 1 Φ χ σ Δ χ σ ε ,
where Φ ( χ ) C ε ( μ , ν ) , and Δ χ σ = χ σ + 1 χ σ with χ 0 = α < χ 1 < < χ N 1 < χ N = β .
The properties of the local fractional integral are as follows:
( S 1 ) 1 Γ ( 1 + ε ) μ x D ( ε ) Φ ( χ ) ( d χ ) ε = φ ( x ) Φ ( μ ) , ( S 2 ) D ( ε ) 1 Γ ( 1 + ε ) α x Φ ( χ ) ( d χ ) ε = Φ ( x ) .

2.3. Basic Operation

It defines the generalized function by [28]
E ε χ ε = κ = 0 χ κ ε Γ ( 1 + κ ε ) , sin ε χ ε = κ = 0 ( 1 ) κ χ ( 2 κ + 1 ) ε Γ ( 1 + ( 2 κ + 1 ) ε ) , cos ε χ ε = κ = 0 ( 1 ) κ χ 2 κ ε Γ ( 1 + 2 κ ε ) , sinh ε χ ε = E ε χ ε E ε χ ε 2 = κ = 0 χ ( 2 κ + 1 ) ε Γ ( 1 + ( 2 κ + 1 ) ε ) , cosh ε χ ε = E ε χ ε + E ε χ ε 2 = κ = 0 χ 2 κ ε Γ ( 1 + 2 κ ε ) .
Local fractional calculus has the following properties [28]
d ε d χ ε E ε C χ ε = C E ε C χ ε , d ε d χ ε sin ε C χ ε = C cos ε C χ ε , d ε d χ ε cos ε C χ ε = C sin ε C χ ε , d ε d χ ε sinh ε C χ ε = C cosh ε C χ ε , d ε d χ ε cosh ε C χ ε = C sinh ε C χ ε , d ε d χ ε χ κ ε Γ ( 1 + κ ε ) = χ ( κ 1 ) ε Γ ( 1 + ( κ 1 ) ε ) ,
and
0 I χ ( ε ) χ κ ε Γ ( 1 + κ ε ) = χ ( κ + 1 ) ε Γ ( 1 + ( κ + 1 ) ε ) .

3. Local Fractional Homotopy Perturbation Method

For a given class of local fractional differential equations,
L σ ( φ σ ) = 0 , φ R ,
where L σ is a local fractional differential operator.
First, let us set up a homotopy mapping [18]:
H σ ( φ , ρ ) = 1 ρ σ L σ φ σ L σ φ 0 σ + ρ σ L σ φ σ , φ R , ρ [ 0 , 1 ] ,
or
H σ ( φ , ρ ) = L σ φ σ L σ φ 0 σ + ρ σ L σ φ 0 σ , φ R , ρ [ 0 , 1 ] ,
where ρ is an imbedding parameter and φ 0 is an initial approximation of Equation (5).
Let H σ ( φ , ρ ) = 0 , then we can obtain, from Equation (6),
H σ ( u , 0 ) = L σ φ σ L σ φ 0 σ = 0 ; H σ ( φ , 1 ) = L σ φ 0 σ = 0 .
According to the homotopy perturbation theory, ρ can be treated as a small parameter. Suppose that the solution φ of Equation (6) can be expressed as a power series of ρ
φ σ = φ 0 σ + ρ σ φ 1 σ + ρ 2 σ φ 2 σ + ρ 2 σ φ 2 σ + = i = 0 n ρ i σ φ i σ .
Substituting Equation (9) into Equation (6), we obtain
H σ i = 0 n ρ i φ i , p = 1 ρ σ L σ i = 0 n ρ i φ i L σ φ 0 + ρ σ L σ i = 0 n ρ i φ i .
An extended form of L σ φ σ [19]
L σ φ σ = L σ φ 0 σ + d σ L σ φ 0 σ d φ σ i = 0 n ρ i φ i φ 0 Γ ( 1 + σ ) + O i = 0 n ρ i φ i φ 0 σ = L σ φ 0 σ + d σ L σ φ 0 σ d φ σ i = 0 n ρ i σ φ i σ φ 0 σ Γ ( 1 + σ ) + O i = 0 n ρ i φ i φ 0 σ ,
such that
H σ ( φ , ρ , σ ) = ( 1 ρ ) σ L σ φ σ L σ φ 0 σ + ρ σ L σ φ σ = ( 1 ρ ) σ L σ φ 0 σ + d σ L σ φ 0 σ d φ σ i = 0 n ρ i σ φ i σ φ 0 σ Γ ( 1 + σ ) + O i = 0 n ρ i φ i φ 0 σ L σ φ 0 σ + ρ σ L σ φ 0 σ + d σ L σ φ 0 σ d φ σ j = 0 n ρ j σ φ i σ φ 0 σ Γ ( 1 + σ ) + O j = 0 n ρ j φ i φ 0 σ ,
which reduces to
H σ ( φ , 0 ) = L σ φ σ L σ φ 0 σ = d σ L σ φ 0 σ d φ σ i = 0 n ρ i φ i φ 0 σ Γ ( 1 + σ ) + d 2 σ L σ φ 0 σ d φ 2 σ i = 0 n ρ i φ i φ 0 2 σ Γ ( 1 + 2 σ ) + d n σ L σ φ 0 σ d φ n σ i = 0 n ρ i φ i φ 0 n σ Γ ( 1 + n σ ) + = 0 ,
and
H σ ( φ , 1 ) = L σ φ σ = L σ φ 0 σ + d σ L σ φ 0 σ d φ σ i = 0 n ρ i φ i φ 0 σ Γ ( 1 + σ ) + d 2 σ L σ φ 0 σ d φ 2 σ i = 0 n ρ i φ i φ 0 2 σ Γ ( 1 + 2 σ ) + d n σ L σ φ 0 σ d φ n σ i = 0 n ρ i φ i φ 0 n σ Γ ( 1 + n σ ) + = 0 .
Using the expression Equation (10), we find that
ρ 0 σ : L σ φ σ L σ φ 0 σ = 0 ,
ρ 1 σ : d σ L σ φ 0 σ d φ σ φ 1 σ Γ ( 1 + σ ) + L σ φ 0 σ = 0 ,
ρ 2 σ : d σ L σ φ 0 σ d φ σ φ 2 σ Γ ( 1 + σ ) + d 2 σ L σ φ 0 σ d φ 2 σ φ 1 2 σ Γ ( 1 + 2 σ ) = 0 .
when ρ 1 , we have the approximate solution of the form:
φ σ = lim ρ 1 i = 0 n ρ i σ φ i σ = i = 0 n φ i σ .

4. Solution of System of Equations

In this section, the solutions of coupled local fractional sine-Gordon equations are discussed.
We consider the following system of local fractional equations:
2 ε ψ t 2 ε 2 ε ψ x 2 ε = δ 2 sin ε ψ ω ε , 2 ε ω t 2 ε c 2 2 ε ω x 2 ε = sin ε ψ ω ε .
In order to solve the coupled local part sine-Gordon equations, we can construct the following nondifferentiable homotopy:
1 ρ ε 2 ε ψ t 2 ε + ρ ε 2 ε ψ t 2 ε 2 ε ψ x 2 ε + δ 2 sin ε ( ψ ω ) = 0 , 1 ρ ε 2 ε ω t 2 ε + ρ ε 2 ε ψ t 2 ε c 2 2 ε ω x 2 ε sin ε ( ψ ω ) = 0 .
Moreover,
ψ = j = 0 ρ j ε ψ j = ψ 0 + ρ ε ψ 1 + ρ 2 ε ψ 2 + , ω = j = 0 ρ j ε ω j = ω 0 + ρ ε ω 1 + ρ 2 ε ω 2 + .
Submitting Equation (13) into Equation (12), we have
2 ε t j = 0 ρ j ε ψ j ( x , t ) 2 ε ψ 0 ( x , t ) t 2 ε = ρ ε 2 ε x 2 ε j = 0 ρ j ε ψ j ( x , t ) δ 2 k = 0 ( 1 ) k ( ψ w ) 2 k + 1 Γ ( 1 + ( 2 k + 1 ) ε ) 2 ε ψ 0 ( x , t ) t 2 ε , 2 ε t 2 ε j = 0 ρ j ε ω j ( x , t ) 2 ε ω 0 ( x , t ) t 2 ε = ρ ε c 2 2 ε x 2 ε j = 0 ρ j ε ω j ( x , t ) + k = 0 ( 1 ) k ( ψ ω ) 2 k + 1 Γ ( 1 + ( 2 k + 1 ) ε ) 2 ε ω 0 ( x , t ) t 2 ε .
Taking k = 1 , 2 , 3 , then rearranging according to the power of ρ -terms, we obtain:
ρ 0 ε : 2 ε ψ 0 ( x , t ) t 2 ε 2 ε ψ 0 ( x , t ) t 2 ε = 0 , ψ 0 ( x , t ) = A cos ε ( k x ε ) , 2 ε ω 0 ( x , t ) t 2 ε ω 0 ( x , t ) t 2 ε = 0 , ω 0 ( x , t ) = 0 .
ρ 1 ε : 2 ε ψ 1 t 2 ε = 2 ε ψ 0 x 2 ε δ 2 ψ 0 ω 0 Γ ( 1 + ε ) + δ 2 ψ 0 3 3 ψ 0 2 ω 0 + 3 ψ 0 ω 0 3 ω 0 3 Γ ( 1 + 3 ε ) δ 2 ψ 0 5 5 ψ 0 4 ω 0 + 10 ψ 0 3 ω 0 2 10 ψ 0 2 ω 0 3 + 5 ψ 0 2 ω 0 3 + 5 ψ 0 ω 0 4 ω 0 5 Γ ( 1 + 5 ε ) 2 ε ψ 0 t 2 ε , 2 ε ω 1 t 2 ε = c 2 2 ω 0 x 2 ε + ψ 0 ω 0 Γ ( 1 + ε ) ψ 0 3 3 ψ 0 2 ω 0 + 3 ψ 0 ω 0 3 ω 0 3 Γ ( 1 + 3 ε ) + ψ 0 5 5 ψ 0 4 ω 0 + 10 ψ 0 3 ω 0 2 10 ψ 0 2 ω 0 3 + 5 ψ 0 2 ω 0 3 + 5 ψ 0 ω 0 4 ω 0 5 Γ ( 1 + 5 ε ) 2 ε ω 0 t 2 ε .
ρ 2 ε : 2 ε ψ 2 t 2 ε = 2 ε ψ 1 x 2 ε δ 2 ψ 1 ω 1 Γ ( 1 + ε ) 1 Γ ( 1 + 3 ε ) 3 ψ 0 2 ψ 1 3 ψ 0 2 ω 1 6 ψ 0 ψ 1 ω 0 + 6 ψ 0 ω 0 ω 1 + 3 ψ 1 ω 0 2 3 ω 0 2 ω 1 + 1 Γ ( 1 + 5 ε ) 2 ψ 0 3 ψ 1 + 2 ψ 0 4 ψ 1 10 ψ 0 ψ 1 ω 0 5 ψ 0 3 ψ 1 ω 0 5 ψ 0 4 ω 1 + 30 ψ 0 2 ψ 1 ω 0 2 + 20 ψ 0 3 ω 0 ω 1 30 ψ 0 2 ω 0 2 ω 1 20 ψ 0 ψ 1 ω 0 3 + 10 ψ 0 ω 0 ω 1 + 5 ψ 0 ω 0 3 ω 1 + 5 ψ 1 ω 0 4 2 ω 0 3 ω 1 2 ω 0 4 ω 1 , 2 ε ω 2 t 2 ε = c 2 2 ε ω 1 x 2 ε + ψ 1 ω 1 Γ ( 1 + ε ) 1 Γ ( 1 + 3 ε ) 3 ψ 0 2 ψ 1 3 ψ 0 2 ω 1 6 ψ 0 ψ 1 ω 0 + 6 ψ 0 ω 0 ω 1 + 3 ψ 1 ω 0 2 3 ω 0 3 ω 1 + 1 Γ ( 1 + 5 ε ) 2 ψ 0 3 ψ 1 + 2 ψ 0 4 ψ 1 10 ψ 0 ψ 1 ω 0 5 ψ 0 3 ψ 1 ω 0 5 ψ 0 4 ω 1 + 30 ψ 0 2 ψ 1 ω 0 2 + 20 ψ 0 3 ω 0 ω 1 30 ψ 0 2 ω 0 2 ω 1 20 ψ 0 ψ 1 ω 0 3 + 10 ψ 0 ω 0 ω 1 + 5 ψ 0 ω 0 3 ω 1 + 5 ψ 1 ω 0 4 2 ω 0 3 ω 1 2 ω 0 4 ω 1 .
Solving Equations (15)–(17), we obtain:
ψ 0 ( x , t ) = A cos k x ε , ω 0 ( x , t ) = 0 .
ψ 1 ( x , t ) = t 2 ε Γ ( 1 + 2 ε ) A k 2 cos ε k x ε δ 2 1 Γ ( 1 + ε ) A cos ε k x ε + δ 2 1 Γ ( 1 + 3 ε ) A 3 cos ε 3 k x ε δ 2 1 Γ ( 1 + 5 ε ) A 5 cos ε 5 k x ε , ω 1 ( x , t ) = t 2 ε Γ ( 1 + 2 ε ) 1 Γ ( 1 + ε ) A cos ε k x ε 1 Γ ( 1 + 3 ε ) A 3 cos ε 3 k x ε + 1 Γ ( 1 + 5 ε ) A 5 cos ε 5 k x ε .
ψ 2 ( x , t ) = t 4 ε Γ ( 1 + 4 ε ) A k 2 + 2 k 2 δ 2 Γ ( 1 + ε ) + δ 4 + δ 2 Γ ( 1 + ε ) Γ ( 1 + ε ) + 6 A 2 k 2 δ 2 Γ ( 1 + 3 ε ) cos ε k x ε A 3 12 k 2 δ 2 Γ ( 1 + 3 ε ) + 20 A 2 k 2 δ 2 Γ ( 1 + 5 ε ) + 4 δ 4 + δ 2 Γ ( 1 + ε ) Γ ( 1 + 3 ε ) cos ε 3 k x ε + A 5 27 k 2 δ 2 Γ ( 1 + 5 ε ) + 3 δ 4 + 6 δ 2 Γ ( 1 + ε ) Γ ( 1 + 5 ε ) + 3 δ 4 + δ 2 Γ ( 1 + 3 ε ) Γ ( 1 + 3 ε ) cos ε 5 k x ε A 7 5 δ 4 + 8 δ 2 Γ ( 1 + 3 ε ) Γ ( 1 + 5 ε ) cos ε 7 k x ε + A 9 2 δ 4 + 5 δ 2 Γ ( 1 + 5 ε ) Γ ( 1 + 5 ε ) cos ε 9 k x ε + A 4 2 k 2 δ 2 Γ ( 1 + 5 ε ) + 2 δ 4 Γ ( 1 + ε ) Γ ( 1 + 5 ε ) cos ε 4 k x ε A 6 2 δ 4 Γ ( 1 + 3 ε ) Γ ( 1 + 5 ε ) cos ε 6 k x ε + A 8 2 δ 4 Γ ( 1 + 5 ε ) Γ ( 1 + 5 ε ) cos ε 8 k x ε ,
ω 2 ( x , t ) = t 4 ε Γ ( 1 + 4 ε ) A k 2 c 2 + 1 Γ ( 1 + ε ) + δ 2 + 1 Γ ( 1 + ε ) Γ ( 1 + ε ) + 6 A 2 k 2 c 2 Γ ( 1 + 3 ε ) cos ε k x ε + A 3 12 k 2 c 2 Γ ( 1 + 3 ε ) + 20 A 2 k 2 c 2 Γ ( 1 + 5 ε ) + 4 δ 2 + 1 Γ ( 1 + ε ) Γ ( 1 + 3 ε ) cos ε 3 k x ε A 5 27 k 2 c 2 Γ ( 1 + 5 ε ) + 3 δ 2 + 6 Γ ( 1 + ε ) Γ ( 1 + 5 ε ) + 3 δ 2 + 1 Γ ( 1 + 3 ε ) Γ ( 1 + 3 ε ) cos ε 5 k x ε + A 7 5 δ 2 + 8 Γ ( 1 + 3 ε ) Γ ( 1 + 5 ε ) cos ε 7 k x ε A 9 2 δ 2 + 5 Γ ( 1 + 5 ε ) Γ ( 1 + 5 ε ) cos ε 9 k x ε A 4 2 k 2 1 Γ ( 1 + 5 ε ) + 2 δ 2 Γ ( 1 + ε ) Γ ( 1 + 5 ε ) cos ε 4 k x ε + A 6 2 δ 2 Γ ( 1 + 3 ε ) Γ ( 1 + 5 ε ) cos ε 6 k x ε A 8 2 δ 2 Γ ( 1 + 5 ε ) Γ ( 1 + 5 ε ) cos ε 8 k x ε .
When ρ 1 , the solutions of the coupled local fractional sine-Gordon equations are
ψ ( x , t ) = ψ 0 ( x , t ) + ψ 1 ( x , t ) + ψ 2 ( x , t ) , ω ( x , t ) = ω 0 ( x , t ) + ω 1 ( x , t ) + ω 2 ( x , t ) .
When the fractal dimension is ε = l n 2 / l n 3 , the corresponding diagrams of ψ ( x , t ) and ω ( x , t ) are as shown in Figure 1 and Figure 2.
The time-dependent behavior of ψ ( x , t ) and ω ( x , t ) obtained from the theory of the local number order homotopy perturbation method is shown in Figure 3. We show the elastic interaction between two solitons. It can be seen that both solitons maintain their initial velocity and shape after interaction. In Figure 4, we compare the local motion diagram of the sine-Gordon equation with integer order and fractal dimension ε = l n 2 / l n 3 . Compared with the general homotopy perturbation method, the behavior graph obtained by this method is more realistic.

5. Conclusions

In this paper, we obtain the exact solution of a coupled local fractional sinusoidal Gordon equation by using the LFHP method, and show the special function graphs defined on the Cantor set when the fractal dimension is ε = l n 2 / l n 3 . The results show that the technique is effective in solving nonlinear partial order equations. The LFHP method solves many nonlinear problems in science and engineering. We consider that the LFHP method can theoretically solve the ac-driven sine-Gordon equation [29]. We will study it in future work. With the help of MATLAB and other mathematical software, this method provides a powerful mathematical tool for more complex nonlinear systems. In a word, the LFHP method provides highly accurate fractal solutions for nonlinear problems.

Author Contributions

Writing—original draft preparation and editing, L.C.; supervision, Q.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No.12062017), Program for Innovative Research Team in Universities of Inner Mongolia Autonomous Region, China (Grant No. NMGIRT2208) and the Natural Science Foundation of Inner Mongolia, China. (Grant No. 2022LHMS01001).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

In this part, the abbreviations used in our article are described below:
HPMhomotopy perturbation method
LFHPlocal fractional homotopy perturbation

References

  1. Charpentier, I. On Higher-order Differentiation in Nonlinear Mechanics. Optim. Methods Softw. 2012, 27, 221–232. [Google Scholar] [CrossRef]
  2. Viana, R.L.; Silva, E.; Kroetz, T.; Caldas, I.T.; Roberto, M.; Sanjuan, A.F. Fractal Structures in Nonlinear Plasma Physics. Philos. Trans. R. Soc. A 2011, 369, 371–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kalashnikov, A.S. The Propagation of Disturbances in Problem of Nonlinear Heat Conduction with Absorption. Comput. Math. Math. Phys. 1974, 14, 70–85. [Google Scholar] [CrossRef]
  4. Chowdhury, A.R.; Pal, S.K. N-fold Backlund Transformation for Deformed Nonlinear Schrödinger Equation. Int. J. Theor. Phys. 1997, 36, 1021–1031. [Google Scholar] [CrossRef]
  5. Liu, W.J.; Tian, B.; Zhang, H.Q.; Li, L.L.; Xue, Y.S. Soliton Interaction in the Higher-order Nonlinear Schrodinger equation Investigated with Hirota’s Bilinear Method. Phys. Rev. E 2008, 77, 066605. [Google Scholar] [CrossRef]
  6. Wazwaz, A.M. A Sine-cosine Method for Handling Nonlinear Wave Equations. Math. Comput. Model. 2004, 40, 499–508. [Google Scholar] [CrossRef]
  7. Ray, S.S.; Bera, R.K. An Approximate Solution of a Nonlinear Fractional Differential Equation by Adomian Decomposition Method. Appli. Math. Comp. 2005, 167, 561–571. [Google Scholar]
  8. Askari, H.; Zhang, D.; Esmailzadeh, E. Periodic Solutions for Nonlinear Oscillations of Nanowires Using Variational Iteration Method. In Proceedings of the 2013 13th IEEE International Conference on Nanotechnology (IEEE-NANO 2013), Beijing, China, 5–8 August 2013. [Google Scholar]
  9. Wu, X.H.; He, J.H. EXP-function Method and its Application to Nonlinear Equations. Chaos Solitons Fractals 2008, 38, 903–910. [Google Scholar]
  10. Rajabi, A.; Ganji, D.D.; Taherian, H. Application of Homotopy Perturbation Method in Nonlinear Heat Conduction and Convection Equations. Phys. Lett. 2007, 360, 570–573. [Google Scholar] [CrossRef]
  11. Asllanaj, F.; Jeandel, G.; Roche, J.R. Numerical Solution of Radiative Transfer Equation Coupled with Nonlinear Heat Conduction Equation. Int. J. Numer. Methods Heat Fluid Flow 2015, 11, 449–473. [Google Scholar] [CrossRef]
  12. He, J.H. Application of Homotopy Perturbation Method to Nonlinear Wave Equations. Chaos Solitons Fractals 2005, 26, 695–700. [Google Scholar] [CrossRef]
  13. Ganji, D.D.; Sadighi, A.; Ganjavi, B. Traveling Wave Solutions of the Sine-gordon and the Coupled Sine-gordon Equations Using the Homotopy-perturbation Method. Sci. Iran. 2009, 16, 189–195. [Google Scholar]
  14. Mustafa Inc. He’s Homotopy Perturbation Method for Solving Korteweg-de Vries Burgers Equation with Initial Condition. Numer. Methods Partial. Differ. Equ. 2010, 26, 1224–1235. [Google Scholar]
  15. Ganji, D.D. The Application of He’s Homotopy Perturbation Method to Nonlinear Equations Arising in Heat Transfer-Science Direct. Phys. Lett. A 2006, 355, 337–341. [Google Scholar] [CrossRef]
  16. Ganji, D.D.; Sadighi, A. Application of He’s Homotopy-perturbation Method to Nonlinear Coupled Systems of Reaction-diffusion Equations. Int. J. Nonlinear Sci. Num. Simul. 2006, 7, 411–418. [Google Scholar] [CrossRef]
  17. Yang, X.J. Advanced Local Fractional Calculus and Its Applications; Orld Science: New York, NY, USA, 2012. [Google Scholar]
  18. Yang, X.J.; Srivastava, H.M.; Cattani, C. Local Fractional Homotopy Perturbation Method for Solving Fractal Partial Differential Equation Arising in Mathematical Physics. Rom. Rep. Phys. 2015, 67, 752–761. [Google Scholar]
  19. Zhang, Y.; Carlo, C. and Yang, X.J. Local Fractional Homotopy Perturbation Method for Solving Non-Homogeneous Heat Conduction Equations in Fractal Domains. Entropy 2015, 17, 6753–6764. [Google Scholar] [CrossRef] [Green Version]
  20. Braun, O.M.; Kivshar, Y.S. Nonlinear Dynamics of the Frenkel-Kontorova Model with Impurities. Phys. Rev. B Condens Matter 1998, 306, 1060–1073. [Google Scholar] [CrossRef]
  21. Saha, S.; Ray. A Numerical Solution of the Coupled Sine-Gordon Equation Using the Modified Decomposition Method. Appl. Math. Comput. 2006, 175, 1046–1054. [Google Scholar]
  22. Hosseini, K.; Mayeli, P.; Kumar, D. New exact solutions of the coupled sine-Gordon equations in nonlinear optics using the modified Kudryashov method. J. Mod. Opt. 2017, 65, 361–364. [Google Scholar] [CrossRef]
  23. Salas, A.H. Exact Solutions of Coupled Sine-Gordon Equations. Nonlinear Anal.-Real 2010, 11, 3930–3935. [Google Scholar] [CrossRef]
  24. Zhao, Y.M.; Liu, H.H.; Yang, Y.J. Exact Solutions for the Coupled Sine-Gordon Equations by a New Hyperbolic Auxiliary Function Method. Appl. Math. Sci. 2011, 5, 1621–1629. [Google Scholar]
  25. Yang, X.J.; Machado, J.A.T.; Hristov, J. Nonlinear Dynamics for Local Fractional Burgers’ equation Arising in Fractal Flow. Nonlinear Dyn. 2016, 84, 3–7. [Google Scholar] [CrossRef] [Green Version]
  26. Ghanbari, B. On Novel Nondifferentiable Exact Solutions to Local Fractional Gardner’s Equation Using an Effective Technique. Math. Method Appl. Sci. 2021, 44, 4673–4685. [Google Scholar] [CrossRef]
  27. Yang, X.J.; Machado, J.T.; Baleanu, D.; Cattani, C. On Exact Traveling-wave Solutions for Local Fractional Korteweg-de Vries Equation. Chaos 2016, 26, 110–118. [Google Scholar] [CrossRef]
  28. Yang, X.J. Local Fractional Integral Transforms and Their Applications; Academic Press: Cambridge, MA, USA, 1995. [Google Scholar]
  29. Quintero, R.N.; Sánchez, A. Ac driven sine-Gordon solitons: Dynamics and stability. Eur. Phys. J. B 1998, 6, 133–142. [Google Scholar] [CrossRef] [Green Version]
Figure 1. In the case of ψ ( x , t ) solution in the local fractional sine-Gordon equations, by choosing A = 1 , k = 1 , c = 1 , δ = 1 . (a) 0 ≤ t ≤ 1; (b) −1 ≤ t ≤ 1.
Figure 1. In the case of ψ ( x , t ) solution in the local fractional sine-Gordon equations, by choosing A = 1 , k = 1 , c = 1 , δ = 1 . (a) 0 ≤ t ≤ 1; (b) −1 ≤ t ≤ 1.
Fractalfract 06 00404 g001
Figure 2. In the case of ω ( x , t ) solution in the local fractional sine-Gordon equations, by choosing A = 1 , k = 1 , c = 1 , δ = 1 . (a) 0 ≤ t ≤ 1; (b) −1 ≤ t ≤ 1.
Figure 2. In the case of ω ( x , t ) solution in the local fractional sine-Gordon equations, by choosing A = 1 , k = 1 , c = 1 , δ = 1 . (a) 0 ≤ t ≤ 1; (b) −1 ≤ t ≤ 1.
Fractalfract 06 00404 g002
Figure 3. In the case of local fractional sine-Gordon equations, by choosing parameters A = 1 , k = 1 , c = 1 , δ = 1 .
Figure 3. In the case of local fractional sine-Gordon equations, by choosing parameters A = 1 , k = 1 , c = 1 , δ = 1 .
Fractalfract 06 00404 g003
Figure 4. Comparison of integer order ε = 1 and fractional dimension ε = l n 2 / l n 3 .
Figure 4. Comparison of integer order ε = 1 and fractional dimension ε = l n 2 / l n 3 .
Fractalfract 06 00404 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Liu, Q. Local Fractional Homotopy Perturbation Method for Solving Coupled Sine-Gordon Equations in Fractal Domain. Fractal Fract. 2022, 6, 404. https://doi.org/10.3390/fractalfract6080404

AMA Style

Chen L, Liu Q. Local Fractional Homotopy Perturbation Method for Solving Coupled Sine-Gordon Equations in Fractal Domain. Fractal and Fractional. 2022; 6(8):404. https://doi.org/10.3390/fractalfract6080404

Chicago/Turabian Style

Chen, Liguo, and Quansheng Liu. 2022. "Local Fractional Homotopy Perturbation Method for Solving Coupled Sine-Gordon Equations in Fractal Domain" Fractal and Fractional 6, no. 8: 404. https://doi.org/10.3390/fractalfract6080404

APA Style

Chen, L., & Liu, Q. (2022). Local Fractional Homotopy Perturbation Method for Solving Coupled Sine-Gordon Equations in Fractal Domain. Fractal and Fractional, 6(8), 404. https://doi.org/10.3390/fractalfract6080404

Article Metrics

Back to TopTop