A species thermodynamic property manager for a phase. More...

#include <MultiSpeciesThermo.h>

Detailed Description

A species thermodynamic property manager for a phase.

This is a general manager that can handle a wide variety of species thermodynamic polynomials for individual species and compute their nondimensional, reference-state thermodynamic properties (that is, as a function of temperature only).

The ThermoPhase object relies on MultiSpeciesThermo to calculate the thermodynamic properties of the reference state for all of the species in the phase, for a range of temperatures. Note, the pressure dependence of the species thermodynamic functions is not handled at this level. Species using the same parameterization are grouped together in order to minimize the operation count and achieve better efficiency.

The most important member function for the MultiSpeciesThermo class is the member function MultiSpeciesThermo::update(). The function calculates the values of Cp/R, H/RT, and S/R for all of the species at once at the specified temperature.

Usually, all of the species in a phase are installed into a MultiSpeciesThermo object. However, there is no requirement that a MultiSpeciesThermo object handles all of the species in a phase. The member function install_STIT() is called to install each species into the MultiSpeciesThermo object.

Definition at line 46 of file MultiSpeciesThermo.h.

Public Member Functions
	MultiSpeciesThermo ()=default
	Constructor.

	MultiSpeciesThermo (const MultiSpeciesThermo &b)=delete

MultiSpeciesThermo &	operator= (const MultiSpeciesThermo &b)=delete

virtual void	install_STIT (size_t index, shared_ptr< SpeciesThermoInterpType > stit)
	Install a new species thermodynamic property parameterization for one species.

virtual void	modifySpecies (size_t index, shared_ptr< SpeciesThermoInterpType > spec)
	Modify the species thermodynamic property parameterization for a species.

virtual void	update_single (size_t k, double T, double cp_R, double h_RT, double *s_R) const
	Get reference-state properties for a single species.

virtual void	update (double T, double cp_R, double h_RT, double *s_R) const
	Compute the reference-state properties for all species.

virtual double	minTemp (size_t k=npos) const
	Minimum temperature.

virtual double	maxTemp (size_t k=npos) const
	Maximum temperature.

virtual double	refPressure (size_t k=npos) const
	The reference-state pressure for species k.

virtual int	reportType (size_t index) const
	This utility function reports the type of parameterization used for the species with index number index.

virtual void	reportParams (size_t index, int &type, double *const c, double &minTemp, double &maxTemp, double &refPressure) const
	This utility function reports back the type of parameterization and all of the parameters for the species with index number index.

virtual double	reportOneHf298 (const size_t k) const
	Report the 298 K Heat of Formation of the standard state of one species (J kmol-1)

virtual void	modifyOneHf298 (const size_t k, const double Hf298New)
	Modify the value of the 298 K Heat of Formation of the standard state of one species in the phase (J kmol-1)

virtual void	resetHf298 (const size_t k)
	Restore the original heat of formation of one or more species.

bool	ready (size_t nSpecies)
	Check if data for all species (0 through nSpecies-1) has been installed.

Protected Types
typedef pair< size_t, shared_ptr< SpeciesThermoInterpType > >	index_STIT

typedef map< int, vector< index_STIT > >	STIT_map

typedef map< int, vector< double > >	tpoly_map

Protected Member Functions
void	markInstalled (size_t k)
	Mark species k as having its thermodynamic data installed.

Protected Attributes
STIT_map	m_sp
	This is the main data structure, which contains the SpeciesThermoInterpType objects, sorted by the parameterization type.

tpoly_map	m_tpoly
	Temperature polynomials for each thermo parameterization.

map< size_t, pair< int, size_t > >	m_speciesLoc
	Map from species index to location within m_sp, such that `m_sp[m_speciesLoc[k].first][m_speciesLoc[k].second]` is the SpeciesThermoInterpType object for species `k`.

double	m_tlow_max = 0.0
	Maximum value of the lowest temperature.

double	m_thigh_min = 1e+30
	Minimum value of the highest temperature.

double	m_p0 = 0.0
	reference pressure (Pa)

vector< bool >	m_installed
	indicates if data for species has been installed

Private Member Functions
SpeciesThermoInterpType *	provideSTIT (size_t k)
	Provide the SpeciesThermoInterpType object.

const SpeciesThermoInterpType *	provideSTIT (size_t k) const

Member Typedef Documentation

◆ index_STIT

typedef pair<size_t, shared_ptr<SpeciesThermoInterpType> > index_STIT

