arxiv: v3 [math-ph] 28 Dec 2016

Transcription

1 The time-fractional radiative transport equation Continuous-time random walk, diffusion approximation, and Legendre-polynomial expansion Manabu Machida a) (Dated: 3 December 1) arxiv:1.53v3 [math-ph] Dec 1 We consider the radiative transport equation in which the time derivative is replaced by the Caputo derivative. Such fractional-order derivatives are related to anomalous transport and anomalous diffusion. In this paper we describe how the time-fractional radiative transport equation is obtained from continuous-time random walk and see how the equation is related to the time-fractional diffusion equation in the asymptotic limit. Then we solve the equation with Legendre-polynomial expansion. PACS numbers: 5..Dd,5..Cd,7.5.+r a) Electronic mail: machida@hama-med.ac.jp; Institute for Medical Photonics Research, Hamamatsu University School of Medicine Hamamatsu, Shizuoka , Japan 1

2 I. INTRODUCTION Anomalous diffusion is often observed in nature 31,33. For example, tracer particles flowing in an aquifer exhibits anomalous diffusion 1. At the macroscopic scale after multiple scattering takes place, such anomalous diffusion is governed by fractional diffusion equations 31,3,39. Considering the fact that the diffusion equation appears in the asymptotic limit of the radiative transport equation or the linear Boltzmann equation 1, one can expect that at the mesoscopic scale there exist anomalous transport phenomena which are described by the fractional radiative transport equation. The use of the radiative transport equation was proposed for predicting the concentration of radionuclides in fractured rock underground,1. If this happens, then its fractional version must appear just like the fractional diffusion equation shows up when the diffusion process takes place in a complex structure. Letα (,1)andσ t,σ s (, )beconstants determined bythemedium under consideration. We suppose σ t > σ s. Let v > be a constant speed. Let u(x,µ,t) (x R, µ [,1], t [, )) be the angular density. We consider the following initial-value problem for the time-fractional radiative transport equation. t α u(x,µ,t)+vµ xu(x,µ,t)+σ t u(x,µ,t) = σ s p(µ,µ )u(x,µ,t)dµ, u(x,µ,) = δ(x)δ(µ µ ), where δ( ) is the Dirac delta function and α t is the Caputo fractional derivative 3, which is defined by 35 α t u(,,t) = 1 t t u(,,t ) Γ(1 α) (t t ) α dt, < α < 1, with Γ( ) the Gamma function. Indeed, u in (1) is the fundamental solution of the timefractional radiative transport equation. We note that recently α t was redefined more generally using fractional Sobolev spaces 1. Compared with the Riemann-Liouville derivative, the Caputo derivative is not singular at t =. Thus we can have the same initial condition in (1) and in the corresponding equation of the first derivative t. The phase function p(µ,µ ) satisfies p(µ,µ )dµ = 1, µ [,1]. Anomalous transport phenomena are in the transport regime when the distance of interest is not too large compared to the scattering mean free path v/σ s, and as is shown (1)

3 below, the time-fractional diffusion equation is obtained from (1) in the asymptotic limit. The time-fractional diffusion equation has been intensively studied. In addition to several examples 31,33, we point out that the behavior of water transport in granite was successfully reproduced by the random walk process with a power-law distribution 11. It is proposed that if there are two porosities, the mass transport in fractured porous aquifer should be governed by the diffusion equation in which both t and t α appear. The Cauchy problem and initial-boundary-value problem, were considered for the time-fractional diffusion equation. The maximum principle was established 3. The technique of eigenfunction expansion was developed 37. Numerical algorithms for the equation have been developed. Moreover the standard time-fractional diffusion equation was generalized to equations with multiple Caputo derivatives 19,5 and distributed-order equations 1,. See the recent review by Jin and Rundell 1. The rest of the paper is organized as follows. In II, we obtain the time-fractional radiative transport equation from continuous-time random walk. In III, we see that the timefractional diffusion equation emerges from the time-fractional radiative transport equation when absorption is small, propagation distance is large, and observation time is long. In IV, we express the solution to the time-fractional radiative transport equation in the form of Legendre polynomial expansion. In V, we numerically compute the solutions of the timefractional radiative transport equation and of the time-fractional diffusion equation. Finally in VI, concluding remarks are made. The subtraction of the ballistic term is considered in Appendix. II. CONTINUOUS-TIME RANDOM WALK We consider the continuous-time random walk whose jump probability density function ϕ(x,t;µ,µ ) (x R, t [, ), µ,µ [,1]) is given by ϕ(x,t;µ,µ ) = [ξ s δ(x)p(µ,µ )+(1 ξ t )δ(x µr)δ(µ µ )]w(t), () where ξ t (,1), ξ s (,ξ t ), and r > are some constants. The first term represents scattering and the second term in the square brackets of () is responsible for transport. The waiting time probability density function w(t) is obtained as (1 ξ a )w(t) = 3 ϕ(x,t;µ,µ )dxdµ,

