Next Article in Journal
Fractals as Julia Sets of Complex Sine Function via Fixed Point Iterations
Next Article in Special Issue
Initial Value Problems of Fuzzy Fractional Coupled Partial Differential Equations with Caputo gH-Type Derivatives
Previous Article in Journal
Fractional-Order Logistic Differential Equation with Mittag–Leffler-Type Kernel
Previous Article in Special Issue
A Fast Preconditioned Semi-Implicit Difference Scheme for Strongly Nonlinear Space-Fractional Diffusion Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mixed Element Algorithm Based on the Modified L1 Crank–Nicolson Scheme for a Nonlinear Fourth-Order Fractional Diffusion-Wave Model

1
School of Statistics and Mathematics, Inner Mongolia University of Finance and Economics, Hohhot 010070, China
2
School of Mathematical Sciences, Inner Mongolia University, Hohhot 010021, China
*
Authors to whom correspondence should be addressed.
Submission received: 15 October 2021 / Revised: 4 December 2021 / Accepted: 10 December 2021 / Published: 14 December 2021
(This article belongs to the Special Issue Novel Numerical Solutions of Fractional PDEs)

Abstract

:
In this article, a new mixed finite element (MFE) algorithm is presented and developed to find the numerical solution of a two-dimensional nonlinear fourth-order Riemann–Liouville fractional diffusion-wave equation. By introducing two auxiliary variables and using a particular technique, a new coupled system with three equations is constructed. Compared to the previous space–time high-order model, the derived system is a lower coupled equation with lower time derivatives and second-order space derivatives, which can be approximated by using many time discrete schemes. Here, the second-order Crank–Nicolson scheme with the modified L 1 -formula is used to approximate the time direction, while the space direction is approximated by the new MFE method. Analyses of the stability and optimal L 2 error estimates are performed and the feasibility is validated by the calculated data.

1. Introduction

Fourth-order fractional partial differential equations (PDEs) including fourth-order fractional subdiffusion models [1,2,3] and fourth-order fractional diffusion-wave models [2,4,5] can be founded in many fields of science and engineering. Thus far, there have been many efficient numerical algorithms for solving linear or nonlinear fourth-order fractional subdiffusion and diffusion-wave models. Liu et al. [6], Liu et al. [7], and Liu et al. [8] considered different mixed element methods to solve fourth-order nonlinear fractional subdiffusion models with the first-order time derivative and developed numerical theories including stability and convergence. Liu et al. [3] introduced a mixed element algorithm with a new approximation of the fractional derivative. Ji et al. [9], Ran et al. [10], Nandal and Pandey [11], Sun et al. [12], and Huang et al. [13] considered some difference schemes for linear or nonlinear fourth-order fractional diffusion or diffusion-wave models. Abbaszadeh and Dehghan [14] studied the direct meshless local Petrov–Galerkin method for solving fourth-order reaction-diffusion problems with a time-fractional derivative. Yang et al. [15] and Zhang et al. [16] found the numerical solutions for a fourth-order fractional model by using the orthogonal spline collocation method. Tariq and Akram [17] considered a quintic spline technique to solve a fourth-order time-fractional subdiffusion model. Guo et al. [18] and Du et al. [19] studied the LDG methods for solving some time-fractional subdiffusion models with fourth-order spatial derivative terms, respectively. In [1], Nikan et al. developed a local radial basis function generated by the finite difference scheme for a time-fractional fourth-order reaction-diffusion model. In [5], Jafari et al. solved a fourth-order fractional diffusion-wave equation by the decomposition method. Hu and Zhang [2] implemented numerical calculations via finite difference methods for fourth-order time fractional subdiffusion and diffusion-wave models. Li and Wong [20] developed an efficient numerical algorithm for a fourth-order time-fractional diffusion-wave model.
Here, we propose a new mixed element algorithm to solve the following nonlinear fourth-order time-fractional diffusion-wave model:
2 u t 2 + R L β u t β + u t + 2 u f ( u ) = g ( x , t ) , ( x , t ) Ω × J , u ( x , t ) = Δ u ( x , t ) = 0 , t J ¯ , u ( x , 0 ) = 0 , u t ( x , 0 ) = u 1 ( x ) , x Ω ¯ ,
where Ω R d ( d 2 ) and J = ( 0 , T ] with 0 < T < are the spatial domain and time interval, respectively. u 1 ( x ) is an initial value function, g ( x , t ) is a given source term, f ( u ) is a polynomial function or bounded function on u satisfying f C 2 ( R ) , and the Riemann–Liouville fractional derivative is defined by
R L β u t β = 1 Γ ( 2 β ) 2 t 2 0 t u ( s ) d s ( t s ) β 1 , 1 < β < 2 ,
where the nonlinear fourth-order fractional diffusion-wave model (1) can be generated by the classical fourth-order hyperbolic wave equation. When β 1 or 2 and f ( u ) = u 3 u , the model (1) can be reduced to an important Cahn–Hilliard equation model [21].
Recently, Zeng and Li [22] developed a new Crank–Nicolson scheme based on a modified L 1 -formula, whose coefficients are different from the famous L 1 -formula (see [23,24] for the fractional parameter α ( 0 , 1 ) ). One should note that this modified L 1 -formula can only approximate the Caputo or Riemann–Liouville fractional derivative with parameter α ( 0 , 1 ) , and it cannot approximate the case β ( 1 , 2 ) . Here, we will develop the modified L 1 -formula for the case of β ( 1 , 2 ) by using some techniques.
In this article, by introducing two auxiliary functions and using some techniques, we propose a new mixed element algorithm. Here, our major contributions are as follows: (1) by the introduction of two auxiliary functions, we reduce the nonlinear fourth-order time-fractional diffusion-wave model to a low-order coupled system; (2) we turn order β ( 1 , 2 ) into order α ( 0 , 1 ) for the Riemann–Liouville fractional derivative; (3) we approximate the derived coupled system with a fractional derivative with order α ( 0 , 1 ) by the modified L 1 Crank–Nicolson scheme with the developed new mixed element method; (4) we derive the stability of the new mixed element scheme and optimal error estimates in the L 2 -norm for three functions.
The structure of this article is as follows: in Section 2, we provide some numerical approximation formulas, propose a new mixed element scheme, and prove the stability of the derived scheme; in Section 3, we derive optimal error estimates for three variables; in Section 4, some numerical data are computed and discussed; Finally, in Section 5, we give some concluding remarks.

