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	c_from_v_e (Real v, Real e, const std::vector< Real > &x) const override
	Speed of sound from specific volume and specific internal energy. More...

virtual void	c_from_v_e (Real v, Real e, const std::vector< Real > &x, Real &c, Real &dc_dv, Real &dc_de, std::vector< Real > &dc_dx) const override
	Speed of sound 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...

virtual Real	e_from_p_rho (Real p, Real rho, const std::vector< Real > &x) const override
	Specific internal energy from pressure and density. More...

virtual void	e_from_p_rho (Real p, Real rho, const std::vector< Real > &x, Real &e, Real &de_dp, Real &de_drho, std::vector< Real > &de_dx) const override
	Specific internal energy and its derivatives from pressure and density. 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	T_from_p_v (Real p, Real v, const std::vector< Real > &x) const
	Temperature from pressure and specific volume. More...

void	T_from_p_v (Real p, Real v, const std::vector< Real > &x, Real &T, Real &dT_dp, Real &dT_dv, std::vector< Real > &dT_dx) const
	Temperature and its derivatives from pressure and specific volume. 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

Static Public Attributes
static const Real	_R = 8.3144598
	Universal gas constant (J/mol/K) More...

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 Real	_T_c2k
	Conversion of temperature from Celsius to Kelvin. More...

const bool	_allow_imperfect_jacobians
	Flag to set unimplemented Jacobian entries to zero. More...

const MooseObject *	_moose_object

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

Static Protected Attributes
constexpr static 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 26 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 33 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 265 of file IdealRealGasMixtureFluidProperties.C.

 {
   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);
 }

Referenced by c_from_p_T().

◆ 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 283 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
   Real dp = p * 1.e-6;
   p_perturbed = p + dp;
   c_perturbed = c_from_p_T(p_perturbed, T, x);
   dc_dp = (c_perturbed - c) / (dp);
  
   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 756 of file IdealRealGasMixtureFluidProperties.C.

 {
   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);
 }

Referenced by c_from_T_v(), and c_from_v_e().

◆ 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 772 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);
  
   Real dv = v * 1.e-6;
   v_perturbed = v + dv;
   c_perturbed = c_from_T_v(T, v_perturbed, x);
   dc_dv = (c_perturbed - c) / (dv);
  
   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);
   }
 }

◆ c_from_v_e() [1/2]

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

overridevirtual

Speed of sound from specific volume and specific internal energy.

Parameters

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

Returns: speed of sound

Implements VaporMixtureFluidProperties.

Definition at line 101 of file IdealRealGasMixtureFluidProperties.C.

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

◆ c_from_v_e() [2/2]

void IdealRealGasMixtureFluidProperties::c_from_v_e	(	Real	v,
		Real	e,
		const std::vector< Real > &	x,
		Real &	c,
		Real &	dc_dv,
		Real &	dc_de,
		std::vector< Real > &	dc_dx
	)		const

overridevirtual

Speed of sound 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]	c	Speed of sound
[out]	dc_dv	derivative of temperature w.r.t. specific volume
[out]	dc_de	derivative of temperature w.r.t. specific internal energy
[out]	dc_dx	derivative of temperature w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 109 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real p, T;
   p_T_from_v_e(v, e, x, p, T);
  
   // sound of speed and derivatives
   Real dc_dT_v, dc_dv_T;
   std::vector<Real> dc_dx_Tv;
   c_from_T_v(T, v, x, c, dc_dT_v, dc_dv_T, dc_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:
   dc_dv = dc_dv_T - dc_dT_v * de_dv_T / de_dT_v;
   dc_de = dc_dT_v / de_dT_v;
  
   // Derivatives with respect to mass fractions:
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     dc_dx[i] = dc_dx_Tv[i] - dc_dT_v * de_dx_Tv[i] / de_dT_v;
   }
 }

◆ 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 325 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 814 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 348 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 836 of file IdealRealGasMixtureFluidProperties.C.

 {
   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;
 }

Referenced by cv_from_p_T().

◆ e_from_p_rho() [1/2]

Real IdealRealGasMixtureFluidProperties::e_from_p_rho	(	Real	p,
		Real	rho,
		const std::vector< Real > &	x
	)		const

overridevirtual

Specific internal energy from pressure and density.

Parameters

[in]	p	pressure
[in]	rho	density
[in]	x	vapor mass fraction values

Returns: specific internal energy

Implements VaporMixtureFluidProperties.

Definition at line 418 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = 1. / rho;
   Real T = T_from_p_v(p, v, x);
  
   return e_from_T_v(T, v, x);
 }

