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

2D differentially heated cavity, which is appropriate solver?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 3, 2010, 14:30
Default 2D differentially heated cavity, which is appropriate solver?
  #1
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 7
vetnav is on a distinguished road
Hello folks,

I am learning OpenFOAM and want to simulate a 2D differentially heated cavity problem for Rayleigh numbers in the range of 10^3 to 10^6, so the flow can be considered laminar for sure but how do I know if it will steady or unsteady? I am unable to decide which solver to use among the three

1) IcoFoam : but I found that the application file icoFoam.c solves only momentum equation so this is not appropriate solver.

2) buoyantBoussinesqPisoFoam : Assuming the flow will be unsteady (even if the flow attains steady state with this solver we observe that after certain number of time steps there will not be much change in the results) this seems appropriate solver but this is written for turbulent case so can I simply comment out the k, epsilon and R parts in fvSolution file to get the results for laminar case?

3) buoyantBoussinesqSimpleFoam : Assuming the flow will reach steady state (in case the flow turns out to be unsteady this is wring choice) this seems appropriate but again this is written for turbulent case.

So, I feel that I should go for option 2 (can I comment out the turbulent part and go ahead??), please tell me your suggestions/opinions.

Thanks a lot
vetnav is offline   Reply With Quote

Old   December 4, 2010, 03:56
Default
  #2
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,

You can very easily switch off turbulence in any solver. You have to mention this in file "$constant/turbulenceProperties". More details can be found in User Guide in Section 7.

As for your specific case, the choice of correct solver can only be advised, if you provide some more details about your case.
nakul is offline   Reply With Quote

Old   December 4, 2010, 15:02
Default
  #3
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 7
vetnav is on a distinguished road
Thanks for the reply Nakul, I will definitely try this. I am also considering to add energey equation to icoFoam (I have found a tutorial by Jassac where he shows how to add a scalar equation in icoFoam, so I am going to follow the same steps to add energy equation) and then run the same problem.

Regarding my case, it is just a two dimensional square cavity with differential heating, top and bottom walls are insulated, left wall is maintained at Th and the right wall is maintained at Tc where Th > Tc. Initially the fluid is at rest, so when this differential heating is applied the motion takes only due to buoyancy effects.

Can you now suggest me how do I decide whether I have to consider a Transient or steady state solver.

Thank you again
vetnav is offline   Reply With Quote

Old   December 5, 2010, 07:41
Default
  #4
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,
I would advice you to use a transient state solver because then you would be able to view many useful changes that happen as the flow tends to steady state.

Although transient solvers take longer time to solve but for your simple case time penalty won't be much.

Since gravity would also play an important role in your solution, hence I think that icoFoam may not be a suitable choice. But I am not sure about this point. It would be better if you first seriously go through icoFoam's code and determine whether it is a suitable solver or not.
nakul is offline   Reply With Quote

Old   December 14, 2010, 14:57
Default
  #5
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 7
vetnav is on a distinguished road
Hello nakul,

So I selected buoynatBoussinesqPimpleFoam for this problem, and in constant/RASProperties I switched the turbulence off, the problem is even though I change the temperature difference between the walls I don't see changes in the flow field.

The fluid properties are in transportProperties. I considered a Courant number of 0.1, and the velocity scale is taken as ratio of thermal diffusivity to characteristic length (0.1 m), I did the calculations for Ra = 10^4 and 10^5, and the flow field appears same for both the values at the end of 3000 time steps.

I am attaching the case, so please suggest me what can go wrong.

Thank you
Attached Files
File Type: gz 2Ddifferentialcavity.tar.gz (3.0 KB, 30 views)
vetnav is offline   Reply With Quote

Old   December 14, 2010, 20:35
Default
  #6
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by vetnav View Post
Hello nakul,

So I selected buoynatBoussinesqPimpleFoam for this problem, and in constant/RASProperties I switched the turbulence off, the problem is even though I change the temperature difference between the walls I don't see changes in the flow field.

