DarcyLaw.h

Go to the documentation of this file.

 
#pragma once
 
#include "BaseLib/StrongType.h"
#include "LiquidDensity.h"
#include "LiquidViscosity.h"
#include "PermeabilityData.h"
#include "SpecificBodyForceData.h"
#include "ThermoOsmosis.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
template <int DisplacementDim>
using DarcyLawData = BaseLib::StrongType<Eigen::Vector<double, DisplacementDim>,
                                         struct DarcyLawDataTag>;
 
constexpr std::string_view ioName(struct DarcyLawDataTag*)
{
    return "velocity";
}
 
template <int DisplacementDim>
struct DarcyLawModel
{
    explicit DarcyLawModel(
        Eigen::Vector<double, DisplacementDim> const& specific_body_force)
        : b_(specific_body_force)
    {
    }
 
    void eval(CapillaryPressureData<DisplacementDim> const& p_cap_data,
              LiquidDensityData const& rho_L_data,
              LiquidViscosityData const& mu_L_data,
              PermeabilityData<DisplacementDim> const& perm_data,
              ThermoOsmosisData<DisplacementDim> const& th_osmosis_data,
              DarcyLawData<DisplacementDim>& out) const;
 
    static DarcyLawModel create(
        SpecificBodyForceData<DisplacementDim> const& specific_body_force_data)
    {
        return DarcyLawModel{specific_body_force_data.specific_body_force};
    }
 
private:
    Eigen::Vector<double, DisplacementDim> const b_;
};
 
extern template struct DarcyLawModel<2>;
extern template struct DarcyLawModel<3>;
}  // namespace ProcessLib::ThermoRichardsMechanics