ENH: lagrangian tutorial updates

This commit is contained in:
andy
2014-01-24 12:38:08 +00:00
parent b6ff378e00
commit b95c4f9ee7
11 changed files with 120 additions and 62 deletions

View File

@ -25,7 +25,7 @@ solvers
p p
{ {
solver GAMG; solver GAMG;
tolerance 0; tolerance 1e-6;
relTol 0.1; relTol 0.1;
smoother GaussSeidel; smoother GaussSeidel;
nPreSweeps 0; nPreSweeps 0;

View File

@ -20,7 +20,7 @@ solvers
p p
{ {
solver GAMG; solver GAMG;
tolerance 0; tolerance 1e-6;
relTol 0.01; relTol 0.01;
smoother GaussSeidel; smoother GaussSeidel;
nPreSweeps 0; nPreSweeps 0;

View File

@ -32,15 +32,16 @@ solvers
relTol 0; relTol 0;
} }
"(U|h)" U
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.1;
maxIter 50;
} }
"(U|h)Final" UFinal
{ {
$U; $U;
relTol 0; relTol 0;
@ -73,9 +74,19 @@ solvers
relTol 0; relTol 0;
} }
"(Yi|O2|N2|H2O)" "(h|Yi|O2|N2|H2O)"
{ {
$hFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
hFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
} }
} }

View File

@ -83,12 +83,12 @@ subModels
((0.6 0.35 1.45) (0 0 -1)) ((0.6 0.35 1.45) (0 0 -1))
); );
massTotal 10; massTotal 5;
parcelsPerInjector 20000; parcelsPerInjector 20000;
parcelsPerSecond 500; parcelsPerSecond 500;
parcelBasisType mass; parcelBasisType mass;
flowRateProfile constant 0.1; flowRateProfile constant 1;
Umag constant 3.0; Umag constant 1;
thetaInner constant 0; thetaInner constant 0;
thetaOuter constant 45; thetaOuter constant 45;

View File

@ -34,9 +34,9 @@ solvers
"(U|h|k|epsilon)" "(U|h|k|epsilon)"
{ {
solver smoothSolver; solver PBiCG;
smoother symGaussSeidel; preconditioner DILU;
tolerance 0; tolerance 1e-5;
relTol 0.1; relTol 0.1;
} }

View File

@ -32,31 +32,27 @@ solvers
relTol 0; relTol 0;
} }
"(U|h)" U
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.1;
maxIter 50;
} }
"(U|h)Final" UFinal
{ {
$U; $U;
relTol 0; relTol 0;
tolerance 1e-05; tolerance 1e-05;
} }
"(k|epsilon)"
{
$UFinal;
}
p_rgh p_rgh
{ {
solver GAMG; solver GAMG;
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.01;
smoother GaussSeidel; smoother GaussSeidel;
nPreSweeps 0; nPreSweeps 0;
nPostSweeps 2; nPostSweeps 2;
@ -64,18 +60,29 @@ solvers
nCellsInCoarsestLevel 10; nCellsInCoarsestLevel 10;
agglomerator faceAreaPair; agglomerator faceAreaPair;
mergeLevels 1; mergeLevels 1;
maxIter 50;
} }
p_rghFinal p_rghFinal
{ {
$p_rgh; $p_rgh;
tolerance 1e-06; tolerance 1e-05;
relTol 0; relTol 0;
} }
"(Yi|O2|N2|H2O)" "(h|Yi|O2|N2|H2O)"
{ {
$hFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
hFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
} }
} }
@ -90,11 +97,11 @@ relaxationFactors
{ {
fields fields
{ {
".*Final" 1; ".*" 1;
} }
equations equations
{ {
".*Final" 1; ".*" 1;
} }
} }

View File

