Class for fluid properties of an arbitrary vapor mixture. More...

#include <IdealRealGasMixtureFluidProperties.h>

Inheritance diagram for IdealRealGasMixtureFluidProperties:

Public Member Functions
	IdealRealGasMixtureFluidProperties (const InputParameters &parameters)

virtual unsigned int	getNumberOfSecondaryVapors () const override
	Number of secondary vapors (non-condensable components) More...

virtual Real	p_from_v_e (Real v, Real e, const std::vector< Real > &x) const override
	Pressure from specific volume and specific internal energy. More...

virtual void	p_from_v_e (Real v, Real e, const std::vector< Real > &x, Real &p, Real &dp_dv, Real &dp_de, std::vector< Real > &dp_dx) const override
	Pressure and its derivatives from specific volume and specific internal energy. More...

virtual Real	T_from_v_e (Real v, Real e, const std::vector< Real > &x) const override
	Temperature from specific volume and specific internal energy. More...

virtual void	T_from_v_e (Real v, Real e, const std::vector< Real > &x, Real &T, Real &dT_dv, Real &dT_de, std::vector< Real > &dT_dx) const override
	Temperature and its derivatives from specific volume and specific internal energy. More...

virtual Real	rho_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Density from pressure and temperature. More...

virtual void	rho_from_p_T (Real p, Real T, const std::vector< Real > &x, Real &rho, Real &drho_dp, Real &drho_dT, std::vector< Real > &drho_dx) const override
	Density and its derivatives from pressure and temperature. More...

virtual Real	e_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Specific internal energy from pressure and temperature. More...

virtual void	e_from_p_T (Real p, Real T, const std::vector< Real > &x, Real &e, Real &de_dp, Real &de_dT, std::vector< Real > &de_dx) const override
	Specific internal energy and its derivatives from pressure and temperature. More...

virtual Real	c_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Speed of sound from pressure and temperature. More...

virtual void	c_from_p_T (Real p, Real T, const std::vector< Real > &x, Real &c, Real &dc_dp, Real &dc_dT, std::vector< Real > &dc_dx) const override
	Speed of sound and its derivatives from pressure and temperature. More...

virtual Real	cp_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Isobaric heat capacity from pressure and temperature. More...

virtual Real	cv_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Isochoric heat capacity from pressure and temperature. More...

virtual Real	mu_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Dynamic viscosity from pressure and temperature. More...

virtual Real	k_from_p_T (Real p, Real T, const std::vector< Real > &x) const override
	Thermal conductivity from pressure and temperature. More...

virtual Real	v_from_p_T (Real p, Real T, const std::vector< Real > &x) const
	Specific volume from pressure and temperature. More...

virtual void	v_from_p_T (Real p, Real T, const std::vector< Real > &x, Real &v, Real &dv_dp, Real &dv_dT, std::vector< Real > &dv_dx) const
	Specific volume and its derivatives from pressure and temperature. More...

void	p_T_from_v_e (Real v, Real e, const std::vector< Real > &x, Real &p, Real &T) const
	Pressure and temperature from specific volume and specific internal energy. More...

void	p_T_from_v_e (Real v, Real e, const std::vector< Real > &x, Real &p, Real &dp_dv, Real &dp_de, std::vector< Real > &dp_dx, Real &T, Real &dT_dv, Real &dT_de, std::vector< Real > &dT_dx) const
	Pressure and temperature from specific volume and specific internal energy. More...

Real	p_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Pressure from temperature and specific volume. More...

void	p_from_T_v (Real T, Real v, const std::vector< Real > &x, Real &p, Real &dp_dT, Real &dp_dv) const
	Pressure and its derivatives from temperature and specific volume. More...

void	p_from_T_v (Real T, Real v, const std::vector< Real > &x, Real &p, Real &dp_dT, Real &dp_dv, std::vector< Real > &dp_dx) const
	Pressure and its derivatives from temperature and specific volume. More...

Real	e_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Specific internal energy from temperature and specific volume. More...

void	e_from_T_v (Real T, Real v, const std::vector< Real > &x, Real &e, Real &de_dT, Real &de_dv, std::vector< Real > &de_dx) const
	Specific internal energy and its derivatives from temperature and specific volume. More...

void	s_from_T_v (Real T, Real v, const std::vector< Real > &x, Real &s, Real &ds_dT, Real &ds_dv) const
	Specific entropy and its derivatives from temperature and specific volume. More...

Real	c_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Speed of sound from temperature and specific volume. More...

void	c_from_T_v (Real T, Real v, const std::vector< Real > &x, Real &c, Real &dc_dT, Real &dc_dv, std::vector< Real > &dc_dx) const
	Speed of sound and its derivatives from temperature and specific volume. More...

Real	cp_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Isobaric heat capacity from temperature and specific volume. More...

Real	cv_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Isochoric heat capacity from temperature and specific volume. More...

Real	mu_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Dynamic viscosity from temperature and specific volume. More...

Real	k_from_T_v (Real T, Real v, const std::vector< Real > &x) const
	Thermal conductivity from temperature and specific volume. More...

Real	xs_prim_from_p_T (Real p, Real T, const std::vector< Real > &x) const
	Mass fraction of primary (condensable) component at saturation from pressure and temperature. More...

Real	primaryMassFraction (const std::vector< Real > &x) const
	Computes the mass fraction of the primary vapor given mass fractions of the secondary vapors. More...

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

Protected Types
enum	NaNMessage { NAN_MESSAGE_NONE = 0, NAN_MESSAGE_WARNING = 1, NAN_MESSAGE_ERROR = 2 }

Protected Member Functions
Real	getNaN () const
	Produces errors, warnings, or just quiet NaNs. More...

template<typename... Args>
Real	getNaN (Args &&... args) const

Protected Attributes
const SinglePhaseFluidProperties *	_fp_primary
	Primary vapor fluid properties. More...

std::vector< const SinglePhaseFluidProperties * >	_fp_secondary
	Secondary vapor fluid properties. More...

const std::vector< UserObjectName >	_fp_secondary_names
	Names of secondary vapor fluid properties. More...

const unsigned int	_n_secondary_vapors
	Number of secondary vapors. More...

const Real	_T_mix_max
	maximum temperature of all components More...

const MooseObject *	_moose_object

enum NaNMessage	_emit_on_nan
	Raise mooseWarning or mooseError? More...

Static Protected Attributes
static constexpr const Real	R_molar = 8.3144598
	molar (or universal) gas constant More...

Detailed Description