2. Numerical Approximation and Stability

Based on the relation between the Riemann–Liouville fractional derivative and Caputo fractional derivative, we take α = β 1 and v = u t to obtain
R L β u t β = 1 Γ ( 2 β ) 0 t 2 u ( s ) s 2 d s ( t s ) β 1 + u ( 0 ) Γ ( 1 β ) t β + u ( 0 ) t Γ ( 2 β ) t 1 β = 1 Γ ( 1 α ) 0 t v s d s ( t s ) α + v ( 0 ) Γ ( 1 α ) t α = R L α v t α , 0 < α < 1 .
Let σ = u f ( u ) ; (1) can be rewritten as the following coupled system:
v = u t , ( x , t ) Ω × J , σ t = v f u ( u ) v , ( x , t ) Ω × J , v t + R L α v t α + v + σ = g ( x , t ) , ( x , t ) Ω × J .
For formulating the fully discrete scheme, we insert the nodes t n = n Δ t ( n = 0 , 1 , 2 , , N ) in time interval [ 0 , T ] , where t n satisfy 0 = t 0 < t 1 < t 2 < < t N = T and the time step length size Δ t = T / N , for some positive integer N. For a smooth function ϕ defined on the time interval [ 0 , T ] , we denote ϕ n = ϕ ( t n ) .
Now, we need to introduce some lemmas on integer and fractional derivatives.
Lemma 1.
At t k + 1 2 , the following relation holds:
ϕ t ( t k + 1 2 ) = ϕ k + 1 ϕ k Δ t + O ( Δ t 2 ) P Δ t ϕ k + 1 2 + O ( Δ t 2 ) .
Lemma 2.
At t k + 1 2 , we have
ϕ ( t k + 1 2 ) = ϕ k + 1 + ϕ k 2 + O ( Δ t 2 ) ϕ k + 1 2 + O ( Δ t 2 ) .
Lemma 3.
([22]). At t k + 1 2 , the Caputo fractional derivative has the following form:
C α ϕ t α ( t k + 1 2 ) = 1 Γ ( 1 α ) 0 t k + 1 2 ϕ s d s ( t k + 1 2 s ) α = Δ t α Γ ( 2 α ) [ a 0 ϕ ( t k + 1 2 ) j = 1 k ( a k j a k j + 1 ) ϕ ( t j 1 2 ) ( a k b k ) ϕ ( t 1 2 ) b k ϕ ( t 0 ) ] + O ( Δ t 2 α ) ,
for k 0 , we have
b k = 2 [ ( k + 1 2 ) 1 α k 1 α ] , a k j = [ ( k j + 1 ) 1 α ( k j ) 1 α ] .
By Lemma 3, we have
Lemma 4.
([22]). At t k + 1 2 , the Riemann–Liouville fractional derivative has the following approximation:
R L α ϕ t α ( t k + 1 2 ) = Δ t α Γ ( 2 α ) [ a 0 ϕ ( t k + 1 2 ) j = 1 k ( a k j a k j + 1 ) ϕ ( t j 1 2 ) ( a k b k ) ϕ ( t 1 2 ) b ^ k ϕ ( t 0 ) ] + O ( Δ t 2 α ) ,
where b ^ k = b k ( 1 α ) ( k + 1 2 ) α .
Remark 1.
In [22], the authors provided the L 1 -formula above, which is different from the usual L 1 -formula and called the modified L 1 -formula.
By the approximation scheme above, we arrive at
( a ) P Δ t u n + 1 2 = v n + 1 2 + R 1 n + 1 2 , ( b ) P Δ t σ n + 1 2 = v n + 1 2 f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 + R 2 n + 1 2 , ( c ) P Δ t v n + 1 2 + Δ t α Γ ( 2 α ) [ a 0 v n + 1 2 j = 1 n ( a n j a n j + 1 ) v j 1 2 ( a n b n ) v 1 2 b ^ n v 0 ] + v n + 1 2 + σ n + 1 2 = g ( x , t n + 1 2 ) + R 3 n + 1 2 ,
where
R 1 n + 1 2 = P Δ t u n + 1 2 u t ( t n + 1 2 ) + ( v ( t n + 1 2 ) v n + 1 2 ) = O ( Δ t 2 ) , R 2 n + 1 2 = P Δ t σ n + 1 2 σ t ( t n + 1 2 ) + ( v ( t n + 1 2 ) v n + 1 2 ) + f u ( u ( t n + 1 2 ) ) v ( t n + 1 2 ) f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 = O ( Δ t 2 ) , R 3 n + 1 2 = P Δ t v n + 1 2 v t ( t n + 1 2 ) + O ( Δ t 2 α ) + ( v n + 1 2 v ( t n + 1 2 ) ) + ( σ n + 1 2 σ ( t n + 1 2 ) ) = O ( Δ t 2 α ) .
For ( φ , ψ , χ ) L 2 × H 0 1 × H 0 1 , we have the following mixed weak formulation:
( a ) ( P Δ t u n + 1 2 , φ ) = ( v n + 1 2 , φ ) + ( R 1 n + 1 2 , φ ) , ( b ) ( P Δ t σ n + 1 2 , ψ ) + ( v n + 1 2 , ψ ) + ( f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 , ψ ) = ( R 2 n + 1 2 , ψ ) , ( c ) ( P Δ t v n + 1 2 , χ ) + ( Δ t α Γ ( 2 α ) [ a 0 v n + 1 2 j = 1 n ( a n j a n j + 1 ) v j 1 2 ( a n b n ) v 1 2 b ^ n v 0 ] , χ ) + ( v n + 1 2 , χ ) ( σ n + 1 2 , χ ) = ( g ( x , t n + 1 2 ) , χ ) + ( R 3 n + 1 2 , χ ) .
For ( φ h , ψ h , χ h ) L h × V h × V h L 2 × H 0 1 × H 0 1 , based on the mixed weak formulation above, we formulate the following new mixed element system:
( a ) ( P Δ t u h n + 1 2 , φ h ) = ( v h n + 1 2 , φ h ) , ( b ) ( P Δ t σ h n + 1 2 , ψ h ) + ( v h n + 1 2 , ψ h ) + ( f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n 2 , ψ h ) = 0 , ( c ) ( P Δ t v h n + 1 2 , χ h ) + ( Δ t α Γ ( 2 α ) [ a 0 v h n + 1 2 j = 1 n ( a n j a n j + 1 ) v h j 1 2 ( a n b n ) v h 1 2 b ^ n v h 0 ] , χ h ) + ( v h n + 1 2 , χ h ) ( σ h n + 1 2 , χ h ) = ( g ( x , t n + 1 2 ) , χ h ) .
Lemma 5.
(See [22]). For b ^ n j , the following important inequality holds:
Δ t 1 α Γ ( 2 α ) j = 0 n 1 b ^ n j C T 1 α Γ ( 2 α ) ,
where C is a positive constant that is independent of space–time step length sizes h and Δ t .
Proof. 
Applying the Taylor formula, we have
Δ t 1 α Γ ( 2 α ) j = 0 n 1 b ^ n j = Δ t 1 α Γ ( 2 α ) j = 0 n 1 [ 2 ( n j + 1 2 ) 1 α 2 ( n j ) 1 α ( 1 α ) ( n j + 1 2 ) α ] = 2 Δ t 1 α Γ ( 2 α ) j = 0 n 1 ( n j ) 1 α [ ( 1 + 1 2 ( n j ) ) 1 α 1 1 α 2 ( n j + 1 2 ) ( 1 + 1 2 ( n j ) ) 1 α ] = Δ t 1 α Γ ( 2 α ) j = 0 n 1 ( n j ) 1 α [ 1 α 2 ( n j ) + ( 1 α ) α 2 ! ( 1 + κ 1 2 ( n j ) ) 1 α 1 4 ( n j ) 2 1 α 2 ( n j + 1 2 ) ( 1 + 1 α 2 ( n j ) + ( 1 α ) α 2 ! ( 1 + κ 1 2 ( n j ) ) 1 α 1 4 ( n j ) 2 ) ] = Δ t 1 α Γ ( 2 α ) j = 0 n 1 ( n j ) 1 α [ 1 α 2 ( n j ) 1 α 2 ( n j + 1 2 ) + O ( 1 ( n j ) 2 ) ] = Δ t 1 α Γ ( 2 α ) j = 0 n 1 ( n j ) 1 α [ 1 α 2 ( n j + 1 2 ) ( n j ) + O ( 1 ( n j ) 2 ) ] .
Noting that Δ t = T N and n j N , we have
Δ t 1 α Γ ( 2 α ) j = 0 n 1 ( n j ) 1 α [ 1 α 2 ( n j + 1 2 ) ( n j ) + O ( 1 ( n j ) 2 ) ] = 1 Γ ( 2 α ) j = 0 n 1 T 1 α ( n j N ) 1 α [ 1 α 2 ( n j + 1 2 ) ( n j ) + O ( 1 ( n j ) 2 ) ] T 1 α Γ ( 2 α ) j = 0 n 1 [ 1 ( n j ) 2 + O ( 1 ( n j ) 2 ) ] C T 1 α Γ ( 2 α ) n = 1 + 1 n 2 C T 1 α Γ ( 2 α ) .
Substitute (15) into (14) to obtain the conclusion. □
Next, we will prove the stability.
Theorem 1.
For n 0 , the stability for the fully discrete system (12) holds:
( a ) . v h n + 1 + σ h n + 1 + ( Δ t 1 α Γ ( 2 α ) j = 1 n + 1 a n j + 1 v h j 1 2 2 ) 1 2 C ( v h 0 + σ h 0 + max 0 j n { g ( x , t j + 1 2 ) } ) , ( b ) . u h n + 1 C ( u h 0 + v h 0 + σ h 0 + max 0 j n { g ( x , t j + 1 2 ) } ) .
Proof. 
In (12) ( a ) , we take φ h = u h n + 1 2 , and use Cauchy–Schwarz inequality as well as Young inequality to obtain
1 2 ( v h n + 1 2 2 + u h n + 1 2 2 ) ( v h n + 1 2 , u h n + 1 2 ) = ( P Δ t u h n + 1 2 , u h n + 1 2 ) u h n + 1 2 u h n 2 2 Δ t .
In (12) ( b ) , set ψ h = σ h n + 1 2 and make use of Cauchy–Schwarz inequality to arrive at
( v h n + 1 2 , σ h n + 1 2 ) = ( P Δ t σ h n + 1 2 , σ h n + 1 2 ) ( f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n 2 , σ h n + 1 2 ) σ h n + 1 2 σ h n 2 2 Δ t + 1 2 ( f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n ) σ h n + 1 2 σ h n + 1 2 σ h n + 1 2 2 Δ t + C ( v h n + 1 2 + v h n 2 + σ h n + 1 2 + σ h n 2 ) .
In (12) ( c ) , set χ h = v h n + 1 2 and use Cauchy–Schwarz inequality to obtain
( σ h n + 1 2 , v h n + 1 2 ) = ( P Δ t v h n + 1 2 , v h n + 1 2 ) ( Δ t α Γ ( 2 α ) [ a 0 v h n + 1 2 j = 1 n ( a n j a n j + 1 ) v h j 1 2 ( a n b n ) v h 1 2 b n v h 0 ] , v h n + 1 2 ) v h n + 1 2 2 + ( g ( x , t n + 1 2 ) , v h n + 1 2 ) v h n + 1 2 v h n 2 2 Δ t ( Δ t α Γ ( 2 α ) [ a 0 v h n + 1 2 j = 1 n ( a n j a n j + 1 ) v h j 1 2 ( a n b n ) v h 1 2 b ^ n v h 0 ] , v h n + 1 2 ) 1 2 v h n + 1 2 2 + 1 2 g ( x , t n + 1 2 ) 2 .
Add (18) and (19) to obtain
v h n + 1 2 v h n 2 2 Δ t + σ h n + 1 2 σ h n 2 2 Δ t + 1 2 v h n + 1 2 2 ( Δ t α Γ ( 2 α ) [ a 0 v h n + 1 2 j = 1 n ( a n j a n j + 1 ) v h j 1 2 ( a n b n ) v h 1 2 b ^ n v h 0 ] , v h n + 1 2 ) + C ( v h n + 1 2 + v h n 2 + σ h n + 1 2 + σ h n 2 + g ( x , t n + 1 2 ) 2 ) .
Refer to Lemma 4.2 in [22] to easily obtain
( Δ t α Γ ( 2 α ) [ a 0 v h n + 1 2 j = 1 n ( a n j a n j + 1 ) v h j 1 2 ( a n b n ) v h 1 2 b n v h 0 ] , v h n + 1 2 ) Δ t α 2 Γ ( 2 α ) ( j = 1 n a n j v h j 1 2 2 j = 1 n + 1 a n j + 1 v h j 1 2 2 + b ^ n v h 0 2 ) .
Combine (20) with (21) to obtain
v h n + 1 2 + σ h n + 1 2 + Δ t v h n + 1 2 2 + Δ t 1 α Γ ( 2 α ) j = 1 n + 1 a n j + 1 v h j 1 2 2 v h n 2 + σ h n 2 + Δ t 1 α Γ ( 2 α ) j = 1 n a n j v h j 1 2 2 + Δ t 1 α Γ ( 2 α ) b ^ n v h 0 2 + C Δ t ( v h n + 1 2 + v h n 2 + σ h n + 1 2 + σ h n 2 + g ( x , t n + 1 2 ) 2 ) .
We denote
Ξ ( v h n + 1 , σ h n + 1 ) = v h n + 1 2 + σ h n + 1 2 + Δ t 1 α Γ ( 2 α ) j = 1 n + 1 a n j + 1 v h j 1 2 2 .
Remove the non-negative term to obtain
Ξ ( v h n + 1 , σ h n + 1 ) ( 1 + Δ t 1 Δ t ) Ξ ( v h n , σ h n ) + 1 1 Δ t Δ t 1 α Γ ( 2 α ) b ^ n v h 0 2 + C Δ t 1 Δ t g ( x , t n + 1 2 ) 2 ( 1 + Δ t 1 Δ t ) 2 Ξ ( v h n 1 , σ h n 1 ) + 1 1 Δ t Δ t 1 α Γ ( 2 α ) v h 0 2 j = 0 1 b ^ n j ( 1 + Δ t 1 Δ t ) j   + C Δ t 1 Δ t j = 0 1 ( 1 + Δ t 1 Δ t ) j g ( x , t n j + 1 2 ) 2 ( 1 + Δ t 1 Δ t ) n Ξ ( v h 1 , σ h 1 ) + 1 1 Δ t Δ t 1 α Γ ( 2 α ) v h 0 2 j = 0 n 1 b ^ n j ( 1 + Δ t 1 Δ t ) j   + C Δ t 1 Δ t j = 0 n 1 ( 1 + Δ t 1 Δ t ) j g ( x , t n j + 1 2 ) 2 .
Noting that ( 1 + Δ t 1 Δ t ) > 1 , Δ t = T / N T / n , we have
( 1 + Δ t 1 Δ t ) n ( 1 + Δ t 1 Δ t ) n + 1 ( 1 + 2 Δ t 1 Δ t ) T Δ t lim Δ t 0 ( 1 + 2 Δ t 1 Δ t ) T ( 1 Δ t ) 2 Δ t 2 1 Δ t = e 2 .
Further, noting that b ^ n j > 0 and using Lemma 5, we have
    1 1 Δ t Δ t 1 α Γ ( 2 α ) v h 0 2 j = 0 n 1 b n j ( 1 + Δ t 1 Δ t ) j + C Δ t 1 Δ t j = 0 n 1 ( 1 + Δ t 1 Δ t ) j g ( x , t n j + 1 2 ) 2 e 2 1 Δ t Δ t 1 α Γ ( 2 α ) v h 0 2 j = 0 n 1 b ^ n j + C Δ t 1 Δ t j = 0 n 1 g ( x , t n j + 1 2 ) 2 C ( T 1 α Γ ( 2 α ) v h 0 2 + max 1 j n { g ( x , t j + 1 2 ) 2 } ) .