4 where ξ a = ξ t ξ s > is the probability for absorption. The left-hand side of the abovementioned equation shows the probability that the test particle is not absorbed in the medium and makes a jump after the time t. Let η(x, µ, t) be the probability density function of just having arrived at position x at time t in direction µ. Let P(x,µ,t) be the probability density function of being at (x,µ,t) R [,1] [, ). We consider the following continuous-time random walk process. η(x,µ,t) = t η(x,µ,t )ϕ(x x,t t ;µ,µ )dx dµ dt +a(x,µ)δ(t), P(x,µ,t) = t η(x,µ,t )Φ(t t )dt, where a(x,µ) is the initial value which is a function of x and µ, Φ(t) is the cumulative probability of not having moved during t, which is given by where Φ(t) = 1 By the Fourier-Laplace transform we have Hence we obtain (LFη)(k,µ,s) = (LFP)(k,µ,s) = [ t w(t )dt. e st e ikx P(x,µ,t)dxdt = (LFη)(k,µ,s)(LΦ)(s), (LΦ)(s) = 1 (Lw)(s). s ξ s p(µ,µ )(LFη)(k,µ,s)dµ + (1 ξ t )(LFη)(k,µ,s)e iµrk ](Lw)(s)+(Fa)(k,µ). We consider small k and use e iµrk 1 iµrk. Thus we arrive at 1 (Lw)(s) [(LP)(x,µ,s) 1s ] (Lw)(s) P(x,µ,) = ξ s p(µ,µ )(LP)(x,µ,s)dµ [ξ t +(1 ξ t )rµ x ](LP)(x,µ,s).

5 Recalling < α < 1, we have 35,3 (L α t f)(s) = sα (Lf)(s) s α f(). Let us assume that the waiting time probability density function behaves as (Lw)(s) 1 (τs) α, < s 1 τ. We introduce σ t = ξ t τ α, σ s = ξ s τ α, v = (1 ξ t)r τ α. We asymptotically obtain This is (1). t α P(x,µ,t)+vµ xp(x,µ,t)+σ t P(x,µ,t) = σ s p(µ,µ )P(x,µ,t)dµ. Remark II.1. In this section we implemented the effect of absorption in our random walk by introducing ξ a. Such extension of the usual continuous-time random walk is done by Hornung, Berkowitz, and Barkai 13, and by Henry, Langlands, and Wearne 1. Indeed, we arrive at the same conclusion by instead writing () as ϕ(x,t;µ,µ ) = [ξ s δ(x)p(µ,µ )+(1 ξ t )δ(x µr)δ(µ µ )] with the waiting time probability density function w(t) introduced as w(t) = We can then give η(x,µ,t) and P(x,µ,t) as η(x,µ,t) = (1 ξ a ) t ϕ(x,t;µ,µ )dxdµ. w(t) 1 ξ a, η(x,µ,t )ϕ(x x,t t ;µ,µ )dx dµ dt +a(x,µ)δ(t), P(x,µ,t) = (1 ξ a ) t η(x,µ,t )Φ(t t )dt. Note that P(x,µ,) = (1 ξ a )a(x,µ). Thus the relation to the past work 1,13 becomes clearer. 5