◆ e_from_p_rho() [2/2]

void IdealRealGasMixtureFluidProperties::e_from_p_rho	(	Real	p,
		Real	rho,
		const std::vector< Real > &	x,
		Real &	e,
		Real &	de_dp,
		Real &	de_drho,
		std::vector< Real > &	de_dx
	)		const

overridevirtual

Specific internal energy and its derivatives from pressure and density.

Parameters

[in]	p	pressure
[in]	rho	density
[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_drho	derivative of specific internal energy w.r.t. density
[out]	de_dx	derivative of specific internal energy w.r.t. vapor mass fraction values

Implements VaporMixtureFluidProperties.

Definition at line 429 of file IdealRealGasMixtureFluidProperties.C.

 {
   Real v = 1. / rho;
   Real T = T_from_p_v(p, v, x);
   Real p_, dp_dT, dp_dv, de_dT, de_dv;
   std::vector<Real> dp_dx_Tv;
   std::vector<Real> de_dx_Tv;
  
   p_from_T_v(T, v, x, p_, dp_dT, dp_dv, dp_dx_Tv);
   e_from_T_v(T, v, x, e, de_dT, de_dv, de_dx_Tv);
  
   de_dp = de_dT / dp_dT;
   de_drho = (-v * v) * (de_dv - de_dT * dp_dv / dp_dT);
  
   // Derivatives with respect to mass fractions:
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     de_dx[i] = de_dx_Tv[i] - de_dT * dp_dx_Tv[i] / dp_dT;
   }
 }

◆ 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 232 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 239 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 666 of file IdealRealGasMixtureFluidProperties.C.

 {
   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;
 }

Referenced by c_from_v_e(), e_from_p_rho(), e_from_p_T(), and p_T_from_v_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 678 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 34 of file FluidProperties.h.

34 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 36 of file FluidProperties.h.

36 {}

◆ getNaN() [1/2]

Real NaNInterface::getNaN ( ) const

inlineprotectedinherited

Produces errors, warnings, or just quiet NaNs.

Definition at line 48 of file NaNInterface.h.

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

Referenced by StiffenedGasFluidProperties::c_from_v_e(), NaNInterface::getNaN(), NaNInterfaceTestFluidProperties::p_from_v_e(), StiffenedGasFluidProperties::rho_from_p_s(), StiffenedGasFluidProperties::rho_from_p_T(), StiffenedGasFluidProperties::s_from_h_p(), StiffenedGasFluidProperties::s_from_p_T(), StiffenedGasFluidProperties::s_from_v_e(), v_from_p_T(), and xs_prim_from_p_T().

◆ getNaN() [2/2]

template<typename... Args>

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

inlineprotectedinherited

Definition at line 51 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 36 of file IdealRealGasMixtureFluidProperties.h.

36 { return _n_secondary_vapors; }

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Definition at line 35 of file FluidProperties.h.

35 {}

◆ 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 390 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 874 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 362 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 848 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 580 of file IdealRealGasMixtureFluidProperties.C.

 {
   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;
 }

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

◆ 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 592 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 616 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 &	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 486 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];
   }
 }

◆ 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 &	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 457 of file IdealRealGasMixtureFluidProperties.C.

 {
   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);
 }

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

◆ 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.

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

Referenced by cp_from_p_T(), cp_from_T_v(), cv_from_p_T(), cv_from_T_v(), e_from_T_v(), 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(), T_from_p_v(), v_from_p_T(), and xs_prim_from_p_T().

◆ 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 142 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 148 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 732 of file IdealRealGasMixtureFluidProperties.C.

 {
   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;
   }
 }

Referenced by c_from_p_T(), and c_from_T_v().

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 39 of file FluidProperties.h.

39 {}

◆ T_from_p_v() [1/2]

Real IdealRealGasMixtureFluidProperties::T_from_p_v	(	Real	p,
		Real	v,
		const std::vector< Real > &	x
	)		const

Temperature from pressure and specific volume.

Parameters

[in]	p	pressure
[in]	v	specific volume
[in]	x	vapor mass fraction values

Returns: temperature

Definition at line 525 of file IdealRealGasMixtureFluidProperties.C.

 {
   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 pressure_diff = [&p, &v, &x, this](Real T) { return this->p_from_T_v(T, v, x) - p; };
  
   BrentsMethod::bracket(pressure_diff, lower_temperature, upper_temperature);
   Real T = BrentsMethod::root(pressure_diff, lower_temperature, upper_temperature);
  
   return T;
 }

Referenced by e_from_p_rho(), and T_from_p_v().

