Green's function number

In mathematical heat conduction, the Green's function number is used to uniquely categorize certain fundamental solutions of the heat equation to make existing solutions easier to identify, store, and retrieve.

Background

Numbers have long been used to identify types of boundary conditions.[1][2][3] The Green's function number system was proposed by Beck and Litkouhi in 1988[4] and has seen increasing use since then.[5][6][7][8][9] The number system has been used to catalog a large collection of Green's functions and related solutions.[10][11] [12]

Although described here for solutions of the heat equation, this number system could also be used for any phenomena described by differential equations such as diffusion, acoustics, electromagnetics, fluid dynamics, etc.

Notation

The Green's function number specifies the coordinate system and the type of boundary conditions that a Green's function satisfies. The Green's function number has two parts, a letter designation followed by a number designation. The letter(s) designate the coordinate system while the numbers designate the type of boundary conditions that are satisfied.

Table 1. Boundary conditions designations for Green's function number system.

Name	Boundary condition	Number
No physical boundary	G is bounded	0
Dirichlet	${\displaystyle G=0}$	1
Neumann	${\displaystyle {\frac {\partial G}{\partial n}}=0}$	2
Robin	${\displaystyle k{\frac {\partial G}{\partial n}}+hG=0}$	3

Some of the designations for the Greens function number system are given next. Coordinate system designations include: X, Y, and Z for Cartesian coordinates; R, Z, ${\displaystyle \phi }$ for cylindrical coordinates; and, RS, ${\displaystyle \phi }$, ${\displaystyle \theta }$ for spherical coordinates. Designations for several boundary conditions are given in Table 1. The zeroth boundary condition is important for identifying the presence of a coordinate boundary where no physical boundary exists, for example, far away in a semi-infinite body or at the center of a cylindrical or spherical body.

Examples in Cartesian coordinates

X11

As an example, number X11 denotes the Green's function that satisfies the heat equation in the domain (0 < x < L) for boundary conditions of type 1 (Dirichlet) at both boundaries x = 0 and x = L. Here X denotes the Cartesian coordinate and 11 denotes the type 1 boundary condition at both sides of the body. The boundary value problem for the X11 Green's function is given by

