Next Article in Journal
Application of Dandelion Optimization Algorithm in Pattern Synthesis of Linear Antenna Arrays
Previous Article in Journal
A Case Study of Accident Analysis and Prevention for Coal Mining Transportation System Based on FTA-BN-PHA in the Context of Smart Mining Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform

School of Advanced Manufacturing, Guangdong University of Technology, Jieyang 522000, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(7), 1110; https://doi.org/10.3390/math12071110
Submission received: 27 February 2024 / Revised: 2 April 2024 / Accepted: 5 April 2024 / Published: 7 April 2024
(This article belongs to the Section Mathematical Physics)

Abstract

:
This paper introduces a novel approach employing the fast cosine transform to tackle the 2-D and 3-D fractional nonlinear Schrödinger equation (fNLSE). The fractional Laplace operator under homogeneous Neumann boundary conditions is first defined through spectral decomposition. The difference matrix Laplace operator is developed by the second-order central finite difference method. Then, we diagonalize the difference matrix based on the properties of Kronecker products. The time discretization employs the Crank–Nicolson method. The conservation of mass and energy is proved for the fully discrete scheme. The advantage of this method is the implementation of the Fast Discrete Cosine Transform (FDCT), which significantly improves computational efficiency. Finally, the accuracy and effectiveness of the method are verified through two-dimensional and three-dimensional numerical experiments, solitons in different dimensions are simulated, and the influence of fractional order on soliton evolution is obtained; that is, the smaller the alpha, the lower the soliton evolution.

1. Introduction

The Nonlinear Schrödinger Equation (NLSE) is a nonlinear partial differential equation that describes the evolution of wave functions in quantum mechanics. It finds extensive applications across various fields of physics, including optics [1], cold atomic physics [2], and plasma physics [3]. In nonlinear optics, the NLSE serves as a fundamental mathematical model to depict the behavior of light waves in nonlinear media. In semiconductor physics, the NLSE describes the propagation of optical solitons in optical fibers made from semiconductor materials [4]. The NLSE is crucial in studying the nonlinear dynamics of semiconductor lasers, particularly in mode-locking and the generation of ultra-short optical pulses [5]. Moreover, it has wide-ranging applications in all-optical communications and all-optical information storage [6]. In recent years, many scholars have introduced fractional-order derivatives into the NLSE, leading to the fractional nonlinear Schrödinger equation (fNLSE). This modification enriches and complicates the dynamical behavior of systems. The fNLSE has been extensively used in various fields such as nonlinear optics [7], Bose–Einstein condensation [8], electromagnetics [9], and quantum mechanics [10]. In nonlinear optics, this equation is employed to describe and predict complex optical phenomena induced by nonlinear effects in optical media, including signal transmission in optical fibers [11], the evolution of wave packets [12], and the interaction among optical solitons [13]. Numerical simulation of solitons through difference equations and exploration of soliton applications in different fields have become hot topics in research. Solitons have been widely applied in various fields such as hydrodynamics [14], dynamics [15,16], ocean engineering [17], etc. Extracting new soliton solutions to study the hidden physical conditions of nonlinear fractional-order partial differential equations has become an important research direction [18,19,20].
We study the following 2-D fNLSE:
i u t ( x , t ) ( Δ u ( x , t ) ) α 2 + V ( x ) u ( x , t ) + β u ( x , t ) 2 u ( x , t ) = 0 , x Ω , t ( 0 , T ]
where Ω represents a 2-D or 3-D with homogeneous Neumann boundary conditions, u ( x , t ) is a complex function, i 2 = 1 is the imaginary unit, V ( x ) is the potential function, and  β is a constant. The symbol ( Δ u ) α 2 denotes the fractional Laplace operator, where α is a positive real number, and  α = 2 stands for the classical Laplace operator. The fractional Laplace operator can be defined in various ways, commonly through Fourier transform methods, integral representations, fractional Sobolev space approaches, and spectral definitions. Since we discuss the homogeneous Neumann boundaries in this paper, the spectral definition method is adopted as follows [21]:
( Δ ) α 2 u x = k = 1 c k μ k α 2 φ k x ,
where μ k is the eigenvalues of ( Δ u ) , and φ k corresponds to the eigenfunction,
c k = Ω u φ k d x
It is well-known that the Nonlinear Schrödinger Equation (NLSE) satisfies the conservation of mass and energy. The development of conservation numerical methods has been a research focus. Numerous scholars have developed various conservation methods. Hendy et al. proposed a method combining finite difference/spectrum and Galerkin–Legendre techniques to solve the coupled nonlinear space–time fractional Schrödinger equation, which exhibits non-smooth solutions in the time domain [22]. Li et al. introduced scalar auxiliary variables to reformulate the Schrödinger equation into a new family of systems. These systems were approximated using the implicit midpoint method, repeated ladder method, and fractional center difference method [23]. Liaqat et al. proposed a novel combinatorial calculation method using the conformable natural transform (CNT) and homotopy perturbation method (HPM) to derive analytical and numerical solutions for the time-fractional suitable Schrödinger equation (TFCSE) [24]. Kaabar et al. defined a new generalized double Laplace transform, coupled with the Adomian decomposition method, to solve the newly formulated nonlinear Schrödinger equation with spatiotemporal dispersion [25]. Zhang [26] investigated the optical soliton solutions of the nonlinear Schrödinger equation with a quintic non-Kerr nonlinear term describing the nonlinear wave state of optical solitons, which is a noteworthy and important model in optical fiber communication. Wang and Huang [27,28] investigated energy-conserving difference schemes combined with Alternating Direction Implicit (ADI) methods. Yang derived a linearized energy-conserving finite difference scheme for a class of nonlinear fractional NLSE equations, studying their energy conservation and convergence properties [29]. Klein et al. proposed a Fourier spectral method for the one-dimensional fractional NLSE [30]. Previous studies mainly focus on Dirichlet boundary conditions and periodic boundary conditions. For the problems with Neumann boundary conditions, to the best of the authors’ knowledge, the relevant papers are limited.
This paper presents a novel, rapid, and conservative method for solving the fractional NLSE under Neumann boundaries in 2-D and 3-D cases. We first develop the difference matrix of the 2-D and 3-D Laplace operators using the finite difference method and then diagonalize the matrix based on the properties of Kronecker products. The difference matrix of the fractional Laplace operator can be obtained by the spectral decomposition. Time discretization employs the Crank–Nicolson method, and the conservation of mass and energy in the fully discrete scheme is proved. The advantage of this method is the implementation of the Fast Discrete Cosine Transform (FDCT). The significant improvement in computational efficiency can be observed in the numerical tests.
This paper is arranged as follows. In Section 2, we discretized 2-D and 3-D Equation (1) by using central difference and the Crank–Nicolson method subsequently. In Section 3, we first introduce the Kronecker operator to simplify the calculation. The conservation of energy and mass of the fully discrete scheme is presented. The implementation of FDCT is discussed in detail in Section 4. In Section 5, we provide numerical examples to verify the efficiency, conservation and accuracy of the method.

