ENH: Updated tutorials

This commit is contained in:
andy
2011-04-15 16:46:35 +01:00
parent 717db62583
commit 154761a3a1
55 changed files with 369 additions and 314 deletions

View File

@ -17,21 +17,30 @@ FoamFile
solvers
{
p
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}
"(p|rho)Final"
{
$p;
tolerance 1e-06;
relTol 0;
}
rho
"(U|b|Su|Xi|h|hu|k)"
{
$p;
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
"(U|ft|fu|b|Xi|Su|h|hu|R|k|epsilon)"
"(U|b|Su|Xi|h|hu|k)Final"
{
solver PBiCG;
preconditioner DILU;
@ -40,12 +49,17 @@ solvers
}
}
PISO
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
"(Xi|Su)" 1;
}
// ************************************************************************* //

View File

@ -17,22 +17,30 @@ FoamFile
solvers
{
p
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
relTol 0.1;
}
rho
"(p|rho)Final"
{
$p;
tolerance 1e-05;
tolerance 1e-06;
relTol 0;
}
"(U|ft|fu|b|Xi|Su|h|hu|R|k|epsilon)"
"(U|b|Su|Xi|h|hu|k)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
"(U|b|Su|Xi|h|hu|k)Final"
{
solver PBiCG;
preconditioner DILU;
@ -41,7 +49,7 @@ solvers
}
}
PISO
PIMPLE
{
nOuterCorrectors 2;
nCorrectors 1;
@ -53,4 +61,5 @@ relaxationFactors
"(Xi|Su)" 1;
}
// ************************************************************************* //

View File

@ -17,22 +17,30 @@ FoamFile
solvers
{
p
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
relTol 0.1;
}
rho
"(p|rho)Final"
{
$p;
tolerance 1e-05;
tolerance 1e-06;
relTol 0;
}
"(U|ft|fu|b|Xi|Su|h|hu|R|k|epsilon)"
"(U|b|Su|Xi|h|hu|k|epsilon)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
"(U|b|Su|Xi|h|hu|k|epsilon)Final"
{
solver PBiCG;
preconditioner DILU;
@ -41,11 +49,11 @@ solvers
}
}
PISO
PIMPLE
{
nCorrectors 2;
nOuterCorrectors 2;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
momentumPredictor yes;
}

View File

@ -17,30 +17,46 @@ FoamFile
solvers
{
p
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 0.1;
}
"(p|rho)Final"
{
$p;
tolerance 1e-09;
relTol 0;
}
rho
{
$p;
tolerance 1e-06;
}
"(U|Yi|hs|k|epsilon)"
"(Yi|k|epsilon)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
"(U|hs)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0.1;
}
"(U|hs)Final"
{
$U;
tolerance 1e-06;
relTol 0;
}
}
PISO
PIMPLE
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;

View File

@ -22,17 +22,48 @@ solvers
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}
rhoFinal
{
$rho;
tolerance 1e-06;
relTol 0;
}
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0.1;
}
pFinal
{
$p;
tolerance 1e-6;
relTol 0.0;
}
"(U|Yi|hs|k|epsilon)"
"(U|hs)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
"(U|hs)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
"(Yi|k|epsilon)"
{
solver PBiCG;
preconditioner DILU;
@ -41,7 +72,7 @@ solvers
}
}
PISO
PIMPLE
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;

View File