Substitute (25) and (26) into (24) to arrive at
Ξ ( v h n + 1 , σ h n + 1 ) C ( Ξ ( v h 1 , σ h 1 ) + T 1 α Γ ( 2 α ) v h 0 2 + max 1 j n { g ( x , t j + 1 2 ) 2 } ) .
Now, we estimate Ξ ( v h 1 , σ h 1 ) . Using (12) ( c ) , taking χ h = v h 1 2 , and using Cauchy–Schwarz inequality, we have
( σ h 1 2 , v h 1 2 ) = ( P Δ t v h 1 2 , v h 1 2 ) ( Δ t α Γ ( 2 α ) [ ( 1 2 ) α v h 1 2 α ( 1 2 ) α v h 0 ] , v h 1 2 ) v h 1 2 2 + ( g ( x , t 1 2 ) , v h 1 2 ) v h 1 2 v h 0 2 2 Δ t ( 1 2 ) α Δ t α Γ ( 2 α ) v h 1 2 2 + Δ t α Γ ( 2 α ) α ( 1 2 ) α ( v h 0 2 + v h 1 2 2 ) 1 2 v h 1 2 2 + 1 2 g ( x , t 1 2 ) 2 .
For n = 0 , we sum for (18) and (28) to obtain
v h 1 2 v h 0 2 2 Δ t + σ h 1 2 σ h 0 2 2 Δ t + ( 1 2 + 2 α Δ t α Γ ( 1 α ) ) v h 1 2 2 α 2 α Δ t α Γ ( 2 α ) v h 0 2 + C ( v h 1 2 + v h 0 2 + σ h 1 2 + σ h 0 2 + g ( x , t 1 2 ) 2 ) .
Noting that 1 α 2 α , ( 0 < α < 1 ) and (23), we have, for sufficiently small Δ t ,
Ξ ( v h 1 , σ h 1 ) = v h 1 2 + σ h 1 2 + Δ t 1 α Γ ( 2 α ) a 0 v h 1 2 2 v h 1 2 + σ h 1 2 + ( 1 + 2 Δ t 1 α Γ ( 1 α ) ) v h 1 2 2 C ( v h 0 2 + σ h 0 2 + g ( x , t 1 2 ) 2 ) .
Substitute (30) into (27) to obtain
Ξ ( v h n + 1 , σ h n + 1 ) C ( v h 0 2 + σ h 0 2 + max 0 j n { g ( x , t j + 1 2 ) 2 } ) , n 0 .
Combine (31) with (17) and use the Gronwall lemma to obtain
u h n + 1 2 C ( u h 0 2 + v h 0 2 + σ h 0 2 + max 0 j n { g ( x , t j + 1 2 ) 2 } ) , n 0 .
Using (31) and (32), we obtain the conclusion. □

