CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

openfoam for heat transfer

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 22, 2010, 16:14
Default openfoam for heat transfer
  #1
New Member
 
Join Date: Jul 2010
Posts: 25
Rep Power: 8
kirankarki is on a distinguished road
Hi everyone,
can anyone tell me if openfoam is capable to handle the heat transfer mechanism through building envelope ? since heat transfer in building should include solar radiation due to outdoor temp, weather that varies with time, conduction through different layers of materials, convection due to natural ventilation air etc that determines the inside temp and heat flux. Which solver is good for this or new solver is to be created myself?

Thanks
kirankarki is offline   Reply With Quote

Old   September 23, 2010, 11:29
Default
  #2
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
You are pretty much in create new solver territory. To treat the wall heat transfer in a transient way, you pretty much have to do some form of conjugate heat transfer as well. If a steady approach to wall heat transfer is ok, then you just need to couple the inside and outside via a layer based resistance.

Let me know if you need a hand, I might have some appropriate code lying around.
eugene is offline   Reply With Quote

Old   September 23, 2010, 11:43
Default For heat transfer thru wall
  #3
New Member
 
Join Date: Jul 2010
Posts: 25
Rep Power: 8
kirankarki is on a distinguished road
Thanks eugene for response,
i m very new to this openfoam and i m construction management major thus no effecient background of programming ,i want to have transient heat transfer through wall and the outside weather data like temp, wind vel changes every hour, if u have some solvers for this case to send me, it will be highly appreciated, also i need to know how many kind of files should i create for system , constant and time directory ,
thanks again
kirankarki is offline   Reply With Quote

Old   September 23, 2010, 19:02
Default
  #4
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
Unfortunately I do not have a ready made solution for you.

There is a however a thermal boundary condition called "wallHeatTransfer" that couples an outside temperature to the inside wall via a fixed transfer coefficient. Unfortunately the BC has several limitations: 1 it wont help you with solar radiation, 2. it assume laminar flow on the inside.

You can use this T boundary in conjunction with the rhoSimpleFoam and other "compressible" solvers.

Good luck,
eugene is offline   Reply With Quote

Old   September 24, 2010, 03:07
Default
  #5
Member
 
Join Date: Nov 2009
Location: Germany
Posts: 96
Rep Power: 9
val46 is on a distinguished road
Hi eugene,

what do you mean by layer based resistance?
Is there a way to set a thermal resistance between 2 regions?

Regards,
Toni
val46 is offline   Reply With Quote

Old   September 24, 2010, 04:53
Default
  #6
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
No, there is no way to currently set layer based resistance between two regions. It isn't that hard to implement if you only want a steady state solution though. If you want a transient solution in the solid, you will need some kind of mesh for the solid.

This is something I would like to to do at some point in the future. If I ever get round to it, I will let you know.
eugene is offline   Reply With Quote

Old   September 24, 2010, 10:30
Default
  #7
New Member
 
Peter Roosen
Join Date: Sep 2010
Posts: 4
Rep Power: 8
petroo is on a distinguished road
Hello,

Quote:
Originally Posted by kirankarki View Post
can anyone tell me if openfoam is capable to handle the heat transfer mechanism through building envelope ? since heat transfer in building should include solar radiation due to outdoor temp, weather that varies with time, conduction through different layers of materials, convection due to natural ventilation air etc that determines the inside temp and heat flux. Which solver is good for this or new solver is to be created myself?
Maybe you might take a look at the EnergyPlus system that is completely dedicated to problems like this. All your requirements are dealt with, including instationary calculation over time. EnergyPlus is capable of individual room-resolved heat load calculations etc., and there is a rather large library of construction materials with their thermodynamic properties.

Regards,

Peter
petroo is offline   Reply With Quote

Old   October 19, 2010, 07:56
Default
  #8
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
hi,
One question:

Can buoyantBousinesqPimpleFoam solver handle a room ventilation simulation with inlet velocity(Fan, Cooler, Heater, ...) and pressure outlet (Door) with temperature (Energy equation)?

i asked this question because i didn't see any tutorial file case same as the case i told and all of them are same as hot room with a hot section without any forced flow like fan.
Best,
Maysam
maysmech is offline   Reply With Quote

Old   October 20, 2010, 06:53
Default
  #9
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
Yes it can. Your boundary conditions will depend on your version of OpenFOAM though.
eugene is offline   Reply With Quote

