×

Numerical solution of a parabolic equation subject to specification of energy. (English) Zbl 1038.65088

Summary: Special attention is devoted to one particular problem, the diffusion equation subject to the specification of mass in a portion of the domain, which has been studied quite extensively, both analytically and numerically, in recent years. Parabolic partial differential equations (PDEs) with a non-local constraint in place of one of the standard boundary specifications feature in the mathematical modeling of many phenomena.
In this paper the application of the method of lines (MOL) to such problems is considered. The MOL semi-discretization approach will be used to transform the model partial differential equation into a system of first-order linear ordinary differential equations (ODEs). The MOL is a method of solving PDEs by discretizing the equation with respect to all but one variable (usually time). The spatial partial derivative is approximated by a finite-difference method.
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. The matrix exponential function is approximated by a rational approximation consisting of four parameters. New parallel algorithms are developed using the resulting approximation. Numerical experiments on two challenging examples are presented to illustrate the performance of the algorithms.

MSC:

65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
35K05 Heat equation
65L12 Finite difference and finite volume methods for ordinary differential equations
65Y05 Parallel numerical computation
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Allegretto, W.; Lin, Y.; Zhou, A., A box scheme for coupled systems resulting from micro-sensor thermistor problems, Dyn. Discret. Contin. Impuls Syst., 5, 209-223 (1999) · Zbl 0979.78023
[2] Cannon, J. R.; van der Hoek, J., Diffusion subject to specification of mass, J. Math. Anal. Appl., 115, 517-529 (1986) · Zbl 0602.35048
[3] Cannon, J. R.; Matheson, A. L., A numerical procedure for diffusion subject to the specification of mass, Int. J. Engng. Sci., 31, 3, 347-355 (1993) · Zbl 0773.65069
[4] Cannon, J. R.; Prez-esteva, S.; van der Hoek, J., A Galerkin procedure for the diffusion equation subject to the specification of mass SIAM, J. Numer. Anal., 24, 499-515 (1987) · Zbl 0677.65108
[5] Cannon, J. R.; van der Hoek, J., Implicit finite difference scheme for the diffusion of mass in porous media, (Noye, J., Numerical Solutions of Partial Differential Equations (1982), North Holland), 527-539
[6] Cannon, J. R.; Yin, H. M., On a class of non-classical parabolic problems, J. Differen. Equat., 79, 266-288 (1989) · Zbl 0702.35120
[7] Cannon, J. R.; Lin, Y.; Wang, S., An implicit finite difference scheme for the diffusion equation subject to the specification of mass, Int. J. Engng. Sci., 28, 7, 573-578 (1990) · Zbl 0721.65046
[8] Capasso, V.; Kunisch, K., A reaction-diffusion system arising in modeling man-environment diseases, Quart. Appl. Math., 46, 431-449 (1988) · Zbl 0704.35069
[9] Cushman, J. H.; Xu, H.; Deng, F., Nonlocal reactive transport with physical and chemical heterogeneity: Localization error, Water Resour. Res., 31, 2219-2237 (1995)
[10] Cushman, J. H.; Ginn, T. R., Nonlocal dispersion in media with continuously evolving scales of heterogeneity, Transp. Porous Med., 13, 123-138 (1993)
[11] Dagan, G., The significance of heterogeneity of evolving scales to transport in porous formations, Water Resour. Res., 13, 3327-3336 (1994)
[12] Day, W. A., Existence of a property of solutions of the heat equation subject to linear thermoelasticity and other theories, Quart. Appl. Math., 40, 319-330 (1982) · Zbl 0502.73007
[13] Day, W. A., A decreasing property of solutions of a parabolic equation with applications to thermoelasticity and other theories, Quart. Appl. Math., 41, 468-475 (1983) · Zbl 0514.35038
[14] Dehghan, M., Fully implicit finite difference methods for two-dimensional diffusion with a non-local boundary condition, Journal of Computational and Applied Mathematics, 106, 255-269 (1999) · Zbl 0931.65091
[15] Dehghan, M., Numerical solution of the three-dimensional parabolic equation with an integral condition, Numer. Meth. Partial Differen. Equat., 18, 2, 193-202 (2002) · Zbl 0996.65079
[16] Deng, K., Comparison principle for some nonlocal problems, Quart. Appl. Math., 50, 3, 517-522 (1992) · Zbl 0777.35006
[17] Ekolin, G., Finite difference methods for a non-local boundary value problem for the heat equation, BIT, 31, 2, 245-255 (1991) · Zbl 0738.65074
[18] Ewing, R. E.; Lazarov, R. D.; Lin, Y., Finite volume element approximations of nonlocal reactive flows in porous media, Numer. Meth. Partial Differen. Equat., 16, 285-311 (2000) · Zbl 0961.76050
[19] Fairweather, G.; Saylor, R. D., The reformulation and numerical solution of certain nonclassical initial-boundary value problems, SIAM. J. Sci. Stat. Comput., 12, 1, 127-144 (1991) · Zbl 0722.65062
[20] Friedman, A., Monotonic decay of solutions of parabolic equation with nonlocal boundary conditions, Quart. Appl. Math., 44, 401-407 (1986) · Zbl 0631.35041
[21] Gourlay, A. R.; Li Morris, J., The extrapolation of first-order methods for parabolic partial differential equations, II, SIAM, J. Numer. Anal., 17, 641-655 (1980) · Zbl 0444.65058
[22] Gumel, A. B., On the numerical solution of the diffusion equation subject to the specification of mass, J. Aust. Math. Soc. Ser. B, 40, 475-483 (1999) · Zbl 0962.65078
[23] Lambert, J. D., Numerical methods for ordinary differential systems: The initial-value problem (1991), John Wiley and Sons: John Wiley and Sons Chichester · Zbl 0745.65049
[24] Murthy, A. S.V.; Verwer, J. G., Solving parabolic integro-differential equations by an explicit integration method, Journal of Computational and Applied Mathematics, 39, 121-132 (1992) · Zbl 0746.65102
[25] Renardy, M.; Hrusa, W.; Nohel, J. A., Mathematical Problems in Viscoelasticity (1987), Longman: Longman England · Zbl 0719.73013
[26] Shelukhin, V. V., A non-local in time model for radionuclides propagation in Stokes fluids, dynamics of fluids with free boundaries, Siberian Russian Acad. Sci., Inst. Hydrodyn., 107, 180-193 (1993) · Zbl 0831.76086
[27] Taj, M. S.A.; Twizell, E. H., A family of fourth-order splitting methods for parabolic partial differential equations, Numer. Meth. Partial Differen. Equat., 13, 357-373 (1997) · Zbl 0881.65085
[28] Voss, D. A.; Khaliq, A. Q.M., Time splitting algorithms for semidiscretized linear parabolic PDEs based on rational approximants with distinct poles, Adv. Comput. Math., 6, 353-363 (1996) · Zbl 0872.65091
[29] Wang, S.; Lin, Y., A numerical method for the diffusion equation with nonlocal boundary specifications, Int. J. Engng. Sci., 28, 6, 543-546 (1990) · Zbl 0718.76096
[30] Dehghan, M., Saulyer’s techniques for solving a parabolic equation with a non linear boundary specification, Intern. J. Computer. Math., 80, 257-265 (2003) · Zbl 1018.65104
[31] Taj, M. S.A.; Twizell, E. H., A family of third-order parallel splitting methods for parabolic partial differential equations, Intern. J. Computer Math., 67, 411-433 (1998) · Zbl 0908.65086
[32] Twizell, E. H., Computational Methods for Partial Differential Equations (1984), Horwood/Wiley: Horwood/Wiley Chichester · Zbl 0546.65062
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.