Advances in Mathematical Physics
Volume 2020 (2020), Article ID 1315426, 7 pages
https://doi.org/10.1155/2020/1315426

Research Article

Stochastic Representation and Monte Carlo Simulation for Multiterm Time-Fractional Diffusion Equation

Longjin Lv and Luna Wang

School of Finance and Information, Ningbo University of Finance and Economics, Ningbo 315000, China

Correspondence should be addressed to Luna Wang; wangluna@nbufe.edu.cn

Received 8 May 2020; Accepted 19 June 2020; Published 14 July 2020

Academic Editor: Laurent Raymond

Copyright © 2020 Longjin Lv and Luna Wang. 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.

Abstract.

In this paper, we mainly study the solution and properties of the multiterm time-fractional diffusion equation. First, we obtained the stochastic representation for this equation, which turns to be a subordinated process. Based on the stochastic representation, we calculated the mean square displacement (MSD) and time average mean square displacement, then proved some properties of this model, including subdiffusion, generalized Einstein relationship, and nonergodicity. Finally, a stochastic simulation algorithm was developed for the visualization of sample path of the abnormal diffusion process. The Monte Carlo method was also employed to show the behavior of the solution of this fractional equation.

1. Introduction

Recently, the diffusion equations that generalize the usual one have received considerable attention due to the broadness of their physical applications, in particular, to the anomalous diffusion. In fact, fractional diffusion equations and the nonlinear fractional diffusion equations have been successfully applied to several physical situations such as percolation of gases through porous media [1], thin saturated regions in porous media [2], standard solid-on-solid model for surface growth [3], thin liquid films spreading under gravity [4], in the transport of fluid in porous media and in viscous fingering [5], modeling of non-Markovian dynamical processes in protein folding [6], relaxation to equilibrium in a system (such as polymer chains and membranes) with long temporal memory [7], and anomalous transport in disordered systems [8], diffusion on fractals [9], and the multiphysical transport in porous media, such as electroosmosis [10, 11]. Moreover, some underlying processes can be more accurately and flexibly modeled by multiterm FPDEs. For example, the multiterm time-fractional diffusion-wave equation is a satisfying mathematical model for viscoelastic damping [12]. In [13], a two-term fractional diffusion equation has been successfully used for distinguishing different states in solute transport.

The multiterm time-fractional advection-diffusion equations are linear integrodifferential equations. They are obtained from their corresponding classical multiterm advection-diffusion equations by replacing the first-order time derivative by fractional derivative, which reads

where with ; here, the is the Caputo fractional derivative of order with respect to as defined as

There are growing interests in studying these equations because of their importance in modeling many physical, biological, medical, chemical, and many other fields. For example, the over diagnostic ultrasound frequencies, acoustic absorption in biological tissue exhibits a power law with a noninteger frequency [14, 15]. Also, in a complex inhomogeneous conducting medium, experimental evidence shows that the sound waves propagate with the power law of noninteger order. For further applications on physics and on real phenomena [16–18], the Caputo time-fractional operator has been widely used instead of the second time derivative to model mathematically such problems in order to discuss the effect of the memory on the studied system [19].

Many analytic and numeric methods are employed to solve this equation. Daftardar-Gejji and Bhalekar considered the multiterm time-fractional diffusion-wave equation using the method of separation of variables [20]. Luchko [21] studied the well-posedness of the multiterm time-fractional diffusion equation based on an approximate maximum principle. Jiang et al. studied the multiterm time-space fractional advection-diffusion equation based on the spectral representation of the fractional Laplacian operator [22]. Meshless analysis based on improved moving least-squares approximation was introduced to solve the two-dimensional two-sided space-fractional wave equation in [23]. Using the method of series expansion, Ye et al. [24] studied the multiterm time-space fractional partial differential equations in 2D and 3D domains. An efficient operational formulation of the spectral tau method for a multiterm time-space fractional differential equation with Dirichlet boundary conditions was proposed in [25]. Liu et al. presented numerical approximations for multiterm time-fractional diffusion equations by using the spectral method in [26] and for multiterm time-fractional wave equations by means of FDMs in [27], respectively. In [28], the authors used finite difference rules to get the approximate solutions of the time-fractional multiterm wave equations. Recently, the stochastic representation method was introduced to solve the fractional diffusion equation. In [29], Kolokoltsov built the relation between the stochastic process and time-fractional diffusion equations with Caputo or Riemann-Liouville derivatives. These generalized Caputo derivatives were further extended to nonmonotone processes, yielding two-sided and even multidimensional extensions. Based on the stochastic representation, the mathematical properties of the related fractional equation were discussed in [30–32]. These papers inspired the research of this paper.