6 III. DIFFUSION APPROXIMATION Let us suppose that the ratio ǫ > of the mean free path to the propagation distance is small. We scale t,x as t ǫ /α t and x ǫx. Furthermore we scale σ a σ a /ǫ assuming σ a is small (recall σ a = σ t σ s ). Although the radiative transport equation (1) has the Caputo derivative, we obtain the time-fractional diffusion equation by following the standard procedure,1,3. In this section we assume that p(µ,µ ) = p(µ,µ). We can write the time-fractional radiative transport equation as ǫ t α u(x,µ,t)+ǫvµ xu(x,µ,t)+ ( ) 1 ǫ σ a +σ s u(x,µ,t) = σs p(µ,µ )u(x,µ,t)dµ. We write u(x,µ,t) = U DA (x,µ,t)+ǫu (1) DA (x,µ,t)+ǫ U () DA (x,µ,t)+. Let us collect terms of order ǫ. We obtain σ s U DA (x,µ,t) = σ s p(µ,µ )U DA (x,µ,t)dµ. The above equation implies that U DA is independent of µ; hereafter we write U DA (x,µ,t) = U DA (x,t). The terms of order ǫ 1 yields We obtain where g (,1) satisfies vµ x U DA (x,t)+σ s U (1) DA (x,µ,t) = σ s U (1) DA (x,µ,t) = v µ x U DA (x,t), (1 g)σ s µg = By collecting terms of order ǫ we have µ p(µ,µ )dµ. p(µ,µ )U (1) DA (x,µ,t)dµ. t α U DA(x,t)+µ x U (1) DA (x,µ,t)+σ au DA (x,t)+σ s U () DA (x,µ,t) 1 = σ s p(µ,µ )U () DA (x,µ,t)dµ. If we integrate the above equation over µ, we obtain α t U DA(x,t) D x U DA(x,t)+σ a U DA (x,t) =, (3) where D = v 3(1 g)σ s. ()

7 Thus the time-fractional diffusion equation is obtained in the asymptotic limit of (1). One remark needs to be made. We have the second derivative for the spatial variable x in (3). In a similar setting, it is known that the space-fractional diffusion equation is obtained if the phase function decays with power-law as a function of the speed of propagating particles 9,3. IV. LEGENDRE-POLYNOMIAL EXPANSION Let us suppose p(µ,µ ) is given by p(µ,µ ) = 1 L β l P l (µ)p l (µ ), l= where L, and β l (l =,1,...,L) are positive constants such as β = 1, β l < l +1 for l 1. Here, P l (µ) are the Legendre polynomials recursively given by (l +1)P l+1 (µ) = (l+1)µp l (µ) lp l (µ), P 1 (µ) = µ, P (µ) = 1, µ [,1]. In the time-independent case, an analytical solution of the space-fractional radiative transport equation was found 15. In this section we solve (1). Let us expand u with Legendre polynomials. (Fu)(k,µ,t) = l+1cl (k,t;µ )P l (µ). (5) We perform the Fourier transform in (1) and substitute (5). We have ( t α +ivkµ+σ t ) l+1cl (k,t;µ )P l (µ) Let us introduce l= l= β l = σ s l+1cl (k,t;µ ) l+1 P l(µ)θ(l l). l= h l = l+1 σ s σ t β l Θ(L l). Let N ( L) be an integer. We take projections with P l (µ) (l =,1,...,N) and obtain ivkl l c l + ivk(l+1) (l+1) c l+1 + t α c l + σ th l l+1 c l =, where we used the recurrence relations and orthogonality relations of Legendre polynomials, µp l (µ) = l+1 l+1 P l+1(µ)+ l l+1 P l(µ), () 7

