Fourier regularization for a final value time-fractional diffusion problem

The evolution process of fractional order describes some phenomenon of anomalous diffusion and transport dynamics in complex system. The equation containing time-fractional derivative provides a suitable mathematical model for describing such a process. The backward problem for this system, which means to recover the initial state for some slow diffusion process from its present status, is very hard to solve due to the nonlocal property of fractional derivative and the irreversibility of time. For this ill-posed problem, we construct a regularizing solution using the Fourier transform method. Both the a-priori choice strategy and the a-posteriori choice strategy for the regularizing parameter are given, with the convergence analysis on the regularizing solution. Numerical implementations are presented to show the validity of the proposed scheme.     

On the uniqueness and reconstruction for an inverse problem of the fractional diffusion process

Consider an inverse problem for the time-fractional diffusion equation in one dimensional spatial space. The aim is to determine the initial status and heat flux on the boundary simultaneously from heat measurement data given on the other boundary. Using the Laplace transform and the unique extension technique, the uniqueness for this inverse problem is proven. Then we construct a regularizing scheme for the reconstruction of boundary flux for known initial status. The convergence rate of the regularizing solution is established under some a priori information about the exact solution. Moreover, the initial distribution can also be recovered approximately from our regularizing scheme. Finally we present some numerical examples, which show the validity of the proposed reconstruction scheme.     

On a terminal value problem for a generalization of the fractional diffusion equation with hyper-Bessel operator

In this paper, we consider an inverse problem of recovering the initial value for a generalization of time-fractional diffusion equation, where the time derivative is replaced by a regularized hyper-Bessel operator. First, we investigate the existence and regularity of our terminal value problem. Then we show that the backward problem is ill-posed, and we propose a regularizing scheme using a fractional Tikhonov regularization method. We also present error estimates between the regularized solution and the exact solution using two parameter choice rules.     

Continuation theorem of fractional order evolutionary integral equations

The fractional order evolutionary integral equations have been considered by the first author in [6], the existence, uniqueness and some other properties of the solution have been proved. Here we study the continuation of the solution and its fractional order derivative. Also we study the generality of this problem and prove that the fractional order diffusion problem, the fractional order wave problem and the initial value problem of the equation of evolution are special cases of it. The abstract diffusion-wave problem will be given also as an application.     

Tikhonov regularization method for a backward problem for the inhomogeneous time-fractional diffusion equation

Fractional (nonlocal) diffusion equations replace the integer-order derivatives in space and time by their fractional-order analogs and they are used to model anomalous diffusion, especially in physics. In this paper, we study a backward problem for an inhomogeneous time-fractional diffusion equation with variable coefficients in a general bounded domain. Such a backward problem is of practically great importance because we often do not know the initial density of substance, but we can observe the density at a positive moment. The backward problem is ill-posed and we propose a regularizing scheme by using Tikhonov regularization method. We also prove the convergence rate for the regularized solution by using an a priori regularization parameter choice rule. Numerical examples illustrate applicability and high accuracy of the proposed method.     

Stable numerical differentiation for the second order derivatives

Consider a numerical differential problem, which aims to compute the second order derivative of a function stably from its given noisy data. For this ill-posed problem, we introduce the Lavrent′ev regularization scheme by reformulating this differentiation problem as an integral equation of the first kind. The advantage of this proposed scheme is that we can give the regularizing solution by an explicit integral expression, therefore it is easy to be implemented. The a-priori and a-posterior choice strategies for the regularization parameter are considered, with convergence analysis and error estimate of the regularizing solution for noisy data based on the integral operator decomposition. The validity of the proposed scheme is shown by several numerical examples.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

THE QUASI-BOUNDARY VALUE METHOD FOR IDENTIFYING THE INITIAL VALUE OF THE SPACE-TIME FRACTIONAL DIFFUSION EQUATION

In this article, we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation. This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem. We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule. Some numerical results in one-dimensional case and two-dimensional case show that our method is effcient and stable.     