Class for fluid properties of an arbitrary vapor mixture.

Definition at line 27 of file IdealRealGasMixtureFluidProperties.h.

Member Enumeration Documentation

◆ NaNMessage

enum NaNInterface::NaNMessage

protectedinherited

Enumerator
NAN_MESSAGE_NONE
NAN_MESSAGE_WARNING
NAN_MESSAGE_ERROR

Definition at line 34 of file NaNInterface.h.

   {
     NAN_MESSAGE_NONE = 0,
     NAN_MESSAGE_WARNING = 1,
     NAN_MESSAGE_ERROR = 2
   };

Constructor & Destructor Documentation

◆ IdealRealGasMixtureFluidProperties()

IdealRealGasMixtureFluidProperties::IdealRealGasMixtureFluidProperties ( const InputParameters & parameters )

Definition at line 39 of file IdealRealGasMixtureFluidProperties.C.

   : VaporMixtureFluidProperties(parameters),
     NaNInterface(this),
     _fp_secondary_names(getParam<std::vector<UserObjectName>>("fp_secondary")),
     _n_secondary_vapors(_fp_secondary_names.size()),
     _T_mix_max(getParam<Real>("_T_mix_max"))
 {
   _fp_primary = &getUserObject<SinglePhaseFluidProperties>("fp_primary");
 
   _fp_secondary.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     _fp_secondary[i] = &getUserObjectByName<SinglePhaseFluidProperties>(_fp_secondary_names[i]);
 }

Member Function Documentation

◆ c_from_p_T() [1/2]

Real IdealRealGasMixtureFluidProperties::c_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Speed of sound from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: speed of sound

Implements VaporMixtureFluidProperties.

Definition at line 224 of file IdealRealGasMixtureFluidProperties.C.

Referenced by c_from_p_T().

 {
   Real v;
   Real p_unused, dp_dT, dp_dv;
   Real s, ds_dT, ds_dv;
 
   v = v_from_p_T(p, T, x);
   p_from_T_v(T, v, x, p_unused, dp_dT, dp_dv);
   s_from_T_v(T, v, x, s, ds_dT, ds_dv);
 
   Real dp_dv_s = dp_dv - dp_dT * ds_dv / ds_dT;
 
   if (dp_dv_s >= 0)
     mooseWarning(name(), ":c_from_p_T(), dp_dv_s = ", dp_dv_s, ". Should be negative.");
   return v * std::sqrt(-dp_dv_s);
 }

◆ c_from_p_T() [2/2]

void IdealRealGasMixtureFluidProperties::c_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x,
		Real &	c,
		Real &	dc_dp,
		Real &	dc_dT,
		std::vector< Real > &	dc_dx
	)		const

overridevirtual

Speed of sound and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values
[out]	c	speed of sound
[out]	dc_dp	derivative of speed of sound w.r.t. pressure
[out]	dc_dT	derivative of speed of sound w.r.t. temperature
[out]	dc_dx	derivative of speed of sound w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 242 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p_perturbed, T_perturbed, c_perturbed;
 
   c = c_from_p_T(p, T, x);
   // For derived properties, we would need higher order derivatives;
   // therefore, numerical derivatives are used here
   p_perturbed = p * 1.e-6;
   c_perturbed = c_from_p_T(p_perturbed, T, x);
   dc_dp = (c_perturbed - c) / (p_perturbed - p);
 
   Real dT = 1.e-6;
   T_perturbed = T + dT;
   c_perturbed = c_from_p_T(p, T_perturbed, x);
   dc_dT = (c_perturbed - c) / (dT);
 
   dc_dx.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real c_perturbed;
     std::vector<Real> x_perturbed(x);
     Real dx_i = 1e-6;
     for (unsigned int j = 0; j < _n_secondary_vapors; j++)
     {
       if (j != i)
         x_perturbed[j] =
             x[j] * (1.0 - (x[i] + dx_i)) / (1.0 - x[i]); // recalculate new mass fractions
     }
     x_perturbed[i] += dx_i;
     c_perturbed = c_from_p_T(p, T, x_perturbed);
     dc_dx[i] = ((c_perturbed - c) / dx_i);
   }
 }

◆ c_from_T_v() [1/2]

Real IdealRealGasMixtureFluidProperties::c_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Speed of sound from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: speed of sound

Definition at line 620 of file IdealRealGasMixtureFluidProperties.C.

Referenced by c_from_T_v().

 {
   Real p, dp_dT, dp_dv;
   Real s, ds_dT, ds_dv;
 
   p_from_T_v(T, v, x, p, dp_dT, dp_dv);
   s_from_T_v(T, v, x, s, ds_dT, ds_dv);
 
   Real dp_dv_s = dp_dv - dp_dT * ds_dv / ds_dT;
 
   if (dp_dv_s >= 0)
     mooseWarning(name(), ":c_from_T_v(), dp_dv_s = ", dp_dv_s, ". Should be negative.");
   return v * std::sqrt(-dp_dv_s);
 }

◆ c_from_T_v() [2/2]

void IdealRealGasMixtureFluidProperties::c_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x,
		Real &	c,
		Real &	dc_dT,
		Real &	dc_dv,
		std::vector< Real > &	dc_dx
	)		const

Speed of sound and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values
[out]	c	speed of sound
[out]	dc_dT	derivative of speed of sound w.r.t. temperature
[out]	dc_dv	derivative of speed of sound w.r.t. specific volume
[out]	dc_dx	derivative of speed of sound w.r.t. vapor mass fraction values

Definition at line 636 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real T_perturbed, v_perturbed, c_perturbed;
 
   c = c_from_T_v(T, v, x);
   // For derived properties, we would need higher order derivatives;
   // therefore, numerical derivatives are used here.
   Real dT = 1.e-6;
   T_perturbed = T + dT;
   c_perturbed = c_from_T_v(T_perturbed, v, x);
   dc_dT = (c_perturbed - c) / (dT);
 
   v_perturbed = v * 1.e-6;
   c_perturbed = c_from_T_v(T, v_perturbed, x);
   dc_dv = (c_perturbed - c) / (v_perturbed - v);
 
   dc_dx.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real c_perturbed;
     std::vector<Real> x_perturbed(x);
     Real dx_i = 1e-6;
     for (unsigned int j = 0; j < _n_secondary_vapors; j++)
     {
       if (j != i)
         x_perturbed[j] =
             x[j] * (1.0 - (x[i] + dx_i)) / (1.0 - x[i]); // recalculate new mass fractions
     }
     x_perturbed[i] += dx_i;
     c_perturbed = c_from_T_v(T, v, x_perturbed);
     dc_dx[i] = ((c_perturbed - c) / dx_i);
   }
 }

