Next Article in Journal
Adaptive Control Design for Euler–Lagrange Systems Using Fixed-Time Fractional Integral Sliding Mode Scheme
Previous Article in Journal
A New Approach to Determining the Network Fractality with Application to Robot-Laser-Hardened Surfaces of Materials
Previous Article in Special Issue
Fractional-Order Total Variation Geiger-Mode Avalanche Photodiode Lidar Range-Image Denoising Algorithm Based on Spatial Kernel Function and Range Kernel Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy-Preserving AVF Methods for Riesz Space-Fractional Nonlinear KGZ and KGS Equations

School of Mathematics and Statistics, Hainan University, Haikou 570228, China
*
Author to whom correspondence should be addressed.
Submission received: 27 June 2023 / Revised: 25 August 2023 / Accepted: 3 September 2023 / Published: 27 September 2023

Abstract

:
The Riesz space-fractional derivative is discretized by the Fourier pseudo-spectral (FPS) method. The Riesz space-fractional nonlinear Klein–Gordon–Zakharov (KGZ) and Klein–Gordon–Schrödinger (KGS) equations are transformed into two infinite-dimensional Hamiltonian systems, which are discretized by the FPS method. Two finite-dimensional Hamiltonian systems are thus obtained and solved by the second-order average vector field (AVF) method. The energy conservation property of these new discrete schemes of the fractional KGZ and KGS equations is proven. These schemes are applied to simulate the evolution of two fractional differential equations. Numerical results show that these schemes can simulate the evolution of these fractional differential equations well and maintain the energy-preserving property.

1. Introduction

Fractional differential equations can better describe the behavior of physical phenomena than integral differential equations. Many scholars have taken great interest in studying fractional differential equations and the theory of the fractional derivative. In general, factional differential equations can not have an exact solution. Numerical simulations for fractional nonlinear differential equations have become very important. Many different numerical methods have been to proposed to solve fractional nonlinear partial differential equations (PDEs) [1,2,3,4]. In this paper, we numerically investigate the following Riesz space fractional KGZ and KGS equations by the energy preserving method.
The space fractional KGZ equation can be written as [5,6,7,8]
2 u ( x , t ) t 2 α u ( x , t ) | x | α + u ( x , t ) + m ( x , t ) u ( x , t ) + | u ( x , t ) | 2 u ( x , t ) = 0 , 2 m ( x , t ) t 2 2 m ( x , t ) x 2 2 ( | u ( x , t ) | 2 ) x 2 = 0 ,
which describes the propagation of Langmuir waves in plasma physics. Suppose the finite domain space Ω = ( a , b ) × ( 0 , T ) with the initial conditions u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 x , m ( x , 0 ) = m 0 ( x ) , m t ( x , 0 ) = m 1 ( x ) , x [ a , b ] and the boundary conditions u ( a , t ) = u ( b , t ) = 0 , m ( a , t ) = m ( b , t ) = 0 , t [ 0 , T ] , where α u ( x , t ) | x | α is the Riesz space fractional derivative with 1 α 2 .
The KGZ equation can have the following invariant energy conservation
E ( t ) = [ | u t | 2 + | α 2 u | x | α 2 | 2 + | u | 2 + m | u | 2 + 1 2 v 2 + 1 2 | u | 4 + 1 2 m 2 ] d x = E ( 0 ) ,
where v ( x , t ) x = m ( x , t ) t . When α = 2 , the Riesz space fractional KGZ equation becomes the classical KGZ equation [9,10,11]. Ma et al. proposed a multi-scale integrator method to solve the KGZ equation [12]. Wang et al. solved the KGZ equation based on the lattice Boltzmann model [13]. Mei et al. proposed the Galerkin finite element method to solve the KGZ equation [14]. A lot of work has also been conducted on the Riesz space fractional KGZ equation [8,15]. Macias et al. proposed an energy conserving scheme to solve the fractional KGZ equation [15,16,17,18,19].
The fractional KGS equation can be written in the following form:
i f t β 2 ( Δ ) α 2 f + α 1 f h = 0 , h t t γ h x x + m 2 h α 1 | f | 2 = 0 ,
where i = 1 . The initial conditions are f ( x , 0 ) = f 0 ( x ) , h t ( x , 0 ) = h 1 ( x ) , x R , and the boundary conditions are f ( x , t ) = h ( x , t ) = 0 , x R Ω , Ω = ( x L , x R ) , t [ 0 , T ] . The complex function f ( x , t ) is the complex neutron field and the real function h ( x , t ) is the meson field. The variables α 1 and β are the coupled constants.
The fractional KGS equation also has invariant energy conservation:
E ˜ ( t ) = + h t 2 + γ h x 2 + m 2 h 2 + β | ( Δ ) α 4 f | 2 2 α 1 | f | 2 h d x = E ˜ ( 0 ) .
When α = 2 and β = γ = 1 , Equation (3) is the classical coupled KGS equation. In quantum field theory, the coupled KGS equation is a mathematical model for the interaction of a conservation complex neutron field and a real meson field. Regarding the coupled KGS equation, many important conclusions have been obtained. Ohta et al. analyzed the stability of the coupled KGS equation [20]. Guo et al. investigated the global well-posedness of the KGS equation [21]. Kong et al. proposed a symplectic method to solve the coupled KGS equation [22]. The fractional KGS equation is an extension of the classical coupled KGS equation. Wang et al. proposed to solve the Riesz space fractional KGS equation using the difference method and the spectral method [23,24].
Recently, energy preserving methods have become important numerical methods in simulating energy conservation nonlinear PDEs, and they are structure-preserving numerical methods. Many different energy preserving methods have been derived, such as the discrete gradient method [25,26], the discrete variational derivative method [27] and the Hamiltonian boundary value method [28]. The AVF method, which is a kind of discrete gradient method, has been widely used to solve energy conservation integral PDEs and has achieved great success [29,30]. However, few people have applied the AVF method to solve energy conservation fractional PDEs. In this paper, we apply the AVF method to solve Riesz space fractional KGZ and KGS equations.
The rest of the paper is organized as follows. In Section 2, the definition and properties of the Riesz fractional derivative are given. The Riesz fractional derivative is discretized by the FPS method. In Section 3, the infinite Hamiltonian symplectic structures of the Riesz space fractional KGZ and KGS equations are obtained. These Hamiltonian systems are discretized by the FPS method and the second-order AVF method. Two energy preserving schemes of the Riesz space fractional KGZ and KGS equations are thus obtained. At last, numerical experiments confirm the advantage of the energy preserving schemes of the Riesz space fractional KGZ and KGS equations in simulating solitary wave behavior and preserving the energy conservation property of the equations. The new energy preserving schemes are superior to the existing second-order energy preserving schemes.

