VapourDiffusionPMQ.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class Phase;
 
class VapourDiffusionPMQ final : public Property
{
public:
    explicit VapourDiffusionPMQ(std::string name,
                                double const base_diffusion_coefficient,
                                double const exponent)
        : base_diffusion_coefficient_(base_diffusion_coefficient),
          exponent_(exponent)
    {
        name_ = std::move(name);
    }
 
    void checkScale() const override
    {
        if (!(std::holds_alternative<Phase*>(scale_) ||
              std::holds_alternative<Component*>(scale_)))
        {
            OGS_FATAL(
                "The property 'VapourDiffusionPMQ' is "
                "implemented on the 'phase' and 'component' scale only.");
        }
    }
 
    PropertyDataType value(VariableArray const& variable_array,
                           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;
 
private:
    double const base_diffusion_coefficient_;
    double const exponent_;
};
 
}  // namespace MaterialPropertyLib