In this paper, we introduce the stochastic representation method to solve this multiterm time-fractional diffusion equation. This paper is organized as follows. In Section 2, we derive a subordinated process, whose PDF is rightly the solution of this equation, where the parent process is a classical diffusion process and the subordinator is the inverse time of the sum of Lévy motions with different parameter. Taking advantage of this result, we study the properties of this multiterm time-fractional diffusion equation in Section 3. We also employ the Monte Carlo method to simulate the solution for this equation in the next section. Section 5 presents our conclusions.

2. Stochastic Representation

In this section, we will give the stochastic representation of the multiterm time-fractional diffusion equation.

Let be the increasing Lévy motion with Laplace transform , then we get the following theorem.

Theorem 1. The subordinated process is the stochastic representation of the multiterm time-fractional advection-diffusion equation (1), where the parent process and subordinator of is defined as

respectively. Here, is independent of and , are also independent with each other for different .

Proof. Following the same procedure shown in [33], we first establish the relation between the PDF of and the PDF of . From the definition of (see Equation (4)), we have [34], therefore

So, the Laplace transform of can be expressed as

Here, the is obtained from the definition of . Then, using the total probability formula and the independence between and , we can get the PDF of , given by

where is the PDF of the parent process . So the Laplace transform of the above equation yields

Thus, the following relation between and holds

Since the process is given by the Itô stochastic differential equation (3), its PDF obeys the classical advection-dispersion equation [35]

The Laplace transform of the above equation with respect to yields

By changing the variable to and using the relation (9), the above equation yields

Since the Laplace transform of the Caputo fractional derivative is given by , by comparing the above equation with the Laplace transform of Equation (3), we get . So the subordinated process is called the stochastic representation of the multiterm time-fractional diffusion equation.

3. Some Properties

The superiority of the stochastic representation approach to the fractional differential equation is that it not only helps us to understand the physical process by providing a description of the dynamical system governed by the fractional differential equation but also provides a way to get the properties of the corresponding equations. For simplicity, we retain two terms for the time-fractional operator, i.e., .

Corollary 2. The subordinated process governed by the multiterm time-fractional diffusion equation (1) is subdiffusive.

Proof. Taking advantage of the relation between and , we can get the following evaluation formula for the mean square displacement

where is the PDF of and its Laplace transformation is given by Equation (4). By inverting the Laplace transform, we can get

Here, , the PDF of inverse time -stable Lévy motion, can be expressed in the form of a Fox function, i.e.,

In order to calculate the mean square displacement, we can first compute its Laplace transform. By using the Laplace transform of Equation (4), and inverting the Laplace transform, we have

where is the Mittag-Leffler function [36]. Here, we have used the equality

From Equation (16), we can know the model resembles a subdiffusion for , and since , the model resembles a subdiffusion for . This result can be got in another way; to see this, note that in distribution and grows faster than for , so the -stable subordinator dominates as and the -stable subordinator dominates as .

Corollary 3. The generalized Einstein relation holds for the multiterm time-fractional diffusion equation (1).

Proof. With the help of the stochastic representation, we can calculate the first moment of of Equation (1) in the presence of a uniform force field ,

Comparing the above result with Equation (16) shows that the generalized Einstein relation holds (the definition of generalized Einstein relation can be found in [37])

To connect to single-particle tracking experiments, we now turn to the time-averaged MSD of the stochastic process, defined by

The time series of length (the measurement time) is thus evaluated in terms of squared differences of the particle position separated by the so-called lag time , which defines the width of the window slid along the time series . Typically, is considered in the limit to obtain good statistics. It is easy to show that for Brownian motion, as long as the measurement is sufficiently long. Therefore, we call the process ergodic: ensemble averages and long-time averages are equivalent in the limit of long measurement times. signifies weak ergodicity breaking [38–40].

