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

heat transfer with OpenFoam

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   March 21, 2011, 03:43
Default heat transfer with OpenFoam
  #1
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Dear all,

I want to model a turbulent flow with heat transfer caused by forced convection and radiation.
How can I model that with OF. I have checked the tutorial files of the heat transfer solvers. I every solver-file there are the following files: alphatwallfunction, epsilon, k, nut, p, p_rgh, T, T.org, U.
- What is the meaning of them (T=temperature, U=velocity.)?
- Is it possible to enter a heat transfer coefficient calculated by my own or does OF calculate everything by itself?

Thank you for your answer.

Best Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Old   March 21, 2011, 17:54
Default
  #2
New Member
 
Fatih
Join Date: Sep 2010
Location: Hamburg
Posts: 12
Rep Power: 7
ang_dipl_ing_fd is on a distinguished road
hi,
may be this will help you:
-alphatWallFunction=turbulent thermal diffusity (or eddy thermal diffusity)
-k=turbulent energy
-epsilon=dissipation
-nut=turbulent viscosity
-p=pressure, exact: pressure/rho for inkompressible Solver
-p_rgh=hydrostatic pressure p_0+rho*g*h (same as above p_0/rho+g*h); if g=(0 0 0) than p=p_rgh
-T=temperature [K]
-T.org=backup for T
-U=velocity (vector)
and the last, which you forget
-kappat=turbulent thermal conductivity

i think you define the heat transfer coefficient over the Pr in transportProperties. but i'm not sure.
palmerlee and mgg like this.
ang_dipl_ing_fd is offline   Reply With Quote

Old   March 22, 2011, 03:03
Default
  #3
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Thank you for your answer!

I have now checked the buoyantSimpleRadiationFoam-solver. There are some more files.
Do you know the meaning of them, too?

G
mut

Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 03:48
Default
  #4
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
I have tried to use the buoyantSimpleRadiationFoam-solver for an easy model.
The folliwing error appears:

Create time
Create mesh for time = 0

Reading g
Reading thermophysical properties
Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThe rmo<hConstThermo<perfectGas>>>>>
#0 Foam::error:rintStack(Foam::Ostream&) in "/data/caehgb.za/studienarbeit-di rk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOA M.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/data/caehgb.za/studienarbeit-dirk/work /OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::calculate() in "/data/caehgb.za/stu dienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPO pt/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7 .1/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ur eMixture<Foam::constTransport<Foam::specieThermo<F oam::hConstThermo<Foam:erfec tGas> > > > > >::New(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk /work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libbasicTher mophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/data/caehgb.za/studienar beit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/lib basicThermophysicalModels.so"
#7 main in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/ OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/buoyantSimpleRadiationFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream: :versionNumber, Foam::IOstream::compressionType) const in "/data/caehgb.za/studi enarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/applications/bin/li nux64GccDPOpt/buoyantSimpleRadiationFoam"
Floating point exception

Can someone help me?
Maybe OF isn't installed correct?

Regards
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 03:53
Default
  #5
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Quote:
Originally Posted by ang_dipl_ing_fd View Post
hi,
may be this will help you:
-alphatWallFunction=turbulent thermal diffusity (or eddy thermal diffusity)
-k=turbulent energy
-epsilon=dissipation
-nut=turbulent viscosity
-p=pressure, exact: pressure/rho for inkompressible Solver
-p_rgh=hydrostatic pressure p_0+rho*g*h (same as above p_0/rho+g*h); if g=(0 0 0) than p=p_rgh
-T=temperature [K]
-T.org=backup for T
-U=velocity (vector)
and the last, which you forget
-kappat=turbulent thermal conductivity

i think you define the heat transfer coefficient over the Pr in transportProperties. but i'm not sure.
Hi ang_dipl_ing_fd,

do you also know the following:
In the thermophysical properties:
mixture air 1 28.9 1000 0 1.8e-05 0.7;

What does that mean?
Does every number have a meaning?

Best Regards
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 04:38
Default
  #6
Member
 
Alpesh Vora
Join Date: Jan 2011
Location: Cottbus, Germany
Posts: 47
Rep Power: 6
alvora is on a distinguished road
Hi

In the thermophysical properties:
mixture air 1 28.9 1000 0 1.8e-05 0.7;

