[BACK]
images Computer Modeling in
Engineering & Sciences
images

DOI: 10.32604/cmes.2021.014950

ARTICLE

Quadratic Finite Volume Element Schemes over Triangular Meshes for a Nonlinear Time-Fractional Rayleigh-Stokes Problem

Yanlong Zhang1, Yanhui Zhou2 and Jiming Wu3,*

1Graduate School of China Academy of Engineering Physics, Beijing, 100088, China
2School of Data and Computer Science, Sun Yat-Sen University, Guangzhou, 510275, China
3Institute of Applied Physics and Computational Mathematics, Beijing, 100088, China
*Corresponding Author: Jiming Wu. Email: wu_jiming@iapcm.ac.cn
Received: 11 November 2020; Accepted: 15 January 2021

Abstract: In this article, we study a 2D nonlinear time-fractional Rayleigh-Stokes problem, which has an anomalous sub-diffusion term, on triangular meshes by quadratic finite volume element schemes. Time-fractional derivative, defined by Caputo fractional derivative, is discretized through L2-1σ formula, and a two step scheme is used to approximate the time first-order derivative at time tn-α/2, where the nonlinear term is approximated by using a matching linearized difference scheme. A family of quadratic finite volume element schemes with two parameters are proposed for the spatial discretization, where the range of values for two parameters are β1(0,1/2), β2(0,2/3). For testing the precision of numerical algorithms, we calculate some numerical examples which have known exact solution or unknown exact solution by several kinds of quadratic finite volume element schemes, and contrast with the results of an existing quadratic finite element scheme by drawing diversified comparison plots and showing the detailed data of L2 error results and convergence orders. Numerical results indicate that, L2 error estimate of one scheme with parameters β1=(3-3)/6, β2=(6+3-21+63)/9 is 𝓞(h3+Δt2), and L2 error estimates of other schemes are 𝓞(h2+Δt2), where h and Δt denote the spatial and temporal discretization parameters, respectively.

Keywords: Quadratic finite volume element schemes; anomalous sub-diffusion term; L2 error estimate; quadratic finite element scheme

1  Introduction

Recently, due to the widespread use of fractional partial differential equations (FPDEs), such as dispersion in a porous medium, statistical mechanics, mathematical biology and so on, numerical solution of FPDEs becomes one of the frontier fields in the research. Fractional partial differential equations can be roughly classified into three categories: Space FPDEs [110], time FPDEs [1129] and space-time FPDEs [3034]. Anomalous sub-diffusion equations, one type of time FPDEs, arise in some physical and biological processes. And the study of FPDEs with anomalous sub-diffusion terms, such as modified anomalous sub-diffusion equations [17,18], fractional Cable equations [11,16] or others, is also meaningful and popular. The problem considered in this article, which belongs to a nonlinear time-fractional Rayleigh-Stokes problem [1922] applied in some non-Newtonian fluids, is a variant of the Stokes’ first problems and Rayleigh-Stokes problems [3538], and it is important in physics and engineering.

At present, numerical simulation is an important and effective way to solve partial differential equations, and the relevant numerical methods can be finite difference methods [7,8,1417], finite element methods (FEMs) [1,2,9,1113,21,29,3133,39], meshless methods [40,41], finite volume methods [36,4254] and so on. Of course, the research for FPDEs by finite volume element methods (FVEMs) [10,2328] has no exception for the local conservation and simple implementation. Sayevand et al. [23] presented a spatially semi-discrete piecewise linear FVEM for the time-fractional sub-diffusion problem and obtained some error estimates of the solution in both FEMs and FVEMs. A linear finite volume element scheme for the 2D time-fractional anomalous sub-diffusion equations was studied and analyzed by Karaa et al. [24], where the convergence rate was of order h2+Δt1+α in the L(L2) norm and the results were improved in [25] for both smooth and nonsmooth initial data. Badr et al. [26] proposed a linear FVEM for the time-fractional advection diffusion problem in one-dimension, and proved that the fully discrete scheme is unconditionally stable. Furthermore, Yazdani et al. [10] solved a space-fractional advection-dispersion problem in one-dimension by using linear FVEM and proved it is stable when the mesh grid size is small enough. Zhao et al. [27] constructed a mixed finite volume element scheme for the time-fractional reaction-diffusion equation, and showed the unconditional stability analysis for it. Moreover, Zhao et al. [28] proposed a linear FVEM for the nonlinear time-fractional mobile/immobile transport equations on triangular grids, and obtained the optimal priori error estimates in L(L2) and L2(H1) norms. To our knowledge, the study of high order finite volume element methods for 2D FPDEs is undiscovered.

There are some research about quadratic finite volume element methods for solving partial differential equations on triangular meshes. Tian et al. [42] presented quadratic element generalized differential methods to solve elliptic equations where two parameters of the quadratic element were β1=β2=1/3 (referring to the definition Eqs. (4), (5)). Liebau [43] solved one type of elliptic boundary value problems by a quadratic element scheme with parameters β1=1/4,β2=1/3, and proved 𝓞(h2) error estimate under some assumption conditions. Xu et al. [44] started to study the structure about two parameters of the quadratic element, and improved some existing coercivity results. Chen et al. [45] established a general framework for construction and analysis of the higher-order finite volume methods. Wang et al. [46] established a unified framework to perform the L2 error analysis for high order finite volume methods on triangular meshes, and proposed a new quadratic scheme with parameters β1=(3-3)/6,β2=(6+3-21+63)/9 to achieve the optimal L2 convergence order. An unconditionally stable quadratic finite volume scheme with parameters β1=β2=(3-3)/6 for elliptic equations was presented by Zou [47], and it had optimal convergence orders under H1 norms. For the quadratic finite volume element schemes with parameters β1=(3-3)/6,β2(0,2/3), Zhou et al. [48] obtained an analytic minimum angle condition and an optimal H1 error estimates under the improved coercivity result. Moreover, a unified framework for the coercivity analysis of a class of quadratic schemes with parameters β1(0,1/2),β2(0,2/3) was established for elliptic boundary value problems [49], which covered all the existing quadratic schemes of Lagrange type, and minimum angle conditions of the existing literatures are improved. All the above papers are mainly confined to the elliptic problems and we have seen some applications to other problems [5053], but none relevant study for FPDEs.