2. Numerical Schemes

In this section, the numerical scheme for the 2-D fractional nonlinear Schrödinger equation is first introduced, and then the schemes are generalized to the 3-D domain.

2.1. 2-D Case

In this subsection, the rectangular computing domain Ω = a , b × c , d is considered, where the boundaries are the zero Neumann boundaries. The domain is divided into N x × N y grids. The coordinates of the grid points as follows: specifically, the domain is divided into N x × N y grids by some lines paralleling the axes, where the step size in the x direction is h x = b a N x and in the y direction is h y = d c N y . Let u j , k be the numerical solution of u ( x , y ) , which are assigned at the center points of the meshes:
T h = x j , y k x j = a + ( j 1 2 ) h x , y k = c + ( k 1 2 ) h y , j = 1 , N x , k = 1 , N y .
Suppose the matrix U for the discrete solution on gird center points, i.e., that U = u j , k N x × N y . Applying the central difference method in the x-direction, we obtain the second-order scheme
u x x j , k u j 1 , k 2 u j 1 , k + u j 1 , k h x 2 , j = 1 , 2 , 3 , N x 1 .
To obtain a second-order accurate method, we might introduce another unknown u 0 , k on the ghost point. For the homogeneous Neumann boundary condition, we can use the centered approximation on the left boundary point to obtain
u 2 , k u 0 , k 2 h x = 0 .
The proceeding two equations result in a difference matrix in the x direction,
K x = 1 1 0 0 0 1 2 1 0 0 0 1 2 1 0 0 0 0 1 2 1 0 0 0 0 1 1 N x × N x ,
and the difference matrix K y in the y direction can be obtained similarly. We can now obtain the discrete scheme of the Laplacian operator as:
Δ u j , k 1 h x 2 K x U + 1 h y 2 U K y j , k .
The tridiagonal matrices K x and K y can be diagonalized as [31]
K x = S x Λ x S x 1 , K y = S y Λ y S y 1
Here, the eigenvector matrices S x and S y are discrete cosine transform matrices, S x j , k = c o s j π N x ( k 1 2 ) , S y j , k = c o s j π N y ( k 1 2 ) , and both are orthogonal matrices. The matrixes Λ x and Λ y are diagonal matrices consisting of eigenvalues, represented in the following form:
Λ x = d i a g λ 1 , λ 2 , , λ N x , λ j = 2 + 2 c o s ( j π / N x )
Λ y = d i a g λ 1 , λ 2 , , λ N y , λ j = 2 + 2 c o s ( j π / N y )
Next, we will rewrite the discrete scheme (5) in matrix-vector multiplication form by the Kronecker product. The Kronecker product A B is defined by multiplying each element of matrix A by the entire matrix B:
A B = a 11 B a 12 B a 1 n B a 21 B a 22 B a 2 n B a m 1 B a m 2 B a m n B ,
where A = a i j m × n and B = b i j p × q . Set arrays U = u j , k N x × N y and vectorize the solution matrix U as
U = V e c U = u 1 , 1 u 2 , 1 u N x , 1 , , u 1 , N y u 2 , N y u N x , N y T .
The elements of the diagonal matrix V are obtained by vectorizing two-dimensional array V.
Set the matrix K = 1 h x 2 K x I y + I x 1 h y 2 K y by the Kronecker product. Then, the discrete scheme (5) can be rewritten into the matrix-vector multiplication form:
K U = 1 h x 2 K x I y + I x 1 h y 2 K y U
Here, the matrixes I x , I y are N x -th and N y -th order identity matrix, respectively. Set the time step τ = T N , t n = n τ , ( n = 0 , 1 , , N 1 ) . The fully discrete scheme of (1) is obtained by using the Crank–Nicolson method as:
i U n + 1 U n τ K α 2 V U n + 1 + U n 2 + β 2 U n + 1 2 + U n 2 U n + 1 + U n 2 = 0
where U n and U n + 1 represent the column vector of U in Equation (7) at different times t, and the computation of the matrix power K α 2 will be presented in Section 3.

