A Note on the Inverse Problem for a Fractional Parabolic Equation
Abdullah Said Erdogan1and Hulya Uygun1
Academic Editor: Ravshan Ashurov
Received15 May 2012
Accepted08 Jul 2012
Published10 Sept 2012
Abstract
For a fractional inverse problem with an unknown time-dependent source term, stability estimates are obtained by using operator theory approach. For the approximate solutions of the problem, the stable difference schemes which have first and second orders of accuracy are presented. The algorithm is tested in a one-dimensional fractional inverse problem.
1. Introduction
Inverse problems arise in many fields of science and engineering such as ion transport problems, chromatography, and heat determination problems with an unknown internal energy source. Different typed of inverse problems have been investigated, and the main results obtained in this field of research were given by many researchers (see [1โ10]). More than three centuries the theory of fractional derivatives developed mainly as a pure theoretical field of mathematics. Fractional integrals and derivatives appear in the theory of control of dynamical systems, when the controlled system or/and the controller is described by a fractional differential equation (see [11]). Recently, many application areas such as bioengineering applications, image and signal processing are also related to fractional calculus. Methods of solutions of problems and theory of fractional calculus have been studied by many researchers [11โ28]. Among them finite difference method is used for solving several fractional differential equations (see [20, 22, 23, 27] and the references therein).
1.1. Statement of the Problem
Many scientists and researchers are trying to enhance mathematical models of real-life cases for investigating and understanding the behavior of them. Therefore, some phenomena have been modeled and investigated as fractional inverse problems (see [29โ33] and the references therein). In this paper, we consider the fractional parabolic inverse problem with the Dirichlet condition
where and are unknown functions, , and is a sufficiently large number. Here, is the standard Riemann-Liouvilleโs derivative of order .
Theorems on the stability of problem (1.1) are analyzed by assuming that is a sufficiently smooth function, and .
2. Main Results
In this section, stability estimates for the solution of (1.1) are investigated. For the mathematical substantiation, we introduce the Banach space , of all continuous functions defined on with satisfying a Hรถlder condition for which the following norm is finite
where is the space of all continuous function defined on with the norm
With the help of a positive operator , we introduce the fractional spaces , consisting of all in a Banach space for which the following norm is finite:
Throughout the paper, positive constants will be indicated by . Here variables are used to focus on the fact that the constant depends only on and the subindex is used to indicate a different constant.
Theorem 2.1. Let , and . Then for the solution of problem (1.1), the following coercive stability estimates
hold.
Proof. Let us search for the solution of inverse problem (1.1) in the following form (see [8]):
where
Using the overdetermined condition, we get
Using identity (2.8) and the triangle inequality, it follows that
for any . Here, is the solution of the following problem:
For simplicity, we assign
where
Note that functions and only contain given functions. Then, we can rewrite problem (2.10) as
So, the end of proof of Theorem 2.1 is based on estimate (2.9) and the following theorem.
Theorem 2.2. For the solution of problem (2.10), the following coercive stability estimate
holds.
Proof. In a Banach space , with the help of the positive operator defined by
with
where is a positive constant, problem (2.10) can be written in the abstract form as an initial-value problem
By the Cauchy formula, the solution can be written as
Applying the formula
we get the following presentation of the solution of abstract problem (2.17):
Changing the order of integration, we obtain that
where
Now, we estimate separately. It is known that [13]
Since operators and commute,
Applying the definition of norm of the spaces and (2.23) and (2.24), we get
for any . Estimation of is as follows:
Let us estimate :
It is proven that (see [28])
Using the definition of norm of the spaces , we can obtain that
Using estimates (2.23) and (2.28), we get
Expanding , estimation of is as follows:
It is known that (see [34])
Since and are known functions, it is easy to see that
Estimation of can be given similar to the estimation of . By (2.23) and (2.32),
Finally combining estimates (2.25), (2.26), (2.30), (2.33), and (2.34), we get
Using the Gronwallโs inequality, we can write
From the last estimate, we can obtain the estimate for by using problem (2.17) and well-posedness of the Cauchy problem in (see [35]). So the following theorem finishes the proof of Theorem 2.2.
Theorem 2.3 (see, [36]). For , the norms of the spaces and are equivalent.
3. Numerical Results
We have not been able to obtain a sharp estimate for the constants figuring in the stability inequalities. So we will provide the following results of numerical experiments of the following problem:
The exact solution of the given problem is and for the control parameter is .
3.1. The First Order of Accuracy Difference Scheme
For the approximate solution of the problem (3.1), the Rothe difference scheme
where denotes greatest integer less than is constructed. Throughout the paper, let us denote
We search the solution of (3.2) in the following form:
where
Moreover for the interior grid point , we have that
From (3.4), (3.5), and the condition , it follows that
where is the solution of the difference scheme
First, applying the first order of accuracy difference scheme (3.10), we obtain system of linear equations and we write them in the matrix form
where
for any , and
Here, for any ,
and is identity matrix. Using (3.11), we can obtain that
To solve the resulting difference equations, we apply the method given in (3.15) step by step for . For the evaluation of is needed. It is obtained in the previous step. Then, the solution pairs are obtained by using the last formulas (3.9) and (3.8).
3.2. The Second Order of Accuracy Difference Scheme
For the approximate solution of the problem (3.1), the Crank-Nicholson difference scheme
is constructed.
Here,
Moreover, applying the second order of approximation formula for
it is obtained (see [27])
Here and throughout the paper,
We search the solution of (3.16) in the following form:
where
We have that
Let us denote
where . Then, one can write
So the values of can be obtained by the following formula:
Let denote
For , one can show that is the solution of the difference scheme
We have the system of linear equations and we write them in the matrix form
where
Here, for any ,
and is identity matrix. Using (3.29), we can obtain that
For ,โโ is the solution of the difference scheme
The system of linear equations given above can be written in the matrix form
where
Here, for any ,
Using (3.34), we can obtain that
For , we can obtain the following difference scheme:
This system can be written in matrix form as
where
for .
Here, for any ,
Finally, from (3.39), it follows that
Applying the last formula step by step, we can reach . Then, using (3.21), (3.25), and (3.26), we reach the approximate solutions of and .
3.3. Error Analysis
In this part, the results of the numerical analysis is given. The numerical solutions are recorded for different values of and and represents the approximate solution of at grid points . Table 1 gives the error analysis between the exact solution and the solutions derived by difference schemes. Table 1 is constructed for , and , respectively. For their comparison, the errors are computed by
Thus, the second order of accuracy difference scheme is more accurate comparing with the first order of accuracy difference scheme.
Acknowledgment
The authors are grateful to Professor Allaberen Ashyralyev (Fatih University, Turkey) for his comments and suggestions to improve the quality of the paper.
References
A. Hasanov, โIdentification of unknown diffusion and convection coefficients in ion transport problems from flux data: an analytical approach,โ Journal of Mathematical Chemistry, vol. 48, no. 2, pp. 413โ423, 2010.
A. Hasanov and S. Tatar, โAn inversion method for identification of elastoplastic properties of a beam from torsional experiment,โ International Journal of Non-Linear Mechanics, vol. 45, pp. 562โ571, 2010.
G. Di Blasio and A. Lorenzi, โIdentification problems for parabolic delay differential equations with measurement on the boundary,โ Journal of Inverse and Ill-Posed Problems, vol. 15, no. 7, pp. 709โ734, 2007.
D. Orlovsky and S. Piskarev, โOn approximation of inverse problems for abstract elliptic problems,โ Journal of Inverse and Ill-Posed Problems, vol. 17, no. 8, pp. 765โ782, 2009.
Y. S. Eidelman, โA boundary value problem for a differential equation with a parameter,โ Differentsia'nye Uravneniya, vol. 14, no. 7, pp. 1335โ1337, 1978.
A. Ashyralyev, โOn a problem of determining the parameter of a parabolic equation,โ Ukranian Mathematical Journal, vol. 62, no. 9, pp. 1200โ1210, 2010.
V. Serov and L. Päivärinta, โInverse scattering problem for two-dimensional Schrödinger operator,โ Journal of Inverse and Ill-Posed Problems, vol. 14, no. 3, pp. 295โ305, 2006.
V. T. Borukhov and P. N. Vabishchevich, โNumerical solution of the inverse problem of reconstructing a distributed right-hand side of a parabolic equation,โ Computer Physics Communications, vol. 126, no. 1-2, pp. 32โ36, 2000.
A. A. Samarskii and P. N. Vabishchevich, Numerical Methods for Solving Inverse Problems of Mathematical Physics, Inverse and Ill-Posed Problems Series, Walter de Gruyter, Berlin, Germany, 2007.
A. I. Prilepko and A. B. Kostin, โSome inverse problems for parabolic equations with final and integral observation,โ Matematicheskiĭ Sbornik, vol. 183, no. 4, pp. 49โ68, 1992.
I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional integrals and derivatives, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993.
A. A. Kilbas, H. M. Sristava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Mathematics Studies, Elsevier Science, 2006.
J.-L. Lavoie, T. J. Osler, and R. Tremblay, โFractional derivatives and special functions,โ SIAM Review, vol. 18, no. 2, pp. 240โ268, 1976.
V. E. Tarasov, โFractional derivative as fractional power of derivative,โ International Journal of Mathematics, vol. 18, no. 3, pp. 281โ299, 2007.
R. Gorenflo and F. Mainardi, โFractional calculus: integral and differential equations of fractional order,โ in Fractals and Fractional Calculus in Continuum Mechanics, A. Carpinteri and F. Mainardi, Eds., vol. 378 of CISM Courses and Lectures, pp. 223โ276, Springer, Vienna, Austria, 1997.
D. Matignon, โStability results for fractional differential equations with applications to control processing,โ in Computational Engineering in System Application, vol. 2, Lille, France, 1996.
F. Mainardi, โFractional calculus: some basic problems in continuum and statistical mechanics,โ in Fractals and Fractional Calculus in Continuum Mechanics, A. Carpinteri and F. Mainardi, Eds., vol. 378 of CISM Courses and Lectures, pp. 291โ348, Springer, New York, NY, USA, 1997.
A. Ashyralyev, F. Dal, and Z. Pinar, โOn the numerical solution of fractional hyperbolic partial differential equations,โ Mathematical Problems in Engineering, vol. 2009, Article ID 730465, 11 pages, 2009.
A. Ashyralyev, โA note on fractional derivatives and fractional powers of operators,โ Journal of Mathematical Analysis and Applications, vol. 357, no. 1, pp. 232โ236, 2009.
E. Cuesta, C. Lubich, and C. Palencia, โConvolution quadrature time discretization of fractional diffusion-wave equations,โ Mathematics of Computation, vol. 75, no. 254, pp. 673โ696, 2006.
P. E. Sobolevskii, โSome properties of the solutions of differential equations in fractional spaces,โ Trudy Naucno-Issledovatel'skogi Instituta Matematiki VGU, vol. 14, pp. 68โ74, 1975 (Russian).
G. Da Prato and P. Grisvard, โSommes d'opérateurs linéaires et équations différentielles opérationnelles,โ Journal de Mathématiques Pures et Appliquées, vol. 54, no. 3, pp. 305โ387, 1975.
A. Ashyralyev and Z. Cakir, โOn the numerical solution of fractional parabolic partial differential equations with the Dirichlet condition,โ in Proceedings of the 2nd International Symposium on Computing in Science and Engineering (ISCSE '11), M. Gunes, Ed., pp. 529โ530, Kusadasi, Aydın, Turkey, June 2011.
A. Ashyralyev, โWell-posedness of the Basset problem in spaces of smooth functions,โ Applied Mathematics Letters, vol. 24, no. 7, pp. 1176โ1180, 2011.
D. A. Murio and C. E. Mejía, โGeneralized time fractional IHCP with Caputo fractional derivatives,โ in Proceedings of the 6th International Conference on Inverse Problems in Engineering: Theory and Practice, vol. 135 of Journal of Physics: Conference Series, pp. 1โ8, Dourdan, Paris, France, 2008.
J. Cheng, J. Nakagawa, M. Yamamoto, and T. Yamazaki, โUniqueness in an inverse problem for a one-dimensional fractional diffusion equation,โ Inverse Problems, vol. 25, no. 11, pp. 1โ16, 2009.
J. Nakagawa, K. Sakamoto, and M. Yamamoto, โOverview to mathematical analysis for fractional diffusion equations—new mathematical aspects motivated by industrial collaboration,โ Journal of Math-for-Industry, vol. 2A, pp. 99โ108, 2010.
K. Sakamoto and M. Yamamoto, โInverse source problem with a final overdetermination for a fractional diffusion equation,โ Mathematical Control and Related Fields, vol. 1, no. 4, pp. 509โ518, 2011.
A. Ashyralyev, โA note on fractional derivatives and fractional powers of operators,โ Journal of Mathematical Analysis and Applications, vol. 357, no. 1, pp. 232โ236, 2009.
A. Ashyralyev and P. E. Sobolevskii, New Difference Schemes for Partial Differential Equations, vol. 148 of Operator Theory: Advances and Applications, Birkhäuser Verlag, Berlin, Germany, 2004.
A. Ashyralyev, โFractional spaces generated by the positive differential and difference operators in a Banach space,โ in Mathematical Methods in Engineering, K. Taş, J. A. Tenreiro Machado, and D. Baleanu, Eds., pp. 13โ22, Springer, Dordrecht, The Netherlands, 2007.