In this article, the quadratic finite volume element method is proposed to solve one class of FPDEs, that is, a 2D nonlinear time-fractional Rayleigh-Stokes problem with the time-fractional derivative defined by Caputo fractional derivative. In spatial direction, this problem is solved by a class of quadratic finite volume element schemes with two parameters β1 and β2. Moreover, this problem is discretized at time tn-α/2 and the time-fractional derivative is discretized through L2-1σ formula in time direction. Numerical experiments indicate the efficiency of the schemes, specifically, the L2 error estimate of one scheme is 𝓞(h3+Δt2), and L2 error estimates of other schemes are 𝓞(h2+Δt2). We find that the new finite volume element schemes are comparable with an existing finite element scheme [29].

The outline of this paper is as follows. In Section 2, we describe in details the specific algorithm steps of the quadratic finite volume element schemes over triangular meshes, and finally obtain the fully discrete schemes. In Section 3, some numerical experiments are performed to investigate the performance of the quadratic finite volume element schemes. The numerical results are also compared with those of an existing quadratic finite element scheme. A brief conclusion ends this article in last section.

2  Quadratic Finite Volume Element Schemes

2.1 Preliminary

In this article, we construct a family of quadratic finite volume element schemes to solve the following 2D nonlinear time-fractional Rayleigh-Stokes problem:

u(x,t)tαΔu(x,t)tαΔu(x,t)=f(u)+g(x,t),(x,t)Ω×(0,T], (1)

u(x,t)=0,(x,t)Ω×(0,T], (2)

u(x,0)=u0(x),xΩ, (3)

where x=(x,y)Ω2 and αΔu(x,t)tα=1Γ(1-α)0tΔu(x,τ)τdτ(t-τ)α, denotes the Caputo fractional derivative of order α(0,1), which is an anomalous sub-diffusion term. f(u) is a nonlinear term subjected to the following conditions, |f(u)|C1|u| and |f(u)|C2 where C1 and C2 are positive constants, and g(x,t) is the source term.

For the numerical solution of Eqs. (1)(3), the space domain Ω is first triangulated to get the so-called primary mesh 𝓒h, see the solid line segments in Fig. 1. The trial function space Uh is then defined with respect to 𝓒h, given by Uh={uhC(Ω̄):uhK2,K𝓒h,uhΩ=0}, where K denotes a generic triangular element and 2 is the set of all polynomials of degree less than or equal to 2. It is easy to see that Uh is the same one as that of the standard quadratic finite element method, but the test function space here is different.

images

Figure 1: The primary mesh and dual mesh

In order to define the test function space, we need to construct a dual mesh associated with 𝓒h, see the dashed line segments in Fig. 1. For any triangular element K=ΔP1P2P3 with P1=(x1,y1),P2=(x2,y2) and P3=(x3,y3), we denote the barycenter of K as Q, and the midpoints of the three sides as M1,M2 and M3, respectively, see Fig. 2. Gk,k+1 and Gk+1,k are the two points on PkPk+1, satisfying

|PkGk,k+1||PkPk+1|=|Gk+1,kPk+1||PkPk+1|=β1(0,12),(4)

while Pk,k+1 is a point on PkMk+1 such that

|PkPk,k+1||PkMk+1|=β2(0,23),(5)

where k is a periodic index with Period 3. Using the above notations, K can be further partitioned into six subcells, i.e., three quadrilaterals and three pentagons, see Fig. 2. Let 𝓝h be the set of all vertices and edge midpoints on the primary mesh 𝓒h. Then, KP*, the dual cell associated with P𝓝h, is defined as the union of the subcells sharing P, and the dual mesh is defined as 𝓒h*={KP*:P𝓝h}. Now, the test function space is chosen as

Vh={vhL2(Ω̄):vhKP*=constant,P𝓝h;vhKP*=0,P𝓝hΩ}, where 𝓝h=𝓝h\Ω. Here we remark that the dual mesh 𝓒h* and the test function space Vh depend on two parameters β1 and β2. Different choices of the pair (β1,β2) lead to different finite volume element schemes.

images

Figure 2: Partition of the triangular element K

2.2 Semi-Discrete Schemes

In this part, we propose the following spatial discrete formulation of Eq. (1), Ω(ut-αΔutα-Δu)vhdxdy=Ω(f(u)+g)vhdxdy,vhVh.

By the divergence theorem and the definition of Vh, we have,

KP*utdxdy-KP*(αutαnP*+unP*)ds=KP*(f(u)+g)dxdy,P𝓝h,(6)

where is the gradient operator and nP* is the unit normal vector outward to KP*. The left hand-side of Eq. (6) can be rewritten as

