VogelsLiquidDynamicViscosity.h

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cmath>
#include <string>
#include <utility>
#include <vector>
 
#include "MaterialLib/Fluid/FluidProperty.h"
 
namespace MaterialLib
{
namespace Fluid
{
template <typename VogelsConstants>
class VogelsLiquidDynamicViscosity final : public FluidProperty
{
public:
    explicit VogelsLiquidDynamicViscosity(VogelsConstants constants)
        : _constants(std::move(constants))
    {
    }
 
    std::string getName() const override
    {
        return "Liquid viscosity by Vogel's equation";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        // Note: the constant of 1.e-3 is for the SI unit conversion.
        return 1.e-3 *
               std::exp(_constants.A + _constants.B / (_constants.C + T));
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        (void)var;
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        const double f_buff = _constants.B / (_constants.C + T);
        // Note: the constant of 1.e-3 is for the SI unit conversion.
        return -1.e-3 * f_buff * std::exp(_constants.A + f_buff) /
               (_constants.C + T);
    }
 
private:
    const VogelsConstants _constants;
};
 
struct VogelsViscosityConstantsWater
{
    VogelsViscosityConstantsWater() = default;
    const double A = -3.7188;
    const double B = 578.919;
    const double C = -137.546;
};
 
struct VogelsViscosityConstantsCO2
{
    VogelsViscosityConstantsCO2() = default;
    const double A = -24.0592;
    const double B = 28535.2;
    const double C = 1037.41;
};
 
struct VogelsViscosityConstantsCH4
{
    VogelsViscosityConstantsCH4() = default;
    const double A = -25.5947;
    const double B = 25392;
    const double C = 969.306;
};
 
}  // namespace Fluid
}  // namespace MaterialLib