1.
Introduction
This paper focuses on the 3D subdiffusion equations with variable coefficients
in which: a1, a2, a3, and b are positive constants; f, φ, and g are smooth functions; and Λ=(0,Υx)×(0,Υy)×(0,Υz) is a rectangular field bounded by ∂Λ. The Caputo fractional derivative can be defined via
The fractional subdiffusion equation with variable coefficients is derived from the classical partial differential equation. Many scholars have studied the existence and uniqueness of the solution of this equation [1,2]. The high-dimensional subdiffusion problem can accurately describe the mechanical and physical processes with non-locality and spatial global correlation, it has become an important tool to describe a variety of complex mechanical and physical behaviors [3,4,5], and many numerical schemes have been proposed in [6,7,8,9,10,11,12,13,14,15]. Among them, the ADI method is favored by many scholars. Tian and Ge [16] offered the fourth-order ADI difference method for 2D unsteady convection diffusion problems. Wang et al. [17,18] considered the 2D fractional integro-differential equation by using an ADI compact difference method. Zhang et al. [19] provided the 2D semilinear multidelay parabolic equations by using Crank-Nicolson ADI compact difference schemes. Zhang et al. [20] considered the 2D space-fractional nonlinear Ginzburg-Landau equation via linearized ADI difference schemes. Wang et al. [21,22] presented an analysis of the convergence of the α-robust H1-norm for the ADI differece scheme for the 2D time fractional diffusion equation (TFDE). Chen et al. [23,24] considered the ADI difference scheme for a class of two-dimensional multi-term TFDEs and gave the point error estimation. Qiao et al. [25,26,27] presented the ADI difference scheme for 2D TFDE with a weakly singular kernel. Zhang et al. [28] presented an efficient ADI difference scheme for the nonlocal evolution problem in 3D space. Zhou et al. [29,30] solved the 3D TFDE via the ADI method.
Along with the previous research, many scholars considered the subdiffusion equation with variable coefficients. Ngondiep [31] discussed a class of convection-diffusion-reaction equations with variable coefficients and gave a two-level fourth-order scheme. Wei et al. [32] gave the anisotropic finite element method (FEM) for 2D time fractional variable coefficient diffusion equations on graded meshes. Zeng and Tan [33] proposed a two-grid FEM for nonlinear TFDEs with variable coefficients. Ma et al. [34] presented the L1 robust analysis of the fourth-order block-centered FDM for 2D TFDEs with variable coefficient. As we all know, the study of the variable coefficient subdiffusion equation is almost based on one or two dimensions. It is not easy to solve three-dimensional fractional subdiffusion equations with variable coefficients due to the existence of the weak singularity near the initial time. In this paper, we put forward an α-robust L1 ADI difference scheme on graded mesh to overcome these difficulties. The L1 scheme on graded meshes is used to discretize the Caputo fractional derivative, and the difference method is used to approximate the spatial derivative. Using the modified Grönwall inequality, we prove the stability and α-robust convergence of the proposed method.
In this article, the main achievements are as follows:
● We construct an ADI difference method for the 3D subdiffusion equation with variable coefficients. Most previous work has focused on 1D or 2D problems.
● We introduce a new norm to prove the stability and convergence of the proposed method.
● The α-robust convergence is proved strictly, which is able to avoid the blow-up phenomenon in the case of α→1−.
The structure of this paper is as follows. In Section 2, the ADI scheme is developed for the 3D variable coefficient subdiffusion problems. In Section 3, we demonstrate the stability and convergence of the method. In Section 4, the three-dimensional numerical example is given to verify the effectiveness and accuracy of the proposed method. We conclude the paper in Section 5.
2.
Derivation of ADI scheme
Set hx=Υxℓ1, hy=Υyℓ2, hz=Υzℓ3 for the positive integers ℓ1, ℓ2, and ℓ3. Let N≥1 be a positive integer number and set tn=T(n/N)γ for n=0,1,⋯,N, γ≥1. Denote time step ιn=tn−tn−1 for 1≤n≤N.
Set ¯Λh={(xi,yr,zk)|0≤i≤ℓ1,0≤r≤ℓ2,0≤k≤ℓ3}, and Λh=¯Λh∩Λ, ∂Λh=Λh∩∂Λ. Let unirk be the approximation solution of (1.1)-(1.3).
For n=1,2,⋯,N, by the L1 scheme [35,36], denote
where
Thus, the problem (1.1)-(1.3) can be approximated by
where
Set Sn=d−1n,1=ιαnΓ(2−α), and to establish the ADI scheme, add a small term
onto the left side of Equation (2.2), which is order O(N−2α). Then, one has
Multiplying both sides of (2.6) by Sn, we get
Then, rewriting Eqs. (2.6)-(2.8) in an ADI form:
in which
Set un−13irk=(I−a3Sn1+bSnδ2z)unirk and un−23irk=(I−a2Sn1+bSnδ2y)unirk. To get the solution {unirk} of problem (2.6)-(2.8), we just need to solve three of the one-dimensional equations.
First, for fixed r∈{1,2,⋯,ℓ2−1} and k∈{1,2,⋯,ℓ3−1}, we can solve {un−23irk} in the x-direction as follows:
Then, for fixed i∈{1,2,⋯,ℓ1−1} and k∈{1,2,⋯,ℓ3−1}, we can solve {un−13irk} in the y-direction as follows:
At last, for i∈{1,2,⋯,ℓ1−1} and r∈{1,2,⋯,ℓ2−1}, we utilize {un−23irk,un−13irk} in the z-direction to solve
where the solution {unijk} is attained.
3.
Analysis of stability and convergence
For grid function u={unirk|0≤i≤ℓ1,0≤r≤ℓ2,0≤k≤ℓ3,0≤n≤N}, we define
where δxui−12,r,k=1hx(uirk−ui−1,r,k). The norms ‖δyun‖, ‖δzun‖, ‖δyδzun‖, ‖δxδzun‖, ‖δxδzδ2yun‖ and ‖δyδzδ2xun‖ can be defined similarly.
Set ¯U={uirk|uirk=0 if (xi,yr,zk)∈∂Λh and (xi,yr,zk)∈¯Λh}, for ∀u, w∈¯U, we define
3.1. Stability analysis
This section mainly discusses the stability analysis of ADI difference scheme (2.6)-(2.8).
Lemma 3.1. For grid functions u,ω∈¯U, one has
where the inequality is equal when ω=u.
Proof. Using the Cauchy-Schwartz inequality, one has
The proof is completed. □
For n=1,2,⋯,N and j=1,2,⋯,n−1, define a positive sequence
Lemma 3.2. [37,38] Assume that the sequences {mn}∞n=1, {gn}∞n=1 are nonnegative, and assume that the grid function {un:n=0,1,⋯,N} satisfies u0≥0 and
Then
Lemma 3.3. [39] When , one has
Next, the stability theorem is presented.
Theorem 3.4. For , the solution of ADI scheme (2.6)-(2.8) satisfies
Proof. Multiplying both sides of , and summing over for , we have
Then, using the Cauchy-Schwartz inequality, Young's inequality, and Lemma 3.1, one has
and
Thus, we get
which obtains
Finally, Lemma 3.2 and Lemma 3.3 are applied in (3.3), and the proof is completed. □
3.2. -norm convergence analysis
Lemma 3.5. [24] Assume that , . Then, for , we have
Lemma 3.6. [40] For grid function , one has
Define
and
According to the definition of the small term (2.5), one has
Using a Taylor expansion, one has
Applying Lemma 3.5 and Eq. (3.8), one has
We now show the -norm convergence analysis.
Theorem 3.7. Suppose that , , . one has
where is an -robust constant that does not blow up as .
Proof. Subtracting Eqs. (1.1)-(1.3) from Eqs. (2.6)-(2.8), the error equations are presented via
Multiplying both sides of (3.11) by and summing over for , one has
Since
and
Then, we get
That is,
Thus
By the Poincaré inequality, Eq. (3.9), and Lemma 3.5, we have
From (3.6), (3.14), and (3.16), using Lemma 3.2 and Lemma 3.3, and noting that , one has
Finally, using Lemma 3.6 and the definition of , the proof is completed. □
4.
Numerical experiment
Example 4.1. For problem (1.1)-(1.3), considering , , , , , and
The exact solution is . Set , , , and . In Tables 1–3, we fix to test the error and convergence in time. In Table 1, choosing different , we give errors, convergence orders, and CPU time in time for . In Table 2, set , , , and , and choosing different we give errors, convergence orders and CPU time in time for . In Table 2, the convergence order goes down for and , which can be caused by the round-off error computing the coefficients of the L1 scheme (2.1), see Remark 2.1 in [41]. In Table 3, set , , , and , and choosing different we give errors, convergence orders, and CPU time in time for . Due to the influence of multiple variable coefficient values such as , , , , the degree of freedom of the equation fitting is higher, and there may be some differences for CPU times. From Tables 1–3, it is obtained that the order of convergence in time is , which is consistent with the theoretical analysis.
In Tables 4–6, we fix to present errors and convergence in space for . In Table 4, set , , , and , and by choosing different we give -norm errors, convergence orders and CPU time in space. In Table 5, choosing different we give -norm errors, convergence orders and CPU time in space for , , , and . In Table 6, choosing different we give -norm errors, convergence orders, and CPU time in space for the case where , , , and . In Tables 4-6, we can find that the spatial convergence order is order 2, which is consistent with our theory. From Tables 1-6 we can notice that our scheme is computationally efficient and consistent with the analysis.
5.
Conclusions
In this paper, an ADI scheme is proposed to solve the 3D variable coefficient subdiffusion problem, and our theoretical analysis shows that the method is unconditionally stable, its temporal convergence order is order , and its spatial convergence order is order 2. We present numerical results when taking different coefficients, which show that our ADI scheme is consistent with the theoretical analysis and is very efficient in solving such problems. The ADI technique proposed can reduce the computational cost from spatial discretization. In the future, we will consider some fast and parallel numerical methods [42,43,44] to improve computational efficiency in time.
Acknowledgments
The work was supported by National Natural Science Foundation of China Mathematics Tianyuan Foundation (12226337, 12226340, 12126321, 12126307), Scientific Research Fund of Hunan Provincial Education Department (21B0550, 22C0323, 23C0193), Hunan Provincial Natural Science Foundation of China (2022JJ50083, 2023JJ50164).
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Conflict of interest
The authors declare there is no conflict of interest.