Transient Conduction in Halfspace: Robin Problem
M.M. Yovanovich
HSROBIN.MWS
One-dimensional conduction in halfspace. The equation is in the halfspace for time .
The initial condition is uniform temperature throughout the halfspace, . The boundary conditions are (i) as
, , and at the free surface , the condition must be satisfied where . The solution is found by means of Laplace Transform methods. The solution is
> restart:
> case1:= (Ti=300, Tf = 350, k=20, alpha=20e-6, h=100);
Solution for Robin condition.
> T:= Ti + (Tf-Ti)*(erfc(x/(2*sqrt(alpha*t))) - exp(h*x/k + h^2*alpha*t/k^2)*erfc(x/(2*sqrt(alpha*t)) + h*sqrt(alpha*t)/k));
Instantaneous surface temperature rise.
> T0:= subs(x=0, T);
Instantaneous surface heat flux.
> q0:= subs(x=0, -k*diff(T,x));
Temperature rise for a particular problem.
> T1:= subs(case1, T);
For case 1, temperature can be calculated for any position given the time, or the time can be calculated given the position and the temperature. Numerical solutions are required.
> evalf(subs(x=0/1000, t=300, T1), 6);
> time1:= fsolve(simplify(subs(x=10/1000, T1)) = 330, t, 0..3000);
>