Abstract

This paper deals with an inverse problem of simultaneously determining the space-dependent diffusion coefficient and the fractional order in the variable-order time fractional diffusion equation by the measurements at one interior point. Numerical solution to the forward problem is given by the finite difference scheme, and the homotopy regularization algorithm is applied to solve the inverse problem utilizing Legendre polynomials as the basis functions of the approximate space. The inversion solutions with noisy data which give good approximations to the exact solution demonstrate effectiveness of the inversion algorithm for the simultaneous inversion problem.

1. Introduction

The fractional calculus and fractional differential equations and their applications have attracted much attention during the last two decades. If the usual integer-order derivative in a classical diffusion equation is replaced with a fractional derivative of any real number in time or space, we get a time or space fractional diffusion equation, which has more and more applications in the fields of physics, mechanics, viscoelastic materials, environmental science, hydrology, nanofluid, etc. On the attempting to describe some real processes and complex dynamical systems using the fractional diffusion models, several researches confronted with the situation that the fractional derivative did not remain constant and changed; say, in the interval from 0 to 1, from 1 to 2 or even from 0 to 2. To study such phenomena in mathematics, an effective way is to employ the fractional derivative of the variable order; i.e., the order of the fractional derivative can change with the time or/and space variables. The variable-order fractional models and the corresponding theory of variable-order operators give a new mathematical framework for dealing with complex dynamical systems, which becomes a new direction in the research of fractional differential equations (see, e.g., [1, 2]). For the concept of variable fractional differentials and integrals and the connection between the mathematical concepts of fractional order and physical process, see [35] and references therein.

There are quite a few of researches on numerical methods for solving the variable-order time/space fractional diffusion equation models. We refer to some work given by F. Liu and his group, see [69] for instance, also see the work given by Sun et al. [10], and the master degree thesis given by Zhang [11] for numerical solutions to the variable-order time fractional diffusion models with variable coefficients. Recently, we refer to the work by Bhrawy and Zaky et al. [1216] for numerical methods for multidimensional space/time variable-order fractional differential equations, where approximate solutions of the forward variable-order fractional differential equations and spectral techniques are discussed based on orthogonal polynomials. On concrete applications of the fractional calculus models, we always encounter with inverse problems which need to determine some unknowns in the model using suitable additional data. As for inverse problems arising in the time/space fractional diffusion equations with constant fractional orders, we refer to [1729], and recently see [3032], etc. However, there are few studies on inverse problems associated with variable-order fractional differential equations to our knowledge. Jia and Peng gave a uniqueness result for the inverse source problem in the variable-order space fractional diffusion equation using a Carleman type estimate [33], and Liu et al. [34] gave numerical inversions for determining the variable fractional order in the variable-order time fractional diffusion equation using the homotopy regularization algorithm.

Since the diffusion coefficient denotes the characteristics of the medium, the variable time fractional order is a key index denoting the correlations in the model, and they are always unknown in real-life problems, it is of much importance to determine these two variable-dependent parameters utilizing the inverse problem method. Although the uniqueness and stability of such a simultaneous inverse problem are very difficult to obtain at present, it is still meaningful to study its numerical inversions from numerics. In this paper, we deal with the simultaneous inversion problem of determining the space-dependent fractional order and diffusion coefficient in the variable-order time fractional diffusion equation utilizing the homotopy regularization algorithm. Such an inversion problem is more difficult than those discussed in the previous work by the following reasons:(i)The model is complicated, which involves the variable-dependent fractional order and diffusion coefficient, and the computational complexity is high.(ii)The inverse problem is to determine two kinds of parameters simultaneously, and both of the diffusion coefficient and the fractional order are spatially dependent functions such that the ill-posedness of the numerical inversion is much severe.(iii)The approximate space for the unknowns is different. We employ the Legendre polynomials as the basis functions instead of the ordinary polynomials to perform the inversion algorithm.