◆ cp_from_p_T()

Real IdealRealGasMixtureFluidProperties::cp_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Isobaric heat capacity from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: isobaric heat capacity

Implements VaporMixtureFluidProperties.

Definition at line 283 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p_unused, dp_dT, dp_dv;
   Real h, dh_dT, dh_dv;
 
   Real v = v_from_p_T(p, T, x);
   p_from_T_v(T, v, x, p_unused, dp_dT, dp_dv);
 
   const Real x_primary = primaryMassFraction(x);
 
   _fp_primary->h_from_T_v(T, v / x_primary, h, dh_dT, dh_dv);
   Real cp = x_primary * (dh_dT - dh_dv * dp_dT / dp_dv);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->h_from_T_v(T, v / x[i], h, dh_dT, dh_dv);
     cp += x[i] * (dh_dT - dh_dv * dp_dT / dp_dv);
   }
 
   return cp;
 }

◆ cp_from_T_v()

Real IdealRealGasMixtureFluidProperties::cp_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Isobaric heat capacity from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: isobaric heat capacity

Definition at line 677 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p, dp_dT, dp_dv;
   Real h, dh_dT, dh_dv;
 
   p_from_T_v(T, v, x, p, dp_dT, dp_dv);
 
   const Real x_primary = primaryMassFraction(x);
 
   _fp_primary->h_from_T_v(T, v / x_primary, h, dh_dT, dh_dv);
   Real cp = x_primary * (dh_dT - dh_dv * dp_dT / dp_dv);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->h_from_T_v(T, v / x[i], h, dh_dT, dh_dv);
     cp += x[i] * (dh_dT - dh_dv * dp_dT / dp_dv);
   }
 
   return cp;
 }

◆ cv_from_p_T()

Real IdealRealGasMixtureFluidProperties::cv_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Isochoric heat capacity from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: isochoric heat capacity

Implements VaporMixtureFluidProperties.

Definition at line 306 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = v_from_p_T(p, T, x);
 
   const Real x_primary = primaryMassFraction(x);
   Real cv = x_primary * _fp_primary->cv_from_T_v(T, v / x_primary);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     cv += x[i] * _fp_secondary[i]->cv_from_T_v(T, v / x[i]);
 
   return cv;
 }

◆ cv_from_T_v()

Real IdealRealGasMixtureFluidProperties::cv_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Isochoric heat capacity from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: isochoric heat capacity

Definition at line 699 of file IdealRealGasMixtureFluidProperties.C.

Referenced by cv_from_p_T().

 {
   const Real x_primary = primaryMassFraction(x);
   Real cv = x_primary * _fp_primary->cv_from_T_v(T, v / x_primary);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     cv += x[i] * _fp_secondary[i]->cv_from_T_v(T, v / x[i]);
 
   return cv;
 }

◆ e_from_p_T() [1/2]

Real IdealRealGasMixtureFluidProperties::e_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Specific internal energy from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: specific internal energy

Implements VaporMixtureFluidProperties.

Definition at line 191 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = v_from_p_T(p, T, x);
   return e_from_T_v(T, v, x);
 }

◆ e_from_p_T() [2/2]

void IdealRealGasMixtureFluidProperties::e_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT,
		std::vector< Real > &	de_dx
	)		const

overridevirtual

Specific internal energy and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values
[out]	e	specific internal energy
[out]	de_dp	derivative of specific internal energy w.r.t. pressure
[out]	de_dT	derivative of specific internal energy w.r.t. temperature
[out]	de_dx	derivative of specific internal energy w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 198 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v, de_dT_v, de_dv_T, p_unused, dp_dT_v, dp_dv_T;
   std::vector<Real> de_dx_Tv, dp_dx_Tv;
   de_dx_Tv.resize(_n_secondary_vapors);
   dp_dx_Tv.resize(_n_secondary_vapors);
 
   v = v_from_p_T(p, T, x);
   e_from_T_v(T, v, x, e, de_dT_v, de_dv_T, de_dx_Tv);
   p_from_T_v(T, v, x, p_unused, dp_dT_v, dp_dv_T, dp_dx_Tv);
 
   de_dp = de_dv_T / dp_dv_T;
   de_dT = de_dT_v - de_dv_T * dp_dT_v / dp_dv_T;
 
   de_dx.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     de_dx[i] = -(de_dv_T - de_dx_Tv[i] * dp_dv_T / dp_dx_Tv[i]) * dp_dx_Tv[i] / dp_dv_T;
 }

◆ e_from_T_v() [1/2]

Real IdealRealGasMixtureFluidProperties::e_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: specific internal energy

Definition at line 530 of file IdealRealGasMixtureFluidProperties.C.

Referenced by e_from_p_T(), and p_T_from_v_e().

 {
   const Real x_primary = primaryMassFraction(x);
   Real e = x_primary * _fp_primary->e_from_T_v(T, v / x_primary);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     e += x[i] * _fp_secondary[i]->e_from_T_v(T, v / x[i]);
 
   return e;
 }

◆ e_from_T_v() [2/2]

void IdealRealGasMixtureFluidProperties::e_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x,
		Real &	e,
		Real &	de_dT,
		Real &	de_dv,
		std::vector< Real > &	de_dx
	)		const

Specific internal energy and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values
[out]	e	specific internal energy
[out]	de_dT	derivative of specific internal energy w.r.t. temperature
[out]	de_dv	derivative of specific internal energy w.r.t. specific volume
[out]	de_dx	derivative of specific internal energy w.r.t. vapor mass fraction values

