Skip to main content Accessibility help

Numerical Solution for a Non-Fickian Diffusion in a Periodic Potential

  • Adérito Araújo (a1), Amal K. Das (a2), Cidália Neves (a1) (a3) and Ercília Sousa (a1)


Numerical solutions of a non-Fickian diffusion equation belonging to a hyperbolic type are presented in one space dimension. The Brownian particle modelled by this diffusion equation is subjected to a symmetric periodic potential whose spatial shape can be varied by a single parameter. We consider a numerical method which consists of applying Laplace transform in time; we then obtain an elliptic diffusion equation which is discretized using a finite difference method. We analyze some aspects of the convergence of the method. Numerical results for particle density, flux and mean-square-displacement (covering both inertial and diffusive regimes) are presented.


Corresponding author



Hide All
[1]Abate, J., Whitt, W., Numerical inversion of Laplace transforms of probability distributions, ORSA Journal on Computing, 7(1) (1995), 3643.
[2]Ahn, J., Kang, S., Kwon, Y., A flexible inverse Laplace transform algorithm and its application, Computing 71(2) (2003), 115131.
[3]Babuska, I.M., Sauter, S.A., Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers?, SIAM Review 42(3) (2000), 451484.
[4]Bao, G., Wei, G.W., Zhao, S., Numerical solution of the Helmholtz equation with high wavenumbers, International Journal for Numerical Methods in Engineering 59 (2004), 389408.
[5]Barbero, G., Macdonald, J. R., Transport process of ions in insulating media in the hyperbolic diffusion regime, Phys. Rev. E 81(5) (2010), 051503.
[6]Branka, A.C., Das, A.K., Heyes, D.M., Overdamped Brownian motion in periodic symmetric potentials, J.Chem.Phys. 113(22) (2000), 99119919.
[7]Crump, K., Numerical inversion of Laplace transforms using a Fourier series approximation, Journal of the Association for Computing Machinery 23(1) (1976), 8996.
[8]Crank, J., Nicolson, P., A practical method for numerical evaluation of solutions of partial differential equations of the heat conduction type, Proc. Cambridge Phil. Soc. 43(1947), 5067.
[9]Das, A.K., A non-Fickian diffusion equation, J.Appl.Phys. 70(3) (1991), 13551358.
[10]Fibich, G., Ilan, B., Tsynkov, S., Backscattering and nonparaxiality arrest collapse of damped nonlinear waves, SIAM J. Appl. Math. 63(5) (2003), 17181736.
[11]Fort, J., Méndez, V., Wavefronts in time-delayed reaction-diffusion systems. Theory and comparison to experiment, Rep. Prog. Phys. 65(2002), 895954.
[12]Huang, al., Direct observation of the full transition from ballistic to diffusive Brownian motion in a liquid, Nature Physics, NPHYS1953, 2011.
[13]Li, T., Kheifets, S., Medellin, D., Raizen, M. G., Measurement of the instantaneous velocity of a Brownian particle, Science 25(328) (2010), 16731675.
[14]Magnasco, M.O., Forced thermal ratchets, Physical Review Letters 71 (1993), 14771481.
[15]Marsden, J.E., Hoffman, M.J., Basic Complex Analysis, W.H. Freeman, 1999.
[16]Neves, C., Araújo, A., Sousa, E., Numerical approximation of a transport equation with a time-dependent dispersion flux, in AIP Conference Proceedings 1048 (2008), 403406.
[17]Oliveira, F.S.B.F., Anastasiou, K., An efficient computational model for water wave propagation in coastal regions, Applied Ocean Research, 20(5) (1998), 263271.
[18]Pusey, P.N., Brownian motion goes ballistic, Science 332 (2011), 802803.
[19]Wang, K., You, C., A note on identifying generalized diagonally dominant matrices, International Journal of Computer Mathematics 84(12) (2007), 18631870.
[20]Varga, R.S., Matrix Iterative Analysis, Springer, 2000.


Related content

Powered by UNSILO

Numerical Solution for a Non-Fickian Diffusion in a Periodic Potential

  • Adérito Araújo (a1), Amal K. Das (a2), Cidália Neves (a1) (a3) and Ercília Sousa (a1)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.