Therefore it is much challenging to deal with the simultaneous inversion problem arising in the variable-order fractional diffusion equation. The most contribution of this paper is devoted to give effective numerical inversions for the simultaneous inversion problem utilizing the homotopy regularization in the framework of Legendre polynomials approximations.

The rest of the paper is organized as follows. In Section 2, a finite difference scheme for solving the forward problem is given. In Section 3, the homotopy regularization algorithm is introduced for solving the inverse problem of determining the fractional order and the diffusion coefficient simultaneously. In Section 4, numerical inversions with noisy data are performed in Legendre polynomial approximate space, and concluding remarks are given in Section 5.

2. The Forward Problem and Finite Difference Scheme

2.1. The Forward Problem

Denote and for and . Consider the variable-order time fractional diffusion equation with a space-dependent diffusion coefficient:where is the state variable, is the fractional differential order depending upon the space and time variables and , is the diffusion coefficient, and is a source term. Here denotes the variable-order fractional derivative in the sense of Caputo, which is defined bywhere is the Gamma function, and there exist two positive constants and such that Given the initial value conditionand the homogeneous boundary conditiona forward problem is composed by (1) with the initial boundary conditions (4)-(5).

2.2. The Finite Difference Scheme

In this subsection, an implicit finite difference scheme is set forth to solve the forward problem (1), (4)-(5) numerically. Although a similar difference scheme has been discussed in [11], we still give it here for the completeness of this paper.

By discretizing the space domain by and the time domain by and the variable-order fractional derivative by the difference discretization, we havewhere and is the space mesh step and is the time mesh step.

Next, denote as the midpoint of the neighboring nodes and . The integer-order diffusion term in (1) is discretized by

Denoting , , , and , substituting (6) and (7) into (1), and omitting high-order terms, we getThe initial boundary value conditions are discretized as follows:and

Finally, denote ; (8) is arranged asfor . In the case of , there isand we have for Denote , and letfor . We rewrite (13) in a matrix form:wherefor , and for and .

Lemma 1. The coefficient matrix () defined by (16) is strictly diagonally dominant, and the difference scheme (15) has only one solution.

Proof. For any , there is . Noting there holdswhich means thatThis shows that the matrix () is strictly diagonally dominant, and the difference scheme (15) has unique solution.
Furthermore, it is not difficult to prove that the scheme (15) is of unconditional stability and convergence with a similar method as used in [34].

Lemma 2 (see [11, 34]). The difference scheme (15) is unconditional stable and convergent to the exact solution of the forward problem for any finite time .

2.3. Numerical Testification

In this subsection, we give a numerical computation for the forward problem utilizing the difference scheme (15). The analytic solution is chosen asand the model parameters are taken as follows:and the source term

The numerical results with different time step and space step are listed in Tables 1 and 2 at given , respectively. The relative error in the solutions is expressed bywhere and denote the exact and numerical solutions, respectively. Moreover, the numerical and exact solutions are plotted in Figure 1.

From Tables 1 and 2 and Figure 1, it can be seen that the error in the solutions gradually decreases as the grid step goes to small demonstrating that the numerical scheme is convergent, and the space step has a significant impact on the numerical scheme for given time step.

3. The Inverse Problem and Inversion Algorithm

3.1. The Inverse Problem

For the variable-order fractional diffusion equation (1), suppose that the variable fractional order and the space-dependent diffusion coefficient are both unknown. Thus we encounter with an inverse problem that is to determine the variable diffusion coefficient and the fractional order simultaneously based on the forward problem and some additional measurements on the solution. Let be the fixed measured point and , , be the measured time. The additional data are given asand the vectoris called the measured data vector.

For simplicity we assume that the fractional order function is time independent on concrete inversions; i.e., the simultaneous inversion problem here is to determine the space-dependent functions and based on the forward problem together with the overposed condition (24).

Together with the forward problem and a priori conditions for the diffusion coefficient and the fractional order, we give the following assumptions:(A1) and for ;(A2), and for .