3. A Priori Error Estimate

Now, we provide two projection operators [25] to derive a priori error estimates of our mixed finite element method.
Lemma 6.
Define the L 2 projection P h : L 2 ( Ω ) L h as
( u P h u , φ h ) = 0 , φ h L h ,
with the estimate inequality
u P h u + u t P h u t C h m + 1 u m + 1 , u L 2 ( Ω ) .
Lemma 7.
Define the elliptic projection Q h : H 0 1 ( Ω ) V h as
( ( v Q h v ) , ϕ h ) = 0 , ϕ h V h ,
with the following inequality:
v Q h v + v t Q h v t + h v Q h v 1 C h k + 1 ( v k + 1 + v t k + 1 ) , v H 0 1 ( Ω ) H k + 1 ( Ω ) .
In what follows, we derive the proof of error estimates in L 2 -norm in detail.
Theorem 2.
For P h u ( 0 ) = u h 0 , Q h v ( 0 ) = v h 0 and Q h σ ( 0 ) = σ h 0 , there exists a positive constant C that is independent of space–time step length sizes ( h , Δ t ) and we have for n 0
    u ( t n + 1 ) u h n + 1 + v ( t n + 1 ) v h n + 1 + σ ( t n + 1 ) σ h n + 1 C [ ( 1 + μ t n + 1 2 1 β ) h k + 1 + Δ t 3 β + h m + 1 ] ,