KP*utdxdy-KP*(αutαnP*+unP*)ds=K𝓒h[KPutdxdy-εK,P(αutα+u)dsnP*], where KP=K KP* and εK,P=K KP*. Based on the above formulation, we formulate the semi-discrete finite volume element solution of the Eq. (6) as: Find uh=uh(x,y,t)Uh with t(0,T], such that

K𝓒h[KPuhtdxdy-εK,P(αuhtα+uh)dsnP*]=K𝓒hKP(f(uh)+g)dxdy,P𝓝h.(7)

For the computation of the terms in Eq. (7), we introduce the following affine mapping that transforms K onto K^ in (λ1,λ2) plane,

{λ1=1JK(|y21y31|x-|x21x31|y+|x2y2x3y3|),λ2=1JK(|y31y11|x-|x31x11|y+|x3y3x1y1|),(8)

where

JK=|x1y11x2y21x3y31|, and K^ is the reference element whose vertices and barycenter are specified by P^1=(1,0),P^2=(0,1),P^3=(0,0) and Q^=(1/3,1/3), respectively, see Fig. 3. Let uh,K be the restriction of uh on K. Then, under the affine mapping Eq. (8), the counterparts of uh and uh,K can be expressed as ûh=ûh(λ1,λ2,t),ûh,K^=j=16uj,Kϕ^j(λ1,λ2), where uj,K=uh(Pj),uj+3,K=uh(Mj),j=1,2,3,

images

Figure 3: The reference element K^ and its associated subcells

ϕ^j(λ1,λ2)(1j6) are the shape functions on K^, given by

{ϕ^1(λ1,λ2)=λ1(2λ1-1),ϕ^2(λ1,λ2)=λ2(2λ2-1),ϕ^3(λ1,λ2)=λ3(2λ3-1),ϕ^4(λ1,λ2)=4λ1λ2,ϕ^5(λ1,λ2)=4λ2λ3,ϕ^6(λ1,λ2)=4λ1λ3, and λ3=1-λ1-λ2. Now we rewrite Eq. (7) to get the following semi-discrete schemes,

K𝓒h[K^PûhtJKdλ1dλ2-ε^K,P(αûhtα+ûh)FK,PdlnP*]=K𝓒hK^P(f(ûh)+ĝ)JKdλ1dλ2,(9)

where P𝓝h and FK,P denotes the ratio of the measures of εK,P and ε^K,P. One can see that ûh,K^=j=16uj,Kϕ^j, where

ϕ^j=(ϕ^jxϕ^jy)=(ϕ^jλ1λ1x+ϕ^jλ2λ2xϕ^jλ1λ1y+ϕ^jλ2λ2y)=(λ1xλ2xλ1yλ2y)^ϕ^j, and ^ is the gradient operator with respect to (λ1,λ2) plane. By straightforward calculations, we have ^ϕ^1=(4λ1-10),^ϕ^2=(04λ2-1),^ϕ^3=(4λ1+4λ2-34λ1+4λ2-3),^ϕ^4=4(λ2λ1),^ϕ^5=4(-λ2-λ1-2λ2+1),^ϕ^6=4(-2λ1-λ2+1-λ1).

2.3 Fully Discrete Schemes

Next we introduce the fully discrete schemes at time tn-α/2. Let 0=t0<t1<<tN=T be a uniform partition of the time interval [0,T] with mesh length Δt=T/N, and tn=nΔt (n=0,1,2,,N). For any function z(t), let zn = z(tn).

Lemma 2.1. ([14], Lemma 2) Suppose z(t)C3[0,T]. Then, we have

zt(tn-α2)={1Δt(z1-z0)+𝓞(Δt),n=1,12Δt[(3-α)zn-(4-2α)zn-1+(1-α)zn-2]+𝓞(Δt2),n2.(10)

Lemma 2.2. ([11], Lemma 2) Assume that f(t)C2[0,T]. Then, the following second-order formula for the approximation of the nonlinear term at time tn-α/2 holds,

f(z(tn-α2))=(2-α2)f(zn-1)-(1-α2)f(zn-2)+𝓞(Δt2),n2.(11)

Lemma 2.3. ([15], Lemma 2) Suppose z(t)C3[0,T]. We have the following L2-1σ formula at time tn-α/2,

αztα(tn-α2)=Δt-αΓ(2-α)[c0nzn-j=1n-1(cn-j-1n-cn-jn)zj-cn-1nz0]+𝓞(Δt3-α),(12)

where c01=a0, for n = 1, and for n2,

cln={a0+b1,l=0,al+bl+1-bl,1ln-2,al-bl,l=n-1, with a0=(1-α2)1-α,al=(l+1-α2)1-α-(l-α2)1-α,l1,bl=12-α[(l+1-α2)2-α-(l-α2)2-α]-12[(l+1-α2)1-α+(l-α2)1-α],l1.

Based on Lemmas 2.1–2.3, we propose the following fully discrete schemes by the Eq. (9) at time tn-α/2:

for n = 1,

KTh{ K^P1Δt(u^h1u^h0)JKdλ1dλ2ε^K,P[ c01ΔtαΓ(2α)(u^h1u^h0)+((1α2)u^h1+α2u^h0) ]FK,PdlnP* }=KThK^P(f(u^h0)+g^1α2)JKdλ1dλ2,PNh°;(13)

for n2,

KTh{K^P12Δt((3α)u^hn(42α)u^hn1+(1α)u^hn2)JKdλ1dλ2ε^K,P[ ΔtαΓ(2α)(c0nu^hnj=1n1(cnj1ncnjn)u^hjcn1nu^h0)+((1α2)u^hn+α2u^hn1) ]×FK,PdlnP*}=KThK^P[ ((2α2)f(u^hn1)(1α2)f(u^hn2))+g^nα2 ]JKdλ1dλ2,PNh°,(14)

where ûhn=ûhn(λ1,λ2,tn), n=0,1,,N.

Let the basis functions of the trial function space Uh be denoted as φk(x,y), (k=1,2,,m) where m is the number of unknowns (i.e., m=#𝓝h), then the numerical solution uhn=(φ1,φ2,,φm)uhnUh, where the vector uhn=(u1n,u2n,,umn)T, (n=0,1,,N). Moreover, we know ûhn=(φ^1,φ^2,,φ^m)uhn, where φ^k(λ1,λ2), (k=1,2,,m) are transformed from φk(x,y), (k=1,2,,m) by the affine mapping Eq. (8). By simplifying and synthesizing the Eqs. (13), (14), we obtain the following matrix form of the fully discrete schemes:

for n = 1,

Muh1-[c01Δt1-αΓ(2-α)+Δt(1-α2)]Auh1=ΔtMf0+Δtg1-α2+Muh0-c01Δt1-αΓ(2-α)Auh0+Δtα2Auh0;(15)

for n2,

(3α)Muhn[ 2c0nΔt1αΓ(2α)+2Δt(1α2) ]Auhn=2ΔtM[ (2α2)fn1(1α2)fn2 ]+2Δtgnα2+(42α)Muhn1(1α)Muhn22Δt1αΓ(2α)A[ j=1n1(cnj1ncnjn)uhj+cn1nuh0 ]+ΔtαAuhn1,(16)

where

fn=(fkn)m×1,fkn=f(uhn(Pk)),PkNh°, (17)

gnα2=(gknα2)m×1,gknα2=KThK^Pkg^nα2JKdλ1dλ2,PkNh°, (18)

mass matrix M=K𝓒hMK and stiffness matrix A=K𝓒hAK. Meanwhile, MK and AK, square matrices of degree m, are expanded by the element matrices MK and AK as follows:

MK=(mij)6×6,mij=K^Piϕ^jJKdλ1dλ2, (19)

AK=(aij)6×6,aij=ε^K,Piϕ^jFK,PdlnPi*, (20)

where KPi=K KPi* and εK,Pi=K KPi* are transformed to K^Pi and ε^K,Pi by the affine mapping Eq. (8) and Pi (i=1,,6) (Pj+3=Mj,j=1,2,3) are vertices or edge midpoints of the element K.

For the above finite volume element schemes, we emphasize that the nonlinear term is approximated by using the linearized difference scheme in Lemma 2.2, and none nonlinear iteration is involved. The whole algorithm is summarized below.

images

3  Numerical Examples

In this section, we use Eqs. (15)(20) to solve four examples on uniform triangular mesh (Mesh I) and random triangular mesh (Mesh II), respectively, see Figs. 4, 5. The first level of Mesh II is constructed from Mesh I by the following random distortion of the interior vertices,

x:=x+ωηxh,y:=y+ωηyh, where ω(0,0.5) is the disturbance coefficient, h is the spatial mesh size, ηx and ηy are two random numbers located in [-1,1]. The subsequent level is refined by the standard bisection procedure from its previous level. In our numerical examples, we choose the disturbance coefficient ω=0.2 for Mesh II (see Fig. 5). Recall that our quadratic finite volume element schemes have two parameters β1 and β2. Here we just investigate the following specific schemes:

•   First scheme (QFVE-1): β1=(3-3)/6, β2=(6+3-21+63)/9;

•   Second scheme (QFVE-2): β1=β2=(3-3)/6;

•   Third scheme (QFVE-3): β1=1/4, β2=1/3;

•   Fourth scheme (QFVE-4): β1=β2=1/3.

images

Figure 4: Uniform triangular mesh (Mesh I)

images

Figure 5: Random triangular mesh (Mesh II)

We remark that the counterparts of the above schemes for elliptic problems have been studied in [42,43,46,47], respectively. The L2 errors Eu and convergence orders Ru1, Ru2 are defined as Eu= max1nN(K𝓒hK(uh(x,y,tn)-u(x,y,tn))2dxdy)12,Ru1=log(Eu(h2)/Eu(h1))log(h2/h1),Ru2=log(Eu(Δt2)/Eu(Δt1))log(Δt2/Δt1), respectively, where h1 and h2 are the spatial mesh sizes of two successive meshes, and Δt1 and Δt2 are the mesh sizes of two successive time levels. Moreover, the results of the quadratic finite element scheme (QFE) [29] are also employed for comparison.

3.1 Example 1

Solve Eqs. (1)(3) with the nonlinear term f(u) = u2 and the source term g(x,t)=[(3+α)t2+α+8π2t3+α+4π2t33Γ(4+α)]sin(2πx)sin(2πy)-t6+2α sin2(2πx)sin2(2πy), where Ω̄=[0,1]2 and T = 1. The exact solution is u(x,t)=t3+α sin(2πx)sin(2πy). We choose Δt=1/2000, h = 1/10 and α=0.1. Fig. 6 shows the exact solution u and the vertex values of the numerical solution uh at x = 0.3 and t = 1 on Mesh I, which implies that the numerical solution can well approximate the exact solution even though the five schemes exhibit difference in accuracy. Tabs. 1, 2 gives some detailed results for the quadratic finite volume element schemes and quadratic finite element scheme. In Figs. 7, 8, we draw log-log plots of the L2 errors vs. with mesh size h on Mesh I and Mesh II when Δt=1/2000. One can see that spatial convergence orders of the QFE and QFVE-1 schemes are close to 3 for different α, while the spatial convergence order of other quadratic finite volume element schemes (QFVE-2, QFVE-3, QFVE-4) are nearly 2 and lower than that of QFVE-1, which agrees with the observations made in [49] for elliptic problems. Moreover, we calculate temporal convergence orders of these schemes in Tabs. 3, 4, and we can find the temporal convergence orders are nearly 2.

images

Figure 6: Comparison of the numerical solutions and the exact solution for Example 1 on Mesh I. (a) A full profile; (b) A local enlarged profile

Table 1: Error results and spatial convergence orders with Δt=1/2000 on Mesh I in Example 1

images

Table 2: Error results and spatial convergence orders with Δt=1/2000 on Mesh II in Example 1

images

images

Figure 7: L2 errors of the numerical solution on Mesh I in Example 1. (a) α=0.1; (b) α=0.2; (c) α=0.5; (d) α=0.9

images

Figure 8: L2 errors of the numerical solution on Mesh II in Example 1. (a) α=0.1; (b) α=0.2; (c) α=0.5; (d) α=0.9

Table 3: Error results and temporal convergence orders with h = 1/160 on Mesh I in Example 1

images

Table 4: Error results and temporal convergence orders with h = 1/160 on Mesh II in Example 1

images

3.2 Example 2

Solve Eqs. (1)(3) with Ω̄=[0,2]2, T = 2, f(u)= sin(u) and g(x,t)=3t2g1(x)-12[t3+6t3-αΓ(4-α)]g2(x)- sin(t3g1(x)), where g1(x)=x(x-0.5)(x-1.5)(x-2)y(y-0.5)(y-1.5)(y-2), g2(x)=x(x-0.5)(x-1.5)(x-2)(y2-2y+1924)+y(y-0.5)(y-1.5)(y-2)(x2-2x+1924).

Table 5: Error results and spatial convergence orders with Δt=1/2000 on Mesh I in Example 2

images

Table 6: Error results and spatial convergence orders with Δt=1/2000 on Mesh II in Example 2

images

images

Figure 9: L2 errors of the numerical solution on Mesh II in Example 2. (a) QFE; (b) QFVE-1; (c) QFVE-2; (d) QFVE-3; (e) QFVE-4

images

Figure 10: Contour plots of |uuh| on Mesh I in Example 3. (a) QFE; (b) QFVE-1; (c) QFVE-2; (d) QFVE-3; (e) QFVE-4

Table 7: Numerical results with Δt=1/2000 on Mesh I in Example 3

images

Table 8: Numerical results with Δt=1/2000 on Mesh II in Example 3

images

Table 9: Numerical results with Δt=1/100 on Mesh I in Example 4

images

The exact solution to this example is:

u(x,t)=t3x(x-0.5)(x-1.5)(x-2)y(y-0.5)(y-1.5)(y-2).

Analogously, we calculate the L2 errors and spatial convergence orders for several kinds of quadratic finite volume element schemes and quadratic finite element scheme, see Tabs. 5, 6. The convergence behavior is similar to that for Example 1. Moreover, when α=0.1, 0.2, 0.5, 0.9 we draw the log–log plots for these five schemes on Mesh II, and find that the converge orders don’t change with α, see Fig. 9.

3.3 Example 3

We take the space-time domain Ω̄×[0,T]=[0,1]2×[0,1], the nonlinear term f(u) = u3u and the source term

g(x,t)=[2et+8π2et+8π2t1-αE1,2-α(t)]sin(2πx)sin(2πy)-e3t sin3(2πx)sin3(2πy), where Mittag–Leffler function E1,2-α(t) is defined by E1,2-α(t)=i=0tiΓ(i+2-α).

In numerical calculation of this example, we use i=01000ti/Γ(i+2-α) to approximate the Mittag–Leffler function E1,2-α(t). The exact solution to the model is u(x,t)=et sin(2πx)sin(2πy). In Fig. 10, we draw contour plots for the absolute value of error between exact solution u and numerical solution uh, i.e., |uuh|, with α=0.5, h = 1/10, Δt=1/2000 at time T = 1 on Mesh I. It is obvious that accuracy of the first three schemes (QFE, QFVE-1, QFVE-2) is better than that of QFVE-3 and QFVE-4, see Tabs. 7, 8 for some detailed data.

3.4 Example 4

In the last example, we choose the nonlinear term f(u) = u3u and the source term g(x, t) = 0 with initial condition u0(x)=x(1-x3)y(1-y3), where Ω̄×[0,T]=[0,1]2×[0,1]. Because of unknown exact solution, we take the numerical solution with h = 1/160, Δt=1/100 as the ‘exact’ solution when computing the errors. The results are given in Tab. 9 where one can see that these schemes still work in this situation.

4  Conclusions

In this article, we study a nonlinear time-fractional Rayleigh-Stokes problem by using the quadratic finite volume element method combined with a specific time discretization. In temporal direction, we use a two step scheme to approximate the equation at time tn-α/2, where L2-1σ formula is used to approximate the time-fractional derivative. The fully discrete schemes of quadratic finite volume element are suggested and we find that only one of these schemes achieves the optimal convergence order in L2 norm in space direction. We calculate some numerical examples by several kinds of quadratic finite volume element schemes and quadratic finite element scheme, space L2 error orders of the QFE and QFVE-1 schemes reach 3. Meanwhile, numerical results of other three quadratic finite volume element schemes (QFVE-2, QFVE-3, QFVE-4) are nearly 2 and lower than the optimal order of QFVE-1. The future work includes the stability analysis and error estimates by following the related results on elliptic problems [46,48,54].

Acknowledgement: The authors would like to thank the editor and the anonymous reviewers for their valuable suggestions.

Funding Statement: This work was partially supported by the National Natural Science Foundation of China (No. 11871009).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

 1.  Ma, J. T., Liu, J. Q., Zhou, Z. Q. (2014). Convergence analysis of moving finite element methods for space fractional differential equations. Journal of Computational and Applied Mathematics, 255(285), 661–670. DOI 10.1016/j.cam.2013.06.021. [Google Scholar] [CrossRef]

 2.  Bu, W. P., Tang, Y. F., Yang, J. Y. (2014). Galerkin finite element method for two-dimensional Riesz space fractional diffusion equations. Journal of Computational Physics, 276, 26–38. DOI 10.1016/j.jcp.2014.07.023. [Google Scholar] [CrossRef]

 3.  Wang, H., Du, N. (2013). A superfast-preconditioned iterative method for steady-state space-fractional diffusion equations. Journal of Computational Physics, 240(E), 49–57. DOI 10.1016/j.jcp.2012.07.045. [Google Scholar] [CrossRef]

 4.  Liu, F., Zhuang, P., Turner, I., Burrage, K., Anh, V. (2014). A new fractional finite volume method for solving the fractional diffusion equation. Applied Mathematical Modelling, 38(15–16), 3871–3878. DOI 10.1016/j.apm.2013.10.007. [Google Scholar] [CrossRef]

 5.  Feng, L. B., Zhuang, P., Liu, F., Turner, I. (2015). Stability and convergence of a new finite volume method for a two-sided space-fractional diffusion equation. Applied Mathematics and Computation, 257(1), 52–65. DOI 10.1016/j.amc.2014.12.060. [Google Scholar] [CrossRef]

 6.  Pan, J. Y., Ng, M. K., Wang, H. (2016). Fast iterative solvers for linear systems arising from time-dependent space-fractional diffusion equations. SIAM Journal on Scientific Computing, 38(5), 2806–2826. DOI 10.1137/15M1030273. [Google Scholar] [CrossRef]

 7.  Li, C., Zhao, S. (2016). Efficient numerical schemes for fractional water wave models. Computers & Mathematics with Applications, 71(1), 238–254. DOI 10.1016/j.camwa.2015.11.018. [Google Scholar] [CrossRef]

 8.  Cheng, X. J., Duan, J. Q., Li, D. F. (2019). A novel compact ADI scheme for two-dimensional Riesz space fractional nonlinear reaction-diffusion equations. Applied Mathematics and Computation, 346(6), 452–464. DOI 10.1016/j.amc.2018.10.065. [Google Scholar] [CrossRef]

 9.  Yin, B. L., Liu, Y., Li, H., He, S. (2019). Fast algorithm based on TT-M FE system for space fractional Allen-Cahn equations with smooth and non-smooth solutions. Journal of Computational Physics, 379(1), 351–372. DOI 10.1016/j.jcp.2018.12.004. [Google Scholar] [CrossRef]

10. Yazdani, A., Mojahed, N., Babaei, A., Cendon, E. V. (2020). Using finite volume-element method for solving space fractional advection-dispersion equation. Progress in Fractional Differentiation and Applications, 6(1), 55–66. DOI 10.18576/pfda/060106. [Google Scholar] [CrossRef]

11. Wang, Y. J., Liu, Y., Li, H., Wang, J. F. (2016). Finite element method combined with second-order time discrete scheme for nonlinear fractional Cable equation. European Physical Journal Plus, 131(61), 1–16. DOI 10.1140/epjp/i2016–16061-3. [Google Scholar] [CrossRef]

12. Liu, Y., Yu, Z. D., Li, H., Liu, F. W., Wang, J. F. (2018). Time two-mesh algorithm combined with finite element method for time fractional water wave model. International Journal of Heat and Mass Transfer, 120(1), 1132–1145. DOI 10.1016/j.ijheatmasstransfer.2017.12.118. [Google Scholar] [CrossRef]

13. Yin, B. L., Liu, Y., Li, H. (2020). A class of shifted high-order numerical methods for the fractional mobile/immobile transport equations. Applied Mathematics and Computation, 368(3), 124799. DOI 10.1016/j.amc.2019.124799. [Google Scholar] [CrossRef]

14. Gao, G. H., Sun, H. W., Sun, Z. Z. (2015). Stability and convergence of finite difference schemes for a class of time-fractional sub-diffusion equations based on certain superconvergence. Journal of Computational Physics, 280, 510–528. DOI 10.1016/j.jcp.2014.09.033. [Google Scholar] [CrossRef]

15. Alikhanov, A. A. (2015). A new difference scheme for the time fractional diffusion equation. Journal of Computational Physics, 280(1), 424–438. DOI 10.1016/j.jcp.2014.09.031. [Google Scholar] [CrossRef]

16. Lin, Y. M., Li, X. J., Xu, C. J. (2011). Finite difference/spectral approximations for the fractional Cable equation. Mathematics of Computation, 80(275), 1369–1396. DOI 10.1090/S0025-5718-2010-02438-X. [Google Scholar] [CrossRef]

17. Liu, F., Yang, C., Burrage, K. (2009). Numerical method and analytical technique of the modified anomalous subdiffusion equation with a nonlinear source term. Journal of Computational and Applied Mathematics, 231(1), 160–176. DOI 10.1016/j.cam.2009.02.013. [Google Scholar] [CrossRef]

18. Cao, X. N., Cao, X. X., Wen, L. P. (2016). The implicit midpoint method for the modified anomalous sub-diffusion equation with a nonlinear source term. Journal of Computational and Applied Mathematics, 318, 199–210. DOI 10.1016/j.cam.2016.10.014. [Google Scholar] [CrossRef]

19. Tuan, N. H., Zhou, Y., Thach, T. N., Can, N. H. (2019). Initial inverse problem for the nonlinear fractional Rayleigh-Stokes equation with random discrete data. Communications in Nonlinear Science and Numerical Simulation, 78(5), 104873. DOI 10.1016/j.cnsns.2019.104873. [Google Scholar] [CrossRef]

20. Zhou, Y., Wang, J. N. (2019). The nonlinear Rayleigh-Stokes problem with Riemann–Liouville fractional derivative. Mathematical Methods in the Applied Sciences, 1–8(3), 2431–2438. DOI 10.1002/mma.5926. [Google Scholar] [CrossRef]

21. Guan, Z., Wang, X. D., Ouyang, J. (2020). An improved finite difference/finite element method for the fractional Rayleigh-Stokes problem with a nonlinear source term. Journal of Applied Mathematics and Computing, 65, 1–29. DOI 10.1007/s12190-020-01399-4. [Google Scholar] [CrossRef]

22. Bao, N. T., Hoang, L. N., van, A. V., Nguyen, H. T., Zhou, Y. (2020). Existence and regularity of inverse problem for the nonlinear fractional Rayleigh-Stokes equations. Mathematical Methods in the Applied Sciences, 44(1), 1–27. DOI 10.1002/mma.6162. [Google Scholar] [CrossRef]

23. Sayevand, K., Arjang, F. (2016). Finite volume element method and its stability analysis for analyzing the behavior of sub-diffusion problems. Applied Mathematics and Computation, 290(2), 224–239. DOI 10.1016/j.amc.2016.06.008. [Google Scholar] [CrossRef]

24. Karaa, S., Mustapha, K., Pani, A. K. (2016). Finite volume element method for two-dimensional fractional subdiffusion problems. IMA Journal of Numerical Analysis, 37, 945–964. DOI 10.1093/imanum/drw010. [Google Scholar] [CrossRef]

25. Karaa, S., Pani, A. K. (2018). Error analysis of a finite volume element method for fractional order evolution equations with nonsmooth initial data. ESAIM Mathematical Modelling and Numerical Analysis, 52(2), 773–801. DOI 10.1051/m2an/2018029. [Google Scholar] [CrossRef]

26. Badr, M., Yazdani, A., Jafari, H. (2018). Stability of a finite volume element method for the time-fractional advection-diffusion equation. Numerical Methods for Partial Differential Equations, 34(5), 1459–1471. DOI 10.1002/num.22243. [Google Scholar] [CrossRef]

27. Zhao, J., Li, H., Fang, Z. C., Liu, Y. (2019). A mixed finite volume element method for time-fractional reaction-diffusion equations on triangular grids. Mathematics, 7(7), 600. DOI 10.3390/math7070600. [Google Scholar] [CrossRef]

28. Zhao, J., Fang, Z. C., Li, H., Liu, Y. (2020). Finite volume element method with the WSGD formula for nonlinear fractional mobile/immobile transport equations. Advances in Difference Equations, 360, 1–20. DOI 10.1186/s13662–020-02786-8. [Google Scholar] [CrossRef]

29. Zhang, Y. L., Yin, B. L., Cao, Y., Liu, Y., Li, H. (2020). A numerical algorithm based on quadratic finite element for two-dimensional nonlinear time fractional thermal diffusion model. Computer Modeling in Engineering & Sciences, 122(3), 1081–1098. DOI 10.32604/cmes.2020.07822. [Google Scholar] [CrossRef]

30. Wu, G. C. (2011). A fractional characteristic method for solving fractional partial differential equations. Applied Mathematics Letters, 24(7), 1046–1050. DOI 10.1016/j.aml.2011.01.020. [Google Scholar] [CrossRef]

31. Li, C. P., Zhao, Z. G., Chen, Y. Q. (2011). Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion. Computers & Mathematics with Applications, 62(3), 855–875. DOI 10.1016/j.camwa.2011.02.045. [Google Scholar] [CrossRef]

32. Feng, L. B., Zhuang, P., Liu, F., Turner, I., Gu, Y. T. (2015). Finite element method for space-time fractional diffusion equation. Numerical Algorithms, 72(3), 749–767. DOI 10.1007/s11075-015-0065-8. [Google Scholar] [CrossRef]

33. Fan, W. P., Liu, F. W., Jiang, X. Y., Turner, I. (2017). A novel unstructured mesh finite element method for solving the time-space fractional wave equation on a two-dimensional irregular convex domain. Fractional Calculus and Applied Analysis, 20(2), 352–383. DOI 10.1515/fca-2017-0019. [Google Scholar] [CrossRef]

34. Zhang, H., Jiang, X. Y., Yang, X. (2018). A time-space spectral method for the time-space fractional Fokker–Planck equation and its inverse problem. Applied Mathematics and Computation, 320(1), 302–318. DOI 10.1016/j.amc.2017.09.040. [Google Scholar] [CrossRef]

35. Fetecau, C., Zierep, J. (2001). On a class of exact solutions of the equations of motion of a second grade fluid. Acta Mechanica, 150(1–2), 135–138. DOI 10.1007/BF01178551. [Google Scholar] [CrossRef]

36. Tan, W. C., Masuoka, T. (2005). Stokes’ first problem for a second grade fluid in a porous half-space with heated boundary. International Journal of Non-Linear Mechanics, 40(4), 515–522. DOI 10.1016/j.ijnonlinmec.2004.07.016. [Google Scholar] [CrossRef]

37. Shen, F., Tan, W. C., Zhao, Y. H., Masuoka, T. (2006). The Rayleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative model. Nonlinear Analysis: Real World Applications, 7(5), 1072–1080. DOI 10.1016/j.nonrwa.2005.09.007. [Google Scholar] [CrossRef]

38. Zierep, J., Fetecau, C. (2007). Energetic balance for the Rayleigh-Stokes problem of a Maxwell fluid. International Journal of Engineering Science, 45(2), 617–627. DOI 10.1016/j.ijengsci.2007.04.015. [Google Scholar] [CrossRef]

39. Parvizi, M., Khodadadian, A., Eslahchi, M. (2020). Analysis of Ciarlet-Raviart mixed finite element methods for solving damped Boussinesq equation. Journal of Computational and Applied Mathematics, 379(1), 112818. DOI 10.1016/j.cam.2020.112818. [Google Scholar] [CrossRef]

40. Abbaszadeh, M., Dehghan, M., Khodadadian, A., Heitzinger, C. (2019). Analysis and application of the interpolating element free Galerkin (IEFG) method to simulate the prevention of groundwater contamination with application in fluid flow. Journal of Computational and Applied Mathematics, 368(155), 112453. DOI 10.1016/j.cam.2019.112453. [Google Scholar] [CrossRef]

41. Abbaszadeh, M., Dehghan, M., Khodadadian, A., Noii, N., Heitzinger, C. et al. (2020). A reduced-order variational multiscale interpolating element free Galerkin technique based on proper orthogonal decomposition for solving Navier–Stokes equations coupled with a heat transfer equation: Nonstationary incompressible Boussinesq equations. Journal of Computational Physics, 426(2), 109875. DOI 10.1016/j.jcp.2020.109875. [Google Scholar] [CrossRef]

42. Tian, M. Z., Chen, Z. Y. (1991). Quadratic element generalized differential methods for elliptic equations. Numerical Mathematics a Journal of Chinese Universities, 13, 99–113. DOI http://dx.doi.org/CNKI: SUN: GDSX.0.1991-02-000. [Google Scholar]

43. Liebau, F. (1996). The finite volume element method with quadratic basis functions. Computing, 57(4), 281–299. DOI 10.1007/BF02252250. [Google Scholar] [CrossRef]

44. Xu, J. C., Zou, Q. S. (2009). Analysis of linear and quadratic simplicial finite volume methods for elliptic equations. Numerische Mathematik, 111(3), 469–492. DOI 10.1007/s00211-008-0189-z. [Google Scholar] [CrossRef]

45. Chen, Z. Y., Wu, J. F., Xu, Y. S. (2012). Higher-order finite volume methods for elliptic boundary value problems. Advances in Computational Mathematics, 37(2), 191–253. DOI 10.1007/s10444-011-9201-8. [Google Scholar] [CrossRef]

46. Wang, X., Li, Y. H. (2016). L2 error estimates for high order finite volume methods on triangular meshes. SIAM Journal on Numerical Analysis, 54(5), 2729–2749. DOI 10.1137/140988486. [Google Scholar] [CrossRef]

47. Zou, Q. S. (2017). An unconditionally stable quadratic finite volume scheme over triangular meshes for elliptic equations. Journal of Scientific Computing, 70(1), 112–124. DOI 10.1007/s10915-016-0244-3. [Google Scholar] [CrossRef]

48. Zhou, Y. H., Wu, J. M. (2020). A family of quadratic finite volume element schemes over triangular meshes for elliptic equations. Computers & Mathematics with Applications, 79(9), 2473–2491. DOI 10.1016/j.camwa.2019.11.017. [Google Scholar] [CrossRef]

49. Zhou, Y. H., Wu, J. M. (2020). A unified analysis of a class of quadratic finite volume element schemes on triangular meshes. Advances in Computational Mathematics, 46(5), 777. DOI 10.1007/s10444-020-09809-8. [Google Scholar] [CrossRef]

50. Wang, P., Zhang, Z. Y. (2010). Quadratic finite volume element method for the air pollution model. International Journal of Computer Mathematics, 87(13), 2925–2944. DOI 10.1080/00207160802680663. [Google Scholar] [CrossRef]

51. Jin, G. H., Li, H. G., Zhang, Q. H., Zou, Q. S. (2016). Linear and quadratic finite volume methods on triangular meshes for elliptic equations with singular solutions. International Journal of Numerical Analysis and Modeling, 13(2), 244–264. [Google Scholar]

52. Xiong, Z. G., Deng, K. (2017). A quadratic triangular finite volume element method for a semilinear elliptic equation. Advances in Applied Mathematics and Mechanics, 9(1), 186–204. DOI 10.4208/aamm.2014.m63. [Google Scholar] [CrossRef]

53. Du, Y. W., Li, Y. H., Sheng, Z. Q. (2019). Quadratic finite volume method for a nonlinear elliptic problem. Advances in Applied Mathematics and Mechanics, 11(4), 838–869. DOI 10.4208/aamm.OA-2017-0231. [Google Scholar] [CrossRef]

54. Zhou, Y. H. (2020). A class of bubble enriched quadratic finite volume element schemes on triangular meshes. International Journal of Numerical Analysis and Modeling, 17(6), 872–899. [Google Scholar]

images This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.