◆ T_from_p_v() [2/2]

void IdealRealGasMixtureFluidProperties::T_from_p_v	(	Real	p,
		Real	v,
		const std::vector< Real > &	x,
		Real &	T,
		Real &	dT_dp,
		Real &	dT_dv,
		std::vector< Real > &	dT_dx
	)		const

Temperature and its derivatives from pressure and specific volume.

Parameters

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

Definition at line 553 of file IdealRealGasMixtureFluidProperties.C.

 {
   T = T_from_p_v(p, v, x);
  
   // 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);
  
   // Compute derivatives using the following rules:
   dT_dp = 1. / dp_dT_v;
   dT_dv = -dp_dv_T / dp_dT_v;
  
   // Derivatives with respect to mass fractions:
   for (unsigned int i = 0; i < _n_secondary_vapors; i++)
   {
     dT_dx[i] = -dp_dx_Tv[i] / dp_dT_v;
   }
 }

◆ 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 38 of file FluidProperties.h.

38 {}

◆ 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 170 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 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;
 }

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().

◆ 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 208 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 900 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

◆ _allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians

protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 48 of file FluidProperties.h.

Referenced by SinglePhaseFluidProperties::fluidPropError().

◆ _emit_on_nan

enum NaNMessage NaNInterface::_emit_on_nan

protectedinherited

Raise mooseWarning or mooseError?

Definition at line 43 of file NaNInterface.h.

Referenced by NaNInterface::getNaN().

◆ _fp_primary

const SinglePhaseFluidProperties* IdealRealGasMixtureFluidProperties::_fp_primary

protected

Primary vapor fluid properties.

Definition at line 542 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(), T_from_p_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 544 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 546 of file IdealRealGasMixtureFluidProperties.h.

Referenced by IdealRealGasMixtureFluidProperties().

◆ _moose_object

const MooseObject* NaNInterface::_moose_object

protectedinherited

Definition at line 40 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 548 of file IdealRealGasMixtureFluidProperties.h.

Referenced by c_from_p_T(), c_from_T_v(), c_from_v_e(), cp_from_p_T(), cp_from_T_v(), cv_from_p_T(), cv_from_T_v(), e_from_p_rho(), 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_p_v(), T_from_v_e(), v_from_p_T(), and xs_prim_from_p_T().

◆ _R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 42 of file FluidProperties.h.

Referenced by HelmholtzFluidProperties::c_from_p_T(), IdealGasFluidProperties::c_from_p_T(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), HelmholtzFluidProperties::p_from_rho_T(), HelmholtzFluidProperties::rho_from_p_T(), IdealGasFluidProperties::rho_from_p_T(), and HelmholtzFluidProperties::s_from_p_T().

◆ _T_c2k

const Real FluidProperties::_T_c2k

protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 46 of file FluidProperties.h.

Referenced by NaClFluidProperties::cp_from_p_T(), BrineFluidProperties::cp_from_p_T_X(), NaClFluidProperties::h_from_p_T(), BrineFluidProperties::h_from_p_T_X(), BrineFluidProperties::haliteSolubility(), NaClFluidProperties::k_from_p_T(), BrineFluidProperties::k_from_p_T_X(), BrineFluidProperties::mu_from_p_T_X(), CO2FluidProperties::partialDensity(), NaClFluidProperties::rho_from_p_T(), and BrineFluidProperties::rho_from_p_T_X().

◆ _T_mix_max

const Real IdealRealGasMixtureFluidProperties::_T_mix_max

protected

maximum temperature of all components

Definition at line 552 of file IdealRealGasMixtureFluidProperties.h.

Referenced by p_T_from_v_e(), and T_from_p_v().

◆ R_molar

constexpr static const Real IdealRealGasMixtureFluidProperties::R_molar = 8.3144598

staticconstexprprotected

molar (or universal) gas constant

Definition at line 550 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

Static Public Attributes

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]

◆ c_from_v_e() [1/2]

◆ c_from_v_e() [2/2]

◆ cp_from_p_T()

◆ cp_from_T_v()

◆ cv_from_p_T()

◆ cv_from_T_v()

◆ e_from_p_rho() [1/2]

◆ e_from_p_rho() [2/2]

◆ 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_p_v() [1/2]

◆ T_from_p_v() [2/2]

◆ 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

◆ _allow_imperfect_jacobians

◆ _emit_on_nan

◆ _fp_primary

◆ _fp_secondary

◆ _fp_secondary_names

◆ _moose_object

◆ _n_secondary_vapors

◆ _R

◆ _T_c2k

◆ _T_mix_max

◆ R_molar