where, for the Caputo fractional derivative, we take μ as 0; for the Riemann–Liouville fractional derivative, we take μ as 1.
Proof. 
For convenience, we write
u ( t n ) u h n = ( u ( t n ) P h u n ) + ( P h u n u h n ) = E n + E n , v ( t n ) v h n = ( v ( t n ) Q h v n ) + ( Q h v n v h n ) = F n + F n , σ ( t n ) σ h n = ( σ ( t n ) Q h σ n ) + ( Q h σ n σ h n ) = H n + H n .
Applying triangle inequality, we have
u ( t n ) u h n E n + E n , v ( t n ) v h n F n + F n , σ ( t n ) σ h n H n + H n .
Using Lemmas 6 and 7, we arrive at the estimates of E n , F n , and H n . Consequently, in the discussion below, we only need to derive the estimates of E n , F n , and H n . Using projections (33) and (35), we have error equations as follows:
( a ) ( P Δ t E n + 1 2 , φ h ) = ( P Δ t E n + 1 2 , φ h ) + ( F n + 1 2 + F n + 1 2 , φ h ) + ( R 1 n + 1 2 , φ h ) , ( b ) ( P Δ t H n + 1 2 , ψ h ) + ( F n + 1 2 , ψ h ) = ( f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n 2 , ψ h ) ( P Δ t H n + 1 2 , ψ h ) + ( R 2 n + 1 2 , ψ h ) , ( c ) ( P Δ t F n + 1 2 , χ h ) + ( Δ t α Γ ( 2 α ) [ a 0 F n + 1 2 j = 1 n ( a n j a n j + 1 ) F j 1 2 ( a n b n ) F 1 2 b ^ n F 0 ] , χ h ) + ( F n + 1 2 + F n + 1 2 , ψ h ) ( H n + 1 2 , χ h ) = ( P Δ t F n + 1 2 , χ h ) ( Δ t α Γ ( 2 α ) [ a 0 F n + 1 2 j = 1 n ( a n j a n j + 1 ) F j 1 2 ( a n b n ) F 1 2 b ^ n F 0 ] , χ h ) + ( R 3 n + 1 2 , χ h ) .
In (39), we set φ h = E n + 1 2 , χ h = F n + 1 2 , and ψ h = H n + 1 2 , and add the resulting equations to obtain
( P Δ t E n + 1 2 , E n + 1 2 ) + ( P Δ t F n + 1 2 , F n + 1 2 ) + ( P Δ t H n + 1 2 , H n + 1 2 ) + ( Δ t α Γ ( 2 α ) [ a 0 F n + 1 2 j = 1 n ( a n j a n j + 1 ) F j 1 2 ( a n b n ) F 1 2 b ^ n F 0 ] , F n + 1 2 ) = ( P Δ t E n + 1 2 , E n + 1 2 ) ( P Δ t F n + 1 2 , F n + 1 2 ) ( P Δ t H n + 1 2 , H n + 1 2 ) + ( F n + 1 2 + F n + 1 2 , E n + 1 2 ) + ( F n + 1 2 + F n + 1 2 , F n + 1 2 ) ( f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n 2 , H n + 1 2 ) ( Δ t α Γ ( 2 α ) [ a 0 F n + 1 2 j = 1 n ( a n j a n j + 1 ) F j 1 2 ( a n b n ) F 1 2 b ^ n F 0 ] , F n + 1 2 ) + ( R 1 n + 1 2 , E n + 1 2 ) + ( R 2 n + 1 2 , H n + 1 2 ) + ( R 3 n + 1 2 , F n + 1 2 ) .
Now, we need to estimate all terms on the right-hand side of (40). Using Cauchy–Schwarz inequality, we have
( P Δ t E n + 1 2 , E n + 1 2 ) ( P Δ t F n + 1 2 , F n + 1 2 ) ( P Δ t H n + 1 2 , H n + 1 2 ) + ( F n + 1 2 + F n + 1 2 , E n + 1 2 ) + ( F n + 1 2 + F n + 1 2 , F n + 1 2 ) C ( P Δ t E n + 1 2 2 + P Δ t F n + 1 2 2 + P Δ t H n + 1 2 2 + F n + 1 2 2 ) + C ( E n + 1 2 2 + F n + 1 2 2 + H n + 1 2 2 ) .
Applying the mean value theorem and Cauchy–Schwarz inequality, we have
( f u ( u n + 1 ) v n + 1 + f u ( u n ) v n 2 f u ( u h n + 1 ) v h n + 1 + f u ( u h n ) v h n 2 , H n + 1 2 ) = 1 2 ( f u ( u n + 1 ) ( v n + 1 v h n + 1 ) + ( f u ( u n + 1 ) f u ( u h n + 1 ) ) v h n + 1 + f u ( u n ) ( v n v h n ) + ( f u ( u n ) f u ( u h n ) ) v h n , H n + 1 2 ) 1 2 ( f u ( u n + 1 ) v n + 1 v h n + 1 + f u u ( θ ¯ n + 1 ) u n + 1 u h n + 1 v h n + 1 + f u ( u n ) v n v h n + f u u ( θ ¯ n ) u n u h n v h n ) H n + 1 2 . C ( E n + 1 2 + F n + 1 2 + E n 2 + F n 2 + E n + 1 2 + F n + 1 2 + H n + 1 2 + E n 2 + F n 2 + H n 2 ) ,
where we use the boundedness of f u ( u n ) and the following bounded inequality:
f u u ( θ ¯ n ) + v h n C ,
where one can apply inverse inequality [25], and use a similar method as the one in [7,26].
Making use of (9), (3), Cauchy–Schwarz inequality, as well as Young inequality, we have
( Δ t α Γ ( 2 α ) [ a 0 F n + 1 2 j = 1 n ( a n j a n j + 1 ) F j 1 2 ( a n b n ) F 1 2 b ^ n F 0 ] , F n + 1 2 ) + ( R 1 n + 1 2 , E n + 1 2 ) + ( R 2 n + 1 2 , H n + 1 2 ) + ( R 3 n + 1 2 , F n + 1 2 ) = ( 1 Γ ( 1 α ) 0 t n + 1 2 F s d s ( t n + 1 2 s ) α + μ F 0 Γ ( 1 α ) t n + 1 2 α + O ( Δ t 2 α ) , F n + 1 2 ) + ( R 1 n + 1 2 , E n + 1 2 ) + ( R 2 n + 1 2 , H n + 1 2 ) + ( R 3 n + 1 2 , F n + 1 2 ) C [ ( 1 + μ t n + 1 2 α ) h 2 k + 2 + Δ t 4 2 α + E n + 1 2 2 + F n + 1 2 2 + H n + 1 2 2 ] .
Making a combination for (41)–(44) and using (18), we have
( E n + 1 2 + F n + 1 2 + H n + 1 2 ) ( E n 2 + F n 2 + H n 2 ) 2 Δ t + Δ t α 2 Γ ( 2 α ) j = 1 n + 1 a n j + 1 F j 1 2 2 = Δ t α 2 Γ ( 2 α ) j = 1 n a n j F j 1 2 2 + Δ t α 2 Γ ( 2 α ) b ^ n F 0 2 + C [ ( 1 + μ t n + 1 2 α ) h 2 k + 2 + Δ t 4 2 α + P Δ t E n + 1 2 2 + P Δ t F n + 1 2 2 + P Δ t H n + 1 2 2 + E n + 1 2 + F n + 1 2 + E n 2 + F n 2 + E n + 1 2 + F n + 1 2 + H n + 1 2 + E n 2 + F n 2 + H n 2 ] .
With given conditions E 0 = 0 , F 0 = 0 , H 0 = 0 , we use (23) to arrive at
Ξ ( F n + 1 , H n + 1 ) + E n + 1 2 Ξ ( F n , H n ) + E n 2 + C Δ t [ ( 1 + μ t n + 1 2 α ) h 2 k + 2 + Δ t 4 2 α + h 2 m + 2 + E n + 1 2 + F n + 1 2 + H n + 1 2 + E n 2 + F n 2 + H n 2 ] .
Sum for (46) with respect to n to arrive at
Ξ ( F n + 1 , H n + 1 ) + E n + 1 2 Ξ ( F 1 , H 1 ) + E 1 2 + C Δ t j = 1 n [ ( 1 + μ t j + 1 2 α ) h 2 k + 2 + Δ t 4 2 α + h 2 m + 2 ] + C Δ t j = 1 n + 1 [ E j 2 + F j 2 + H j 2 ] .
For n = 0 , we use a similar derivation to the one of n 1 and apply triangle inequality to arrive at
Ξ ( F 1 , H 1 ) + E 1 2 C [ ( 1 + μ t 1 2 α ) h 2 k + 2 + Δ t 4 2 α + h 2 m + 2 ] .
Substitute (48) into (47) and use the Gronwall lemma to obtain
Ξ ( F n + 1 , H n + 1 ) + E n + 1 2 C [ ( 1 + μ t n + 1 2 α ) h 2 k + 2 + Δ t 4 2 α + h 2 m + 2 ] , n 0 .
Combining (49), (34), and (36) with (38) and noting that α = β 1 , we complete the proof of the theorem. □
Remark 2.
Compared with the classical mixed element method for fourth-order partial differential equations, our method can approximate simultaneously three variables with optimal error estimates in L 2 -norm. More importantly, we can obtain directly optimal error estimates in L 2 -norm for auxiliary variables in solving fourth-order PDEs, which are difficult to achieve by using classical mixed element methods [6,7,8].

