Next Article in Journal
A Family of Derivative Free Algorithms for Multiple-Roots of Van Der Waals Problem
Next Article in Special Issue
Symmetric Spectral Collocation Method for a Kind of Nonlinear Volterra Integral Equation
Previous Article in Journal
Symmetries in the Superposition Model of Extensive Air Shower Development
Previous Article in Special Issue
A Quadruple Integral Involving Product of the Struve Hv(βt) and Parabolic Cylinder Du(αx) Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Estimation for a Type of Fractional Diffusion Equation Based on Compact Difference Scheme

1
School of Statistics and Mathematics, Zhongnan University of Economics and Law, Wuhan 430073, China
2
School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(3), 560; https://doi.org/10.3390/sym14030560
Submission received: 10 February 2022 / Revised: 4 March 2022 / Accepted: 7 March 2022 / Published: 11 March 2022

Abstract

:
Numerical solution and parameter estimation for a type of fractional diffusion equation are considered. Firstly, the symmetrical compact difference scheme is applied to solve the forward problem of the fractional diffusion equation. The stability and convergence of the symmetrical difference scheme are presented. Then, the Bayesian method is considered to estimate the unknown fractional order α of the fractional diffusion equation model. To validate the efficiency of the symmetrical numerical scheme and the estimation method, some simulation tests are considered. The simulation results demonstrate the accuracy of the compact difference scheme and show that the proposed estimation algorithm can provide effective statistical characteristics of the parameter.

1. Introduction

In recent years, the studies of fractional calculus have attracted many researchers [1,2,3], and fractional differential equations have many applications in economics [4], finance [5], physics [6,7,8], bioecology [9,10], porous media dynamics [11,12,13], and so on. However, only few exact solutions can be founded for some linear fractional differential equations. Many researchers turned to develop efficient numerical schemes. For example, the authors in [14,15,16] considered the difference schemes. Li et al. [17,18] proposed a novel algorithm based on linear and quadratic approximations. Zhou et al. [19] considered a nonuniform Alikhanov linearized Galerkin finite element methods. The authors in [20,21] considered the Chebyshev collocation method.
Recently, the research about the inverse problem has attracted more and more interest, one can refer to [22,23,24] and the references therein. Most of the literature on the inverse problem of the fractional differential equations exploited deterministic techniques, such as exact matching, least squares optimizations, without considering measurement error and numerical error. In general, very roughly speaking, one may split the corresponding approaches of recovering the order into the following two categories: solving the corresponding inverse problems analytically or using numerical–analytical methods [25,26,27,28,29,30,31], or using some more soft, metaheuristic, or statistical optimization and regularization techniques [32,33,34,35,36,37]. However, the observations are usually contaminated with measurement error, and the forward problem of the models will bring the numerical error, so the uncertainties are non-ignorable. Some studies for the inverse problem under uncertainties have been provided, such as the extended maximum likelihood method [38], sensitivity analysis method [39], and the Bayesian method [40,41].
The work of this paper was mainly motivated by the study in [41]. However, in the paper, Fan et al. considered homogeneous fractional diffusion equation, and the center Box difference method was used to obtain the numerical solution of the equation. An inhomogeneous fractional diffusion equation was considered in this paper, and a high-order compact difference scheme instead of the second-order center Box difference method was used to obtain the corresponding numerical solution. We considered the Bayesian method to estimate the parameter of the fractional diffusion equation model based on a compact difference scheme. Firstly, a high-order compact difference scheme is employed to solve the forward problem of the fractional diffusion equation. The convergence and stability of the difference scheme are proved, where the difference scheme has 2 α order in the temporal direction and 4 order in the spatial direction. Then, the Bayesian method is considered to estimate the unknown fractional order α of the model, where some additional observations contaminated with measurement errors must be provided. To validate the efficiency of the numerical scheme and the estimation method, some simulation tests are considered. The simulation results demonstrate the accuracy of the compact difference scheme and show that the proposed estimation algorithm can provide effective statistical characteristics of the parameters.
The outline of this paper is as follows. In Section 2, mathematical model of the problem is introduced, and a compact difference scheme is constructed to solve the forward problem. Section 3 introduces the Bayesian method for the inverse problem. In Section 4, numerical tests are presented to validate the performance of the estimation algorithm. Section 5 gives some discussion and remarks.

2. Mathematical Model of the Problem and the Compact Difference Scheme

This section described the mathematical model of the problem, and a compact difference scheme was constructed to solve the forward problem.

2.1. Mathematical Model of the Problem

