interfacialarea.hh

Go to the documentation of this file.

// -*- 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/>.   *
 *****************************************************************************/
#ifndef DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_HH
#define DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_HH
 
#include <dumux/common/parameters.hh>
#include <dumux/material/fluidmatrixinteractions/fluidmatrixinteraction.hh>
#include <dumux/material/fluidmatrixinteractions/2p/efftoabsdefaultpolicy.hh>
#include <dumux/material/fluidmatrixinteractions/2p/noregularization.hh>
 
namespace Dumux::FluidMatrix {
 
template<class ScalarType,
         class BaseLaw,
         template<class> class InterfaceType,
         class Regularization = NoRegularization,
         class EffToAbsPolicy = TwoPEffToAbsDefaultPolicy>
class InterfacialArea : public Adapter<InterfacialArea<ScalarType, BaseLaw, InterfaceType, Regularization, EffToAbsPolicy>, InterfaceType>
{
public:
 
    using Scalar = ScalarType;
 
    using BasicParams = typename BaseLaw::template Params<Scalar>;
    using EffToAbsParams = typename EffToAbsPolicy::template Params<Scalar>;
    using RegularizationParams = typename Regularization::template Params<Scalar>;
 
    using EffToAbs = EffToAbsPolicy;
 
    static constexpr bool isRegularized()
    { return !std::is_same<Regularization, NoRegularization>::value; }
 
    InterfacialArea() = delete;
 
    explicit InterfacialArea(const std::string& paramGroup)
    : basicParams_(makeBasicParams(paramGroup))
    , effToAbsParams_(makeEffToAbsParams(paramGroup))
    {
        if constexpr (isRegularized())
            regularization_.init(this, paramGroup);
    }
 
    InterfacialArea(const BasicParams& baseParams,
                    const EffToAbsParams& effToAbsParams = {},
                    const RegularizationParams& regParams = {})
    : basicParams_(baseParams)
    , effToAbsParams_(effToAbsParams)
    {
        if constexpr (isRegularized())
            regularization_.init(this, baseParams, effToAbsParams, regParams);
    }
 
    template<bool enableRegularization = isRegularized()>
    Scalar area(const Scalar sw, const Scalar pc) const
    {
        const auto swe = EffToAbs::swToSwe(sw, effToAbsParams_);
        if constexpr (enableRegularization)
        {
            const auto regularized = regularization_.area(swe, pc);
            if (regularized)
                return regularized.value();
        }
 
        return BaseLaw::area(swe, pc, basicParams_);
    }
 
    template<bool enableRegularization = isRegularized()>
    Scalar darea_dpc(const Scalar sw, const Scalar pc) const
    {
        if constexpr (enableRegularization)
        {
            const auto regularized = regularization_.darea_dpc(sw, pc);
            if (regularized)
                return regularized.value();
        }
 
        return BaseLaw::darea_dpc(sw, pc, basicParams_);
    }
 
    template<bool enableRegularization = isRegularized()>
    Scalar darea_dsw(const Scalar sw, const Scalar pc) const
    {
        if constexpr (enableRegularization)
        {
            const auto regularized = regularization_.darea_dsw(sw, pc);
            if (regularized)
                return regularized.value();
        }
 
        return BaseLaw::darea_dsw(sw, pc, basicParams_);
    }
 
    bool operator== (const InterfacialArea& o) const
    {
        return basicParams_ == o.basicParams_
               && effToAbsParams_ == o.effToAbsParams_
               && regularization_ == o.regularization_;
    }
 
    static BasicParams makeBasicParams(const std::string& paramGroup)
    {
        return BaseLaw::template makeParams<Scalar>(paramGroup);
    }
 
    const BasicParams& basicParams() const
    { return basicParams_; }
 
    static EffToAbsParams makeEffToAbsParams(const std::string& paramGroup)
    {
        return EffToAbs::template makeParams<Scalar>(paramGroup);
    }
 
    const EffToAbsParams& effToAbsParams() const
    { return effToAbsParams_; }
 
private:
    BasicParams basicParams_;
    EffToAbsParams effToAbsParams_;
    Regularization regularization_;
};
 
} // end namespace Dumux::FluidMatrix
 
#endif