ENH: Updated LTSReactingParcelFoam tutorials

This commit is contained in:
andy
2011-03-25 14:35:50 +00:00
parent 3bb159794f
commit 28cd7fc92f
7 changed files with 68 additions and 49 deletions

View File

@ -51,8 +51,6 @@ maxCo 0.5;
maxDeltaT 1; maxDeltaT 1;
alphaTauSpecie 1e10;
alphaTauTemp 0.005; alphaTauTemp 0.005;
// ************************************************************************* // // ************************************************************************* //

View File

@ -17,7 +17,7 @@ FoamFile
ddtSchemes ddtSchemes
{ {
default localEuler invTauFlow; default localEuler invTau;
} }
gradSchemes gradSchemes

View File

@ -24,6 +24,13 @@ solvers
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.1;
} }
rhoFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-4;
relTol 0;
}
p p
{ {
solver GAMG; solver GAMG;
@ -44,6 +51,14 @@ solvers
maxIter 50; maxIter 50;
}; };
pFinal
{
$p;
tolerance 1e-4;
relTol 0;
};
"(U|Yi|hs|k|epsilon)" "(U|Yi|hs|k|epsilon)"
{ {
solver smoothSolver; solver smoothSolver;
@ -51,9 +66,17 @@ solvers
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.1;
} }
"(U|Yi|hs|k|epsilon)Final"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-4;
relTol 0;
}
} }
PISO PIMPLE
{ {
nCorrectors 2; nCorrectors 2;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
@ -64,8 +87,8 @@ PISO
additional additional
{ {
eWork true; pressureWork true;
hWork true; pressureWorkTimeDerivative true;
solveSpecies true; solveSpecies true;
} }

View File

@ -122,7 +122,7 @@ subModels
massFlowRate 0.8e-03; massFlowRate 0.8e-03;
parcelBasisType mass; parcelBasisType mass;
patchName inletCentral; patchName inletCentral;
parcelsPerSecond 400; parcelsPerSecond 100;
duration 1; // NOTE: set to 1 for steady state duration 1; // NOTE: set to 1 for steady state
U0 (0 40 0); U0 (0 40 0);
flowRateProfile constant 1; flowRateProfile constant 1;

View File

@ -47,15 +47,11 @@ runTimeModifiable yes;
maxCo 5; maxCo 5;
alphaTauSmooth 1e10; alphaTauSmooth 1;
alphaTauTemp 1e10; alphaTauTemp 1;
alphaTauRho 1e10; maxDeltaT 1;
alphaTauSpecie 1e10;
maxDeltaT 1e10;
functions functions
{ {

View File

@ -17,7 +17,7 @@ FoamFile
ddtSchemes ddtSchemes
{ {
default localEuler invTauFlow; default localEuler invTau;
} }
gradSchemes gradSchemes

View File

@ -1,14 +1,14 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev | | \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.com | | \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
FoamFile FoamFile
{ {
version 2.0; version 2.0;
format ascii; format binary;
class dictionary; class dictionary;
location "system"; location "system";
object fvSolution; object fvSolution;
@ -17,11 +17,21 @@ FoamFile
solvers solvers
{ {
"(rho|.*Final)"
{
solver PCG;
preconditioner DIC;
tolerance 0;
relTol 0.1;
}
// "(p|.*Final)"
p p
{ {
solver GAMG; solver GAMG;
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.05;
smoother DICGaussSeidel; smoother DICGaussSeidel;
nPreSweeps 0; nPreSweeps 0;
@ -33,47 +43,39 @@ solvers
agglomerator faceAreaPair; agglomerator faceAreaPair;
mergeLevels 1; mergeLevels 1;
maxIter 20; maxIter 50;
}; };
"(rho|G)"
pFinal
{ {
solver PCG; $p;
preconditioner DIC; tolerance 1e-4;
tolerance 0; }
relTol 0.1;
maxIter 20; "(U|Yi|hs|k|omega|.*Final)"
};
"(Yi|hs)"
{ {
solver smoothSolver; solver smoothSolver;
smoother GaussSeidel; smoother GaussSeidel;
tolerance 0; tolerance 0;
relTol 0.1; relTol 0.1;
maxIter 20; maxIter 20;
}; }
"(U|k|omega)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 0;
relTol 0.1;
maxIter 20;
};
} }
PISO PIMPLE
{ {
nOuterCorrectors 1;
nCorrectors 2; nCorrectors 2;
nNonOrthogonalCorrectors 0; // 10; nNonOrthogonalCorrectors 0;
momentumPredictor yes; momentumPredictor yes;
rhoMin rhoMin [ 1 -3 0 0 0 ] 0; // 0.5; rhoMin rhoMin [1 -3 0 0 0] 0.1;
rhoMax rhoMax [ 1 -3 0 0 0 ] 100; // 1.5; rhoMax rhoMax [1 -3 0 0 0] 1.5;
} }
additional additional
{ {
eWork true; pressureWork true;
hWork true; pressureWorkTimeDerivative false; // true;
solveSpecies true; solveSpecies true;
} }