Sphere heating by convection - boundary condition(numerical)
Sphere heating by convection - boundary condition(numerical)
(OP)
Howdy,
I've been fighting with coding a numerical solution to the unsteady temperature distribution of a cool solid sphere being heated by forced convection of hot gas surrounding it. The problem for me lies in the surface boundary condition.
I must specify heat flux by equating the heat transfered away from the surface to the heat transferred in by convection:
-k(dT/dr) = h(Tsurface - Tgas)
Now, in finite difference form,
dT/dr = (Ts1 - Ts2)/dr
dr = spacial step size now, Ts1 = surface temp, Ts2 is the temperature one step into the sphere from the surface.
Solving for surface temp:
Ts = ((k/dr/h)*Ts2 + Tgas)/(k/dr/h + 1)
This just doesn't work somehow, and I can't figure out why.. I'm dealing with small spheres, and the above equation gives me a surface temp about equal to that of the next layer in, and the damn sphere heats up very very slowly.. I've also compared it to the analytical solution which confirms it's not right.
Does anybody know of a numerical method for handling this boundary which works?
many thanks,
eric
(sorry for the lengthy post)
I've been fighting with coding a numerical solution to the unsteady temperature distribution of a cool solid sphere being heated by forced convection of hot gas surrounding it. The problem for me lies in the surface boundary condition.
I must specify heat flux by equating the heat transfered away from the surface to the heat transferred in by convection:
-k(dT/dr) = h(Tsurface - Tgas)
Now, in finite difference form,
dT/dr = (Ts1 - Ts2)/dr
dr = spacial step size now, Ts1 = surface temp, Ts2 is the temperature one step into the sphere from the surface.
Solving for surface temp:
Ts = ((k/dr/h)*Ts2 + Tgas)/(k/dr/h + 1)
This just doesn't work somehow, and I can't figure out why.. I'm dealing with small spheres, and the above equation gives me a surface temp about equal to that of the next layer in, and the damn sphere heats up very very slowly.. I've also compared it to the analytical solution which confirms it's not right.
Does anybody know of a numerical method for handling this boundary which works?
many thanks,
eric
(sorry for the lengthy post)





RE: Sphere heating by convection - boundary condition(numerical)
If you consider quasi-steady flow:
Your heat transfer coefficient, h, will vary over the surface. Consider where the flow hits the sphere dead on and where it leaves,180 degrees later. The velocity distribution and therefore its effect on h can be significant.
Further, I believe you are implying with the boundary condition dT/dr = (Ts1 - Ts2)/dr, that the outer surface of the sphere is only a function of time and not coordinate. If that is true, then you are using 1D equation for temp distribution in the sphere which could lead to error.
I would look for quasi-steady heat transfer coef as a function of geometry and model with more than one dimension.
RE: Sphere heating by convection - boundary condition(numerical)
The basis of the standard heat flow equation starts with a source of heat that travels through a series of layers and ends up in a heat sink, which is the air.
The convection coefficient of the air layer presumes an infinite sink at tambient.
Tsurface is a function of the heatflow across the interface and the interior temperature of the sphere, which is not constant. Your equations need to account for that.
TTFN
RE: Sphere heating by convection - boundary condition(numerical)
2) do you have sufficient mesh resolution to solve for the boundary layer?
3) is the flow laminar or turbulent?
4) the heat transfer coeff is not a constant, but a function of space for steady state.
5) at the surface, there's only balance of conduction heat transfer: heat flux of conduction of gas to surface equals to that of solid. another condition is temperature continuity. depend on your gas temperature, you can neglect radiation.
RE: Sphere heating by convection - boundary condition(numerical)
Hope this has helped.
RE: Sphere heating by convection - boundary condition(numerical)
corus
RE: Sphere heating by convection - boundary condition(numerical)
The problem isn't 2D
Clearly, if you assume the sphere temp constant, the convection coef and surrounding ambient constant, then the problem is one dimensional.
Corus--What is your basis for the above assumptions?
RE: Sphere heating by convection - boundary condition(numerical)
You could also consider it as 1 dimensional with spherical co-ordinates, or in 2 dimensions when, as modelled as an arc, it would be axisymmetric. There'd be no need for 2D though with a single heat transfer coefficient and either 1D spherical co-ordinates or my approximate analytical method would do.
corus
RE: Sphere heating by convection - boundary condition(numerical)
I suggest this approach if you compare closed form solutions of a sphere with the approx film coef, same diam and thermal diffusivity. If temp at center of sphere compares with lumped analysis, then lumped is satisfactory way to go.
RE: Sphere heating by convection - boundary condition(numerical)
Sinilar to what you describe, I have seen a formula for the 'biot modulus' based upon the epxression hL/k for a cylinder of length L. What it is for a sphere I have no idea. Howver, the test of whether the approximation leads to negigible error is if this biot modulus is small. What small is, is another matter.
corus
RE: Sphere heating by convection - boundary condition(numerical)