4. Numerical Tests

Here, we will verify the theoretical results by numerical computing. In (1), we take space domain Ω ¯ = [ 0 , 1 ] 2 , time interval J ¯ = [ 0 , 1 ] , nonlinear term f ( u ) = u 3 u , initial conditions with u ( x , y , 0 ) = 0 , u 1 ( x , y ) = 0 , and exact solution u = t 3 sin ( 2 π x ) sin ( 2 π y ) ; we can obtain the source term g ( x , y , t ) and two auxiliary variables v = 3 t 2 sin ( 2 π x ) sin ( 2 π y ) and σ = t 3 sin ( 2 π x ) sin ( 2 π y ) ( 8 π 2 + t 6 sin ( 2 π x ) 2 sin ( 2 π y ) 2 1 ) . In the following numerical calculations, the order of convergence in space is calculated by the following formula with a sufficiently small time step size Δ t
Order = log h 1 h 2 ϕ ϕ h 1 ϕ ϕ h 2 ,
where h k ( k = 1 , 2 ) represents different space mesh step lengths.
For implementing the new mixed element algorithm, we approximate the spatial direction by the finite element method with the basis function P ( x , y ) = a + b x + c y + d x y and discretize the time direction by using the modified L 1 Crank–Nicolson scheme. In Table 1, by taking the fixed time mesh parameter Δ t = 1 / 200 , changed spatial step length sizes h = 2 / 9 , 2 / 16 and 2 / 25 , and different parameters β = 1.1 , 1.5 , 1.9 , we show the L 2 -norm error estimates and second-order convergence data in space. In Table 2 and Table 3, we compute the convergence results v and σ , respectively. From Table 1, Table 2 and Table 3, one can see that the numerical method is effective for solving nonlinear fourth-order fractional diffusion-wave equation models with a smooth solution.

