applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/porousExplicitSourceReactingParcelFoam.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Application
  25     porousExplicitSourceReactingParcelFoam
  26
  27 Description
  28     Transient PISO solver for compressible, laminar or turbulent flow with
  29     reacting multiphase Lagrangian parcels for porous media, including explicit
  30     sources for mass, momentum and energy
  31
  32     The solver includes:
  33     - reacting multiphase parcel cloud
  34     - porous media
  35     - mass, momentum and energy sources
  36
  37     Note: ddtPhiCorr not used here when porous zones are active
  38     - not well defined for porous calculations
  39
  40 \*---------------------------------------------------------------------------*/
  41
  42 #include "fvCFD.H"
  43 #include "hReactionThermo.H"
  44 #include "turbulenceModel.H"
  45 #include "basicReactingMultiphaseCloud.H"
  46 #include "rhoChemistryModel.H"
  47 #include "chemistrySolver.H"
  48 #include "radiationModel.H"
  49 #include "porousZones.H"
  50 #include "timeActivatedExplicitSource.H"
  51 #include "SLGThermo.H"
  52 #include "pimpleControl.H"
  53
  54 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  55
  56 int main(int argc, char *argv[])
  57 {
  58     #include "setRootCase.H"
  59
  60     #include "createTime.H"
  61     #include "createMesh.H"
  62     #include "readChemistryProperties.H"
  63     #include "readGravitationalAcceleration.H"
  64     #include "createFields.H"
  65     #include "createRadiationModel.H"
  66     #include "createClouds.H"
  67     #include "createExplicitSources.H"
  68     #include "createPorousZones.H"
  69     #include "initContinuityErrs.H"
  70     #include "readTimeControls.H"
  71     #include "compressibleCourantNo.H"
  72     #include "setInitialDeltaT.H"
  73
  74     pimpleControl pimple(mesh);
  75
  76     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  77
  78     Info<< "\nStarting time loop\n" << endl;
  79
  80     while (runTime.run())
  81     {
  82         #include "readTimeControls.H"
  83         #include "readAdditionalSolutionControls.H"
  84         #include "compressibleCourantNo.H"
  85         #include "setDeltaT.H"
  86
  87         runTime++;
  88
  89         Info<< "Time = " << runTime.timeName() << nl << endl;
  90
  91         parcels.evolve();
  92
  93         #include "chemistry.H"
  94         #include "rhoEqn.H"
  95
  96         // --- Pressure-velocity PIMPLE corrector loop
  97         for (pimple.start(); pimple.loop(); pimple++)
  98         {
  99             #include "UEqn.H"
 100             #include "YEqn.H"
 101             #include "hsEqn.H"
 102
 103             // --- PISO loop
 104             for (int corr=0; corr<pimple.nCorr(); corr++)
 105             {
 106                 #include "pEqn.H"
 107             }
 108
 109             if (pimple.turbCorr())
 110             {
 111                 turbulence->correct();
 112             }
 113         }
 114
 115         rho = thermo.rho();
 116
 117         if (runTime.write())
 118         {
 119             chemistry.dQ()().write();
 120         }
 121
 122         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 123             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 124             << nl << endl;
 125     }
 126
 127     Info<< "End\n" << endl;
 128
 129     return(0);
 130 }
 131
 132
 133 // ************************************************************************* //