Function SemiInf8(T_i, q``_dot_s, alpha, k, a, r, time) returns the temperature within a semi-infinite material that is located outside of a spherical surface and is subjected to a step change in the surface heat flux. Solution is from Carslaw and Jaeger as reported in the Handbook of Heat Transfer.

The calling protocol is:

T=SemiInf8(T_i, q``_dot_s, alpha, k, a, r, time)

T_i = initial temperature of the material [C] or [K]

q``_dot_s = surface heat flux at r = a [W/m^2]

alpha = thermal diffusivity [m^2/s]

k = thermal conductivity [W/m-K]

a = radius of sphere [m]

r = radial position [m]

time = time relative to beginning of surface disturbance [s]

This function can be used with English units set in EES. In this case, T_i is in [F] or [R], q``_dot_s is in [Btu/hr-ft^2], alpha is in [ft^2/hr], k is in [Btu/hr-ft], a is in [ft], r is in [ft], and time is in [s].

T is the temperature in [C] or [K] (or [F] or [R] in English units)

Rohsenow, W.M., J.P. Hartnett, and Y.I. Cho,

$UnitSystem SI Mass J K Pa

$VarInfo T units=K

T_i=300 [K]

q``_dot_s=1e4 [W/m^2]

alpha=1e-4 [m^2/s]

k=1.5 [W/m-K]

a=0.01 [m]

r=0.02 [m]

time=1 [s]

T=

{

T=307.6 [K]}