2.2. 3-D Case

Next, we consider the 3-D case. The computation domain Ω = a , b × c , d × e , f is divided into N x × N y × N z grids along the x, y, and z directions, respectively. Use x j , y k , z l to represent the center of the meshes, which is defined similarly as in (4). The time step remains consistent with the 2-D case. Let the 3-D array U ˜ = u j , k , l N x × N y × N z represent the numerical solution at x j , y k , z l , then vectorize the solution array as:
U ˜ = V e c U ˜ = u 1 , 1 , 1 , , u 1 , 1 , N z , u 1 , 2 , 1 , , u 1 , 2 , N z , , u 1 , N y , N z , u 2 , 1 , 1 , , u N x , N y , N z T
We apply the second-order central difference scheme in space and the Crank–Nicolson in time for Equation (1) in 3-D. The numerical scheme in the matrix-vector multiplication form can be obtained as:
i U ˜ n + 1 U ˜ n τ K ˜ α 2 V ˜ U ˜ n + 1 + U ˜ n 2 + β 2 U ˜ n + 1 2 + U ˜ n 2 U ˜ n + 1 + U ˜ n 2 = 0
The difference matrix is defined as:
K ˜ = I z I y K x + I z K y I x + K z I y I x
where I z is the identity matrix of order N z .

3. Conservation

Before proving the conservative properties of the difference scheme, some notations are first introduced. Let u , v C N x × N y be a complex grid function on grids T h in (4). The inner product is defined as: ( u , v ) = h x h y j = 1 N x k = 1 N y u j , k v j , k ¯ , where v j , k ¯ is the conjugate of v j , k . The definition of the l 2 -norm of grid function v is v = v , v . We also define the l p -norm as v p = h x h y j = 1 N x k = 1 N y v j , k p 1 p , 0 < p < . When p = , we obtain the l -norm v = max 1 j N 1 | u j , k | . The inner product and l p -norm in the 3-D case can be similarly defined and we omit it here.
The conservation of mass and energy will be proved under the fully discrete scheme (8) and (10). Before proceeding, it is necessary to decompose K in (8). To make the paper self-contained, we first present some important properties of the Kronecker product:
i.
( A B ) ( C D ) = A C B D
ii.
( A B ) 1 = A 1 B 1
iii.
( A B ) T = A T B T
iv.
B T A v e c X = v e c A X B
Based on the properties of the Kronecker product, we give the following Lemmas.
 Lemma 1. 