Definition at line 542 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real e_primary, de_dT_primary, de_dv_primary, de_dx_primary, dxi_dx_primary;
   Real de_dT_sec, dxj_dxi, dx_primary_dxi;
   std::vector<Real> e_sec, de_dv_sec;
 
   const Real x_primary = primaryMassFraction(x);
   de_dx.resize(_n_secondary_vapors);
   e_sec.resize(_n_secondary_vapors);
   de_dv_sec.resize(_n_secondary_vapors);
 
   _fp_primary->e_from_T_v(T, v / x_primary, e_primary, de_dT_primary, de_dv_primary);
   e = x_primary * e_primary;
   de_dT = x_primary * de_dT_primary;
   de_dv = de_dv_primary;
   de_dx_primary = e_primary - x_primary * de_dv_primary * v / (x_primary * x_primary);
 
   // get the partial pressures and their derivatives first
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->e_from_T_v(T, v / x[i], e_sec[i], de_dT_sec, de_dv_sec[i]);
 
     e += x[i] * e_sec[i];
     de_dT += x[i] * de_dT_sec;
     de_dv += de_dv_sec[i];
     dxi_dx_primary = -x[i] / (1. - x_primary);
     de_dx_primary +=
         dxi_dx_primary * e_sec[i] - x[i] * de_dv_sec[i] * v / (x[i] * x[i]) * dxi_dx_primary;
   }
 
   // get the composition dependent derivatives of the secondary vapors
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     de_dx[i] = e_sec[i] - x[i] * de_dv_sec[i] * v / (x[i] * x[i]);
     for (unsigned int j = 0; j < _n_secondary_vapors; j++)
     {
       if (j == i)
         continue;
       dxj_dxi = -x[j] / (1. - x[i]);
       de_dx[i] += dxj_dxi * e_sec[j] - x[j] * de_dv_sec[j] * v / (x[j] * x[j]) * dxj_dxi;
     }
     dx_primary_dxi = -x_primary / (1. - x[i]);
     de_dx[i] += dx_primary_dxi * e_primary -
                 x_primary * de_dv_primary * v / (x_primary * x_primary) * dx_primary_dxi;
   }
 }

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Definition at line 27 of file FluidProperties.h.

27 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 29 of file FluidProperties.h.

29 {}

◆ getNaN() [1/2]

Real NaNInterface::getNaN ( ) const

inlineprotectedinherited

Produces errors, warnings, or just quiet NaNs.

Definition at line 49 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), NaNInterfaceTestFluidProperties::p_from_v_e(), v_from_p_T(), and xs_prim_from_p_T().

49 { return getNaN("A NaN was produced."); }

NaNInterface::getNaN

Real getNaN() const

Produces errors, warnings, or just quiet NaNs.

Definition: NaNInterface.h:49

◆ getNaN() [2/2]

template<typename... Args>

Real NaNInterface::getNaN ( Args &&... args ) const

inlineprotectedinherited

Definition at line 52 of file NaNInterface.h.

   {
     switch (_emit_on_nan)
     {
       case (NAN_MESSAGE_WARNING):
         mooseWarning(_moose_object->name(), ": ", std::forward<Args>(args)...);
         break;
       case (NAN_MESSAGE_ERROR):
         mooseError(_moose_object->name(), ": ", std::forward<Args>(args)...);
         break;
       default:
         break;
     }
     // return a quiet NaN
     return std::nan("");
   }

◆ getNumberOfSecondaryVapors()

virtual unsigned int IdealRealGasMixtureFluidProperties::getNumberOfSecondaryVapors ( ) const

inlineoverridevirtual

Number of secondary vapors (non-condensable components)

Returns: number of secondary vapors

Implements VaporMixtureFluidProperties.

Definition at line 37 of file IdealRealGasMixtureFluidProperties.h.

37 { return _n_secondary_vapors; }

IdealRealGasMixtureFluidProperties::_n_secondary_vapors

const unsigned int _n_secondary_vapors

Number of secondary vapors.

Definition: IdealRealGasMixtureFluidProperties.h:462

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Definition at line 28 of file FluidProperties.h.

28 {}

◆ k_from_p_T()

Real IdealRealGasMixtureFluidProperties::k_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Thermal conductivity from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: thermal conductivity

Implements VaporMixtureFluidProperties.

Definition at line 348 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = v_from_p_T(p, T, x);
 
   const Real x_primary = primaryMassFraction(x);
   Real M_primary = _fp_primary->molarMass();
 
   Real sum = x_primary / M_primary;
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     sum += x[i] / _fp_secondary[i]->molarMass();
   Real M_star = 1. / sum;
 
   Real vp = v / x_primary;
   Real ep = _fp_primary->e_from_T_v(T, vp);
   Real k = x_primary * M_star / M_primary * _fp_primary->k_from_v_e(vp, ep);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real vi = v / x[i];
     Real ei = _fp_secondary[i]->e_from_T_v(T, vp);
     Real Mi = _fp_secondary[i]->molarMass();
     k += x[i] * M_star / Mi * _fp_secondary[i]->k_from_v_e(vi, ei);
   }
 
   return k;
 }

◆ k_from_T_v()

Real IdealRealGasMixtureFluidProperties::k_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Thermal conductivity from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: thermal conductivity

Definition at line 737 of file IdealRealGasMixtureFluidProperties.C.

 {
   const Real x_primary = primaryMassFraction(x);
   Real M_primary = _fp_primary->molarMass();
 
   Real sum = x_primary / M_primary;
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     sum += x[i] / _fp_secondary[i]->molarMass();
   Real M_star = 1. / sum;
 
   Real vp = v / x_primary;
   Real ep = _fp_primary->e_from_T_v(T, vp);
   Real k = x_primary * M_star / M_primary * _fp_primary->k_from_v_e(vp, ep);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real vi = v / x[i];
     Real ei = _fp_secondary[i]->e_from_T_v(T, vp);
     Real Mi = _fp_secondary[i]->molarMass();
     k += x[i] * M_star / Mi * _fp_secondary[i]->k_from_v_e(vi, ei);
   }
 
   return k;
 }

◆ mu_from_p_T()

Real IdealRealGasMixtureFluidProperties::mu_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Dynamic viscosity from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: dynamic viscosity

Implements VaporMixtureFluidProperties.

Definition at line 320 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = v_from_p_T(p, T, x);
 
   const Real x_primary = primaryMassFraction(x);
   Real M_primary = _fp_primary->molarMass();
 
   Real sum = x_primary / M_primary;
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     sum += x[i] / _fp_secondary[i]->molarMass();
   Real M_star = 1. / sum;
 
   Real vp = v / x_primary;
   Real ep = _fp_primary->e_from_T_v(T, vp);
   Real mu = x_primary * M_star / M_primary * _fp_primary->mu_from_v_e(vp, ep);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real vi = v / x[i];
     Real ei = _fp_secondary[i]->e_from_T_v(T, vp);
     Real Mi = _fp_secondary[i]->molarMass();
     mu += x[i] * M_star / Mi * _fp_secondary[i]->mu_from_v_e(vi, ei);
   }
 
   return mu;
 }

◆ mu_from_T_v()

Real IdealRealGasMixtureFluidProperties::mu_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Dynamic viscosity from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: dynamic viscosity

