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]} The number system has been used to catalog a large collection of Green's functions and related solutions. ^{[9] [10] [11]}

Although the examples given below are for the heat equation, this number system applies to 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, φ for cylindrical coordinates; and, RS, φ, θ 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

Equation	${\displaystyle {\frac {\partial ^{2}G}{\partial x^{2}}}+{\frac {1}{\alpha }}\delta (t-\tau )\delta (x-x')={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<x<L}$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(x=0)=0}$, ${\displaystyle G(x=L)=0}$ ${\displaystyle G(t<\tau )=0}$

Here ${\displaystyle \alpha }$ is the thermal diffusivity (m²/s) and ${\displaystyle \delta }$ is the Dirac delta function. This GF is developed elsewhere. ^{[12] [13]}

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

Equation	${\displaystyle {\frac {\partial ^{2}G}{\partial x^{2}}}+{\frac {1}{\alpha }}\delta (t-\tau )\delta (x-x')={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<x<\infty }$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle \left.{\frac {\partial G}{\partial n}}\right|_{x=0}=0}$, ${\displaystyle \left.G\right|_{x\to \infty }}$ is bounded, ${\displaystyle G(t<\tau )=0}$

This GF is published elsewhere. ^{[14] [15]}

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

Equation	${\displaystyle {\frac {\partial ^{2}G}{\partial x^{2}}}+{\frac {\partial ^{2}G}{\partial y^{2}}}+{\frac {1}{\alpha }}\delta (t-\tau )\delta (x-x')\delta (y-y')={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<x<\infty }$ ${\displaystyle 0<y<\infty }$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(x=0)=0}$, ${\displaystyle \left.G\right|_{x\to \infty }}$ is bounded, ${\displaystyle \left.{\frac {\partial G}{\partial n}}\right|_{y=0}=0}$, ${\displaystyle G|_{y\to \infty }}$ is bounded, ${\displaystyle G(t<\tau )=0}$

Applications of related half-space and quarter-space GF are available. ^[16]

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

Equation	${\displaystyle {\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial G}{\partial r}}\right)+{\frac {1}{\alpha }}\delta (t-\tau ){\frac {\delta (r-r')}{2\pi r'}}={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<r<a}$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(r=0)}$ is bounded, ${\displaystyle \left.k{\frac {\partial G}{\partial n}}\right|_{r=a}+\left.hG\right|_{r=a}=0}$, ${\displaystyle G(t<\tau )=0}$

Here ${\displaystyle k}$ is thermal conductivity (W/(m K)) and ${\displaystyle h}$ is the heat transfer coefficient (W/(m² K)). See Carslaw & Jaeger (1959 , p. 369), Cole et al. (2011 , p. 543) for this GF.

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

Equation	${\displaystyle {\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial G}{\partial r}}\right)+{\frac {1}{\alpha }}\delta (t-\tau ){\frac {\delta (r-r')}{2\pi r'}}={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle a<r<\infty }$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(r=a)=0}$, ${\displaystyle G|_{r\to \infty }}$ is bounded, ${\displaystyle G(t<\tau )=0}$

This GF is available elsewhere. ^{[17] [18]}

R01φ00

As a two dimensional example, number R01φ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 φ denotes the angular (azimuthal) 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φ00 Green's function is given by

Equation	${\displaystyle {\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial G}{\partial r}}\right)+{\frac {1}{r^{2}}}{\frac {\partial ^{2}G}{\partial \phi ^{2}}}+{\frac {1}{\alpha }}\delta (t-\tau ){\frac {\delta (r-r')}{2\pi r'}}\delta (\phi -\phi ')={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<r<a}$ ${\displaystyle 0<\phi <2\pi }$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(r=0)}$ is bounded, ${\displaystyle G(r=a)=0}$, ${\displaystyle G(\phi =0)=G(\phi =2\pi )}$, ${\displaystyle \left.{\frac {\partial G}{\partial \phi }}\right|_{\phi =0}=\left.{\frac {\partial G}{\partial \phi }}\right|_{\phi =2\pi }}$ ${\displaystyle G(t<\tau )=0}$

Both a transient ^[19] and steady form ^[20] of this GF are available.

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

Equation	${\displaystyle {\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial G}{\partial r}}\right)+{\frac {1}{\alpha }}\delta (t-\tau ){\frac {\delta (r-r')}{4\pi r^{2}}}={\frac {1}{\alpha }}{\frac {\partial G}{\partial t}}}$
Domain	${\displaystyle 0<r<b}$ ${\displaystyle t>0}$
Boundary conditions	${\displaystyle G(r=0)}$ is bounded, ${\displaystyle \left.{\frac {\partial G}{\partial n}}\right|_{r=a}=0}$, ${\displaystyle G(t<\tau )=0}$

This GF is available elsewhere. ^[21]

See also

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

References

1. Luikov 1968
2. Özışık, M. Necati (1980). Heat conduction (1st ed.). New York: Wiley. ISBN   047105481X.
3. Nowak, A.; Białecki, R.; Kurpisz, K. (February 1987). "Evaluating eigenvalues for boundary value problems of heat conduction in rectangular and cylindrical co-ordinate systems". International Journal for Numerical Methods in Engineering. 24 (2): 419–445. doi:10.1002/nme.1620240210.
4. Beck, James V.; Litkouhi, Bahman (March 1988). "Heat conduction numbering system for basic geometries". International Journal of Heat and Mass Transfer. 31 (3): 505–515. doi:10.1016/0017-9310(88)90032-4.
5. Al-Nimr, M. A.; Alkam, M. K. (19 September 1997). "A generalized thermal boundary condition". Heat and Mass Transfer. 33 (1–2): 157–161. doi:10.1007/s002310050173. S2CID   119549322.
6. de Monte, Filippo (September 2006). "Multi-layer transient heat conduction using transition time scales". International Journal of Thermal Sciences. 45 (9): 882–892. doi:10.1016/j.ijthermalsci.2005.11.006.
7. Lefebvre, G. (December 2010). "A general modal-based numerical simulation of transient heat conduction in a one-dimensional homogeneous slab". Energy and Buildings. 42 (12): 2309–2322. doi:10.1016/j.enbuild.2010.07.024.
8. Toptan, A.; Porter, N. W.; Hales, J. D. (2020). "Construction of a code verification matrix for heat conduction with finite element code applications". Journal of Verification, Validation and Uncertainty Quantification. 5 (4): 041002. doi:10.1115/1.4049037.
9. Cole et al. 2011
10. "Green's Function Library" . Retrieved November 19, 2020.
11. "Exact Analytical Conduction Toolbox" . Retrieved March 4, 2021.
12. Luikov 1968 , p. 388
13. Cole et al. 2011
14. Luikov 1968 , p. 387
15. Carslaw & Jaeger 1959 , p. 276
16. Beck, J. V.; Wright, N.; Haji-Sheikh, A.; Cole, K. D; Amos. D. (2008). "Conduction in rectangular plates with boundary temperatures specified". International Journal of Heat and Mass Transfer. 52 (19–20): 4676–4690. doi:10.1016/j.ijheatmasstransfer.2008.02.020. S2CID   12677235.
17. Carslaw & Jaeger 1959 , p. 378
18. Thambynayagam, R. K. M. (2011). The Diffusion Handbook. McGraw-Hill. p. 432. ISBN   9780071751841.
19. Cole et al. 2011 , p. 554
20. Melnikov, Y. A. (1999). Influence Functions and Matrices. New York: Marcel Dekker. p. 223. ISBN   9780824719418.
21. Cole et al. 2011 , p. 309

• Carslaw, H. S.; Jaeger, J. C. (1959). Conduction of Heat in Solids. Oxford University Press. ISBN   9780198533689.
• Cole, Kevin D.; Beck, James; Haji-Sheikh, A.; Litkouhi, Bahman (2011). Heat Conduction Using Greens Functions (2nd ed.). doi:10.1201/9781439895214. ISBN   9781439813546.
• Luikov, A. V. (1968). Analytical Heat Diffusion Theory. Academic Press. p. 388. ISBN   0124597564.