Brownian subordinators and fractional Cauchy problems

Fractional Cauchy problems replace the usual first-order time derivative by a fractional derivative. This paper develops classical solutions and stochastic analogues for fractional Cauchy problems in a bounded domain DR d with Dirichlet boundary conditions. Stochastic solutions are constructed via an inverse stable subordi- nator whose scaling index corresponds to the order of the fractional time derivative. Dirichlet problems corresponding to iterated Brow- nian motion in a bounded domain are then solved by establishing a correspondence with the case of a half-derivative in time. 1. Introduction. In this paper, we extend the approach of Meerschaert and Scheffler ( 23) and Meerschaert et al. (24) to fractional Cauchy problems on bounded domains. Our methods involve eigenfunction expansions, killed Markov processes and inverse stable subordinators. In a recent related paper (7), we establish a connection between fractional Cauchy problems with index β = 1/2 on an unbounded domain, and iterated Brownian motion (IBM), defined as Zt = B(|Yt|), where B is a Brownian motion with values in R d and Y is an independent one-dimensional Brownian motion. Since IBM is also the stochastic solution to a Cauchy problem involving a fourth-order derivative in space (2, 14), that paper also establishes a connection between certain higher-order Cauchy problems and their time-fractional analogues. More generally, Baeumer, Meerschaert and Nane (7) shows a connection between fractional Cauchy problems with β = 1/2 and higher-order Cauchy problems that involve the square of the generator. In the present paper, we

Over the last two decades, anomalous diffusion processes in which the mean squares variance grows slower or faster than that in a Gaussian process have found many applications. At a macroscopic level, these processes are adequately described by fractional differential equations, which involves fractional derivatives in time or/and space. The fractional derivatives describe either history mechanism or long range interactions of particle motions at a microscopic level. The new physics can change dramatically the behavior of the forward problems. For example, the solution operator of the time fractional diffusion diffusion equation has only limited smoothing property, whereas the solution for the space fractional diffusion equation may contain weak singularity. Naturally one expects that the new physics will impact related inverse problems in terms of uniqueness, stability, and degree of ill-posedness. The last aspect is especially important from a practical point of view, i.e., stably reconstructing the quantities of interest. In this paper, we employ a formal analytic and numerical way, especially the two-parameter Mittag-Leffler function and singular value decomposition, to examine the degree of ill-posedness of several ‘classical’ inverse problems for fractional differential equations involving a Djrbashian–Caputo fractional derivative in either time or space, which represent the fractional analogues of that for classical integral order differential equations. We discuss four inverse problems, i.e., backward fractional diffusion, sideways problem, inverse source problem and inverse potential problem for time fractional diffusion, and inverse Sturm–Liouville problem, Cauchy problem, backward fractional diffusion and sideways problem for space fractional diffusion. It is found that contrary to the wide belief, the influence of anomalous diffusion on the degree of ill-posedness is not definitive: it can either significantly improve or worsen the conditioning of related inverse problems, depending crucially on the specific type of given data and quantity of interest. Further, the study exhibits distinct new features of ‘fractional’ inverse problems, and a partial list of surprising observations is given below. (a) Classical backward diffusion is exponentially ill-posed, whereas time fractional backward diffusion is only mildly ill-posed in the sense of norms on the domain and range spaces. However, this does not imply that the latter always allows a more effective reconstruction. (b) Theoretically, the time fractional sideways problem is severely ill-posed like its classical counterpart, but numerically can be nearly well-posed. (c) The classical Sturm–Liouville problem requires two pieces of spectral data to uniquely determine a general potential, but in the fractional case, one single Dirichlet spectrum may suffice. (d) The space fractional sideways problem can be far more or far less ill-posed than the classical counterpart, depending on the location of the lateral Cauchy data. In many cases, the precise mechanism of these surprising observations is unclear, and awaits further analytical and numerical exploration, which requires new mathematical tools and ingenuities. Further, our findings indicate fractional diffusion inverse problems also provide an excellent case study in the differences between theoretical ill-conditioning involving domain and range norms and the numerical analysis of a finite-dimensional reconstruction procedure. Throughout we will also describe known analytical and numerical results in the literature.

We introduce a class of interesting stochastic processes based on Brownian-time processes. These are obtained by taking Markov processes and replacing the time parameter with the modulus of Brownian motion. They generalize the iterated Brownian motion (IBM) of Burdzy and the Markov snake of Le Gall, and they introduce new interesting examples. After defining Brownian-time processes, we relate them to fourth order parabolic partial differential equations (PDE’s). We then study their exit problem as they exit nice domains in $\mathbb{R}^d$ , and connect it to elliptic PDE’s. We show that these processes have the peculiar property that they solve fourth order parabolic PDE’s, but their exit distribution—at least in the standard Brownian time process case—solves the usual second order Dirichlet problem. We recover fourth order PDE’s in the elliptic setting by encoding the iterative nature of the Brownian-time process, through its exit time, in a standard Brownian motion. We also show that it is possible to assign a formal generator to these non-Markovian processes by giving such a generator in the half-derivative sense.

