CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Energy equation with convective heat transfer.

Register Blogs Community New Posts Updated Threads Search

Like Tree6Likes
  • 1 Post By arungovindneelan
  • 1 Post By FMDenaro
  • 1 Post By FMDenaro
  • 1 Post By sbaffini
  • 1 Post By Lucas_Gasparino
  • 1 Post By Lucas_Gasparino

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 23, 2017, 22:52
Default Energy equation with convective heat transfer.
  #1
Member
 
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14
arungovindneelan is on a distinguished road
Hi,
In books, while driving conservative equations we use control-volume approach and in the energy equation, we consider only conduction equation (Fourier’s law). Please tell me how can we incorporate convection into it.

http://www.bakker.org/dartmouth06/engs150/03-cnsrv.pdf

I think we can't directly put q_{convection}=h*A_surface*(T1-T2) into it because convection is surface phenomena, not a volumetric phenomenon. For example, let us assume we have one quad element at the interface between the fluid and solid, where the cross-sectional area ( direction of grad(T)) of the face of the element is not equal to surface area so assuming

q=q_{convective} +q_{conductive} could introduce some error. Am I right?

Please provide some link where they incorporate all heat transfer mode in the energy equation.
arungovindneelan is offline   Reply With Quote

Old   August 24, 2017, 04:39
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
If you are considering convection as that appearing in conjugate heat transfer that is a BC.s applied to the surface that ensures that q is the same at the interface. The equation for the total energy has only a general expression for the surface integral of the normal compononent of the q vector.
FMDenaro is offline   Reply With Quote

Old   August 24, 2017, 05:52
Default
  #3
Member
 
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14
arungovindneelan is on a distinguished road
In my understanding direction of conductive heat transfer is along gradient(T), since we are not having complete Temperature field we can't calculate grad(T) but at least we can derive that locally.



In the figure "s" stands for solid and "f" stands for fluid.

In figure (a) we have split the solid block into two elements(mesh) using slant line, heat transfer between those elements should be same as the figure (b), so the area to calculate heat transfer rate is the cross-sectional area.

Figure (c) and (d) having interface between fluid and solid, where conduction area should be same as figure(a) and (b) but convection area is wetted area that will differ from figure (c) and (d) so convection heat transfer will be higher in figure (c) then figure (d)

In figure(c) and (d), I think conduction heat transfer area is lower than convection but we know B.C we use never take care of that for example

This is one of the examples I could think of. My intention is: Is there any pre-ready Energy equation available that take care of all heat transfer mode in it.

You may consider a fluid element, that can transfer heat to another element by all mode of heat transfer (convection and radiation could be dominant). How can I account for all heat transfer mode in an element?
granzer likes this.
arungovindneelan is offline   Reply With Quote

Old   August 24, 2017, 06:57
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Generally the continuity of the flux is prescribed at the interface solid/fluid.
If you want to have an idea of general heat transfer models that can be practically applied in CFD I suggest a reading of
http://www.afs.enea.it/fluent/Public.../PDF/chp11.pdf

Then, you can focus on theoretical aspects on texbooks of heat transfer.
arungovindneelan likes this.
FMDenaro is offline   Reply With Quote

Old   August 24, 2017, 07:13
Default
  #5
