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

Heat Transfer possibility with DnsFoam ?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 5, 2010, 09:23
Default Heat Transfer possibility with DnsFoam ?
  #1
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Hi Foamers !

I am trying to simulate a very simple case : a turbulent flow into a duct.
I have to "see" the turbulent structures in the boundary layer (maybe not the hairpins, but ideally it would be nice) and to analyze heat transfer at this boundary layer.

I am running a case with DnsFoam : 1024*1024*2 mesh for a duct with 10 millimeter high. Inlet is initialized with a LES model for turbulences.
Normally I should have good results with that (still running at present) ^^

But now I wonder if there is a possibility, with DnsFoam, to run the same case and to "see" the heat transfer at the walls. I mean : is someone has already run a heat transfer case with DnsFoam, and how ?
Any kind of help would be greatly appreciated !
NorbertB is offline   Reply With Quote

Old   July 6, 2010, 09:42
Default
  #2
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Actually I got another problem with dnsFoam : it seems to be impossible to run a parallel computing case.
Kmesh always outputs a fatal error about the number of cells (I use a regular grid, but while running dnsFoam after decomposePar it gives me that error). I think I understand that error, but since the "-parallel" option is available for dnsFoam I wonder if there is a way to avoid that error.

Thank you in advance
NorbertB is offline   Reply With Quote

Old   July 6, 2010, 10:10
Default
  #3
Member
 
Francois Gallard
Join Date: Mar 2010
Location: Edinburgh
Posts: 39
Rep Power: 7
fgal is on a distinguished road
Hi Norbert,

Just a remark, it seems that you are using a 2D grid, well, with 2cells depth. I does not make sense to do DNS or LES in 2D because all the turbulent energy cascade is wrong in 2d as turbulence structures are fundamentally 3D.

For the heat transfer, you can have a look at this boundary conditions :
Conductive Heat Transfer in OpenFoam?

You may have made a mistake in the decomposeparDict file.

Cheers,

Francois
fgal is offline   Reply With Quote

Old   July 7, 2010, 04:38
Default
  #4
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Hi Francois,

I know a pure 2D DNS is impossible, because the physics is all wrong, but I thought using 2 cells in the z direction would be enough since I am not interested in the behavior in that direction. Am I wrong ?

Thank you for the link, I will check this out.

Does anyone have an experience with a parallel DNS case and could help me ?
My decomposePar file is the same as in the UserGuide, but I use a 15mm*10mm*1mm (1024*1024*2) box and not a cube, could it put any problem for decomposePar or Kmesh ?

Thank you in advance
NorbertB is offline   Reply With Quote

Old   July 7, 2010, 07:07
Default
  #5
Member
 
Francois Gallard
Join Date: Mar 2010
Location: Edinburgh
Posts: 39
Rep Power: 7
fgal is on a distinguished road
Hello Norbert,

Having in mind that turbulence is the presence in the flow chaotic eddies, existing on a wide range of length scales, and all these scales exchanging energy via the stretching and viscous phenomena, the conclusion is that the mesh must be able to hold (I didn't find a better word for that) all this range of spacial wave length. But with only 2 points in the z direction, you can't discretize any signal. Basically, it is the same thing as discretizing a sinus function with 2 points. The signal is completely destroyed. You would need around 10 points per period to discretize correctly the signal.

As a wide range of eddies exist, typically, when you have a look at the Kolmogorov spectrum you will see that often more than one decade of eddies exist in frequency. So a fast conclusion would be that you need at least 100 points in all directions, whatever your problem is. Of course it is more complicated than that. Actually, in papers, we most often see DNS meshes with the same number of cells in all directions (people, tell me if I am wrong), because of that completely 3D structure of the turbulence.

I suggest you to read these pages :
http://www.cfd-online.com/Wiki/Intro..._of_turbulence

http://www.cfd-online.com/Wiki/Direct_numerical_simulation_(DNS)

The solver that you use should not be linked with parallelization issues (again, tell me if I am wrong), because in OF solvers source codes, you will notice that there is no parallelization instruction as this is completely in charge of the OF low level classes. And this is a strength of OF because parallelization is not easy to program.

Cheers,

Francois.

Last edited by fgal; July 7, 2010 at 08:02.
fgal is offline   Reply With Quote

Old   July 7, 2010, 08:14
Default
  #6
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Hi Francois,

Indeed most cases in papers have the same number of cells in all directions, you may be right about that : I will enlarge my number of cells in the z direction.

However I red on this forum (do not find out the link ...) that a parallelization is not possible at present with dnsFoam.
Is there any expert who could come here and give us a definitive answer plz

Thanks for your time
NorbertB is offline   Reply With Quote

Old   July 7, 2010, 08:39
Default
  #7
Member
 
Francois Gallard
Join Date: Mar 2010
Location: Edinburgh
Posts: 39
Rep Power: 7
fgal is on a distinguished road
Hi,

You seem to be right, I get the same error when I launched the tutorial after decomposition of the domain. This is very surprising to me ! Also doing DNS without parallelism is a bit restrictive...

Cheers and good luck,

Francois.
fgal 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 coefficient - what is waht Stan FLUENT 20 March 11, 2014 11:42
Convective / Conductive Heat Transfer in Hypersonic flows enigma Main CFD Forum 2 November 1, 2009 23:53
How can I increase Heat Transfer at Domain Interf? B.Simon CFX 3 October 28, 2008 19:53
Question on heat transfer coefficient!!! Benny FLUENT 7 June 7, 2005 09:25
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


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