Definition at line 711 of file IdealRealGasMixtureFluidProperties.C.

 {
   const Real x_primary = primaryMassFraction(x);
   Real M_primary = _fp_primary->molarMass();
 
   Real sum = x_primary / M_primary;
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     sum += x[i] / _fp_secondary[i]->molarMass();
   Real M_star = 1. / sum;
 
   Real vp = v / x_primary;
   Real ep = _fp_primary->e_from_T_v(T, vp);
   Real mu = x_primary * M_star / M_primary * _fp_primary->mu_from_v_e(vp, ep);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     Real vi = v / x[i];
     Real ei = _fp_secondary[i]->e_from_T_v(T, vp);
     Real Mi = _fp_secondary[i]->molarMass();
     mu += x[i] * M_star / Mi * _fp_secondary[i]->mu_from_v_e(vi, ei);
   }
 
   return mu;
 }

◆ p_from_T_v() [1/3]

Real IdealRealGasMixtureFluidProperties::p_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x
	)		const

Pressure from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: pressure

Definition at line 444 of file IdealRealGasMixtureFluidProperties.C.

Referenced by c_from_p_T(), c_from_T_v(), cp_from_p_T(), cp_from_T_v(), e_from_p_T(), p_T_from_v_e(), and v_from_p_T().

 {
   const Real x_primary = primaryMassFraction(x);
   Real p = _fp_primary->p_from_T_v(T, v / x_primary);
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     p += _fp_secondary[i]->p_from_T_v(T, v / x[i]);
 
   return p;
 }

◆ p_from_T_v() [2/3]

void IdealRealGasMixtureFluidProperties::p_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x,
		Real &	p,
		Real &	dp_dT,
		Real &	dp_dv
	)		const

Pressure and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values
[out]	p	pressure
[out]	dp_dT	derivative of pressure w.r.t. temperature
[out]	dp_dv	derivative of pressure w.r.t. specific volume

Definition at line 456 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p_primary, dp_dT_primary, dp_dv_primary;
   Real p_sec, dp_dT_sec, dp_dv_sec;
 
   const Real x_primary = primaryMassFraction(x);
 
   _fp_primary->p_from_T_v(T, v / x_primary, p_primary, dp_dT_primary, dp_dv_primary);
   p = p_primary;
   dp_dT = dp_dT_primary;
   dp_dv = dp_dv_primary / x_primary;
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->p_from_T_v(T, v / x[i], p_sec, dp_dT_sec, dp_dv_sec);
 
     p += p_sec;
     dp_dT += dp_dT_sec;
     dp_dv += dp_dv_sec / x[i];
   }
 }

◆ p_from_T_v() [3/3]

void IdealRealGasMixtureFluidProperties::p_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x,
		Real &	p,
		Real &	dp_dT,
		Real &	dp_dv,
		std::vector< Real > &	dp_dx
	)		const

Pressure and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values
[out]	p	pressure
[out]	dp_dT	derivative of pressure w.r.t. temperature
[out]	dp_dv	derivative of pressure w.r.t. specific volume
[out]	dp_dx	derivative of pressure w.r.t. vapor mass fraction values

Definition at line 480 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p_primary, dp_dT_primary, dp_dv_primary, dp_dx_primary, dxi_dx_primary;
   Real p_sec, dp_dT_sec, dxj_dxi;
   std::vector<Real> dp_dv_sec;
 
   const Real x_primary = primaryMassFraction(x);
   dp_dx.resize(_n_secondary_vapors);
   dp_dv_sec.resize(_n_secondary_vapors);
 
   _fp_primary->p_from_T_v(T, v / x_primary, p_primary, dp_dT_primary, dp_dv_primary);
   p = p_primary;
   dp_dT = dp_dT_primary;
   dp_dv = dp_dv_primary / x_primary;
   dp_dx_primary = -dp_dv_primary * v / (x_primary * x_primary);
 
   // get the partial pressures and their derivatives first
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->p_from_T_v(T, v / x[i], p_sec, dp_dT_sec, dp_dv_sec[i]);
 
     p += p_sec;
     dp_dT += dp_dT_sec;
     dp_dv += dp_dv_sec[i] / x[i];
     dxi_dx_primary = -x[i] / (1. - x_primary);
     dp_dx_primary += -dp_dv_sec[i] * v / (x[i] * x[i]) * dxi_dx_primary;
   }
 
   // get the composition dependent derivatives of the secondary vapors
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     dp_dx[i] = -dp_dv_sec[i] * v / (x[i] * x[i]);
     for (unsigned int j = 0; j < _n_secondary_vapors; j++)
     {
       if (j == i)
         continue;
       dxj_dxi = -x[j] / (1. - x[i]);
       dp_dx[i] += -dp_dv_sec[j] * v / (x[j] * x[j]) * dxj_dxi;
     }
     dp_dx[i] += -dp_dv_primary * v / (x_primary * x_primary) * (-x_primary / (1. - x[i]));
   }
 }

◆ p_from_v_e() [1/2]

Real IdealRealGasMixtureFluidProperties::p_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x
	)		const

overridevirtual

Pressure from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values

Returns: pressure

Implements VaporMixtureFluidProperties.

Definition at line 55 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p, T;
   p_T_from_v_e(v, e, x, p, T);
 
   return p;
 }

◆ p_from_v_e() [2/2]

void IdealRealGasMixtureFluidProperties::p_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x,
		Real &	p,
		Real &	dp_dv,
		Real &	dp_de,
		std::vector< Real > &	dp_dx
	)		const

overridevirtual

Pressure and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values
[out]	p	pressure
[out]	dp_dv	derivative of pressure w.r.t. specific volume
[out]	dp_de	derivative of pressure w.r.t. specific internal energy
[out]	dp_dx	derivative of pressure w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 64 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real T, dT_dv, dT_de;
   std::vector<Real> dT_dx(_n_secondary_vapors);
   p_T_from_v_e(v, e, x, p, dp_dv, dp_de, dp_dx, T, dT_dv, dT_de, dT_dx);
 }

◆ p_T_from_v_e() [1/2]

void IdealRealGasMixtureFluidProperties::p_T_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x,
		Real &	p,
		Real &	T
	)		const

Pressure and temperature from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values
[out]	p	pressure
[out]	T	temperature

Definition at line 376 of file IdealRealGasMixtureFluidProperties.C.