protected

Definition at line 202 of file MultiSpeciesThermo.h.

◆ STIT_map

typedef map<int, vector<index_STIT> > STIT_map

protected

Definition at line 203 of file MultiSpeciesThermo.h.

◆ tpoly_map

typedef map<int, vector<double> > tpoly_map

protected

Definition at line 204 of file MultiSpeciesThermo.h.

Constructor & Destructor Documentation

◆ MultiSpeciesThermo()

MultiSpeciesThermo ( )

default

Constructor.

Member Function Documentation

◆ install_STIT()

void install_STIT	(	size_t	index,
		shared_ptr< SpeciesThermoInterpType >	stit
	)

virtual

Install a new species thermodynamic property parameterization for one species.

Parameters

index	Index of the species being installed
stit	Pointer to the SpeciesThermoInterpType object This will set up the thermo for one species

Definition at line 21 of file MultiSpeciesThermo.cpp.

◆ modifySpecies()

void modifySpecies	(	size_t	index,
		shared_ptr< SpeciesThermoInterpType >	spec
	)

virtual

Modify the species thermodynamic property parameterization for a species.

Parameters

index	Index of the species being installed
spec	Pointer to the SpeciesThermoInterpType object

Definition at line 53 of file MultiSpeciesThermo.cpp.

◆ update_single()

void update_single	(	size_t	k,
		double	T,
		double *	cp_R,
		double *	h_RT,
		double *	s_R
	)		const

virtual

Get reference-state properties for a single species.

Parameters

k	species index
T	Temperature (Kelvin)
cp_R	Dimensionless heat capacity
h_RT	Dimensionless enthalpy
s_R	Dimensionless entropy

Definition at line 85 of file MultiSpeciesThermo.cpp.

◆ update()

void update	(	double	T,
		double *	cp_R,
		double *	h_RT,
		double *	s_R
	)		const

virtual

Compute the reference-state properties for all species.

Given temperature T in K, this method updates the values of the non- dimensional heat capacity at constant pressure, enthalpy, and entropy, at the reference pressure, Pref of each of the standard states.

Parameters

T	Temperature (Kelvin)
cp_R	Vector of Dimensionless heat capacities. (length m_kk).
h_RT	Vector of Dimensionless enthalpies. (length m_kk).
s_R	Vector of Dimensionless entropies. (length m_kk).

Definition at line 94 of file MultiSpeciesThermo.cpp.

◆ minTemp()

double minTemp ( size_t k = npos ) const

virtual

Minimum temperature.

If no argument is supplied, this method returns the minimum temperature for which all parameterizations are valid. If an integer index k is supplied, then the value returned is the minimum temperature for species k in the phase.

Parameters

k	Species index

Definition at line 130 of file MultiSpeciesThermo.cpp.

◆ maxTemp()

double maxTemp ( size_t k = npos ) const

virtual

Maximum temperature.

If no argument is supplied, this method returns the maximum temperature for which all parameterizations are valid. If an integer index k is supplied, then the value returned is the maximum temperature for parameterization k.

Parameters

k	Species Index

Definition at line 141 of file MultiSpeciesThermo.cpp.

◆ refPressure()

double refPressure ( size_t k = npos ) const

virtual

The reference-state pressure for species k.

Returns the reference state pressure in Pascals for species k. If k is left out of the argument list, it returns the reference state pressure for the first species.

Parameters

k	Species Index

Deprecated:: The species index parameter is deprecated and will be removed after Cantera 3.0. All species in a phase must have the same reference pressure.

Definition at line 152 of file MultiSpeciesThermo.cpp.

◆ reportType()

int reportType ( size_t index ) const

virtual

This utility function reports the type of parameterization used for the species with index number index.

Parameters

index Species index

Definition at line 108 of file MultiSpeciesThermo.cpp.

◆ reportParams()

void reportParams	(	size_t	index,
		int &	type,
		double *const	c,
		double &	minTemp,
		double &	maxTemp,
		double &	refPressure
	)		const

virtual

This utility function reports back the type of parameterization and all of the parameters for the species with index number index.

Parameters

index	Species index
type	Integer type of the standard type
c	Vector of coefficients used to set the parameters for the standard state.
minTemp	output - Minimum temperature
maxTemp	output - Maximum temperature
refPressure	output - reference pressure (Pa).

Definition at line 117 of file MultiSpeciesThermo.cpp.

◆ reportOneHf298()

double reportOneHf298 ( const size_t k ) const