Member
 
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14
arungovindneelan is on a distinguished road
What I have understood so for is, there is no need to put explicit convection heat transfer in energy equation because that would be taken care by convective part of the energy equation (but I didn't understand the explicit relation between convection part and h*(T1-T2)) and radiation is treated as a source term.

Thank you.
arungovindneelan is offline   Reply With Quote

Old   August 24, 2017, 07:13
Default
  #6
Member
 
Ali Noaman Ibrahim
Join Date: Sep 2015
Location: US_Chicago
Posts: 97
Rep Power: 11
alinuman15 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Generally the continuity of the flux is prescribed at the interface solid/fluid.
If you want to have an idea of general heat transfer models that can be practically applied in CFD I suggest a reading of
http://www.afs.enea.it/fluent/Public.../PDF/chp11.pdf

Then, you can focus on theoretical aspects on texbooks of heat transfer.


Thank you prof. Filippo for that nice reference for H.T. Do you know any reference relates heat transfer to phase change phenomena in connection with OpenFoam capabilities?
Thank you in advance!


Sent from my iPhone using CFD Online Forum mobile app
alinuman15 is offline   Reply With Quote

Old   August 24, 2017, 07:15
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by arungovindneelan View Post
Sir link is broken, I google that link and read that. That may be useful to others. You could edit it.
What I have understood so for is, there is no need to put explicit convection heat transfer in energy equation because that would be taken care by convective part of the energy equation (but I didn't understand the explicit relation between convection part and h*(T1-T2)) and radiation is treated as a source term.

Thank you.

The link is working on my PC ...
FMDenaro is offline   Reply With Quote

Old   August 24, 2017, 07:28
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by alinuman15 View Post
Thank you prof. Filippo for that nice reference for H.T. Do you know any reference relates heat transfer to phase change phenomena in connection with OpenFoam capabilities?
Thank you in advance!


Sent from my iPhone using CFD Online Forum mobile app
No, sorry...
FMDenaro is offline   Reply With Quote

Old   August 24, 2017, 07:39
Default
  #9
Member
 
Ali Noaman Ibrahim
Join Date: Sep 2015
Location: US_Chicago
Posts: 97
Rep Power: 11
alinuman15 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
No, sorry...


Thanks for replying!
Actually, I am asking because we have very nice experimental setup for two phase flow in our lab. Therefore, it would be great if we integrate that with OpenFoam numerical capabilities.


Sent from my iPhone using CFD Online Forum mobile app
alinuman15 is offline   Reply With Quote

Old   August 24, 2017, 07:43
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Try to ask in the OF forum
alinuman15 likes this.
FMDenaro is offline   Reply With Quote

Old   August 24, 2017, 08:26
Default
  #11
Member
 
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14
arungovindneelan is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The link is working on my PC ...
Probably my internet might be disconnected at that time .
arungovindneelan is offline   Reply With Quote

Old   August 25, 2017, 06:22
Default
  #12
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,190
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
As already pointed out by Filippo, the convective part in the energy equation just contains all the convection. Expressing it as h(T-T_inf) is just a mean to handle overall convection effects at the system level or, for numerical simulations, when you don't actually want to simulate convection (i.e., your domain consists of the solid part only and you approximate convection in the fluid as a bc). Note however that, for actual simulations with such fluid-solid interfaces, you typically need an actual grid interface to separate them. Then the only difference among the different cases ends up being in the different numerical accuracies (typically lower for distorted interfaces).
granzer likes this.
sbaffini is offline   Reply With Quote

Old   August 25, 2017, 06:28
Default
  #13
New Member
 
Lucas Gasparino
Join Date: Jul 2017
Location: Swansea, UK
Posts: 23
Rep Power: 9
Lucas_Gasparino is on a distinguished road
Simply put, convection goes into energy equation as a flux term qh= hdT, where dT is simply temperature variation. Now, if you're implying natural convection and Boussinesq assumption for buoyancy, then this is a whole new matter; still, the heat part is as I described. Frank Incropera has a good chapter in convection, and O. C. Zienkiewicz The FE set vol3 states the full energy equation and how to discretise it.

Sent from my XT1635-02 using CFD Online Forum mobile app
granzer likes this.
Lucas_Gasparino is offline   Reply With Quote

Old   August 25, 2017, 06:47
Default
  #14
New Member
 
Lucas Gasparino
Join Date: Jul 2017
Location: Swansea, UK
Posts: 23
Rep Power: 9
Lucas_Gasparino is on a distinguished road
Well, I was trying to post the picture, but it's not happening... So here it is:

Convection is a heat flux, just like radiation (in the surface sense); thus, it crosses the volume boundary as a flux qh, and is defined as is, i.e., don't need to change the classical definition. Now, of course the area is dA, meaning that this is simply a boundary term. When a weak form is obtained for the energy equation, this term adds to other Neumann BCs. Then, the difference lies in that one temperature is unknown; this is a mixed, or Cauchy condition, in which both a Dirichlet and a Neumann term are collected into the surface. So, when a discrete form Ku=f is achieved, the known term hTref is a known flux and moves to the rhs. The unknown term modifies part of the element matrix in which it occurs. Now, the catch is that h depends on the flow: so a coupling relation with the momentum equation is introduced to interact h, T and flow properties (velocity and buoyancy). In a thermal FE software, you assume h known, and the coupling disappears, as the flux is well defined. But it is always driven by temperature difference, obeying zeroth law of thermodynamics. Hope it helps!

Sent from my XT1635-02 using CFD Online Forum mobile app
granzer likes this.
Lucas_Gasparino is offline   Reply With Quote

Reply

Tags
conservative form, energy balance


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
Multiphase flow - incorrect velocity on inlet Mike_Tom CFX 6 September 29, 2016 02:27
Separate Convective and Radiative heat transfer in CFD post using fluent as a solver Lemanes FLUENT 1 July 6, 2015 11:31
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 07:28
Mass transfer Term in Mixture Energy Equation CFDtoy Main CFD Forum 0 July 3, 2008 17:14


All times are GMT -4. The time now is 12:40.