Next Article in Journal
Special Types of Locally Conformal Closed G2-Structures
Next Article in Special Issue
Periodic Solution and Asymptotic Stability for the Magnetohydrodynamic Equations with Inhomogeneous Boundary Condition
Previous Article in Journal
Exponentially Harmonic Maps into Spheres
Previous Article in Special Issue
On the Shape Differentiability of Objectives: A Lagrangian Approach and the Brinkman Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem

by
Manuel D. Echeverry
and
Carlos E. Mejía
*,†
Escuela de Matemáticas, Universidad Nacional de Colombia, Carrera 65 # 59 A-110, Medellín 050034, Colombia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2018, 7(4), 89; https://doi.org/10.3390/axioms7040089
Submission received: 16 September 2018 / Revised: 13 November 2018 / Accepted: 20 November 2018 / Published: 23 November 2018
(This article belongs to the Special Issue Applications of Differential Equations and Dynamical Systems)

Abstract

:
We consider a two-dimensional time fractional diffusion equation and address the important inverse problem consisting of the identification of an ingredient in the source term. The fractional derivative is in the sense of Caputo. The necessary regularization procedure is provided by a two-dimensional discrete mollification operator. Convergence results and illustrative numerical examples are included.

1. Introduction

The mollification method has been known since the 1980s as a regularization procedure for ill-posed problems [1,2]. Later on, it was implemented as an accelerator of explicit schemes [2,3,4] and there are a variety of implementations and variations. Among the most recent, Reference [5] presents a new and very close method, called the convolution regularization method.
Applications of one-dimensional mollification methods are found, for instance, in References [6,7,8]. Mollification in two dimensions is less common and it is worth mentioning that it serves as a smoothing operator in Reference [9], appears as an accelerator of explicit schemes in Reference [4], and is a regularization tool for an identification problem in Reference [10].
In this paper, we deal with the two-dimensional discrete mollification operator defined in Section 3.2 of Reference [4]. This operator, which we call a 2DDM operator, will be defined shortly. Based on this operator, we propose the stable solution of a challenging inverse problem consisting of the identification of an ingredient in a source term for a time fractional diffusion equation, where the fractional derivative is defined in the sense of Caputo.
This inverse problem was solved in References [11,12,13]. The first two papers implemented the Tikhonov method as the required regularization procedure. The third one introduced a modified quasi-boundary regularization method. Our main original contribution is the implementation of two-dimensional discrete mollification as a regularization method for the numerical solution of this inverse source problem. It is important to have alternatives and this is the case with regularization. We offer an alternative that has not been considered before for the stable solution of this inverse problem.
The solution of the inverse problem requires the solution of the corresponding direct problem. In Reference [12], there are no details on the numerical solution, although the authors provide abundant numerical results. In Reference [11], the numerical solution is obtained by a finite difference scheme, but the only examples are one-dimensional. In our numerical solution of the direct problem, we make use of the implicit scheme defined in Reference [14] for a more general equation.
Problems based on time fractional diffusion equations are useful in the description of anomalous diffusion in highly heterogeneous porous media. There are a variety of time fractional inverse source problems. According to the number of fractional differentiation orders, they are classified as single-term and multi-term problems. This paper deals with a single-term problem that is ill-posed in the Hadamard sense. However, there are other situations, for instance in Reference [15], where the problem is single-term and depends on a parameter. For some values of this parameter, it is a well-posed problem in the Hadamard sense. Finally, we mention that in Reference [16], there is a theoretical analysis and some numerical results for a fairly general multiterm time fractional inverse source problem.
The paper is organized as follows: Section 2 is devoted to the 2DDM operator. Section 3 deals with the convergence analysis for the solution of the inverse problem, and Section 4 includes the numerical experiments. After this section, there are some final comments and conclusions.

2. The 2D Discrete Mollification Operator (2DDM Operator)

2.1. One-Dimensional Discrete Mollification

The one-dimensional mollification of a discrete function [1,2,4] is defined as follows: Let y = { y j } j Z be a discrete function, which can, for example, be the evaluations or cell averages of a real function y = y ( x ) at equidistant grid points on:
X = x j : x j = x 0 + j h , j Z .
The discrete mollification of y , denoted J η y , is:
( J η y ) j : = i = - η η w i y j + i ,
where η is the integer support parameter.

