12−√36 | 14 | 14−√36 |
12+√36 | 14+√36 | 14 |
12 | 12 |
Citation: Júlia Novaes Matias, Gyovanna Sorrentino dos Santos Campanari, Gabriela Achete de Souza, Vinícius Marinho Lima, Ricardo José Tofano, Claudia Rucco Penteado Detregiachi, Sandra M. Barbalho. Metabolic syndrome and COVID-19[J]. AIMS Bioengineering, 2020, 7(4): 242-253. doi: 10.3934/bioeng.2020021
[1] | Junjie Wang, Yaping Zhang, Liangliang Zhai . Structure-preserving scheme for one dimension and two dimension fractional KGS equations. Networks and Heterogeneous Media, 2023, 18(1): 463-493. doi: 10.3934/nhm.2023019 |
[2] | Tingting Ma, Yayun Fu, Yuehua He, Wenjie Yang . A linearly implicit energy-preserving exponential time differencing scheme for the fractional nonlinear Schrödinger equation. Networks and Heterogeneous Media, 2023, 18(3): 1105-1117. doi: 10.3934/nhm.2023048 |
[3] | Min Li, Ju Ming, Tingting Qin, Boya Zhou . Convergence of an energy-preserving finite difference method for the nonlinear coupled space-fractional Klein-Gordon equations. Networks and Heterogeneous Media, 2023, 18(3): 957-981. doi: 10.3934/nhm.2023042 |
[4] | Jin Cui, Yayun Fu . A high-order linearly implicit energy-preserving Partitioned Runge-Kutta scheme for a class of nonlinear dispersive equations. Networks and Heterogeneous Media, 2023, 18(1): 399-411. doi: 10.3934/nhm.2023016 |
[5] | Leqiang Zou, Yanzi Zhang . Efficient numerical schemes for variable-order mobile-immobile advection-dispersion equation. Networks and Heterogeneous Media, 2025, 20(2): 387-405. doi: 10.3934/nhm.2025018 |
[6] | Jie Gu, Lijuan Nong, Qian Yi, An Chen . Two high-order compact difference schemes with temporal graded meshes for time-fractional Black-Scholes equation. Networks and Heterogeneous Media, 2023, 18(4): 1692-1712. doi: 10.3934/nhm.2023074 |
[7] | Huanhuan Li, Meiling Ding, Xianbing Luo, Shuwen Xiang . Convergence analysis of finite element approximations for a nonlinear second order hyperbolic optimal control problems. Networks and Heterogeneous Media, 2024, 19(2): 842-866. doi: 10.3934/nhm.2024038 |
[8] | Hyungjin Huh . Remarks on the Schrödinger-Lohe model. Networks and Heterogeneous Media, 2019, 14(4): 759-769. doi: 10.3934/nhm.2019030 |
[9] | Jinhu Zhao . Natural convection flow and heat transfer of generalized Maxwell fluid with distributed order time fractional derivatives embedded in the porous medium. Networks and Heterogeneous Media, 2024, 19(2): 753-770. doi: 10.3934/nhm.2024034 |
[10] | Zhe Pu, Maohua Ran, Hong Luo . Effective difference methods for solving the variable coefficient fourth-order fractional sub-diffusion equations. Networks and Heterogeneous Media, 2023, 18(1): 291-309. doi: 10.3934/nhm.2023011 |
In 2000, Laskin established the nonlinear Schrödinger (NLS) model [18] with fractional Laplacian operator by extending the Feynmann path integral to the Lévy path [17]. The derived equation is more accurate than the traditional NLS equation in describing the variation law of quantum states of nonlocal physical systems with time [16]. In recent, scholars have done a lot of research on the fractional NLS equation in theoretical analysis and numerical algorithms, and achieved fruitful results [11,13,34].
In this paper, we aim to present and analyze high-order schemes with conservation properties for the fractional CNLS equation with periodic boundary conditions that can be written in the form [28]
iφt−ϑLs2φ+ζ(|φ|2+ε|ϕ|2)φ=0, x∈Ω, 0≤t≤T, | (1.1) |
iϕt−ϑLs2ϕ+ζ(|ϕ|2+ε|φ|2)ϕ=0, x∈Ω, 0≤t≤T, | (1.2) |
where i is the imaginary unit root, Ω⊂Rd(d=1,2), ϑ, ζ and ε are real constants, φ(x,t) and ϕ(x,t) are complex valued functions, 1<s≤2, and Ls2 denotes the fractional Laplacian operator and can be expressed by the Fourier transform, namely
^Ls2φ(ξ)=|ξ|sˆφ(ξ), | (1.3) |
where ˆφ(ξ)=∫Ωφ(x)e−iξxdx represents the Fourier transform for φ(x). The fractional CNLS equation will reduce to the classical CNLS equation when s=2. Therefore, similar to the classical system, the fractional CNLS equation with periodic boundary conditions possesses the following physical invariants that do not change with time [30]
M(t):=∫Ω(|φ|2+|ϕ|2)dx≡M(0), | (1.4) |
H(t):=ϑ2∫Ω[(Ls4φ)2+(Ls4ϕ)2]dx−ζ4∫Ω(|φ|4+|ϕ|4+2ε|φ|2|ϕ|2)dx≡H(0), | (1.5) |
where M is the fractional mass, and we refer to H as the fractional Hamiltonian energy.
By setting φ=˜p+i˜q, ϕ=ˆp+iˆq, the original system (1.1), (1.2) can be rewritten as the following real system
˜pt−ϑLs2˜q+ζ[˜p2+˜q2+ε(ˆp2+ˆq2)]˜q=0, | (1.6) |
˜qt+ϑLs2˜p+ζ[˜p2+˜q2+ε(ˆp2+ˆq2)]˜p=0, | (1.7) |
ˆpt−ϑLs2ˆq+ζ[ˆp2+ˆq2+ε(˜p2+˜q2)]ˆq=0, | (1.8) |
ˆqt+ϑLs2ˆp+ζ[ˆp2+ˆq2+ε(˜p2+˜q2)]ˆp=0. | (1.9) |
According to the variational derivative formula [33], the original fractional CNLS system (1.6)–(1.9) can be expressed by an infinite-dimensional Hamiltonian system
dydt=SδH(y)δy, S=(OI2−I20), | (1.10) |
where δH(y)δy is the vector of variational derivatives for y, y=(˜q,ˆq,˜p,ˆp)T, I2 is the second order unit matrix, and
H(y)=ϑ2∫Ω[(Ls4˜p)2+(Ls4˜q)2+(Ls4ˆp)2+(Ls4ˆq)2]dx−ζ4∫Ω[(ˆp2+ˆq2)2+(˜p2+˜q2)2+2ε(ˆp2+ˆq2)(˜p2+˜q2)]dx. | (1.11) |
In recent years, some numerical algorithms [15,23,24,37] have been given for solving the fractional NLS equation as explicit analytical solutions of fractional differential equations can not be obtained [10]. However, these traditional algorithms can not inherit the conservative laws of the equation and can not be implemented in long-time numerical calculations. Inspired by the idea of the structure-preserving algorithm of the classical Hamiltonian system, scholars have developed many conservative schemes for the fractional CNLS system. For example, Wang and Xiao [29,30,31] first derived invariants of system (1.1) and (1.2) and constructed finite difference schemes to conserve these invariants; the existence and uniqueness of solutions and convergence of schemes are discussed in these works. In [22,26], Li et al. developed an efficient conservative difference and finite element scheme for solving the strongly fractional CNLS equation. There are many corresponding structure algorithms for fractional NLS-type equations, and the readers can refer to the literature [6,32,35].
Unfortunately, these numerical schemes mentioned above only have second-order accuracy temporally or require small time step sizes to obtain satisfactory numerical solutions for NLS-type equations. Therefore, it is necessary to construct and analyze high-accuracy conservative numerical schemes for the equation. Over the past decade, some numerical methods have been proposed to construct high-accuracy conservative methods for Hamiltonian systems [1,21]. These methods can also be applied to develop high order conservative schemes for fractional NLS-type equations. But, these schemes can both preserve not the mass and energy. Fortunately, some energy quadratic methods exist, including the invariant energy quadratization (IEQ) [9,36] and the scalar auxiliary variable (SAV) approaches [3,25,27], originally developed for gradient flow models. These methods have been applied to construct high-accuracy modified energy-preserving schemes for conservative equations [7,8,14,19,20] by combining with the symplectic Runge-Kutta method [12]. The quadratic auxiliary variable (QAV) approach is newly proposed to construct high-accuracy schemes for the KdV equation [2], the derived schemes can conserve the original energy of the system. This method provides us with a way to construct some numerical schemes that can preserve the original energy of conservative systems. But, the QAV method is not applicable for some equations with complex energy functionals.
As is shown in Eq (1.5), the Hamiltonian energy functional is very complicated and has high degree polynomials, and cannot be handled by introducing a single auxiliary variable. Therefore, the original QAV approach can not be used to construct numerical schemes to inherit the Hamiltonian energy for the equation. Influenced by the idea of the QAV method, we developed a new method which we called the multiple quadratic auxiliary variables (MQAV) approach by introducing more quadratic auxiliary variables and construct conservative schemes for system (1.1), (1.2) based on the symplectic Runge-Kutta method. The constructed schemes can preserve the mass and Hamiltonian energy of the system and have high-accuracy in the temporal direction. In addition, the MQAV approach also can be applied to other conservative systems with complex energy functionals and high degree polynomials.
The following is the outline of this paper. Section 2 gives an equivalent system for the fractional CNLS equation by introducing two quadratic auxiliary variables. Section 3 applies the symplectic Runge-Kutta method to the equivalent system in the time direction can derive a conservative semi-discrete scheme, and discusses the scheme's conservation properties. Section 4 uses the Fourier pseudo-spectral method to approximate the semi-discrete system to obtain a fully-discrete conservative scheme and gives the proposed schemes' stability. In Section 5, numerical results are displayed to illustrate the theoretical results. The final section contains some conclusions.
We first derive an equivalent system for the fractional CNLS system based on the idea of the MQAV method.
By setting u=|φ|2, v=|ϕ|2, the original Hamiltonian energy can be rewritten as
E(t)=ϑ2∫Ω[(Ls4φ)2+(Ls4ϕ)2]dx−ζ4∫Ω(u2+v2+2εuv)dx. | (2.1) |
Considering the following initial conditions
φ0=φ(x,0), ϕ=ϕ(x,0), u0=|φ0|2, v0=|ϕ0|2, |
and taking the variational derivative for the new energy Eq (2.1), the original fractional CNLS system (1.1), (1.2) can be rewritten as the following equivalent MQAV system
iφt−ϑLs2φ+ζ(u+εv)φ=0, | (2.2) |
iϕt−ϑLs2ϕ+ζ(v+εu)ϕ=0, | (2.3) |
ut=2Re(φt,φ), | (2.4) |
vt=2Re(ϕt,ϕ), | (2.5) |
here (∗,∙):=∫Ω∗ˉ∙dx, ˉ∙ is the conjugate of ∙, Re(∗) represents taking the real part of ∗.
Theorem 2.1. The newly proposed MQAV system can inherit some invariants, namely
(i).I(t)≡I(0); (ii).M(t)≡M(0); (iii).E(t)≡E(0), |
where
I(t):=u−|φ|2+v−|ϕ|2, |
and M(t), E(t) are given in Eqs (1.4) and (2.1), respectively.
Proof. First, system (2.4), (2.5) can be rewritten as
∂t(u−|φ|2)=0, |
∂t(v−|ϕ|2)=0, |
and adding the above two formulas we can derive
ddtI(t)=0. |
Then, we compute the inner products of equation (2.2), (2.3) with φ and ϕ, and take the imaginary part of the derivation, leading to
ddtM(t)=0. |
Further, by computing the inner products of system (2.2)–(2.5) with φt,ϕt,u,v, respectively, we derive
ddtE(t)=0. |
Remark 2.1. It is worth noting that the equivalent MQAV system not only preserves the modified energy and mass but also conserves a new invariant I(t), which implies the modified energy E(t) is equivalent to the original Hamiltonian energy H(t) in the continuous case. For the constructed new system, the auxiliary variables relation can be regarded as weak properties, while the original energy reduces to a weak invariant [12].
We all know that the construction of the high-accuracy conservative schemes for the fractional CNLS equation is challenging. Fortunately, the MQAV system (2.2)–(2.5) can preserve the quadratic conserved quantities. This allows us to use the symplectic RK method that can conserve the system's quadratic invariants to construct a high-order accuracy conservative scheme for the fractional CNLS equation. First, we outline some notations of the symplectic RK method.
We consider the following ordinary differential equation
dw(t)dt=f(t), w(t0)=w0. | (3.1) |
We define τ=TN as the time step, and set tn=nτ, n=0,1,⋯,N. An s-stage RK method for system (3.1) in [tn,tn+1] is given by
win=wn+τs∑j=1aijkj, ki=f(tn+ciτ), i=1,2,⋯,s, | (3.2) |
wn+1=wn+τs∑i=1biki, | (3.3) |
where aij,bi, ci=s∑i=1aij (i,j=1,2,⋯,s) are real coefficients. If these coefficients have the following relationship
aijbi+ajibj−bibj=0, | (3.4) |
the proposed RK method is called as the symplectic RK method [12].
To develop a high-order conservative scheme for the MQAV system that preserves some quadratic invariants in the time direction, a natural idea is to use the symplectic RK method to approximate system (2.2)–(2.5). The resulting semi-discrete RK scheme is given as follows
{φin=φn+τs∑j=1aijkjφ, kiφ=−i(ϑLs2φin−ζ(uin+εvin)φin),ϕin=ϕn+τs∑j=1aijkjϕ, kiϕ=−i(ϑLs2ϕin−ζ(vin+εuin)ϕin),uin=un+τs∑j=1aijkju, kiu=2Re(ˉφinkiφ),vin=vn+τs∑j=1aijkjv, kiv=2Re(ˉϕinkiϕ), | (3.5) |
Then (φn+1,ϕn+1,un+1,vn+1)T can be updated by
{φn+1=φn+τs∑i=1bikiφ,ϕn+1=ϕn+τs∑i=1bikiϕ,un+1=un+τs∑i=1bikiu,vn+1=vn+τs∑i=1bikiv, | (3.6) |
where zn represents the numerical approximation to z(⋅,t) and tn for any function z. For the semi-discrete scheme, we can derive the following theorem for the structure-preserving properties.
Theorem 3.1. The system (3.5), (3.6) can inherit the following invariants, i.e.,
Mn=M0, In=I0, En=E0, |
where Mn is the mass function and is defined by
Mn=(φn,φn)+(ϕn,ϕn), | (3.7) |
and
In=(un,1)+(vn,1)−(φn,φn)−(ϕn,ϕn), | (3.8) |
En is the energy function, which has the form
En=ϑ2[(Ls2φn,φn)+(Ls2ϕn,ϕn)]−ζ4[(un,un)+(vn,vn)+2ε(un,vn)]. | (3.9) |
Proof. According to system (3.6) and the definition of Mn, we have
Mn+1−Mn=(φn+1,φn+1)−(φn,φn)+(ϕn+1,ϕn+1)−(ϕn,ϕn)=τs∑i=1bi(φin,kiφ)+τs∑i=1bi(kiφ,φin)+τs∑i=1bi(ϕin,kiϕ)+τs∑i=1bi(kiϕ,ϕin)=2τs∑i=1biRe(kiφ,φin)+2τs∑i=1biRe(kiϕ,ϕin), | (3.10) |
Plugging kiφ=−i(ϑLs2φin−ζ(uin+εvin)φin), kiϕ=−i(ϑLs2ϕin−ζ(vin+εuin)ϕin) into above system, we derive
Mn+1−Mn=2τs∑i=1biRe(−i(ϑLs2φin−ζ(uin+εvin)φin),φin)+2τs∑i=1biRe(−i(ϑLs2ϕin−ζ(vin+εuin)ϕin),ϕin)=2τs∑i=1biRe(−i(ϑ|Ls4φin|2−ζ(uin+εvin)|φin|2))+2τs∑i=1biRe(−i(ϑ|Ls4ϕin|2−ζ(vin+εuin)|ϕin|2))=0, | (3.11) |
where symplectic condition (3.4) and s∑i,j=1aijbikikj = s∑i,j=1ajibjkikj were used. Noticing that
(un+1,1)−(un,1)+(vn+1,1)−(vn,1)=τs∑i=1bikiu+τs∑i=1bikiv=2τs∑i=1biRe(ˉφinkiφ)+2τs∑i=1biRe(ˉϕinkiϕ)=2τs∑i=1biRe(ˉφinkiφ)+2τs∑i=1biRe(ˉϕinkiϕ)=2τs∑i=1biRe(kiφ,φin)+2τs∑i=1biRe(kiϕ,ϕin)=τs∑i=1bi(φin,kiφ)+τs∑i=1bi(kiφ,φin)+τs∑i=1bi(ϕin,kiϕ)+τs∑i=1bi(kiϕ,ϕin), | (3.12) |
and
(φn+1,φn+1)−(φn,φn)+(ϕn+1,ϕn+1)−(ϕn,ϕn)=τs∑i=1bi(φin,kiφ)+τs∑i=1bi(kiφ,φin)+τs∑i=1bi(ϕin,kiϕ)+τs∑i=1bi(kiϕ,ϕin)=2τs∑i=1biRe(kiφ,φin)+2τs∑i=1biRe(kiϕ,ϕin). | (3.13) |
This together with Eq (3.12) can derive
In+1−In=0. |
Then, we further obtain
(Ls2φn+1,φn+1)−(Ls2φn,φn)=(Ls2φn+Ls2τs∑i=1bikiφ,φn+τs∑i=1bikiφ)−(Ls2φn,φn)=(Ls2φn,τs∑i=1bikiφ)+(Ls2τs∑i=1bikiφ,φn)+(Ls2τs∑i=1bikiφ,τs∑i=1bikiφ)=2τs∑i=1Re(Ls2φn,bikiφ)+τ2s∑i,j=1(Ls2bikiφ,bjkjφ)=2τs∑i=1Re(Ls2φin−τs∑j=1aijkjφ,bikiφ)+τ2s∑i,j=1(Ls2bikiφ,bjkjφ)=2τs∑i=1Re(Ls2φin,bikiφ)−2τ2s∑i,j=1(Ls2aijkjp,bikiφ)+τ2s∑i,j=1(Ls2bikiφ,bjkjφ)=2τs∑i=1Re(Ls2φin,bikiφ), | (3.14) |
where the symplectic condition (3.4) and s∑i,j=1aijbikikj = s∑i,j=1ajibjkikj were used. Furthermore, similar discussions, we obtain
(Ls2ϕn+1,ϕn+1)−(Ls2ϕn,ϕn)=2τs∑i=1Re(Ls2ϕin,bikiϕ), | (3.15) |
(un+1,un+1)−(un,un)=2τs∑i=1(uin,bikiu), | (3.16) |
(vn+1,vn+1)−(vn,vn)=2τs∑i=1(vin,bikiv). | (3.17) |
Therefore, we derive
(Ls2φn+1,φn+1)−(Ls2φn,φn)+(Ls2ϕn+1,ϕn+1)−(Ls2ϕn,ϕn)=2τs∑i=1Re(Ls2φin,bikiφ)+2τs∑i=1Re(Ls2ϕin,bikiϕ), | (3.18) |
(un+1,un+1)−(un,un)+(vn+1,vn+1)−(vn,vn)=2τs∑i=1(uin,bikiu)+2τs∑i=1(vin,bikiv), | (3.19) |
(un+1,vn+1)−(un,vn)=τs∑i=1(uin,bikiv)+τs∑i=1(vin,bikiu), | (3.20) |
Based on the definition of En, we can deduce
En+1−En=ϑ2[(Ls2φn+1,φn+1)−(Ls2φn,φn)+(Ls2ϕn+1,ϕn+1)−(Ls2ϕn,ϕn)]−ζ4[(un+1,un+1)−(un,un)+(vn+1,vn+1)−(vn,vn)+ε(un+1,vn+1)−ε(un,vn)]=τs∑i=1biRe[(ϑLs2φin−ζφinuin−ζεφinvin,kiφ)+(ϑLs2ϕin−ζϕinvin−ζεϕinuin,kiϕ)]=τs∑i=1biRe[(−ikiφ,kiφ)+(−ikiϕ,kiϕ)]=0. | (3.21) |
The proof is completed.
Then, under the consistent initial condition u0=|φ0|, v0=|ϕ0|, we have the following original energy conservation theorem for the semi-discrete MQAV-RK system.
Theorem 3.2. The MQAV-RK system (3.5), (3.6) can also inherit the original semi-discrete Hamiltonian energy of the fractional CNLS system, namely
Hn=H0, |
where Hn is the original semi-discrete energy function and given by
Hn=ϑ2[(Ls2φn,φn)+(Ls2ϕn,ϕn)]−ζ4[((φn)2,(φn)2)+((ϕn)2,(ϕn)2)+2ε((φn)2,(ϕn)2)]. | (3.22) |
Proof. Based on Theorem 3.1 and the two introduced quadratic auxiliary variable u=|φ|2, v=|ϕ|2, one can obtain un=|φn|2, vn=|ϕn|2 for 0≤n≤N, which leads to
En=ϑ2[(Ls2φn,φn)+(Ls2ϕn,ϕn)]−ζ4[((φn)2,(φn)2)+((ϕn)2,(ϕn)2)+2ε((φn)2,(ϕn)2)]=Hn, | (3.23) |
and we can immediately obtain Hn=H0.
The proposed conservative scheme has fourth-order accuracy in time, and the Fourier pseudo-spectral method is chosen to approximate the semi-discrete system in space so that the accuracy of time direction is not affected. We consider system (1.1), (1.2) in a periodic region Ω=[−L,L], and define xj=−L+jh=−L+j2LN,0≤j≤N−1, N is a positive even integer, and h:=2LN. Furthermore, we denote
Vh={v|v=(v0,v1,⋯,vN−1)}, |
be a vector space of grid functions with vj=v(xj). Then defining the interpolation space SpN:=span{gj(x)} with
gj(x)=1NN/2∑l=−N/21aleilμ(x−xj), |
where μ=πL, al={1,|l|<N/2,2,|l|=N/2,. Then, we set IN:L2(Ω)→SpN as the interpolation operator
INv(x)=N−1∑j=0vjgj(x)=N/2∑l=−N/2ˆvleilμ(x+L), |
with
ˆvl=1NalN−1∑j=0v(xj)e−ilμ(xj+L), | (4.1) |
where ˆv−N/2=ˆvN/2. Applying the Laplacian operator Ls2 to INv(x) leads to
Ls2INv(xj)=Ds2v:=N/2∑l=−N/2|μl|sˆvleilμ(xj+L)=(F−1NΛFN)v, | (4.2) |
where v∈Vh, and
Λ=diag[0s,|μ|s,⋯,|μN2|s,|μ(−N2+1)|s,⋯,|−2μ|s,|−μ|s]. |
and the discrete Fourier transform and its inverse are given by
(FNv)k=1NN−1∑j=0vje−2πijkN,(F−1Nˆv)k=N/2−1∑k=−N/2ˆvke2πijkN. | (4.3) |
It is can be computed by using the FFT technique in long time simulations.
This subsection aims to give a high-order fully-discrete scheme for the Eqs (1.1), (1.2), and the scheme is derived by using the Fourier pseudo-spectral method to approximate equation (3.5), (3.6) in space, namely
{φin=φn+τs∑j=1aijkjφ, kiφ=−i(ϑDs2φin−ζ(uin+εvin)φin),ϕin=ϕn+τs∑j=1aijkjϕ, kiϕ=−i(ϑDs2ϕin−ζ(vin+εuin)ϕin),uin=un+τs∑j=1aijkju, kiu=2Re(ˉφinkiφ),vin=vn+τs∑j=1aijkjv, kiv=2Re(ˉϕinkiϕ). | (4.4) |
Then (φn+1,ϕn+1,un+1,vn+1)T can be updated by
{φn+1=φn+τs∑i=1bikiφ,ϕn+1=ϕn+τs∑i=1bikiϕ,un+1=un+τs∑i=1bikiu,vn+1=vn+τs∑i=1bikiv. | (4.5) |
For a fully-discrete scheme, we can derive the following theorem for the structure-preserving properties.
Theorem 4.1. System (4.4), (4.5) can conserve some conservation laws, i.e.,
Mn=M0, In=I0, En=E0, |
where Mn is the mass function and defined by
Mn=⟨φn,φn⟩+⟨ϕn,ϕn⟩, | (4.6) |
and
In=⟨un,1⟩+⟨vn,1⟩−⟨φn,φn⟩−⟨ϕn,ϕn⟩, | (4.7) |
and En is the energy function, which has the form
En=ϑ2[⟨Ds2φn,φn⟩+⟨Ds2ϕn,ϕn⟩]−ζ4[⟨un,un⟩+⟨vn,vn⟩+2ε⟨un,vn⟩], | (4.8) |
where ⟨⋆,∙⟩=hN−1∑j=0⋆jˉ∙j.
Proof. The proof process is the same as the Theorem 3.1.
In addition, we can also prove the developed schemes can inherit the original Hamiltonian energy.
Theorem 4.2. System (4.4), (4.5) can conserve the original Hamiltonian energy of the fractional CNLS system, namely
Hn=H0, |
where
Hn=ϑ2[⟨Ds2φn,φn⟩+⟨Ds2ϕn,ϕn⟩]−ζ4[⟨|φn|2,|φn|2⟩+⟨|ϕn|2,|ϕn|2⟩+2ε⟨|φn|2,|ϕn|2⟩]. | (4.9) |
Remark 4.1. Actually, the QAV method is similar to the IEQ and SAV methods, both of which are energy quadratic methods. Therefore, they can also be used to develop high-order conservative schemes for conservative equations. However, the derived schemes based on the IEQ and SAV methods can only inherit the modified energy [5,8]. Theorem 4.1 and Theorem 4.2 indicate that the constructed high-order MQAV-RK scheme can inherit both the modified energy and the Hamiltonian energy.
According to the mass conservation law of Theorem 4.1, we can conclude that the proposed schemes are stable.
We give some examples to confirm the conservation properties and the high-accuracy of the developed schemes. The exact solution of the system is not given, and thus we compute the numerical errors by using the formula
error=‖ϕN(h,τ)−ϕ2N(h,τ/2)‖∞. | (5.1) |
The accuracy of the constructed schemes can be computed by
Rate=ln(error1/error2)/ln(τ1/τ2), | (5.2) |
where τj, errorj,(j=1,2) are the time step and the maximum-norm errors with τj, respectively. Furthermore, to test the conservation properties, we define the relative errors of some invariants as
RMn=|Mn−M0M0|, RHn=|Hn−H0H0|, REn=|En−E0E0|, | (5.3) |
where Mn, Hn and En denote the invariant at t=nτ. Without losing generality, the paper will consider the following 2s-4th order symplectic RK method [12]
To compare the advantages of the constructed scheme with some conservative schemes, we define
● MQAV-RK4: The constructed fourth-order conservative MQAV-RK scheme in this paper;
● SAV-RK4: The proposed fourth-order conservative RK scheme based on the SAV method [8].
● CN2: The proposed second-order Crank-Nicolson conservative scheme for the Eqs (1.1) and (1.2) [30].
Example 5.1. First, consider a decoupled fractional NLS equation in Ω=(−π,π) with the form
iφt−ϑLs2φ+ζ|φ|2φ=0, |
with the exact solution
φ(x,t)=ηexp(i(γx−δt)), δ=|γ|s−θ|η|2. |
In the simulation, we take η=2, γ=1, θ=2 to exhibit the theoretical results.
First, we set s=1.6 and T=1 to test the convergence of the given three numerical schemes. The maximum-norm errors and convergence rates are shown in Table 2. The listed data in the table indicate that the MQAV-RK4 and SAV-RK4 both have fourth-order convergence rates in time, and the CN2 scheme only gets second-order accuracy. Furthermore, we investigate the maximum-norm errors of the numerical solution versus the CPU time using the three selected schemes. The results are summarized in Figure 1. For a given global error, we observe that (ⅰ) the cost of the CN2 scheme is the most expensive because of the low-order accuracy in time; (ⅱ) the cost of SAV-RK4 is much cheaper than the CN2 scheme; (ⅲ) the cost of our proposed MQAV-RK4 scheme is the cheapest. We test the conservative properties of four conservative schemes at T=100 by using the denoted formulas in (5.3). Figure 2 reveals relative errors of conservation laws. As is shown above, the SAV-RK4 scheme can only conserve the modified energy and mass of the equation. But, our proposed methods can both inherit the original energy and mass.
12−√36 | 14 | 14−√36 |
12+√36 | 14+√36 | 14 |
12 | 12 |
Scheme | τ=150 | τ=1100 | τ=1200 | τ=1400 | |
CN2[30] | error | 5.78e-05 | 1.42e-05 | 3.54e-06 | 8.82e-07 |
Rate | * | 2.01 | 2.01 | 2.00 | |
SAV-RK4 [8] | error | 9.00e-07 | 7.14e-08 | 4.98e-09 | 3.28e-10 |
Rate | * | 3.65 | 3.84 | 3.92 | |
MQAV-RK4 | error | 5.31e-07 | 3.32e-08 | 2.07e-09 | 1.29e-10 |
Rate | * | 3.99 | 3.99 | 4.00 |
Example 5.2. This example considers the fractional CNLS equation with initial conditions [30]
u(x,0)=sech(x+δ0)⋅exp(iω0x), |
v(x,0)=sech(x−δ0)⋅exp(−iω0x). |
In this example, x∈(−20,20), and we choose δ0=5,ω0=3,ϑ=1,ζ=2, ε=1.
First, Eqs (5.1) and (5.2) are used to compute the L∞−normal and corresponding convergence rate of time for some conservative schemes at T=1. The obtained numerical results are given in Table 3. From Table 3, we can find the MQAV-RK4 scheme has fourth-order convergence. By using the denoted formulas in (5.3), we test the conservative properties of four conservative schemes at T=50. Figure 3 reveals relative errors of conservation laws. As is shown above, the proposed method can both inherit the mass and original energy exactly. Finally, we define Ψ=|φn|+|ϕn| and study the relationship between the evolution of numerical solution and the fractional s. As is shown in Figure 4, the s affects the evolution of solitons, and the shape of the solitons change quickly as s changes.
s | τ=120 | τ=140 | τ=180 | τ=1160 | |
1.4 | error | 1.19e-04 | 7.62e-06 | 4.79e-07 | 2.99e-08 |
Rate | * | 3.97 | 3.99 | 3.99 | |
1.6 | error | 3.38e-04 | 2.17e-05 | 1.36e-06 | 8.56e-08 |
Rate | * | 3.96 | 3.98 | 3.99 | |
2 | error | 2.90e-03 | 1.83e-04 | 1.18e-05 | 7.46e-07 |
Rate | * | 3.98 | 3.95 | 3.98 |
Example 5.3. In this example, we choose the initial condition
φ(x,y,0)=sech(0.5x+y)exp(i(−0.5x−y)), |
ϕ(x,y,0)=sech(0.5x+y)exp(i(−0.5x−y)), |
where (x,y)∈[−8,8]2, and, some parameters are given as follows
ϑ=12,ζ=1,ε=1. |
In Table 4, we give the accuracy of the proposed scheme at s=1.3,1.7,2. As is shown in the table, the scheme has fourth-order accuracy in the time direction for different s. In addition, the numerical errors also increase with the increase of s. Figure 5 implies that the scheme can inherit the original energy and mass of the system for different s.
s | τ=110 | τ=120 | τ=140 | τ=180 | |
s=1.3 | error | 5.22e-07 | 3.26e-08 | 2.04e-09 | 1.27e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=1.7 | error | 3.33e-06 | 2.08e-07 | 1.30e-08 | 8.14e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=2 | error | 1.01e-05 | 6.36e-07 | 3.97e-08 | 2.48e-09 |
Rate | * | 3.99 | 3.99 | 3.99 |
In this paper, we first propose the Hamiltonian structure of the fractional CNLS equation, and then develop the MQAV approach for the fractional CNLS equation with a complex energy function that can not be handled by using the original QAV approach. Then, a class of MQAV-RK schemes is proposed by combining with the symplectic Runge-Kutta method. The constructed schemes have high-accuracy and inherit the mass and Hamiltonian energy. Numerical results also demonstrate that the schemes have good numerical stability. In addition, the approaches presented in this work can be used to solve other conservative differential equations.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The research is supported by the Science Foundation of Henan (No. 212300410323).
The authors declare there is no conflict of interest.
[1] |
Harapan H, Itoh N, Yufika A, et al. (2020) Coronavirus disease 2019 (COVID-19): A literature review. J Infect Public Health 13: 667-673. doi: 10.1016/j.jiph.2020.03.019
![]() |
[2] |
Jiang F, Deng L, Zhang L, et al. (2020) Review of the clinical characteristics of coronavirus disease 2019 (COVID-19). J Gen Intern Med 35: 1545-1549. doi: 10.1007/s11606-020-05762-w
![]() |
[3] |
Shi Y, Wang G, Cai XP, et al. (2020) An overview of COVID-19. J Zhejiang Univ Sci B 21: 343-360. doi: 10.1631/jzus.B2000083
![]() |
[4] | Wang L, Wang Y, Ye D, et al. (2020) Review of the 2019 novel coronavirus (SARS-CoV-2) based on current evidence. Int J Antimicrob Agents 55: 105948. |
[5] |
Ahn DG, Shin HJ, Kim MH, et al. (2020) Current status of epidemiology, diagnosis, therapeutics, and vaccines for novel coronavirus disease 2019 (COVID-19). J Microbiol Biotechnol 30: 313-324. doi: 10.4014/jmb.2003.03011
![]() |
[6] |
Sohrabi C, Alsafi Z, O'Neill N, et al. (2020) World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19). Int J Surg 76: 71-76. doi: 10.1016/j.ijsu.2020.02.034
![]() |
[7] | World Health Organization, Coronavirus Disease 2019 (COVID-19), Situation Report–83, 2000. Available from: https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200412-sitrep-83-covid-19.pdf?sfvrsn=697ce98d_4. |
[8] |
Da Silva C S, Monteiro CRA, Da Silva GHF, et al. (2020) Assessing the metabolic impact of ground chia seed in overweight and obese prepubescent children: results of a double-blind randomized clinical trial. J Med Food 23: 224-232. doi: 10.1089/jmf.2019.0055
![]() |
[9] |
Zhai P, Ding Y, Wu X, et al. (2020) The epidemiology, diagnosis and treatment of COVID-19. Int J Antimicrob Agents 55: 105955. doi: 10.1016/j.ijantimicag.2020.105955
![]() |
[10] |
Tian S, Hu N, Lou J, et al. (2020) Characteristics of COVID-19 infection in Beijing. J Infect 80: 401-406. doi: 10.1016/j.jinf.2020.02.018
![]() |
[11] |
Ge HP, Wang XF, Yuan XN, et al. (2020) The epidemiology and clinical information about COVID-19. Eur J Clin Microbiol Infect Dis 39: 1011-1019. doi: 10.1007/s10096-020-03874-z
![]() |
[12] |
Madabhavi I, Sarkar M, Kadakol N (2020) COVID-19: a review. doi: 10.4081/monaldi.2020.1298
![]() |
[13] |
Tufan A, Avanoğlu Güler A, Matucci-Cerinic M (2020) COVID-19, immune system response, hyperinflammation and repurposing antirheumatic drugs. Turk J Med Sci 50: 620-632. doi: 10.3906/sag-2004-168
![]() |
[14] |
Velavan TP, Meyer CG (2020) The COVID-19 epidemic. Trop Med Int Health 25: 278-280. doi: 10.1111/tmi.13383
![]() |
[15] |
Muniyappa R, Gubbi S (2020) COVID-19 pandemic, coronaviruses, and diabetes mellitus. Am J Physiol Endocrinol Metab 318: E736-E741. doi: 10.1152/ajpendo.00124.2020
![]() |
[16] |
Schiffrin EL, Flack JM, Ito S, et al. (2020) Hypertension and COVID-19. Am J Hypertens 33: 373-374. doi: 10.1093/ajh/hpaa057
![]() |
[17] |
Afshin A, Forouzanfar MH, Reitsma MB, et al. (2017) Health effects of overweight and obesity in 195 countries over 25 years. N Engl J Med 377: 13-27. doi: 10.1056/NEJMoa1614362
![]() |
[18] |
Ogurtsova K, Da Rocha Fernandes JD, Huang Y, et al. (2017) IDF Diabetes Atlas: Global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract 128: 40-50. doi: 10.1016/j.diabres.2017.03.024
![]() |
[19] |
Saklayen MG (2018) The global epidemic of the metabolic syndrome. Curr Hypertens Rep 20: 12. doi: 10.1007/s11906-018-0812-z
![]() |
[20] | Jiang X, Coffee M, Bari A, et al. (2020) Towards an artificial intelligence framework for data-driven prediction of coronavirus clinical severity. Comput Mater Con 63: 537-551. |
[21] |
McCracken E, Monaghan M, Sreenivasan S (2018) Pathophysiology of the metabolic syndrome. Clin Dermatol 36: 14-20. doi: 10.1016/j.clindermatol.2017.09.004
![]() |
[22] |
Barbalho SM, Tofano RJ, De Campos AL, et al. (2018) Association between vitamin D status and metabolic syndrome risk factors. Diabetes Metab Syndr 12: 501-507. doi: 10.1016/j.dsx.2018.03.011
![]() |
[23] |
Marhl M, Grubelnik V, Magdič M, et al. (2020) Diabetes and metabolic syndrome as risk factors for COVID-19. Diabetes Metab Syndr 14: 671-677. doi: 10.1016/j.dsx.2020.05.013
![]() |
[24] |
Singhal T (2020) A review of coronavirus disease-2019 (COVID-19). Indian J Pediatr 87: 281-286. doi: 10.1007/s12098-020-03263-6
![]() |
[25] |
Soy M, Keser G, Atagündüz P, et al. (2020) Cytokine storm in COVID-19: pathogenesis and overview of anti-inflammatory agents used in treatment. Clin Rheumatol 39: 2085-2094. doi: 10.1007/s10067-020-05190-5
![]() |
[26] |
Ye Q, Wang B, Mao J (2020) The pathogenesis and treatment of the ‘Cytokine Storm’ in COVID-19. J Infect 80: 607-613. doi: 10.1016/j.jinf.2020.03.037
![]() |
[27] |
Butler MJ, Barrientos RM (2020) The impact of nutrition on COVID-19 susceptibility and long-term consequences. Brain Behav Immun 87: 53-54. doi: 10.1016/j.bbi.2020.04.040
![]() |
[28] |
Petrakis D, Margină D, Tsarouhas K, et al. (2020) Obesity—a risk factor for increased COVID-19 prevalence, severity and lethality. Mol Med Rep 22: 9-19. doi: 10.3892/mmr.2020.11127
![]() |
[29] |
Tadic M, Cuspidi C, Sala C (2020) COVID-19 and diabetes: Is there enough evidence?. doi: 10.1111/jch.13912
![]() |
[30] |
Engin AB, Engin ED, Engin A (2020) Two important controversial risk factors in SARS-CoV-2 infection: obesity and smoking. Environ Toxicol Pharmacol 78: 103411. doi: 10.1016/j.etap.2020.103411
![]() |
[31] |
Ryan PM, Caplice NM (2020) Is adipose tissue a reservoir for viral spread, immune activation, and cytokine amplification in coronavirus disease 2019?. doi: 10.1002/oby.22843
![]() |
[32] |
Zabetakis I, Lordan R, Norton C, et al. (2020) COVID-19: The inflammation link and the role of nutrition in potential mitigation. Nutrients 12: 1466. doi: 10.3390/nu12051466
![]() |
[33] |
Dhar D, Mohanty A (2020) Gut microbiota and Covid-19- possible link and implications. Virus Res 285: 198018. doi: 10.1016/j.virusres.2020.198018
![]() |
[34] |
Araújo JR, Tomas J, Brenner C, et al. (2017) Impact of high-fat diet on the intestinal microbiota and small intestinal physiology before and after the onset of obesity. Biochimie 141: 97-106. doi: 10.1016/j.biochi.2017.05.019
![]() |
[35] |
Cox AJ, West NP, Cripps AW (2015) Obesity, inflammation, and the gut microbiota. Lancet Diabetes Endocrinol 3: 207-215. doi: 10.1016/S2213-8587(14)70134-2
![]() |
[36] |
Gupta R, Hussain A, Misra A (2020) Diabetes and COVID-19: evidence, current status and unanswered research questions. Eur J Clin Nutr 74: 864-870. doi: 10.1038/s41430-020-0652-1
![]() |
[37] |
Hussain A, Bhowmik B, Do Vale Moreira NC (2020) COVID-19 and diabetes: Knowledge in progress. Diabetes Res Clin Pract 162: 108142. doi: 10.1016/j.diabres.2020.108142
![]() |
[38] |
Schofield J, Leelarathna L, Thabit H (2020) COVID-19: Impact of and on diabetes. Diabetes Ther 11: 1429-1435. doi: 10.1007/s13300-020-00847-5
![]() |
[39] |
Cristelo C, Azevedo C, Marques JM, et al. (2020) SARS-CoV-2 and diabetes: New challenges for the disease. Diabetes Res Clin Pract 164: 108228. doi: 10.1016/j.diabres.2020.108228
![]() |
[40] |
Gupta R, Ghosh A, Singh AK, et al. (2020) Clinical considerations for patients with diabetes in times of COVID-19 epidemic. Diabetes Metab Syndr 14: 211-212. doi: 10.1016/j.dsx.2020.03.002
![]() |
[41] |
Brufsky A (2020) Hyperglycemia, hydroxychloroquine, and the COVID-19 pandemic. J Med Virol 92: 770-775. doi: 10.1002/jmv.25887
![]() |
[42] |
Drucker DJ (2020) Coronavirus infections and type 2 diabetes—shared pathways with therapeutic implications. Endocr Rev 41: 457-470. doi: 10.1210/endrev/bnaa011
![]() |
[43] |
Orioli L, Hermans MP, Thissen JP, et al. (2020) COVID-19 in diabetic patients: Related risks and specifics of management. Ann Endocrinol 81: 101-109. doi: 10.1016/j.ando.2020.05.001
![]() |
[44] |
Sardu C, Gambardella J, Morelli MB, et al. (2020) Hypertension, thrombosis, kidney failure, and diabetes: Is COVID-19 an endothelial disease? A comprehensive evaluation of clinical and basic evidence. J Clin Med 9: 1417. doi: 10.3390/jcm9051417
![]() |
[45] |
Zheng Z, Peng F, Xu B, et al. (2020) Risk factors of critical & mortal COVID-19 cases: A systematic literature review and meta-analysis. doi: 10.1016/j.jinf.2020.04.021
![]() |
[46] |
Bouhanick B, Cracowski JL, Faillie JL (2020) Diabetes and COVID-19. doi: 10.1016/j.therap.2020.05.006
![]() |
[47] |
Zheng YY, Ma YT, Zhang JY, et al. (2020) COVID-19 and the cardiovascular system. Nat Rev Cardiol 17: 259-260. doi: 10.1038/s41569-020-0360-5
![]() |
[48] |
Li B, Yang J, Zhao F, et al. (2020) Prevalence and impact of cardiovascular metabolic diseases on COVID-19 in China. Clin Res Cardiol 109: 531-538. doi: 10.1007/s00392-020-01626-9
![]() |
[49] |
Ke C, Zhu X, Zhang Y, et al. (2018) Metabolomic characterization of hypertension and dyslipidemia. Metabolomics 14: 117. doi: 10.1007/s11306-018-1408-y
![]() |
[50] |
Chobanian AV (2017) Guidelines for the management of hypertension. Med Clin N Am 101: 219-227. doi: 10.1016/j.mcna.2016.08.016
![]() |
[51] |
Gupta R, Misra A (2020) Contentious issues and evolving concepts in the clinical presentation and management of patients with COVID-19 infectionwith reference to use of therapeutic and other drugs used in Co-morbid diseases (Hypertension, diabetes etc). Diabetes Metab Syndr 14: 251-254. doi: 10.1016/j.dsx.2020.03.012
![]() |
[52] |
Cheng H, Wang Y, Wang GQ (2020) Organ-protective effect of angiotensin-converting enzyme 2 and its effect on the prognosis of COVID-19. J Med Virol 92: 726-730. doi: 10.1002/jmv.25785
![]() |
[53] |
South AM, Diz DI, Chappell MC (2020) COVID-19, ACE2, and the cardiovascular consequences. Am J Physiol Heart Circ Physiol 318: H1084-H1090. doi: 10.1152/ajpheart.00217.2020
![]() |
[54] |
D'Ardes D, Boccatonda A, Rossi I, et al. (2020) COVID-19 and RAS: unravelling an unclear relationship. Int J Mol Sci 21: 3003. doi: 10.3390/ijms21083003
![]() |
[55] | Angel-Korman A, Brosh T, Glick K, et al. (2020) COVID-19, the kidney and hypertension. Harefuah 159: 231-234. |
[56] |
Wang D, Hu B, Hu C, et al. (2020) Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus–infected pneumonia in Wuhan, China. Jama 323: 1061-1069. doi: 10.1001/jama.2020.1585
![]() |
[57] |
Meng J, Xiao G, Zhang J, et al. (2020) Renin-angiotensin system inhibitors improve the clinical outcomes of COVID-19 patients with hypertension. Emerg Microbes Infect 9: 757-760. doi: 10.1080/22221751.2020.1746200
![]() |
[58] |
Cao X, Yin R, Albrecht H, et al. (2020) Cholesterol: A new game player accelerating endothelial injuries caused by SARS-CoV-2? Am J Physiol Endocrinol Metab 319: E197-E202. doi: 10.1152/ajpendo.00255.2020
![]() |
[59] |
Silva AME, Aguiar C, Duarte J S, et al. (2019) CODAP: A multidisciplinary consensus among Portuguese experts on the definition, detection and management of atherogenic dyslipidemia. Rev Port Cardiol 38: 531-542. doi: 10.1016/j.repc.2019.03.005
![]() |
[60] |
Helkin A, Stein JJ, Lin S, et al. (2016) Dyslipidemia part 1—review of lipid metabolism and vascular cell physiology. Vasc Endovascular Surg 50: 107-118. doi: 10.1177/1538574416628654
![]() |
[61] |
Zheng KI, Gao F, Wang XB, et al. (2020) Letter to the Editor: Obesity as a risk factor for greater severity of COVID-19 in patients with metabolic associated fatty liver disease. Metabolism 108: 154244. doi: 10.1016/j.metabol.2020.154244
![]() |
[62] |
Vekic J, Zeljkovic A, Stefanovic A, et al. (2019) Obesity and dyslipidemia. Metabolism 92: 71-81. doi: 10.1016/j.metabol.2018.11.005
![]() |
[63] |
Wei X, Zeng W, Su J, et al. (2020) Hypolipidemia is associated with the severity of COVID-19. J Clin Lipidol 14: 297-304. doi: 10.1016/j.jacl.2020.04.008
![]() |
[64] |
Vuorio A, Watts GF, Kovanen PT (2020) Familial hypercholesterolaemia and COVID-19: triggering of increased sustained cardiovascular risk. J Intern Med 287: 746-747. doi: 10.1111/joim.13070
![]() |
12−√36 | 14 | 14−√36 |
12+√36 | 14+√36 | 14 |
12 | 12 |
Scheme | τ=150 | τ=1100 | τ=1200 | τ=1400 | |
CN2[30] | error | 5.78e-05 | 1.42e-05 | 3.54e-06 | 8.82e-07 |
Rate | * | 2.01 | 2.01 | 2.00 | |
SAV-RK4 [8] | error | 9.00e-07 | 7.14e-08 | 4.98e-09 | 3.28e-10 |
Rate | * | 3.65 | 3.84 | 3.92 | |
MQAV-RK4 | error | 5.31e-07 | 3.32e-08 | 2.07e-09 | 1.29e-10 |
Rate | * | 3.99 | 3.99 | 4.00 |
s | τ=120 | τ=140 | τ=180 | τ=1160 | |
1.4 | error | 1.19e-04 | 7.62e-06 | 4.79e-07 | 2.99e-08 |
Rate | * | 3.97 | 3.99 | 3.99 | |
1.6 | error | 3.38e-04 | 2.17e-05 | 1.36e-06 | 8.56e-08 |
Rate | * | 3.96 | 3.98 | 3.99 | |
2 | error | 2.90e-03 | 1.83e-04 | 1.18e-05 | 7.46e-07 |
Rate | * | 3.98 | 3.95 | 3.98 |
s | τ=110 | τ=120 | τ=140 | τ=180 | |
s=1.3 | error | 5.22e-07 | 3.26e-08 | 2.04e-09 | 1.27e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=1.7 | error | 3.33e-06 | 2.08e-07 | 1.30e-08 | 8.14e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=2 | error | 1.01e-05 | 6.36e-07 | 3.97e-08 | 2.48e-09 |
Rate | * | 3.99 | 3.99 | 3.99 |
12−√36 | 14 | 14−√36 |
12+√36 | 14+√36 | 14 |
12 | 12 |
Scheme | τ=150 | τ=1100 | τ=1200 | τ=1400 | |
CN2[30] | error | 5.78e-05 | 1.42e-05 | 3.54e-06 | 8.82e-07 |
Rate | * | 2.01 | 2.01 | 2.00 | |
SAV-RK4 [8] | error | 9.00e-07 | 7.14e-08 | 4.98e-09 | 3.28e-10 |
Rate | * | 3.65 | 3.84 | 3.92 | |
MQAV-RK4 | error | 5.31e-07 | 3.32e-08 | 2.07e-09 | 1.29e-10 |
Rate | * | 3.99 | 3.99 | 4.00 |
s | τ=120 | τ=140 | τ=180 | τ=1160 | |
1.4 | error | 1.19e-04 | 7.62e-06 | 4.79e-07 | 2.99e-08 |
Rate | * | 3.97 | 3.99 | 3.99 | |
1.6 | error | 3.38e-04 | 2.17e-05 | 1.36e-06 | 8.56e-08 |
Rate | * | 3.96 | 3.98 | 3.99 | |
2 | error | 2.90e-03 | 1.83e-04 | 1.18e-05 | 7.46e-07 |
Rate | * | 3.98 | 3.95 | 3.98 |
s | τ=110 | τ=120 | τ=140 | τ=180 | |
s=1.3 | error | 5.22e-07 | 3.26e-08 | 2.04e-09 | 1.27e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=1.7 | error | 3.33e-06 | 2.08e-07 | 1.30e-08 | 8.14e-10 |
Rate | * | 3.99 | 3.99 | 4.00 | |
s=2 | error | 1.01e-05 | 6.36e-07 | 3.97e-08 | 2.48e-09 |
Rate | * | 3.99 | 3.99 | 3.99 |