Numerical solution of the space fractional Fokker-Planck equation

The traditional second-order Fokker-Planck equation may not adequately describe the movement of solute in an aquifer because of large deviation from the dynamics of Brownian motion. Densities of α-stable type have been used to describe the probability distribution of these motions. The resulting governing equation of these motions is similar to the traditional Fokker-Planck equation except that the order α of the highest derivative is fractional.In this paper, a space fractional Fokker-Planck equation (SFFPE) with instantaneous source is considered. A numerical scheme for solving SFFPE is presented. Using the Riemann-Liouville and Grunwald-Letnikov definitions of fractional derivatives, the SFFPE is transformed into a system of ordinary differential equations (ODE). Then the ODE system is solved by a method of lines. Numerical results for SFFPE with a constant diffusion coefficient are evaluated for comparison with the known analytical solution. The numerical approximation of SFFPE with a time-dependent diffusion coefficient is also used to simulate Levy motion with α-stable densities. We will show that the numerical method of SFFPE is able to more accurately model these heavy-tailed motions.

[1]  Massimiliano Giona,et al.  A theory of transport phenomena in disordered systems , 1992 .

[2]  W. Woyczynski,et al.  Fractal Burgers Equations , 1998 .

[3]  Rudolf Hilfer,et al.  EXACT SOLUTIONS FOR A CLASS OF FRACTAL TIME RANDOM WALKS , 1995 .

[4]  S. Bhatia,et al.  Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions , 2001 .

[5]  D. Benson,et al.  The fractional‐order governing equation of Lévy Motion , 2000 .

[6]  S. Bhatia,et al.  Effect of pore blockage on adsorption isotherms and dynamics: anomalous adsorption of iodine on activated carbon , 2000 .

[7]  V. Anh,et al.  Fractional diffusion and fractional heat equation , 2000, Advances in Applied Probability.

[8]  J. Bouchaud,et al.  Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications , 1990 .

[9]  W. Wyss,et al.  THE FRACTIONAL BLACK-SCHOLES EQUATION , 2000 .

[10]  C. W. Gear,et al.  The automatic integration of ordinary differential equations , 1971, Commun. ACM.

[11]  W. Wyss.  The fractional diffusion equation , 1986 .

[12]  Socio Michele Caputo.  The Green function of the diffusion of fluids in porous media with memory , 1996 .

[13]  Fawang Liu,et al.  An unstructured mesh finite volume method for modelling saltwater intrusion into coastal aquifers , 2002 .

[14]  Roberto Monaco,et al.  Waves and Stability in Continuous Media , 1991 .

[15]  John H. Cushman,et al.  Dynamics of fluids in hierarchical porous media , 1990 .

[16]  N. Leonenko,et al.  Spectral Analysis of Fractional Kinetic Equations with Random Data , 2001 .

[17]  Numerical simulation of the power density distribution generated in a multimode cavity by using the method of lines technique to solve directly for the electric field , 1996 .

[18]  M. Celia,et al.  Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts, 2, Analysis of spatial moments for a nonreactive tracer , 1991 .

[19]  K. Miller,et al.  An Introduction to the Fractional Calculus and Fractional Differential Equations , 1993 .

[20]  D. Benson,et al.  Application of a fractional advection‐dispersion equation , 2000 .

[21]  S. Bhatia,et al.  Application of Petrov–Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids , 2001 .

[22]  周建和,et al.  Anomalous Diffusion in Disordered Media , 2000 .

[23]  S. Bhatia,et al.  Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles , 1999 .

[24]  Cass T. Miller,et al.  Accurate and economical solution of the pressure-head form of Richards' equation by the method of lines , 1997 .

[25]  O. Marichev,et al.  Fractional Integrals and Derivatives: Theory and Applications , 1993 .

[26]  Vo Anh,et al.  Scaling laws for fractional diffusion-wave equations with singular data , 2000 .

[27]  I. Podlubny.  Fractional differential equations , 1998 .

[28]  W. Schneider,et al.  Fractional diffusion and wave equations , 1989 .

[29]  R. Gorenflo,et al.  Wright functions as scale-invariant solutions of the diffusion-wave equation , 2000 .

[30]  A two-dimensional finite volume method for transient simulation of time- and scale-dependent transport in heterogeneous aquifer systems , 2003 .

[31]  F. R. Norwood,et al.  IUTAM symposium : nonlinear waves in solids , 1995 .

[32]  J. Klafter,et al.  The random walk's guide to anomalous diffusion: a fractional dynamics approach , 2000 .

[33]  W. W. Wood,et al.  Large-Scale Natural Gradient Tracer Test in Sand and Gravel, , 1991 .

[34]  H. Risken.  Fokker-Planck Equation , 1984 .

[35]  S. Bhatia,et al.  Numerical solution of hyperbolic models of transport in bidisperse solids , 2000 .

[36]  I. Turner,et al.  Time fractional advection-dispersion equation , 2003 .

[37]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.