CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

Boundary condition for natural convection

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 26, 2012, 02:28
Default Boundary condition for natural convection
  #1
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
Hey

I am new to openfoam, I want to simulate some tubes inside a cylinder having steam, air is entering from bottom of the cylinder and leaving from the top due to buoyancy. I am using the zeroGradient BC for the inlet(bottom) and the outlet(top). But the solution is diverging, can anybody help me what boundary condition should i give for pressure and velocity in this case.

Thank you
Ank is offline   Reply With Quote

Old   July 26, 2012, 02:36
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
What boundary condition did you give? Which solver do you use? What is diverging exactly? Please provide more information, because there are so many option. If I remember correctly, for p_rgh you have a buoyantPressure bc that may suit your needs.
Bernhard is offline   Reply With Quote

Old   July 26, 2012, 03:13
Default
  #3
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
I have used buoyantBoussinesqPisoFoam with Quick scheme, actually I have used uniform 0 for pressure at the inlet and the outlet, zeroGradient for velocity at the inlet and the outlet. After that the courant number is going very high, final residual for pressure is also very high after 3rd iteration only, hence the solution is getting diverged.I am having openfoam 1.6

Last edited by Ank; July 26, 2012 at 03:13. Reason: additional information
Ank is offline   Reply With Quote

Old   July 26, 2012, 17:21
Default
  #4
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
you trap pressure in your field so it is normal to diverge.
and zero gradient for both side is not good idea.
can you show your test case clearly?
niaz is offline   Reply With Quote

Old   July 27, 2012, 01:09
Default
  #5
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
hey thanks for the reply people,

My case is to do simulation of one cylinder having some 20 tubes inside, so I am simulating a quarter of the cylinder having 5 tubes. Steam is flowing from top to bottom at 413 K inside the tubes, and cylinder is open from both the side so that there can be a natural flow of air due to buoyancy.

Thanks..
Ank is offline   Reply With Quote

Old   July 27, 2012, 04:00
Default
  #6
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
Dear Ank
1- zeroGradient affects the solution of test cases. is your top and bottom limited?
I mean " can you increase the length of the main cylinder?"
2- change one of your pressure to zeroGradient.
tell me: does it work properly?
niaz is offline   Reply With Quote

Old   July 27, 2012, 05:01
Default
  #7
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
Hey Niaz

Thank you for the reply, I did what you said to me but it is stil now working, I will try some other BC as well..if you know some other Bc then please let me know
Ank is offline   Reply With Quote

Old   July 27, 2012, 05:36
Default
  #8
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
add these lines to your controldict
adjustTimeStep yes;

maxCo 0.2;
you may use a big times step. it helps to control courant number. and tiemstep.
try and tell me what is the result.
niaz is offline   Reply With Quote

Old   July 27, 2012, 06:00
Default
  #9
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
I did it, but courant number has already gone to 4, after 4 iterations...so again there is no result. I dont know what is the problem with the pressure..it is fluctuating now..I will let you know if it gets stable..thanks..
Ank is offline   Reply With Quote

Old   July 27, 2012, 06:03
Default
  #10
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Did you try decreasing your initial time-step? Please post a log-file, it will contains a lot of information.
Bernhard is offline   Reply With Quote

Old   July 27, 2012, 06:16
Default
  #11
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
I am using 1e-06 as initial time step, I am posting my case files, please have a look at it..solution has diverged again.
Attached Files
File Type: zip casefiles.zip (7.1 KB, 7 views)
Ank is offline   Reply With Quote

Old   July 27, 2012, 07:38
Default
  #12
Member
 
Primoz Ternik
Join Date: Apr 2009
Location: Maribor, Slovenia
Posts: 65
Rep Power: 7
ternik is on a distinguished road
Quote:
Originally Posted by Ank View Post
Hey

I am new to openfoam, I want to simulate some tubes inside a cylinder having steam, air is entering from bottom of the cylinder and leaving from the top due to buoyancy. I am using the zeroGradient BC for the inlet(bottom) and the outlet(top). But the solution is diverging, can anybody help me what boundary condition should i give for pressure and velocity in this case.

Thank you
Dear all,

I am having some difficulties with the natural convection as well;details were posted here:

buoyantBoussinesqSimpleFoam

After quite some time spent I have concluded that the problem might be with the boundary conditions for p and p_rgh (zero velocity and zero gradient/constant temperature on solid walls are pretty straightforward).

It is true that I have found some (similar) posts regarding the natural convection problem, but no final (and validated) results are available.

Is there someone that can give me usefull hints how to set the correct pressure boundary conditions (not only from numerical but from physical point of veiw as well)?!

Best,
Primoz
ternik is offline   Reply With Quote

Old   July 27, 2012, 10:46
Default
  #13
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
Quote:
Originally Posted by Ank View Post
I am using 1e-06 as initial time step, I am posting my case files, please have a look at it..solution has diverged again.
Dear Ank
at first, change your pressure BC. it is not correct to use uniform 0 for inlet and outlet. this BC means that no flow can in or out and it forces to circulate in the tube.
secondly, I think courant control is not work in of1.6 which you used
niaz is offline   Reply With Quote

Old   July 28, 2012, 08:59
Default
  #14
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
Hey Niaz

I tried other BCs as well, but didnt work out. Can you suggest me some specific BC for this kind of case, which you have used?
Thanks
Ank is offline   Reply With Quote

Old   July 28, 2012, 15:33
Default
  #15
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
Dear Ank
are you check that your case is not turbulent?
niaz is offline   Reply With Quote

Old   July 30, 2012, 01:40
Default
  #16
Ank
Member
 
ankur
Join Date: May 2012
Location: India
Posts: 50
Rep Power: 4
Ank is on a distinguished road
Hey Niaz,
If the flow is turbulent then what difference does it make to the boundary conditions?
Ank is offline   Reply With Quote

Old   July 30, 2012, 05:26
Default
  #17
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
dear ank
firstly, you have a problem in convergency of your problem.
it can affect the solution.
niaz 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
Convection Boundary condition tomcatbobby FLUENT 2 April 30, 2012 14:50
Domain Imbalance HMR CFX 3 March 6, 2011 21:10
asking for Boundary condition in FLUENT Destry FLUENT 0 July 27, 2010 01:55
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
How to set boundary condition in Fluent for the fo Peiyong FLUENT 1 November 10, 2006 12:44


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