CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Newton's cooling BC

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   August 20, 2004, 17:14
Default Newton's cooling BC
Posts: n/a
Hi, I am trying to solve a 2-d steady state heat transfer problem with finite difference (successive under-relaxation). I am trying to apply a convective boundary condition on a horizontal surface with a variable convection heat transfer coeffient. i.e. Heat is to convect outside from a horizontal surface with known but variable (in x, of course) convection heat transfer coefficient, h. I expect the temperature on the surface be variable. But it is not. I use this BC:

|K. grad(T)|=h(T-T_infinity)

where || means "magnitude," and K is the conductivity. Then I discretize it: backward in x and forward in y (since the convective surface is a bottom surface) and find T and use it as BC in my equations. It really seemed like an easy problem at first. When I use an scheme O(dx^2, dy^2), backward for x and forward in y, it never converges. But O(dx, dy^2) converges however the temperature is constant on the surface which it should not be.

I really appreciate any help/idea. Thanks, Mohsen
  Reply With Quote

Old   August 22, 2004, 03:30
Default Re: Newton's cooling BC
Posts: n/a

I think the BC should relate the NORMAL flux due to conduction with the convection, i.e.,

K grad(T)*n = h(T-T_infinity) (where *n is the scalar product with the normal unit vector), rather than the BC you use ( |K. grad(T)|=h(T-T_infinity) ).

  Reply With Quote

Old   August 22, 2004, 11:47
Default Re: Newton's cooling BC
Posts: n/a
thanks Rami. I tried that one first. And I think you are right. In that case, I get a uniform temprature on the surface whose heat convection coefficient is non-uniform; which is not correct. It could be a mere coding problem though.
  Reply With Quote

Old   August 23, 2004, 02:59
Default Re: Newton's cooling BC
Posts: n/a

<head> <meta http-equiv=Content-Type content="text/html; charset=windows-1255"> <meta name=Generator content="Microsoft Word 10 (filtered)"> <title>Mohsen,</title>




<body lang=EN-US>

<div class=Section1>

<span style='color:#333333'>Mohsen,</span>

<span style='color:#333333'>Let's apply the BC</span>

<span style='color:#333333'>-k grad(T) n = h(T-T<sub>&infin;</sub>)</span>

<p style='text-align:justify'><span style='color:#333333'>at y=0. For simplicity I use vertex-centered temperature and 1<sup>st</sup> order discretization for the derivative. After some simple manipulation, the discrete boundary temperature, T<sub>i,0</sub>, is</span>

<span style='color:#333333'>T<sub>i,0</sub> <sub>*</sub>= ( T<sub>i,1 </sub>*- Nu T<sub>&infin;</sub>) / (1-Nu)</span>

<p class=MsoNormal><span style='color:#333333'>where Nu, the local Nusselt number is defined as</span>

<p class=MsoNormal><span style='color:#333333'>Nu = h(x) ∆y / k(x,0) .</span>

<p class=MsoNormal><span style='color:#333333'>I hope the derivation is correct (please check) and this may help you to validate your code.</span>



  Reply With Quote

Old   August 23, 2004, 12:11
Default Re: Newton's cooling BC
Posts: n/a
Thanks Rami. Last night I figured all this time the code was giving me correct results! My geometry is super small (micro) therefore the temperature x-gradient is very small. I multiplied the dimensions by 1000 and observed there IS a non-zero x-gradient. I used the same equation as yours for floor and this:

k grad(T) n = h(T-T)

fot ceiling. I used second order backward for x. For ceiling I used second order forward and for floor I used second order backward. And it worked!
  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
About Cooling channel alefem FLOW-3D 1 May 28, 2010 11:15
Regenerative cooling nozzle. Brian NUMECA 3 September 1, 2005 08:53
Regenerative cooling rocket engine Brian Main CFD Forum 0 November 12, 2003 07:59
Regenerative cooling nozzle. Brian Main CFD Forum 0 November 10, 2003 02:11
Passive cooling system for batch reactor Darek Main CFD Forum 0 September 30, 2003 05:06

All times are GMT -4. The time now is 17:03.