The 2-D difference matrix K in (8) is symmetric positive definite matrices, and it follows that K α 2 = L T L . And the 3-D difference matrix in (11) can also be decomposed as K ˜ α 2 = L ˜ T L ˜ .
Proof. 
Firstly, we prove the decomposition of the 2-D difference matrix K . Since the discrete cosine transform matrices S x , S y are orthogonal, the matrices K x , K y in (5) can be written as K x = S x Λ x S x T , K y = S y Λ y S y T . Furthermore, by the properties (ii) and (iii) of the Kronecker product, we have:   
1 h x 2 K x I y = S x 1 h x 2 Λ x S x T I y = S x 1 h x 2 Λ x S x T S y I y S y T = S x S y 1 h x 2 Λ x S x T I y S y T = S x S y 1 h x 2 Λ x I y S x T S y T = S x S y 1 h x 2 Λ x I y S x S y T
similarly, we obtain:
I x 1 h y 2 K y = S x S y I x 1 h y 2 Λ y S x S y T .
From the above two equations, we have:
K α 2 = S x S y 1 h x 2 Λ x I y I x 1 h y 2 Λ y α 2 S x S y T
It is obvious that Λ = 1 h x 2 Λ x I y I x 1 h y 2 Λ y is a diagonal matrix with positive diagonal elements. Therefore, we can compute the fractional power of the matrix as K α 2 = S x S y Λ α 4 Λ α 4 S x S y T = L T L , where L = Λ α 4 S x S y T .
Next, we consider the decomposition of K ˜ in (10). Based on the properties (i) and (ii) of the Kronecker product, we have:
I z I y G x = S z I z S z T S y I y S y T S x Λ x S x T       = S z S y S x I z I y Λ x S z T S y T S x T
A similar computation can be applied to I z G y I x , G z I y I x . Then, we can rewrite K ˜ as:
K ˜ = S z S y S x Λ ˜ S z T S y T S x T
where Λ ˜ = I z I y 1 h x 2 Λ x + I z 1 h y 2 Λ y I x + 1 h z 2 Λ z I y I x is a diagonal matrix with positive diagonal elements. The fractional power of the matrix K ˜ can be decomposed as K ˜ α 2 = S z S y S x Λ ˜ α 4 Λ ˜ α 4 S z T S y T S x T T = L ˜ T L ˜ .    □
Now, we complete Lemma 1. The following proofs in 3-D have analogs for the 2-D case, but the formulas are different, i.e., the formulas in the 3-D case are accompanied by a tilde above them. The proof below does not need to be written for both cases, so we will only present the proofs in the 2-D case. By Lemma 1, Lemma 2 can be readily obtained.
 Lemma 2. 
L T L U , U = L U , L U = L U 2 .
With the l 2 - and l p -norm, the definitions of mass and energy conservation are given as Q n = U n 2 and E n = L U n 2 V U n 2 β 2 U n 4 , respectively.
 Lemma 3. 
For a complex vector U , the following equation holds:
(1) 
Im ( K α 2 U n + 1 2 , U n + 1 2 ) = 0
(2) 
Im U n + 1 2 + U n 2 U n + 1 2 , U n + 1 2 = 0
(3) 
Re K α 2 U n + 1 2 , U n + 1 U n = 1 2 L U n + 1 2 L U n 2
(4) 
Re β 2 U n + 1 2 + U n 2 U n + 1 2 , U n + 1 U n = β 4 U n + 1 4 U n 4 ,
where Re denotes the real part and Im denotes the imaginary part. The proof is straightforward by using Lemma 1 and Lemma 2. For more details, the reader is referred to the existing works in [27].
 Theorem 1. 
The fully discrete format (5) satisfies discrete mass conservation:
Q n + 1 = Q n = = Q 0
Proof. 
Taking the inner product of both sides of Equation (10) with U n + 1 + U n and taking the imaginary part, using Lemma 3 (1) and (2), we have:
U n + 1 2 U n 2 = 0
U n + 1 2 = U n 2
Thus, Q n + 1 = Q n , and the proof is complete by iteration.    □
 Theorem 2. 
The fully discrete format (9) satisfies discrete energy conservation.
E n + 1 = E n = = E 0
Proof. 
Taking the inner product of both sides of Equation (9) with U n + 1 U n and taking the real part, using Lemma 3 (3) and (4), we have:
L U n + 1 2 U n + 1 2 β 2 U n + 1 4 = L U n 2 V U n 2 β 2 U n 4
Therefore, we obtain E n + 1 = E n . This immediately implies discrete energy conservation.    □

4. Fast Implementation

Next, we will give the implementation of the solution of Equation (10). The implementation in 3-D is similar to the 2-D case, so we provide the computation with the 2-D case as an example.
i τ I K α 2 V 2 U n + 1 = i τ I + K α 2 V 2 U n + F n + 1 2
where F n + 1 2 = β 4 U n + 1 2 + U n 2 U n + 1 + U n and I is the N x N y -th order identity matrix. Equation (17) becomes:
U n + 1 = A 1 B U n + A 1 F n + 1 2 ,
where the matrices A and B are defined as:
A = S x S y i τ I 1 2 1 h x 2 Λ x I y I x 1 h y 2 Λ y α 2 + V 2 S x S y 1 B = S x S y i τ I + 1 2 1 h x 2 Λ x I y I x 1 h y 2 Λ y α 2 V 2 S x S y 1
according to Equation (17). We solve the above nonlinear system (18) by Picard iteration. The specific iteration algorithm is as following Algorithm 1 [32]:
Algorithm 1: The Picard iteration for the nonlinear system (18)
U 0 = U 0 x j , y k f o r n = 0 , 1 , N U 0 n + 1 = U n f o r p = 1 , 2 U n + 1 p + 1 = A 1 B U n + A 1 F ( U n ) + F ( U n + 1 p ) 2 w h i l e | | U n + 1 p + 1 U n + 1 p | | < ε = 10 13 c o n t i n u e e n d U n = U n + 1 p + 1 e n d
Next, we will use a fast discrete cosine transform to solve Equation (18). Let
Λ 1 = i τ I 1 h x 2 Λ x I y I x 1 h y 2 Λ y α 2 + V 2 , Λ 2 = i τ I + 1 h x 2 Λ x I y I x 1 h y 2 Λ y α 2 V 2
we rewrite Equation (18) as the following equations:
U n + 1 = S x S y Λ 1 1 Λ 2 S x S y 1 U n + S x S y Λ 1 1 S x S y 1 F n + 1 2
Due to the similarity in computation steps between the two terms on the right-hand side of Equation (19), we take the second term as an example. The computational procedure is divided into three steps:
Step 1:
According to the properties of the Kronecker product (ii) and (iv), the matrix-vector multiplication of S x S y 1 F in Equation (17) can be achieved as follows:
S x S y 1 v e c F = S x 1 S y 1 v e c F = v e c S y 1 F S x
It can be achieved using the fast discrete cosine transform.
Step 2:
Since it is a diagonal matrix, Λ 1 1 v e c S y 1 F S x = v e c S y 1 F S x . / Λ , where Λ 1 is represented by . / Λ .
Step 3:
Since S x = S x T , then S x S y v e c S y 1 F S x . / Λ = v e c ( S y S y 1 F S x . / Λ S x ) . This part can be achieved using the fast inverse discrete cosine transform.