Referenced by p_from_v_e(), p_T_from_v_e(), and T_from_v_e().

 {
   Real v_primary = v / primaryMassFraction(x);
   static const Real vc = 1. / _fp_primary->criticalDensity();
   static const Real ec = _fp_primary->criticalInternalEnergy();
 
   // Initial estimate of a bracketing interval for the temperature
   Real lower_temperature, upper_temperature;
   if (v_primary > vc)
   {
     Real e_sat_primary = _fp_primary->e_spndl_from_v(v_primary);
     lower_temperature = _fp_primary->T_from_v_e(v_primary, e_sat_primary);
   }
   else
     lower_temperature = _fp_primary->T_from_v_e(v_primary, ec);
 
   upper_temperature = _T_mix_max;
 
   // Use BrentsMethod to find temperature
   auto energy_diff = [&v, &e, &x, this](Real T) { return this->e_from_T_v(T, v, x) - e; };
 
   BrentsMethod::bracket(energy_diff, lower_temperature, upper_temperature);
   T = BrentsMethod::root(energy_diff, lower_temperature, upper_temperature);
 
   p = p_from_T_v(T, v, x);
 }

◆ p_T_from_v_e() [2/2]

void IdealRealGasMixtureFluidProperties::p_T_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x,
		Real &	p,
		Real &	dp_dv,
		Real &	dp_de,
		std::vector< Real > &	dp_dx,
		Real &	T,
		Real &	dT_dv,
		Real &	dT_de,
		std::vector< Real > &	dT_dx
	)		const

Pressure and temperature from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values
[out]	p	pressure
[out]	dp_dv	derivative of pressure w.r.t. specific volume
[out]	dp_de	derivative of pressure w.r.t. specific internal energy
[out]	dp_dx	derivative of pressure w.r.t. vapor mass fraction values
[out]	T	temperature
[out]	dT_dv	derivative of temperature w.r.t. specific volume
[out]	dT_de	derivative of temperature w.r.t. specific internal energy
[out]	dT_dx	derivative of temperature w.r.t. vapor mass fraction values

Definition at line 405 of file IdealRealGasMixtureFluidProperties.C.

 {
   p_T_from_v_e(v, e, x, p, T);
 
   // pressure and derivatives
   Real p_unused, dp_dT_v, dp_dv_T;
   std::vector<Real> dp_dx_Tv;
   p_from_T_v(T, v, x, p_unused, dp_dT_v, dp_dv_T, dp_dx_Tv);
 
   // internal energy and derivatives
   Real e_unused, de_dT_v, de_dv_T;
   std::vector<Real> de_dx_Tv;
   e_from_T_v(T, v, x, e_unused, de_dT_v, de_dv_T, de_dx_Tv);
 
   // Compute derivatives using the following rules:
   dp_dv = dp_dv_T - dp_dT_v * de_dv_T / de_dT_v;
   dp_de = dp_dT_v / de_dT_v;
   dT_dv = -de_dv_T / de_dT_v;
   dT_de = 1. / de_dT_v;
 
   // Derivatives with respect to mass fractions:
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     dT_dx[i] = -de_dx_Tv[i] / de_dT_v;
     dp_dx[i] = dp_dx_Tv[i] + dp_dT_v * dT_dx[i];
   }
 }

◆ primaryMassFraction()

Real VaporMixtureFluidProperties::primaryMassFraction ( const std::vector< Real > & x ) const

inherited

Computes the mass fraction of the primary vapor given mass fractions of the secondary vapors.

This uses the relation

\[ \sum\limits_i^N x_i = 1 , \]

where the mass fractions \(x_i, i=2\ldots N\) correspond to the secondary vapors.

Definition at line 29 of file VaporMixtureFluidProperties.C.

Referenced by GeneralVaporMixtureFluidProperties::c_from_p_T(), GeneralVaporMixtureFluidProperties::cp_from_p_T(), cp_from_p_T(), cp_from_T_v(), GeneralVaporMixtureFluidProperties::cv_from_p_T(), cv_from_p_T(), cv_from_T_v(), GeneralVaporMixtureFluidProperties::e_from_p_T(), e_from_T_v(), GeneralVaporMixtureFluidProperties::k_from_p_T(), k_from_p_T(), k_from_T_v(), GeneralVaporMixtureFluidProperties::mu_from_p_T(), mu_from_p_T(), mu_from_T_v(), p_from_T_v(), p_T_from_v_e(), s_from_T_v(), GeneralVaporMixtureFluidProperties::v_from_p_T(), v_from_p_T(), and xs_prim_from_p_T().

 {
   return 1 - std::accumulate(x.begin(), x.end(), 0.0);
 }

◆ rho_from_p_T() [1/2]

Real IdealRealGasMixtureFluidProperties::rho_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

overridevirtual

Density from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: density

Implements VaporMixtureFluidProperties.

Definition at line 101 of file IdealRealGasMixtureFluidProperties.C.

 {
   return 1.0 / v_from_p_T(p, T, x);
 }

◆ rho_from_p_T() [2/2]

void IdealRealGasMixtureFluidProperties::rho_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		std::vector< Real > &	drho_dx
	)		const

overridevirtual

Density and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values
[out]	rho	density
[out]	drho_dp	derivative of density w.r.t. pressure
[out]	drho_dT	derivative of density w.r.t. temperature
[out]	drho_dx	derivative of density w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 107 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v, dv_dp, dv_dT;
   std::vector<Real> dv_dx;
   v_from_p_T(p, T, x, v, dv_dp, dv_dT, dv_dx);
 
   rho = 1.0 / v;
   const Real drho_dv = -1.0 / (v * v);
   drho_dp = drho_dv * dv_dp;
   drho_dT = drho_dv * dv_dT;
   drho_dx.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     drho_dx[i] = drho_dv * dv_dx[i];
 }

◆ s_from_T_v()

void IdealRealGasMixtureFluidProperties::s_from_T_v	(	Real	T,
		Real	v,
		const std::vector< Real > &	x,
		Real &	s,
		Real &	ds_dT,
		Real &	ds_dv
	)		const

Specific entropy and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[in]	x	vapor mass fraction values
[out]	s	specific entropy
[out]	ds_dT	derivative of specific entropy w.r.t. temperature
[out]	ds_dv	derivative of specific entropy w.r.t. specific volume
[out]	ds_dx	derivative of specific entropy w.r.t. vapor mass fraction values

Definition at line 596 of file IdealRealGasMixtureFluidProperties.C.

