Boundary condition setting for adding gravity in simpleFoam
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;
2 I make change to UEqn.H
eqnResidual = solve
UEqn() == -fvc::grad(p)+g
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;
value uniform 0;
I would suggest you to look into solvers where gravity is implemented, such as:
multiphase and heat transfer solvers.
Thanks for your kind attention. I have looked into the heat transfer solvers and I use buoyantSimpleFoam now it works
|All times are GMT -4. The time now is 10:42.|