On the numerical solution of the diffusion equation with a nonlocal boundary condition.

Parabolic partial differential equations with nonlocal boundary specifications feature in the mathematical modeling of many phenomena. In this paper, numerical schemes are developed for obtaining approximate solutions to the initial boundary value problem for one-dimensional diffusion equation with a nonlocal constraint in place of one of the standard boundary conditions. The method of lines (MOL) semidiscretization approach is used to transform the model partial differential equation into a system of first-order linear ordinary differential equations (ODEs). The partial derivative with respect to the space variable is approximated by a second-order finite-difference approximation. The solution of the resulting system of first-order ODEs satisfies a recurrence relation which involves a matrix exponential function. Numerical techniques are developed by approximating the exponential matrix function in this recurrence relation. We use a partial fraction expansion to compute the matrix exponential function via Pade approximations, which is particularly useful in parallel processing. The algorithm is tested on a model problem from the literature.

[1]  Mehdi Dehghan,et al.  Fully explicit finite-difference methods for two-dimensional diffusion with an integral condition , 2002 .

[2]  John R. Cannon,et al.  A Galerkin procedure for the diffusion equation subject to the specification of mass , 1987 .

[3]  M. Dehghan.  Numerical solution of a non‐local boundary value problem with Neumann's boundary conditions , 2002 .

[4]  Karl Kunisch,et al.  A reaction-diffusion system arising in modelling man-environment diseases , 1988 .

[5]  A. Matheson,et al.  A numerical procedure for diffusion subject to the specification of mass , 1993 .

[6]  Gedeon Dagan,et al.  The significance of heterogeneity of evolving scales to transport in porous formations , 1994 .

[7]  John H. Cushman,et al.  Nonlocal Reactive Transport with Physical and Chemical Heterogeneity: Localization Errors , 1995 .

[8]  Abdul-Qayyum M. Khaliq,et al.  Time-stepping algorithms for semidiscretized linear parabolic PDEs based on rational approximants with distinct real poles , 1996, Adv. Comput. Math..

[9]  Gunnar Ekolin,et al.  Finite difference methods for a nonlocal boundary value problem for the heat equation , 1991 .

[10]  Lin Yanping,et al.  A numerical method for the diffusion equation with nonlocal boundary specifications , 1990 .

[11]  John van der Hoek,et al.  Diffusion subject to the specification of mass , 1986 .

[12]  R. Lazarov,et al.  Finite volume element approximations of nonlocal reactive flows in porous media , 2000 .

[13]  John R. Cannon,et al.  A class of non-linear non-classical parabolic equations , 1989 .

[14]  W. A. Day.  A decreasing property of solutions of parabolic equations with applications to thermoelasticity , 1983 .

[15]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[16]  E. H. Twizell,et al.  A family of fourth‐order parallel splitting methods for parabolic partial differential equations , 1997 .

[17]  Michael Renardy,et al.  Mathematical problems in viscoelasticity , 1987 .

[18]  Graeme Fairweather,et al.  The Reformulation and Numerical Solution of Certain Nonclassical Initial-Boundary Value Problems , 1991, SIAM J. Sci. Comput..

[19]  Avner Friedman,et al.  Monotonic decay of solutions of parabolic equations with nonlocal boundary conditions , 1986 .

[20]  A. R. Gourlay,et al.  The Extrapolation of First Order Methods for Parabolic Partial Differential Equations, II , 1978 .

[21]  J. G. Verwer,et al.  Solving parabolic integro-differential equations by an explicit integration method , 1992 .

[22]  Timothy R. Ginn,et al.  Nonlocal dispersion in media with continuously evolving scales of heterogeneity , 1993 .

[23]  Lin Yanping,et al.  An implicit finite difference scheme for the diffusion equation subject to mass specification , 1990 .

[24]  Jan ter Maten,et al.  Splitting methods for fourth order parabolic partial differential equations , 1986, Computing.