Old   October 21, 2010, 05:57
Default
  #10
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
Quote:
Originally Posted by eugene View Post
Yes it can. Your boundary conditions will depend on your version of OpenFOAM though.
Thanks Eugene,
i need a solver to simulate a room ventilation with inlet and outlet. the most similar solver i found was Heat Transfer solvers like buoyantBousinesqPimpleFoam. but i don't know why it doesn't understand the flow temperature for example i run it for a room with inlet fan with colder T than room but the flow effects on temperature was not real physically i think it set temperature of fan as a fixed wall instead of setting inlet patch flow temperature as desired T. Also velocity contour shows velocity to top in below of fan and it couldn't be true.
maysmech is offline   Reply With Quote

Old   October 21, 2010, 12:41
Default
  #11
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
The solver you used should give ok answers for the problem you described. Not sure why its not working.
eugene is offline   Reply With Quote

Old   October 21, 2010, 17:36
Question
  #12
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
Quote:
Originally Posted by eugene View Post
The solver you used should give ok answers for the problem you described. Not sure why its not working.

Here is geometry of the problem.
It is a room with 6 inlet fans at ceiling with T=290K and U=0.1 m/s. A window with 305 K is shown and walls and door have zeroGradient temperature.



After solving for 120 sec with time step of 0.1 the solution is:


--


and in real it is not acceptable to have T over all other temperatures. also contour of U_Z shows velocity to top instead of down of room in below of fans
I don't know where does the problem come form? the first thing i think is the solver may be not suitable for simulating flow with inlet and outlet with temperature like modeling of a room ventilation.
maysmech is offline   Reply With Quote

Old   October 22, 2010, 04:51
Default
  #13
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
Are you sure your fans are blowing into and not out of the room?
eugene is offline   Reply With Quote

Old   October 22, 2010, 05:57
Default
  #14
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
Yes, in boundary cells of fans u_z is -0.1 but below of them u_z is positive to up and shows a circulating flow to both sides of fan.
maysmech is offline   Reply With Quote

Old   October 23, 2010, 21:48
Default
  #15
New Member
 
Join Date: May 2009
Posts: 14
Rep Power: 9
Hasselhoff is on a distinguished road
Did you check your gravity settings? Tell us if you solved that problem, I was dealing with quite a similar problem some time ago.
Hasselhoff is offline   Reply With Quote

Old   October 24, 2010, 01:57
Default
  #16
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
----------

Last edited by maysmech; October 24, 2010 at 02:23.
maysmech is offline   Reply With Quote

Old   October 24, 2010, 02:22
Default
  #17
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
hi,
Gravity is -9.81 m/s2 as default.
when i increase the fans inlet velocity from 0.1 to 1 the velocity contour along z is more acceptable. i mean in u_z=-0.1 velocity near ceiling fans is to up but for u_z=-1 is to down and i think it is because of changing dominant between buoyancy (U to up) and inlet fan velocity (to down). I am not sure about this. Any more comment will be appreciated.
maysmech is offline   Reply With Quote

Old   October 24, 2010, 04:17
Default
  #18
New Member
 
Join Date: May 2009
Posts: 14
Rep Power: 9
Hasselhoff is on a distinguished road
In my case, the cold air from the radiators sinks down to the floor as it is supposed to be. Maybe the 0.1m/s is just too weak of a flow in your case. Or maybe your domain bounding box is scaled too big?
Good luck..
/M
Hasselhoff is offline   Reply With Quote

Old   October 24, 2010, 06:02
Default
  #19
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 10
maysmech is on a distinguished road
The other problem is how can we define T B.C for windows.

I think fixedValue is not a good assumption for windows temperature.
JuanRodriguez likes this.
maysmech is offline   Reply With Quote

Old   October 24, 2010, 08:45
Default
  #20
New Member
 
Join Date: May 2009
Posts: 14
Rep Power: 9
Hasselhoff is on a distinguished road
For windows I use heat flux boundaries. What kind of turbulence model are you using?
Hasselhoff is offline   Reply With Quote

Reply

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
Conductive Heat Transfer in OpenFoam? FloK OpenFOAM 8 March 18, 2011 14:51
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 06:25
Modified OpenFOAM Forum Structure and New Mailing-List pete Site News & Announcements 0 June 29, 2009 05:56
64bitrhel5 OF installation instructions mirko OpenFOAM Installation 2 August 12, 2008 18:07
Adventure of fisrst openfoam installation on Ubuntu 710 jussi OpenFOAM Installation 0 April 24, 2008 14:25


All times are GMT -4. The time now is 00:58.