Consider the following fractional parabolic equation,
α u t α = Q 2 u ( x , t ) x 2 + f ( x , t ) , ( x , t ) ( 0 , 1 ) × ( 0 , 1 ] ,
u ( x , 0 ) = ϕ ( x ) , x [ 0 , 1 ] ,
u ( 0 , t ) = ψ 1 ( t ) , u ( 1 , t ) = ψ 2 ( t ) , t [ 0 , 1 ] ,
where 0 < α < 1 and Q > 0 are constants, ϕ , ψ 1 , and ψ 2 are given functions. The time fractional partial derivative α u t α is defined in the Caputo sense as the following
α u t α = 1 Γ ( 1 α ) 0 t ( t ξ ) α u ( x , ξ ) ξ d ξ ,
where Γ ( ) is the gamma function.
If α is a known parameter, Equations (1)–(3) represent a forward problem; extensive studies have been provided. However, in some practical situations, we cannot know the exact values of α . Under this situation, additional data must be provided to estimate the unknown parameter, which can be called the inverse problem. The additional data may be observed at some measured time for a fixed point x f [ 0 , 1 ] . For simplicity, let x f = 1 in this paper, and we have the following observations,
Y ( t ) = u ( 1 , t ) , 0 < t 1 .
Let F be an operator, F : θ u ( 1 , t ; θ ) , where θ = α ; then, the inverse problem can be seen as the solution of the following equation,
F ( θ ) = Y .
In fact, the observations may be contaminated with measurement errors inevitably. We will consider the Bayesian approach to deal with the inverse problem.

2.2. The Compact Difference Scheme for the Forward Problem

