MaterialLib::Fluid::LiquidDensity Class Reference

Detailed Description

Class of a generic density model for liquid fluids varying with temperature or pressure.

The formula is given on engineeringtoolbox which reads

$$\rho_l = \rho_0/(1+\beta(T-T_0))/(1-(p-p_0)/E)$$

where

$$\begin{eqnarray*} &\rho_l:& \mbox{liquid density,}\\ &\rho_0:& \mbox{initial liquid density,}\\ &\beta: & \mbox{volumetric temperature expansion coefficient,}\\ &T:& \mbox{temperature,}\\ &T_0:& \mbox{initial temperature,}\\ &p:& \mbox{pressure,}\\ &p_0:& \mbox{initial pressure,}\\ &E:& \mbox{bulk modulus fluid elasticity.}\\ \end{eqnarray*}$$

Definition at line 51 of file LiquidDensity.h.

#include <LiquidDensity.h>

Inheritance diagram for MaterialLib::Fluid::LiquidDensity:

Collaboration diagram for MaterialLib::Fluid::LiquidDensity:

Public Member Functions
	LiquidDensity (const double beta, const double rho0, const double T0, const double p0, const double E)
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) const override
Public Member Functions inherited from MaterialLib::Fluid::FluidProperty
virtual	~FluidProperty ()=default

Private Member Functions
double	dLiquidDensity_dT (const double T, const double p) const
double	dLiquidDensity_dp (const double T, const double p) const

Private Attributes
const double	_beta
	Volumetric temperature expansion coefficient. More...
const double	_rho0
	Initial density. More...
const double	_temperature0
	Initial temperature. More...
const double	_p0
	Initial pressure. More...
const double	_bulk_modulus
	Bulk modulus. More...

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

Constructor & Destructor Documentation

LiquidDensity()

MaterialLib::Fluid::LiquidDensity::LiquidDensity ( const double  beta, const double  rho0, const double  T0, const double  p0, const double  E  )

inlineexplicit

Definition at line 54 of file LiquidDensity.h.

        : _beta(beta), _rho0(rho0), _temperature0(T0), _p0(p0), _bulk_modulus(E)
    {
    }

Member Function Documentation

dLiquidDensity_dp()

double MaterialLib::Fluid::LiquidDensity::dLiquidDensity_dp ( const double  T, const double  p  ) const

inlineprivate

Calculate the derivative of density of fluids with respect to pressure.

Parameters

T	Temperature (K).
p	Pressure (Pa).

Definition at line 132 of file LiquidDensity.h.

    {
        const double fac_p = 1. - (p - _p0) / _bulk_modulus;
        return _rho0 / (1. + _beta * (T - _temperature0)) /
               (fac_p * fac_p * _bulk_modulus);
    }

References _beta, _bulk_modulus, _p0, _rho0, _temperature0, MaterialLib::Fluid::p, and MaterialLib::Fluid::T.

Referenced by getdValue().

dLiquidDensity_dT()

double MaterialLib::Fluid::LiquidDensity::dLiquidDensity_dT ( const double  T, const double  p  ) const

inlineprivate

Calculate the derivative of density of fluids with respect to temperature.

Parameters

T	Temperature (K).
p	Pressure (Pa).

Definition at line 120 of file LiquidDensity.h.

    {
        const double fac_T = 1. + _beta * (T - _temperature0);
        return -_beta * _rho0 / (fac_T * fac_T) /
               (1. - (p - _p0) / _bulk_modulus);
    }

References _beta, _bulk_modulus, _p0, _rho0, _temperature0, MaterialLib::Fluid::p, and MaterialLib::Fluid::T.

Referenced by getdValue().

getdValue()

double MaterialLib::Fluid::LiquidDensity::getdValue ( const ArrayType &  var_vals, const PropertyVariableType  var  ) const

inlineoverridevirtual

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

Parameters

var_vals	Variable values in an array. The order of its elements is given in enum class PropertyVariableType.
var	Variable type.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 82 of file LiquidDensity.h.

    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
        switch (var)
        {
            case PropertyVariableType::T:
                return dLiquidDensity_dT(T, p);
            case PropertyVariableType::p:
                return dLiquidDensity_dp(T, p);
            default:
                return 0.;
        }
    }

References dLiquidDensity_dp(), dLiquidDensity_dT(), MaterialLib::Fluid::p, and MaterialLib::Fluid::T.

getName()

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

inlineoverridevirtual

Get density model name.

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 61 of file LiquidDensity.h.

    {
        return "Temperature or pressure dependent liquid density";
    }

getValue()

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

inlineoverridevirtual

Get density value.

Parameters

var_vals

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

Implements MaterialLib::Fluid::FluidProperty.

Definition at line 69 of file LiquidDensity.h.

    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
        return _rho0 / (1. + _beta * (T - _temperature0)) /
               (1. - (p - _p0) / _bulk_modulus);
    }

References _beta, _bulk_modulus, _p0, _rho0, _temperature0, MaterialLib::Fluid::p, and MaterialLib::Fluid::T.

Member Data Documentation

_beta

const double MaterialLib::Fluid::LiquidDensity::_beta

private

Volumetric temperature expansion coefficient.

Definition at line 100 of file LiquidDensity.h.

Referenced by dLiquidDensity_dp(), dLiquidDensity_dT(), and getValue().

_bulk_modulus

const double MaterialLib::Fluid::LiquidDensity::_bulk_modulus

private

Bulk modulus.

Definition at line 112 of file LiquidDensity.h.

Referenced by dLiquidDensity_dp(), dLiquidDensity_dT(), and getValue().

_p0

const double MaterialLib::Fluid::LiquidDensity::_p0

private

Initial pressure.

Definition at line 109 of file LiquidDensity.h.

Referenced by dLiquidDensity_dp(), dLiquidDensity_dT(), and getValue().

_rho0

const double MaterialLib::Fluid::LiquidDensity::_rho0

private

Initial density.

Definition at line 103 of file LiquidDensity.h.

Referenced by dLiquidDensity_dp(), dLiquidDensity_dT(), and getValue().

_temperature0

const double MaterialLib::Fluid::LiquidDensity::_temperature0

private

Initial temperature.

Definition at line 106 of file LiquidDensity.h.

Referenced by dLiquidDensity_dp(), dLiquidDensity_dT(), and getValue().

---

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

• MaterialLib/Fluid/Density/LiquidDensity.h