
[Sponsors] 
June 16, 2006, 07:02 
Hello everybody,
Did anybod

#1 
New Member
Alexander Jarosch
Join Date: Mar 2009
Location: Iceland
Posts: 6
Rep Power: 10 
Hello everybody,
Did anybody already modify the simpleFoam solver to include a body force term, i.e. gravity. I would like to add that and look for some hints. Thanks for all advice 

June 16, 2006, 08:27 
look at buoyantFoam, it is in

#2 
Senior Member
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 10 
look at buoyantFoam, it is in there
markus 

June 16, 2006, 09:02 
thanks, right I just look at t

#3 
New Member
Alexander Jarosch
Join Date: Mar 2009
Location: Iceland
Posts: 6
Rep Power: 10 
thanks, right I just look at that. I'm new to OpenFoam so it might take me a while to get into the code. buoyantFoam uses a variable density as I understand it, whereas I treat a constant density case.
I want to let the fluid move driven by gravity. Is it maybe simpler to just define the pressure gradient instead of including it in the solver? I guess this is possible in OpenFoam ;) 

June 16, 2006, 09:58 
twoLiquidMixingFoam is incompr

#4 
Senior Member
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 10 
twoLiquidMixingFoam is incompressible, try that


January 3, 2012, 10:11 

#5 
Member
Join Date: May 2010
Posts: 38
Rep Power: 9 
Is this the last word in this case?
Is it not possible to model steadystate, single phase, incompressible, laminar flow with gravity in OpenFoam? 

May 10, 2013, 03:27 

#6 
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 754
Rep Power: 10 

May 14, 2013, 01:03 

#7  
Senior Member

Quote:
Instead of variable density suppose I want constant density. Do you think the answers in both cases will differ? I too want to know the answer.. Last edited by Tushar@cfd; May 14, 2013 at 02:10. 

May 14, 2013, 11:44 

#8  
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 10 
Quote:
I guess you might have already traced out solution.... Anyways if your rho is fixed you just have to add it as source term in UEqn.h() with g as vector. u may define them in create Fields.h. So the efftect will be seen as source term in weight of fluid in only 1 momentum eq but rest would remaiin the same as before... Correct me if i am wrong ... Sachin 

May 15, 2013, 00:32 

#9  
Senior Member

Quote:
I have already done that thing with buoyantSimpleFoam... But later dropped the idea as it didn't suit my requirements and worked with Boussinesq approximation. If you see the code here variable density is needed for buoyancy.. I think I should say in detail, since I am working with buoyantBoussinesqSimpleFoam... I have solution but these differs from the numerical approach of the papers. For a very simple case the results matches exactly with the journal papers. But for the case of complex geometry results are showing deviations from paper. Which I don't think is permissible, with the use of higher order schemes and since the results are reaching convergence(I have set below 1e05 for all variable is my convergence criteria). When I happen to see the Isotherms and streamlines the dissipations seems to match with the paper but the gradT value is differs. It seems convective dissipation of heat in OF is more. Anyways, I am trying to figure it out. Is there anyone who has solved a complex case of natural convection (especially with buoyantBoussinesqSimpleFoam)? Last edited by Tushar@cfd; May 15, 2013 at 07:28. 

May 15, 2013, 00:56 

#10 
Senior Member

Well,
Anyone successful in solving natural convection of complex geometry (like heat transfer inside a room with a source object). If yes please, try to share your views.. Are you getting the exact results (as compared with numerical results of papers)? Thanks in Advance.. Last edited by Tushar@cfd; May 15, 2013 at 07:28. 

May 15, 2013, 09:13 

#11 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 10 
Sounds interesting ... Can you share the data ... what is the complexity...
CAn you share the difference ... Points you can check  1. What is effect of variation in Tref in boussinesq approximation 2. What is your Rayliegh number ... is it very high .. in which case it would turbulent flow 3. Changes in mesh Best Luck Sachin 

May 15, 2013, 23:32 

#12  
Senior Member

Quote:
I didn't understand this. What do you mean by variation in Tref in boussinesq approximation? Can you please explain, what kind of variation you are talking about? Condition 2: I already shared the data check my post... http://www.cfdonline.com/Forums/ope...implefoam.html Condition 3: Well, if I am right then mesh size won't have much impact on results in case of SIMPLE algorithm. Last edited by Tushar@cfd; May 16, 2013 at 00:49. 

May 16, 2013, 07:47 

#13 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 10 
I was refering to your ref temperature/ambient temperature...
I understand that these solvers might have some problems with complex geometries but the case described doesnt really seems to be a complex geometry ... For mesh, I dont know why it is not important. your concern from what understand is not something about no convergence but about dissipation of heat. So, if the mesh does capture the heated boundary condition, the results might improve. To check if the OF is giving some problem, you can compare it with some other s/w. With no real case results, we can keep guessing what you can try... If I could still thing of any more ideas, I would keep you posted. If you do solve your problem, let me know. Best Luck, Sachin 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
DPM Body Force  Sandilya Garimella  FLUENT  1  April 8, 2008 03:30 
Flow under gravity force  danvica  OpenFOAM Running, Solving & CFD  0  December 12, 2007 11:22 
Body force  Jacob  FLUENT  0  August 4, 2003 16:38 
how to treat the gravity force  Angela  Main CFD Forum  0  January 17, 2003 04:39 
UDFBody Force  Fluent_User  FLUENT  0  August 9, 2002 09:07 