8 and P l (µ)p l (µ)dµ = l +1 δ ll. The above equation is expressed as A(k)c(k,t;µ )+ α t c(k,t;µ ) =, where A(k) is an (N +1) (N +1) matrix and c(k,t;µ ) is an N +1 dimensional vector defined by {A(k)} ll = ivkl l δ l,l + σ th l l+1 δ ivk(l +1) l,l + (l+1) δ l+1,l, (7) {c(k,t;µ )} l = c l (k,t;µ ). () When the Legendre polynomial expansion is used, tridiagonal matrices such as A(k) appear due to the three-term recurrence relation () 7,,1,3. By taking the Laplace transform we have where we used (Lc)(k,s;µ ) = (A(k)+s α ) s α c(k,;µ ), (L α t c)(k,s;µ ) = s α (Lc)(k,s;µ ) s α c(k,;µ ), < α 1. Let us recall that the Mittag-Leffler function is given by 35 E α (z) := n= and the Laplace transform is obtained as z n, z,α C, Rα >, Γ(αn+1) L{E α (zt α ); s} = sα s α z, z,s,α C, Rs,Rα >, z s α < 1. Thus we find c(k,t;µ ) = E α ( A(k)t α )c(k,;µ ). Since δ(µ µ ) = l+1 l= P l(µ)p l (µ ), we obtain {c(k,;µ )} l = l+1 P l (µ ). Let λ n (k) and v n (k) be the nth eigenvalue and eigenvector of the matrix A(k). We can write A(k) as A(k) = Q(k)D(k)Q(k),

9 where Q(k) = (v (k) v 1 (k) v N (k)), D(k) = diag(λ (k),λ 1 (k),...,λ N (k)). We have {A(k)} ij = { Q(k)D(k)Q(k) } ij = N n= λ n (k)v (i) n (k)v(j) n (k), where v (i) n (k) is the ith component of v n (k). Therefore we can write {c(k,t;µ )} l = N j +1 P j (µ ) j= N n= v (l) n (k)v(j) n (k)e α ( λ n (k)t α ). Noting (), Eq. (5) yields u(x,µ,t) u(x,µ,t;n) := 1 π Since k appears always as ik, we see We obtain u(x,µ,t;n) = N e ikx l+1cl (k,t;µ )P l (µ)dk. (9) l= c l ( k,t;µ ) = c l (k,t;µ ). N l+1 P l (µ) π l= [cos(kx)rc l (k,t;µ ) sin(kx)ic l (k,t;µ )] dk. (1) Remark IV.1. Although in this section we directly calculated u in (1), indeed, it is possible to directly relate u(x,µ,t) to u 1 (x,µ,t) which is the solution of (1) with α = 1. Let f α (t) be a function such that (Lf α )(s) = e sα. For example, we have If we introduce we have f 1/ (t) = t 3/ π e/(t). ϕ(τ,t) = ( ) t t ατ 1+1/αf α, τ 1/α (Lϕ)(τ,s) = s α e τsα. 9

10 Let us consider the Laplace transform of u with respect to s and u 1 with respect to s α. Assuming u(x,µ,) = u 1 (x,µ,), we obtain s α (Lu)(x,µ,s) s α u(x,µ,)+vµ x (Lu)(x,µ,s)+σ t (Lu)(x,µ,x) The above equations imply = σ s p(µ,µ )(Lu)(x,µ,s)dµ, s α (Lu 1 )(x,µ,s α ) u 1 (x,µ,)+vµ x (Lu 1 )(x,µ,s α )+σ t (Lu 1 )(x,µ,s α ) (Lu)(x,µ,s) = s α (Lu 1 )(x,µ,s α ) = Therefore u and u 1 are related as u(x,µ,t) = = σ s p(µ,µ )(Lu 1 )(x,µ,s α )dµ. u 1 (x,µ,τ)ϕ(τ,t)dτ. u 1 (x,µ,t)s α e τsα dτ. This means that we can obtain u by integrating u 1, which is the solution of the first-order equation. The solution u is subordinated to the solution u V. NUMERICAL CALCULATION The energy density U(x, t) is introduced as U(x,t) = Each N gives an approximated value of U(x,t) as u(x,µ,t)dµ. U(x,t) U(x,t;N), where U(x,t;N) = u(x,µ,t;n)dµ. We note that U(x,t) = U(x,t; ). Let us calculate U(x,t;N) for the initial condition U(x,;N) = δ(x). 1

11 From (1) we obtain U(x,t;N) = = 1 π N e ikx n= N In this section we set v () v () n (k) u(x,µ,t;n)dµdµ n (k) E α ( λ n (k)t α ) dk = π n= ( cos(kx)re α ( λ n (k)t α ) sin(kx)ie α ( λ n (k)t α ) ) dk. v = 1, σ a =, L = N = 1, and The matrix A(k) in (7) is given by where we introduced σ s = 1, g = β 1 3 =.9. A(k) = 1 ik 3 ik k c k c := 3 σ s(1 g). Its eigenvalues and eigenvectors are obtained as λ(k) = k ( ) c k 1± 1, 3 k c and Thus we have. ( k v(k) = 1 ik c k 1± 1 3, N = 3 k c N λ(k) v () (k) = ) ( ) 1 k 1 1, k k c, k c 11 1, k > k c., k k c, 3 k, k > k c.

