Rationalized heat release rate functions

Combined 'dQ()' and 'Sh()' into 'Qdot()' which returns the heat-release rate in the normal units [kg/m/s3] and used as the heat release rate source term in the energy equations, to set the field 'Qdot' in several combustion solvers and for the evaluation of the local time-step when running LTS.
2025-11-28 03:28:01 +00:00 · 2016-12-15 17:10:21 +00:00
parent 8b5516abd9
commit b99817d924
54 changed files with 142 additions and 366 deletions
--- a/applications/solvers/combustion/chemFoam/output.H
+++ b/applications/solvers/combustion/chemFoam/output.H
@ -1,6 +1,6 @@
    runTime.write();

-    Info<< "Sh = " << Sh
+    Info<< "Qdot = " << Qdot
        << ", T = " << thermo.T()[0]
        << ", p = " << thermo.p()[0]
        << ", " << Y[0].name() << " = " << Y[0][0]
@ -8,4 +8,3 @@

    post<< runTime.value() << token::TAB << thermo.T()[0] << token::TAB
        << thermo.p()[0] << endl;
-
--- a/applications/solvers/combustion/chemFoam/solveChemistry.H
+++ b/applications/solvers/combustion/chemFoam/solveChemistry.H
@ -1,3 +1,3 @@
-    dtChem = chemistry.solve(runTime.deltaT().value());
-    scalar Sh = chemistry.Sh()()[0]/rho[0];
-    integratedHeat += Sh*runTime.deltaT().value();
+dtChem = chemistry.solve(runTime.deltaT().value());
+scalar Qdot = chemistry.Qdot()()[0]/rho[0];
+integratedHeat += Qdot*runTime.deltaT().value();
--- a/applications/solvers/combustion/fireFoam/YEEqn.H
+++ b/applications/solvers/combustion/fireFoam/YEEqn.H
@ -11,7 +11,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection
 {
    radiation->correct();
    combustion->correct();
-    dQ = combustion->dQ();
+    Qdot = combustion->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
@ -66,7 +66,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
-        combustion->Sh()
+        Qdot
      + radiation->Sh(thermo)
      + parcels.Sh(he)
      + surfaceFilm.Sh()
--- a/applications/solvers/combustion/fireFoam/createFields.H
+++ b/applications/solvers/combustion/fireFoam/createFields.H
@ -131,18 +131,18 @@ Switch solvePyrolysisRegion
    additionalControlsDict.lookupOrDefault<bool>("solvePyrolysisRegion", true)
 );

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );


--- a/applications/solvers/combustion/reactingFoam/EEqn.H
+++ b/applications/solvers/combustion/reactingFoam/EEqn.H
@ -17,7 +17,7 @@
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
-        reaction->Sh()
+        Qdot
      + fvOptions(rho, he)
    );

