Research Article  Open Access
ZhiLiang Deng, XiaoMei Yang, XiaoLi Feng, "A Mollification Regularization Method for a FractionalDiffusion Inverse Heat Conduction Problem", Mathematical Problems in Engineering, vol. 2013, Article ID 109340, 9 pages, 2013. https://doi.org/10.1155/2013/109340
A Mollification Regularization Method for a FractionalDiffusion Inverse Heat Conduction Problem
Abstract
The illposed problem of attempting to recover the temperature functions from one measured transient data temperature at some interior point of a onedimensional semiinfinite conductor when the governing linear diffusion equation is of fractional type is discussed. A simple regularization method based on Dirichlet kernel mollification techniques is introduced. We also propose a priori and a posteriori parameter choice rules and get the corresponding error estimate between the exact solution and its regularized approximation. Moreover, a numerical example is provided to verify our theoretical results.
1. Introduction
It is well known that the classical diffusion equation requires second spatial derivatives and first time derivatives. However, people are shifting their partial focus to fractionalorder differential equations with the realization that the use of fractionalorder derivatives and integrals leads to formulas of certain physical processes which is more economical and useful than the classical approach in terms of Fick's laws of diffusion [1–4]. Some time fractional diffusion equations involving only a firstorder spatial derivative and a halforder time derivative are successfully used for modeling some anomalous diffusion physical phenomena (see, e.g., [5–7] and the references therein). One can also find some applications of such fractional diffusion equations in several scenarios in [2, 7], such as relaxation to equilibrium in systems (such as polymers chains and membranes) with long temporal memory, anomalous transport in disordered systems, diffusion on fractals, and to model nonMarkovian dynamical processes in protein folding. Fractional diffusion equations have been extensively investigated both in theory itself and numerical computation based on the broad applications in many application fields, especially in describing phenomena related to anomalous diffusion processes. The following is a partial list of articles which contain theoretical results and numerical tests. Some fundamental solutions and Green functions of fractional differential equations are given in [8, 9]. Finite difference techniques [10–12] and finite element methods [13, 14] provide some efficiency numerical results for several kinds of fractional differential equations. In [15], the decomposition method is used to construct analytical approximate solutions of timefractional wave equation subject to specified boundary conditions. An analytical solution of a fractional diffusion equation by Adomian decomposition method is presented in [16]. One can also refer [17–19] to the regularity of the solution of fractional diffusion equations and some a priori estimates.
In this paper, we consider the following onedimensional fractional diffusion problem on a semiinfinite slab: suppose the temperature at some interior point is approximately measurable. For convenience, we set . The temperature at the is desired and unknown. The mathematical description of the fractional diffusion problem is listed in the following. The unknown temperature satisfies where is the constant diffusivity coefficient, , a constant. The halftime differentiation indicates the RiemannLiouville fractional derivative with order which is given by the convolution integral where is the Gamma function. More detailed information on fractional derivatives and a general historical perspective may be found in [7, 20, 21]. The situations we will treat are those in which the system is initially at equilibrium, so that At a perturbation of the system commences by some unspecified process occurring at the boundary. During times of interest , this perturbation does not affect regions remote from the boundary, so that the relationship applies. In [22], Murio discussed the similar problems using mollification method with Gauss kernel. The idea used in this current work is a development of the ideas in [22]. However, there is only formal stability and discrete error discussion in the cited paper. Here, we give some analysis of error estimates under a priori and a posteriori regularization parameter, and a comparison of those two choice methods in our regularization method. In the following, the regularization method is outlined.
In order to simplify the Fourier analysis of fractional diffusion problem (1), and in the rest of the paper, we assume without loss of generality, , . We also assume that all the functions involved are and use the corresponding norm, as defined in the following, to measure errors: If the Fourier transform of a function is written as then the Fourier inversion formula reads Fourier transforming (1) with respect to variable , it follows that according to the fact that [7] The firstorder ordinary differential equation (8) has the general solution or equivalently, Now, we can write the problem (1) as where is a multiplication operator.
Since we require and , we see from (4) and (10) that the real part of is positive. Thereby, the multiplication operator is a bounded linear operator for . Then its inverse operator is unbounded. Therefore problem (12) is linear ill posed. For the general theory of linear ill posed problems, we can refer to [23]. It is worth pointing out that ill posed problems of a large number of diffusion equations, both fractionalorder as well as integral order, have been discussed by many authors. Yang et al. [24–28] discuss the identification of source terms for some integralorder diffusion equations using some regularization strategies. Hon et al. [29, 30] apply some meshless methods to the ill posed problems of heat conduction equations. In [17, 22, 29, 31, 32], some uniqueness results and numerical methods are given for some fractional diffusion ill posed problems. Here, we apply a simple stabilizing method, namely, the mollification method with Dirichlet kernel [33], to stabilize the problem (1). Suppose that the measured data function satisfies where is noise level. Take the Dirichlet function as the mollifier kernel, where is a positive constant. Define operator as Then, we have the following associated problem: for some , find that satisfies In fact, plays the role of regularization parameter.
In Section 2 an a priori parameter choice rule and the corresponding error estimate are discussed. In Section 3, we propose an a posteriori parameter choice rule and get the error bound. Numerical tests are given in Section 4.
2. The Error Estimate with a Priori Parameter Choice
In this section, the error estimate of the mollification regularization method will be derived under the a priori parameter choice rule. Suppose that the following source condition holds: The source condition (17) is indispensable, otherwise, there can be no uniform convergence rates for any regularization method, that is, the convergence is arbitrarily slow (see Proposition 3.11 in [23]). For , by Parseval formula and triangle inequality, we know that Since , for , we get by (13) For , we use the source condition (17) and obtain Therefore, Minimizing the righthand side of (21), we set and get which also implies Therefore, we get the following theorem.
Theorem 1. Assume conditions (13), (17) hold. If the regularization parameter is taken by (22), then for , there holds the error estimate (23).
Remark 2. The error estimate (23) is order optimal in the sense of Tautenhahn [34]. In our application is usually not known, therefore we have no exact a priori bound and cannot choose the parameter according to (22). However, if selecting , we can obtain the convergence rate
Theorem 1 provides no information about the convergence and convergence rates of at . The question is settled by our next result. We now give the error estimate at under a stronger a priori assumption
where denotes the norm on Sobolev space defined by
for . We only need to reestimate for . Under the stronger a priori bound (25), it is not hard to get
Then, (18), (19), and (27) lead to the following error bound:
The error bound (28) does not provide the convergence as obviously. Hence, we need to choose a proper parameter . For this reason, some proper lower bound of the righthand side of (28) should be given. In order to minimize the righthand side of (28), it is necessary to introduce the following lemma and its proof can be found in [34].
Lemma 3 (see [34]). Let the function be given by with a constant and positive constants , and , then for the inverse function one has
Minimizing the righthand side of (28), we let . Denote . Simple computation shows . By using (30) in Lemma 3, we obtain which also implies that Therefore, the following estimate holds where is a constant, which also shows that the convergence rate at is logarithmic.
Theorem 4. Assume conditions (13), (25) hold. If the regularization parameter is taken as given by (32), then the error estimate (33) holds.
Remark 5. From the error estimate (23), as , we see that the accuracy of regularization solution becomes progressively lower and even cannot get convergence for . This is common in the theory of ill posed problems. Nevertheless, if a stronger a priori assumption (25) is imposed, the regularization solution converges to the exact solution at , but only in a slower way.
3. The Error Estimate with A Posteriori Parameter Choice
In this section, we consider the a posteriori regularization parameter choice rule. Choose the regularization parameter as the solution of the equation where is a constant. To establish existence and uniqueness of solution of (34), we need the following lemma.
Lemma 6. If , then there hold the following:(a) is a continuous function;(b)(c);(d) is a strictly decreasing function.
The proof is very easy and we omit it here.
Denote We give the main result of this section as follows.
Theorem 7. Assume the conditions (11) and (17) hold and . Take the solution of (34) as the regularization parameter, then there holds the error estimate for : where is a constant depending on .
Proof. By (35), the Parseval formula, triangle inequality, and (17), we have By virtue of (34), we know that which means Thus it is obvious that Moreover, inserting (40) into (37), we get the following inequality: In addition, It is easy to see Denote . By using the Hölder inequality, we know that Therefore, we obtain which also implies that Therefore, we complete the proof.
4. Numerical Examples
For linear heat diffusion, analytic solutions for the temperature distribution , , in a semiinfinite solid with zero initial temperature, and at the surface, are obtained using the integral equation where the kernel function is the temperature distribution corresponding to a unit step boundary temperature, . Here, denotes the complementary error function defined by We take the example from [22]. As an interesting and challenging test for the numerical method, we proposed the surface temperature function , which is 1 between 0.2 and 0.6 and zero otherwise. The exact solutions for the FICHP, in this example, are the functions Consequently, the exact interior data temperature is given by Suppose the vector represents samples from the function . The noisy discrete data function is generated by adding random errors to the exact data function, at every grid point, that is, where the 's are independent random Gaussian variables with variance . The absolute and relative weighted errors for the recovered interior temperatures are calculated as respectively.
It is easy to implement the algorithm described in (16). Using an available fast Fourier transform (FFT) subroutine, a simple program was written to test the algorithm for the above example in Matlab. Firstly, some tests are implemented to verify the effectiveness of the regularization method under a priori parameter choice rule at several interior points. In addition, we would like to compare the a posteriori parameter choice rule (34) with the a priori parameter choice rule (22). The a priori bound in (22) can be calculated easily in this example. The Newton's bisection is used to solve (34), where we choose .
Figures 1 and 2 provide the comparisons between the exact solutions and regularization solutions with error bound and at interior point using a priori parameter choice rule (22), respectively. Figure 3 gives the comparison between the exact solutions and regularization solutions with error bound at boundary point using a priori parameter choice rule (32). And intuitively, it seems that there are better numerical effects for closer distance from . We list the error in Tables 1 and 2 to verify our result. Figure 4 demonstrates the comparisons between a priori and a posteriori parameter choice rules and Table 3 tells us that it has better effect using a posteriori parameter choice rule than a priori parameter choice rule. We also see that the difference of the numerical results between a priori choice rule and a posteriori choice rule is slight, which agrees with our theoretical results.



