Detailed Description

Water density model base on the IAPWS Industrial Formulation 1997 IF97-Rev

Definition at line 28 of file WaterDensityIAPWSIF97Region1.h.

#include <WaterDensityIAPWSIF97Region1.h>

Inheritance diagram for MaterialLib::Fluid::WaterDensityIAPWSIF97Region1:

Collaboration diagram for MaterialLib::Fluid::WaterDensityIAPWSIF97Region1:

Public Member Functions
	WaterDensityIAPWSIF97Region1 ()

std::string	getName () const override
	Get density 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 DimensionLessGibbsFreeEnergyRegion1	_gibbs_free_energy

const double	_ref_T = 1386
	reference temperature in K. More...

const double	_ref_p = 1.653e7

const double	_sR = 461.526
	Specific water vapour gas constant in J/(kgK). More...

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

Constructor & Destructor Documentation

◆ WaterDensityIAPWSIF97Region1()

MaterialLib::Fluid::WaterDensityIAPWSIF97Region1::WaterDensityIAPWSIF97Region1 ( )

inline

Definition at line 31 of file WaterDensityIAPWSIF97Region1.h.

32 : _gibbs_free_energy(DimensionLessGibbsFreeEnergyRegion1()){};

MaterialLib::Fluid::WaterDensityIAPWSIF97Region1::_gibbs_free_energy

const DimensionLessGibbsFreeEnergyRegion1 _gibbs_free_energy

Definition: WaterDensityIAPWSIF97Region1.h:85

Member Function Documentation

◆ getdValue()

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

inlineoverridevirtual

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

Parameters

var_vals	Variable values in an array of temperature and pressure.
var_type	Variable type.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 58 of file WaterDensityIAPWSIF97Region1.h.

     {
         const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
         const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
  
         const double tau = _ref_T / T;
         const double pi = p / _ref_p;
  
         const double dgamma_dpi = _gibbs_free_energy.get_dgamma_dpi(tau, pi);
         switch (var_type)
         {
             case PropertyVariableType::T:
                 return -(_ref_p -
                          tau * _ref_p *
                              _gibbs_free_energy.get_dgamma_dtau_dpi(tau, pi) /
                              dgamma_dpi) /
                        (_sR * T * T * dgamma_dpi);
             case PropertyVariableType::p:
                 return -_gibbs_free_energy.get_dgamma_dpi_dpi(tau, pi) /
                        (_sR * T * dgamma_dpi * dgamma_dpi);
             default:
                 return 0.0;
         }
     }

References _gibbs_free_energy, _ref_p, _ref_T, _sR, MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1::get_dgamma_dpi(), MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1::get_dgamma_dpi_dpi(), MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1::get_dgamma_dtau_dpi(), MaterialLib::Fluid::p, and MaterialLib::Fluid::T.

◆ getName()

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

inlineoverridevirtual

Get density model name.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 35 of file WaterDensityIAPWSIF97Region1.h.

     {
         return "Water density model based on IAPWSIF97";
     }

◆ getValue()

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

inlineoverridevirtual

Get density value.

Parameters

var_vals Variable values in an array of temperature and pressure.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 42 of file WaterDensityIAPWSIF97Region1.h.

     {
         const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
         const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
         const double tau = _ref_T / T;
         const double pi = p / _ref_p;
  
         return _ref_p / (_sR * T * _gibbs_free_energy.get_dgamma_dpi(tau, pi));
     }