Quenching phenomenon of a time‐fractional diffusion equation with singular source term

In this paper, a time‐fractional diffusion equation with singular source term is considered. The Caputo fractional derivative with order 0<α ?1 is applied to the temporal variable. Under specific initial and boundary conditions, we find that the time‐fractional diffusion equation presents quenching solution that is not globally well‐defined as time goes to infinity. The quenching time is estimated by using the eigenfunction of linear fractional diffusion equation. Moreover, by implementing a finite difference scheme, we give some numerical simulations to demonstrate the theoretical analysis. Copyright © 2017 John Wiley & Sons, Ltd.     

Analysis of a Multi-Term Variable-Order Time-Fractional Diffusion Equation and Its Galerkin Finite Element Approximation

In this paper, we study the well-posedness and solution regularity of a multi-term variable-order time-fractional diffusion equation, and then develop an optimal Galerkin finite element scheme without any regularity assumption on its true solution. We show that the solution regularity of the considered problem can be affected by the maximum value of variable-order at initial time $t = 0$. More precisely, we prove that the solution to the multi-term variable-order time-fractional diffusion equation belongs to $C^2([0,T])$ in time provided that the maximum value has an integer limit near the initial time and the data has sufficient smoothness, otherwise the solution exhibits the same singular behavior like its constant-order counterpart. Based on these regularity results, we prove optimal-order convergence rate of the Galerkin finite element scheme. Furthermore, we develop an efficient parallel-in-time algorithm to reduce the computational costs of the evaluation of multi-term variable-order fractional derivatives. Numerical experiments are put forward to verify the theoretical findings and to demonstrate the efficiency of the proposed scheme.     

The dependence on fractional orders of mild solutions to the fractional diffusion equation with memory

In this paper we investigate the Cauchy problem for a fractional diffusion equation and the time-fractional derivative is taken in the Caputo type sense. We give a representation of solutions under Fourier series and analyze initial value problems for the semi-linear fractional diffusion equation with a memory term. We also discuss the stability of the fractional derivative order for the time under some assumptions on the input data. Our key idea is to use Mittag-Leffler functions, the Banach fixed point theorem, and some Sobolev embeddings.     

Regularity results for fractional diffusion equations involving fractional derivative with Mittag–Leffler kernel

This paper studies partial differential equation model with the new general fractional derivatives involving the kernels of the extended Mittag–Leffler type functions. An initial boundary value problem for the anomalous diffusion of fractional order is analyzed and considered. The fractional derivative with Mittag–Leffler kernel or also called Atangana and Baleanu fractional derivative in time is taken in the Caputo sense. We obtain results on the existence, uniqueness, and regularity of the solution.     

High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

In this work, an effective and fast finite element numerical method with high-order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two-level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H¹-norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.     

分数阶微分方程边值问题研究简介

分数阶微积分是一个古老而又新颖的课题,近30年来,由于在包括分形现象在内的物理、工程等诸多应用学科领域应用的拓展,激发了科研人员对分数阶微积分的巨大热情。分数阶微分方程现在已应用于分数物理学、混沌与湍流、粘弹性力学与非牛顿流体力学、高分子材料的解链、自动控制理论、化学物理、随机过程和反常扩散等许多科学领域。分数阶微分方程边值问题是非线性常微分方程理论研究中一个活跃而成果丰硕的领域。本文讨论了分数阶微分方程边值问题的一些理论,介绍了作者的著作《分数阶微分方程边值问题理论及应用》的基本内容。     

A pseudo-spectral scheme for the approximate solution of a family of fractional differential equations

In this paper, the pseudo-spectral method is generalized for solving fractional differential equations with initial conditions. For this purpose, an appropriate representation of the solution is presented and the pseudo-spectral differentiation matrix of fractional order is derived. Then, by using pseudo-spectral scheme, the problem is reduced to the solution of a system of algebraic equations. Through several numerical examples, we evaluate the accuracy and performance of our proposed method.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