--- a/applications/solvers/combustion/reactingFoam/YEqn.H
+++ b/applications/solvers/combustion/reactingFoam/YEqn.H
@ -11,7 +11,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection

 {
    reaction->correct();
-    dQ = reaction->dQ();
+    Qdot = reaction->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
--- a/applications/solvers/combustion/reactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/createFields.H
@ -117,18 +117,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H
@ -117,18 +117,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
@ -96,18 +96,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/combustion/reactingFoam/setRDeltaT.H
+++ b/applications/solvers/combustion/reactingFoam/setRDeltaT.H
@ -77,7 +77,7 @@ License
    {
        volScalarField::Internal rDeltaTT
        (
-            mag(reaction->Sh())/(alphaTemp*rho*thermo.Cp()*T)
+            mag(Qdot)/(alphaTemp*rho*thermo.Cp()*T)
        );

        Info<< "    Temperature = "
--- a/applications/solvers/lagrangian/coalChemistryFoam/EEqn.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/EEqn.H
@ -18,7 +18,7 @@
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
        rho*(U&g)
-      + combustion->Sh()
+      + Qdot
      + coalParcels.Sh(he)
      + limestoneParcels.Sh(he)
      + radiation->Sh(thermo)
--- a/applications/solvers/lagrangian/coalChemistryFoam/YEqn.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/YEqn.H
@ -12,7 +12,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection

 {
    combustion->correct();
-    dQ = combustion->dQ();
+    Qdot = combustion->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
--- a/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
@ -131,18 +131,18 @@ volScalarField dpdt
 Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/lagrangian/coalChemistryFoam/setRDeltaT.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/setRDeltaT.H
@ -80,7 +80,7 @@ License
            (
               (coalParcels.hsTrans() + limestoneParcels.hsTrans())
              /(mesh.V()*runTime.deltaT())
-             + combustion->Sh()()
+             + Qdot
            )
           /(
               alphaTemp
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/EEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/EEqn.H
@ -21,7 +21,7 @@
      + parcels.Sh(he)
      + surfaceFilm.Sh()
      + radiation->Sh(thermo)
-      + combustion->Sh()
+      + Qdot
      + fvOptions(rho, he)
    );

--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/YEqn.H
@ -12,7 +12,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection

 {
    combustion->correct();
-    dQ = combustion->dQ();
+    Qdot = combustion->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
--- a/applications/solvers/lagrangian/reactingParcelFilmFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFilmFoam/createFields.H
@ -134,18 +134,18 @@ Switch solvePrimaryRegion
    additionalControlsDict.lookup("solvePrimaryRegion")
 );

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/lagrangian/reactingParcelFoam/EEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/EEqn.H
@ -20,7 +20,7 @@
        rho*(U&g)
      + parcels.Sh(he)
      + radiation->Sh(thermo)
-      + combustion->Sh()
+      + Qdot
      + fvOptions(rho, he)
    );

--- a/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H
@ -11,7 +11,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection

 {
    combustion->correct();
-    dQ = combustion->dQ();
+    Qdot = combustion->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
--- a/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
@ -121,18 +121,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/lagrangian/reactingParcelFoam/setRDeltaT.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/setRDeltaT.H
@ -79,7 +79,7 @@ License
            mag
            (
               parcels.hsTrans()/(mesh.V()*runTime.deltaT())
-             + combustion->Sh()()
+             + Qdot
            )
           /(
               alphaTemp
--- a/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/EEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/EEqn.H
@ -14,7 +14,7 @@
        rho*(U&g)
      + parcels.Sh(he)
      + radiation->Sh(thermo)
-      + combustion->Sh()
+      + Qdot
      + fvOptions(rho, he)
    );

--- a/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/YEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/YEqn.H
@ -11,7 +11,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection

 {
    combustion->correct();
-    dQ = combustion->dQ();
+    Qdot = combustion->Qdot();
    volScalarField Yt(0.0*Y[0]);

    forAll(Y, i)
--- a/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/createFields.H
@ -103,18 +103,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
-        IOobject::NO_READ,
+        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/lagrangian/sprayFoam/createFields.H
+++ b/applications/solvers/lagrangian/sprayFoam/createFields.H
@ -117,18 +117,18 @@ forAll(Y, i)
 }
 fields.add(thermo.he());

-volScalarField dQ
+volScalarField Qdot
 (
    IOobject
    (
-        "dQ",
+        "Qdot",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
-    dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
+    dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
 );

 #include "createMRF.H"
--- a/applications/solvers/lagrangian/sprayFoam/sprayDyMFoam/sprayDyMFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/sprayDyMFoam/sprayDyMFoam.C
@ -135,7 +135,7 @@ int main(int argc, char *argv[])

        if (runTime.write())
        {
-            combustion->dQ()().write();
+            combustion->Qdot()().write();
        }

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
--- a/applications/solvers/lagrangian/sprayFoam/sprayEngineFoam/sprayEngineFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/sprayEngineFoam/sprayEngineFoam.C
@ -110,7 +110,7 @@ int main(int argc, char *argv[])

        if (runTime.write())
        {
-            combustion->dQ()().write();
+            combustion->Qdot()().write();
        }

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
--- a/applications/solvers/lagrangian/sprayFoam/sprayFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/sprayFoam.C
@ -102,7 +102,7 @@ int main(int argc, char *argv[])

            if (runTime.write())
            {
-                combustion->dQ()().write();
+                combustion->Qdot()().write();
            }
        }
        else
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C
@ -135,7 +135,7 @@ Foam::AnisothermalPhaseModel<BasePhaseModel>::heEqn()
            he
        )
     ==
-        this->Sh()
+        this->Qdot()
    );

    // Add the appropriate pressure-work term
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/InertPhaseModel/InertPhaseModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/InertPhaseModel/InertPhaseModel.C
@ -65,7 +65,7 @@ Foam::InertPhaseModel<BasePhaseModel>::R

 template<class BasePhaseModel>
 Foam::tmp<Foam::volScalarField>
-Foam::InertPhaseModel<BasePhaseModel>::Sh() const
+Foam::InertPhaseModel<BasePhaseModel>::Qdot() const
 {
    return tmp<volScalarField>
    (
@ -73,7 +73,7 @@ Foam::InertPhaseModel<BasePhaseModel>::Sh() const
        (
            IOobject
            (
-                IOobject::groupName("Sh", this->name()),
+                IOobject::groupName("Qdot", this->name()),
                this->mesh().time().timeName(),
                this->mesh()
            ),
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/InertPhaseModel/InertPhaseModel.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/InertPhaseModel/InertPhaseModel.H
@ -78,8 +78,8 @@ public:
                volScalarField& Yi
            ) const;

-            //- Return the reaction heat source
-            virtual tmp<volScalarField> Sh() const;
+            //- Return the heat release rate
+            virtual tmp<volScalarField> Qdot() const;
 };


--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/ReactingPhaseModel/ReactingPhaseModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/ReactingPhaseModel/ReactingPhaseModel.C
@ -105,9 +105,9 @@ Foam::ReactingPhaseModel<BasePhaseModel, ReactionType>::R

 template<class BasePhaseModel, class ReactionType>
 Foam::tmp<Foam::volScalarField>
-Foam::ReactingPhaseModel<BasePhaseModel, ReactionType>::Sh() const
+Foam::ReactingPhaseModel<BasePhaseModel, ReactionType>::Qdot() const
 {
-    return reaction_->Sh();
+    return reaction_->Qdot();
 }


--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/ReactingPhaseModel/ReactingPhaseModel.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/ReactingPhaseModel/ReactingPhaseModel.H
@ -87,8 +87,8 @@ public:
            volScalarField& Yi
        ) const;

-        //- Return the reacton heat source
-        virtual tmp<volScalarField> Sh() const;
+        //- Return heat release rate
+        virtual tmp<volScalarField> Qdot() const;
 };