CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Radiative heat transfer with natural convection inside a square cavity (http://www.cfd-online.com/Forums/openfoam/71654-radiative-heat-transfer-natural-convection-inside-square-cavity.html)

msarkar January 11, 2010 08:16

Radiative heat transfer with natural convection inside a square cavity
 
I am performing a 2D simulation of radiative heat transfer with natural convection inside a square cavity using buoyantSimpleRadiationFoam. The top and bottom walls are treated as adiabatic. The left wall is hot wall and right wall is cold wall. Inside air temperature is same as the cold wall temperature. All wall assumed to have same emissivity.

Mesh has been created using blockMesh. The patches are created using the following:

patches
(
wall topAndBottom
(
(3 7 6 2)
(1 5 4 0)
)
wall leftWall
(
(0 4 7 3)
)
wall rightWall
(
(2 6 5 1)
)
empty frontAndBack
(
(0 3 2 1)
(4 5 6 7)
)

Initial conditions for pressure, velocity and temperature are as follows:

Pressure:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 100000;

boundaryField
{
topAndBottom
{
type zeroGradient;
}

leftWall
{
type zeroGradient;
}

rightWall
{
type zeroGradient;
}

frontAndBack
{
type empty;
}
}

Velocity:
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 0 0);

boundaryField
{
floor
{
type fixedValue;
value uniform (0 0 0);
}

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

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

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

Temperature:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 323;

boundaryField
{
topAndBottom
{
type zeroGradient;
}

leftWall
{
type fixedValue;
value uniform 373.0;
}

rightWall
{
type fixedValue;
value uniform 323.0;
}

frontAndBack
{
type empty;
}
}

After running, it created uniform pressure inside the whole computational domain that is not correct. I am not able to find the error as I am new in OpenFOAM. Can anyone please help me to solve this problem?

I would appreciate any help...

msarkar January 11, 2010 23:21

Initial velocity condition I posted is wrong. Correct velocity condition is below:

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 0 0);

boundaryField
{
topAndBottom
{
type fixedValue;
value uniform (0 0 0);
}

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

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

frontAndBack
{
type empty;
}
}


All times are GMT -4. The time now is 22:10.