LinearPressureDependentViscosity.h

Go to the documentation of this file.

 
#pragma once
 
#include <string>
 
#include "MaterialLib/Fluid/FluidProperty.h"
 
namespace MaterialLib
{
namespace Fluid
{
class LinearPressureDependentViscosity final : public FluidProperty
{
public:
    explicit LinearPressureDependentViscosity(const double mu0,
                                              const double p0,
                                              const double gamma)
        : _mu0(mu0), _p0(p0), _gamma(gamma)
    {
    }
 
    std::string getName() const override
    {
        return "Linear pressure dependent viscosity";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
        return _mu0 * (1 + _gamma * (p - _p0));
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        (void)var_vals;
        (void)var;
        return _mu0 * _gamma;
    }
 
private:
    const double _mu0;    
    const double _p0;     
    const double _gamma;  
};
 
}  // namespace Fluid
}  // namespace MaterialLib