@ -47,25 +47,5 @@ graphFormat raw;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -59,6 +59,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5;
rhoMax rhoMax [ 1 -3 0 0 0 ] 1.5;
residualControl
{
p 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -47,25 +47,5 @@ graphFormat raw;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p 1e-3;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -53,6 +53,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5;
rhoMax rhoMAx [ 1 -3 0 0 0 ] 2.0;
residualControl
{
p 1e-3;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -45,25 +45,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -39,6 +39,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -45,25 +45,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -39,6 +39,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-4;
T 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -44,25 +44,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -46,6 +46,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -45,25 +45,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -39,6 +39,16 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -45,26 +45,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
G 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -46,6 +46,17 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
G 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -45,26 +45,5 @@ timePrecision 6;
runTimeModifiable true;
functions
{
residualControl1
{
type residualControl;
functionObjectLibs ( "libjobControl.so" );
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
G 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
}
// ************************************************************************* //

View File

@ -43,6 +43,17 @@ SIMPLE
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
residualControl
{
p_rgh 1e-2;
U 1e-3;
h 1e-3;
G 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -26,8 +26,7 @@ solvers
hFinal
{
solver PCG;
preconditioner DIC;
$h;
tolerance 1e-8;
relTol 0;
}
@ -37,11 +36,18 @@ solvers
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0.1;
}
hUFinal
{
$hU;
tolerance 1e-6;
relTol 0;
}
}
PISO
PIMPLE
{
nOuterCorrectors 3;
nCorrectors 1;

View File

@ -47,21 +47,6 @@ runTimeModifiable true;
functions
{
convergenceChecks
{
type residualControl;
functionObjectLibs ("libjobControl.so");
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
}
}
streamLines
{
type streamLine;

View File

@ -69,6 +69,13 @@ solvers
SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -47,21 +47,6 @@ runTimeModifiable true;
functions
{
convergenceChecks
{
type residualControl;
functionObjectLibs ("libjobControl.so");
outputControl timeStep;
outputInterval 1;
maxResiduals
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
}
}
streamLines
{
type streamLine;

View File

@ -69,6 +69,13 @@ solvers
SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
}
}
relaxationFactors

View File

@ -11,9 +11,4 @@
// Add resistance on the actuation disks
actuationDisks.addSu(UEqn());
eqnResidual = solve
(
UEqn() == -fvc::grad(p)
).initialResidual();
maxResidual = max(eqnResidual, maxResidual);
solve(UEqn() == -fvc::grad(p));

View File

@ -1,5 +0,0 @@
// initialize values for convergence checks
scalar eqnResidual = 1, maxResidual = 0;
scalar convergenceCriterion = 0;
simple.readIfPresent("convergence", convergenceCriterion);

View File

@ -7,7 +7,7 @@
adjustPhi(phi, U, p);
// Non-orthogonal pressure corrector loop
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
for (int nonOrth=0; nonOrth<=simple.nNonOrthCorr(); nonOrth++)
{
fvScalarMatrix pEqn
(
@ -15,23 +15,15 @@
);
pEqn.setReference(pRefCell, pRefValue);
// retain the residual from the first iteration
if (nonOrth == 0)
{
pEqn.solve();
}
else
{
pEqn.solve();
}
pEqn.solve();
if (nonOrth == nNonOrthCorr)
if (nonOrth == simple.nNonOrthCorr())
{
phi -= pEqn.flux();
}
}
# include "continuityErrs.H"
#include "continuityErrs.H"
// Explicitly relax pressure for momentum corrector
p.relax();

View File

@ -34,6 +34,7 @@ Description
#include "singlePhaseTransportModel.H"
#include "RASModel.H"
#include "IObasicSourceList.H"
#include "simpleControl.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -45,17 +46,16 @@ int main(int argc, char *argv[])
#include "createFields.H"
#include "initContinuityErrs.H"
simpleControl simple(mesh);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Info<< "\nStarting time loop\n" << endl;
while (runTime.loop())
while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;
#include "readSIMPLEControls.H"
#include "initConvergenceCheck.H"
p.storePrevIter();
// Pressure-velocity SIMPLE corrector
@ -71,9 +71,6 @@ int main(int argc, char *argv[])
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
#include "convergenceCheck.H"
}
Info<< "End\n" << endl;

View File

@ -60,7 +60,7 @@ solvers
relTol 0;
}
"(U|k|omega)"
"(U|k|omega).*"
{
solver smoothSolver;
@ -72,7 +72,7 @@ solvers
};
}
PISO
PIMPLE
{
momentumPredictor yes;

View File

@ -50,7 +50,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -50,7 +50,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 4;

View File

@ -17,7 +17,15 @@ FoamFile
solvers
{
"(rho|U|k|omega)"
rho
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
rhoFinal
{
solver PBiCG;
preconditioner DILU;
@ -25,6 +33,7 @@ solvers
relTol 0;
}
p
{
solver GAMG;
@ -45,17 +54,37 @@ solvers
tolerance 1e-08;
relTol 0;
}
"(k|omega)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
UFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}
}
PISO
PIMPLE
{
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}
SIMPLE
{
}
// ************************************************************************* //

View File

@ -17,11 +17,18 @@ FoamFile
solvers
{
"(rho|U|k|omega)"
rho
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
rhoFinal
{
$rho;
tolerance 1e-08;
relTol 0;
}
@ -45,17 +52,36 @@ solvers
tolerance 1e-08;
relTol 0;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
UFinal
{
$U;
tolerance 1e-08;
relTol 0;
}
"(k|omega)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}
}
PISO
PIMPLE
{
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}
SIMPLE
{
}
// ************************************************************************* //

View File

@ -17,7 +17,15 @@ FoamFile
solvers
{
"(rho|U|k|omega)"
rho
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
rhoFinal
{
solver PBiCG;
preconditioner DILU;
@ -25,6 +33,7 @@ solvers
relTol 0;
}
p
{
solver GAMG;
@ -45,16 +54,36 @@ solvers
tolerance 1e-08;
relTol 0;
}
"(k|omega)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}
UFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}
}
PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 1;
}
SIMPLE
PIMPLE
{
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}

View File

@ -39,7 +39,7 @@ solvers
maxIter 100;
}
rho
"(rho|rhoFinal)"
{
solver diagonal;
}
@ -100,7 +100,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
transSonic no;

View File

@ -39,7 +39,7 @@ solvers
maxIter 100;
}
rho
"(rho|rhoFinal)"
{
solver diagonal;
}
@ -100,7 +100,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
transSonic no;

View File

@ -97,7 +97,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -109,7 +109,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -89,7 +89,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -49,7 +49,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -49,7 +49,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -95,7 +95,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 2;

View File

@ -58,7 +58,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -59,7 +59,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 3;

View File

@ -82,12 +82,19 @@ solvers
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0;
relTol 0.1;
nSweeps 1;
}
UFinal
{
$U;
tolerance 1e-08;
relTol 0;
}
}
PISO
PIMPLE
{
nCorrectors 4;
nNonOrthogonalCorrectors 0;

View File

@ -87,7 +87,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor no;
nCorrectors 4;

View File

@ -22,29 +22,27 @@ solvers
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
relTol 0.1;
}
"U.*"
p_rghFinal
{
solver PBiCG;
preconditioner DILU;
$p_rgh;
tolerance 1e-07;
relTol 0;
}
k
"(U|Alpha|k|epsilon)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
relTol 0.1;
}
epsilon
"(U|Alpha|k|epsilon)Final"
{
solver PBiCG;
preconditioner DILU;
$k;
tolerance 1e-07;
relTol 0;
}
@ -54,19 +52,18 @@ solvers
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
relTol 0.1;
}
Alpha
rhoFinal
{
solver PBiCG;
preconditioner DILU;
$rho;
tolerance 1e-07;
relTol 0;
}
}
PISO
PIMPLE
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;

View File

@ -66,7 +66,7 @@ solvers
}
}
PISO
PIMPLE
{
momentumPredictor yes;
nOuterCorrectors 1;