The fluid properties are in transportProperties. I considered a Courant number of 0.1, and the velocity scale is taken as ratio of thermal diffusivity to characteristic length (0.1 m), I did the calculations for Ra = 10^4 and 10^5, and the flow field appears same for both the values at the end of 3000 time steps.

I am attaching the case, so please suggest me what can go wrong.

Thank you

Hi, dear vetnav, cheering.

in the private message you sent to me, you said that you tried two cases Rayleighnumber 10^4 (temperature difference between the walls ~ 4 degrees) and 10^5 ( temperature difference between the walls ~ 45 degrees). Did you calculate the cases with only different Rayleigh numbers, or with both different Rayleigh numbers and temperatures?

If you calculated them with only different Rayleigh numbers, I guess maybe the results are not large different, for the reason that both them are under or over the critical limit. You know, Rayleigh number is a parameter that determines when convection may occur, and convection occurs when the Rayleigh number exceeds a critical limit.

Best regards,
Chiven
chiven is offline   Reply With Quote

Old   December 14, 2010, 20:43
Default
  #7
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 7
vetnav is on a distinguished road
Hello Chiven,

Thank you for the reply, I used different temperature differences also, as shown in the paper by De Vahl Davis titled "Natural convection of air in a square cavity: A benchmark numerical solution", Intl. J. for Nume. methods in fluids, 3, 249 (1983), if everything goes fine I should see the different patterns for stream function (figure 3 in the paper), but with my OpenFOAM case I can not see any differences.

Reagrds
vetnav is offline   Reply With Quote

Old   December 14, 2010, 21:08
Default
  #8
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by vetnav View Post
Hello Chiven,

Thank you for the reply, I used different temperature differences also, as shown in the paper by De Vahl Davis titled "Natural convection of air in a square cavity: A benchmark numerical solution", Intl. J. for Nume. methods in fluids, 3, 249 (1983), if everything goes fine I should see the different patterns for stream function (figure 3 in the paper), but with my OpenFOAM case I can not see any differences.

Reagrds
Hi, vetnav,

I look through your case and find nothing wrong this time. would you please send the paper to me?

Best regards,
Chiven
chiven is offline   Reply With Quote

Old   December 14, 2010, 21:15
Default
  #9
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 7
vetnav is on a distinguished road
Thank you for your time Chiven, the pdf file size is 1.3 MB so I can not attach it here, but I am sending the figure from the paper which shows the stream function.

I can send you the paper if you can give me your email id.

Thank you
Attached Files
File Type: doc DeVahlDavis.doc (76.0 KB, 48 views)
vetnav is offline   Reply With Quote

Old   December 14, 2010, 21:22
Default
  #10
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by vetnav View Post
Thank you for your time Chiven, the pdf file size is 1.3 MB so I can not attach it here, but I am sending the figure from the paper which shows the stream function.

I can send you the paper if you can give me your email id.

Thank you

Hi, vetnav, you are welcome anytime.

I am not sure I can help you in the end, but I shall try the best. I used buoyantBoussinesq-Foam last year, and I am still working in OF. My email is chiven77@hotmail.com

Best regards,
Chiven
chiven is offline   Reply With Quote

Old   December 22, 2010, 08:32
Default
  #11
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi Naveen,

Can you please also send me the paper. I am not sure, but if I get time I will try to model the case and would inform you of results.

My email id is nakulthemaster@gmail.com
nakul is offline   Reply With Quote

Old   May 26, 2014, 12:06
Default
  #12
Senior Member
 
Srivathsan N
Join Date: Jan 2013
Location: India
Posts: 101
Rep Power: 4
Sherlock_1812 is on a distinguished road
Hi all

I know this is a very old post. But could you mail the paper to me too? DeVahl Davis? My ID is: srivathsan.iitm@gmail.com

Many thanks!
__________________
Regards,

Srivaths
Sherlock_1812 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
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
cavity in flat plate and drag prediction Far FLUENT 0 May 19, 2010 14:47
drag of flat plate with cavity Far FLUENT 0 May 18, 2010 04:57
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
Error during Solver cfd guy CFX 4 May 8, 2001 06:04


All times are GMT -4. The time now is 23:32.