2.2. Mollification Weights

The weights are defined as integrals over cells. Given positive real numbers p and h, compute:
δ = η + 1 / 2 h p .
This is the continuous mollification parameter. The weights are obtained by numerical integration of the truncated Gaussian kernel:
κ p δ ( t ) = A p δ - 1 exp ( - t 2 / δ 2 ) , t p δ 0 , t > p δ ,
where the normalization constant:
A p = - p p exp ( - s 2 ) d s - 1
is chosen in such a way that R κ p δ = 1 . Thus, the mollification weights are obtained by integrations of this kernel over cells, namely,
w i = t i t i + 1 κ p δ ( s ) d s ,
where:
t j = j - 1 / 2 h , j Z .
We usually take p = 3 . Chosen in this way, the weights w i satisfy:
w i = w - i , 0 w i w i - 1 , i = 1 , , η ; i = - η η w i = 1 , i = - η η i w i = 0 .
The set of weights has a feature worth mentioning: Independence of the step size h.
w i = t i t i + 1 κ δ p ( s ) d s = A p δ - 1 t i t i + 1 exp ( - s 2 δ 2 ) d s = A p t i / δ t i + 1 / δ exp ( - s 2 ) d s
and
t i δ = p i - 1 / 2 h η + 1 / 2 h = p i - 1 / 2 η + 1 / 2 ,
then the w i s depend only on p and η . Some weights are shown in Table 1.

2.3. Two-Dimensional Discrete Mollification

Likewise [4], two-dimensional discrete mollification requires a uniform bidimensional grid of the form ( x i , y j ) where x i = x 0 + i h x and y j = y 0 + j h y . Let Y be a discrete function defined by Y ( x i , y j ) = Y i j . For each positive integer η , the two dimensional discrete mollification of Y is:
J η 2 Y i j = k = - η η m = - η η w k w m Y i + m , j + k .
The weights of this two dimensional convolution are denoted by w k m = w k w m . If h x = h y = h and δ is given by Equation (2), the shape of a weight is shown in Figure 1.
The main relation with one dimensional mollification is [4]:
J η 2 Y i j = k = - η η m = - η η w k w m Y i + m , j + k = k = - η η w k m = - η η w m Y i + m , j + k = k = - η η w k J η x Y i , j + k = J η y J η x Y i j .

2.4. Abstract Setting

The discrete convolutions that define mollification in one and two dimensions have their counterparts in infinite dimensions. It is true that in numerical computations, one deals with discrete functions, but what about the error estimates? They relate the abstract model based on differential equations with the discretized problem.
Convolutions are troublesome near borders; thus, whenever a convolution is mentioned, some kind of domain restriction is in order. Our presentation of the abstract setting follows Reference [17]. Let I = [ 0 , 1 ] × [ 0 , 1 ] , I p = [ - p , p ] × [ - p , p ] , I p δ = [ - p δ , p δ ] × [ - p δ , p δ ] and I δ = [ p δ , 1 - p δ ] × [ p δ , 1 - p δ ] . If s = ( s 1 , s 2 ) , the two-dimensional truncated Gaussian kernel is defined by multiplication of the known one dimensional kernels, that is:
ϕ δ p ( s 1 , s 2 ) = κ δ p x ( s 1 ) κ δ p y ( s 2 ) = A p δ - 2 exp ( - ( s 1 2 + s 2 2 ) / δ 2 ) , s I p δ 0 , o t h e r w i s e
where A p = I p exp ( - s 2 2 ) d s - 1 . This kernel satisfies:
  • ϕ 0 ,
  • ϕ 0 outside I p δ ,
  • R 2 ϕ δ p = 1 .
