Adding in HKFT parameterization
This commit is contained in:
parent
f36a5c4863
commit
2c29a41bcf
3 changed files with 15 additions and 3 deletions
|
|
@ -17,6 +17,8 @@ namespace Cantera {
|
|||
|
||||
|
||||
|
||||
|
||||
|
||||
HKFT_PDSS::HKFT_PDSS(ThermoPhase *tp, int spindex) :
|
||||
PDSS(tp, spindex)
|
||||
{
|
||||
|
|
@ -342,6 +344,13 @@ namespace Cantera {
|
|||
throw CanteraError("HKFT_PDSS::satPressure()", "unimplemented");
|
||||
return (0.0);
|
||||
}
|
||||
|
||||
|
||||
//! Internal formula for the calculation of a_g()
|
||||
double HKFT_PDSS::ag(const double temp) const {
|
||||
static double ag_coeff[3] = { -2.037662, 5.747000E-3, -6.557892E-6};
|
||||
double t2 = temp * temp;
|
||||
double val = ag_coeff[0] + ag_coeff[1] * temp * ag_coeff[2] * t2;
|
||||
return val;
|
||||
}
|
||||
|
||||
}
|
||||
|
|
|
|||
|
|
@ -134,6 +134,9 @@ namespace Cantera {
|
|||
virtual void initThermo();
|
||||
virtual void setParametersFromXML(const XML_Node& eosdata);
|
||||
|
||||
private:
|
||||
double ag(const double temp) const;
|
||||
|
||||
protected:
|
||||
|
||||
|
||||
|
|
|
|||
|
|
@ -58,13 +58,13 @@ ELECTRO_OBJ = MolalityVPSSTP.o VPStandardStateTP.o \
|
|||
WaterPropsIAPWSphi.o WaterPropsIAPWS.o WaterProps.o \
|
||||
PDSS.o WaterPDSS.o \
|
||||
HMWSoln.o HMWSoln_input.o DebyeHuckel.o \
|
||||
IdealGasPDSS.o WaterSSTP.o
|
||||
IdealGasPDSS.o WaterSSTP.o HKFT_PDSS.o
|
||||
|
||||
ELECTRO_H = MolalityVPSSTP.h VPStandardStateTP.h \
|
||||
IdealMolalSoln.h \
|
||||
WaterPropsIAPWSphi.h WaterPropsIAPWS.h WaterProps.h \
|
||||
PDSS.h WaterPDSS.h HMWSoln.h electrolytes.h \
|
||||
DebyeHuckel.h IdealGasPDSS.h WaterSSTP.h
|
||||
DebyeHuckel.h IdealGasPDSS.h WaterSSTP.h HKFT_PDSS.h
|
||||
endif
|
||||
ifeq ($(do_issp),1)
|
||||
ISSP_OBJ = IdealSolidSolnPhase.o StoichSubstanceSSTP.o SingleSpeciesTP.o
|
||||
|
|
|
|||
Loading…
Add table
Reference in a new issue