
[Sponsors] 
December 3, 2010, 14:30 
2D differentially heated cavity, which is appropriate solver?

#1 
Member
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 9 
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 

December 4, 2010, 03:56 

#2 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 8 
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. 

December 4, 2010, 15:02 

#3 
Member
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 9 
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 

December 5, 2010, 07:41 

#4 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 8 
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. 

December 14, 2010, 14:57 

#5 
Member
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 9 
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 

December 14, 2010, 20:35 

#6  
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 9 
Quote:
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 

December 14, 2010, 20:43 

#7 
Member
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 9 
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 

December 14, 2010, 21:08 

#8  
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 9 
Quote:
I look through your case and find nothing wrong this time. would you please send the paper to me? Best regards, Chiven 

December 14, 2010, 21:15 

#9 
Member
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 9 
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 

December 14, 2010, 21:22 

#10  
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 9 
Quote:
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 buoyantBoussinesqFoam last year, and I am still working in OF. My email is chiven77@hotmail.com Best regards, Chiven 

December 22, 2010, 08:32 

#11 
Senior Member
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 8 
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 

May 26, 2014, 12:06 

#12 
Senior Member
Srivathsan N
Join Date: Jan 2013
Location: India
Posts: 101
Rep Power: 6 
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 

Thread Tools  
Display Modes  


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 