[fluidmatrixinteraction] remove porosityreactivebed

This is identical to PorosityPrecipitation with a unit reference porosity.

dumux/material/fluidmatrixinteractions/CMakeLists.txt

@@ -12,5 +12,4 @@ diffusivitymillingtonquirk.hh
 permeabilitykozenycarman.hh
 permeabilityrutqvisttsang.hh
 porosityprecipitation.hh
-porosityreactivebed.hh
 DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/material/fluidmatrixinteractions)

dumux/material/fluidmatrixinteractions/porosityreactivebed.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 2 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 Fluidmatrixinteractions
- * \brief Class for the evaluation of the porosity subject to precipitation.
- */
-#ifndef DUMUX_POROSITY_REACTIVE_BED_HH
-#define DUMUX_POROSITY_REACTIVE_BED_HH
-
-#include <dumux/common/properties.hh>
-#include <dumux/discretization/evalsolution.hh>
-
-namespace Dumux {
-
-/*!
- * \ingroup Fluidmatrixinteractions
- * \brief Calculates the porosity depeding on the volume fractions of different solid species.
- */
-template<class TypeTag>
-class PorosityReactiveBed
-{
-    using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar);
-    using Problem = typename GET_PROP_TYPE(TypeTag, Problem);
-    using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
-    using SpatialParams = typename GET_PROP_TYPE(TypeTag, SpatialParams);
-    using FVElementGeometry = typename GET_PROP_TYPE(TypeTag, FVGridGeometry)::LocalView;
-    using SubControlVolume = typename FVElementGeometry::SubControlVolume;
-
-    static const int dim = GridView::dimension;
-    static const int dimWorld = GridView::dimensionworld;
-    static const int numComponents = GET_PROP_TYPE(TypeTag, ModelTraits)::numComponents();
-    static const int numSolidPhases = GET_PROP_TYPE(TypeTag, ModelTraits)::numSPhases();
-
-    using Element = typename GridView::template Codim<0>:: Entity;
-
-public:
-    void init(const SpatialParams& spatialParams)
-    {
-        spatialParamsPtr_ = &spatialParams;
-    }
-
-    /*!
-     * \brief Calculates the porosity in a sub-control volume.
-     *
-     * \param element element
-     * \param elemSol the element solution
-     * \param scv sub control volume
-     */
-    template<class ElementSolution>
-    Scalar evaluatePorosity(const Element& element,
-                            const SubControlVolume& scv,
-                            const ElementSolution& elemSol) const
-    {
-        auto priVars = evalSolution(element, element.geometry(), elemSol, scv.center());
-
-        Scalar sumPrecipitates = 0.0;
-        for (unsigned int solidPhaseIdx = 0; solidPhaseIdx < numSolidPhases; ++solidPhaseIdx)
-            sumPrecipitates += priVars[numComponents + solidPhaseIdx];
-
-        return (1 - sumPrecipitates);
-    }
-
-private:
-    const SpatialParams& spatialParams_() const
-    { return *spatialParamsPtr_; }
-
-    const SpatialParams* spatialParamsPtr_;
-};
-
-} // namespace Dumux
-
-#endif