looking for help w/ heat equation + surface convection
looking for help w/ heat equation + surface convection
(OP)
Hi,
I have an application in which I have a semi-infinite solid subjected to surface convection at the x=0 boundary and a time varying fluid temperature. I am familiar with the heat equation (d2Tdx2 = 1/alpha dTdt) and I've been able to solve my problem by discretizing the spatial domain and time marching with an rk4 solver, but I am trying a method which does not require tracking as many state variables as I potentially need to solve for many different semi-infinite bodies at once.
I came across an approximate integral method (ref. Ozisik, Heat Conduction) in which an approximate spatial distribution is enforced over a thermal boundary distance (delta) and the heat equation is reduced to solving a differential equation describing delta.
My problem is actually integrating the equation. At time zero the surface temperature rate of change is infinite (also the thermal boundary layer has zero thickness), which I've confirmed by looking at available exact solutions for my problem which have a constant fluid temperature. This equates to a thermal boundary rate initially at infinity. At any time after t=0 there will be a positive boundary layer thickness and the equations won't blow up.
How do I handle time zero and getting the whole thing started? The text just leave it at "Equation (9-36) can be integrated numerically".
Thanks for any insight anyone can offer.
I have an application in which I have a semi-infinite solid subjected to surface convection at the x=0 boundary and a time varying fluid temperature. I am familiar with the heat equation (d2Tdx2 = 1/alpha dTdt) and I've been able to solve my problem by discretizing the spatial domain and time marching with an rk4 solver, but I am trying a method which does not require tracking as many state variables as I potentially need to solve for many different semi-infinite bodies at once.
I came across an approximate integral method (ref. Ozisik, Heat Conduction) in which an approximate spatial distribution is enforced over a thermal boundary distance (delta) and the heat equation is reduced to solving a differential equation describing delta.
My problem is actually integrating the equation. At time zero the surface temperature rate of change is infinite (also the thermal boundary layer has zero thickness), which I've confirmed by looking at available exact solutions for my problem which have a constant fluid temperature. This equates to a thermal boundary rate initially at infinity. At any time after t=0 there will be a positive boundary layer thickness and the equations won't blow up.
How do I handle time zero and getting the whole thing started? The text just leave it at "Equation (9-36) can be integrated numerically".
Thanks for any insight anyone can offer.





RE: looking for help w/ heat equation + surface convection
Heat transfer requires approximations and assumptions for everything else.
RE: looking for help w/ heat equation + surface convection
RE: looking for help w/ heat equation + surface convection
Your time interval for heat transfer in a real-world large application (unless it is an explosion or internal combustion engine fuel combustion) is on the order of minutes and tens of minutes not milliseconds. So I'd try to stay away from the very, very tiny increments that will distort your starting conditions.