Generalized Trapezoidal Formulas For The Symmetric Heat Equation In Polar Coordinates II. The Case Of Point Sources

• Published in: International journal of computer mathematics 2002-01, Vol.79 (11), p.1187-1200
• Main Authors: Chawla, M. M., Al-Zanaidi, M. A., Evans, D. J.
• Format: Article
• Language: English
• Subjects: Continuous Point Sources; Exact sciences and technology; Generalized Finite Hankel Transforms; Generalized Trapezoidal Formulas; Mathematics; Numerical analysis; Numerical analysis. Scientific computation; Partial differential equations, miscellaneous problems; Polar Coordinates; Sciences and techniques of general use; Symmetric Heat Equation
• ISSN: 0020-7160; 1029-0265
• DOI: 10.1080/00207160213940

The present paper is in continuation of our previous paper Chawla et al. [4]. An important class of applications of the radially symmetric heat equation in polar coordinates: u_{t}=v ( u_{rr}\,+\; ( a/r ) u_{r} ), involve the presence of a continuous point source of heat at the centre of the sphere ( a=2 ) or on the axis of the cylinder ( a=1 ). This necessitates a modification of the radial grid used in [4]; our modification of the radial grid in the present paper accommodates a point source of heat at r=0 in a natural way. We then extend generalized trapezoidal formulas GTF( \alpha ) for the one-step time integration of these problems. Again, with the help of the generalized finite Hankel transforms introduced in [4] we are able to obtain, in a natural way, analytical solutions of the heat equation in the presence of point sources of heat for both the cases a=1 and a=2 . Numerical experiments are provided to compare the performance of the GTF( \alpha ) time integration scheme with the schemes based on the backward Euler and the classical arithmetic-mean trapezoidal formula.