Under the above assumptions for the unknowns we give the inversion algorithm in the following.

3.2. The Homotopy Regularization Algorithm

Let denote the admissible space for the unknowns satisfying conditions (A1) and (A2), andfor some constant . For further assume that they have the generalized Fourier expansionsandrespectively, where is a group of basis functions in the continuous space and and are the expansion coefficients of and , respectively. In numerical approximations, letandbe the approximations to the exact functions and and and be the dimensions of the approximate spaces, respectively. On concrete computations, it is convenient to set a S-dimensional vector corresponding to and a K-dimensional vector corresponding to , respectively. Therefore, to get an approximate solution is equivalent to finding a vector . In what follows, we set and equip with Euclid norm.

For given , the forward problem is solved numerically whose solution is denoted by . Then we get for correspondingly to the measured time and position, which is denoted by a parameter-dependent vectorCombining with the additional conditions (24)-(25), we are to look for an optimal solution to the inverse problem from numerics by minimizing a cost functional between the measured data and the computational output data:here and in the follows denotes the Euclid norm. Furthermore by the homotopy idea and regularization method, we get the following minimization problem:where is the homotopy parameter. The inversion algorithm is proceeded based on the minimization problem (33) by iteration, linearization, differentiation, and homotopy approximations as done in [27, 34, 35].

In fact, by linearization and numerical differentiation approximations, solving (33) is transformed to solve a normal equation combining with an iteration process:where is just the homotopy parameter which taking values from near decreasingly approximating to , is a perturbation vector for any given , denotes the iterative number and is an initial iteration, andwhere is the numerical differential step, () is the unit basis vector in , and , are the computational output and the measured data vectors, respectively. The algorithm can be terminated as long as an optimal perturbation satisfying the condition: , where eps is the given convergent precision. On concrete inversions we employ the homotopy parameter bywhere is the number of iterations, is the preestimated number, and is the adjust parameter.

In the follows, we utilize the above inversion algorithm to perform numerical inversions for the simultaneous inverse problem. We always set the initial iteration , the differential step , the preestimated number , and the adjust parameter in the homotopy parameter given by (36) if there is no specification. Moreover, we choose and in computation of the solution to the forward problem using the finite scheme (15). The step-by-step procedures of the inversion algorithm are given as follows.

Step 1. Give an approximate space composed by orthogonal polynomials on , choose , as the dimensions of the approximate spaces for and , respectively, choose an initial iteration vector (), numerical differentiation step , and convergent precision , and give the additional data vector .

Step 2. Solve the forward problem with the difference scheme (15) to get and the vector and then to obtain the matrix by formula (35).

Step 3. Choose the homotopy parameter via (36) to work out an optimal perturbation by formula (34).

Step 4. If there holds , then the inversion algorithm can be terminated, and is taken as the inversion solution vector. Otherwise, go to Step 2 by replacing with .

4. Numerical Inversion

Consider the forward problem for and We will perform numerical inversions for determining and on simultaneously by the additional data utilizing the homotopy regularization algorithm introduced in Section 3.2.

4.1. Inversion with Accurate Data

Let the exact solution of the forward problem be , and the exact diffusion coefficient and the fractional order be and on , respectively. The additional data are obtained by computing the forward problem using the exact diffusion coefficient and the fractional order, and we utilize Legendre polynomials as the basis functions in the approximate spaces. In the process of numerical implementation, let the approximate diffusion coefficient , and the approximate fractional order . In other words, the exact solution of the simultaneous inversion problem in the finite-dimensional approximate space is expressed via

By utilizing the inversion algorithm with accurate data, the inversion results are obtained which listed in Table 3, where denotes the dimension of the space of the additional data, which is the number of the additional data measured at and (), denotes the relative error in the exact and inversion solutions, and denotes the vector of the inverse solution.

It can be seen from Table 3 that the inversion results give good approximations to the exact solutions even using a few of measurements, and the solutions error goes to small rapidly as the number of the additional data increasing.