5. Numerical Experiments and Discussion

5.1. Numerical Experiments

In this section, we report some numerical results of the 2-D and 3-D fNLS Equation (1) to support our theoretical analysis.
 Example 1. 
We consider the problem (1) in Ω = π , π × π , π , with the potential function V ( x , y ) = 1 4 cos 2 2 x sin 2 2 y . When α = 2 , the problem collapses to the classical cubic nonlinear Schrödinger equation, and the exact solution is given by u ( x , y , t ) = e 9 i t cos 2 x sin 2 y . In this example, we compute the l -norm errors e ( τ , h ) = u ( T ) U N t of the numerical solution at T = 1. The time step is fixed with a small time step, τ = 1 × 10 4 .
Firstly, we test the discrete mass conservation law. Figure 1 gives the evolution of mass Qn and energy En for 1.6, 1.7, 1.8, 2, respectively, with τ = h = 0.01.
In Table 1, we compute the maximum norm errors of the energy and mass of the numerical solution with τ = h = 0.01 at t = 1 , which helps to better evaluate the accuracy and stability of the numerical method. Then, the accuracy of the example is demonstrated by computing the error between the analytical solution and the numerical solution, with the results presented in Table 2. Table 2 shows that the scheme has second-order accuracy in space. It is worth noting that the scheme with DCT performs the computation in CPU time, which is several orders of magnitude shorter with respect to the no DCT method. Table 3 is obtained by scheme (9) with different time steps. The spatial step is also chosen to be relatively small (h = 0.05). The results from Table 3 indicate that our scheme is second-order in time.
 Example 2. 
Now, we investigate the impact of the collision of two solitons brought by fractional order. In this example, we compute the interaction between two solitons. Take V ( x , y ) = 0 , β = 2 and Ω = [ 30 , 60 ] × [ 30 , 60 ] in (1). The initial conditions are given as follows: u ( x , y , 0 ) = e 2 i ( x + y ) sech x sech y + e 0.05 i ( x + y 15 ) sech ( x 15 ) sech ( y 15 ) . In order to illustrate the fact clearly, we first specify the elastic collisions with α = 2 . 0 . It can be observed that the waves retain their shape and velocity after interaction (see Figure 2). Then, we choose different α to test the collisions of two solitons. Figure 3 shows the evolution of the modulus of the collision at t = 4 for α = 1.6 , 1.7 , 1.8 , respectively. We can find that the order α will greatly affect the collision time. The smaller α becomes, the longer the collision time. Generally, a reduction in the order may indicate a stronger memory effect in the system, leading to a slower interaction between isolated solitons. This feature is consistent with the fact that the fractional order introduces more historical information.
 Example 3. 
The solutions of the nonlinear Schrödinger equation can exhibit a phenomenon known as ‘blow-up’. In this example, we show the singular solutions for the FNLS equation. We choose V ( x , y ) = 0 , β = 1 and Ω = [ 10 , 10 ] × [ 10 , 10 ] in (1). The initial condition is u ( x , y , 0 ) = 6 2 exp ( ( x 2 + y 2 ) ) . Figure 4 shows the modulus of the solution u ( x , y ) with different fractional order, α = 1.3 , α = 1.5 , α = 1.8 , α = 2 . 0 . The blow-up effect is obtained in finite time with different α. These plots show that the blow-up time becomes smaller for progressively increasing α. This is because the fractional-order derivatives introduce a dependence on the system’s past states, causing the evolution of the system to be influenced by previous states over longer distances.
 Example 4. 