Corollary 4. The subordinated process governed by the multiterm time-fractional diffusion equation (1) is weak ergodicity breaking.

Proof. From Equation (16), we have

According to the definition of time-averaged MSD, we calculate the of .

The disparity between the ensemble and the -averaged MSD exists. Even in the limit of long measurement times , £¬, and therefore, the disparity still exists, which ends our proof.

4. Stochastic Simulation

The stochastic representation provides two ways to get the solution of the multiterm time-fractional advection-diffusion equation (1). One way is to get the analytical solution by substituting and into Equation (7). The other way is to simulate the stochastic representation, then use Monte Carlo to simulate the solution. The Monte Carlo method is firstly proposed to simulate the solution of fractional order equation [41]. Here, we mainly introduce how to simulate the sample path of the stochastic representation and get the simulated solution of the multiterm time-fractional diffusion equation.

The algorithm of simulation of the subordinated process is divided into two steps. is the horizon and .

Step 1. This step is aimed at simulating the subordinator (see Equation (4)). Since the is the strictly increasing -stable Levy motion with independent increments, then the process on the mesh () can be simulated as follows:

where is the i.i.d. strictly increasing -stable Levy noise [42, 43], generated by

Here, is uniform distribution on , and is exponential distribution with mean 1. Then, we can get the simulation of the process . From the definition of (4), we know the subordinator is the first passage time of . So, for every element , we only need to find the element such that , then . Since is a pure-jump process. For every jump of , there is a corresponding flat period of its inverse These heavy-tailed flat periods of represent long waiting times in which the subdiffusive particle gets immobilized in the trap. The sample path of can be found in Figure 1, and the corresponding waiting time in Figure 2. From the figures, we find the subordinator stop at the same time, which leads the waiting time to be fluctuant.

Step 2. This step is aimed at simulating the process . Since the parent process is driven by Brown motion, then we employ the Euler scheme to simulate the process , given by

where is the i.i.d. standard normal noise, . The sample path of can be found in Figure 3. Then, the Monte Carlo method can be employed to estimate the solution of Equation (1) (see Figure 4). From the figure, we find that the solution of Equation (1) has sharp peak and heavy tails, in contrast with normal distribution, which is called the stretched Gaussian distribution. The figure for is plotted to show these results more clearly (see Figure 4). Here, we remark that all the numerical results are obtained by the software Matlab.

Figure 1: The sample paths of the stochastic process , where and .

Figure 2: The corresponding waiting times with and .

Figure 3: The sample paths of the stochastic process , where , , , and .

Figure 4: The approximated solutions (up) of Equation (1) with different and are obtained for . Noting that we estimate 5000 times to get these solutions. Below is the corresponding Logarithmic solution.

5. Conclusions

In this paper, an advection-diffusion equation with multiterm time-fractional derivatives is employed. We obtained its stochastic representation, which is driven by the Brown motion and the inverse time of the sum of Lévy motions with different parameters. Then, the mean square displacement indicates the model is subdiffusive and the generalized Einstein relation is also retained, but weak ergodicity is breaking. At last, an algorithm is constructed to simulate the sample paths of the stochastic process. With the help of stochastic representation, the Monte Carlo method is employed to approximate the solution of the corresponding equation. We find that the solution is heavy tailed and sharp peaked, which is common in statistical physics and finance. So, we expect that the results obtained here may be useful for the discussion of the anomalous diffusion systems.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that no conflicts of interest exist regarding this manuscript.

Acknowledgments

We are grateful to Professor Zhongdi Cen for useful comments and suggestions. This work is supported by the National Natural Science Foundation of China No. 11801288, the Natural Science Foundation of Ningbo No. 2017A610130, and the Natural Science Foundation of Zhejiang Province No. LY17A010020.

References

M. Muskat, The Flow of Homogeneous Fluid through Porous Media, McGraw-Hill, New York, NY, USA, 1937.

P. I. Polubarinova-Koch, Theory of Ground Water Movement, Princeton University Press, 2015.

