PorosityFromMassBalance.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;

class PorosityFromMassBalance final : public Property
{
private:
    ParameterLib::Parameter<double> const& phi0_;
    double const phi_min_;  //< Lower limit for the porosity.
    double const phi_max_;  //< Upper limit for the porosity.
 
public:
    PorosityFromMassBalance(
        std::string name,
        ParameterLib::Parameter<double> const& initial_porosity,
        double const phi_min, double const phi_max)
        : phi0_(initial_porosity), phi_min_(phi_min), phi_max_(phi_max)
    {
        name_ = std::move(name);
    }
 
    void checkScale() const override;
 
    PropertyDataType initialValue(ParameterLib::SpatialPosition const& pos,
                                  double const t) const override
    {
        return fromVector(phi0_(t, pos));
    }
 
    PropertyDataType value(VariableArray const& variable_array,
                           ParameterLib::SpatialPosition const& pos,
                           double const t, double const dt) const override;
    PropertyDataType value(VariableArray const& variable_array,
                           VariableArray const& variable_array_prev,
                           ParameterLib::SpatialPosition const& pos,
                           double const t, double const dt) const override;
    PropertyDataType dValue(VariableArray const& variable_array,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t, double const dt) const override;
};
}  // namespace MaterialPropertyLib