The thermalPhsical Properties is calculated based which model you specified, like
if model is hPsiThermo, calculation based on enthalpy 'h' and compressibility 'psi'.
and if model is hRhoThermo, calculation based on enthalpy or sensible enthalpy.
and if model is pureMixture, calculation based on passive gas mixtures.

parameters :

air is mixture fluid name
1 is for pure Mixture (mixture property) (number of molecules of species)
28.9 for molecule weight (kg/kmol)
1000 for heat capacity at constant pressure [J/kmol K]
0 for reference enthalpy [J/kmol]
1.8e-05 for Dynamics viscosity (kg/ms)
0.7 is for Prandtl number.

you can change properties according to your problem (fluid).

best luck

Alpesh



alvora is offline   Reply With Quote

Old   March 22, 2011, 04:58
Default
  #7
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 6
sixwp is on a distinguished road
Quote:
Originally Posted by tH3f0rC3 View Post
Hi ang_dipl_ing_fd,

do you also know the following:
In the thermophysical properties:
mixture air 1 28.9 1000 0 1.8e-05 0.7;

What does that mean?
Does every number have a meaning?

Best Regards
Hi tH3f0rc3 (hard to say :-p)

As far as I know:

1 => Number of moles of this specie
28,9 => Molecular weight
1000 => Cp (Specific heat capacity)
0 => heat of fusion (useless, I guess in your case)
1.8e-05 => dynamic viscosity
0.7 => Prantl number

Quote:
Originally Posted by tH3f0rC3 View Post
Reading g
Reading thermophysical properties
Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThe rmo<hConstThermo<perfectGas>>>>>
#0 Foam::error:rintStack(Foam::Ostream&) in "/data/caehgb.za/studienarbeit-di rk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOA M.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/data/caehgb.za/studienarbeit-dirk/work /OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::calculate() in "/data/caehgb.za/stu dienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPO pt/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7 .1/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ur eMixture<Foam::constTransport<Foam::specieThermo<F oam::hConstThermo<Foam:erfec tGas> > > > > >::New(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk /work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libbasicTher mophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/data/caehgb.za/studienar beit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/lib basicThermophysicalModels.so"
#7 main in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/ OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/buoyantSimpleRadiationFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream: :versionNumber, Foam::IOstream::compressionType) const in "/data/caehgb.za/studi enarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/applications/bin/li nux64GccDPOpt/buoyantSimpleRadiationFoam"
Floating point exception
I'm also interested in this matter (I've got the same issue but with chtMultiRegionSimpleFoam). Someone clever than us? ^^

Cheers guys

[Edit]: I was too slow, apparently (damn coffee break ! :-D), sorry for this double answer
__________________
Sam
sixwp is offline   Reply With Quote

Old   March 22, 2011, 05:26
Default
  #8
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Hey,
thanks for your answer.
Do you also know, where I can type in the properties of the walls, for example that a wall consists of nickel?

Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 06:28
Default
  #9
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 6
sixwp is on a distinguished road
It depends what you want to do:

for a fixed temperature wall, you don't need to input any properties. (only create a radiationProperties file => check the tutorial hotRadiationRoom). In 0/T, fixe the wall temperature.

If you're really eager to use the wall properties, include it in the mesh and use a MultiRegion solver (much more complicated, for my level of knowledge obviously).

cheers




__________________
Sam
sixwp is offline   Reply With Quote

Old   March 22, 2011, 06:44
Default
  #10
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
OK, but I think i have to gibe the solver an information of the material which should be heaten up. For example the density or cp...
Where can I do that?

Best Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 06:47
Default
  #11
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 6
sixwp is on a distinguished road
I think I got your case wrong:
Can you explain briefly what your heat transfer consist in? And add your directories files also, if you don't mind.
__________________
Sam
sixwp is offline   Reply With Quote

Old   March 22, 2011, 06:58
Default
  #12
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
I want to simulate a heat transfer from a turbulent (I guess so) flow to a surface, that should be heaten up.
In addition to that there is also radiation from the hot wall to the surface.

If this description is too short, please let me know.

Which files do you want to see?

Regards,
tH3f0rC3

Last edited by tH3f0rC3; March 22, 2011 at 07:13.
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 07:18
Default
  #13
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Quote:
Originally Posted by sixwp View Post
Hi tH3f0rc3 (hard to say :-p)
I'm also interested in this matter (I've got the same issue but with chtMultiRegionSimpleFoam). Someone clever than us? ^^
Hi sixwp,

it is exactly the same issue just with other file dorectory?