2. Discretization of the Riesz Space Fractional Derivative

Definition 1. 
When n 1 α n , n is a positive integer, the Riesz space fractional derivative with α order can be defined as [2]
α u ( x , t ) | x | α = 1 2 c o s ( π α 2 ) Γ ( n α ) n x n + u ( ξ , t ) | x ξ | α + 1 n d ξ ,
where Γ ( . ) is the Gamma function. | x | α is denoted as the α order derivative of u ( x , t ) at ( x , t ) .
Lemma 1. 
In the infinite domain interval ( < x < ) , the function u ( x , t ) is equivalent to
α u ( x , t ) | x | α = 1 2 c o s ( π α 2 ) Γ ( n α ) n x n + u ( ξ , t ) | x ξ | α + 1 n d ξ = ( ) α 2 u ( x , t ) ,
where n 1 < α < n .
In the infinite domain interval ( < x < ) , the fractional Laplace operator is defined as
( ) α 2 u ( x , t ) = F 1 | x | α F u ( x , t ) ,
where F and F 1 are represented as the Fourier transformation and the Fourier inverse transform of u ( x , t ) , and we can obtain
( ) α 2 u ( x , t ) = 1 2 π + e i x ξ | x | α + e i ξ η u ( η , t ) d η d ξ .
On the other hand, in the finite domain interval Ω = ( a , b ) , the Fourier series can be defined as
( ) α 2 u ( x , t ) = l Z | v l | α u ^ l e i v l ( x a ) ,
where ν l = 2 l π b a , and the coefficients of the Fourier series are
u ^ l = 1 b a Ω u ( x , t ) e i v l ( x a ) d x .
We apply the FPS method to discrete the α order Riesz space fractional derivative. Suppose the space integral interval is Ω = [ a , b ] , the interval Ω is divided into N equal parts. N is an even number. The space step length is h 1 = b a N . Take x j = a + j h 1 , where j = 0 , 1 , , N 1 are the space Fourier collocation points. I N u ( x , t ) is the approximation of u ( x , t ) , we have
( I N u ) ( x , t ) = u N ( x , t ) = k = N / 2 N / 2 u ˜ k e i k μ ( x a ) ,
where u ˜ k , μ = 2 π b a and
u ˜ k = 1 N c k j = 0 N 1 u ( x j , t ) e i k μ ( x j a ) , μ = 2 π b a , k < N / 2 , c k = 1 , k = ± N / 2 c k = 2 .
When | k | < N / 2 , c k = 1 , and when k = ± N / 2 , c k = 2 .
We can obtain
( ) α 2 I N u ( x j , t ) = k = N / 2 N / 2 k μ α u ˜ k e i k μ ( x j a ) .
The α order derivative of the approximation function I N u ( x , t ) can be denoted as
α I N u ( x j , t ) | x | α = ( ) α 2 I N u ( x j , t ) = k = N / 2 N / 2 | k μ | α ( 1 N c k l = 0 N 1 u l e i k μ ( x l a ) ) e i k μ ( x j a ) = l = 0 N 1 u l ( k = N / 2 N / 2 1 N c k | k μ | α e i k μ ( x j x l ) = ( D 2 α U ) j ,
where D 2 α is N × N matrix, U = ( U 0 , U 1 , , U N 1 ) and the coefficients of the matrix D 2 α are
( D 2 α ) j , l = k = N / 2 N / 2 1 N c k | k μ | α e i k μ ( x j x l ) .

3. Energy Preserving Method of Fractional KGZ and KGS Equations

In this section, we first give the Hamiltonian structures of the Riesz space fractional KGZ and KGS equations. Then, the space fractional derivatives of these two equations are discretized by the FPS method. The AVF method is applied to solve the semi-discrete fractional KGZ and KGS equations.

3.1. Energy Preserving Method of the Fractional KGZ Equation

Let w ( x , t ) = u t ( x , t ) , 2 q x x ( x , t ) = m t ( x , t ) , then the space fractional KGZ Equation (1) is equivalent to
u t = w , m t = 2 q x x , w t = α u | x | α u m u | u | 2 u , q t = 1 2 m 1 2 | u | 2 .
Equation (16) can be expressed by the following infinite dimensional Hamiltonian system
d z d t = J δ H ( z ) δ z , J = O I I O , I = 1 0 0 1 ,
where z = ( u , m , w , q ) T , O is 2 × 2 zero matrix, I is 2 × 2 unite matrix and the corresponding Hamiltonian energy function is
H ( z ) = [ 1 2 w 2 + ( q x ) 2 + 1 2 | α 2 u | x | α 2 | 2 + 1 2 | u | 2 + 1 2 m | u | 2 + 1 4 | u | 4 + 1 4 m 2 ] d x .
The second-order partial derivative q x x of Equation (16) can be approximated by the FPS method [31,32]. Suppose I N q ( x , t ) is the FPS approximation of the function q ( x , t ) . The values of the derivatives d d x I N q ( x , t ) and d 2 d x 2 I N p ( x , t ) at the collocation points x j are obtained in terms of the value of q j , i.e.,
d d x I N q ( x , t ) | x = x j = l = 0 N 1 q l d g l ( x j ) d x = ( D 1 Q ) j ,
d 2 d x 2 I N q ( x , t ) | x = x j = l = 0 N 1 q l d 2 g l ( x j ) d x 2 = ( D 2 Q ) j ,
The function g l ( x ) is the trigonometric polynomial explicitly given by
g l ( x ) = 1 N k = N / 2 N / 2 1 c k exp i k μ ( x x l ) ,
where c k = 1 ( | k | N / 2 ) , c N / 2 = c N / 2 = 2 , μ = 2 π L and
( D 2 ) i , j = 1 2 μ 2 ( 1 ) i + j + 1 1 sin 2 ( μ x i x j 2 ) , i j , μ 2 N 2 + 2 12 , i = j .
Based on the above FPS method, we can obtain the semi-discrete system of the space fractional KGZ equation
d d t u j = w j , d d t m j = 2 ( D 2 Q ) j , d d t w j = ( D 2 α U ) j u j m j u j | u j | 2 u j , d d t q j = 1 2 m j 1 2 | u j | 2 ,
where j = 0 , 1 , , N 1 .
Equation (23) can be written as the following semi-discrete Hamiltonian system
d Z d t = J ^ Z H ( Z ) , J ^ = 0 0 I N 0 0 0 0 I N I N 0 0 0 0 I N 0 0 ,
where Z = ( U T , M T , W T , Q T ) , M = ( m 0 , m 1 , , m N 1 ) T , W = ( w 0 , w 1 , , w N 1 ) T , Q = ( q 0 , q 1 , , q N 1 ) T . 0 and I N are the N × N zeros matrix and unite matrix. The corresponding Hamiltonian function is
H ( Z ) = Q T D 2 Q 1 2 U T D 2 α U + j = 0 N 1 [ 1 2 w j 2 + 1 2 | u j | 2 + 1 2 m j | u j | 2 + 1 4 | u j | 4 + 1 4 m j 2 ] .
The semi-discrete Hamiltonian system (24) is solved by the following second-order AVF method
Z n + 1 Z n τ = J ^ 0 1 H ( ( 1 ξ ) Z n + ξ Z n + 1 ) d ξ .
Equation (26) is equivalent to the following equations
u j n + 1 u j n τ = 0 1 ( ( 1 ξ ) w j n + ξ w j n + 1 ) d ξ ,
m j n + 1 m j n τ = 2 0 1 ( ( 1 ξ ) ( D 2 Q n ) j + ξ ( D 2 Q n + 1 ) j ) d ξ ,
w j n + 1 w j n τ = 0 1 [ ( ( 1 ξ ) D 2 α U n j + ξ D 2 α U n + 1 j ) ( ( 1 ξ ) u j n + ξ u j n + 1 ) ] d ξ 0 1 ( 1 ξ ) m j n + ξ m j n + 1 ( 1 ξ ) u j n + ξ u j n + 1 d ξ 0 1 ( 1 ξ ) u j n + ξ u j n + 1 2 ( 1 ξ ) u j n + ξ u j n + 1 d ξ ,
q j n + 1 q j n τ = 1 2 0 1 ( ( ( 1 ξ ) m j n + ξ m j n + 1 ) + ( ( 1 ξ ) u j n + ξ u j n + 1 ) 2 ) d ξ .
From Equations (27)–(30), the auxiliary variables w, q can be deleted. We can obtain
u j n + 1 2 u j n + u j n 1 τ 2 = ( D 2 α U n + 1 + U n 4 ) j u j n + 1 + u j n 4 m j n + 1 + m j n 4 u j n 2 m j n + 1 + m j n 12 ( u j n + 1 u j n ) 1 6 ( u j n + 1 ) 2 + 1 6 ( u j n ) 2 + 1 6 u j n + 1 u j n u j n 1 8 ( u j n + 1 ) 2 + 1 24 ( u j n ) 2 + 1 12 u j n + 1 u j n ( u j n + 1 u j n ) + ( D 2 α U n + U n 1 4 ) j u j n + u j n 1 4 m j n + m j n 1 4 u j n 1 2 m j n + m j n 1 12 ( u j n u j n 1 ) 1 6 ( u j n ) 2 + 1 6 ( u j n 1 ) 2 + 1 6 u j n u j n 1 u j n 1 1 8 ( u j n ) 2 + 1 24 ( u j n 1 ) 2 + 1 12 u j n u j n 1 ( u j n u j n 1 ) ,
m j n + 1 2 m j n + m j n 1 τ 2 = ( D 2 M n + 1 + 2 M n + M n 1 4 ) j + l = 1 N d j , l ( 1 6 ( u j n + 1 ) 2 + 1 6 ( u j n ) 2 + 1 6 u j n + 1 u j n + 1 6 ( u j n ) 2 + 1 6 ( u j n 1 ) 2 + 1 6 u j n u j n 1 ) .
Theorem 1. 
The semi-discrete scheme (26) can preserve the discrete energy conservation of the finite dimensional Hamiltonian system
H ( Z n + 1 ) = H ( Z n ) .
Proof. 
From Equation (26), we can obtain that
( 0 1 H ( ( 1 ξ ) Z n + ξ Z n + 1 ) d ξ ) T Z n + 1 Z n τ = ( 0 1 H ( ( 1 ξ ) Z n + ξ Z n + 1 ) d ξ ) T J 0 1 H ( ( 1 ξ ) Z n + ξ Z n + 1 ) d ξ ,
where J is the skew symmetric matrix. We can find that
1 τ 0 1 ( Z n Z n + 1 ) T H ( ( 1 ξ ) Z n + ξ Z n + 1 ) d ξ = 0 .
From the basic theorem of the integral theory, we can get
1 τ ( H ( Z n + 1 ) H ( Z n ) ) = 0 .
The proof of the theorem is completed. □

3.2. Energy Preserving Method of the Fractional KGS Equation

Suppose f ( x , t ) = p ( x , t ) + r ( x , t ) i , p ( x , t ) = p , r ( x , t ) = r , then Equation (3) can be written as the following differential equations
r t = β 2 ( ) α 2 p + α 1 h p , p t = β 2 ( ) α 2 r α 1 h r , v t = γ 2 h x x m 2 2 h + α 1 2 ( p 2 + r 2 ) , h t = 2 v .
Equation (37) can be transformed into the following infinite dimensional Hamiltonian system
d z ˜ d t = J ˜ δ H ( z ˜ ) δ z ˜ ,
where z ˜ = ( r , p , v , h ) T ,
J ˜ = 0 1 0 0 1 0 0 0 0 0 0 1 0 0 1 0 ,
and the corresponding Hamiltonian energy function is
H ( z ˜ ) = β 4 ( ( ) α 4 p ) 2 ( ) α 4 r ) 2 ) γ 4 ( h x ) 2 1 4 m 2 h 2 + α 1 2 h ( p 2 + r 2 ) v 2 d x .
Equation (37) is discretized by the FPS method, and we can obtain
d d t r j = β 2 ( D 2 α p ) j + α 1 h j p j , d d t p j = β 2 ( D 2 α r ) j α 1 h j r j , d d t v j = γ 2 ( D 2 h ) j m 2 2 h j + α 1 2 ( p j 2 + r j 2 ) , d d t h j = 2 v j .
Equation (41) can be transformed into the finite dimensional Hamiltonian system
d Z ˜ d t = S z H ( Z ˜ ) ,
where Z ˜ = ( r 0 , , r N 1 , p 0 , , p N 1 , v 0 , , v N 1 , h 0 , , h N 1 ) T and
S = 0 I N 0 0 I N 0 0 0 0 0 0 I N 0 0 I N 0 .
H ( Z ˜ ) = β 4 ( p T D 2 α p + r T D 2 α r ) + γ 4 ( h T D 2 h + j = 0 N 1 [ α 1 2 h j ( p j 2 + r j 2 ) 1 4 m 2 h j 2 v j 2 ] .
The semi-discrete Hamiltonian system (42) is solved by the second order AVF method, and we can get
Z ˜ n + 1 Z ˜ n τ = S 0 1 Z ˜ H ( ( 1 ξ ) Z ˜ n + ξ Z ˜ n + 1 ) d ξ .
Equation (45) is equivalent to the following schemes
r j n + 1 r j n τ = β 2 ( D 2 α p n + 1 + p n 2 ) j + α 1 ( 1 3 h j n + 1 p j n + 1 + 1 6 h j n + 1 p j n + 1 6 h j n p j n + 1 + 1 3 h j n p j n ) ,
p j n + 1 p j n τ = β 2 ( D 2 α r n + 1 + r n 2 ) j α 1 ( 1 3 h j n + 1 r j n + 1 + 1 6 h j n + 1 r j n + 1 6 h j n r j n + 1 + 1 3 h j n r j n ) ,
v j n + 1 v j n τ = γ 2 ( D 2 h n + 1 + h n 2 ) j m 2 2 ( h j n + 1 + h j n 2 ) + α 1 6 ( p j n + 1 p j n + 1 + p j n + 1 p j n + p j n p j n + r j n + 1 r j n + 1 + r j n + 1 r j n + r j n r j n ) ,
h j n + 1 h j n τ = 2 v j n + 1 + v j n 2 .
According to the above Theorem 1, schemes (46)–(49) can also preserve the discrete Hamiltonian energy of the Riesz space fractional KGS equation exactly.

4. Numerical Simulation

4.1. Numerical Simulation of the Fractional KGZ Equation

In this section, we test the following discrete Hamiltonian energy errors
E r r o r H n = | H ( Z n ) H ( Z 0 ) H ( Z 0 ) | ,
by schemes (31) and (32), where
H ( Z n ) = ( Q n ) T D 2 Q n 1 2 ( U n ) T D 2 α U n + j = 0 N 1 [ 1 2 ( w j n ) 2 + 1 2 | u j n | 2 + 1 2 m j n | u j 2 | 2 + 1 4 | u j 2 | 4 + 1 4 ( m j n ) 2 ] ,
and H ( Z 0 ) is the discrete Hamiltonian energy at t = 0 .
The two level initial conditions can be taken as follows [16]:
u 0 = 10 2 2 sech ( 1 + 5 2 x ) exp ( i 2 1 + 5 x ) , m 0 = 2 sech 2 ( 1 + 5 2 x ) ,
and
u 1 = 10 2 2 ( tanh x 1 ) sech ( 1 + 5 2 x ) exp ( i 2 1 + 5 x ) , m 1 = 4 sech 2 ( 1 + 5 2 x ) tanh ( 1 + 5 2 x ) .
Figure 1 gives the numerical solution of the KGZ equation with α = 2 , t [ 0 , 10 ] . Figure 2 gives the numerical solution of the Riesz space fractional KGZ equation with α = 1.7 , t [ 0 , 16 ] . These numerical results are consistent with the existing numerical results [16]. From Figure 1 and Figure 2, we can see that the new scheme can simulate the evolution of solitary waves of the Riesz space fractional KGZ equation well. The numerical solutions and the numerical schemes are stable. Figure 3 shows the energy error of the Riesz space fractional KGZ equation with ( a ) : α = 2 , ( b ) : α = 1.7 at t [ 0 , 20 ] . The energy error is up to 10 13 , which can neglected. It is obvious that the new scheme can well preserve the discrete energy of the fractional KGZ equation.

4.2. Numerical Simulation of the Fractional KGS Equation

We apply schemes (46)–(49) to simulate the Riesz space fractional KGS equation. The discrete Hamiltonian energy errors can be defined as
E r r o r H n = | H ( Z ˜ n ) H ( Z ˜ 0 ) H ( Z ˜ 0 ) | ,
where
H ( Z ˜ n ) = β 4 ( ( p n ) T D 2 α p n + ( r n ) T D 2 α r n ) + γ 4 ( h n ) T D 2 h n + j = 0 N 1 [ α 1 2 h j n ( ( p j n ) 2 + ( r j n ) 2 ) 1 4 m 2 ( h j n ) 2 ( v j n ) 2 ]
and H ( Z ˜ 0 ) is the discrete Hamiltonian energy at t = 0 .
The parameter variables are taken as υ = 0.8 , x 0 = 10 . First, we consider the evolution of single solitary waves. The initial condition of single solitary wave is taken as follows
f 0 = f ( x x 0 , 0 , v ) = 3 2 4 ( 1 υ 2 ) sech 2 ( 1 2 1 υ 2 ( x x 0 ) exp ( i υ x ) , h 0 = h ( x x 0 , 0 , v ) = 3 4 ( 1 υ 2 ) sech 2 ( 1 2 1 υ 2 ( x x 0 ) .
Figure 4 shows the numerical solution of a single solitary wave of the KGS equation at α = 2 , t [ 0 , 20 ] . Figure 5 shows the numerical solution of solitary waves f ( x , t ) and h ( x , t ) of the Riesz space fractional KGS equation at α = 1.2 , t [ 0 , 20 ] . From Figure 4 and Figure 5, we can see that the new scheme of the Riesz space fractional KGS equation can also simulate the evolution of a solitary wave of the equation well. The numerical solutions and the numerical schemes are stable. Figure 6 shows the energy error of the Riesz space fractional KGS equation. The error can be neglected. The new scheme of the equation can also well preserve the discrete energy of the equation.
Then, we consider the evolution of two solitary waves with the following initial conditions
f 0 = f ( x x 0 , 0 , v ) + f ( x + x 0 , 0 , v ) , h 0 = h ( x x 0 , 0 , v ) + h ( x + x 0 , 0 , v ) .
Figure 7 shows the numerical solutions of two solitary waves of the KGS equation at α = 2 , t [ 0 , 30 ] . Figure 8 shows numerical solutions of two solitary waves of the Riesz space fractional KGS equation at α = 1.6 , t [ 0 , 30 ] . The numerical results are consistent with the existing results [23]. From Figure 7 and Figure 8, we can see that the new scheme of the Riesz space fractional KGS equation can also simulate the evolution of multiple solitary waves of the equation well. The numerical solutions and the numerical schemes are stable. Figure 9 shows the energy error of the Riesz space fractional KGS equation. The error can be neglected. The new scheme of the Riesz space fractional KGS equation can also preserve the discrete energy of the equation well.

5. Conclusions

In this paper, two new energy preserving schemes for the Riesz space fractional KGZ and KGS equations are proposed based on the FPS method and the AVF method. The energy conservation property is proven. Numerical results show that these new schemes can simulate the evolution behavior of the solitary waves of these fractional differential equations well and preserve the discrete energy conservation property. The existing energy conserving schemes of the space fractional KGZ and KGS equations are in general of second-order accuracy. However, these energy conservation schemes based on the second-order AVF method can be easily extended to high-order energy preserving schemes. In future work, we will construct high-order energy preserving schemes for fractional differential equations based on the high-order AVF method and analyze the convergence and stability of these new energy-preserving schemes.

Author Contributions

Conceptualization, J.S., S.Y. and L.Z.; methodology, J.S., S.Y. and L.Z.; data collection, S.Y. and L.Z.; writing—original draft and preparation, J.S., S.Y. and L.Z.; data analysis, L.Z. and S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Hainan Provincial Natural Science Foundation of China (Grant No. 120RC450) and the Natural Science Foundation of China (Grant Nos. 11961020, 41974114).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
KGSKlein–Gordon–Schrödinger
KGZKlein–Gordon–Zakharov
AVFaverage vector field
PDEspartial differential equations
FPSFourier pseudo-spectral

References

  1. Ray, S.S. A new analytical modelling for nonlocal generalized Riesz fractional sine-Gordon equation. J. King Saud Univ. Sci. 2016, 28, 48–54. [Google Scholar]
  2. Wang, P.D.; Huang, C.M. Structure-preserving numerical methods for the fractional Schrödinger equation. Appl. Numer. Math. 2018, 129, 137–158. [Google Scholar] [CrossRef]
  3. Sun, Z.; Gao, G. Finite Difference Methods for Fractional Differential Equations; Science Press: Beijing, China, 2015. [Google Scholar]
  4. Li, C.P.; Zeng, F.H. Numerical Methods for Fractional Calculus; CRC Press: New York, NY, USA, 2015; pp. 1–17. [Google Scholar]
  5. Ding, H.; Li, C. Fractional-compact numerical algorithms for Riesz spatial fractional reaction-dispersion equations. Fract. Calc. Appl. Anal. 2017, 20, 722–764. [Google Scholar] [CrossRef]
  6. Denghan, M.; Nikpour, A. The solitary wave solution of coupled Klein-Gordon-Zakharov equations via two different numerical methods. Comput. Phys. Commun. 2013, 184, 2145–2158. [Google Scholar] [CrossRef]
  7. Wang, J.J. Solitary wave propagation and interactions for the Klein-Gordon-Zakharov equations in plasma physics. J. Phys. A Math. Theory 2009, 42, 085205. [Google Scholar] [CrossRef]
  8. Macias-Diaz, J.E.; Hemdy, A.S.; Staelen, R.D. A compact fourth-order in space energy-preserving method for Riesz space-fractional nonlinear wave equations. Appl. Math. Comput. 2018, 325, 1–14. [Google Scholar] [CrossRef]
  9. Su, C.M.; Zhao, X.F. A uniformly first-order accurate method for Klein-Gordon-Zakharov system in simultaneous high-plasma-frequency and subsonic limit regime. J. Comput. Phys. 2021, 428, 110064. [Google Scholar] [CrossRef]
  10. Bao, W.; Su, C. Uniform error bounds of a finite difference method for the Klein-Gordon-Zakharov system in the subsonic limit regime. Math. Comput. 2018, 87, 2133–2158. [Google Scholar] [CrossRef]
  11. Akbulut, A.; Tascn, F. Application of conservation theorem and modified extended tanh-function method to (1+1) dimensional nonlinear coupled Klein-Gordon-Zakharov equation. Chaos Solitons Fractals 2017, 104, 33–40. [Google Scholar] [CrossRef]
  12. Ma, Y.; Su, C.M. A uniformly and optimally accurate multiscale time integrator method for the Klein-Gordon-Zakharov system in the subsonic limit regime. Comput. Math. Appl. 2018, 76, 602–619. [Google Scholar] [CrossRef]
  13. Wang, H.M. Numerical simulation for solitary wave of Klein-Gordon-Zakharov equation based on the lattice Boltzmann model. Comput. Math. Appl. 2019, 78, 3941–3955. [Google Scholar] [CrossRef]
  14. Gao, Y.L.; Mei, L.Q.; Li, R. Galerkin finite element methods for the generalized Klein-Gordon-Zakharov equations. Comput. Math. Appl. 2017, 74, 2466–2484. [Google Scholar] [CrossRef]
  15. Macias-Diaz, J.E.; Hendy, A.S.; Staelen, R.D. A pseudo energy-invariant method for relativistic wave equations with Riesz space fractional derivatives. Comput. Phys. Commun. 2017, 224, 98–107. [Google Scholar] [CrossRef]
  16. Hendy, A.S.; Macias-Diaz, J.E. A numerically efficient and conservative model for a Riesz space-fractional Klein-Gordon-Zakharov system. Commun. Nonlinear Sci. Numer. Simul. 2019, 71, 22–37. [Google Scholar] [CrossRef]
  17. Martinez, R.; Macias-Diaz, J.E.; Hendy, A.S. Theoretical analysis of an explicit energy-conserving scheme for a fractional Klein-Gordon-Zakharov system. Appl. Numer. Math. 2019, 146, 245–259. [Google Scholar] [CrossRef]
  18. Martinez, R.; Macias-Diaz, J.E. An energy-preserving and efficient scheme for a double-fractional conservative Klein-Gordon-Zakharov system. Appl. Numer. Math. 2020, 158, 292–313. [Google Scholar] [CrossRef]
  19. Macias-Diaz, J.E. Existence of solutions of an explicit energy-conserving scheme for a fractional Klein-Gordon-Zakharov system. Appl. Numer. Math. 2020, 151, 40–43. [Google Scholar] [CrossRef]
  20. Ohta, M. Stability of stationary states for coupled Klein-Gordon-Schrödinger equations. Nonlinear Anal. Theory, Methods Appl. 1996, 27, 455–461. [Google Scholar] [CrossRef]
  21. Huang, C.; Guo, B.; Huang, D.; Li, Q. Global well-posedness of the fractional Klein-Gordon-Schrödinger system with rough initial data. Sci. China Math. 2016, 59, 1345–1366. [Google Scholar] [CrossRef]
  22. Kong, L.H.; Chen, M.; Yin, X.L. A novel kind of efficient symplectic scheme for Klein-Gordon-Schrödinger equation. Appl. Numer. Math. 2019, 135, 481–496. [Google Scholar] [CrossRef]
  23. Wang, J.J.; Xiao, A.G. A efficient conservation difference scheme for fractional Klein-Gordon-Schrödinger equations. Appl. Math. Comput. 2018, 320, 691–709. [Google Scholar] [CrossRef]
  24. Wang, J.J.; Xiao, A.G. Conservative fourier spectral method and numerical investigation of space fractional Klein-Gordon-Schrödinger equations. Appl. Math. Comput. 2019, 350, 348–365. [Google Scholar] [CrossRef]
  25. McLachlan, R.I.; Quispel, G.R.W.; Robidoux, N. Geometric integration using discrete gradients. Philos. Trans. R Biol. Sci. 1999, 357, 1021–1045. [Google Scholar] [CrossRef]
  26. Quispel, G.R.W.; Mclaren, D.I. A new class of energy-preserving numerical integration methods. J. Phys. A Math. Theor. 2008, 41, 045206. [Google Scholar] [CrossRef]
  27. Furihata, D.; Matsuo, T. Discrete Variational Derivative Method, a Structure-Preserving Numerical Method for Partial Differential Equations; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  28. Brugnano, L.; Iavernarivo, F.; Trigiante, D. Hamiltonian boundary value methods (Energy preserving discrete line integral methods). J. Numer. Anal. Appl. Math. 2010, 5, 17–37. [Google Scholar]
  29. Celledonia, E.; Grimmb, V.; McLachlanc, R.I.; McLarend, D.I.; O’Nealed, D.; Owrena, B.; Quispel, G.R.W. Preserving energy resp. dissipation in numerical PDEs using the “Average Vector Field” method. J. Comput. Phys. 2012, 231, 6770–6789. [Google Scholar] [CrossRef]
  30. Jiang, C.L.; Sun, J.Q.; He, X.F. High order energy-preserving method of the Good Boussinesq equation. Numer. Methods Appl. 2016, 9, 111–122. [Google Scholar] [CrossRef]
  31. Bridges, T.J. Multi-symplectic structures and wave propagation. Math. Proc. Camb. Phil. Soc. 1997, 121, 147–190. [Google Scholar] [CrossRef]
  32. Chen, J.B. Symplectic and multisymplectic Fourier pseudospectral discretizations for the Klein-Gordon equation. Lett. Math. Phys. 2006, 75, 293–305. [Google Scholar] [CrossRef]
Figure 1. Evolution of solitary waves (a) u ( x , t ) and (b) m ( x , t ) at α = 2 , t [ 0 , 10 ] .
Figure 1. Evolution of solitary waves (a) u ( x , t ) and (b) m ( x , t ) at α = 2 , t [ 0 , 10 ] .
Fractalfract 07 00711 g001
Figure 2. Evolution of solitary waves (a) u ( x , t ) and (b) m ( x , t ) at α = 1.7 , t [ 0 , 16 ] .
Figure 2. Evolution of solitary waves (a) u ( x , t ) and (b) m ( x , t ) at α = 1.7 , t [ 0 , 16 ] .
Fractalfract 07 00711 g002
Figure 3. Energy error of fractional KGZ equation: (a) α = 2 , (b) α = 1.7 at t [ 0 , 20 ] .
Figure 3. Energy error of fractional KGZ equation: (a) α = 2 , (b) α = 1.7 at t [ 0 , 20 ] .
Fractalfract 07 00711 g003
Figure 4. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 2 , t [ 0 , 20 ] .
Figure 4. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 2 , t [ 0 , 20 ] .
Fractalfract 07 00711 g004
Figure 5. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 1.2 , t [ 0 , 20 ] .
Figure 5. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 1.2 , t [ 0 , 20 ] .
Fractalfract 07 00711 g005
Figure 6. Energy error of the fractional KGS equation: (a) α = 2 , (b) α = 1.2 at t [ 0 , 20 ] .
Figure 6. Energy error of the fractional KGS equation: (a) α = 2 , (b) α = 1.2 at t [ 0 , 20 ] .
Fractalfract 07 00711 g006
Figure 7. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 2 , t [ 0 , 30 ] .
Figure 7. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 2 , t [ 0 , 30 ] .
Fractalfract 07 00711 g007
Figure 8. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 1.6 , t [ 0 , 30 ] .
Figure 8. Evolution of solitary waves (a) f ( x , t ) and (b) h ( x , t ) at α = 1.6 , t [ 0 , 30 ] .
Fractalfract 07 00711 g008
Figure 9. Energy errors of the fractional KGS equations: (a) α = 2 , (b) α = 1.6 at t [ 0 , 30 ] .
Figure 9. Energy errors of the fractional KGS equations: (a) α = 2 , (b) α = 1.6 at t [ 0 , 30 ] .
Fractalfract 07 00711 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, J.; Yang, S.; Zhang, L. Energy-Preserving AVF Methods for Riesz Space-Fractional Nonlinear KGZ and KGS Equations. Fractal Fract. 2023, 7, 711. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7100711

AMA Style

Sun J, Yang S, Zhang L. Energy-Preserving AVF Methods for Riesz Space-Fractional Nonlinear KGZ and KGS Equations. Fractal and Fractional. 2023; 7(10):711. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7100711

Chicago/Turabian Style

Sun, Jianqiang, Siqi Yang, and Lijuan Zhang. 2023. "Energy-Preserving AVF Methods for Riesz Space-Fractional Nonlinear KGZ and KGS Equations" Fractal and Fractional 7, no. 10: 711. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7100711

Article Metrics

Back to TopTop