@ -19,15 +19,15 @@ solvers
{ {
"(Uf|hf|deltaf\*rhof)" "(Uf|hf|deltaf\*rhof)"
{ {
solver smoothSolver; solver PBiCG;
smoother symGaussSeidel; preconditioner DILU;
tolerance 1e-10; tolerance 1e-10;
relTol 0; relTol 0;
} }
deltaf deltaf
{ {
solver smoothSolver; solver PBiCG;
smoother symGaussSeidel; preconditioner DILU;
tolerance 1e-10; tolerance 1e-10;
relTol 0; relTol 0;
} }
@ -46,9 +46,11 @@ relaxationFactors
{ {
fields fields
{ {
".*" 1;
} }
equations equations
{ {
".*" 1;
} }
} }

View File

@ -32,7 +32,7 @@ solvers
relTol 0; relTol 0;
} }
"(U|h|k|epsilon)" "(U|k|epsilon)"
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
@ -40,7 +40,7 @@ solvers
relTol 0.1; relTol 0.1;
} }
"(U|h|k|epsilon)Final" "(U|k|epsilon)Final"
{ {
$U; $U;
tolerance 1e-06; tolerance 1e-06;
@ -70,7 +70,21 @@ solvers
"(Yi|O2|N2|H2O)" "(Yi|O2|N2|H2O)"
{ {
$hFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
h
{
$Yi;
relTol 0.1;
}
hFinal
{
$Yi;
} }
} }

View File

@ -32,7 +32,7 @@ solvers
relTol 0; relTol 0;
} }
"(U|h|k|epsilon)" "(U|k|epsilon)"
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
@ -40,7 +40,7 @@ solvers
relTol 0.1; relTol 0.1;
} }
"(U|h|k|epsilon)Final" "(U|k|epsilon)Final"
{ {
$U; $U;
tolerance 1e-05; tolerance 1e-05;
@ -71,7 +71,19 @@ solvers
"(Yi|O2|N2|H2O)" "(Yi|O2|N2|H2O)"
{ {
$hFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
h
{
$Yi;
relTol 0.1;
}
hFinal
{
$Yi;
} }
} }

View File

@ -32,7 +32,7 @@ solvers
relTol 0; relTol 0;
} }
"(U|h|k|omega)" "(U|k|omega)"
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
@ -40,7 +40,7 @@ solvers
relTol 0.1; relTol 0.1;
} }
"(U|h|k|omega)Final" "(U|k|omega)Final"
{ {
$U; $U;
tolerance 1e-06; tolerance 1e-06;
@ -59,7 +59,6 @@ solvers
nCellsInCoarsestLevel 10; nCellsInCoarsestLevel 10;
agglomerator faceAreaPair; agglomerator faceAreaPair;
mergeLevels 1; mergeLevels 1;
maxIter 50;
} }
pFinal pFinal
@ -71,23 +70,31 @@ solvers
"(Yi|O2|N2|H2O)" "(Yi|O2|N2|H2O)"
{ {
$hFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
h
{
$Yi;
relTol 0.1;
}
hFinal
{
$Yi;
} }
} }
PIMPLE PIMPLE
{ {
nCorrectors 3; transonic no;
nOuterCorrectors 1;
nCorrectors 2;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
momentumPredictor yes; momentumPredictor yes;
} }
potentialFlow
{
// used for potentialFoam initialisation
nNonOrthogonalCorrectors 5;
}
relaxationFactors relaxationFactors
{ {
fields fields
@ -95,8 +102,9 @@ relaxationFactors
} }
equations equations
{ {
".*" 1; ".*Final" 1;
} }
} }
// ************************************************************************* // // ************************************************************************* //

View File

@ -32,11 +32,11 @@ solvers
relTol 0; relTol 0;
} }
"(U|h)" "(U|k|epsilon)"
{ {
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
tolerance 1e-05; tolerance 1e-06;
relTol 0.1; relTol 0.1;
} }
@ -61,26 +61,30 @@ solvers
relTol 0; relTol 0;
} }
"(U|h)Final" "(U|k|epsilon)Final"
{ {
$U; $U;
tolerance 1e-05; tolerance 1e-06;
relTol 0; relTol 0;
} }
"(k|epsilon)"
{
$U;
}
"(k|epsilon)Final"
{
$UFinal;
}
"(Yi|O2|N2|H2O)" "(Yi|O2|N2|H2O)"
{ {
$UFinal; solver PBiCG;
preconditioner DILU;
tolerance 1e-6;
relTol 0;
}
h
{
$Yi;
relTol 0.1;
}
hFinal
{
$Yi;
} }
} }