CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Selecting numerical schemes for temperature equation (https://www.cfd-online.com/Forums/openfoam-pre-processing/62196-selecting-numerical-schemes-temperature-equation.html)

 maka May 21, 2006 13:35

I'm solving a case with transi

I'm solving a case with transient coupled (trough buoyancy term) momentum and energy equation (temperature eq, boussinesq assumption). I need some comment on my selection of fvSchemes. I need to apply central difference for convective and viscous terms together with Crank Nicholson for time.

http://www.opencfd.co.uk/openfoam/do...#x20-1050004.4

but I have two questions:

a) how does fluxRequired dictionary in fvSchemes affect the calculation? Should I put T there also? Note: I'm familiar with PISO except for the part of interpolation to surface.

b) Should I put a limiter on T?

you are welcome to give any comments on the following selections:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default CrankNicholson 1;
}

{
default none; //Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,T) Gauss linear; // energy equation. Gauss limitedLinear 1
}

laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian(1|A(U),p) Gauss linear corrected;
laplacian(alpha,T) Gauss linear corrected; // energy equation.
//alpha: thermal diffussivity.
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

{
default corrected;
}

fluxRequired
{
default no;
p;
T;
}

// ************************************************** *********************** //

best regards,
Maka.

 All times are GMT -4. The time now is 18:06.