H. Spohn, “Surface dynamics below the roughening transition,” Journal de Physique, vol. 3, no. 1, pp. 69–81, 1993.

J. Buckmaster, “Viscous sheets advancing over dry beds,” Journal of Fluid Mechanics, vol. 81, no. 4, pp. 735–756, 1977.

P. Grosfils and J. P. Boon, “Nonextensive statistics in viscous fingering,” Physica A, vol. 362, no. 1, pp. 168–173, 2006.

S. S. Plotkin and P. G. Wolynes, “Non-Markovian configurational diffusion and reaction coordinates for protein folding,” Physical Review Letters, vol. 80, no. 22, pp. 5015–5018, 1998.

D. S. F. Crothers, D. Holland, Y. P. Kalmykov, and W. T. Coffey, “The role of Mittag-Leffler functions in anomalous relaxation,” Journal of Molecular Liquids, vol. 114, no. 1-3, pp. 27–34, 2004.

R. Metzler, E. Barkai, and J. Klafter, “Anomalous transport in disordered systems under the influence of external fields,” Physica A, vol. 266, no. 1-4, pp. 343–350, 1999.

D. Campos, V. Mendez, and J. Fort, “Description of diffusive and propagative behavior on fractals,” Physical Review E, vol. 69, no. 3, article 031115, 2004.

M. Wang and S. Chen, “Electroosmosis in homogeneously charged micro- and nanoscale random porous media,” Journal of Colloid and Interface Science, vol. 314, no. 1, pp. 264–273, 2007.

M. Wang, J. Wang, S. Chen, and N. Pan, “Electrokinetic pumping effects of charged porous media in microchannels using the lattice Poisson-Boltzmann method,” Journal of Colloid and Interface Science, vol. 304, no. 1, pp. 246–253, 2006.

P. J. Torvik and R. L. Bagley, “On the appearance of the fractional derivative in the Behavior of real materials,” Journal of Applied Mechanics, vol. 51, no. 2, pp. 294–298, 1984.

R. Schumer, D. A. Benson, M. M. Meerschaert, and B. Baeumer, “Fractal mobile/immobile solute transport,” Water Resources Research, vol. 39, no. 10, p. 1296, 2003.

F. A. Duck, Physicl Properties of Tissue: A Comperhensive Reference Book, Academic Press, Boston, MA, USA, 1990.

W. Chen and S. Holm, “Modified Szabo’s wave equation models for lossy media obeying frequency power law,” The Journal of the Acoustical Society of America, vol. 114, no. 5, pp. 2570–2574, 2003.

J. K. Kelly, R. J. McGough, and M. M. Meerschaert, “Analytical time-domain Green's functions for power-law media,” The Journal of the Acoustical Society of America, vol. 124, no. 5, pp. 2861–2872, 2008.

M. Liebler, S. Ginter, T. Dreyer, and R. E. Riedlinger, “Full wave modeling of therapeutic ultrasound: efficient time-domain implementation of the frequency power-law attenuation,” The Journal of the Acoustical Society of America, vol. 116, no. 5, pp. 2742–2750, 2004.

T. L. Szabo, “Time domain wave equations for lossy media obeying a frequency power law,” The Journal of the Acoustical Society of America, vol. 96, no. 1, pp. 491–500, 1994.

M. Caputo, “Linear models of dissipation whose Q is almost frequency independent–II,” Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967.

V. Daftardar-Gejji and S. Bhalekar, “Boundary value problems for multi-term fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 345, no. 2, pp. 754–765, 2008.

Y. Luchko, “Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation,” Journal of Mathematical Analysis and Applications, vol. 374, no. 2, pp. 538–548, 2011.

H. Jiang, F. Liu, I. Turner, and K. Burrage, “Analytical solutions for the multi-term time-space Caputo-Riesz fractional advection-diffusion equations on a finite domain,” Journal of Mathematical Analysis and Applications, vol. 389, no. 2, pp. 1117–1127, 2012.

R. Cheng, F. Sun, and J. Wang, “Meshless analysis of two-dimensional two-sided space-fractional wave equation based on improved moving least-squares approximation,” International Journal of Computer Mathematics, vol. 95, no. 3, pp. 540–560, 2017.