Assuming θ in the operator F to be known in this subsection, we construct the high-order compact difference scheme to solve (6).
Let x i = i h , i = 0 , 1 , , M , and t k = k τ , k = 0 , 1 , , N , where h = 1 / N , τ = 1 / M . Introducing Ω h τ = Ω h × Ω τ , where Ω h = { x i | 0 i M } , Ω τ = { t k | 0 k N } . Denote U i k = u ( x i , t k ) , 0 i M , 0 k N . Let
V = { v i k | 0 i M , 0 k N , }
be the grid function space defined on Ω h τ . Introducing the following notations
δ x v i + 1 / 2 k = v i + 1 k v i k h , δ x 2 v i k = v i + 1 k 2 v i k + v i 1 k h 2 , H v i k = 1 12 ( v i 1 k + 10 v i k + v i + 1 k ) ,
δ t α v k = τ α Γ ( 2 α ) [ d 0 v k j = 1 k 1 ( d k j 1 d k j ) v j d k 1 v 0 ] ,
where d j = ( j + 1 ) 1 α j 1 α , j 0 .
One can use Lemma 1 to deal with (4).
Lemma 1
([42]). Suppose 0 < α < 1 , y C 2 [ 0 , t k ] , it holds that
| 1 Γ ( 1 α ) 0 t k y ( s ) d s ( t k s ) α τ α Γ ( 2 α ) [ d 0 y ( t k ) j = 1 k 1 ( d k j 1 d k j ) y ( t j ) d k 1 y ( t 0 ) ] | 1 Γ ( 2 α ) [ 1 α 12 + 2 2 α 2 α ( 1 + 2 α ) ] max 0 t t k | y ( t ) | τ 2 α .
And d j satisfies Lemma 2.
Lemma 2
([43]). Assume 0 < α < 1 , then it holds that
(1) 
d j decreases monotonically as j increases, and 0 < d j 1 ;
(2) 
d 0 = 1 , j = 1 k 1 ( d k j 1 d k j ) = d 0 d k 1 .
The following Lemma 3 will play the key role to construct the high-order compact difference scheme.
Lemma 3
([44]). Suppose y ( x ) C 6 [ x i 1 , x i + 1 ] , then we have
1 12 [ y ( x i 1 ) + 10 y ( x i ) + y ( x i + 1 ) ] 1 h 2 [ y ( x i 1 ) 2 y ( x i ) + y ( x i + 1 ) ] = h 4 240 y ( 6 ) ( ι i ) ,
where ι i ( x i 1 , x i + 1 ) .
Considering (1) at the point ( x i , t k ) , we have
α u t α ( x i , t k ) = Q 2 u x 2 ( x i , t k ) + f ( x i , t k ) , 0 i M , 0 k N .
By Lemma 1, we have
α u t α ( x i , t k ) = δ t α U i k + O ( τ 2 α ) .
Substituting (8) into (7), then acting H on (7), we have
H δ t α U i k = Q ( H 2 u x 2 ( x i , t k ) ) + H f ( x i , t k ) + O ( τ 2 α ) .
From Lemma 3 and Taylor expansion, we have
H 2 u x 2 ( x i , t k ) = δ x 2 U i k + O ( h 4 ) .
Substituting (10) into (9), we have
H δ t α U i k = Q δ x 2 U i k + H f ( x i , t k ) + O ( τ 2 α + h 4 ) , 1 i M 1 , 1 k N .
Discretizing (2) and (3), we obtain
U i 0 = ϕ ( x i ) , 0 i M ,
U 0 k = ψ 1 ( t k ) , U M k = ψ 2 ( t k ) , 1 k N .
Replacing U i k by u i k in (11)–(13), the following compact difference scheme solving for (1)–(3) can be obtained,
H δ t α u i k = Q δ x 2 u i k + H f ( x i , t k ) , 1 i M 1 , 1 k N ,
u i 0 = ϕ ( x i ) , 0 i M ,
u 0 k = ψ 1 ( t k ) , u M k = ψ 2 ( t k ) , 1 k N .
There are some results on schemes (14)–(16).
Theorem 1.
The difference scheme in Equations (14)–(16) has a unique solution.
The proof of Theorem 1 can be referred to the corresponding proof in [45].
Denote e i k = U i k u i k , 0 i M , 0 k N and e = max 1 i M | e i | ; then, we can obtain the following error equations
H δ t α e i k = Q δ x 2 e i k + R i k , 1 i M 1 , 1 k N ,
e i 0 = 0 , 0 i M ,
e 0 k = 0 , e M k = 0 , 1 k N ,
where R i k = O ( τ 2 α + h 4 ) .
To prove the following theorems, introducing the following notations
( u , v ) = h i = 1 M 1 u i v i , u = ( u , u ) , u = max 1 i M 1 | u i | , δ x u , δ x v = h i = 0 M 1 ( δ x u i + 1 2 ) ( δ x v i + 1 2 ) , | δ x u | 1 = δ x u , δ x u .
Theorem 2.
Assume u ( x , t ) C 6 , 2 ( [ 0 , 1 ] × ( 0 , 1 ] ) is the solution of (1)–(3), and { u i k | 0 i M , 0 k N } is the solution of (14)–(16). Then, we have
e k C ( τ 2 α + h 4 ) , 0 k N ,
where C is a positive constant independent of h and τ.
Proof. 
Multiplying (17) with h δ t α e i k and summing up for i from 1 to M-1, one can obtain
h i = 1 M 1 ( H δ t α e i k ) ( δ t α e i k ) = Q h i = 1 M 1 ( δ x 2 e i k ) ( δ t α e i k ) + h i = 1 M 1 ( R i k ) ( δ t α e i k ) , 1 i M 1 , 1 k N .
Then, we will estimate each term of (21). From the discrete Green formula, we have
h i = 1 M 1 ( H δ t α e i k ) ( δ t α e i k ) = δ t α e i k 2 + h 2 12 h i = 1 M 1 ( δ x 2 ( δ t α e i k ) ) ( δ t α e i k ) = δ t α e i k 2 h 2 12 δ x δ t α e i k 2 2 3 δ t α e i k 2 ,
where inequality ( a + b ) 2 2 ( a 2 + b 2 ) is used.
By the discrete Green formula, one has
Q h i = 1 M 1 ( δ x 2 e i k ) ( δ t α e i k ) = Q h i = 0 M 1 ( δ x e i + 1 / 2 k ) ( δ t α δ x e i + 1 / 2 k ) = Q λ h i = 0 M 1 ( δ x e i + 1 / 2 k ) [ δ x e i + 1 / 2 k j = 1 k 1 ( d k j 1 d k j ) δ x e i + 1 / 2 j d k 1 δ x e i + 1 / 2 0 ] = Q λ { | δ x e k | 1 2 j = 1 k 1 ( d k j 1 d k j ) δ x e k , δ x e j d k 1 δ x e k , δ x e 0 } Q λ { | δ x e k | 1 2 j = 1 k 1 ( d k j 1 d k j ) | δ x e k | 1 2 + | δ x e j | 1 2 2 d k 1 | δ x e k | 1 2 + | δ x e 0 | 1 2 2 } = Q λ { | δ x e k | 1 2 2 j = 1 k 1 ( d k j 1 d k j ) | δ x e j | 1 2 2 d k 1 | δ x e 0 | 1 2 2 } .
From Cauchy–Schwarz inequality, we have
h i = 1 M 1 ( R i k ) ( δ t α e i k ) h i = 1 M 1 ( ( R i k ) 2 2 ε + ε 2 ( δ t α e i k ) 2 ) = 1 2 ε R k 2 + ε 2 δ t α e k 2 .
Inserting Equations (22)–(24) into (21), and take ε = 4 / 3 , we obtain
| δ x e k | 1 2 j = 1 k 1 ( d k j 1 d k j ) | δ x e j | 1 2 + d k 1 | δ x e 0 | 1 2 + 3 λ 4 Q R k 2 .
Noticing R i k = O ( τ 2 α + h 4 ) , from Lemma 6 in [45], we have
| δ x e k | 1 2 C 1 ( τ 2 α + h 4 ) 2 ,
where C 1 is a positive constant independent of h and τ . Then, we have
e k C 1 ( τ 2 α + h 4 ) 2 / 2 C ( τ 2 α + h 4 ) .
The proof is completed. □
Similar to the proof of Theorem 2, we have the following result about the stability of the difference scheme in Equations (14)–(16).
Theorem 3.
Under the same condition of Theorem 2, the difference scheme in Equations (14)–(16) is stable.