Definition 1.
A function g : R 2 R is called Lipschitz continuous, if there exists a constant C L i p such that:
| g ( x 1 ) - g ( x 2 ) | C L i p x 1 - x 2
The next lemma establishes the consistency and stability of the two-dimensional discrete mollification.
Lemma 1.
Suppose g is a Lipschitz continuous function with Lipschitz constant L defined on I and G = G i j is the associated discrete function defined on a uniform grid covering I δ with step size h in both directions. Then, there exists a constant C such that:
1. 
Consistency
J η 2 G - g , I δ C δ + h .
Proof. 
This is Lemma 3.1 Part (1) of Reference [17]. □
Remark 1.
Notice that due to Equation (2), we may say that the right hand side bound is of the form C δ , where C depends on L , p , and η.
2. 
Stability: If g ϵ is a Lipschitz continuous approximation of g so that g ϵ - g , I ϵ , then:
J η 2 g ϵ - g 2 , I δ C δ + ϵ .
Remark 2.
Actually, the right hand side bound is of the form C ( h + δ + ϵ ) , according to the previous remark.
Proof. 
It follows from Part 1 that since:
J η 2 g ϵ - g 2 2 = I δ | J η 2 g ϵ - g | 2 d x d y J η 2 g ϵ - g , I δ 2 μ ( I δ ) J η 2 g ϵ - g ϵ , I δ + g ϵ - g , I δ 2 C δ + ϵ 2 ,
where μ stands for measure. Notice that this estimate is valid in any other square subset of R 2 . Only the constant may change. □

3. The Inverse Source Problem

Fractional derivatives have been known since the seventeenth century, but only recently have they become the subject of attention due in part to a variety of interesting applications. Podlubny [18] says: “Fractional derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes.” We mention only two examples: For materials, we chose the damage evolution of brain tissue modeled in Reference [19] by a model in which the time derivative is a Caputo–Almeida fractional derivative. For processes, we selected the anomalous diffusion of a contaminant in fractured porous rocks modeled in Reference [20]. The model is based on fractional differentiation in the sense of Caputo.

3.1. Direct Problem

Our concern is a time fractional diffusion equation. More precisely, we are interested in an initial/boundary value problem, in which the time derivative is a Caputo fractional derivative of order α , 0 < α < 1 . The Caputo fractional derivative of order α of a differentiable function v is:
D t α v ( t ) = 1 Γ ( 1 - α ) 0 t v ( ξ ) ( t - ξ ) α d ξ , 0 < α < 1 , v ( t ) , α = 1 ,
where Γ is the Gamma function given by:
Γ ( z ) = 0 x z - 1 e - x d x .
We are interested in the following initial/boundary value problem:
D t α u ( z , t ) - ( L u ) ( z , t ) = p ( t ) f ( z ) , z Ω t ( 0 , T ) u ( z , t ) = 0 , z Ω t ( 0 , T ) u ( z , 0 ) = 0 , z Ω ¯
where Ω R d and:
L u ( z ) = i = 1 d z i j = 1 d a i j ( z ) z j u ( z ) + c ( z ) u ( z ) , z Ω .
The coefficients satisfy a i j = a j i C 1 ( Ω ¯ ) , i , j = 1 d a i j ζ i ζ j > θ i = 1 d | ζ i | 2 ( θ > 0 ), c ( z ) 0 , c C ( Ω ¯ ) . The theoretical analysis that follows is from References [11,13,21]. The operator - L is a symmetric uniformly elliptic operator defined on D ( - L ) = H 2 ( Ω ) H 0 1 ( Ω ) .
The direct problem is to find u that safisfies Equation (13), with p and f known functions and α a known fractional order of differentiation. If we think of this problem as the model for a contaminant diffusion in groundwater, then f ( z ) plays the role of contaminant discharge intensity and p ( t ) is an attenuation coefficient.
Let 0 < λ 1 λ 2 λ n be the eigenvalues of - L , with associated eigenfunctions { χ n } n N . It is known that l i m n λ n = and { χ n } n N is an orthonormal basis of L 2 ( Ω ) .
The first step is the consideration of a fractional power operator ( - L ) γ , γ R and the space:
D ( - L ) γ = ψ L 2 ( Ω ) : n = 1 λ n 2 γ | ( ψ , χ n ) | 2 < ,
where ( · , · ) is the inner product in L 2 ( Ω ) . It is a Hilbert space with the norm:
ψ D ( - L ) γ = n = 1 λ n 2 γ | ( ψ , χ n ) | 2 1 2 .
The Mittag–Leffler function is a generalization of the exponential function and is very useful when dealing with fractional derivatives.
Definition 2.
A two-parameter function of the Mittag–Leffler type is defined by:
E α , β ( z ) = k = 0 z k Γ ( α k + β ) ,
where α > 0 , β > 0 are real constants and Γ is the Gamma function defined by Equation (12).
The formal solution of the direct problem (13) can be written in terms of a Mittag–Leffler function. It is:
u ( z , t ) = n = 1 f n 0 t p ( τ ) t - τ α - 1 E α , α - λ n t - τ α d τ χ n ( z ) ,
where f n = ( f , χ n ) .
From Equation (18):
u ( z , T ) = n = 1 f n 0 T p ( τ ) T - τ α - 1 E α , α - λ n T - τ α d τ χ n ( z ) , = q ( z ) .
Denote:
Q n ( t ) = 0 t p ( τ ) t - τ α - 1 E α , α - λ n t - τ α d τ .
Thus:
q ( z ) = n = 1 f n Q n ( T ) χ n ( z ) ,
and its Fourier coefficients are q n = f n Q n ( T ) . We conclude f n = q n Q n ( T ) .

