velikaDjoka |
March 26, 2021 13:39 |
Pressure driven flow
2 Attachment(s)
Hello everybody,
I am trying to do a pressure derived flow in a quenching chamber working at 5 bars. With the following geometry.
Attachment 83539
The tubes are the inlets and the blue circle at the back is the outlet.
Obviously I want to have a flow perpendicular to the inlet patches . To assure this I wrote the boundary conditions for the velocity : pressureNormalInletOutletVelocity.
But how come my results of the velocity field are all parallel to the X axis ?
Attachment 83540
And how the U magnitude is not the same for the bottom inlets ? This might be because I am imposing p_rgh ( so there is some gravitational effect that enters. But then how can I impose a same inlet velocity for all the inlets with a pressure driven flow? ( I am using chtMultiRegionFoam, I don't know if it is important to mention )
Does any body know how to get the flow I want to have. That is a flow coming from the inlet patches ( perpendicularly ) and a flow leaving the chamber through the outlet. With a pressure driven flow.
Here are my boundary conditions are the following for p_rgh and U:
p_rgh
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0/vaccum";
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 1 -1 -2 0 0 0 0 ];
internalField uniform 500000;
boundaryField
{
outlet
{
type prghTotalPressure;
value uniform 499900;
p0 uniform 499900;
}
walls
{
type fixedFluxPressure;
value $internalField;
}
inletN
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
inletNE
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
inletSE
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
inletS
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
inletSW
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
inletNW
{
type prghTotalPressure;
value uniform 500050;
p0 uniform 500050;
}
vaccum_to_block
{
type fixedFluxPressure;
value $internalField;
}
}
// ************************************************************************* //
U
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0/vaccum";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 1 -1 0 0 0 0 ];
internalField uniform ( 0 0 0 );
boundaryField
{
outlet
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
walls
{
type noSlip;
value uniform ( 0.01 0 0 );
}
inletN
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
inletNE
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
inletSE
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
inletS
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
inletSW
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
inletNW
{
type pressureNormalInletOutletVelocity;
value uniform ( 0 0 0 );
}
vaccum_to_block
{
type noSlip;
value uniform ( 0 0 0 );
}
}
// ************************************************************************* //
|