CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   fvDOM validation with analytical results (https://www.cfd-online.com/Forums/openfoam-solving/216315-fvdom-validation-analytical-results.html)

MuriloHenriqueMoreira April 3, 2019 07:16

fvDOM validation with analytical results
 
Dear all,
I am trying to reproduce the work of Amouzandeh et al. He uses analytical results to benchmark both P1 and fvDOM. I am in doubt of how to set up the fvDOM case for infinite parallel black walls separated by a gray media. I have a square domain of 1x1x0.1 with a discretization of 40x40x1. For the frontaAndBack faces I have an empty boundary condition. For top and bottom I have a symmetry plane and as reported on the paper, I am using greyDiffusiveRadiation for the IDefault file as follows:


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

internalField uniform 0;

boundaryField
{

top
{
type symmetryPlane;
}

bottom
{
type symmetryPlane;
}

right
{
type greyDiffusiveRadiation;
value uniform 400;
}

left
{
type greyDiffusiveRadiation;
value uniform 600;
}

frontAndBack
{
type empty;
}
}

The problem is that my temperature calculated by buoyantSimpleFoam is converging to zero.

Does anybody knows how should I set this up?
For the P1 case all I had to do was to set G as:

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

internalField uniform 0;

boundaryField
{
top
{
type symmetryPlane;
}

bottom
{
type symmetryPlane;
}

right
{
type MarshakRadiationFixedTemperature;
Trad uniform 400; // radiation temperature field
value uniform 0; // place holder
}

left
{
type MarshakRadiationFixedTemperature;
Trad uniform 600; // radiation temperature field
value uniform 0; // place holder
}

frontAndBack
{
type empty;
}
}

I am felling that I am losing something important
Thank you very much!
All the best,
Murilo

mabinty April 23, 2019 15:26

2 Attachment(s)
Hallo Murilo,


please find attached the case set-ups of the mentioned "validation" cases.




Cheers,
Aram

MuriloHenriqueMoreira April 23, 2019 16:07

Thenk you very much
 
Dear Aram,
Thank you very much!
It worked!
All the best,
Murilo


All times are GMT -4. The time now is 21:16.