3.2. Inverse Problem

We are interested in the inverse problem of determining the unknown source ingredient f ( z ) in Equation (13) based on the overposed data u ( z , T ) = q ( z ) . Time T is a final time and the data q ( t ) are not known exactly; what is known is a measurement q ϵ so that q - q ϵ < ϵ .
As in Reference [11], the identification of f can be formulated as the operator equation:
( K f ) ( z ) = n = 1 f n 0 T p ( τ ) T - τ α - 1 E α , α - λ n T - τ α d τ χ n ( z ) = q ( z ) .
The operator K is a linear self-adjoint compact operator and Problem (19) is ill-posed.
Our aim is to obtain a function f δ ϵ that minimizes the functional:
H f = q - J η 2 q ϵ 2 2 ,
where J η 2 q ϵ is the mollification of a discrete version of the noisy data q ϵ . Since K is a self-adjoint compact operator, there exists a minimizer of H, f δ ϵ , given by:
f δ ϵ ( z ) = n = 1 J η 2 q ϵ n Q n ( T ) χ n ( z ) , w h e r e J η 2 q ϵ n = J η 2 q ϵ , χ n .
Theorem 1.
([13] Remark 3.1) If p ( t ) C [ 0 , T ] satisfying p ( t ) p 0 > 0 for all t [ 0 , T ] , set p C [ 0 , T ] = sup t [ t [ 0 , T ] ] | p ( t ) | . Then, we have:
p 0 C ̲ λ n Q n ( T ) = 0 T p ( τ ) ( T - τ ) α - 1 E α , α ( - λ n ( T - τ ) α ) d τ p C [ 0 , T ] λ n ,
where C ̲ = 1 - E α , 1 ( - λ 1 T α ) .
Theorem 2.
Suppose p C [ 0 , T ] satisfies p ( t ) p 0 > 0 , t [ 0 , T ] . Moreover, suppose q ϵ - q < ϵ and there exist m R + and E, a constant such that f D ( ( - L ) m 2 ) E , then:
f δ ϵ - f 2 C 1 δ + ϵ m / ( m + 2 ) ,
where C 1 depends on α , T , λ 1 , μ ( Ω ) .
Proof. 
f δ ϵ D ( ( - L ) m 2 ) 2 = n = 1 λ n m ( f δ ϵ , χ n ( z ) ) 2 = n = 1 λ n m ( f δ ϵ ) n 2 = n = 1 λ n m J η 2 q ϵ n 2 Q n ( T ) 2 n = 1 λ n m + 2 C ̲ 2 p 0 2 ( J η 2 q ϵ , χ n ( z ) ) 2 = 1 C ̲ 2 p 0 2 n = 1 λ n m + 2 ( J η 2 q ϵ , χ n ( x , y ) ) 2 = 1 C ̲ 2 p 0 2 J η 2 q ϵ D ( ( - L ) m 2 + 1 ) 2 .
The inequality follows from Theorem 1. However, J η 2 q ϵ C ( Ω ) , which implies f δ ϵ D ( ( - L ) m 2 ) . Taking M = max { J η 2 q ϵ D ( ( - L ) m 2 + 1 ) , E } , we obtain:
K f δ ϵ - q 2 K f δ ϵ - J η 2 q ϵ 2 + J η 2 q ϵ - q 2 K f - J η 2 q ϵ 2 + J η 2 q ϵ - q 2 = 2 J η 2 q ϵ - q 2 C δ + ϵ . f δ ϵ - f D ( ( - L ) m / 2 ) | f δ ϵ D ( ( - L ) m / 2 ) + f D ( ( - L ) m / 2 ) M 1 + 1 C ̲ p 0 .
Following Reference [11] and taking into consideration Lemma 1, Part 1, we conclude:
f δ ϵ - f 2 ( C ̲ p 0 ) - m m + 2 K f δ ϵ - q 2 m / ( m + 2 ) f δ ϵ - f D ( ( - L ) m / 2 ) 2 / ( m + 2 ) ( C ̲ p 0 ) - m m + 2 M 1 + 1 C ̲ p 0 2 / ( m + 2 ) C δ + ϵ m / ( m + 2 )