(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
Acknowledgments
The authors are indebted to the referees of this paper for their most helpful comments and suggestions, which helped to improve the presentation greatly. This work was supported by the Fundamental Research Funds for the Central Universities ZYGX2011J104 and SWJTU11BR078 and the NSF of China (Nos. 11126102, 11226040, and 11126187).
References
 V. V. Kulish and J. L. Lage, “Fractionaldiffusion solution for transient local temperature and heat flux,” Transactions of ASME, vol. 122, pp. 372–376, 2000. View at: Google Scholar
 K. B. Oldham and J. Spanier, The Fractional Calculus: Theory and Application of Differential and Integration to Arbitrary Order, Academic Press, 1974. View at: MathSciNet
 K. B. Oldham and J. Spanier, “A general solution of the diffusion equation for semiinfinite geometries,” Journal of Mathematical Analysis and Applications, vol. 39, pp. 655–669, 1972. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. B. Oldham and J. Spanier, “The replacement of Fick’s law by a formulation involving semidifferentiation,” Journal of Electroanalytical Chemistry, vol. 26, pp. 331–341, 1970. View at: Google Scholar
 R. Gorenflo, F. Mainardi, D. Moretti, and P. Paradisi, “Time fractional diffusion: a discrete random walk approach,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 129–143, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. V. Anh and N. N. Leonenko, “NonGaussian scenarios for the heat equation with singular initial conditions,” Stochastic Processes and their Applications, vol. 84, no. 1, pp. 91–114, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 I. Podlubny, Fractional Differential Equations, Academic Press, 1999.
 M. M. ElBorai, “The fundamental solutions for fractional evolution equations of parabolic type,” Boletín de la Asociación Matemática Venezolana, vol. 6, no. 1, pp. 29–43, 2004. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 F. Mainardi, Y. Luchko, and G. Pagnini, “The fundamental solution of the spacetime fractional diffusion equation,” Fractional Calculus & Applied Analysis, vol. 4, no. 2, pp. 153–192, 2001. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 F. Liu, V. Anhb, and I. Turnerb, “Numerical solution of the space fractional FokkerPlanck equation,” Journal of Computational and Applied Mathematics, vol. 166, pp. 209–219, 2004. View at: Google Scholar
 M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for fractional advectiondispersion flow equations,” Journal of Computational and Applied Mathematics, vol. 172, no. 1, pp. 65–77, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. B. Yuste, “Weighted average finite difference methods for fractional diffusion equations,” Journal of Computational Physics, vol. 216, no. 1, pp. 264–274, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. J. Fix and J. P. Roop, “Least squares finiteelement solution of a fractional order twopoint boundary value problem,” Computers & Mathematics with Applications, vol. 48, no. 78, pp. 1017–1033, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. P. Roop, “Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in ${R}^{2}$,” Journal of Computational and Applied Mathematics, vol. 193, no. 1, pp. 243–268, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. M. Odibat and S. Momani, “Approximate solutions for boundary value problems of timefractional wave equation,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 767–774, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. S. Ray and R. K. Bera, “Analytical solution of a fractional diffusion equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 174, no. 1, pp. 329–336, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. J. Liu and M. Yamamoto, “A backward problem for the timefractional diffusion equation,” Applicable Analysis, vol. 89, no. 11, pp. 1769–1788, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. McLean, “Regularity of solutions to a timefractional diffusion equation,” The ANZIAM Journal, vol. 52, no. 2, pp. 123–138, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Prüss, Evolutionary Integral Equations and Applications, vol. 87, Birkhäuser, Basel, Switzerland, 1993. View at: Publisher Site  MathSciNet
 A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, 2006. View at: MathSciNet
 S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Drivatives, Gordon and Breach Science, Yverdon, Switzerland, 1993. View at: MathSciNet
 D. A. Murio, “Stable numerical solution of a fractionaldiffusion inverse heat conduction problem,” Computers & Mathematics with Applications, vol. 53, no. 10, pp. 1492–1501, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. W. Engl, M. Hanke, and A. Neubauer, Regularization of Inverse Problems, Kluwer Academic, 1996. View at: Publisher Site  MathSciNet
 L. Yang, Z.C. Deng, J.N. Yu, and G.W. Luo, “Two regularization strategies for an evolutional type inverse heat source problem,” Journal of Physics A, vol. 42, no. 36, Article ID 365203, 16 pages, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Yang, Z.C. Deng, J.N. Yu, and G.W. Luo, “Optimization method for the inverse problem of reconstructing the source term in a parabolic equation,” Mathematics and Computers in Simulation, vol. 80, no. 2, pp. 314–326, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Yang, M. Dehghan, J.N. Yu, and G.W. Luo, “Inverse problem of timedependent heat sources numerical reconstruction,” Mathematics and Computers in Simulation, vol. 81, no. 8, pp. 1656–1672, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Yang, J. N. Yu, G. W. Luo, and Z. C. Deng, “Reconstruction of a space and time dependent heat source from finite measurement data,” International Journal of Heat and Mass Transfer, vol. 55, pp. 6573–6581, 2012. View at: Google Scholar
 L. Yang, J. N. Yu, G. W. Luo, and Z. C. Deng, “Numerical identification of source terms for a two dimensional heat conduction problem in polar coordinate system,” Applied Mathematical Modelling, vol. 37, pp. 939–957, 2013. View at: Google Scholar
 F. F. Dou and Y. C. Hon, “Kernelbased approximation for Cauchy problem of the timefractional diffusion equation,” Engineering Analysis with Boundary Elements, vol. 36, no. 9, pp. 1344–1352, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Y. C. Hon and T. Wei, “A fundamental solution method for inverse heat conduction problem,” Engineering Analysis With Boundary Elements, vol. 28, no. 5, pp. 489–495, 2004. View at: Google Scholar
 J. Cheng, J. Nakagawa, M. Yamamoto, and T. Yamazaki, “Uniqueness in an inverse problem for a onedimensional fractional diffusion equation,” Inverse Problems, vol. 25, no. 11, pp. 1–16, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Jin and W. Rundell, “An inverse problem for a onedimensional timefractional diffusion problem,” Inverse Problems, vol. 28, no. 7, Article ID 075010, 19 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. N. Hào, “A mollification method for illposed problems,” Numerische Mathematik, vol. 68, no. 4, pp. 469–506, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 U. Tautenhahn, “Optimality for illposed problems under general source conditions,” Numerical Functional Analysis and Optimization, vol. 19, no. 34, pp. 377–398, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2013 ZhiLiang Deng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.