Referenced by c_from_p_T(), and c_from_T_v().

 {
   Real s_primary, ds_dT_primary, ds_dv_primary;
   Real s_sec, ds_dT_sec, ds_dv_sec;
 
   const Real x_primary = primaryMassFraction(x);
 
   _fp_primary->s_from_T_v(T, v / x_primary, s_primary, ds_dT_primary, ds_dv_primary);
   s = x_primary * s_primary;
   ds_dT = x_primary * ds_dT_primary;
   ds_dv = ds_dv_primary;
 
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     _fp_secondary[i]->s_from_T_v(T, v / x[i], s_sec, ds_dT_sec, ds_dv_sec);
 
     s += x[i] * s_sec;
     ds_dT += x[i] * ds_dT_sec;
     ds_dv += ds_dv_sec;
   }
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 32 of file FluidProperties.h.

32 {}

◆ T_from_v_e() [1/2]

Real IdealRealGasMixtureFluidProperties::T_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x
	)		const

overridevirtual

Temperature from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values

Returns: temperature

Implements VaporMixtureFluidProperties.

Definition at line 78 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p, T;
   p_T_from_v_e(v, e, x, p, T);
 
   return T;
 }

◆ T_from_v_e() [2/2]

void IdealRealGasMixtureFluidProperties::T_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x,
		Real &	T,
		Real &	dT_dv,
		Real &	dT_de,
		std::vector< Real > &	dT_dx
	)		const

overridevirtual

Temperature and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[in]	x	vapor mass fraction values
[out]	T	temperature
[out]	dT_dv	derivative of temperature w.r.t. specific volume
[out]	dT_de	derivative of temperature w.r.t. specific internal energy
[out]	dT_dx	derivative of temperature w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 87 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p, dp_dv, dp_de;
   std::vector<Real> dp_dx(_n_secondary_vapors);
   p_T_from_v_e(v, e, x, p, dp_dv, dp_de, dp_dx, T, dT_dv, dT_de, dT_dx);
 }

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Definition at line 31 of file FluidProperties.h.

31 {}

◆ v_from_p_T() [1/2]

Real IdealRealGasMixtureFluidProperties::v_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

virtual

Specific volume from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values

Returns: specific volume

Definition at line 129 of file IdealRealGasMixtureFluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), k_from_p_T(), mu_from_p_T(), rho_from_p_T(), and v_from_p_T().

 {
   const Real x_primary = primaryMassFraction(x);
   Real M_primary = _fp_primary->molarMass();
 
   Real sum = x_primary / M_primary;
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     sum += x[i] / _fp_secondary[i]->molarMass();
   Real M_star = 1. / sum;
   Real v_ideal = R_molar * T / (M_star * p);
 
   // check range of validity for primary (condensable) component
   static const Real Tc = _fp_primary->criticalTemperature();
   static const Real vc = 1. / _fp_primary->criticalDensity();
   Real v_spndl, e_spndl;
   if (T < Tc)
     _fp_primary->v_e_spndl_from_T(T, v_spndl, e_spndl);
   else
     v_spndl = vc;
 
   Real lower_spec_volume = v_spndl * x_primary;
   Real upper_spec_volume = v_ideal; // p*v/(RT) <= 1
 
   // Initial estimate of a bracketing interval for the temperature
   Real p_max = p_from_T_v(T, lower_spec_volume, x);
   if (p > p_max || upper_spec_volume < lower_spec_volume)
     return getNaN();
 
   // Use BrentsMethod to find temperature
   auto pressure_diff = [&T, &p, &x, this](Real v) { return this->p_from_T_v(T, v, x) - p; };
 
   BrentsMethod::bracket(pressure_diff, lower_spec_volume, upper_spec_volume);
   Real v = BrentsMethod::root(pressure_diff, lower_spec_volume, upper_spec_volume);
 
   return v;
 }

◆ v_from_p_T() [2/2]

void IdealRealGasMixtureFluidProperties::v_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x,
		Real &	v,
		Real &	dv_dp,
		Real &	dv_dT,
		std::vector< Real > &	dv_dx
	)		const

virtual

Specific volume and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[in]	x	vapor mass fraction values
[out]	v	specific volume
[out]	dv_dp	derivative of specific volume w.r.t. pressure
[out]	dv_dT	derivative of specific volume w.r.t. temperature
[out]	dv_dx	derivative of specific volume w.r.t. vapor mass fraction values

Definition at line 167 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p_unused, dp_dT, dp_dv;
   std::vector<Real> dp_dx;
   dp_dx.resize(_n_secondary_vapors);
 
   v = v_from_p_T(p, T, x);
   p_from_T_v(T, v, x, p_unused, dp_dT, dp_dv, dp_dx);
 
   dv_dp = 1. / dp_dv;
   dv_dT = -dp_dT / dp_dv;
 
   dv_dx.resize(_n_secondary_vapors);
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
     dv_dx[i] = -dp_dx[i] / dp_dv;
 }

◆ xs_prim_from_p_T()

Real IdealRealGasMixtureFluidProperties::xs_prim_from_p_T	(	Real	p,
		Real	T,
		const std::vector< Real > &	x
	)		const

Mass fraction of primary (condensable) component at saturation from pressure and temperature.

Parameters

[in]	T	temperature
[in]	p	pressure

Returns: mass fraction of primary (condensable) component at saturation

Definition at line 763 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real T_c = _fp_primary->criticalTemperature();
   Real xs;
   if (T > T_c)
     // return 1. to indicate that no water would condense for
     // given (p,T)
     xs = 1.;
   else
   {
     Real pp_sat = _fp_primary->pp_sat_from_p_T(p, T);
     if (pp_sat < 0.)
     {
       // return 1. to indicate that no water would condense for
       // given (p,T)
       xs = 1.;
       return xs;
     }
     Real v_primary = _fp_primary->v_from_p_T(pp_sat, T);
     Real pp_sat_secondary = p - pp_sat;
 
     Real v_secondary;
     if (_n_secondary_vapors == 1)
       v_secondary = _fp_secondary[0]->v_from_p_T(pp_sat_secondary, T);
     else
     {
       std::vector<Real> x_sec(_n_secondary_vapors);
       const Real x_primary = primaryMassFraction(x);
       Real sum = 0.;
       for (unsigned int i = 0; i < _n_secondary_vapors; i++)
       {
         x_sec[i] = x[i] / (1. - x_primary);
         sum += x_sec[i] / _fp_secondary[i]->molarMass();
       }
       Real M_star = 1. / sum;
       v_secondary = R_molar * T / (M_star * pp_sat_secondary);
       int it = 0;
       double f = 1., df_dvs, pp_sec, p_sec, dp_dT_sec, dp_dv_sec, dp_dT, dp_dv;
       double tol_p = 1.e-8;
       while (std::fabs(f / pp_sat_secondary) > tol_p)
       {
         pp_sec = 0.;
         dp_dT = 0.;
         dp_dv = 0.;
         for (unsigned int i = 0; i < _n_secondary_vapors; i++)
         {
           _fp_secondary[i]->p_from_T_v(T, v_secondary / x_sec[i], p_sec, dp_dT_sec, dp_dv_sec);
           pp_sec += p_sec;
           dp_dT += dp_dT_sec;
           dp_dv += dp_dv_sec / x_sec[i];
         }
         f = pp_sec - pp_sat_secondary;
         df_dvs = dp_dv;
         v_secondary -= f / df_dvs;
         if (it++ > 15)
           return getNaN();
       }
     }
 
     xs = v_secondary / (v_primary + v_secondary);
   }
 
   return xs;
 }

