MichiB |
March 18, 2010 09:52 |
turbDyMFoam settings
Hi everybody!
we are trying to calculate the transient turbulent flow through a pump, yet with little success. We where using turbDymFoam, with a uniform (p and U) initial calculation domain. The calculation somehow diverges, and we can't figure out why.
The fvSchemes we set look as follows:
Quote:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind Gauss;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,omega) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DomegaEff,omega) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
laplacian(rAU,p) Gauss linear corrected;
laplacian(rAU,pcorr) Gauss linear corrected;
laplacian(nu,U) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
interpolate(HbyA) linear;
interpolate(1|A) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
pcorr;
}
// ************************************************** *********************** //
|
We now wanted to initialize the field with a steady state solution using MRFSimpleFoam. But again -> ocillating residuals.
Does anybody have an idea what we could have done wrong?
|