4. Numerical Results

There are many alternatives for the numerical solution of a time fractional diffusion problem. If we restrict our attention to methods based on finite difference discretizations, we mention Reference [22], in which Caputo fractional derivatives are computed by the so-called Caputo fractional derivative rule, which in turn depends on a modified trapezoidal rule. This method is successfully implemented in Reference [23] for plane strain/stress elasticity in the presence of nonlocal elasticity.
Our finite difference method for the solution of the time fractional diffusion equation is taken from Reference [14], which in turn is inspired by Reference [24]. The method is implicit and its properties about stability and convergence are discussed in Reference [14].
We consider the following time fractional diffusion problem:
D t α u ( z , t ) - ( L u ) ( z , t ) = p ( t ) f ( z ) , z Ω t ( 0 , T ) u ( z , t ) = 0 , z Ω t ( 0 , T ) u ( z , 0 ) = 0 , z Ω ¯ t ( 0 , T )
where:
L u ( x , y ) = x a 11 ( x , y ) x u ( x , y ) + y a 22 ( x , y ) y u ( x , y ) .
After discretization, operator L leads to a matrix l, whose eigenvectors and eigenvalues play the roles of the eigenfunctions and eigenvalues of the differential operator.
The Q n s are approximated according to Reference [11] by:
Q n = - 1 λ n E α , 1 ( - λ n τ α ) p ( T - τ ) | τ = T - p ( T - t ) | τ = 0 - 0 T E α , 1 ( - λ n τ α ) p ( T - τ ) d τ .
The integrals are approximated by composite trapezoidal rule.
The errors are measured by the root mean square (rms) norm given by:
f - f δ ϵ r m s = 1 N M n = 1 N m = 1 M ( f n m - ( f δ ϵ ) n m ) 2 1 / 2 .
The particular examples considered here are for:
f ( x , y ) = s i n ( 3 π x ) s i n ( 3 π y )
p ( t ) = 1 + e x p ( - ( t + α ) ) .
The space step size h = 1 / 32 is used for both coordinates, x and y. Final time is T = 1 and Δ t = 1 / 10 .

Example 1

Let a 11 ( x , y ) = x 2 + 1 , a 22 ( x , y ) = y 2 + 1 . The relative l 2 norms of the identification are shown in Table 2 and Table 3.

Example 2

Let a 11 ( x , y ) 1 and a 22 ( x , y ) 1 . Relative l 2 errors for the identification of the source term are shown in Table 4 and Table 5. Figure 2 gives an idea of the quality of the identification.

5. Discussion and Final Comments

The 2DDM operator is a useful regularization procedure when dealing with inverse source problems for time fractional diffusion equations. It is also appropriate for surface and gradient fitting as shown in Reference [17], but we do not pursue this objective here. There are more applications of this operator and we hope to master some of them; for instance, we would like to modify our identification problem and deal with the factor p ( t ) as unknown and the factor f ( x , y ) as data.

Author Contributions

Conceptualization, M.E. and C.M.; Investigation, M.E. and C.M.; Writing—Original Draft Preparation, M.E. and C.M.; Writing—Review & Editing, M.E. and C.M.

Funding