12 The energy density is written as kc U(x,t;1) = 1 cos(kx) [ π ( k π ( k 1 k c k c k c cos(kx)re α ) E α ( ) E α ( ( k c + ) kc k t α 3 k c k c k 3 t α )]dk k c i k kc t )dk. α (11) 3 In the diffusion approximation the energy density is given as follows. If the initial condition is given by U DA (x,) = δ(x), we have 7, U DA (x,t) = 1 cos(kx)e α ( D k t α )dk (1) π = 1 ( ) x t α Mα/, D D t α/ where M α (z) is the M-Wright function defined by () n z n M α (z) := n!γ( α(n+1)+1). n= Equations (11) and (1) are implemented in Fortran. The numerical implementation of the Mittag-Lifter function relies on the algorithm by Gorenflo, Loutchko, and Luchko 9. Although we saw in III that U(x,t) asymptotically becomes U DA (x,t), they are different in general. In Figs. 1 through 3, we plot U(x,t;1) and U DA (x,t) for α =.5,.5, and.75, respectively. For all the cases, we see that U(x,t;1) stays near the source at x = for a relatively long time whereas U DA (x,t) broadens quickly. When α =.75 we can see that U(x,t;1) has two peaks. Such a double-peak structure shows up for α > 1 in the case of the fractional diffusion equation. This behavior can be understood from the relation 5 E α (z)+e α ( z) = E α (z ), z C. For sufficiently large k, which corresponds to small x, we asymptotically have 35 ( E α k c i ) ( k kc t α 1 3 α exp k c i ) 1/α k kc t α. 3 1

13 Hence in (11) we have RE α ( k c i ) [ ( k kc t α 1 3 α exp i k ) ] [ 1/α ( 3 tα + 1 α exp i k 1 E α (i 3 k ) t α + 1 E α ( i 3 k ) t α ( = E α 1 ) 3 k t α. 3 tα ) 1/α ] The above calculation implies that the double-peak behavior for the fractional diffusion equationwithα > 1canbeseenforthefractionalradiativetransportequationwithα > 1/ x U U DA x U U DA x U U DA FIG. 1. Comparison of U(x,t) and U DA (x,t) as a function of x, from the left, for t =.1,.5, and.1, respectively when α = x U U DA x U U DA x U U DA FIG.. Comparison of U(x,t) and U DA (x,t) as a function of x, from the left, for t =.1,.5, and.1, respectively when α =.5. 13

14 x U U DA x U U DA x U U DA FIG. 3. Comparison of U(x,t) and U DA (x,t) as a function of x, from the left, for t =.5,.1, and., respectively when α =.75. VI. CONCLUDING REMARKS One of the purposes of the present paper is to see the connection between the timefractional radiative transport equation and the time-fractional diffusion equation. Roughly speaking, the time-fractional radiative transport equation of α t behaves as the timefractional diffusion equation of α t for large x and behaves as the time-fractional diffusion equation of α t near x = as is investigated in III and V. When u(x,µ,t)in(1) is expressed inthe formof the collision expansion, the ballistic term is singular. If u(x, µ, t) itself is numerically computed, it is desirable to subtract the ballistic term. In a straightforward manner, we can extend the calculation in IV. This calculation is summarized in Appendix. Appendix A: Subtraction of the ballistic term Let us split u(x,µ,t) in (1) into the ballistic and scattered parts as u(x,µ,t) = u b (x,µ,t)+u s (x,µ,t), where u b (x,µ,t) and u s (x,µ,t) respectively satisfy t α u b (x,µ,t)+µ x u b (x,µ,t)+σ t u b (x,µ,t) =, u b (x,µ,) = δ(x)δ(µ µ ), 1

