BUG: solidIsothermalReactionRate: redefine the latent energy (fixes #2646)

src/thermophysicalModels/solidSpecie/reaction/reactionRate/solidIsothermalReactionRate/solidIsothermalReactionRate.H

@@ -5,7 +5,7 @@
     \\  /    A nd           | www.openfoam.com
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
-    Copyright (C) 2021 OpenCFD Ltd.
+    Copyright (C) 2021-2022 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -29,6 +29,32 @@ Class
 Description
     Isothermal reaction rate for solids
 
+Usage
+    Minimal example:
+    \verbatim
+    {
+        // Mandatory entries
+        C       <scalar>;
+        Cp      <scalar>;
+        Tpc     <scalar>;
+        Elat    <scalar>;
+    }
+    \endverbatim
+
+    where the entries mean:
+    \table
+      Property  | Description             | Type   | Reqd | Deflt
+      C         | Model constant          | scalar | yes  | -
+      Cp | Specific heat of solid material (assumed constant) | scalar | yes | -
+      Tpc       | Phase change temperature | scalar | yes | -
+      Elat      | Latent heat of phase change | scalar | yes | -
+    \endtable
+
+Note
+  - \c Elat, the latent energy is defined in the present
+    implementation as a negative value (using the convention
+    endothermic process has negative enthalpy.)
+
 SourceFiles
     solidIsothermalReactionRateI.H
 

src/thermophysicalModels/solidSpecie/reaction/reactionRate/solidIsothermalReactionRate/solidIsothermalReactionRateI.H

@@ -5,7 +5,7 @@
     \\  /    A nd           | www.openfoam.com
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
-    Copyright (C) 2021 OpenCFD Ltd.
+    Copyright (C) 2021-2022 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -51,7 +51,7 @@ inline Foam::solidIsothermalReactionRate::solidIsothermalReactionRate
     C_(dict.get<scalar>("C")),
     Cp_(dict.get<scalar>("Cp")),
     Tpc_(dict.get<scalar>("Tpc")),
-    Elat_(dict.getCheck<scalar>("Elat", scalarMinMax::ge(SMALL)))
+    Elat_(dict.get<scalar>("Elat"))
 {}