4.2. Inversion with Noisy Data

In this subsection, we implement the inversion algorithm using noisy data. Suppose the disturbed additional data vector is expressed bywhere is the noise level and is an random number ranged in .

4.2.1. Inversion for Polynomial Coefficient

Take the model and parameters as utilized in Section 4.1, i.e., the exact solution to the inverse problem is also set to be the vector given by (38), and the inversion is performed using noisy data given by (39). Noting that the computational result for each iteration is different due to the random noise in the data, we repeat the inversion processes ten times with each noise level and take an average of gained inversion results. The average inversion results are listed in Table 4, where denotes the average of the inversion solutions and is the relative error in the solutions. Moreover, the average inversion results and the exact solutions for the diffusion coefficient and the fractional order with different noise levels are plotted in Figures 2(a)2(d) and Figures 3(a)3(d), respectively.

It can be seen from Table 4 and Figures 2 and 3 that the inversion algorithm is stable from numerics for the simultaneous inversion problem of determining the space-dependence diffusion coefficient and variable fractional order; however, the inversion becomes difficult when the noise level is greater than 5% showing that the joint inversion here is of high ill-posedness.

Remark 3. If utilizing the ordinary polynomials as the basic functions of the approximate space to perform the inversion algorithm with noisy data, the average inversion results are listed in Table 5.
It can be seen from Table 5 as compared with Table 4 that the inversion results with ordinary polynomials in the approximate spaces are worse than those of using Legendre orthogonal polynomials.

4.2.2. Inversion for Exponential Coefficient

Also let the exact solution to the forward problem (37) be , but let the exact diffusion coefficient be and the exact fraction order be . Since the exponential function has the truncated error in the polynomials space, we choose to approximate the diffusion coefficient and also to approximate fractional order on the concrete inversions. Therefore, the exact solution vector to the inverse problem here in the approximate space is expressed by

Like done in Section 4.2.1, the inversion is performed utilizing the same parameters with noisy data. The average inversion results are listed in Table 6, and the exact and inversion solutions with different noise levels are plotted in Figures 4(a)4(d) and Figures 5(a)5(d), respectively. From Table 6 and Figures 4 and 5, we can be see that the inversion has a similar result as shown in Section 4.2.1, and the inversion algorithm is testified again to be efficiency for the simultaneous inversion problem.

4.2.3. Inversion for Nonlinear Order Function

In this subsection we perform numerical inversions for the fractional order taking nonlinear functions. Let the exact solution of the forward problem be and the exact diffusion coefficient also be but the fractional order be for . On the finite-dimensional Legendre polynomials approximation we choose to approximate the diffusion coefficient and to approximate the fractional order, and the exact solution vector of the inverse problem is expressed via

Like done in the above, the inversion is performed utilizing the same parameters with noisy data. The average inversion results are listed in Table 7, and the exact and inversion solutions with different noise levels are plotted in Figures 6(a)6(d) and Figures 7(a)7(d), respectively. By Table 7 and Figures 6 and 7, we can see again that the inversion solutions approximate to the exact solutions as the noise level goes to small demonstrating the efficiency of the simultaneous inversion algorithm.

5. Conclusion

The simultaneous inverse problem of determining space-dependent diffusion coefficient and fractional order in the variable-order time fractional diffusion equation is investigated by using measurements at one interior point. Generally speaking, an inverse problem of determining multiple types of coefficients in a PDE is more difficult than that of determining a single type of coefficient, especially for the variable-order fractional differential equations. However, numerical inversions for such kinds of inverse problems can be performed successfully by the homotopy regularization algorithm using Legendre polynomials as the basis functions of the approximate space. The inversion becomes very difficult when the fractional order is time-space-dependent and the noise level is greater than 1%. We will focus our attention on the uniqueness of such inverse problems and theoretical analysis for the inversion algorithm utilizing limited additional data in the near future.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (nos. 11371231 and 11871313).