H. Ye, F. Liu, I. Turner, V. Anh, and K. Burrage, “Series expansion solutions for the multi-term time and space fractional partial differential equations in two- and three-dimensions,” The European Physical Journal Special Topics, vol. 222, no. 8, pp. 1901–1914, 2013.

A. H. Bhrawy and M. A. Zaky, “A method based on the Jacobi tau approximation for solving multi-term time-space fractional partial differential equations,” Journal of Computational Physics, vol. 281, pp. 876–895, 2015.

M. Zheng, F. Liu, V. Anh, and I. Turner, “A high-order spectral method for the multi-term time-fractional diffusion equations,” Applied Mathematical Modelling, vol. 40, no. 7-8, pp. 4970–4985, 2016.

F. Liu, M. M. Meerschaert, R. McGough, P. Zhuang, and Q. Liu, “Numerical methods for solving the multi-term time fractional wave-diffusion equations,” Fractional Calculus and Applied Analysis, vol. 16, no. 1, pp. 9–25, 2013.

E. A. Abdel-Rehim, A. M. A. El-Sayed, and A. S. Hashem, “Simulation of the approximate solutions of the time-fractional multi-term wave equations,” Computers Mathematics with Applications, vol. 73, no. 6, pp. 1134–1154, 2017.

V. Kolokoltsov, “On fully mixed and multidimensional extensions of the Caputo and Riemann-Liouville derivatives, related Markov processes and fractional differential equations,” Fractional Calculus and Applied Analysis, vol. 18, no. 4, pp. 1039–1073, 2015.

M. E. Hernandez-Hernandez, V. N. Kolokoltsov, and L. Toniazzi, “Generalised fractional evolution equations of Caputo type,” Chaos Solitons Fractals, vol. 102, pp. 184–196, 2017.

V. N. Kolokoltsov, “The probabilistic point of view on the generalized fractional partial differential equations,” Fractional Calculus and Applied Analysis, vol. 22, no. 3, pp. 543–600, 2019.

A. N. Kochubei and Y. Kondratiev, “Fractional kinetic hierarchies and intermittency,” Kinetic and Related Models, vol. 10, no. 3, pp. 725–740, 2017.

L. Lv and L. Wang, “Option pricing based on modified advection-dispersion equation: stochastic representation and applications,” Discrete Dynamics in Nature and Society, vol. 2020, Article ID 7168571, 8 pages, 2020.

B. P. Becker-Kern, M. M. Meerschaert, and H. Scheffler, “Limit theorems for coupled continuous time random walks,” The Annals of Probability, vol. 32, pp. 730–756, 2004.

H. Risken, The Fokker-Planck Equation, Springer, Berlin, 1996.

I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Academic Press, 1999.

R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,” Physics Reports, vol. 339, no. 1, pp. 1–77, 2000.

R. Metzler, J. H. Jeon, A. G. Cherstvy, and E. Barkai, “Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking,” Physical Chemistry Chemical Physics, vol. 16, no. 44, pp. 24128–24164, 2014.

L. Lv, F. Y. Ren, J. Wang, and J. Xiao, “Correlated continuous time random walk with time averaged waiting time,” Physica A: Statistical Mechanics and its Applications, vol. 422, pp. 101–106, 2015.

C. Merrigan, S. K. Birwa, S. Tewari, and B. Chakraborty, “Ergodicity breaking dynamics of arch collapse,” Physical Review E, vol. 97, no. 4, article 040901, 2018.

V. V. Uchaikin and V. V. Saenko, “Stochastic solution to partial differential equations of fractional orders,” Сuбupcκuй жypнал вычucлumeльнoй мameмamuκu, vol. 6, no. 2, pp. 197–203, 2003.

A. Janicki and A. Weron, Simulation and Chaotic Behaviour of α-Stable Stochastic Processes, Marcel Dekker, New York, NY, USA, 1994.

M. Magdziarz, “Stochastic representation of subdiffusion processes with time-dependent drift,” Stochastic Processes and their Applications, vol. 119, no. 10, pp. 3238–3252, 2009.