[EshitaPal]

Hello,

I am using Openfoam 2.2 and trying to simulate natural circulation in a tank with multiple inlets at the bottom and multiple outlets at the top. As I run the transient buoyantBoussinesqPimpleFoam solver, the flow reverses very quickly, and the water starts flowing out of the inlet. Are my boundary conditions appropriate? Thanks.

U file

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{


pipe
{
type fixedValue;
value uniform (0 0 0);
}

inlet1
{
type pressureInletOutletVelocity;
value uniform (-0.1414 0 0);
}
inlet2
{
type pressureInletOutletVelocity;
value uniform (-0.130636565 -0.054111437 0);
}
inlet3
{
type pressureInletOutletVelocity;
value uniform (-0.099984898 -0.099984898 0);
}
outlet
{
type zeroGradient;
}
symmetry
{
type symmetryPlane;
}

wall
{
type fixedValue;
value uniform (0 0 0);
}


}



p_rgh file

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
pipe
{
type buoyantPressure;
rho rhok;
value uniform 0;
}
inlet1
{
type totalPressure;
p0 uniform 51.993; // this is equal to g*(height of the tank)
U U;
phi phi;
rho rhok;
psi none;
gamma 1;
value uniform 0;
}
inlet2
{
type totalPressure;
p0 uniform 51.993;
U U;
phi phi;
rho rhok;
psi none;
gamma 1;
value uniform 0;
}
inlet3
{
type totalPressure;
p0 uniform 51.993;
U U;
phi phi;
rho rhok;
psi none;
gamma 1;
value uniform 0;
}
outlet
{
type buoyantPressure;
rho rhok;
value uniform 0;
}
symmetry
{
type symmetryPlane;
}


wall
{
type buoyantPressure;
rho rhok;
value uniform 0;
}



}