MaterialLib::Fluid::WaterViscosityIAPWS Class Reference

Detailed Description

A class for viscosity model that is defined by The International Association for the Properties of Water and Steam IAPWS

With the definition, the viscosity is a function of temperature and water density

Attention

The critical enhancement, \(\bar{\mu}_2\), which is significant only within the boundaries specified by

\[ T (\mbox{in K}) \in (645.91, 650.77) \]

and

\[ \rho (\mbox{in kg m}^{-3}) \in (245.8, 405.3)\]

, is not considered.

Definition at line 42 of file WaterViscosityIAPWS.h.

#include <WaterViscosityIAPWS.h>

Inheritance diagram for MaterialLib::Fluid::WaterViscosityIAPWS:

Collaboration diagram for MaterialLib::Fluid::WaterViscosityIAPWS:

Public Member Functions
	WaterViscosityIAPWS ()=default
std::string	getName () const override
	Get model name. More...
double	getValue (const ArrayType &var_vals) const override
double	getdValue (const ArrayType &var_vals, const PropertyVariableType var_type) const override
Public Member Functions inherited from MaterialLib::Fluid::FluidProperty
virtual	~FluidProperty ()=default

Private Attributes
const double	_ref_T = 647.096
	reference temperature in K More...
const double	_ref_rho = 322.0
	reference density in kg/m^3 More...
const double	_ref_mu = 1.0e-6
	reference viscosity in Pa.s More...

Additional Inherited Members
Public Types inherited from MaterialLib::Fluid::FluidProperty
using	ArrayType = std::array< double, PropertyVariableNumber >

Constructor & Destructor Documentation

WaterViscosityIAPWS()

MaterialLib::Fluid::WaterViscosityIAPWS::WaterViscosityIAPWS ( )

default

Member Function Documentation

getdValue()

double MaterialLib::Fluid::WaterViscosityIAPWS::getdValue ( const ArrayType &  var_vals, const PropertyVariableType  var_type  ) const

overridevirtual

Get the partial differential of the density with respect to temperature or liquid pressure.

Parameters

var_vals

Variable values of temperature and water density in an array. The order of its elements is given in enum class PropertyVariableType.

Attention

The variable values must follow the SI standard, e.g temperature unit is K.

Parameters

var_type

Variable type.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 59 of file WaterViscosityIAPWS.cpp.

{
    const double bar_T =
        var_vals[static_cast<unsigned>(PropertyVariableType::T)] / _ref_T;
    const double bar_rho =
        var_vals[static_cast<unsigned>(PropertyVariableType::rho)] / _ref_rho;
 
    switch (var_type)
    {
        case PropertyVariableType::T:
            return _ref_mu * computedBarMu_dbarT(bar_T, bar_rho) / _ref_T;
        case PropertyVariableType::rho:
            return _ref_mu * computedBarMu_dbarRho(bar_T, bar_rho) / _ref_rho;
        default:
            return 0.;
    }
}

References _ref_mu, _ref_rho, _ref_T, MaterialLib::Fluid::computedBarMu_dbarRho(), MaterialLib::Fluid::computedBarMu_dbarT(), MaterialLib::Fluid::rho, and MaterialLib::Fluid::T.

getName()

std::string MaterialLib::Fluid::WaterViscosityIAPWS::getName ( ) const

inlineoverridevirtual

Get model name.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 48 of file WaterViscosityIAPWS.h.

    {
        return "IAPWS temperature-density dependent viscosity model";
    }

getValue()

double MaterialLib::Fluid::WaterViscosityIAPWS::getValue ( const ArrayType &  var_vals ) const

overridevirtual

Get density value.

Parameters

var_vals

Variable values of temperature and water density in an array. The order of its elements is given in enum class PropertyVariableType.

Attention

The variable values must follow the SI standard, e.g temperature unit is K.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 42 of file WaterViscosityIAPWS.cpp.

{
    const double bar_T =
        var_vals[static_cast<unsigned>(PropertyVariableType::T)] / _ref_T;
    const double bar_rho =
        var_vals[static_cast<unsigned>(PropertyVariableType::rho)] / _ref_rho;
 
    const double mu0 = 100. * std::sqrt(bar_T) / computeBarMu0Factor(bar_T);
 
    const auto& series_factorT = computeSeriesFactorTForMu1(bar_T);
    const auto& series_factorRho = computeSeriesFactorRhoForMu1(bar_rho);
    const double mu1 = std::exp(
        bar_rho * computeBarMu1Factor(series_factorT, series_factorRho));
 
    return mu0 * mu1 * _ref_mu;
}

References _ref_mu, _ref_rho, _ref_T, MaterialLib::Fluid::computeBarMu0Factor(), MaterialLib::Fluid::computeBarMu1Factor(), MaterialLib::Fluid::computeSeriesFactorRhoForMu1(), MaterialLib::Fluid::computeSeriesFactorTForMu1(), MaterialLib::Fluid::rho, and MaterialLib::Fluid::T.

Member Data Documentation

_ref_mu

const double MaterialLib::Fluid::WaterViscosityIAPWS::_ref_mu = 1.0e-6

private

reference viscosity in Pa.s

Definition at line 80 of file WaterViscosityIAPWS.h.

Referenced by getdValue(), and getValue().

_ref_rho

const double MaterialLib::Fluid::WaterViscosityIAPWS::_ref_rho = 322.0

private

reference density in kg/m^3

Definition at line 79 of file WaterViscosityIAPWS.h.

Referenced by getdValue(), and getValue().

_ref_T

const double MaterialLib::Fluid::WaterViscosityIAPWS::_ref_T = 647.096

private

reference temperature in K

Definition at line 78 of file WaterViscosityIAPWS.h.

Referenced by getdValue(), and getValue().

---

The documentation for this class was generated from the following files:

• MaterialLib/Fluid/Viscosity/WaterViscosityIAPWS.h
• MaterialLib/Fluid/Viscosity/WaterViscosityIAPWS.cpp