ENH: Tutorial update

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T

@@ -36,7 +36,6 @@ boundaryField
         type            inletOutlet;
         inletValue      uniform 293;
         value           uniform 293;
-
     }
     frontAndBack
     {

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U

@@ -33,7 +33,8 @@ boundaryField
     }
     outlet
     {
-        type            zeroGradient;
+        type            pressureInletOutletVelocity;
+        value           $internalField;
     }
     frontAndBack
     {

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon

@@ -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;
-    }
-}
-
-
-// ************************************************************************* //

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k

@@ -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;
-    }
-}
-
-
-// ************************************************************************* //

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p

@@ -29,14 +29,16 @@ boundaryField
     {
         type            zeroGradient;
     }
-    inlet2
-    {
-        type            zeroGradient;
-    }
     outlet
     {
-        type            fixedValue;
+        type            totalPressure;
-        value           uniform 1e5;
+        p0              $internalField;
+        U               U;
+        phi             phi;
+        rho             none;
+        psi             none;
+        gamma           1;
+        value           $internalField;
     }
     frontAndBack
     {

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;
+}
 
 // ************************************************************************* //

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties

@@ -15,15 +15,12 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-combustionModel  PaSR<rhoChemistryCombustion>;
+combustionModel  laminar<rhoChemistryCombustion>;
 
 active  true;
 
-PaSRCoeffs
+laminarCoeffs
 {
-    Cmix                0.1;
-    turbulentReaction   off;
-    useReactionRate     true;
 }
 
 

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary

@@ -9,37 +9,39 @@ FoamFile
 {
     version     2.0;
     format      ascii;
-    class       volScalarField;
+    class       polyBoundaryMesh;
-    location    "0";
+    location    "constant/polyMesh";
-    object      alphat;
+    object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-dimensions      [1 -1 -1 0 0 0 0];
+4
-
+(
-internalField   uniform 0;
-
-boundaryField
-{
     fuel
     {
-        type            fixedValue;
+        type            patch;
-        value           uniform 0;
+        nFaces          40;
+        startFace       7860;
     }
     air
     {
-        type            fixedValue;
+        type            patch;
-        value           uniform 0;
+        nFaces          40;
+        startFace       7900;
     }
     outlet
     {
-        type            zeroGradient;
+        type            patch;
+        nFaces          200;
+        startFace       7940;
     }
     frontAndBack
     {
         type            empty;
+        inGroups        1(empty);
+        nFaces          8000;
+        startFace       8140;
     }
-}
+)
-
 
 // ************************************************************************* //

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties

@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  RASModel;
+simulationType  laminar;
 
 
 // ************************************************************************* //

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict

@@ -23,7 +23,7 @@ startTime       0;
 
 stopAt          endTime;
 
-endTime         10000;
+endTime         2000;
 
 deltaT          1;
 

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution

@@ -17,79 +17,69 @@ FoamFile
 
 solvers
 {
-    rho
+    "rho.*"
     {
-        solver          PCG;
+        solver          diagonal;
-        preconditioner  DIC;
-        tolerance       0;
-        relTol          0.1;
-    }
-
-    rhoFinal
-    {
-        solver          PCG;
-        preconditioner  DIC;
-        tolerance       1e-4;
-        relTol          0;
     }
 
     p
     {
-        solver           GAMG;
+        solver           PCG;
-
+        preconditioner   DIC;
-        tolerance        0;
+        tolerance        1e-6;
-        relTol           0.05;
+        relTol           0.1;
-
+    }
-        smoother         DICGaussSeidel;
-        nPreSweeps       0;
-        nPostSweeps      2;
-
-        cacheAgglomeration true;
-
-        nCellsInCoarsestLevel 10;
-        agglomerator     faceAreaPair;
-        mergeLevels      1;
-
-        maxIter          50;
-    };
 
     pFinal
     {
         $p;
-
-        tolerance        1e-4;
         relTol           0;
-    };
+    }
 
-    "(U|Yi|h|k|epsilon)"
+    "(U|h|k|epsilon)"
     {
-        solver          smoothSolver;
+        solver          PBiCG;
-        smoother        GaussSeidel;
+        preconditioner  DILU;
-        tolerance       0;
+        tolerance       1e-6;
         relTol          0.1;
     }
 
-    "(U|Yi|h|k|epsilon)Final"
+    "(U|h|k|epsilon)Final"
     {
-        solver          smoothSolver;
+        $U;
-        smoother        GaussSeidel;
+        relTol          0.1;
-        tolerance       1e-4;
+    }
-        relTol          0;
+
+    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;
+    maxDeltaT       1e-2;
-    rDeltaTSmoothingCoeff 0.1;
+    maxCo           1;
-    maxDeltaT       2e-4;
+    alphaTemp       0.05;
-    alphaTemp       0.005;
+    rDeltaTSmoothingCoeff 1;
+    rDeltaTDampingCoeff 1;
+}
+
+relaxationFactors
+{
+    fields
+    {
+    }
+    equations
+    {
+        ".*" 1;
+    }
 }