From 84e3fb4db2893895580e17776ae81c7fa35a2823 Mon Sep 17 00:00:00 2001 From: andy Date: Tue, 19 Nov 2013 16:38:44 +0000 Subject: [PATCH] ENH: Tutorial update --- .../counterFlowFlame2D/0/T | 1 - .../counterFlowFlame2D/0/U | 3 +- .../counterFlowFlame2D/0/epsilon | 45 ---------- .../counterFlowFlame2D/0/k | 45 ---------- .../counterFlowFlame2D/0/p | 14 +-- .../constant/chemistryProperties | 19 +++- .../constant/combustionProperties | 7 +- .../{0/alphat => constant/polyMesh/boundary} | 34 +++---- .../constant/turbulenceProperties | 2 +- .../counterFlowFlame2D/system/controlDict | 2 +- .../counterFlowFlame2D/system/fvSolution | 90 +++++++++---------- 11 files changed, 87 insertions(+), 175 deletions(-) delete mode 100644 tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon delete mode 100644 tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k rename tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/{0/alphat => constant/polyMesh/boundary} (66%) diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T index d5d567628e..fcf963a775 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T @@ -36,7 +36,6 @@ boundaryField type inletOutlet; inletValue uniform 293; value uniform 293; - } frontAndBack { diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U index 27773aea2b..c235ba3d65 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U @@ -33,7 +33,8 @@ boundaryField } outlet { - type zeroGradient; + type pressureInletOutletVelocity; + value $internalField; } frontAndBack { diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon deleted file mode 100644 index ad65f7ea0f..0000000000 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon +++ /dev/null @@ -1,45 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ -| ========= | | -| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | -| \\ / A nd | Web: www.OpenFOAM.org | -| \\/ M anipulation | | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - location "0"; - object epsilon; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 2 -3 0 0 0 0]; - -internalField uniform 50; - -boundaryField -{ - fuel - { - type fixedValue; - value uniform 1e-05; - } - air - { - type fixedValue; - value uniform 1e-05; - } - outlet - { - type zeroGradient; - } - frontAndBack - { - type empty; - } -} - - -// ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k deleted file mode 100644 index 3dcedaa2da..0000000000 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k +++ /dev/null @@ -1,45 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ -| ========= | | -| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | -| \\ / A nd | Web: www.OpenFOAM.org | -| \\/ M anipulation | | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - location "0"; - object k; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 2 -2 0 0 0 0]; - -internalField uniform 3.75e-05; - -boundaryField -{ - fuel - { - type fixedValue; - value uniform 3.75e-05; - } - air - { - type fixedValue; - value uniform 3.75e-05; - } - outlet - { - type zeroGradient; - } - frontAndBack - { - type empty; - } -} - - -// ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p index 79f071ce31..deca9ba60d 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p @@ -29,14 +29,16 @@ boundaryField { type zeroGradient; } - inlet2 - { - type zeroGradient; - } outlet { - type fixedValue; - value uniform 1e5; + type totalPressure; + p0 $internalField; + U U; + phi phi; + rho none; + psi none; + gamma 1; + value $internalField; } frontAndBack { diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties index 254f5c4633..ebdde7591f 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties @@ -17,14 +17,25 @@ FoamFile chemistryType { - chemistrySolver noChemistrySolver; + chemistrySolver EulerImplicit; chemistryThermo rho; } -chemistry on; // off; +chemistry on; -chemCalcFreq 1; +initialChemicalTimeStep 1e-07; -initialChemicalTimeStep 1e-8; // NOT USED +EulerImplicitCoeffs +{ + cTauChem 1; + equilibriumRateLimiter off; +} + +odeCoeffs +{ + solver Rosenbrock43; + absTol 1e-12; + relTol 0.01; +} // ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties index 06715e90f1..bad68a398b 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties @@ -15,15 +15,12 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -combustionModel PaSR; +combustionModel laminar; active true; -PaSRCoeffs +laminarCoeffs { - Cmix 0.1; - turbulentReaction off; - useReactionRate true; } diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/alphat b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary similarity index 66% rename from tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/alphat rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary index c74018bbdb..4ba3243474 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/alphat +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary @@ -9,37 +9,39 @@ FoamFile { version 2.0; format ascii; - class volScalarField; - location "0"; - object alphat; + class polyBoundaryMesh; + location "constant/polyMesh"; + object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [1 -1 -1 0 0 0 0]; - -internalField uniform 0; - -boundaryField -{ +4 +( fuel { - type fixedValue; - value uniform 0; + type patch; + nFaces 40; + startFace 7860; } air { - type fixedValue; - value uniform 0; + type patch; + nFaces 40; + startFace 7900; } outlet { - type zeroGradient; + type patch; + nFaces 200; + startFace 7940; } frontAndBack { type empty; + inGroups 1(empty); + nFaces 8000; + startFace 8140; } -} - +) // ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties index 3721a46a2e..c2c3b28a1b 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties @@ -15,7 +15,7 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -simulationType RASModel; +simulationType laminar; // ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict index bb20dffbf2..e27ca26539 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict @@ -23,7 +23,7 @@ startTime 0; stopAt endTime; -endTime 10000; +endTime 2000; deltaT 1; diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution index 552f948a84..cf6feec247 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution @@ -17,79 +17,69 @@ FoamFile solvers { - rho + "rho.*" { - solver PCG; - preconditioner DIC; - tolerance 0; - relTol 0.1; - } - - rhoFinal - { - solver PCG; - preconditioner DIC; - tolerance 1e-4; - relTol 0; + solver diagonal; } p { - solver GAMG; - - tolerance 0; - relTol 0.05; - - smoother DICGaussSeidel; - nPreSweeps 0; - nPostSweeps 2; - - cacheAgglomeration true; - - nCellsInCoarsestLevel 10; - agglomerator faceAreaPair; - mergeLevels 1; - - maxIter 50; - }; + solver PCG; + preconditioner DIC; + tolerance 1e-6; + relTol 0.1; + } pFinal { $p; - - tolerance 1e-4; relTol 0; - }; + } - "(U|Yi|h|k|epsilon)" + "(U|h|k|epsilon)" { - solver smoothSolver; - smoother GaussSeidel; - tolerance 0; + solver PBiCG; + preconditioner DILU; + tolerance 1e-6; relTol 0.1; } - "(U|Yi|h|k|epsilon)Final" + "(U|h|k|epsilon)Final" { - solver smoothSolver; - smoother GaussSeidel; - tolerance 1e-4; - relTol 0; + $U; + relTol 0.1; + } + + Yi + { + $U; + relTol 0.1; } } PIMPLE { - nCorrectors 2; + momentumPredictor no; + nOuterCorrectors 1; + nCorrectors 1; nNonOrthogonalCorrectors 0; - momentumPredictor yes; - rhoMin 0.1; - rhoMax 1.5; - maxCo 0.1; - rDeltaTSmoothingCoeff 0.1; - maxDeltaT 2e-4; - alphaTemp 0.005; + maxDeltaT 1e-2; + maxCo 1; + alphaTemp 0.05; + rDeltaTSmoothingCoeff 1; + rDeltaTDampingCoeff 1; +} + +relaxationFactors +{ + fields + { + } + equations + { + ".*" 1; + } }