[nonisothermal] add fv-type operators

dumux/porousmediumflow/nonisothermal/CMakeLists.txt

@@ -3,5 +3,6 @@ indices.hh
 iofields.hh
 localresidual.hh
 model.hh
+operators.hh
 volumevariables.hh
 DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/porousmediumflow/nonisothermal)

dumux/porousmediumflow/nonisothermal/operators.hh

+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ *   See the file COPYING for full copying permissions.                      *
+ *                                                                           *
+ *   This program is free software: you can redistribute it and/or modify    *
+ *   it under the terms of the GNU General Public License as published by    *
+ *   the Free Software Foundation, either version 3 of the License, or       *
+ *   (at your option) any later version.                                     *
+ *                                                                           *
+ *   This program is distributed in the hope that it will be useful,         *
+ *   but WITHOUT ANY WARRANTY; without even the implied warranty of          *
+ *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the            *
+ *   GNU General Public License for more details.                            *
+ *                                                                           *
+ *   You should have received a copy of the GNU General Public License       *
+ *   along with this program.  If not, see <http://www.gnu.org/licenses/>.   *
+ *****************************************************************************/
+/*!
+ * \file
+ * \ingroup NIModel
+ * \brief Sub-control entity-local evaluation of the operators
+ *        involved in the system of equations non-isothermal models
+ *        that assume thermal equilibrium between all phases.
+ */
+#ifndef DUMUX_FV_NON_ISOTHERMAL_OPERATORS_HH
+#define DUMUX_FV_NON_ISOTHERMAL_OPERATORS_HH
+
+namespace Dumux::Experimental {
+
+/*!
+ * \ingroup NIModel
+ * \brief Sub-control entity-local evaluation of the operators
+ *        involved in the system of equations of non-isothermal models
+ *        that assume thermal equilibrium between all phases.
+ * \tparam ModelTraits Model-specific traits.
+ * \tparam LocalContext Element-local context (geometry & primary/secondary variables)
+ */
+template<class ModelTraits, class LocalContext>
+class FVNonIsothermalOperators
+{
+    // The variables required for the evaluation of the equation
+    using ElementVariables = typename LocalContext::ElementVariables;
+
+    // The grid geometry on which the scheme operates
+    using GridGeometry = typename LocalContext::ElementGridGeometry::GridGeometry;
+    using SubControlVolume = typename GridGeometry::SubControlVolume;
+
+public:
+
+    /*!
+     * \brief The energy storage in the fluid phase with index phaseIdx
+     *
+     * \param storage The mass of the component within the sub-control volume
+     * \param scv The sub-control volume
+     * \param context The element-local context (primary/secondary variables)
+     * \param phaseIdx The phase index
+     */
+    template<class NumEqVector>
+    static void addFluidPhaseStorage(NumEqVector& storage,
+                                     const SubControlVolume& scv,
+                                     const LocalContext& context,
+                                     int phaseIdx)
+    {
+        // do no-op in case energy balance is not active
+        if constexpr (ModelTraits::enableEnergyBalance())
+        {
+            const auto& volVars = context.elementVariables().elemVolVars()[scv];
+            storage[ModelTraits::Indices::energyEqIdx] += volVars.porosity()
+                                                          *volVars.density(phaseIdx)
+                                                          *volVars.internalEnergy(phaseIdx)
+                                                          *volVars.saturation(phaseIdx);
+        }
+    }
+
+    /*!
+     * \brief The energy storage in the solid matrix
+     *
+     * \param storage The mass of the component within the sub-control volume
+     * \param scv The sub-control volume
+     * \param context The element-local context (primary/secondary variables)
+     */
+    template<class NumEqVector>
+    static void addSolidPhaseStorage(NumEqVector& storage,
+                                     const SubControlVolume& scv,
+                                     const LocalContext& context)
+    {
+        // do no-op in case energy balance is not active
+        if constexpr (ModelTraits::enableEnergyBalance())
+        {
+            const auto& volVars = context.elementVariables().elemVolVars()[scv];
+            storage[ModelTraits::Indices::energyEqIdx] += volVars.temperature()
+                                                          *volVars.solidHeatCapacity()
+                                                          *volVars.solidDensity()
+                                                          *(1.0 - volVars.porosity());
+        }
+    }
+
+    /*!
+     * \brief The advective phase energy fluxes
+     *
+     * \param flux The flux
+     * \param fluxVars The flux variables.
+     * \param phaseIdx The phase index
+     */
+    template<class NumEqVector, class FluxVariables>
+    static void addHeatConvectionFlux(NumEqVector& flux,
+                                      FluxVariables& fluxVars,
+                                      int phaseIdx)
+    {
+        // do no-op in case energy balance is not active
+        if constexpr (ModelTraits::enableEnergyBalance())
+        {
+            auto upwindTerm = [phaseIdx](const auto& volVars)
+            { return volVars.density(phaseIdx)*volVars.mobility(phaseIdx)*volVars.enthalpy(phaseIdx); };
+
+            flux[ModelTraits::Indices::energyEqIdx] += fluxVars.advectiveFlux(phaseIdx, upwindTerm);
+        }
+    }
+
+    /*!
+     * \brief The diffusive energy fluxes
+     *
+     * \param flux The flux
+     * \param fluxVars The flux variables.
+     */
+    template<class NumEqVector, class FluxVariables>
+    static void addHeatConductionFlux(NumEqVector& flux,
+                                      FluxVariables& fluxVars)
+    {
+        // do no-op in case energy balance is not active
+        if constexpr (ModelTraits::enableEnergyBalance())
+            flux[ModelTraits::Indices::energyEqIdx] += fluxVars.heatConductionFlux();
+    }
+
+    /*!
+     * \brief heat transfer between the phases for nonequilibrium models
+     *
+     * \param source The source which ought to be simulated
+     * \param element An element which contains part of the control volume
+     * \param fvGeometry The finite-volume geometry
+     * \param context The element-local context (primary/secondary variables)
+     * \param scv The sub-control volume over which we integrate the source term
+     */
+    template<class NumEqVector>
+    static void addEnergySource(NumEqVector& source,
+                                const LocalContext& context,
+                                const SubControlVolume &scv)
+    {}
+};
+
+} // end namespace Dumux::Experimental
+
+#endif