From d5cbbe0c1ea5ee98ac63b664d9eddd026f71daef Mon Sep 17 00:00:00 2001 From: Ian Bell Date: Sun, 26 Apr 2015 17:36:27 -0600 Subject: [PATCH] Don't evaluate exactly at the critical point to help out with Water and CO2 reference state setting --- src/Backends/Helmholtz/Fluids/FluidLibrary.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/Backends/Helmholtz/Fluids/FluidLibrary.cpp b/src/Backends/Helmholtz/Fluids/FluidLibrary.cpp index 0cdf5a2a..931acc38 100644 --- a/src/Backends/Helmholtz/Fluids/FluidLibrary.cpp +++ b/src/Backends/Helmholtz/Fluids/FluidLibrary.cpp @@ -39,8 +39,8 @@ void JSONFluidLibrary::set_fluid_enthalpy_entropy_offset(const std::string &flui it2->second.EOS().hs_anchor.hmolar = HEOS->hmolar(); it2->second.EOS().hs_anchor.smolar = HEOS->smolar(); - // Calculate the new enthalpy and entropy values at the reducing stat - HEOS->update(DmolarT_INPUTS, it2->second.EOS().reduce.rhomolar, it2->second.EOS().reduce.T); + // Calculate the new enthalpy and entropy values at the reducing state + HEOS->update(DmolarT_INPUTS, it2->second.EOS().reduce.rhomolar+1e-12, it2->second.EOS().reduce.T+1e-12); it2->second.EOS().reduce.hmolar = HEOS->hmolar(); it2->second.EOS().reduce.smolar = HEOS->smolar(); }