5. Concluding Remarks

From the calculated results in Table 1, Table 2 and Table 3, one can see that our method for solving fourth-order fractional diffusion-wave equations in this article can obtain optimal error estimates in L 2 -norm for three variables, which is in agreement with the derived theoretical results. These results for auxiliary variables are difficult to achieve directly by using classical mixed element methods [6,7,8].
In the future, we will improve our mixed element method by combining other techniques [7,27,28] with high-order time approximate schemes and develop their optimal numerical theories.

Author Contributions

Conceptualization, J.W.; methodology, J.W., Y.L. and H.L.; software, B.Y.; validation, Y.L., H.L. and Z.F.; formal analysis, J.W.; writing—original draft preparation, J.W.; writing—review and editing, Y.L., H.L. and Z.F.; funding acquisition, J.W., Y.L., H.L. and Z.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (12061053, 12161063), Natural Science Foundation of Inner Mongolia (2020MS01003, 2021MS01018), Young Innovative Talents Project of Grassland Talents Project, Scientific Research Projection of Higher Schools of Inner Mongolia (NJZY21266), and Program for Innovative Research Team in Universities of Inner Mongolia Autonomous Region (NMGIRT2207).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are computed by using our mixed element method.

Conflicts of Interest

The authors declare no conflict of interest.The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
MDPIMultidisciplinary Digital Publishing Institute
DOAJDirectory of Open Access Journals
TLAThree-letter acronym
LDLinear dichroism

References

  1. Nikan, O.; Tenreiro Machado, J.A.; Golbabai, A. Numerical solution of time-fractional fourth-order reaction-diffusion model arising in composite environments. Appl. Math. Model. 2021, 89, 819–836. [Google Scholar] [CrossRef]
  2. Hu, X.L.; Zhang, L.M. On finite difference methods for fourth-order fractional diffusion-wave and subdiffusion systems. Appl. Math. Comput. 2012, 218, 5019–5034. [Google Scholar] [CrossRef]
  3. Liu, Y.; Fang, Z.C.; Li, H.; He, S. A mixed finite element method for a time-fractional fourth-order partial differential equation. Appl. Math. Comput. 2014, 243, 703–717. [Google Scholar] [CrossRef]
  4. Agrawal, O.P. A general solution for a fourth-order fractional diffusion-wave equation defined in a bounded domain. Comput. Struct. 2001, 79, 1497–1501. [Google Scholar] [CrossRef]
  5. Jafari, H.; Dehghan, M.; Sayevand, K. Solving a fourth-order fractional diffusion-wave equation in a bounded domain by decomposition method. Numer. Meth. Part. Differ. Equ. 2008, 24, 1115–1126. [Google Scholar] [CrossRef]
  6. Liu, N.; Liu, Y.; Li, H.; Wang, J.F. Time second-order finite difference/finite element algorithm for nonlinear time-fractional diffusion problem with fourth-order derivative term. Comput. Math. Appl. 2018, 75, 3521–3536. [Google Scholar]
  7. Liu, Y.; Du, Y.W.; Li, H.; Li, J.C.; He, S. A two-grid mixed finite element method for a nonlinear fourth-order reaction-diffusion problem with time-fractional derivative. Comput. Math. Appl. 2015, 70, 2474–2492. [Google Scholar] [CrossRef]
  8. Liu, Y.; Du, Y.W.; Li, H.; He, S.; Gao, W. Finite difference/finite element method for a nonlinear time-fractional fourth-order reaction-diffusion problem. Comput. Math. Appl. 2015, 70, 573–591. [Google Scholar] [CrossRef]
  9. Ji, C.C.; Sun, Z.Z.; Hao, Z.P. Numerical algorithms with high spatial accuracy for the fourth-order fractional sub-diffusion equations with the first Dirichlet boundary conditions. J. Sci. Comput. 2016, 66, 1148–1174. [Google Scholar] [CrossRef]
  10. Ran, M.H.; Zhang, C.J. New compact difference scheme for solving the fourth-order time fractional sub-diffusion equation of the distributed order. Appl. Numer. Math. 2018, 129, 58–70. [Google Scholar] [CrossRef]
  11. Nandal, S.; Pandey, D.N. Numerical solution of non-linear fourth order fractional sub-diffusion wave equation with time delay. Appl. Math. Comput. 2020, 369, 124900. [Google Scholar] [CrossRef]
  12. Sun, H.; Sun, Z.Z.; Du, R. A linearized second-order difference scheme for the nonlinear time-fractional fourth-order reaction-diffusion equation. Numer. Math. Theory Meth. Appl. 2019, 12, 1168–1190. [Google Scholar]
  13. Huang, J.F.; Qiao, Z.; Zhang, J.N.; Arshad, S.; Tang, Y.F. Two linearized schemes for time fractional nonlinear wave equations with fourth-order derivative. J. Appl. Math. Comput. 2021, 66, 561–579. [Google Scholar] [CrossRef]
  14. Abbaszadeh, M.; Dehghan, M. Direct meshless local Petrov–Galerkin (DMLPG) method for time-fractional fourth-order reaction–diffusion problem on complex domains. Comput. Math. Appl. 2020, 79, 876–888. [Google Scholar] [CrossRef]
  15. Yang, X.H.; Zhang, H.X.; Xu, D. Orthogonal spline collocation method for the fourth-order diffusion system. Comput. Math. Appl. 2018, 75, 3172–3185. [Google Scholar] [CrossRef]
  16. Zhang, H.X.; Yang, X.H.; Xu, D. A high-order numerical method for solving the 2D fourth-order reaction-diffusion equation. Numer. Algor. 2019, 80, 849–877. [Google Scholar] [CrossRef]
  17. Tariq, H.; Akram, G. Quintic spline technique for time fractional fourth-order partial differential equation. Numer. Meth. Part Differ. Equ. 2017, 33, 445–466. [Google Scholar] [CrossRef]
  18. Guo, L.; Wang, Z.B.; Vong, S.K. Fully discrete local discontinuous Galerkin methods for some time-fractional fourth-order problems. Int. J. Comput. Math. 2016, 93, 1665–1682. [Google Scholar] [CrossRef]
  19. Du, Y.W.; Liu, Y.; Li, H.; Fang, Z.C.; He, S. Local discontinuous Galerkin method for a nonlinear time-fractional fourth-order partial differential equation. J. Comput. Phys. 2017, 344, 108–126. [Google Scholar] [CrossRef]
  20. Li, X.H.; Wong, P.J.Y. An efficient numerical treatment of fourth-order fractional diffusion-wave problems. Numer. Meth. Part. Differ. Equ. 2018, 34, 1324–1347. [Google Scholar] [CrossRef]
  21. Grasselli, M.; Pierre, M. A splitting method for the Cahn-Hilliard equation with inertial term. Math. Model. Meth. Appl. Sci. 2010, 20, 1363–1390. [Google Scholar] [CrossRef] [Green Version]
  22. Zeng, F.H.; Li, C.P. A new Crank-Nicolson finite element method for the time-fractional subdiffusion equation. Appl. Numer. Math. 2017, 121, 82–95. [Google Scholar] [CrossRef]
  23. Sun, Z.Z.; Wu, X.N. A fully discrete difference scheme for a diffusion-wave system. Appl. Numer. Math. 2006, 56, 193–209. [Google Scholar] [CrossRef]
  24. Lin, Y.M.; Xu, C.J. Finite difference/spectral approximations for the time-fractional diffusion equation. J. Comput. Phys. 2007, 225, 1533–1552. [Google Scholar] [CrossRef]
  25. Ciarlet, P.G. The Finite Element Method for Elliptic Problems; Elsevier: Amsterdam, The Netherlands, 1978. [Google Scholar]
  26. Li, D.F.; Zhang, J.W.; Zhang, Z.M. Unconditionally optimal error estimates of a linearized galerkin method for nonlinear time fractional reaction-subdiffusion equations. J. Sci. Comput. 2018, 76, 848–866. [Google Scholar] [CrossRef]
  27. Bu, W.P.; Liu, X.T.; Tang, Y.F.; Yang, J.Y. Finite element multigrid method for multi-term time fractional advection diffusion equations. Int. J. Model. Simulat. Sci. Comput. 2015, 6, 1540001. [Google Scholar] [CrossRef]
  28. Yue, X.Q.; Shu, S.; Xu, X.; Bu, W.; Pan, K. Parallel-in-time multigrid for space–time finite element approximations of two-dimensional space-fractional diffusion equations. Comput. Math. Appl. 2019, 78, 3471–3484. [Google Scholar] [CrossRef]
