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

Boundary condition setting for adding gravity in simpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 8, 2013, 07:51
Unhappy Boundary condition setting for adding gravity in simpleFoam
  #1
New Member
 
Bin Xu
Join Date: Apr 2012
Location: Singapore
Posts: 23
Rep Power: 5
norkistar is on a distinguished road
Hi Foamers,

I am just starting simulation water flow in a horizontal circular pipe. Without adding gravity, it works well. With adding gravity, it can't converge. The result diverse and become unbelievable. I think problem might due to the boundary condition setting for p. I will list the step that I did for adding the gravity. Your kindly attention would be much appreciated.

1 I add g in the createFields.H as follows
Info<< "Reading g\n" << endl;
uniformDimensionedVectorField g
(
IOobject
(
"g",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);

2 I make change to UEqn.H
eqnResidual = solve
(
UEqn() == -fvc::grad(p)+g
).initialResidual();

3 I make change to pEqn.H
U -= (fvc::grad(p)-g)/AU;

After that, I wmake to make these changes take effect.

In the constant folder, the g file is listed as follow
dimensions [0 1 -2 0 0 0 0];
value ( 0 0 -9.81 );

In the 0 folder, the p is setting as previous simulation which is without gravity
dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
tubewall
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
plane
{
type symmetryPlane;
}
}
norkistar is offline   Reply With Quote

Old   February 14, 2013, 12:59
Default
  #2
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 129
Rep Power: 8
ziemowitzima is on a distinguished road
I would suggest you to look into solvers where gravity is implemented, such as:
multiphase and heat transfer solvers.
ziemowitzima is offline   Reply With Quote

Old   February 15, 2013, 20:06
Default
  #3
New Member
 
Bin Xu
Join Date: Apr 2012
Location: Singapore
Posts: 23
Rep Power: 5
norkistar is on a distinguished road
Thanks for your kind attention. I have looked into the heat transfer solvers and I use buoyantSimpleFoam now it works
norkistar is offline   Reply With Quote

Reply

Tags
boundary condition, convergence, gravity, simplefoam

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
setting Boundary condition for additional variable fek66 CFX 2 May 19, 2011 04:44
UDF setting wall boundary condition with a DEFINE_PROFILE NLao FLUENT 1 April 17, 2011 05:23
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 02:53
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Warning 097- AB CD-adapco 6 November 15, 2004 05:41


All times are GMT -4. The time now is 07:28.