Finally, we consider the 3-D problem (1) in Ω = [ 10 , 10 ] 3 , with V ( x , y ) = 0 , β = 2 . We set the final computation time as t = 1. The initial conditions are u ( x , y , z , 0 ) = sech ( x ) sech ( y ) sech ( z ) exp ( i ( x + y + z ) ) . To better illustrate the results, Figure 5 describes the time evolutions of solitons. Figure 5a shows the initial solution with the isosurfaces u = 0.25 . Figure 5b–d demonstrate the solution of u = 0.25 , where the fractional orders are chosen as α = 1.2 , 1.6 , and 2. It can be seen from the figure that the order α will greatly affect the soliton. The soliton propagates at a slower speed when α decreases.

5.2. Discussion

In Figure 1, the conservation of mass E and energy Q is clearly depicted, alongside the discernible trend that, under varying α , the mass remains nearly constant while the energy Q decreases with α . Figure 2 presents simulation diagrams illustrating optical solitons with integer orders at different time intervals. These diagrams vividly show the loss and diffraction phenomena occurring during the propagation of optical solitons. By comparing with Figure 2, Figure 3 elucidates the influence of α on the motion of optical solitons: as α decreases, the propagation speed of optical solitons decelerates. The occurrence of blasting phenomena significantly impacts the stability of optical soliton propagation. In Figure 4, the blasting time for different α when solving the Schrödinger equation is demonstrated. To accentuate the blasting phenomenon, a three-dimensional image is provided for better observation, revealing that smaller α results in slower blasting phenomena. Figure 5 presents the isosurface diagram concerning u, offering a theoretical foundation for the study of optical solitons.

6. Conclusions

This paper solves the two-dimensional and three-dimensional fractional nonlinear Schrödinger equation under the Neumann boundaries. Our scheme employs Kronecker products to represent the differential matrix of the fully discrete equation. This is more helpful for our proof of the conservation of mass and energy. In addition, it obviously improves the computing speed. In the implementation process, the application of the Fast Discrete Cosine Transform effectively reduces computation time. Furthermore, extensive numerical experiments in the 2-D case and the 3-D case are given to confirm the method’s accuracy, efficiency, and stability. By varying fractional order, we could observe the impact of fractional order. Specifically, we noted that the smaller the fractional order, the slower the movement and formation of optical solitons. Furthermore, the approach delineated in this article extends its applicability beyond the specific differential equation under consideration. It not only serves to verify the conservation properties of various other differential equations but also introduces an innovative methodology for their solutions. This provides a theoretical basis for the study of optical soliton propagation in optical fibers and nonlinear optics. This highlights the need for further study in the future.

7. Future Work

In the future, we will focus on the study of nonlinear optical solitons depicted by different differential equations and extend the methodology to other nonlinear models.

Author Contributions

