TemperatureDependentViscosity.h

Go to the documentation of this file.

 
#pragma once
 
#include <string>
#include <cmath>
 
#include "MaterialLib/Fluid/FluidProperty.h"
 
namespace MaterialLib
{
namespace Fluid
{
class TemperatureDependentViscosity final : public FluidProperty
{
public:
    explicit TemperatureDependentViscosity(const double mu0,
                                           const double T_c,
                                           const double T_v)
        : _mu0(mu0), _temperature_c(T_c), _temperature_v(T_v)
    {
    }
 
    std::string getName() const override
    {
        return "Temperature dependent viscosity";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        return _mu0 * std::exp(-(T - _temperature_c) / _temperature_v);
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        (void)var;
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        return -_mu0 / _temperature_v *
               std::exp(-(T - _temperature_c) / _temperature_v);
    }
 
private:
    const double _mu0;            
    const double _temperature_c;  
    const double _temperature_v;  
};
 
}  // namespace Fluid
}  // namespace MaterialLib