Sedullo |
January 15, 2020 09:58 |
Natural convection boundary condition problem
1 Attachment(s)
I need to study the heat transfer on a wall by natural convection, I am simulating a problem where the central part of the right side wall is heated, the left side wall is set to slip condition. I am using buoyantBoussinesqPimpleFoam
As you can see in the picture, a kind of vortex is created, instead, I expected that the flux should have gone upward entering from below, I attach here the bc I set, I hope you can find a way to fix the problem, because I have been struggling for some days on it and I cannot go ahead.
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
sideWallheated
{
type fixedFluxPressure;
}
rightsideWallsNotheated
{
type fixedFluxPressure;
}
leftsideWallsNotheated
{
type fixedFluxPressure;
}
inlet
{
type totalPressure;
p0 uniform 0;
}
outlet
{
type totalPressure;
p0 uniform 0;
}
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
sideWallheated
{
type fixedValue;
value uniform 340;
}
rightsideWallsNotheated
{
type zeroGradient;
}
leftsideWallsNotheated
{
type zeroGradient;
}
inlet
{
type inletOutlet;
inletValue uniform 300;
value uniform 300;
}
outlet
{
type inletOutlet;
inletValue uniform 300;
value uniform 300;
}
#includeEtc "caseDicts/setConstraintTypes"
}
Code:
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
sideWallheated
{
type noSlip;
}
rightsideWallsNotheated
{
type noSlip;
}
leftsideWallsNotheated
{
type slip;
}
inlet
{
type outletInlet;
outletValue uniform (0 0 0);
value uniform (0 0 0);
}
outlet
{
type zeroGradient;
}
#includeEtc "caseDicts/setConstraintTypes"
}
|