• The Laplacian $\Delta$ is the infinitesimal generator of isotropic Brownian motion, being the limit process of normal diffusion, while the fractional Laplacian $\Delta^{\beta/2}$ serves as the infinitesimal generator of the limit process of isotropic L\'{e}vy process. Taking limit, in some sense, means that the operators can approximate the physical process well after sufficient long time. We introduce the nonlocal operators (being effective from the starting time), which describe the general processes undergoing normal diffusion. For anomalous diffusion, we extend to the anisotropic fractional Laplacian $\Delta_m^{\beta/2}$ and the tempered one $\Delta_m^{\beta/2,\lambda}$ in $\mathbb{R}^n$. Their definitions are proved to be equivalent to an alternative one in Fourier space. Based on these new nonlocal diffusion operators, we further derive the deterministic governing equations of some interesting statistical observables of the very general jump processes with multiple internal states. Finally, we consider the associated initial and boundary value problems and prove their well-posedness of the Galerkin weak formulation in $\mathbb{R}^n$. To obtain the coercivity, we claim that the probability density function $m(Y)$ should be nondegenerate.
  • In this paper, we study the numerical schemes for the two-dimensional Fokker-Planck equation governing the probability density function of the tempered fractional Brownian motion. The main challenges of the numerical schemes come from the singularity in the time direction. When $0<H<0.5$, a change of variables $\partial \left(t^{2H}\right)=2Ht^{2H-1}\partial t$ avoids the singularity of numerical computation at $t=0$, which naturally results in nonuniform time discretization and greatly improves the computational efficiency. For $0.5<H<1$, the time span dependent numerical scheme and nonuniform time discretization are introduced to ensure the effectiveness of the calculation and the computational efficiency. By numerically solving the corresponding Fokker-Planck equation, we obtain the mean squared displacement of stochastic processes, which conforms to the characteristics of the tempered fractional Brownian motion.
  • Topological defects resulted from boundary constraints in confined liquid crystals have attracted extensive research interests. In this paper, we use numerical simulation to study the phase transition dynamics in the context of stochastic resonance in a bistable liquid crystal device containing defects. This device is made of nematic liquid crystals confined in a shallow square well, and is described by the planar Lebwohl-Lasher model. The stochastic phase transition processes of the system in the presence of a weak oscillating potential is simulated using an over-damped Langevin dynamics. Our simulation results reveal that, depending on system size, the phase transition may follow two distinct pathways: in small systems the pre-existing defect structures at the corners hold until the last stage and there is no newly formed defect point in the bulk during the phase transition; In large systems new defect points appear spontaneously in the bulk and eventually merge with the pre-existing defects at the corners. For both transition pathways stochastic resonance can be observed, but show dramatic difference in their responses to the boundary anchoring strength. In small systems we observe a "sticky-boundary" effect for a certain range of anchoring strength in which the phase transition gets stuck and stochastic resonance becomes de-activated. Our work demonstrates the dynamical interplay among defects, noises, and boundary conditions in confined liquid crystals.
  • Applying proper orthogonal decomposition to a usual finite element (FE) formulation for space fractional partial differential equation, we get a reduced FE model, which greatly reduces the complexity of computation. Then, the stability analysis and error estimate for the reduced model are presented. Finally, we verify the effectiveness of the algorithm by numerical experiments.
  • Tempered fractional Laplacian is the generator of the tempered isotropic L\'evy process [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., 16(1), 125-149, 2018]. This paper provides the finite difference discretization for the two dimensional tempered fractional Laplacian $(\Delta+\lambda)^{\frac{\beta}{2}}$. Then we use it to solve the tempered fractional Poisson equation with Dirichlet boundary conditions and derive the error estimates. Numerical experiments verify the convergence rates and effectiveness of the schemes.
  • Based on the theory of continuous time random walks (CTRW), we build the models of characterizing the transitions among anomalous diffusions with different diffusion exponents, often observed in natural world. In the CTRW framework, we take the waiting time probability density function (PDF) as an infinite series in three parameter Mittag-Leffler functions. According to the models, the mean squared displacement of the process is analytically obtained and numerically verified, in particular, the trend of its transition is shown; furthermore the stochastic representation of the process is presented and the positiveness of the PDF of the position of the particles is strictly proved. Finally, the fractional moments of the model are calculated, and the analytical solutions of the model with external harmonic potential are obtained and some applications are proposed.
  • The generalized Langevin equation describes anomalous dynamics. Noise is not only the origin of uncertainty but also plays a positive role in helping to detect signal with information, termed stochastic resonance (SR). This paper analyzes the anomalous resonant behaviors of the generalized Langevin system with a multiplicative dichotomous noise and an internal tempered Mittag-Leffler noise. For the system with fluctuating harmonic potential, we obtain the exact expressions of several SR, such as, the first moment, the amplitude and the autocorrelation function for the output signal as well as the signal-noise ratio. We analyze the influence of the tempering parameter and memory exponent on the bona fide SR and the general SR. Moreover, it is detected that the critical memory exponent changes regularly with the increase of tempering parameter. Almost all the theoretical results are validated by numerical simulations.
  • Levy walk is a fundamental model with applications ranging from quantum physics to paths of animal foraging. Taking animal foraging as an example, a natural idea that comes to one's mind is to introduce the multiple internal states for dealing with the dependence of the PDF of waiting time on the energy of the animal and richness of the food at a particular location, etc; the framework can also be used to model the moving trajectories of smart animals without returning to the directions or locations which they come from immediately. After building the Levy walk model with multiple internal states and deriving the governing equation of the distribution of the positions of the particles, some applications are discussed with specific transition matrices. The type of diffusion for non-immediately-repeating L\'{e}vy walk is uncovered, and the distribution and average of first passage time are numerically simulated.
  • The fractional Laplacian $\Delta^{\beta/2}$ is the generator of $\beta$-stable L\'evy process, which is the scaling limit of the L\'evy fight. Due to the divergence of the second moment of the jump length of the L\'evy fight it is not appropriate as a physical model in many practical applications. However, using a parameter $\lambda$ to exponentially temper the isotropic power law measure of the jump length leads to the tempered L\'evy fight, which has finite second moment. For short time the tempered L\'evy fight exhibits the dynamics of L\'evy fight while after sufficiently long time it turns to normal diffusion. The generator of tempered $\beta$-stable L\'evy process is the tempered fractional Laplacian $(\Delta+\lambda)^{\beta/2}$ [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., in press, 2017]. In the current work, we present new computational methods for the tempered fractional Laplacian equation, including the cases with the homogeneous and nonhomogeneous generalized Dirichlet type boundary conditions. We prove the well-posedness of the Galerkin weak formulation and provide convergence analysis of the single scaling B-spline and multiscale Riesz bases finite element methods. We propose a technique for efficiently generating the entries of the dense stiffness matrix and for solving the resulting algebraic equation by preconditioning. We also present several numerical experiments to verify the theoretical results.
  • Aging, the process of growing old or maturing, is one of the most widely seen natural phenomena in the world. For the stochastic processes, sometimes the influence of aging can not be ignored. For example, in this paper, by analyzing the functional distribution of the trajectories of aging particles performing anomalous diffusion, we reveal that for the fraction of the occupation time $T_+/t$ of strong aging particles, $\langle (T^+(t)^2)\rangle=\frac{1}{2}t^2$ with coefficient $\frac{1}{2}$, having no relation with the aging time $t_a$ and $\alpha$ and being completely different from the case of weak (none) aging. In fact, we first build the models governing the corresponding functional distributions, i.e., the aging forward and backward Feynman-Kac equations; the above result is one of the applications of the models. Another application of the models is to solve the asymptotic behaviors of the distribution of the first passage time, $g(t_a,t)$. The striking discovery is that for weakly aging systems, $g(t_a,t)\sim t_a^{\frac{\alpha}{2}}t^{-1-\frac{\alpha}{2}}$, while for strongly aging systems, $g(t_a,t)$ behaves as $ t_a^{\alpha-1}t^{-\alpha}$.
  • This paper focuses on the adaptive discontinuous Galerkin (DG) methods for the tempered fractional (convection) diffusion equations. The DG schemes with interior penalty for the diffusion term and numerical flux for the convection term are used to solve the equations, and the detailed stability and convergence analyses are provided. Based on the derived posteriori error estimates, the local error indicator is designed. The theoretical results and the effectiveness of the adaptive DG methods are respectively verified and displayed by the extensive numerical experiments. The strategy of designing adaptive schemes presented in this paper works for the general PDEs with fractional operators.
  • In this paper we derive new uniform convergence estimates for the V-cycle MGM applied to symmetric positive definite Toeplitz block tridiagonal matrices, by also discussing few connections with previous results. More concretely, the contributions of this paper are as follows: (1) It tackles the Toeplitz systems directly for the elliptic PDEs. (2) Simple (traditional) restriction operator and prolongation operator are employed in order to handle general Toeplitz systems at each level of the recursion. Such a technique is then applied to systems of algebraic equations generated by the difference scheme of the two-dimensional fractional Feynman-Kac equation, which describes the joint probability density function of non-Brownian motion. In particular, we consider the two coarsening strategies, i.e., doubling the mesh size (geometric MGM) and Galerkin approach (algebraic MGM), which lead to the distinct coarsening stiffness matrices in the general case: however, several numerical experiments show that the two algorithms produce almost the same error behaviour.
  • This paper provides a theoretical framework of deriving the forward and backward Feynman-Kac equations for the distribution of functionals of the path of a particle undergoing both diffusion and chemical reaction. Very general forms of the equations are obtained. Once given the diffusion type and reaction rate, a specific forward or backward Feynman-Kac equation can be obtained. The listed in the paper include the ones for normal/anomalous diffusions and reactions with linear/nonlinear rates. Using the derived equations, we also study the occupation time in half-space, the first passage time to a fixed boundary, and the occupation time in half-space with absorbing or reflecting boundary conditions.
  • This paper further discusses the tempered fractional Brownian motion, its ergodicity, and the derivation of the corresponding Fokker-Planck equation. Then we introduce the generalized Langevin equation with the tempered fractional Gaussian noise for a free particle, called tempered fractional Langevin equation (tfLe). While the tempered fractional Brownian motion displays localization diffusion for the long time limit and for the short time its mean squared displacement has the asymptotic form $t^{2H}$, we show that the asymptotic form of the mean squared displacement of the tfLe transits from $t^2$ (ballistic diffusion for short time) to $t^{2-2H}$, and then to $t^2$ (again ballistic diffusion for long time). On the other hand, the overdamped tfLe has the transition of the diffusion type from $t^{2-2H}$ to $t^2$ (ballistic diffusion). The tfLe with harmonic potential is also considered.
  • We provide and analyze the high order algorithms for the model describing the functional distributions of particles performing anomalous motion with power-law jump length and tempered power-law waiting time. The model is derived in [Wu, Deng, and Barkai, Phys. Rev. E., 84 (2016), 032151], being called the time-tempered fractional Feynman-Kac equation. The key step of designing the algorithms is to discretize the time tempered fractional substantial derivative, being defined as $${^S\!}D_t^{\gamma,\widetilde{\lambda}} G(x,p,t)\!=\!D_t^{\gamma,\widetilde{\lambda}} G(x,p,t)\!-\!\lambda^\gamma G(x,p,t) ~{\rm with}~\widetilde{\lambda}=\lambda+ pU(x),\, p=\rho+J\eta,\, J=\sqrt{-1},$$ where $$D_t^{\gamma,\widetilde{\lambda}} G(x,p,t) =\frac{1}{\Gamma(1-\gamma)} \left[\frac{\partial}{\partial t}+\widetilde{\lambda} \right] \int_{0}^t{\left(t-z\right)^{-\gamma}}e^{-\widetilde{\lambda}\cdot(t-z)}{G(x,p,z)}dz,$$ and $\lambda \ge 0$, $0<\gamma<1$, $\rho>0$, and $\eta$ is a real number. The designed schemes are unconditionally stable and have the global truncation error $\mathcal{O}(\tau^2+h^2)$, being theoretically proved and numerically verified in {\em complex} space. Moreover, some simulations for the distributions of the first passage time are performed, and the second order convergence is also obtained for solving the `physical' equation (without artificial source term).
  • For the particles undergoing the anomalous diffusion with different waiting time distributions for different internal states, we derive the Fokker-Planck and Feymann-Kac equations, respectively, describing positions of the particles and functional distributions of the trajectories of particles; in particular, the equations governing the functional distribution of internal states are also obtained. The dynamics of the stochastic processes are analyzed and the applications, calculating the distribution of the first passage time and the distribution of the fraction of the occupation time, of the equations are given.
  • For characterizing the Brownian motion in a bounded domain: $\Omega$, it is well-known that the boundary conditions of the classical diffusion equation just rely on the given information of the solution along the boundary of a domain; on the contrary, for the L\'evy flights or tempered L\'evy flights in a bounded domain, it involves the information of a solution in the complementary set of $\Omega$, i.e., $\mathbb{R}^n\backslash \Omega$, with the potential reason that paths of the corresponding stochastic process are discontinuous. Guided by probability intuitions and the stochastic perspectives of anomalous diffusion, we show the reasonable ways, ensuring the clear physical meaning and well-posedness of the partial differential equations (PDEs), of specifying `boundary' conditions for space fractional PDEs modeling the anomalous diffusion. Some properties of the operators are discussed, and the well-posednesses of the PDEs with generalized boundary conditions are proved.
  • The mean first exit (passage) time characterizes the average time of a stochastic process never leaving a fixed region in the state space, while the escape probability describes the likelihood of a transition from one region to another for a stochastic system driven by discontinuous (with jumps) L\'evy motion. This paper discusses the two deterministic quantities, mean first exit time and escape probability, for the anomalous processes having the tempered L\'{e}vy stable waiting times with the tempering index $\lambda>0$ and the stability index $0<\alpha \le 1$; as for the distribution of jump lengths or the type of the noises driving the system, two cases are considered, i.e., Gaussian white noise and non-Gaussian (tempered) $\beta$-stable ($0<\beta<2$) L\'{e}vy noise. Firstly, we derive the nonlocal elliptic partial differential equations (PDEs) governing the mean first exit time and escape probability. Based on the derived PDEs, it is observed that the mean first exit time depends strongly on the domain size and the values of $\alpha$, $\beta$ and $\lambda$; when $\lambda$ is close to zero, the mean first exit time tends to $\infty$. In particular, we also find an interesting result that the escape probability of a particle with (tempered) power-law jumping length distribution has no relation with the distribution of waiting times for the model considered in this paper. For the solutions of the derived PDEs, the boundary layer phenomena are observed, which inspires the motivation for developing the boundary layer theory for nonlocal PDEs.
  • This paper focuses on providing the high order algorithms for the space-time tempered fractional diffusion-wave equation. The designed schemes are unconditionally stable and have the global truncation error $\mathcal{O}(\tau^2+h^2)$, being theoretically proved and numerically verified.
  • For describing the probability distribution of the positions and times of particles performing anomalous motion, fractional PDEs are derived from the continuous time random walk models with waiting time distribution having divergent first order moment and/or jump length distribution which has divergent second order moment. It can be noted that the fractional PDEs are essentially dealing with the multiscale issues. Generally the regularity of the solutions for fractional PDEs is weak at the areas close to boundary and initial time. This paper focuses on developing the applications of wavelet bases to numerically solving fractional PDEs and digging out the potential benefits of wavelet methods comparing with other numerical methods, especially in the aspects of realizing preconditioning, adaptivity, and keeping the Toeplitz structure. More specifically, the contributions of this paper are as follows: 1. the techniques of efficiently generating stiffness matrix with computational cost $\mathcal{O}(2^J)$ are provided for first, second, and any order bases; 2. theoretically and numerically discuss the effective preconditioner for time-independent equation and multigrid method for time-dependent equation, respectively; 3. the wavelet adaptivity is detailedly discussed and numerically applied to solving the time-dependent (independent) equations. In fact, having reliable, simple, and local regularity indicators is the striking benefit of the wavelet in adaptively solving fractional PDEs (it seems hard to give a local posteriori error estimate for the adaptive finite element method because of the global property of the operator).
  • Recently, nonlocal models attract the wide interests of scientist. They mainly come from two applied scientific fields: peridyanmics and anomalous diffusion. Even though the matrices of the algebraic equation corresponding the nonlocal models are usually Toeplitz (denote a0 as the principal diagonal element, a1 as the trailing diagonal element, etc). There are still some differences for the models in these two fields. For the model of anomalous diffusion, a0/a1 is uniformly bounded; most of the time, a0/a1 of the model for peridyanmics is unbounded as the stepsize h tends to zero. Based on the uniform boundedness of a0/a1, the convergence of the two-grid method is well established [Chan, Chang, and Sun, SIAM J. Sci. Comput., 19 (1998), pp. 516--529; Pang and Sun, J. Comput. Phys., 231 (2012), pp. 693--703; Chen, Wang, Cheng, and Deng, BIT, 54 (2014), pp. 623--647]. This paper provides the detailed proof of the convergence of the two-grid method for the nonlocal model of peridynamics. Some special cases of the full multigrid and the V-cycle multigrid are also discussed. The numerical experiments are performed to verify the convergence.
  • This paper focuses on providing the computation methods for the backward time tempered fractional Feynman-Kac equation, being one of the models recently proposed in [Wu, Deng, and Barkai, Phys. Rev. E, 84 (2016) 032151]. The discretization for the tempered fractional substantial derivative is derived, and the corresponding finite difference and finite element schemes are designed with well established stability and convergence. The performed numerical experiments show the effectiveness of the presented schemes.
  • Because of the finiteness of the life span and boundedness of the physical space, the more reasonable or physical choice is the tempered power-law instead of pure power-law for the CTRW model in characterizing the waiting time and jump length of the motion of particles. This paper focuses on providing the variational formulation and efficient implementation for solving the corresponding deterministic/macroscopic models, including the space tempered fractional equation and time tempered fractional equation. The convergence, numerical stability, and a series of variational equalities are theoretically proved. And the theoretical results are confirmed by numerical experiments.
  • Fractional operators are widely used in mathematical models describing abnormal and nonlocal phenomena. Although there are extensive numerical methods for solving the corresponding model problems, theoretical analysis such as the regularity result, or the relationship between the left-side and right-side fractional operators are seldom mentioned. In stead of considering the fractional derivative spaces, this paper starts from discussing the image spaces of Riemann-Liouville fractional integrals of $L_p(\Omega)$ functions, since the fractional derivative operators that often used are all pseudo-differential. Then high regularity situation---the image spaces of Riemann-Liouville fractional integral operators on $W^{m,p}(\Omega)$ space are considered. Equivalent characterizations of the defined spaces, as well as of the intersection of the left-side and right-side spaces are given. The behavior of the functions in the defined spaces at both the nearby boundary point/ponits and the points in the domain are demonstrated in a clear way. Besides, tempered fractional operators show to be reciprocal to the corresponding Riemann-Liouville fractional operators, which is expected to make some efforts on theoretical support for relevant numerical methods. Last, we also provide some instructions on how to take advantage of the introduced spaces when numerically solving fractional equations.
  • The traditional wave equation models wave propagation in an ideal conducting medium. For characterizing the wave propagation in inhomogeneous media with frequency dependent power-law attenuation, the space-time fractional wave equation appears; further incorporating the additive white Gaussian noise coming from many natural sources leads to the stochastic space-time fractional wave equation. This paper discusses the Galerkin finite element approximations for the stochastic space-time fractional wave equation forced by an additive space-time white noise. We firstly discretize the space-time additive noise, which introduces a modeling error and results in a regularized stochastic space-time fractional wave equation; then the regularity of the regularized equation is analyzed. For the discretization in space, the finite element approximation is used and the definition of the discrete fractional Laplacian is introduced. We derive the mean-squared $L^2$-norm priori estimates for the modeling error and for the approximation error to the solution of the regularized problem; and the numerical experiments are performed to confirm the estimates. For the time-stepping, we calculate the analytically obtained Mittag-Leffler type function.