Changed calcDensity to be a virtual function started at VPStandardStateTP().
It's a protected function as well.
At the VPStandardStateTP level, calcDensity() produces an error when called, because the EOS isn't specified. Derived classes must define this function, supplying the EOS.
WaterPropsIAPW: I added documentation to this routine, defining what functions changes the state of the object.
IdealMolalSoln: I had forgotten to update this object to employ the newer setState_TP() strategy that I had used in VPStandardStateTP().
system so that I can do kinetics with the ThermoPhase object.
For molten salts, we're going with a standard concentration
of 1.0, until there is an argument otherwise. This means
that the kinetics objects, homogeneous and heterogeneous,
will treat the molten salts as unitless.
This is a vector that is a baseline "safe" mole fraction
vector that can be applied anywhere. I was running into problems
with specifying this for molten salts, because just setting
the first component to a mole fraction of 1 causes
charge imbalance.
M thermo/ThermoPhase.h
M thermo/ThermoPhase.cpp
IonsFromNeutralVPSSTP : This is initially used for molten salts.
Molten salt thermo is treated via the neutral molecule
molten salt expressions. However, ions migrate by themselves.
therefore, you must treat ions within the transport equations.
This thermophase object treats the ions as species, and is
made to be consistent with the neutral molecules thermophase
object that is used to express the non-ideal thermodynamics
of the molten salt solution. Currently, MargulesVPSSTP
is used for this as the LiKCl system is parameterized using
this excess Gibbs free energy formulation.
MargulesVPSSTP
Various upgrades and fillers. It now reads and is initiated
from the input file.
setState_TP(), setTemperature(), and setPressure() within
VPStandardStateTP() formulations
There was some circular logic, previously. Now setTemperature() and setPressure()
are subserviant to setState_TP(). Basically, it doesn't make
sense to calculate the standard state without first getting
both the temperature and pressure (not just one and then the other)
before you calculate the standard state for VPStandardStateTP derivatives.
Think about the water EOS, which is a real fluid. You don't want
to be trying to calculate the standard state within the enveloppe where
neither gas nor liquid is stable.
PseudoBinaryVPSSTP
Deprecated.