3. Parameter Estimation for the Inverse Problem

Considering the measurement error and numerical error in the observations, then (6) has the following form,
Y = F ( θ ) + ε ,
where ε N ( 0 , σ 2 ) , and σ 2 is assumed to be a known parameter. θ = α is the parameter to be estimated from the observations governed by (26).
Assume α U ( 0 , 1 ) is the prior probability distributions. We have the prior probability density function of θ ,
p ( θ ) = p ( α ) = 1 , 0 < α < 1 , 0 , o t h e r w i s e ,
From (26), we have the following likelihood function,
p ( Y | θ ) = 1 ( 2 π σ 2 ) n / 2 exp ( ( F ( θ ) Y ) ( F ( θ ) Y ) 2 σ 2 ) ,
where Y is the observation vector, n is the length of Y.
By the Bayes theorem, the posterior probability density function is
p ( θ | Y ) exp ( ( F ( θ ) Y ) ( F ( θ ) Y ) 2 σ 2 ) .
The Bayesian method for estimating θ = α is as follows,
(1)
Draw an initial value of θ 0 U ( 0 , 1 ) ,
(2)
for i = 0: N s i m 1
draw u U ( 0 , 1 ) ,
draw θ q ( θ | θ i ) ,
if u < ς , then θ i + 1 = θ ,
else θ i + 1 = θ i ,
where
ς = m i n { 1 , p ( θ | Y ) q ( θ i | θ ) p ( θ i | Y ) q ( θ | θ i ) } ,
and
q ( · | · ) U ( 0 , 1 ) ,
where N s i m is the number of iterations.
In the Bayesian method, we will use the posterior mean as the estimation of θ .
Remark 1.
F ( θ ) in (29) can be obtained through (5), where a higher compact difference scheme is used. However, the traditional Crank–Nicolson difference scheme and other algorithms such as finite element method and finite volume method can also be used to solve (5). The advantage of the compact difference is that it is highly accurate, easy to program, and easily parallel. The advantage of the Bayesian estimation algorithm is that it can deal with contaminated observations and inhomogeneous fractional diffusion equation, and the most disadvantage of the Bayesian estimation algorithm is that it is time consuming.

4. Simulation

This section considers the following example. The configuration of the computer is Intel(R) Core(TM) i5-8265U CPU @ 1.60 GHz 1.80 GHz, Memory 4 GB (Lenovo, Beijing, China). All the codes are compiled in MATLAB 2009(a).
Example 1.
α u t α = Q 2 u ( x , t ) x 2 + f ( x , t ) , ( x , t ) ( 0 , 1 ) × ( 0 , 1 ] , u ( x , 0 ) = 0 , x [ 0 , 1 ] , u ( 0 , t ) = 0 , u ( 1 , t ) = 0 , t ( 0 , 1 ] ,
the exact solution of (30) is u ( x , t ) = e x x 2 ( 1 x ) 2 t 2 + α , and
f ( x , t ) = Γ ( 3 + α ) 2 t 2 e x x 2 ( 1 x ) 2 Q e x ( 2 8 x + x 2 + 6 x 3 + x 4 ) t 2 + α .

4.1. Simulation for the Forward Problem