Manuel D. Echeverry would like to thank COLCIENCIAS convocatoria-647 for financial support through a Ph.D. study grant.

Acknowledgments

The authors acknowledge the support provided by Universidad Nacional de Colombia through the research project Nonlocal operators and operators with memory, Hermes code number 33154.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Acosta, C.D.; Mejía, C.E. Stable Computations by Discrete Mollification; Universidad Nacional de Colombia: Bogotá, Colombia, 2014. [Google Scholar]
  2. Murio, D.A. Mollification and Space Marching. In Inverse Engineering Handbook; Woodbury, K., Ed.; CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar]
  3. Acosta, C.D.; Mejía, C.E. Stabilization of explicit methods for convection diffusion equations by discrete mollification. Comput. Math. Appl. 2008, 55, 368–380. [Google Scholar] [CrossRef]
  4. Acosta, C.D.; Buerger, R. Difference schemes stabilized by discrete mollification for degenerate parabolic equations in two space dimensions. IMA J. Numer. Anal. 2012, 32, 1509–1540. [Google Scholar] [CrossRef]
  5. Shi, C.; Wang, C.; Wei, T. Convolution regularization method for backward problems of linear parabolic equations. Appl. Numer. Math. 2016, 108, 143–156. [Google Scholar] [CrossRef]
  6. Garshasbi, M.; Dastour, H. Estimation of unknown boundary functions in an inverse heat conduction problem using a mollified marching scheme. Numer. Algorithms 2015, 68, 769–790. [Google Scholar] [CrossRef]
  7. Mejía, C.E.; Piedrahita, A. Solution of a time fractional inverse advection-dispersion problem by discrete mollification. Rev. Colomb. Mat. 2017, 51, 83–102. [Google Scholar] [CrossRef]
  8. Yang, F.; Fu, C.; Li, X. A mollification regularization method for unknown source in time-fractional diffusion equation. Int. J. Comput. Mat. 2014, 91, 1516–1534. [Google Scholar] [CrossRef]
  9. Beatson, R.; Bui, H.Q. Mollification formulas and implicit smoothing. Adv. Comput. Mat. 2007, 27, 125–149. [Google Scholar] [CrossRef]
  10. Yi, Z.; Murio, D.A. Identification of Source Terms in 2-D IHCP. Comput. Math. Appl. 2004, 47, 1517–1533. [Google Scholar] [CrossRef]
  11. Ma, Y.; Prakash, P.; Deiveegan, A. Generalized Tikhonov methods for an inverse source problem of the time-fractional diffusion equation. Chaos Solitons Fract. 2018, 108, 39–48. [Google Scholar] [CrossRef]
  12. Nguyen, H.T.; Le, D.L.; Nguyen, V.T. Regularized solution of an inverse source problem for a time fractional diffusion equation. Appl. Math. Model. 2016, 40, 8244–8264. [Google Scholar] [CrossRef]
  13. Wei, T.; Wang, J. A modified quasi-boundary value method for an inverse source problem of the time-fractional diffusion equation. Appl. Numer. Math. 2014, 78, 95–111. [Google Scholar] [CrossRef]
  14. Mejía, C.E.; Piedrahita, A. A finite difference approximation of a two dimensional time fractional advection-dispersion problem. arXiv 2018, arXiv:1807.07393v1. [Google Scholar]
  15. Sakamoto, K.; Yamamoto, M. Inverse source problem with a final overdetermination for a fractional diffusion equation. Math. Control Relat. Fields 2011, 1, 509–518. [Google Scholar] [CrossRef] [Green Version]
  16. Jiang, D.; Li, Z.; Liu, Y.; Yamamoto, M. Weak unique continuation property and a related inverse source problem for time-fractional diffusion-advection equations. Inverse Probl. 2017, 33, 055013. [Google Scholar] [CrossRef] [Green Version]
  17. Zhan, S.; Murio, D.A. Surface fitting and numerical gradient computation by discrete mollification. Comput. Math. Appl. 1999, 37, 145–159. [Google Scholar] [CrossRef]
  18. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  19. Voyiadjis, G.Z.; Sumelka, W. Brain modelling in the framework of anisotropic hyperelasticity with time fractional damage evolution governed by the Caputo-Almeida fractional derivative. J. Mech. Behav. Biomed. Mater. 2019, 89, 209–216. [Google Scholar] [CrossRef] [PubMed]
  20. Fomin, S.; Chugunov, V.; Hashida, T. Application of fractional differential equations for modeling the anomalous diffusion of contaminant from fracture into porous rock matrix with bordering alteration zone. Transp. Porous Med. 2010, 81, 187–205. [Google Scholar] [CrossRef]
  21. Sakamoto, K.; Yamamoto, M. Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems. J. Math. Anal. Appl. 2011, 382, 426–447. [Google Scholar] [CrossRef]
  22. Odibat, Z. Approximations of fractional integrals and Caputo fractional derivatives. Appl. Math. Comput. 2006, 178, 527–533. [Google Scholar] [CrossRef]
  23. Sumelka, W.; Szajek, K.; Lodygowski, T. Plain strain and plane stress elasticity under fractional continuum mechanics. Arch. Appl. Mech. 2015, 85, 1527–1544. [Google Scholar] [CrossRef]
  24. Zhuang, P.; Liu, F. Finite difference approximation for two-dimensional time fractional diffusion equation. J. Algorithms Comput. Technol. 2007, 1, 1–16. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A weight for δ = 0 . 015 and η = 4 .