virtual

Report the 298 K Heat of Formation of the standard state of one species (J kmol-1)

The 298K Heat of Formation is defined as the enthalpy change to create the standard state of the species from its constituent elements in their standard states at 298 K and 1 bar.

Parameters

k	species index

Returns: the current value of the Heat of Formation at 298K and 1 bar

Definition at line 186 of file MultiSpeciesThermo.cpp.

◆ modifyOneHf298()

void modifyOneHf298	(	const size_t	k,
		const double	Hf298New
	)

virtual

Modify the value of the 298 K Heat of Formation of the standard state of one species in the phase (J kmol-1)

The 298K heat of formation is defined as the enthalpy change to create the standard state of the species from its constituent elements in their standard states at 298 K and 1 bar.

Parameters

k	Index of the species
Hf298New	Specify the new value of the Heat of Formation at 298K and 1 bar. units = J/kmol.

Definition at line 196 of file MultiSpeciesThermo.cpp.

◆ resetHf298()

void resetHf298 ( const size_t k )

virtual

Restore the original heat of formation of one or more species.

Resets changes made by modifyOneHf298(). If the species index is not specified, the heats of formation for all species are restored.

Definition at line 204 of file MultiSpeciesThermo.cpp.

◆ ready()

bool ready ( size_t nSpecies )

Check if data for all species (0 through nSpecies-1) has been installed.

Definition at line 212 of file MultiSpeciesThermo.cpp.

◆ provideSTIT() [1/2]

SpeciesThermoInterpType * provideSTIT ( size_t k )

private

Provide the SpeciesThermoInterpType object.

Parameters

k	species index

Returns: pointer to the SpeciesThermoInterpType object.

Definition at line 166 of file MultiSpeciesThermo.cpp.

◆ provideSTIT() [2/2]

const SpeciesThermoInterpType * provideSTIT ( size_t k ) const

private

Definition at line 176 of file MultiSpeciesThermo.cpp.

◆ markInstalled()

void markInstalled ( size_t k )

protected

Mark species k as having its thermodynamic data installed.

Definition at line 224 of file MultiSpeciesThermo.cpp.

Member Data Documentation

◆ m_sp

STIT_map m_sp

protected

This is the main data structure, which contains the SpeciesThermoInterpType objects, sorted by the parameterization type.

m_sp[i] is the vector of [species index, STIT] pairs which use parameterization i.

Definition at line 210 of file MultiSpeciesThermo.h.

◆ m_tpoly

tpoly_map m_tpoly

mutableprotected

Temperature polynomials for each thermo parameterization.

Definition at line 213 of file MultiSpeciesThermo.h.

◆ m_speciesLoc

map<size_t, pair<int, size_t> > m_speciesLoc

protected

Map from species index to location within m_sp, such that m_sp[m_speciesLoc[k].first][m_speciesLoc[k].second] is the SpeciesThermoInterpType object for species k.

Definition at line 218 of file MultiSpeciesThermo.h.

◆ m_tlow_max

double m_tlow_max = 0.0

protected

Maximum value of the lowest temperature.

Definition at line 221 of file MultiSpeciesThermo.h.

◆ m_thigh_min

double m_thigh_min = 1e+30

protected

Minimum value of the highest temperature.

Definition at line 224 of file MultiSpeciesThermo.h.

◆ m_p0

double m_p0 = 0.0

protected

reference pressure (Pa)

Definition at line 227 of file MultiSpeciesThermo.h.

◆ m_installed

vector<bool> m_installed

protected

indicates if data for species has been installed

Definition at line 230 of file MultiSpeciesThermo.h.

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

Detailed Description

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Member Typedef Documentation

◆ index_STIT

◆ STIT_map

◆ tpoly_map

Constructor & Destructor Documentation

◆ MultiSpeciesThermo()

Member Function Documentation

◆ install_STIT()

◆ modifySpecies()

◆ update_single()

◆ update()

◆ minTemp()

◆ maxTemp()

◆ refPressure()

◆ reportType()

◆ reportParams()

◆ reportOneHf298()

◆ modifyOneHf298()

◆ resetHf298()

◆ ready()

◆ provideSTIT() [1/2]

◆ provideSTIT() [2/2]

◆ markInstalled()

Member Data Documentation

◆ m_sp

◆ m_tpoly

◆ m_speciesLoc

◆ m_tlow_max

◆ m_thigh_min

◆ m_p0

◆ m_installed