This paper explicitly computes the transition densities of a spectrally negative stable process with index greater than one, reflected at its infimum. First we derive the forward equation using the theory of sun-dual semigroups. The resulting forward equation is a boundary value problem on the positive half-line that involves a negative Riemann-Liouville fractional derivative in space, and a fractional reflecting boundary condition at the origin. Then we apply numerical methods to explicitly compute the transition density of this space-inhomogeneous Markov process, for any starting point, to any desired degree of accuracy. Finally, we discuss an application to fractional Cauchy problems, which involve a positive Caputo fractional derivative in time.

Coupling between hyperbolic and diffusive systems: A port-Hamiltonian formulation

In this article, homotopy analysis method is successfully applied to obtain the approximate analytical solutions of the characteristic Cauchy reaction-diffusion equation with fractional time derivative. The beauty of the article is the wonderful application of Caputo fractional order time derivative. The linear interactions of the merging populations are examined using perturbation theory and the method of matched asymptotic expansions. The solutions of the problem for different particular cases are presented graphically.

In this paper, we analyze the stability of the equilibrium point and Hopf bifurcation point in the three-component time-fractional differential equation, which describes the predator-prey interaction between different species. In the dynamics, the classical first-order derivative in time is modelled by either the Caputo or the Atangana-Baleanu fractional derivative of order α,0<α<1. We utilized a fractional version of the Adams-Bashforth formula to discretize these fractional derivatives in time. The results of the linear stability analysis presented are confirmed by computer simulation results.

Suppose that S is a subordinator with a nonzero drift and W is an independent 1-dimensional Brownian motion. We study the subordinate Brownian motion X defined by X t = W(S t ). We give sharp bounds for the Green function of the process X killed upon exiting a bounded open interval and prove a boundary Harnack principle. In the case when S is a stable subordinator with a positive drift, we prove sharp bounds for the Green function of X in (0, ∞ ), and sharp bounds for the Poisson kernel of X in a bounded open interval.

Numerical approximation of nonlinear fractional parabolic differential equations with Caputo–Fabrizio derivative in Riemann–Liouville sense

We describe an extension of previous results on degenerate abstract equations described by Favini and Yagi in their monograph cited in the references.This allows us to study degenerate differential equations with fractional derivative in time. Both direct and relative inverse problems are studied. Finally, various applications to degenerate partial differential equations are described.

Fractional partial differential equations with more than one fractional derivative in time describe some important physical phenomena, such as the telegraph equation, the power law wave equation, or the Szabo wave equation. In this paper, we consider two- and three-dimensional multi-term time and space fractional partial differential equations. The multi-term time-fractional derivative is defined in the Caputo sense, whose order belongs to the interval (1,2],(2,3],(3,4] or (0,m], and the space-fractional derivative is referred to as the fractional Laplacian form. We derive series expansion solutions based on a spectral representation of the Laplacian operator on a bounded region. Some applications are given for the two- and three-dimensional telegraph equation, power law wave equation and Szabo wave equation.

The paper presents a concise summary of main properties of generalized fractional derivatives, so-called convolution type derivatives with respect to Bernstein functions. Applications are considered to modeling time dependent random fields on the sphere as solutions to partial differential equations with the generalized fractional derivative in time and random initial condition.

To model wave propagation in inhomogeneous media with frequency-dependent power-law attenuation, it is needed to use the fractional powers of symmetric coercive elliptic operators in space and the Caputo tempered fractional derivative in time. The model studied in this paper is semilinear stochastic space-time fractional wave equations driven by infinite dimensional multiplicative Gaussian noise and fractional Gaussian noise, because of the potential fluctuations of the external sources. We first give a representation of the mild solution and some stability estimates for the homogeneous problem, and then prove the existence and uniqueness of the mild solution by using fixed point theorem. Finally, the decay and regularity theory of the solution are provided.

We consider initial/boundary value problems for the subdiffusion and diffusion-wave equations involving a Caputo fractional derivative in time. We develop two fully discrete schemes based on the piecewise linear Galerkin finite element method in space and convolution quadrature in time with the generating function given by the backward Euler method/second-order backward difference method, and establish error estimates optimal with respect to the regularity of problem data. These two schemes are first- and second-order accurate in time for both smooth and nonsmooth data. Extensive numerical experiments for two-dimensional problems confirm the convergence analysis and robustness of the schemes with respect to data regularity.

Some existence and regularity results for a non-local transport-diffusion equation with fractional derivatives in time and space