Figure 1. A weight for δ = 0 . 015 and η = 4 .
Axioms 07 00089 g001
Figure 2. Identification of f for h = 1 / 32 , ϵ = 0 . 025 and η = 3 .
Figure 2. Identification of f for h = 1 / 32 , ϵ = 0 . 025 and η = 3 .
Axioms 07 00089 g002
Table 1. Mollification weights.
Table 1. Mollification weights.
η w 0 w 1 w 2 w 3 w 4 w 5
18.4272 × 10 - 1 7.8640 × 10 - 2
26.0387 × 10 - 1 1.9262 × 10 - 1 5.4438 × 10 - 3
34.5556 × 10 - 1 2.3772 × 10 - 1 3.3291 × 10 - 2 1.2099 × 10 - 3
43.6266 × 10 - 1 2.4003 × 10 - 1 6.9440 × 10 - 2 8.7275 × 10 - 3 4.7268 × 10 - 4
53.0028 × 10 - 1 2.2625 × 10 - 1 9.6723 × 10 - 2 2.3430 × 10 - 2 3.2095 × 10 - 3 2.4798 × 10 - 4
Table 2. Error norms, ϵ = 0 . 025 and η = 5 .
Table 2. Error norms, ϵ = 0 . 025 and η = 5 .
α 0.10.50.70.9
f - f δ ϵ r m s 0.112372580.113415730.114829340.11338589
Table 3. Error norms, α = 0 . 7 .
Table 3. Error norms, α = 0 . 7 .
ϵ 0.0050.010.0250.05
η 3456
f - f δ ϵ r m s 0.056343420.111415610.176831410.24996681
Table 4. Error norms, ϵ = 0 . 025 and η = 5 .
Table 4. Error norms, ϵ = 0 . 025 and η = 5 .
α 0.10.50.70.9
f - f δ ϵ r m s 0.186798460.183303550.191290970.19253567
Table 5. Error norms, α = 0 . 7 .
Table 5. Error norms, α = 0 . 7 .
ϵ 0.0050.010.0250.05
η 3456
f - f δ ϵ r m s 0.07158180.120908870.192954880.26167402

Share and Cite

MDPI and ACS Style

Echeverry, M.D.; Mejía, C.E. A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem. Axioms 2018, 7, 89. https://doi.org/10.3390/axioms7040089

AMA Style

Echeverry MD, Mejía CE. A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem. Axioms. 2018; 7(4):89. https://doi.org/10.3390/axioms7040089

Chicago/Turabian Style

Echeverry, Manuel D., and Carlos E. Mejía. 2018. "A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem" Axioms 7, no. 4: 89. https://doi.org/10.3390/axioms7040089

APA Style

Echeverry, M. D., & Mejía, C. E. (2018). A Two Dimensional Discrete Mollification Operator and the Numerical Solution of an Inverse Source Problem. Axioms, 7(4), 89. https://doi.org/10.3390/axioms7040089

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