${\displaystyle {\begin{aligned}{\dfrac {\partial ^{2}G}{\partial x^{2}}}+{\dfrac {1}{\alpha }}\delta (t-\tau )\delta (x-x')={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};&\;\;0<x<L;\;\;\;t>0\\G|_{x=0}=0;\;\;G|_{x=L}=0;\;\;G|_{t<\tau }=0.&\end{aligned}}}$

Here ${\displaystyle \alpha }$ is the thermal diffusivity (m²/s) and ${\displaystyle \delta }$ is the Dirac delta function.

X20

As another Cartesian example, number X20 denotes the Green's function in the semi-infinite body (${\displaystyle 0<x<\infty }$) with a Neumann (type 2) boundary at x = 0. Here X denotes the Cartesian coordinate, 2 denotes the type 2 boundary condition at x = 0 and 0 denotes the zeroth type boundary condition (boundedness) at ${\displaystyle x=\infty }$. The boundary value problem for the X20 Green's function is given by

${\displaystyle {\begin{aligned}{\dfrac {\partial ^{2}G}{\partial x^{2}}}+{\dfrac {1}{\alpha }}\delta (t-\tau )\delta (x-x')={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};&\;\;0<x<\infty$ ;\;\;\;t>0\\{\dfrac {\partial G}{\partial n}}|_{x=0}=0;\;\;G|_{x\rightarrow \infty }{\mbox{ is bounded}};\;\;G|_{t<\tau }=0.&\end{aligned}}}

X10Y20

As a two-dimensional example, number X10Y20 denotes the Green's function in the quarter-infinite body (${\displaystyle 0<x<\infty }$, ${\displaystyle 0<y<\infty }$) with a Dirichlet (type 1) boundary at x = 0 and a Neumann (type 2) boundary at y = 0. The boundary value problem for the X10Y20 Green's function is given by

${\displaystyle {\begin{aligned}{\dfrac {\partial ^{2}G}{\partial x^{2}}}+{\dfrac {\partial ^{2}G}{\partial y^{2}}}+{\dfrac {1}{\alpha }}\delta (t-\tau )\delta (x-x')\delta (y-y')={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};&\;\;0<x<\infty$ ;\;\;0<y<\infty ;\;\;\;t>0\\G|_{x=0}=0;\;\;{\dfrac {\partial G}{\partial n}}|_{y=0}=0;\;\;&\\G|_{x\rightarrow \infty }{\mbox{ is bounded}};\;\;G|_{y\rightarrow \infty }{\mbox{ is bounded}};\;\;G|_{t<\tau }=0.&\end{aligned}}}

Examples in cylindrical coordinates

R03

As an example in the cylindrical coordinate system, number R03 denotes the Green's function that satisfies the heat equation in the solid cylinder (0 < r < a) with a boundary condition of type 3 (Robin) at r = a. Here letter R denotes the cylindrical coordinate system, number 0 denotes the zeroth boundary condition (boundedness) at the center of the cylinder (r = 0), and number 3 denotes the type 3 (Robin) boundary condition at r = a. The boundary value problem for R03 Green's function is given by

${\displaystyle {\dfrac {1}{r}}{\dfrac {\partial }{\partial r}}\left(r{\dfrac {\partial G}{\partial r}}\right)+{\dfrac {1}{\alpha }}\delta (t-\tau ){\dfrac {\delta (r-r')}{2\pi r'}}={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};\;\;0<r<a;\;\;\;t>0}$

${\displaystyle G|_{r=0}{\mbox{ is bounded}};\;\;k{\dfrac {\partial G}{\partial n}}|_{r=a}+hG|_{r=a}=0;\;\;G|_{t<\tau }=0.}$

Here ${\displaystyle k}$ is thermal conductivity (W/(m K)) and ${\displaystyle h}$ is the heat transfer coefficient (W/(m² K)).

R10

As another example, number R10 denotes the Green's function in a large body containing a cylindrical void (a < r < ${\displaystyle \infty }$) with a type 1 (Dirichlet) boundary condition at r = a. Again letter R denotes the cylindrical coordinate system, number 1 denotes the type 1 boundary at r = a, and number 0 denotes the type zero boundary (boundedness) at large values of r. The boundary value problem for the R10 Green's function is given by

${\displaystyle {\begin{aligned}{\dfrac {1}{r}}{\dfrac {\partial }{\partial r}}\left(r{\dfrac {\partial G}{\partial r}}\right)+{\dfrac {1}{\alpha }}\delta (t-\tau ){\dfrac {\delta (r-r')}{2\pi r'}}={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};&\;\;a<r<\infty$ ;\;\;\;t>0\\G|_{r=a}=0;\;\;G|_{r\rightarrow \infty }{\mbox{ is bounded}};\;\;G|_{t<\tau }=0.&\end{aligned}}}

R01${\displaystyle \phi }$00

As a two dimensional example, number R01${\displaystyle \phi }$00 denotes the Green's function in a solid cylinder with angular dependence, with a type 1 (Dirichlet) boundary condition at r = a. Here letter ${\displaystyle \phi }$ denotes the angular coordinate, and numbers 00 denote the type zero boundaries for angle; here no physical boundary takes the form of the periodic boundary condition. The boundary value problem for the R01${\displaystyle \phi }$00 Green's function is given by

${\displaystyle {\begin{aligned}{\dfrac {1}{r}}{\dfrac {\partial }{\partial r}}\left(r{\dfrac {\partial G}{\partial r}}\right)+{\dfrac {1}{r^{2}}}{\dfrac {\partial ^{2}G}{\partial \phi ^{2}}}+{\dfrac {1}{\alpha }}\delta (t-\tau ){\dfrac {\delta (r-r')}{2\pi r'}}\delta (\phi -\phi ')={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};&\;\;0<r<a;\;\;\;0<\phi <2\pi$ ;\;\;\;t>0\\G|_{r=0}{\mbox{ is bounded,}}\;\;G|_{r=a}=0;&\\G|_{\phi =0}=G|_{\phi =2\pi };{\dfrac {\partial G}{\partial \phi }}|_{\phi =0}={\dfrac {\partial G}{\partial \phi }}|_{\phi =2\pi };\;\;G|_{t<\tau }=0.&\end{aligned}}}

Example in spherical coordinates

RS02

As an example in the spherical coordinate system, number RS02 denotes the Green's function for a solid sphere (0 < r < b ) with a type 2 (Neumann) boundary condition at r = b. Here letters RS denote the radial-spherical coordinate system, number 0 denotes the zeroth boundary condition (boundedness) at r=0, and number 2 denotes the type 2 boundary at r = b. The boundary value problem for the RS02 Green's function is given by

${\displaystyle {\dfrac {1}{r^{2}}}{\dfrac {\partial }{\partial r}}\left(r^{2}{\dfrac {\partial G}{\partial r}}\right)+{\dfrac {1}{\alpha }}\delta (t-\tau ){\dfrac {\delta (r-r')}{4\pi r^{2}}}={\dfrac {1}{\alpha }}{\dfrac {\partial G}{\partial t}};\;\;0<r<b;\;\;\;t>0}$

${\displaystyle G|_{r=0}{\mbox{ is bounded}};\;\;{\dfrac {\partial G}{\partial n}}|_{r=b}=0;\;\;G|_{t<\tau }=0.}$

See also

• Fundamental solution
• Dirichlet boundary condition
• Neumann boundary condition
• Robin boundary condition
• Heat equation

References

1. Luikov, A. V. (1968) Analytical Heat Diffusion Theory, Academic Press, ISBN 0124597564.
2. Ozisik, M. N. (1980) Heat Conduction, John Wiley, p. 13, ISBN 047105481X.
3. Nowak, A., Bialecki R., and Kurpisz, K. (1987) Evaluating eigenvalues for boundary value problems of heat conduction in rectangular and cylindrical coordinates, Int. J for Numerical Methods in Engineering, 24, 419 - 445.
4. Beck, J. V. and Litkouhi, B, (1988) Heat conduction number system, International Journal of Heat and Mass Transfer, 31, 505-515.
5. Al-Nimr, M. A. and Alkam, M. K. (1997) A generalized thermal boundary condition, Heat and Mass Transfer, v. 33, pp. 157 – 161.
6. De Monte, F. (2006) Multi-layer transient heat conduction using transition time scales, Int. Journal Thermal Sciences, v. 45, pp. 882 – 892.
7. Lefebvre, G. (2010) A general modal-based numerical simulation of transient heat conduction in a one dimensional homogeneous slab, Energy and Buildings, v. 42, no. 12, pp. 2309 – 2322.
8. Sarkar, D. and Haji-Sheikh, A. (2012) A view of the thermal wave behaviors in thin plates, International Communications in Heat and Mass Transfer, v. 39, No. 8, pp. 1009-1017.
9. Zhou, Y. (2012) Heat conduction in a semi-infinite body with power-type initial and boundary conditions, International Journal of Thermophysics, v. 33, No. 12, pp. 2390-2406.
10. Cole, K.D., Beck, J. V., Haji-Sheikh, A., and Litkouhi, B. (2011), Heat Conduction Using Green's Functions,Taylor and Francis, (2nd ed.) ISBN 9781439813546.
11. Green's Function Library, http://www.greensfunction.unl.edu/
12. Exact Analytical Conduction Toolbox, http://exact.unl.edu/