Conceptualization, R.Z.; methodology, P.W.; software, P.W.; validation, S.P.; writing—original draft, Y.C.; writing—review and editing, R.Z.; supervision, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ibarra-Villalon, H.E.; Pottiez, O.; Gómez-Vieyra, A.; Lauterio-Cruz, J.P.; Bracamontes-Rodriguez, Y.E. Numerical approaches for solving the nonlinear Schrödinger equation in the nonlinear fiber optics formalism. J. Opt. 2020, 22, 043501. [Google Scholar] [CrossRef]
  2. Vowe, S.; Lämmerzahl, C.; Krutzik, M. Detecting a logarithmic nonlinearity in the Schrödinger equation using Bose-Einstein condensates. Phys. Rev. A 2020, 101, 043617. [Google Scholar] [CrossRef]
  3. Sultana, S. Review of heavy-nucleus-acoustic nonlinear structures in cold degenerate plasmas. Rev. Mod. Plasma Phys. 2022, 6, 6. [Google Scholar] [CrossRef]
  4. Rao, J.-G.; Chen, S.-A.; Wu, Z.-J.; He, J.-S. General higher-order rogue waves in the space-shifted symmetric nonlocal nonlinear Schrödinger equation. Acta Phys. Sin. 2023, 72, 104204-1–104204-9. [Google Scholar] [CrossRef]
  5. Li, M.; Wang, B.T.; Xu, T.; Shui, J.J. Study on the generation mechanism of bright and dark solitary waves and rogue wave for a fourth-order dispersive nonlinear Schrödinger equation. Acta Phys. Sin. 2020, 69, 010502-1–010502-10. [Google Scholar] [CrossRef]
  6. Wen, J.-M.; Bo, W.-B.; Wen, X.-K.; Dai, C.-Q. Multipole vector solitons in coupled nonlinear Schrödinger equation with saturable nonlinearity. Acta Phys. Sin. 2023, 72, 100502-1–100502-7. [Google Scholar] [CrossRef]
  7. Ahmad, J.; Akram, S.; Noor, K.; Nadeem, M.; Bucur, A.; Alsayaad, Y. Soliton solutions of fractional extended nonlinear Schrödinger equation arising in plasma physics and nonlinear optical fiber. Sci. Rep. 2023, 13, 10877. [Google Scholar] [CrossRef] [PubMed]
  8. Jiang, T.; Huang, J.-J.; Lu, L.-G.; Ren, J.-L. Numerical study of nonlinear Schrödinger equation with high-order split-step corrected smoothed particle hydrodynamics method. Acta Phys. Sin. 2019, 68, 090203-1–090203-14. [Google Scholar] [CrossRef]
  9. Qureshi, S.; Chang, M.M.; Shaikh, A.A. Analysis of series RL and RC circuits with time-invariant source using truncated M, Atangana beta and conformable derivatives. J. Ocean Eng. Sci. 2021, 6, 217–227. [Google Scholar] [CrossRef]
  10. Stephanovich, V.A.; Olchawa, W.; Kirichenko, E.V.; Dugaev, V.K. 1D solitons in cubic-quintic fractional nonlinear Schrödinger model. Sci. Rep. 2022, 12, 15031. [Google Scholar] [CrossRef] [PubMed]
  11. Islam, Z.; Abdeljabbar, A.; Sheikh, M.A.N.; Taher, M.A. Optical solitons to the fractional order nonlinear complex model for wave packet envelope. Results Phys. 2022, 43, 106095. [Google Scholar] [CrossRef]
  12. Xie, P.; Zhu, Y. Wave Packets in the Fractional Nonlinear Schrödinger Equation with a Honeycomb Potential. Multiscale Model. Simul. 2021, 19, 951–979. [Google Scholar] [CrossRef]
  13. Riaz, M.B.; Atangana, A.; Jahngeer, A.; Jarad, F.; Awrejcewicz, J. New optical solitons of fractional nonlinear Schrodinger equation with the oscillating nonlinear coefficient: A comparative study. Results Phys. 2022, 37, 105471. [Google Scholar] [CrossRef]
  14. Shahen, N.H.M.; Rahman, M.M. Dispersive solitary wave structures with MI Analysis to the unidirectional DGH equation via the unified method. Partial Differ. Equ. Appl. Math. 2022, 6, 100444. [Google Scholar]
  15. An, T.; Shahen, N.H.M.; Ananna, S.N.; Hossain, M.F.; Muazu, T. Exact and explicit travelling-wave solutions to the family of new 3D fractional WBBM equations in mathematical physics. Results Phys. 2020, 19, 103517. [Google Scholar]
  16. Shahen, N.H.M.; Rahman, M.M.; Alshomrani, A.S.; Inc, M. On fractional order computational solutions of low-pass electrical transmission line model with the sense of conformable derivative. Alex. Eng. J. 2023, 81, 87–100. [Google Scholar]
  17. Iqbal, M.A.; Miah, M.M.; Ali, H.S.; Shahen, N.H.M.; Deifalla, A. New applications of the fractional derivative to extract abundant soliton solutions of the fractional order PDEs in mathematics physics. Partial Differ. Equ. Appl. Math. 2024, 9, 100597. [Google Scholar] [CrossRef]
  18. Justin, M.; David, V.; Shahen, N.H.M.; Sylvere, A.S. Sundry optical solitons and modulational instability in Sasa-Satsuma model. Opt. Quantum Electron. 2022, 54, 1–15. [Google Scholar] [CrossRef]
  19. Shahen, N.H.M.; Ali, M.S.; Rahman, M.M. Interaction among lump, periodic, and kink solutions with dynamical analysis to the conformable time-fractional Phi-four equation. Partial Differ. Equ. Appl. Math. 2021, 4, 100038. [Google Scholar] [CrossRef]
  20. Shahen, N.H.M.; Foyjonnesa Bashar, M.H.; Tahseen, T.; Hossain, S. Solitary and rogue wave solutions to the conformable time fractional modified kawahara equation in mathematical physics. Adv. Math. Phys. 2021, 2021, 6668092. [Google Scholar] [CrossRef]
  21. Lischke, A.; Pang, G.; Gulian, M.; Song, F.; Glusa, C.; Zheng, X.; Mao, Z.; Cai, W.; Meerschaert, M.M.; Ainsworth, M.; et al. What is the fractional Laplacian? A comparative review with new results. J. Comput. Phys. 2020, 404, 109009. [Google Scholar] [CrossRef]
  22. Hendy, A.S.; Zaky, M.A. Combined Galerkin spectral/finite difference method over graded meshes for the generalized nonlinear fractional Schrödinger equation. Nonlinear Dyn. 2021, 103, 2493–2507. [Google Scholar] [CrossRef]
  23. Li, X.; Wen, J.; Li, D. Mass-and energy-conserving difference schemes for nonlinear fractional Schrödinger equations. Appl. Math. Lett. 2021, 111, 106686. [Google Scholar] [CrossRef]
  24. Liaqat, M.I.; Akgül, A. A novel approach for solving linear and nonlinear time-fractional Schrödinger equations. Chaos Solitons Fractals 2022, 162, 112487. [Google Scholar] [CrossRef]
  25. Kaabar, M.K.; Martínez, F.; Gómez-Aguilar, J.F.; Ghanbari, B.; Kaplan, M.; Günerhan, H. New approximate analytical solutions for the nonlinear frac-tional Schrödinger equation with second-order spatio-temporal dispersion via double Laplace transform method. Math. Methods Appl. Sci. 2021, 44, 11138–11156. [Google Scholar] [CrossRef]
  26. Zhang, K.; Han, T. The optical soliton solutions of nonlinear Schrödinger equation with quintic non-Kerr nonlinear term. Results Phys. 2023, 48, 106397. [Google Scholar] [CrossRef]
  27. Wang, P.; Huang, C. An energy conservative difference scheme for the nonlinear fractional Schrödinger equations. J. Comput. Phys. 2015, 293, 238–251. [Google Scholar] [CrossRef]
  28. Wang, P.; Huang, C. Split-step alternating direction implicit difference scheme for the fractional Schrödinger equation in two dimensions. Comput. Math. Appl. 2016, 71, 1114–1128. [Google Scholar] [CrossRef]
  29. Yang, Z. A class of linearized energy-conserved finite difference schemes for nonlinear space-fractional Schrödinger equations. Int. J. Comput. Math. 2016, 93, 609–626. [Google Scholar] [CrossRef]
  30. Klein, C.; Sparber, C.; Markowich, P. Numerical study of fractional nonlinear Schrödinger equations. Proc. Math. Phys. Eng. Sci. 2014, 470, 20140364. [Google Scholar] [CrossRef]
  31. LeVeque, R.J. Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems; SIAM: Philadelphia, PA, USA, 2007. [Google Scholar]
  32. Chen, Z.; Gou, Q. Piecewise Picard iteration method for solving nonlinear fractional differential equation with proportional delays. Appl. Math. Comput. 2019, 348, 465–478. [Google Scholar] [CrossRef]