As for the forward problem, we will validate the performance of the scheme (14)–(16) for solving (30). Introducing the following maximum error,
e ( h , τ ) = max 1 k N U k u k ,
and we define the following convergence order as
R τ = l o g ( e ( h , τ 1 ) / e ( h , τ 2 ) ) l o g ( τ 1 / τ 2 ) , R h = l o g ( e ( h 1 , τ ) / e ( h 2 , τ ) ) l o g ( h 1 / h 2 ) .
For calculating R τ , h needs to be fixed and small enough. While for R h , τ needs to be fixed and small enough.
To compare the errors between the numerical solution and the exact solution of (30), Figure 1 was provided. Taking temporal and spatial step sizes as τ = h = 1 / 10 , 1 / 20 , 1 / 40 , 1 / 80 respectively, let α = 0.2 , and we can see that the value of the error is within [ 0 , 1.6 × 10 5 ] . The numerical solution can approximate the exact solution well enough.
Table 1 shows the maximum errors in the temporal directions for α = 0.2 , 0.6 , and 0.8, respectively, where we fixed the spatial step to be h = 1 / 2000 . We can see that the temporal convergence order is 2 α .
Table 2 shows the maximum errors in the spatial directions for α = 0.2 , 0.6 , and 0.8, respectively, where we fixed the temporal step τ = 1/10,000. In most of the cases, we can conclude that the spatial convergence order is 4. However, we can see that the convergence of the considered numerical scheme is not steady under the condition that α = 0.8 .

4.2. Simulation for the Inverse Problem

As for the inverse problem, we will testify to the validity of the Bayesian method provided. To obtain the observations Y governed by (26), let θ = α = 0.2 in Equation (1) be the given value (note: θ = α = 0.6 and 0.8 will take the same conclusion, only θ = α = 0.2 is considered in this subsection for the sake of brevity). Firstly, we applied the compact difference scheme (14)–(16) to discrete (30), where M = N = 100 . A series of points { u M k } , k = 0 , 1 , , N were generated. Then, we picked up every fourth (two or one) point from the generated points, finally add the noise on each data, and we get the observations with length n = 25 ( n = 50 , n = 100 ).
However, the length of observations n and σ 2 may affect the estimation results of θ . In the following tests, we will evaluate such effects caused by different values of n and σ 2 . In the following tests, the length of the Markov chain is taken to be N s i m = 5000 , and the burn in number is taken to be 1000.
Let n = 100 be fixed, and assume σ [ 1 × 10 8 , 5 × 10 7 ] . For different σ , we will inspect the posterior mean, standard error, and quantile of α . Table 3 was obtained, and we can see that the measurement error with bigger σ may take more uncertainty of the estimation results.
To inspect the effect of different n for α , let σ = 2 × 10 7 be fixed, assume n = 25 , n = 50 , and n = 100 . For different n, we also inspect the posterior mean, standard error, and quantile of α . Table 4 was obtained, when n = 50 and n = 100 , the different values of n make little impact on the estimation results. While when n = 25 , there is more uncertainty in the estimation results of α . We can draw the conclusion that there is little influence for estimating α under bigger values of different n. From the computational time, it can be seen that we need more estimation time when a bigger n is given.
Remark 2.
In the tests of simulation for the inverse problem, we set the seed as ’sum(100*clock)’ in the ’rand’ function to specify the settings of the random number generator. Then, the estimation results can always be the same.

5. Conclusions

This paper considered both the forward problem and inverse problem of a type of fractional diffusion equation. Firstly, a compact difference scheme was provided to solve the forward problem of the fractional diffusion equation. The convergence and stability of the difference scheme were proved, and the convergence order is O ( τ 2 α + h 4 ) . Secondly, as for the inverse problem, a Bayesian algorithm was provided to estimate the fractional order α of the diffusion equation under the condition that some additional observations contaminated with measurement errors were provided. Numerical experiments were provided to testify the convergence results both in the temporal direction and in the spatial direction and the estimation efficiency of the Bayesian algorithm.

Author Contributions