Table 1. The convergence results for u with Δ t = 1 / 200 .
Table 1. The convergence results for u with Δ t = 1 / 200 .
β h u u h OrderCPU-Time (s)
1.1 2 / 9 4.1181  × 10 2 1.13
2 / 16 1.3341  × 10 2 1.95904.04
2 / 25 5.5001  × 10 3 1.985519.03
1.5 2 / 9 4.1175  × 10 2 1.10
2 / 16 1.3339  × 10 2 1.95904.02
2 / 25 5.4991  × 10 3 1.985519.40
1.9 2 / 9 4.1169  × 10 2 1.14
2 / 16 1.3336  × 10 2 1.95914.07
2 / 25 5.4977  × 10 3 1.985719.33
Table 2. The convergence results for v with Δ t = 1 / 200 .
Table 2. The convergence results for v with Δ t = 1 / 200 .
β h v v h OrderCPU-Time (s)
1.1 2 / 9 1.2293  × 10 1 1.13
2 / 16 3.9815  × 10 2 1.95944.04
2 / 25 1.6417  × 10 2 1.985119.03
1.5 2 / 9 1.2292  × 10 1 1.10
2 / 16 3.9816  × 10 2 1.95934.02
2 / 25 1.6417  × 10 2 1.985219.40
1.9 2 / 9 1.2293  × 10 1 1.14
2 / 16 3.9819  × 10 2 1.95924.07
2 / 25 1.6418  × 10 2 1.985219.33
Table 3. The convergence results for σ with Δ t = 1 / 200 .
Table 3. The convergence results for σ with Δ t = 1 / 200 .
β h σ σ h OrderCPU-Time (s)
1.1 2 / 9 1.8858  × 10 + 0 1.13
2 / 16 5.9584  × 10 1 2.00244.04
2 / 25 2.4398  × 10 1 2.000719.03
1.5 2 / 9 1.8853  × 10 + 0 1.10
2 / 16 5.9568  × 10 1 2.00254.02
2 / 25 2.4391  × 10 1 2.000719.40
1.9 2 / 9 1.8849  × 10 + 0 1.14
2 / 16 5.9550  × 10 1 2.00264.07
2 / 25 2.4381  × 10 1 2.001019.33
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, J.; Yin, B.; Liu, Y.; Li, H.; Fang, Z. A Mixed Element Algorithm Based on the Modified L1 Crank–Nicolson Scheme for a Nonlinear Fourth-Order Fractional Diffusion-Wave Model. Fractal Fract. 2021, 5, 274. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5040274

AMA Style

Wang J, Yin B, Liu Y, Li H, Fang Z. A Mixed Element Algorithm Based on the Modified L1 Crank–Nicolson Scheme for a Nonlinear Fourth-Order Fractional Diffusion-Wave Model. Fractal and Fractional. 2021; 5(4):274. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5040274

Chicago/Turabian Style

Wang, Jinfeng, Baoli Yin, Yang Liu, Hong Li, and Zhichao Fang. 2021. "A Mixed Element Algorithm Based on the Modified L1 Crank–Nicolson Scheme for a Nonlinear Fourth-Order Fractional Diffusion-Wave Model" Fractal and Fractional 5, no. 4: 274. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5040274

Article Metrics

Back to TopTop