Regards
tH3f0rC3 is offline   Reply With Quote

Old   March 22, 2011, 09:38
Default
  #14
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 6
sixwp is on a distinguished road
Quote:
Originally Posted by tH3f0rC3 View Post
I want to simulate a heat transfer from a turbulent (I guess so) flow to a surface, that should be heaten up.
In addition to that there is also radiation from the hot wall to the surface.

If this description is too short, please let me know.

Which files do you want to see?

Regards,
tH3f0rC3
I'm sorry, I must be tired, I don't really see what's your geometry's like. Can you show me a raw scheme?

Concerning the directories files: can you post your 0, constant and system folders? Maybe (not likely but, who knows), I'll be able to find out your problem and therefore, mine :-)
__________________
Sam
sixwp is offline   Reply With Quote

Old   March 22, 2011, 09:42
Default
  #15
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Yes, of course.
But you have to give me an email adress of yours, because I can't post fotos.

Regards
tH3f0rC3 is offline   Reply With Quote

Old   March 23, 2011, 06:57
Default
  #16
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Hi,

is there someone else, who understands the following error message?

Create time
Create mesh for time = 0

Reading g
Reading thermophysical properties
Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThe rmo<hConstThermo<perfectGas>>>>>
#0 Foam::error:rintStack(Foam::Ostream&) in "/data/caehgb.za/studienarbeit-di rk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOA M.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/data/caehgb.za/studienarbeit-dirk/work /OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::calculate() in "/data/caehgb.za/stu dienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPO pt/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<F oam::hConstThermo<Foam:erfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7 .1/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ur eMixture<Foam::constTransport<Foam::specieThermo<F oam::hConstThermo<Foam:erfec tGas> > > > > >::New(Foam::fvMesh const&) in "/data/caehgb.za/studienarbeit-dirk /work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libbasicTher mophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/data/caehgb.za/studienar beit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/lib basicThermophysicalModels.so"
#7 main in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/ OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/buoyantSimpleRadiationFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream: :versionNumber, Foam::IOstream::compressionType) const in "/data/caehgb.za/studi enarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/applications/bin/li nux64GccDPOpt/buoyantSimpleRadiationFoam"
Floating point exception
tH3f0rC3 is offline   Reply With Quote

Old   March 23, 2011, 08:40
Default
  #17
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
I know now, that this error message appears every time I want to involve the heat transfer.
Maybe this helps.
tH3f0rC3 is offline   Reply With Quote

Old   March 27, 2011, 15:57
Default
  #18
New Member
 
Fatih
Join Date: Sep 2010
Location: Hamburg
Posts: 12
Rep Power: 7
ang_dipl_ing_fd is on a distinguished road
hi tH3f0rC3,
did you run the tutorials without any change?
ang_dipl_ing_fd is offline   Reply With Quote

Old   March 29, 2011, 02:46
Default
  #19
Member
 
Join Date: Nov 2009
Location: Germany
Posts: 96
Rep Power: 7
val46 is on a distinguished road
Quote:
Originally Posted by tH3f0rC3 View Post
OK, but I think i have to gibe the solver an information of the material which should be heaten up. For example the density or cp...
Where can I do that?

Best Regards,
tH3f0rC3
Hi,

if I get you right you just have one fluid domain.
If you want to know the interaction of the fluid with a connected solid than you need a multiregion solver.
(density and cp are values connected to a volume and not to a surface)

Regards,
Toni
val46 is offline   Reply With Quote

Old   March 30, 2011, 02:24
Default
  #20
Senior Member
 
Join Date: Mar 2011
Posts: 155
Rep Power: 6
tH3f0rC3 is on a distinguished road
Hello,

I have another question to
-k=turbulent energy
-epsilon=dissipation
-nut=turbulent viscosity

I have to set up these files in the 0-file as boundary cinditions.
k and epsilon must have a value for the inlet surfaces. But how can I calculate these values?
And what is meant with the values
tH3f0rC3 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
Heat transfer problem seojaho CFX 6 May 6, 2010 00:32
Heat Transfer Coefficient los OpenFOAM Running, Solving & CFD 5 January 31, 2010 18:44
Which Heat transfer coeffcient to use? tengra FLUENT 1 May 1, 2009 13:49
No results for solid domain Gary Holland CFX 10 March 13, 2009 04:30
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


All times are GMT -4. The time now is 21:14.