All the authors contributed equally to this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by National Natural Science Foundation of P.R. China (No. 71974204).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bazhlekova, E. Fractional Evolution Equations in Banach Spaces. Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2001. [Google Scholar]
  2. Kostic, M. Abstract Volterra Integro-Differential Equations; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  3. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA; London, UK, 1974. [Google Scholar]
  4. Xie, C.; Xia, X.; Aghdam, Y.E.; Farnam, B.; Jafari, H.; Wang, S. The numerical strategy of tempered fractional derivative in European double barrier option. Fractals 2022, 30, 2240049. [Google Scholar] [CrossRef]
  5. Mesgarani, H.; Aghdam, Y.E.; Tavakoli, H. Numerical Simulation to Solve Two-Dimensional Temporal–Space Fractional Bloch–Torrey Equation Taken of the Spin Magnetic Moment Diffusion. Int. J. Appl. Comput. Math. 2021, 7, 94. [Google Scholar] [CrossRef]
  6. Bouchaud, J.; Georges, A. Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications. Phys. Rep. 1990, 195, 127–293. [Google Scholar] [CrossRef]
  7. Li, D.; Wang, J.; Zhang, J. Unconditionally convergent L1-Galerkin FEMs for nonlinear time-fractional schrödinger equations. SIAM J. Sci. Comput. 2017, 39, A3067–A3088. [Google Scholar] [CrossRef]
  8. Li, L.; Li, D. Exact solutions and numerical study of time fractional Burgers’ equations. Appl. Math. Lett. 2020, 100, 106011. [Google Scholar] [CrossRef]
  9. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  10. Li, D.; Zhang, C. Long time numerical behaviors of fractional pantograph equations. Math. Comput. Simul. 2020, 172, 244–257. [Google Scholar] [CrossRef]
  11. Berardi, M.; Difonzo, F.V. A quadrature-based scheme for numerical solutions to Kirchhoff transformed Richards’ equation. J. Comput. Dyn. 2022. [Google Scholar] [CrossRef]
  12. Liu, F.; Fukumoto, Y.; Zhao, X. A multi level linearized CrankCNicolson scheme for Richards equation under variable flux boundary conditions. Appl. Anal. 2021. [Google Scholar] [CrossRef]
  13. Garrappa, R. Trapezoidal methods for fractional differential equations: Theoretical and computational aspects. Math. Comput. Simul. 2015, 110, 96–112. [Google Scholar] [CrossRef] [Green Version]
  14. Alikhanov, A. A new difference scheme for the time fractional diffusion equation. J. Comput. Phys. 2015, 280, 424–438. [Google Scholar] [CrossRef] [Green Version]
  15. Langlands, T.; Henry, B. The accuracy and stability of an implicit solution method for the fractional diffusion equation. J. Comput. Phys. 2005, 205, 719–736. [Google Scholar] [CrossRef]
  16. Chen, X.; Di, Y.; Duan, J.; Li, D. Linearized compact ADI schemes for nonlinear time-fractional Schrodinger equations. Appl. Math. Lett. 2018, 84, 160–167. [Google Scholar] [CrossRef]
  17. Li, D.; Sun, W.; Wu, C. A novel numerical approach to time-fractional parabolic equations with nonsmooth solutions. Numer. Math. Theory Meth. Appl. 2021, 14, 355–376. [Google Scholar]
  18. Li, D.; Liao, H.; Sun, W.; Wang, J.; Zhang, J. Analysis of L1-Galerkin FEMs for time-fractional nonlinear parabolic problems. Commun. Comput. Phys. 2018, 24, 86–103. [Google Scholar] [CrossRef] [Green Version]
  19. Zhou, B.; Chen, X.; Li, D. Nonuniform Alikhanov linearized Galerkin finite element methods for nonlinear time-fractional parabolic equations. J. Sci. Comput. 2020, 85, 39. [Google Scholar] [CrossRef]
  20. Aghdam, Y.E.; Mesgarani, H.; Moremedi, G.M.; Khoshkhahtinat, M. High-accuracy numerical scheme for solving the space-time fractional advection-diffusion equation with convergence analysis. Alex. Eng. J. 2022, 61, 217–225. [Google Scholar] [CrossRef]
  21. Mesgarani, H.; Beiranvand, A.; Aghdam, Y.E. The impact of the Chebyshev collocation method on solutions of the time-fractional Black-Scholes. Math. Sci. 2021, 15, 137–143. [Google Scholar] [CrossRef]
  22. Janno, J.; Kinash, N. Reconstruction of an order of derivative and a source term in a diffusion equation from final measurements. Inverse Probl. 2018, 34, 025007. [Google Scholar] [CrossRef]
  23. Li, G.; Zhang, D.; Jia, X.; Yamamoto, M. Simultaneous inversion for the space-dependent diffusion coefficient and the order in the time-diffusion equation. Inverse Probl. 2013, 29, 65014. [Google Scholar] [CrossRef]
  24. Miller, L.; Yamamoto, M. Coefficient inverse problem for a diffusion equation. Inverse Probl. 2013, 29, 75013. [Google Scholar] [CrossRef]
  25. Hatano, Y.; Nakagawa, J.; Wang, S.; Yamamoto, M. Determination of order in fractional diffusion equation. J. Math-for-Ind. 2013, 5, 118. [Google Scholar]
  26. Jin, B.; Kian, Y. Recovery of the Order of Derivation for Fractional Diffusion Equations in an Unknown Medium. arXiv 2021, arXiv:2103.15166. [Google Scholar]
  27. Jin, B.; Kian, Y. Recovering multiple fractional orders in time-fractional diffusion in an unknown medium. Proc. R. Math. Phys. Eng. Sci. 2021, 477, 20210468. [Google Scholar] [CrossRef]
  28. Jin, B.; Zhou, Z. Recovering the potential and order in one-dimensional time-fractional diffusion with unknown initial condition and source. Inverse Probl. 2021, 37, 105009. [Google Scholar] [CrossRef]
  29. Janno, J. Determination of the order of fractional derivative and a kernel in an inverse problem for a generalized time fractional diffusion equation. Electron. J. Differ. Equ. 2016, 2016, 1–28. [Google Scholar]
  30. Yamamoto, M. Uniqueness in determining fractional orders of derivatives and initial values. Inverse Probl. 2021, 37, 095006. [Google Scholar] [CrossRef]
  31. El Hamidi, A.; Tfayli, A. Identification of the derivative order in fractional differential equations. Math. Methods Appl. Sci. 2021, 44, 8397–8413. [Google Scholar] [CrossRef]
  32. Krasnoschok, M.; Pereverzyev, S.; Siryk, S.V.; Vasylyeva, N. Regularized reconstruction of the order in semilinear subdiffusion with memory. In Springer Proceedings in Mathematics and Statistics; Springer: Singapore, 2020; Volume 310, pp. 205–236. [Google Scholar]
  33. Krasnoschok, M.; Pereverzyev, S.; Siryk, S.V.; Vasylyeva, N. Determination of the fractional order in semilinear subdiffusion equations. Fract. Calc. Appl. Anal. 2020, 23, 694–722. [Google Scholar] [CrossRef]
  34. Gao, F.; Lee, X.J.; Tong, H.Q.; Fei, F.X.; Zhao, H.L. Identification of Unknown Parameters and Orders via Cuckoo Search Oriented Statistically by Differential Evolution for Noncommensurate Fractional-Order Chaotic Systems. Abstr. Appl. Anal. 2013, 2013, 382834. [Google Scholar] [CrossRef] [Green Version]
  35. Yuan, L.G.; Yang, Q.G. Parameter identification and synchronization of fractional-order chaotic systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 305–316. [Google Scholar] [CrossRef]
  36. Bohaienko, V.; Gladky, A. Parameters identification for fractional-fractal model of filtration-consolidation using GPU. In Proceedings of the CMIS 2021: 4th International Workshop on Computer Modeling and Intelligent Systems, Zaporizhzhia, Ukraine, 27 April 2021; pp. 409–418. [Google Scholar]
  37. Bohaienko, V.; Gladky, A.; Romashchenko, M.; Matiash, T. Identification of fractional water transport model with ψ-Caputo derivatives using particle swarm optimization algorithm. Appl. Math. Comput. 2021, 390, 125665. [Google Scholar] [CrossRef]
  38. Fadale, T.D.; Nenarokomov, A.V.; Emery, A.F. Uncertainties in parameter estimation: The inverse problem. Int. J. Heat Mass Transf. 1995, 38, 511–518. [Google Scholar] [CrossRef]
  39. Sun, N.Z.; Yeh, W.G. A stochastic inverse solution for transient groundwater flow: Parameter identification and reliability analysis. Water Resour. Res. 1992, 28, 3269–3280. [Google Scholar] [CrossRef]
  40. Kaipio, J.; Somersalo, E. Statistical and Computational Inverse Problems; Springer: New York, NY, USA, 2005. [Google Scholar]
  41. Fan, W.; Jiang, X.; Chen, S. Parameter estimation for the fractional fractal diffusion model based on its numerical solution. Comput. Math. Appl. 2016, 71, 642–651. [Google Scholar] [CrossRef]
  42. Sun, Z.; Wu, X. A fully discrete difference scheme for a diffusion-wave system. Appl. Numer. Math. 2006, 56, 193–209. [Google Scholar] [CrossRef]
  43. Chen, S.; Liu, F.; Zhuang, P.; Anh, V. Finite difference approximations for the fractional Fokker-Planck equation. Appl. Math. Model. 2009, 33, 256–273. [Google Scholar] [CrossRef] [Green Version]
  44. Sun, Z. The Numerical Methods for Partial Equations; Science Press: Beijing, China, 2005. (In Chinese) [Google Scholar]
  45. Gu, W.; Zhou, Y.; Ge, X. Compact Difference Scheme for Solving Fractional Neutral Parabolic Differential Equation with Proportional Delay. J. Funct. Spaces 2017, 2017, 3679526. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The error curves of | u ( x , 1 ) U ( x , 1 ) | for different temporal and spatial step sizes τ = h = 1 / 10 , 1 / 20 , 1 / 40 , 1 / 80 for solving (30), where α = 0.2 , u ( x , 1 ) , and U ( x , 1 ) denote the numerical solution and exact solution, respectively.
