Function SemiInf14(T_i, alpha, a, DELTAT, x, time) returns the temperature within a semi-infinite body that has an initial temperature that is non-uniform. The region immediately adjacent to the adiabatic surface (0 < x < a) is at an elevated temperature: T = T_i+DELTAT*(1-x/a) while the remainder of the material (x>a) is at T_i.

The calling protocol is:

T = SemiInf14(T_i, alpha, a, DELTAT, x, time)

T_i = temperature far from the surface [C] or [K]

alpha = thermal diffusivity [m^2/s]

a = extent of the elevated temperature region [m]

DELTAT = temperature elevation at x=0 [C] or [K]

x = position relative to surface [m]

time = time relative to beginning of heat flux cycle [s]

This function can be used with English units set in EES. In this case, T_i is in [F] or [R], alpha is in [ft^2/hr], a is in [ft], DELTAK is in [F] or [R], x is in [ft], and time is in [s].

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

Carslaw, H.S. and J.C. Jaeger,

$UnitSystem SI Mass J K Pa Radian

$VarInfo T units=K

T_i=300 [K]

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

x=0.002 [m]

time=0.1 [s]

a=0.005 [m]

DELTAT=100 [K]

T=

{

T=313.8 [K]

}