CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

turbDyMFoam settings

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 18, 2010, 10:52
Default turbDyMFoam settings
  #1
New Member
 
Michael Buchmayr
Join Date: Mar 2010
Posts: 16
Rep Power: 7
MichiB is on a distinguished road
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?
MichiB is offline   Reply With Quote

Old   March 18, 2010, 10:54
Default
  #2
New Member
 
Michael Buchmayr
Join Date: Mar 2010
Posts: 16
Rep Power: 7
MichiB is on a distinguished road
Forgot the solution schemes

Quote:
solvers {
p PCG {
preconditioner DIC;
tolerance 1e-05;
relTol 0.0002;
minIter 2; maxIter 5; };
pFinal PCG {
preconditioner DIC;
tolerance 1e-05;
relTol 0.0001;
minIter 2; maxIter 5; };
U BiCGStab {
preconditioner { type DILU; }
minIter 1; maxIter 10;
tolerance 1e-07;
relTol 0; };
k BiCGStab { // and epsilon
preconditioner { type DILU; }
minIter 1; maxIter 5;
tolerance 1e-07;
relTol 0; };
epsilon BiCGStab {
preconditioner { type DILU; }
minIter 1; maxIter 5;
tolerance 1e-07;
relTol 0; };
omega BiCGStab {
preconditioner { type DILU; }
minIter 1; maxIter 5;
tolerance 1e-07;
relTol 0; };
R GAMG
{
preconditioner GAMG;
tolerance 1e-06;
relTol 0.1;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
smoother GaussSeidel;
minIter 0;
maxIter 100;
cacheAgglomeration true;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
};
nuTilda GAMG
{
preconditioner GAMG;
tolerance 1e-06;
relTol 0.1;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
smoother GaussSeidel;
minIter 0;
maxIter 100;
cacheAgglomeration true;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
};
}
SIMPLE
{
nNonOrthogonalCorrectors 2;
}
PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 2;
nGammaCorr 1;
nGammaSubCycles 4;
cGamma 1;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
omega 0.7;
R 0.7;
nuTilda 0.7;
}

// ************************************************** *********************** //
MichiB is offline   Reply With Quote

Old   August 4, 2010, 08:14
Default
  #3
Member
 
Nick Gardiner
Join Date: Apr 2009
Location: Chichester, UK
Posts: 93
Rep Power: 8
NickG is on a distinguished road
Hi

Try running turbFoam on it first to get reasonable values for k and epsilon.

I have momentumPredictor yes; in my PISO loop. Don't know if that's going to make a difference but the initial turbulence model variables certainly affect the stability.

Good luck
NickG is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
BC settings to expand pressure on atmosphere - simpleFoam / totalPressure sErik OpenFOAM Running, Solving & CFD 1 June 15, 2011 02:49
update turbDyMFoam antonio_ing OpenFOAM Installation 12 June 2, 2010 08:16
Export settings in the material panel sebastian FLUENT 1 November 18, 2009 06:11
Default settings for 3D boundary layer-Gambit 2.2 Ramine FLUENT 1 July 31, 2006 12:36
help on memory settings in TASCflow Dr C.Bhasker CFX 1 September 5, 2001 05:57


All times are GMT -4. The time now is 07:53.