The space-time fractional diffusion equation with Caputo derivatives

We deal with the Cauchy problem for the space-time fractional diffusion equation, which is obtained from standard diffusion equation by replacing the second-order space derivative with a Caputo (or Riemann-Liouville) derivative of order β∈(0, 2] and the first-order time derivative with Caputo derivative of order α∈(0, 1]. The fundamental solution (Green function) for the Cauchy problem is investigated with respect to its scaling and similarity properties, starting from its Fourier-Laplace representation. We derive explicit expression of the Green function. The Green function also can be interpreted as a spatial probability density function evolving in time. We further explain the similarity property by discussing the scale-invariance of the space-time fractional diffusion equation.     

A fractional partial differential equation based multiscale denoising model for texture image

Traditional integer‐order partial differential equation based image denoising approach can easily lead edge and complex texture detail blur, thus its denoising effect for texture image is always not well. To solve the problem, we propose to implement a fractional partial differential equation (FPDE) based denoising model for texture image by applying a novel mathematical method—fractional calculus to image processing from the view of system evolution. Previous studies show that fractional calculus has some unique properties that it can nonlinearly enhance complex texture detail in digital image processing, which is obvious different with integer‐order differential calculus. The goal of the modeling is to overcome the problems of the existed denoising approaches by utilizing the aforementioned properties of fractional differential calculus. Using classic definition and property of fractional differential calculus, we extend integer‐order steepest descent approach to fractional field to implement fractional steepest descent approach. Then, based on the earlier fractional formulas, a FPDE based multiscale denoising model for texture image is proposed and further analyze optimal parameters value for FPDE based denoising model. The experimental results prove that the ability for preserving high‐frequency edge and complex texture information of the proposed fractional denoising model are obviously superior to traditional integral based algorithms, as for texture detail rich images. Copyright © 2013 John Wiley & Sons, Ltd.     

Singular fractional differential equations

Fractional differential equations have received increasing attention during recent years since the behavior of many physical systems can be properly described using the fractional order system theory. By fractional analog for Duhamel principle we give the existence-uniqueness result for linear and nonlinear time fractional evolution equations with singularities in corresponding norm in extended Colombeau algebra of generalized functions. In order to find the explicit solutions we use integral representation of the solution obtained via Laplace and Fourier transforms in succession and their inverses. We deal with some nonlinear models with singularities appearing in viscoelasticity and in anomalous processes, extending the results in viscoelasticity, continuum random walk, seismology, continuum mechanics and many other branches of life and science. The main task is finding existence-uniqueness results like in the case of evolution equations with entire derivatives. By examining the fractional evolution equations it turns out that they lead to till now known results from the evolution equations with entire derivatives in limiting case. They give more, behavior of the solution when order of derivatives are inside the intervals of entire points. In this way we can follow the influence of the operators generated by entire derivative in many fractional time evolution PDEs especially with singular initial data, and non-Lipschitz's nonlinear term. Apart from evolution equations we prove also an existence-uniqueness result for an initial value problem with singularities for linear and nonlinear fractional elliptic equation of Helmholtz type and fractional order α, where 1 < Re(α) ≤ 2, with respect to the one variable from R ₊. As a framework, we employ also Colombeau algebra of generalized functions containing fractional derivatives and operations among them in order to deal with the fractional equations with singularities. We apply the same techniques to the fractional Laplace and Poisson equation linear and nonlinear ones. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非光滑时间分数阶齐次扩散方程的修正

当初值不光滑时,时间分数阶齐次扩散方程数值方法的精度会下降.为了得到高阶时间收敛格式,提出加权移位的Grünwald-Letnikov的修正格式,运用Lubich的修正方法,得到非光滑时间分数阶齐次扩散方程的收敛阶仍为O(k²).最后,通过数值算例验证了数值计算结果与理论计算结果一致.