This paper is concerned with the numerical approximations for the variable coefficient fourth-order fractional sub-diffusion equations subject to the second Dirichlet boundary conditions. We construct two effective difference schemes with second order accuracy in time by applying the second order approximation to the time Caputo derivative and the sum-of-exponentials approximation. By combining the discrete energy method and the mathematical induction method, the proposed methods proved to be unconditional stable and convergent. In order to overcome the possible singularity of the solution near the initial stage, a difference scheme based on non-uniform mesh is also given. Some numerical experiments are carried out to support our theoretical results. The results indicate that the our two main schemes has the almost same accuracy and the fast scheme can reduce the storage and computational cost significantly.
Citation: Zhe Pu, Maohua Ran, Hong Luo. Effective difference methods for solving the variable coefficient fourth-order fractional sub-diffusion equations[J]. Networks and Heterogeneous Media, 2023, 18(1): 291-309. doi: 10.3934/nhm.2023011
[1] | Khaled Kefi, Nasser S. Albalawi . Three weak solutions for degenerate weighted quasilinear elliptic equations with indefinite weights and variable exponents. AIMS Mathematics, 2025, 10(2): 4492-4503. doi: 10.3934/math.2025207 |
[2] | Haohao Jia, Feiyao Ma, Weifeng Wo . Large positive solutions to an elliptic system of competitive type with nonhomogeneous terms. AIMS Mathematics, 2021, 6(8): 8191-8204. doi: 10.3934/math.2021474 |
[3] | Mengyu Wang, Xinmin Qu, Huiqin Lu . Ground state sign-changing solutions for fractional Laplacian equations with critical nonlinearity. AIMS Mathematics, 2021, 6(5): 5028-5039. doi: 10.3934/math.2021297 |
[4] | Labudan Suonan, Yonglin Xu . Existence of solutions to mixed local and nonlocal anisotropic quasilinear singular elliptic equations. AIMS Mathematics, 2023, 8(10): 24862-24887. doi: 10.3934/math.20231268 |
[5] | Diane Denny . Existence of a solution to a semilinear elliptic equation. AIMS Mathematics, 2016, 1(3): 208-211. doi: 10.3934/Math.2016.3.208 |
[6] | Dan Wang, Yongxiang Li . Existence and uniqueness of radial solution for the elliptic equation system in an annulus. AIMS Mathematics, 2023, 8(9): 21929-21942. doi: 10.3934/math.20231118 |
[7] | Zhiying Deng, Yisheng Huang . Existence and multiplicity results for a singular fourth-order elliptic system involving critical homogeneous nonlinearities. AIMS Mathematics, 2023, 8(4): 9054-9073. doi: 10.3934/math.2023453 |
[8] | Khaled Kefi . Existence and multiplicity of triple weak solutions for a nonlinear elliptic problem with fourth-order operator and Hardy potential. AIMS Mathematics, 2024, 9(7): 17758-17773. doi: 10.3934/math.2024863 |
[9] | Nafissa Toureche Trouba, Mohamed E. M. Alngar, Reham M. A. Shohib, Haitham A. Mahmoud, Yakup Yildirim, Huiying Xu, Xinzhong Zhu . Novel solitary wave solutions of the (3+1)–dimensional nonlinear Schrödinger equation with generalized Kudryashov self–phase modulation. AIMS Mathematics, 2025, 10(2): 4374-4411. doi: 10.3934/math.2025202 |
[10] | Ninghe Yang . Exact wave patterns and chaotic dynamical behaviors of the extended (3+1)-dimensional NLSE. AIMS Mathematics, 2024, 9(11): 31274-31294. doi: 10.3934/math.20241508 |
This paper is concerned with the numerical approximations for the variable coefficient fourth-order fractional sub-diffusion equations subject to the second Dirichlet boundary conditions. We construct two effective difference schemes with second order accuracy in time by applying the second order approximation to the time Caputo derivative and the sum-of-exponentials approximation. By combining the discrete energy method and the mathematical induction method, the proposed methods proved to be unconditional stable and convergent. In order to overcome the possible singularity of the solution near the initial stage, a difference scheme based on non-uniform mesh is also given. Some numerical experiments are carried out to support our theoretical results. The results indicate that the our two main schemes has the almost same accuracy and the fast scheme can reduce the storage and computational cost significantly.
Let Ω be a bounded domain of RN(N≥3) with smooth boundary ∂Ω. In this paper, we consider the existence of W1,10(Ω)solutions to the following elliptic problem
{−div(M(x,u)∇u)+|∇u|2uθ=f,x∈Ω,u=0,x∈∂Ω, | (1.1) |
where NN−1≤θ<2, M:Ω×R→RN2 is a symmetric Carathéodory matrix function, which satisfies the following assumptions: for some real constants γ>0, α>0,β>0,
|M(x,s)|≤β,M(x,s)ξ⋅ξ≥α(a(x)+|s|)γ|ξ|2, | (1.2) |
for almost every x∈Ω, (s,ξ)∈R×RN, where a(x) is a measurable function, such that
0<ζ≤a(x)≤ρ, | (1.3) |
for some positive constants ζ,ρ.
We note that there are two difficulties in dealing with (1.1), the first one is the fact that, due to hypothesis (1.2), the differential operator A(u)=−div(M(x,u)∇u) is well defined in H10(Ω), but it not coercive on H10(Ω) when u is large enough. Therefore, the classical Leray-Lions theorem cannot be applied even if f is sufficiently regular. The second difficulty is dealing with lower order term which singular natural growth with respect to the gradient. In order to overcome these difficulties, we approximate problem (1.1) by means of truncations in M(x,s) to get a coercive differential operator on H10(Ω).
The existence of W1,10(Ω) solution to elliptic problem has been studied by many authors. Boccardo and Croce [4] proved the existence of W1,10(Ω) solutions to problem
{−div(a(x)∇u(1+|u|)γ)=f,x∈Ω,u=0,x∈∂Ω, |
where a:Ω→R is a measurable function which satisfies (1.3), f∈Lm(Ω) with
m=NN+1−γ(N−1),1N−1<γ<1. |
In the literature [6], the authors considered the existence and regularity of solutions to the following elliptic equation with noncoercivity
{−div(a(x,u)∇u)=f,x∈Ω,u=0,x∈∂Ω, | (1.4) |
where Ω is an open bounded subset of RN(N≥3), f∈Lm(Ω) and a(x,s):Ω×R→R is a Carathéodory function which satisfies
α(1+|s|)γ≤a(x,s)≤β, |
where 0≤γ<1.The existence results of solutions to problem (1.4)are as following:
● There exists a weak solution u∈H10(Ω)∩L∞(Ω) to (1.4) if m>N2.
● There exists a weak solution u∈H10(Ω)∩Lr(Ω) to (1.4) with r=Nm(1−γ)N−2m if
2NN+2−γ(N−2)≤m<N2. |
● There exists a distributional solution u∈W1,q0(Ω) to (1.4) with q=Nm(1−γ)N−m(1+γ)<2 if
NN+1−γ(N−1)<m<2NN+2−γ(N−2). |
In [15], the under the assumption (1.2)-(1.3), Souilah proved the existence results of solutions to problem
{−div(M(x,u)∇u)+|∇u|2uθ=f+λur,x∈Ω,u=0,x∈∂Ω, | (1.5) |
where 0<θ<1,0<r<2−θ,λ>0,f∈Lm(Ω)(m≥1). There exists at least a solution to problem (1.5):
● If 2N2N−θ(N−2)≤m<N2, then u∈H10(Ω)∩L∞(Ω).
● If 1<m<2N2N−θ(N−2), then u∈W1,q0(Ω) with q=Nm(2−θ)N−mθ.
● If m≥N2, then u∈H10(Ω)∩L∞(Ω).
Moreover, the existence of solutions u∈H10(Ω) to problem (1.5) with λ=0 have been obtained in [9]. Some other related results see [1,3,5,7,10,11,12,14,16].
Based on the above research results, the aim of this article is to study the existence of W1,10(Ω) solution to problem (1.1).
In order to state the main results of this paper, the following definition need to be introduced. We use the following notion of distributional solution to problem (1.1).
Definition 1.1. We say that u∈W1,10(Ω) is a distributional solution to problem (1.1) if u>0 in Ω, |∇u|2uθ∈L1(Ω) and
∫ΩM(x,u)∇u⋅∇φ+∫Ω|∇u|2uθφ=∫Ωfφ, |
for every φ∈C∞0(Ω).
Our main results are following:
Theorem 1.2. Assume that (1.2)-(1.3) hold, f∈Lm(Ω) is a nonnegative function with
m=N2N−θ(N−1),NN−1<θ<2. | (1.6) |
Then there exists a distributional solution u∈W1,10(Ω) to problem (1.1).
Remark 1.3. Notice that the result of previous theorem do not depend on γ.
Remark 1.4. Observe that, m>1 if and only if θ>NN−1.
For f∈L1(Ω), we have the following theorem.
Theorem 1.5. Assume (1.2)-(1.3) hold, f∈L1(Ω) is a nonnegative function and θ=NN−1. Then there exists a distributional solution u∈W1,10(Ω) to problem (1.1).
The paper is organized as follows. In section 2, we collect some definitions and useful tools. The proof of Theorem 1.2 and 1.5 be given in section 3.
In order to prove our main results, we need to introduce a basic definition and some lemmas.
Definition 2.1. For all k≥0, the truncation function defined by
Tk(s)=max{−k,min{k,s}},Gk(s)=s−Tk(s). |
Let 0<ε<1, we approximate problem (1.1) by the following non-singular problem
{−div(M(x,T1ε(uε))∇uε)+uε|∇uε|2(|uε|+ε)θ+1=fε,x∈Ω,uε=0,x∈∂Ω, | (2.1) |
where fε=T1ε(f). Problem (2.1) admits at least a solution uε∈H10(Ω)∩L∞(Ω) by Theorem 2 of [8]. Due to the fact that fε≥0 and quadratic lower order term has the same sign of the solution, it is easy to prove that uε≥0 by taking u−ε as a test function in (2.1).
Lemma 2.2. Let uε be the solutions to problem (2.1). Then
∫Ωuε|∇uε|2(uε+ε)θ+1≤∫Ωf. | (2.2) |
Proof. For fixed h>0, taking Th(uε)h as a test function in (2.1). Dropping the first term, we obtain
∫Ωuε|∇uε|2(uε+ε)θ+1Th(uε)h≤∫ΩfεTh(uε)h. |
Using the fact that fε≤f and Th(uε)h≤1, then
∫Ωuε|∇uε|2(uε+ε)θ+1Th(uε)h≤∫Ωf. |
Letting h→0, we deduce (2.2) by the Fatou Lemma.
Lemma 2.3. Let δ>0 and 0<ε<1. Then there exists C>0, such that
αδ(t+ε)θ−2(ρ+t)γ+tt+ε≥C. |
for every t≥0.
Proof. Clearly, if t≥ε, we have tt+ε≥12, while if t<ε, we have
αδ(t+ε)θ−2(ρ+t)γ≥αδ(ρ+t)γ(2ε)2−θ≥αδ22−θ(ρ+1)γ, |
since ε<1. Therefore, Lemma 2.3 is proved.
In this section, C denotes a generic constant whose value might change from line to line. We prove the existence results of Theorems 1.2 and 1.5 by considering the following approximate problem
{−div(M(x,T1ε(uε))∇uε)+uε|∇uε|2(uε+ε)θ+1=fε,x∈Ω,uε=0,x∈∂Ω. | (3.1) |
Proof of Theorem 1.2. Step 1: Let δ=θ−NN−1, then δ>0 by (1.6). Choosing (uε+ε)δ−(uε+ε)δ−1 as a test function in the approximate problem (3.1), we find
∫ΩM(x,T1ε(uε))∇uε⋅∇uε[δ(uε+ε)δ−1+(1−δ)(uε+ε)δ−2]+∫Ωuε(uε+ε)δ|∇uε|2(uε+ε)θ+1=∫Ωuε|∇uε|2(uε+ε)θ+1(uε+ε)δ−1+∫Ωfε[(uε+ε)δ−(uε+ε)δ−1]. |
Combining (1.2)-(1.3) and dropping the positive term, we obtain
∫Ω|∇uε|2(uε+ε)δ−θ[α(1−δ)(uε+ε)θ−2(ρ+uε)γ+uεuε+ε]≤∫Ωuε|∇uε|2(uε+ε)θ+1(uε+ε)δ−1+∫Ωfε(uε+ε)δ. |
Since 1−δ>0, according to Lemma 2.3, we have
C∫Ω|∇uε|2(uε+ε)δ−θ≤∫Ωuε|∇uε|2(uε+ε)θ+1(uε+ε)δ−1+∫Ωfε(uε+ε)δ. |
Using the fact that uε≥0,fε≤f and (2.2), we obtain
C∫Ω|∇uε|2(uε+ε)δ−θ≤εδ−1∫Ωuε|∇uε|2(uε+ε)θ+1+∫Ωf(uε+ε)δ≤εδ−1∫Ωf+∫Ωf(uε+ε)δ. | (3.2) |
Observe that the left hand side of (3.2) can be rewritten as
C∫Ω|∇[(uε+ε)δ−θ+22−εδ−θ+22]|2. | (3.3) |
Then, (3.2) and (3.3) imply
C∫Ω|∇[(uε+ε)δ−θ+22−εδ−θ+22]|2≤εδ−1∫Ωf+∫Ωf(uε+ε)δ. | (3.4) |
By the Sobolev inequality, satisfy
[∫Ω|(uε+ε)δ−θ+22−εδ−θ+22|2∗]22∗≤C∫Ω|∇[(uε+ε)δ−θ+22−εδ−θ+22]|2. | (3.5) |
Using the Hölder inequality and (3.4)-(3.5), we get
[∫Ω|(uε+ε)δ−θ+22−εδ−θ+22|2∗]22∗≤C‖f‖Lm(Ω)+C‖f‖Lm(Ω)[∫Ω(uε+ε)δm′]1m′. |
Since |(t+ε)s−εs|2∗≥C[(t+ε)2∗s−1] for every t≥0 and for suitable constant C independent on ε, then we find
(∫Ω[(uε+ε)2∗(δ−θ+2)2−1])22∗≤C‖f‖Lm(Ω)+C‖f‖Lm(Ω)[∫Ω(uε+ε)δm′]1m′. | (3.6) |
Thanks to the choice of δ, we have
2∗(δ−θ+2)2=δm′=NN−1. |
Moreover 22∗>1m′ since m<N2. Then (3.6) implies that
∫ΩuNN−1ε≤C. | (3.7) |
Observe that δ−θ=−NN−1, then, (3.2), (3.7) follow
∫Ω|∇uε|2(ε+uε)NN−1≤C. | (3.8) |
Combining (3.7)-(3.8) with the Hölder inequality, we obtain
∫Ω|∇uε|=∫Ω∇uε(ε+uε)N2N−2(ε+uε)N2N−2≤[∫Ω|∇uε|2(ε+uε)NN−1]12[∫Ω(ε+uε)NN−1]12≤C. |
Then we get that {uε} is bounded in W1,10(Ω). Hence, there exists a subsequence {uε}, which converges to a measurable function u a.e. in Lr(Ω) with 1≤r<NN−1.
Step 2: First, we are going to estimate ∫{uε≥k}|∇uε|. Choosing [(uε+ε)δ−(k+ε)δ]+ as a test function in (3.1). By (1.2)-(1.3) and Lemma 2.3, we have
∫{uε≥k}|∇uε|2(ε+uε)NN−1≤(∫{uε≥k}|f|m)1m(∫{uε≥k}(ε+uε)NN−1)1m′≤C(∫{uε≥k}|f|m)1m. |
Using the Hölder inequality and (3.7), we find
∫{uε≥k}|∇uε|=∫{uε≥k}∇uε(ε+uε)N2N−2(ε+uε)N2N−2≤C(∫{uε≥k}|f|m)12m. | (3.9) |
Choosing Tk(uε) as a test function in (3.1). Dropping the nonnegative lower order term, by (1.2)-(1.3) and the boundedness of uε in LNN−1(Ω), we get
∫Ω|∇Tk(uε)|2≤k(ρ+k)γα‖f‖L1(Ω). | (3.10) |
This implies that Tk(uε)⇀Tk(u) weakly in W1,20(Ω).
Let E be a measurable subset of Ω, and i=1,⋯,N. By the Hölder inequality and (3.9)-(3.10), we obtain
∫E|∂uε∂xi|≤∫E|∇uε|≤∫E|∇Tk(uε)|+∫{uε≥k}|∇uε|≤meas(E)12(∫E|∇Tk(uε)|2)12+C(∫{uε≥k}|f|m)12m. | (3.11) |
The estimates (3.7) and (3.11) shows that the sequence {∂uε∂xi} is equi-integrable. Thus, by the Dunford–Pettis theorem, there exists a subsequence {uε} and Vi in L1(Ω), such that ∂un∂xi⇀Vi in L1(Ω). Since ∂uε∂xi is the distributional partial derivative of uε, then we have
∫Ω∂uε∂xiφ=−∫Ωuε∂φ∂xi,∀φ∈C∞0(Ω), |
for every ε>0.
Since ∂uε∂xi⇀Vi in L1(Ω) and uε→u in L1(Ω), we find
∫ΩViφ=−∫Ωu∂φ∂xi,∀φ∈C∞0(Ω). |
This implies that Vi=∂u∂xi for every i.
Step 3: We prove that uε|∇uε|2(uε+ε)θ+1 is equi-integrable. Let E⊂⊂Ω, then
∫Euε|∇uε|2(uε+ε)θ+1≤∫E∩{uε≤k}uε|∇uε|2(uε+ε)θ+1+∫E∩{uε≥k}uε|∇uε|2(uε+ε)θ+1. |
For every subset E⊂⊂Ω,
∫E∩{uε≤k}uε|∇uε|2(uε+ε)θ+1≤∫E∩{uε≤k}1uθε|∇Tk(uε)|2≤C∫E∩{uε≤k}|∇Tk(uε)|2, |
since uε≥C>0 in E by Proposition 2 of [9]. Moreover, since Tk(uε)⇀Tk(u) weakly in W1,20(Ω), then there exists εn,δ>0, such that
∫E∩{uε≤k}|∇Tk(uε)|dx≤ϵ2,∀ε≥εn, | (3.12) |
for every ϵ>0 if μ(E)<δ.
Choosing T1(uε−Tk−1(uε)) as a test function in the approximate problem (3.1), dropping the nonnegative term, we have
∫{uε≥k}uε|∇uε|2(uε+ε)θ+1≤∫{uε≥k−1}f. | (3.13) |
Observe there exists a constant C>0, such that μ(uε≥k−1)≤Ck−1. As uε are uniformly bounded in LNN−1(Ω). This implies the right hand side of (3.13) converges to 0 as k→∞. Thus, we deduce there exists k0>1, such that
∫{uε≥k}uε|∇uε|2(uε+ε)θ+1≤ϵ2,∀k>k0, | (3.14) |
for every ϵ>0. The (3.12), (3.14) imply that uε|∇uε|2(uε+ε)θ+1 is equi-integrable and converges a.e. to |∇u|2uθ.
Let u the weak limit of the sequence of approximated solutions uε. Thanks to (2.2), we have
∫Ωuε|∇uε|2(uε+ε)θ+1≤∫Ωf. |
Using the Fatou lemma, that uε convergence to u a.e, ∇uε convergence to ∇u a.e and the strict positivity of uε imply
∫Ω|∇u|2uθ≤∫Ωf≤C. |
This show that |∇u|2uθ∈L1(Ω).
Since uε is bounded and ∇uε→∇u a.e, it follow M(x,T1ε(uε)∇uε→M(x,u)∇u a.e. Hence, we can pass to the limit in (3.1). Thus prove that u∈W1,10(Ω) is a distributional solution of (1.1) and yields the conclusion of the proof of Theorem 1.1.
Proof of Theorem 1.5. Step 1: For 0<ε<1, according to Lemma 2.2, we have
12θ+1∫{uε≥1}|∇uε|2uθε≤∫{uε≥1}uε|∇uε|2(uε+ε)θ+1≤‖f‖L1(Ω). | (3.15) |
By the Sobolev inequality, (3.15) lead to
[∫Ω|u2−θ2ε−1|2∗]22∗≤C‖f‖L1(Ω), | (3.16) |
which implies that
[∫Ωu(2−θ)2∗2ε]22∗≤C+C‖f‖L1(Ω). | (3.17) |
Observe that θ=(2−θ)2∗2=NN−1. Then (3.17) shows that
∫ΩuNN−1ε≤C. | (3.18) |
Using the Hölder inequality and (3.15), (3.18), we obtain
∫Ω|∇G1(uε)|=∫{uε≥1}|∇G1(uε)|uθ2εuθ2ε≤[∫{uε≥1}|∇uε|2uθε]12[∫{uε≥1}uθε]12≤C‖f‖L1(Ω). |
This fact show that G1(uε) is bounded in W1,10(Ω).
Choosing T1(uε) as a test function in (3.1), it is easy to prove that T1(uε) is bounded in H10(Ω)), hence in W1,10(Ω). Since uε=G1(uε)+T1(uε), we deduce that uε is bounded in W1,10(Ω).
Moreover, due to (3.15) and the Hölder inequality, we have
∫{uε≥k}|∇uε|=∫{uε≥k}|∇uε|uθ2εuθ2ε≤C‖f‖12L1(Ω). | (3.19) |
That (3.10), (3.19) implies, for every measurable subset E, we have
∫E|∂uε∂xi|≤∫E|∇uε|≤∫E|∇Tk(uε)|+∫{uε≥k}|∇uε|≤meas(E)12[k(ρ+k)γα‖f‖L1(Ω)]12+C‖f‖12L1(Ω). |
Thus, we prove that uε⇀u in W1,10(Ω). Then pass to the limit in problem (3.1), as in the proof of Theorem 1.2, it is sufficient to observe that u∈W1,10(Ω) is a distributional solution of (1.1). This concludes the proof the Theorem 1.5.
In this paper, we main consider the existence of W1,10(Ω) solutions to a elliptic equation with principal part having noncoercivity. The main results show that the singular quadratic term has an important impact on this existence.
This research was partially supported by the National Natural Science Foundation of China (No. 11761059), Program for Yong Talent of State Ethnic Affairs Commission of China (No. XBMU-2019-AB-34), Fundamental Research Funds for the Central Universities (No.31920200036) and First-rate Discipline of Northwest Minzu University.
The authors declare that there is no conflict of interests regarding the publication of this article.
[1] |
V. Srivastava, K. N. Rai, A multi-term fractional diffusion equation for oxygen delivery through a capillary to tissues, Math. Comput. Model., 51 (2010), 616–624. https://doi.org/10.1016/j.mcm.2009.11.002 doi: 10.1016/j.mcm.2009.11.002
![]() |
[2] |
Y. Z. Povstenko, Fractional Cattaneo-type equations and generalized thermoelasticity, J. Therm. Stresses, 34 (2011), 97–114. https://doi.org/10.1080/01495739.2010.511931 doi: 10.1080/01495739.2010.511931
![]() |
[3] |
F. Mainardi, M. Reberto, R. Gorenflo, E. Scalas, Fractional caculus and continuous-time finance II: the waiting-time distribution, Physica A, 287 (2000), 468–481. https://doi.org/10.1016/S0378-4371(00)00386-1 doi: 10.1016/S0378-4371(00)00386-1
![]() |
[4] | D. Benson, R. Schumer, M. Meerschaert, S. Wheatcraft, Fractional dispersion, Lévy motion, and the MADE tracer tests, Transport Porous Med., 42 (2001), 211–240. https://doi.org/10.1023/A:1006733002131 |
[5] | I. Podlubny, Fractional Differential Equations, New York: Academic press, 1999. |
[6] |
J. B. Greer, A. L. Bertozzi, G. Sapiro, Fourth order partial differential equations on general geometries, J. Comput. Phys., 216 (2006), 216–246. https://doi.org/10.1016/j.jcp.2005.11.031 doi: 10.1016/j.jcp.2005.11.031
![]() |
[7] | F. Memoli, G. Sapiro, P. Thompson, Implicit brain imaging, NeuroImage, 23 (2004), 179–188. https://doi.org/10.1016/j.neuroimage.2004.07.072 |
[8] |
M. Hofer, H. Pottmann, Energy-minimizing splines in manifolds, ACM Trans. Graph, 23 (2004), 284–293. https://doi.org/10.1145/1015706.1015716 doi: 10.1145/1015706.1015716
![]() |
[9] |
D. Halpern, O. E. Jensen, J. B. Grotberg, A theoretical study of surfactant and liquid delivery into the lung, J. Appl. Physiol., 85 (1998), 333–352. https://doi.org/10.1152/jappl.1998.85.1.333 doi: 10.1152/jappl.1998.85.1.333
![]() |
[10] | T. G. Myers, J. P. F. Charpin, S. J. Chapman, The flow and solidification of a thin fluid film on an arbitrary three-dimensional surface, Phys. Fluids, 14 (8) (2002), 2788–2803. https://doi.org/10.1063/1.1488599 |
[11] |
T. G. Myers, J. P. F. Charpin, A mathematical model for atmospheric ice accretion and water flow on a cold surface, Int. J. Heat Mass Transf., 47 (2004), 5483–5500. https://doi.org/10.1016/j.ijheatmasstransfer.2004.06.037 doi: 10.1016/j.ijheatmasstransfer.2004.06.037
![]() |
[12] |
X. Hu, L. Zhang, On finite difference methods for fourth-order fractional diffusion-wave and subdiffusion systems, Appl. Math. Comput., 218 (2012), 5019–5034. https://doi.org/10.1016/j.amc.2011.10.069 doi: 10.1016/j.amc.2011.10.069
![]() |
[13] |
X. Hu, L. Zhang, A compact finite difference scheme for the fourth-order fractional diffusion-wave system, Comput. Phys. Commun., 182 (2011), 1645–1650. https://doi.org/10.1016/j.cpc.2011.04.013 doi: 10.1016/j.cpc.2011.04.013
![]() |
[14] |
C. Ji, Z. Sun, Z. Hao, Numerical algorithms with high spatial accuracy for the fourth-order fractional sub-diffusion equations with the first Dirichlet boundary conditions, J. Sci. Comput., 66 (2015), 1148–1174. https://doi.org/10.1007/s10915-015-0059-7 doi: 10.1007/s10915-015-0059-7
![]() |
[15] |
P. Zhang, H. Pu, A second-order compact difference scheme for the fourth-order fractional sub-diffusion equation, Numer Algor, 76 (2017), 573–598. https://doi.org/10.1007/s11075-017-0271-7 doi: 10.1007/s11075-017-0271-7
![]() |
[16] |
A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, J. Comput. Phys., 280 (2015), 424–438. https://doi.org/10.1016/j.jcp.2014.09.031 doi: 10.1016/j.jcp.2014.09.031
![]() |
[17] |
M. Ran, C. Zhang, New compact difference scheme for solving the fourth-order time fractional sub-diffusion equation of the distributed order, Appl. Numer. Math, 129 (2018), 58–70. https://doi.org/10.1016/j.apnum.2018.03.005 doi: 10.1016/j.apnum.2018.03.005
![]() |
[18] |
X. Zhao, Q. Xu, Efficient numerical schemes for fractional sub-diffusion equation with the spatially variable coefficient, Appl. Math. Model., 38 (2014), 3848–3859. https://doi.org/10.1016/j.apm.2013.10.037 doi: 10.1016/j.apm.2013.10.037
![]() |
[19] |
S. Vong, P. Lyu, Z. Wang, A compact difference scheme for fractional sub-diffusion equations with the spatially variable coefficient under Neumann boundary conditions, J. Sci. Comput., 66 (2016), 725–739. https://doi.org/10.1007/s10915-015-0040-5 doi: 10.1007/s10915-015-0040-5
![]() |
[20] | A. A. Samarskii, V. V. Andreev, Difference methods for elliptic equations (in Russian), Moscow: Izdatel'stvo Nauka, 1976. |
[21] | Z. Sun, Numerical Methods of Partial Differential Equations (in Chinese), Beijing: Science Press, 2012. |
[22] |
T. Yan, Z. Sun, J. Zhang, Fast evalution of the Caputo fractional derivative and its applications to fractional diffusion equations: A second-order scheme, Commun. Comput. Phys., 22 (2017), 1028–1048. https://doi.org/10.4208/cicp.OA-2016-0136 doi: 10.4208/cicp.OA-2016-0136
![]() |
[23] |
H. Liao, D. Li, J. Zhang, Sharp error estimate of the nonuniform L1 formula for linear reaction-subdiffusion equations, SIAM J. Numer. Anal., 56 (2018), 1112–1133. https://doi.org/10.1137/17M1131829 doi: 10.1137/17M1131829
![]() |
[24] |
B. Jin, R. Lazarov, Z. Zhou, Numerical methods for time-fractional evolution equations with nonsmooth data: A concise overview, Comput. Methods Appl. Mech. Engrg., 346 (2019), 332–358. https://doi.org/10.1016/j.cma.2018.12.011 doi: 10.1016/j.cma.2018.12.011
![]() |
[25] |
M. Stynes, E. O'Riordan, J. L. Gracia, Error analysis of a finite difference method on graded meshes for a time-fractional diffusion equation, SIAM J. Numer. Anal., 55 (2017), 1057–1079. https://doi.org/10.1137/16M1082329 doi: 10.1137/16M1082329
![]() |
[26] |
M. Cui, A compact difference scheme for time-fractional Dirichlet biharmonic equation on temporal graded meshe, East Asian J. Appl. Math., 11 (2021), 164–180. https://doi.org/10.4208/eajam.270520.210920 doi: 10.4208/eajam.270520.210920
![]() |
1. | Mounim El Ouardy, Youssef El Hadfi, Some Nonlinear Parabolic Problems with Singular Natural Growth Term, 2022, 77, 1422-6383, 10.1007/s00025-022-01631-6 | |
2. | Abdelkrim Moussaoui, Kamel Saoudi, Existence and Location of Nodal Solutions for Quasilinear Convection–Absorption Neumann Problems, 2024, 47, 0126-6705, 10.1007/s40840-024-01669-5 | |
3. | Hocine Ayadi, Rezak Souilah, The impact of a singular first-order term in some degenerate elliptic equations involving Hardy potential, 2024, 9, 2662-2009, 10.1007/s43036-024-00324-x |