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

Error in solver? (buoyantPisoFoam)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2010, 18:51
Default Error in solver? (buoyantPisoFoam)
  #1
Member
 
James Baker
Join Date: Dec 2009
Posts: 35
Rep Power: 7
fijinx is on a distinguished road
This is for you smart ones out there (or at least smarter than me)! I have a simplified problem where I have a 1 meter cubed box with an inflow of 1kg/sec (trying to model filling a tank with gas). after 15 seconds I check the values, and the pressure is very much off (factor of 12ish). I have check, re-checked, and re-re-checked my values in all the files, but they all look accurate to me. Based off hand calcs, I'm expecting a pressure of around 1.724 MPa after 15 seconds, but it only gets to 135 kPa. Anyone out there that can crack this? I am using buoyantPisoFoam to solve the problem.

http://rapidshare.com/files/343125087/meterBox.zip.html
fijinx is offline   Reply With Quote

Old   January 29, 2010, 18:56
Default
  #2
Member
 
James Baker
Join Date: Dec 2009
Posts: 35
Rep Power: 7
fijinx is on a distinguished road
The link above is the full case file.
fijinx is offline   Reply With Quote

Old   February 1, 2010, 14:27
Default
  #3
Member
 
James Baker
Join Date: Dec 2009
Posts: 35
Rep Power: 7
fijinx is on a distinguished road
Anyone have a chance to look at this yet? I'm still pulling my hair out.
fijinx is offline   Reply With Quote

Old   February 2, 2010, 03:43
Default
  #4
Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 92
Rep Power: 8
matejfor is on a distinguished road
To save some of your boldness for other cases, just without much compressible experience and source checking, I would say, that this buoyant solver, although compressible, is only weakly compressible. It means it does compute the buoyancy force from density difference and not the Boussinesq approximation, but is not suited for real compressible case like yours.

I would suggest to try rhoPisoFoam or some other real compressible solver.

I hope someone will correct me if I'm wrong.

Good luck

matej
matejfor is offline   Reply With Quote

Old   February 2, 2010, 07:00
Default
  #5
Member
 
Join Date: Dec 2009
Posts: 57
Rep Power: 7
Graham81 is on a distinguished road
Hey James,

Since Im waiting on a simulation, Im trying to think with you....though Im rather newbish to openFOAM. I downloaded your case, but was wondering where you prescribed the inflow of gass.

Pieter
Graham81 is offline   Reply With Quote

Old   February 2, 2010, 07:05
Default
  #6
Member
 
Join Date: Dec 2009
Posts: 57
Rep Power: 7
Graham81 is on a distinguished road
My bad, I see it now in the BC for U.
Graham81 is offline   Reply With Quote

Old   February 2, 2010, 15:19
Default
  #7
Member
 
James Baker
Join Date: Dec 2009
Posts: 35
Rep Power: 7
fijinx is on a distinguished road
Thank you matejfor! You were correct, using rhoPisoFoam gave me almost exact results to my hand calculations. I do however have one delima now.

Does rhoPisoFoam support gravity, and how would it be implimented?

--James
fijinx is offline   Reply With Quote

Old   February 2, 2010, 15:58
Default
  #8
Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 92
Rep Power: 8
matejfor is on a distinguished road
Good to hear!

now the gravitation. Before turning most of my brain on, do you really need it?
You are filling a volume quite fast and compress it a lot. Does the gravity really plays important role in your analysis.
Is yes, then let's have a look how the gravity is implemented.
Let's have a look into: $FOAM/applications/solvers/heatTransfer/buoyantPisoFoam/buoyantPisoFoam.C

line 49 is important ( in my 1.6.x installation) with
#include "readGravitationalAcceleration.H"

then the gravitational force is taken into account as a source of momentum in UEqn.H
and pEqn.H

especially with a pressure equation it is not going to be copy&paste as I bet the pressure equation looks different in the fully compressible case. But with a pen and pencil, DoNotDisturb sign and mug of tea, you should found out how to do it.


hope this helps

matej
matejfor 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
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 08:45
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25


All times are GMT -4. The time now is 13:56.