ENH: Updated steadyReactingParcelFoam -> LTSReactingParcelFoam

applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2008-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2010-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -22,17 +22,13 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Application
-    porousExplicitSourceReactingParcelFoam
+    LTSReactingParcelFoam
 
 Description
-    Transient PISO solver for compressible, laminar or turbulent flow with
+    Local time stepping (LTS) solver for steady, compressible, laminar or
-    reacting multiphase Lagrangian parcels for porous media, including explicit
+    turbulent reacting and non-reacting flow with multiphase Lagrangian
-    sources for mass, momentum and energy
+    parcels and porous media, including explicit sources for mass, momentum
-
+    and energy
-    The solver includes:
-    - reacting multiphase parcel cloud
-    - porous media
-    - mass, momentum and energy sources
 
     Note: ddtPhiCorr not used here when porous zones are active
     - not well defined for porous calculations
@@ -75,7 +71,7 @@ int main(int argc, char *argv[])
 
     while (runTime.run())
     {
-        #include "readSIMPLEControls.H"
+        #include "readPISOControls.H"
         #include "readChemistryProperties.H"
         #include "readAdditionalSolutionControls.H"
         #include "readTimeControls.H"
@@ -84,14 +80,9 @@ int main(int argc, char *argv[])
 
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
-        p.storePrevIter();
-
-        // --- Pressure-velocity corrector
-        {
         parcels.evolve();
 
         #include "chemistry.H"
-
         #include "timeScales.H"
 
         #include "rhoEqn.H"
@@ -99,10 +90,13 @@ int main(int argc, char *argv[])
         #include "YEqn.H"
         #include "hsEqn.H"
 
+        // --- PISO loop
+        for (int corr=0; corr<nCorr; corr++)
+        {
             #include "pEqn.H"
+        }
 
         turbulence->correct();
-        }
 
         if (runTime.write())
         {

applications/solvers/lagrangian/LTSReactingParcelFoam/Make/files

@@ -0,0 +1,3 @@
+LTSReactingParcelFoam.C
+
+EXE = $(FOAM_APPBIN)/LTSReactingParcelFoam

applications/solvers/lagrangian/LTSReactingParcelFoam/createFields.H

@@ -73,12 +73,12 @@
 
     dimensionedScalar rhoMax
     (
-        mesh.solutionDict().subDict("SIMPLE").lookup("rhoMax")
+        mesh.solutionDict().subDict("PISO").lookup("rhoMax")
     );
 
     dimensionedScalar rhoMin
     (
-        mesh.solutionDict().subDict("SIMPLE").lookup("rhoMin")
+        mesh.solutionDict().subDict("PISO").lookup("rhoMin")
     );
 
     Info<< "Creating turbulence model\n" << endl;

applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H

@@ -48,9 +48,6 @@
         }
     }
 
-    // Explicitly relax pressure for momentum corrector
-    p.relax();
-
     Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
 
     // Second part of thermodynamic density update

applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2010-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2010-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License

applications/solvers/lagrangian/LTSReactingParcelFoam/rhoEqn.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2008-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2008-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License

applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2010-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2010-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License

applications/solvers/lagrangian/steadyReactingParcelFoam/Make/files

@@ -1,3 +0,0 @@
-steadyReactingParcelFoam.C
-
-EXE = $(FOAM_APPBIN)/steadyReactingParcelFoam

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun

@@ -8,4 +8,4 @@ cd ${0%/*} || exit 1    # run from this directory
 runApplication blockMesh
 
 # run the solver
-runApplication steadyReactingParcelFoam
+runApplication LTSReactingParcelFoam

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict

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

tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution

@@ -53,8 +53,9 @@ solvers
     }
 }
 
-SIMPLE
+PISO
 {
+    nCorrectors         2;
     nNonOrthogonalCorrectors 0;
     momentumPredictor yes;
     rhoMin rhoMin [1 -3 0 0 0] 0.1;

tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allrun

@@ -8,6 +8,6 @@ cd ${0%/*} || exit 1    # run from this directory
 runApplication blockMesh
 
 # run the solver
-runApplication steadyReactingParcelFoam
+runApplication LTSReactingParcelFoam
 
 # ----------------------------------------------------------------- end-of-file

tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution

@@ -61,8 +61,9 @@ solvers
     };
 }
 
-SIMPLE
+PISO
 {
+    nCorrectors          2;
     nNonOrthogonalCorrectors 0; // 10;
     momentumPredictor yes;
     rhoMin rhoMin [ 1 -3 0 0 0 ] 0; // 0.5;