15 and t α u s(x,µ,t)+µ x u s (x,µ,t)+σ t u s (x,µ,t) = σ s p(µ,µ )u s (x,µ,t)dµ +S(x,µ,t), u s (x,µ,) =. Here the source term for u s (x,µ,t) is given by Noting that we obtain and S(x,µ,t;µ ) = σ s p(µ,µ )u b (x,µ,t)dµ. (LFu b )(k,µ,s) = s α s α +ikµ+σ t δ(µ µ ), u b (x,µ,t) = 1 π δ(µ µ ) e ikx E α [ (ikµ +σ t )t α ] dk, s α (LFS)(k,µ,s;µ ) = σ s p(µ,µ ). s α +ikµ +σ t Let us expand u s with Legendre polynomials. For l N we obtain (Fu s )(k,µ,t) = l +1cl (k,t;µ )P l (µ). l= (A1) A(k)c(k,t;µ )+ α t c(k,t;µ ) = w(k,t;µ ), where w(k,t;µ ) is an N +1 dimensional vector defined by {w(k,t;µ )} l = l+1 P l (µ)(fs)(k,µ,t;µ )dµ. By taking the Laplace transform we have (Lc)(k,s;µ ) = (A(k)+s α ) [ s α c(k,;µ )+(Lw)(k,s;µ ) ]. Let us express the Laplace transform of w(k,t;µ ) as (Lw)(k,s;µ ) = s α s α +ikµ +σ t b(µ ), 15

16 where Using the relation {b(µ )} l = σ sβ l l+1 Θ(L l)p l(µ ). we find = (A(k)+ikµ +σ t ) ( (A(k) s α ) (Lw)(k,s;µ ) ) s α sα b(µ s α +ikµ +σ t s α ), A(k) c(k,t;µ ) = E α ( A(k)t α )c(k,;µ ) +(A(k)+ikµ +σ t ) [E α ( (ikµ +σ t )t α ) E α (A(k)t α )]b(µ ). REFERENCES 1 Adams, E. E. and Gelhar, L. W., Field study of dispersion in a heterogeneous aquifer. Spatial moments analysis, Water Res. Res., (199). Arridge, S. R. and Schotland, J. C., Optical tomography: forward and inverse problems, Inverse Problems 5, 131 (9). 3 Caputo, M., Linear model of dissipation whose Q is almost frequency independent-ii, Geophys. J. R. Astr. Soc. 13, (197). Eidelman, S. D. and Kochubei, A. N., Cauchy problem for fractional diffusion equations, J. Diff. Eq. 199, (). 5 Erdélyi, A., Magnus, W., Oberhettinger, F., and Tricomi, F. G., Higher Transcendental Functions Vol. 3 (McGraw-Hill, 1955). Fomin, S.A., Chugunov, V.A., andhashida, T., Non-Fickianmasstransportinfractured porous media, Adv. Water Resour. 3, 5 1 (11). 7 Garcia, R. D. M. and Siewert C. E., On discrete spectrum calculations in radiative transfer, J. Quant. Spec. Rad. Trans., (199). Gershenson, M., Time-dependent equation for the intensity in the diffusion limit using a higher-order angular expansion, Phys. Rev. E 59, (1999). 9 Gorenflo, R., Loutchko, J., and Luchko, Y., Computation of the Mittag-Leffler function E α,β (z) and its derivative, Fract. Calc. Appl. Anal. 5, (). 1 Gorenflo, R., Luchko, Y., and Yamamoto, M., Time-fractional diffusion equation in the fractional Sobolev spaces, Fract. Calc. Appl. Anal. 1, 799 (15). 1