Figure 1. The evolution of discrete mass and energy for different values of α .
Figure 1. The evolution of discrete mass and energy for different values of α .
Mathematics 12 01110 g001
Figure 2. The interaction between two solitons at different times with α = 2 , τ = 0.01 .
Figure 2. The interaction between two solitons at different times with α = 2 , τ = 0.01 .
Mathematics 12 01110 g002
Figure 3. The interactions between two solitons at different times with α = 1.6 , α = 1.7 , α = 1.8 .
Figure 3. The interactions between two solitons at different times with α = 1.6 , α = 1.7 , α = 1.8 .
Mathematics 12 01110 g003
Figure 4. The modulus of solution in Example 3 with different α and blow-up time.
Figure 4. The modulus of solution in Example 3 with different α and blow-up time.
Mathematics 12 01110 g004
Figure 5. The evolution of the soliton in 3-D with different fractional order.
Figure 5. The evolution of the soliton in 3-D with different fractional order.
Mathematics 12 01110 g005
Table 1. The maximal errors of the energy and mass with τ = h = 0.01 at t = 1 .
Table 1. The maximal errors of the energy and mass with τ = h = 0.01 at t = 1 .
α E(n)Q(n)
1.62.46 × 10 12 2.27 × 10 12
1.72.52 × 10 12 3.40 × 10 12
1.82.61 × 10 12 4.82 × 10 12
22.98 × 10 12 8.79 × 10 12
Table 2. The error and CPU with T = 1 and τ = 1 × 10 4 .
Table 2. The error and CPU with T = 1 and τ = 1 × 10 4 .
DCTN3264128256
Error2.6 × 10 3 6.42 × 10 4 1.61 × 10 4 4.02 × 10 5
CPU(s)0.20.41.634.12
NO DCTN3264128256
Error2.6 × 10 3 6.42 × 10 4 6.42 × 10 4 -
CPU(s)1.9752.361815.55-
Table 3. Order of convergence in terms of time.
Table 3. Order of convergence in terms of time.
τ a = 1.8a = 1.9a = 2.0
e( τ ) e( τ ) e( τ )
1/1001.49 × 10 3 8.19 × 10 4 7.67 × 10 4
1/2003.73 × 10 4 2.06 × 10 4 1.92 × 10 4
1/4009.32 × 10 5 5.12 × 10 5 4.82 × 10 5
1/8002.26 × 10 5 1.28 × 10 5 1.20 × 10 5
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

Wang, P.; Peng, S.; Cao, Y.; Zhang, R. The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform. Mathematics 2024, 12, 1110. https://doi.org/10.3390/math12071110

AMA Style

Wang P, Peng S, Cao Y, Zhang R. The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform. Mathematics. 2024; 12(7):1110. https://doi.org/10.3390/math12071110

Chicago/Turabian Style

Wang, Peiyao, Shangwen Peng, Yihao Cao, and Rongpei Zhang. 2024. "The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform" Mathematics 12, no. 7: 1110. https://doi.org/10.3390/math12071110

APA Style

Wang, P., Peng, S., Cao, Y., & Zhang, R. (2024). The Conservative and Efficient Numerical Method of 2-D and 3-D Fractional Nonlinear Schrödinger Equation Using Fast Cosine Transform. Mathematics, 12(7), 1110. https://doi.org/10.3390/math12071110

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