Member Data Documentation

◆ _emit_on_nan

enum NaNMessage NaNInterface::_emit_on_nan

protectedinherited

Raise mooseWarning or mooseError?

Definition at line 44 of file NaNInterface.h.

Referenced by NaNInterface::getNaN().

◆ _fp_primary

const SinglePhaseFluidProperties* IdealRealGasMixtureFluidProperties::_fp_primary

protected

Primary vapor fluid properties.

Definition at line 456 of file IdealRealGasMixtureFluidProperties.h.

Referenced by cp_from_p_T(), cp_from_T_v(), cv_from_p_T(), cv_from_T_v(), e_from_T_v(), IdealRealGasMixtureFluidProperties(), k_from_p_T(), k_from_T_v(), mu_from_p_T(), mu_from_T_v(), p_from_T_v(), p_T_from_v_e(), s_from_T_v(), v_from_p_T(), and xs_prim_from_p_T().

◆ _fp_secondary

std::vector<const SinglePhaseFluidProperties *> IdealRealGasMixtureFluidProperties::_fp_secondary

protected

Secondary vapor fluid properties.

Definition at line 458 of file IdealRealGasMixtureFluidProperties.h.

Referenced by cp_from_p_T(), cp_from_T_v(), cv_from_p_T(), cv_from_T_v(), e_from_T_v(), IdealRealGasMixtureFluidProperties(), k_from_p_T(), k_from_T_v(), mu_from_p_T(), mu_from_T_v(), p_from_T_v(), s_from_T_v(), v_from_p_T(), and xs_prim_from_p_T().

◆ _fp_secondary_names

const std::vector<UserObjectName> IdealRealGasMixtureFluidProperties::_fp_secondary_names

protected

Names of secondary vapor fluid properties.

Definition at line 460 of file IdealRealGasMixtureFluidProperties.h.

Referenced by IdealRealGasMixtureFluidProperties().

◆ _moose_object

const MooseObject* NaNInterface::_moose_object

protectedinherited

Definition at line 41 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), and NaNInterface::NaNInterface().

◆ _n_secondary_vapors

const unsigned int IdealRealGasMixtureFluidProperties::_n_secondary_vapors

protected

Number of secondary vapors.

Definition at line 462 of file IdealRealGasMixtureFluidProperties.h.

Referenced by c_from_p_T(), c_from_T_v(), cp_from_p_T(), cp_from_T_v(), cv_from_p_T(), cv_from_T_v(), e_from_p_T(), e_from_T_v(), getNumberOfSecondaryVapors(), IdealRealGasMixtureFluidProperties(), k_from_p_T(), k_from_T_v(), mu_from_p_T(), mu_from_T_v(), p_from_T_v(), p_from_v_e(), p_T_from_v_e(), rho_from_p_T(), s_from_T_v(), T_from_v_e(), v_from_p_T(), and xs_prim_from_p_T().

◆ _T_mix_max

const Real IdealRealGasMixtureFluidProperties::_T_mix_max

protected

maximum temperature of all components

Definition at line 466 of file IdealRealGasMixtureFluidProperties.h.

Referenced by p_T_from_v_e().

◆ R_molar

constexpr const Real IdealRealGasMixtureFluidProperties::R_molar = 8.3144598

staticprotected

molar (or universal) gas constant

Definition at line 464 of file IdealRealGasMixtureFluidProperties.h.

Referenced by v_from_p_T(), and xs_prim_from_p_T().

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

fluid_properties/include/userobjects/IdealRealGasMixtureFluidProperties.h
fluid_properties/src/userobjects/IdealRealGasMixtureFluidProperties.C

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ NaNMessage

Constructor & Destructor Documentation

◆ IdealRealGasMixtureFluidProperties()

Member Function Documentation

◆ c_from_p_T() [1/2]

◆ c_from_p_T() [2/2]

◆ c_from_T_v() [1/2]

◆ c_from_T_v() [2/2]

◆ cp_from_p_T()

◆ cp_from_T_v()

◆ cv_from_p_T()

◆ cv_from_T_v()

◆ e_from_p_T() [1/2]

◆ e_from_p_T() [2/2]

◆ e_from_T_v() [1/2]

◆ e_from_T_v() [2/2]

◆ execute()

◆ finalize()

◆ getNaN() [1/2]

◆ getNaN() [2/2]

◆ getNumberOfSecondaryVapors()

◆ initialize()

◆ k_from_p_T()

◆ k_from_T_v()

◆ mu_from_p_T()

◆ mu_from_T_v()

◆ p_from_T_v() [1/3]

◆ p_from_T_v() [2/3]

◆ p_from_T_v() [3/3]

◆ p_from_v_e() [1/2]

◆ p_from_v_e() [2/2]

◆ p_T_from_v_e() [1/2]

◆ p_T_from_v_e() [2/2]

◆ primaryMassFraction()

◆ rho_from_p_T() [1/2]

◆ rho_from_p_T() [2/2]

◆ s_from_T_v()

◆ subdomainSetup()

◆ T_from_v_e() [1/2]

◆ T_from_v_e() [2/2]

◆ threadJoin()

◆ v_from_p_T() [1/2]

◆ v_from_p_T() [2/2]

◆ xs_prim_from_p_T()

Member Data Documentation

◆ _emit_on_nan

◆ _fp_primary

◆ _fp_secondary

◆ _fp_secondary_names

◆ _moose_object

◆ _n_secondary_vapors

◆ _T_mix_max

◆ R_molar