17 11 Hatano, Y. and Hatano, N., Dispersive transport of ions in column experiments: An explanation of long-tailed profiles, Water Resour. Res. 3, (199). 1 Henry, B. I., Langlands, T. A. M., and Wearne, S. L., Anomalous diffusion with linear reaction dynamics: From continuous time random walks to fractional reaction-diffusion equations, Phys. Rev. E 7, 3111 () 13 Hornung, G., Berkowitz, B., and Barkai, N., Morphogen gradient formation in a complex environment: An anomalous diffusion model, Phys. Rev. E 7, 191 (5). 1 Jin, B. and Rundell, W., A tutorial on inverse problems for anomalous diffusion processes, Inverse Problems 31, 353 (15). 15 Kadem, A., Luchko, Y., and Baleanu, D., Spectral method for solution of the fractional transport equation, Rep. Math. Phys., (1). 1 Kochubei, A. N., Distributed order calculus and equations of ultraslow diffusion, J. Math. Anal. Appl. 3, 5 1 (). 17 Langlands, T. A. M., Henry, B. I., and Wearne, S. L., Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions, J. Math. Biol. 59, 71 (9). 1 Larsen, E. W. and Keller, J. B., Asymptotic solution of neutron transport problems for small mean free paths, J. Math. Phys. 15, 75 1 (197). 19 Li, Z., Liu, Y., and Yamamoto, M., Initial-boundary value problems for multi-term timefractional diffusion equations with positive constant coefficients, Appl. Math. Comp. 57, (15). Li, Z., Luchko, Y., and Yamamoto, M., Asymptotic estimates of solutions to initialboundary-value problems for distributed order time-fractional diffusion equations, Fract. Cal. Appl. Anal. 17, (1). 1 Liemert, A. and Kienle, A., Infinite space Green s function of the time-dependent radiative transfer equation, Biomed. Opt. Exp. 3, (1). Lin, Y. and Xu, C., Finite difference/spectral approximations for the time-fractional diffusion equation, J. Comp. Phys. 5, (7). 3 Luchko, Y., Maximum principle for the generalized time-fractional diffusion equation, J. Math. Anal. Appl. 351, 1 3 (9). Luchko, Y., Some uniqueness and existence results for the initial-boundary-value problems for the generalized time-fractional diffusion equation, Comp. Math. Appl. 59, 17 17

18 177 (1). 5 Luchko, Y., Initial-boundary-value problems for the generalized multi-term timefractional diffusion equation, J. Math. Anal. Appl. 37, 53 5 (11). Luchko, Y., Initial-boundary-value problems for the one-dimensional time-fractional diffusion equation, Fract. Cal. Appl. Anal. 15, 11 1 (1). 7 Mainardi, F., Thefundamentalsolutionsforthefractionaldiffusion-waveequation, Appl. Math. Lett. 9, 3 (199). Mainardi, F., Luchko, Y., Pagnini, G., The fundamental solution of the space-time fractional diffusion equation, Fract. Cal. Appl. Anal., (1). 9 Mellet, A., Fractional diffusion limit for collisional kinetic equations: A moments method, Indiana Univ. Math. J. 59, (1). 3 Mellet, A., Mischler, S., and Mouhot, C., Fractional diffusion limit for collisional kinetic equations, Arch. Rational Mech. Anal. 199, (11). 31 Metzler, R. and Klafter, J., The random walk s guide to anomalous diffusion: a fractional dynamics approach, Phys. Rep. 339, 1 77 (). 3 Metzler, R. and Klafter, J., The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics, J. Phys. A: Math. Gen. 37, R11 R (). 33 Metzler, R., Jeon, J.-H., Cherstvya, A. G., and Barkaid, E., Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking, Phys. Chem. Chem. Phys. 1, 1 1 (1). 3 Panasyuk, G., Schotland, J. C., and Markel, V. A., Radiative transport equation in rotated reference frames, J. Phys. A: Math. Gen. 39, (). 35 Podlubny, I., Fractional Differential Equations (Academic Press, 1999). 3 Ryzhik, L., Papanicolaou, G., and Keller, J. B., Transport equations for elastic and other waves in random media, Wave Motion, (199). 37 Sakamoto, K. and Yamamoto, M., Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems, J. Math. Anal. Appl. 3, 7 (11). 3 Samko, S. G., Kilbas, A. A., and Marichev, O. I., Fractional integrals and derivatives: theory and applications (Gordon and Breach Science, 1993). 39 Sokolov, I., Klafter, J., and Blumen, A., Fractional Kinetics, Physics Today 55, 5 1

19 (). Williams, M. M. R., Stochastic problems in the transport of radioactive nuclides in fractured rock, Nucl. Sci. Eng. 11, 15 3 (199). 1 Williams, M. M. R., Radionuclide transport in fractured rock a new model: application and discussion, Ann. Nucl. Energy, (1993). 19