Figure 1. The error curves of | u ( x , 1 ) U ( x , 1 ) | for different temporal and spatial step sizes τ = h = 1 / 10 , 1 / 20 , 1 / 40 , 1 / 80 for solving (30), where α = 0.2 , u ( x , 1 ) , and U ( x , 1 ) denote the numerical solution and exact solution, respectively.
Symmetry 14 00560 g001
Table 1. Maximum errors and convergence order for solving (30) in the temporal direction, where Q = 1.8 , h = 1 / 2000 , and α = 0.2 , 0.6 , and 0.8. * means no value.
Table 1. Maximum errors and convergence order for solving (30) in the temporal direction, where Q = 1.8 , h = 1 / 2000 , and α = 0.2 , 0.6 , and 0.8. * means no value.
τ α = 0.2 α = 0.6 α = 0.8
e ( h , τ ) Rate τ e ( h , τ ) Rate τ e ( h , τ ) Rate τ
1/50 1.121 × 10 6 * 2.388 × 10 5 * 8.765 × 10 5 *
1/100 3.373 × 10 7 1.733 9.139 × 10 6 1.386 3.836 × 10 5 1.192
1/150 1.664 × 10 7 1.742 5.202 × 10 6 1.390 2.363 × 10 5 1.195
1/200 1.007 × 10 7 1.747 3.485 × 10 6 1.392 1.675 × 10 5 1.196
Table 2. Maximum errors and convergence order for solving (30) in the spatial direction, where Q = 1.8 , τ = 1/10,000, and α = 0.2 , 0.6 , and 0.8. * means no value.
Table 2. Maximum errors and convergence order for solving (30) in the spatial direction, where Q = 1.8 , τ = 1/10,000, and α = 0.2 , 0.6 , and 0.8. * means no value.
h α = 0.2 α = 0.6 α = 0.8
e ( h , τ ) Rate h e ( h , τ ) Rate h e ( h , τ ) Rate h
1/10 3.166 × 10 5 * 3.059 × 10 5 * 2.968 × 10 5 *
1/20 1.981 × 10 6 3.998 1.901 × 10 6 4.008 1.718 × 10 6 4.111
1/40 1.238 × 10 7 4.001 1.056 × 10 7 4.171 7.446 × 10 8 4.528
1/80 7.646 × 10 9 4.017 8.946 × 10 9 3.561 1.466 × 10 7 −0.977
Table 3. Impact of different σ on estimation for α of (30), where Q = 1.8 , the true value of α is 0.2 .
Table 3. Impact of different σ on estimation for α of (30), where Q = 1.8 , the true value of α is 0.2 .
σ α mean α std 95 % α quantile Time (s)
5 × 10 7 0.20490.0227[0.1701, 0.2367]819.1719
1 × 10 7 0.19540.0042[0.1884, 0.2014]865.3281
1 × 10 8 0.19980.0003[0.1991, 0.2004]769.8125
Table 4. Impact of different n on estimation for α of (30), where Q = 1.8 , and the true value of α is 0.2 .
Table 4. Impact of different n on estimation for α of (30), where Q = 1.8 , and the true value of α is 0.2 .
n α mean α std 95 % α quantile Time (s)
250.21180.0145[0.1878, 0.2357]827.8750
500.17780.0120[0.1572, 0.1984]836.0469
1000.18730.0082[0.1744, 0.2023]844.9219
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gu, W.; Wei, F.; Li, M. Parameter Estimation for a Type of Fractional Diffusion Equation Based on Compact Difference Scheme. Symmetry 2022, 14, 560. https://doi.org/10.3390/sym14030560

AMA Style

Gu W, Wei F, Li M. Parameter Estimation for a Type of Fractional Diffusion Equation Based on Compact Difference Scheme. Symmetry. 2022; 14(3):560. https://doi.org/10.3390/sym14030560

Chicago/Turabian Style

Gu, Wei, Fang Wei, and Min Li. 2022. "Parameter Estimation for a Type of Fractional Diffusion Equation Based on Compact Difference Scheme" Symmetry 14, no. 3: 560. https://doi.org/10.3390/sym14030560

APA Style

Gu, W., Wei, F., & Li, M. (2022). Parameter Estimation for a Type of Fractional Diffusion Equation Based on Compact Difference Scheme. Symmetry, 14(3), 560. https://doi.org/10.3390/sym14030560

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop