iapws-if97 2.0.2 → 2.1.1

package/dist/saturation/common.js

@@ -1,11 +1,12 @@
 import { Region } from '../types.js';
-import { Pc, R2_T_MIN, RHOc } from '../constants.js';
+import { Pc, R2_T_MIN, RHOc, Tt } from '../constants.js';
 import { region1 } from '../regions/region1.js';
 import { region2 } from '../regions/region2.js';
 import { region3ByRhoT } from '../regions/region3.js';
 import { saturationPressure, saturationTemperature } from '../regions/region4.js';
 import { region3SatVolume } from '../regions/region3-subregions.js';
 import { newtonRaphson } from '../solvers/newton-raphson.js';
+import { normalizeRegion4Pressure, normalizeRegion4Temperature, } from './region4-boundaries.js';
 /** Tag a single-phase state as Region 4 with the given vapour quality. */
 function withRegion4Metadata(state, quality) {
     return { ...state, region: Region.Region4, quality };
@@ -38,24 +39,25 @@ function clampQuality(x) {
  * @param p - Saturation pressure [MPa]
  */
 export function saturationEndpointsAtPressure(p) {
-    const temperature = saturationTemperature(p);
-    if (Math.abs(p - Pc) < 1e-10) {
+    const pressure = normalizeRegion4Pressure(p);
+    const temperature = normalizeRegion4Temperature(saturationTemperature(pressure));
+    if (pressure === Pc) {
         const state = region3ByRhoT(RHOc, temperature);
-        return { pressure: p, temperature, liquid: state, vapor: state };
+        return { pressure, temperature, liquid: state, vapor: state };
     }
     if (temperature > R2_T_MIN) {
         return {
-            pressure: p,
+            pressure,
             temperature,
-            liquid: region3ByRhoT(solveR3Density(p, temperature, 0), temperature),
-            vapor: region3ByRhoT(solveR3Density(p, temperature, 1), temperature),
+            liquid: region3ByRhoT(solveR3Density(pressure, temperature, 0), temperature),
+            vapor: region3ByRhoT(solveR3Density(pressure, temperature, 1), temperature),
         };
     }
     return {
-        pressure: p,
+        pressure,
         temperature,
-        liquid: region1(p, temperature),
-        vapor: region2(p, temperature),
+        liquid: region1(pressure, temperature),
+        vapor: region2(pressure, temperature),
     };
 }
 /**
@@ -63,7 +65,17 @@ export function saturationEndpointsAtPressure(p) {
  * @param T - Saturation temperature [K]
  */
 export function saturationEndpointsAtTemperature(T) {
-    return saturationEndpointsAtPressure(saturationPressure(T));
+    const temperature = normalizeRegion4Temperature(T);
+    const pressure = normalizeRegion4Pressure(saturationPressure(temperature));
+    if (temperature < Tt) {
+        return {
+            pressure,
+            temperature,
+            liquid: region1(pressure, temperature),
+            vapor: region2(pressure, temperature),
+        };
+    }
+    return saturationEndpointsAtPressure(pressure);
 }
 /**
  * Compute vapour quality from any extensive saturation property.

package/dist/saturation/common.js.map

@@ -1 +1 @@
-{"version":3,"file":"common.js","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,EAAE,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,iBAAiB,CAAC;AACrD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AACtD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,uBAAuB,CAAC;AAClF,OAAO,EAAE,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAU7D,0EAA0E;AAC1E,SAAS,mBAAmB,CAAC,KAAsB,EAAE,OAAe;IAClE,OAAO,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC;AACvD,CAAC;AAED;;;;;;GAMG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAQ;IACpD,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,OAAO,aAAa,CAClB,CAAC,GAAG,EAAE,EAAE,CAAC,aAAa,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,EAC3C,CAAC,GAAG,EAAE,CACP,CAAC;AACJ,CAAC;AAED,qFAAqF;AACrF,SAAS,YAAY,CAAC,CAAS;IAC7B,IAAI,CAAC,IAAI,KAAK;QAAE,OAAO,CAAC,CAAC;IACzB,IAAI,CAAC,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC;IAC7B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,6BAA6B,CAAC,CAAS;IACrD,MAAM,WAAW,GAAG,qBAAqB,CAAC,CAAC,CAAC,CAAC;IAE7C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC;QAC7B,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,EAAE,WAAW,CAAC,CAAC;QAC/C,OAAO,EAAE,QAAQ,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IACnE,CAAC;IAED,IAAI,WAAW,GAAG,QAAQ,EAAE,CAAC;QAC3B,OAAO;YACL,QAAQ,EAAE,CAAC;YACX,WAAW;YACX,MAAM,EAAE,aAAa,CAAC,cAAc,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;YACrE,KAAK,EAAE,aAAa,CAAC,cAAc,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;SACrE,CAAC;IACJ,CAAC;IAED,OAAO;QACL,QAAQ,EAAE,CAAC;QACX,WAAW;QACX,MAAM,EAAE,OAAO,CAAC,CAAC,EAAE,WAAW,CAAC;QAC/B,KAAK,EAAE,OAAO,CAAC,CAAC,EAAE,WAAW,CAAC;KAC/B,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gCAAgC,CAAC,CAAS;IACxD,OAAO,6BAA6B,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,6BAA6B,CAC3C,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,OAAO,YAAY,CAAC,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAChC,SAA8B,EAC9B,YAAoB;IAEpB,MAAM,OAAO,GAAG,YAAY,CAAC,YAAY,CAAC,CAAC;IAC3C,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,SAAS,CAAC;IAE3D,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC;IAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACvC,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ;QACR,WAAW;QACX,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,OAAO,EAAE,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,OAAO,GAAG,OAAO;QACjE,QAAQ,EAAE,MAAM,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,OAAO;QACpE,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,IAAI;QACR,YAAY,EAAE,IAAI;QAClB,OAAO;QACP,iBAAiB,EAAE,IAAI;QACvB,yBAAyB,EAAE,IAAI;KAChC,CAAC;AACJ,CAAC"}
+{"version":3,"file":"common.js","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,EAAE,EAAE,QAAQ,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,iBAAiB,CAAC;AACzD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AACtD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,uBAAuB,CAAC;AAClF,OAAO,EAAE,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAC7D,OAAO,EACL,wBAAwB,EACxB,2BAA2B,GAC5B,MAAM,yBAAyB,CAAC;AAUjC,0EAA0E;AAC1E,SAAS,mBAAmB,CAAC,KAAsB,EAAE,OAAe;IAClE,OAAO,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC;AACvD,CAAC;AAED;;;;;;GAMG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAQ;IACpD,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,OAAO,aAAa,CAClB,CAAC,GAAG,EAAE,EAAE,CAAC,aAAa,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,EAC3C,CAAC,GAAG,EAAE,CACP,CAAC;AACJ,CAAC;AAED,qFAAqF;AACrF,SAAS,YAAY,CAAC,CAAS;IAC7B,IAAI,CAAC,IAAI,KAAK;QAAE,OAAO,CAAC,CAAC;IACzB,IAAI,CAAC,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC;IAC7B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,6BAA6B,CAAC,CAAS;IACrD,MAAM,QAAQ,GAAG,wBAAwB,CAAC,CAAC,CAAC,CAAC;IAC7C,MAAM,WAAW,GAAG,2BAA2B,CAAC,qBAAqB,CAAC,QAAQ,CAAC,CAAC,CAAC;IAEjF,IAAI,QAAQ,KAAK,EAAE,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,EAAE,WAAW,CAAC,CAAC;QAC/C,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,IAAI,WAAW,GAAG,QAAQ,EAAE,CAAC;QAC3B,OAAO;YACL,QAAQ;YACR,WAAW;YACX,MAAM,EAAE,aAAa,CAAC,cAAc,CAAC,QAAQ,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;YAC5E,KAAK,EAAE,aAAa,CAAC,cAAc,CAAC,QAAQ,EAAE,WAAW,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC;SAC5E,CAAC;IACJ,CAAC;IAED,OAAO;QACL,QAAQ;QACR,WAAW;QACX,MAAM,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;QACtC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;KACtC,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gCAAgC,CAAC,CAAS;IACxD,MAAM,WAAW,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,wBAAwB,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC;IAE3E,IAAI,WAAW,GAAG,EAAE,EAAE,CAAC;QACrB,OAAO;YACL,QAAQ;YACR,WAAW;YACX,MAAM,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;YACtC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,WAAW,CAAC;SACtC,CAAC;IACJ,CAAC;IAED,OAAO,6BAA6B,CAAC,QAAQ,CAAC,CAAC;AACjD,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,6BAA6B,CAC3C,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,OAAO,YAAY,CAAC,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAChC,SAA8B,EAC9B,YAAoB;IAEpB,MAAM,OAAO,GAAG,YAAY,CAAC,YAAY,CAAC,CAAC;IAC3C,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,SAAS,CAAC;IAE3D,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC;IAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACvC,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ;QACR,WAAW;QACX,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,cAAc,EAAE,MAAM,CAAC,cAAc,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,cAAc,GAAG,OAAO;QACtF,OAAO,EAAE,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,OAAO,GAAG,OAAO;QACjE,QAAQ,EAAE,MAAM,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,OAAO;QACpE,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,IAAI;QACR,YAAY,EAAE,IAAI;QAClB,OAAO;QACP,iBAAiB,EAAE,IAAI;QACvB,yBAAyB,EAAE,IAAI;KAChC,CAAC;AACJ,CAAC"}

package/dist/saturation/region4-boundaries.d.ts

@@ -0,0 +1,14 @@
+export declare function normalizeRegion4Pressure(p: number): number;
+export declare function normalizeRegion4Temperature(T: number): number;
+export declare function clampRegion4TemperatureBelowCritical(T: number): number;
+export declare function isRegion4CriticalPressure(p: number): boolean;
+export declare function isRegion4CriticalTemperature(T: number): boolean;
+export declare function assertRegion4PressureStateInput(p: number, solverName: string): number;
+export declare function assertRegion4TemperatureStateInput(T: number, solverName: string): number;
+export declare function assertRegion4StateAllowed(pressure: number, temperature: number, solverName: string): void;
+export declare function isCriticalRegion4Enthalpy(h: number): boolean;
+export declare function isCriticalRegion4Entropy(s: number): boolean;
+export declare function assertCriticalRegion4PHInput(p: number, h: number, solverName: string): void;
+export declare function assertCriticalRegion4PSInput(p: number, s: number, solverName: string): void;
+export declare function assertCriticalRegion4HSInput(h: number, s: number, solverName: string): void;
+//# sourceMappingURL=region4-boundaries.d.ts.map

package/dist/saturation/region4-boundaries.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-boundaries.d.ts","sourceRoot":"","sources":["../../src/saturation/region4-boundaries.ts"],"names":[],"mappings":"AAGA,wBAAgB,wBAAwB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK1D;AAED,wBAAgB,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK7D;AAED,wBAAgB,oCAAoC,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtE;AAED,wBAAgB,yBAAyB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE5D;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE/D;AAMD,wBAAgB,+BAA+B,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM,CASrF;AAED,wBAAgB,kCAAkC,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM,CASxF;AAED,wBAAgB,yBAAyB,CACvC,QAAQ,EAAE,MAAM,EAChB,WAAW,EAAE,MAAM,EACnB,UAAU,EAAE,MAAM,GACjB,IAAI,CAIN;AAED,wBAAgB,yBAAyB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE5D;AAED,wBAAgB,wBAAwB,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO,CAE3D;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F;AAED,wBAAgB,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,IAAI,CAI3F"}

package/dist/saturation/region4-boundaries.js

@@ -0,0 +1,79 @@
+import * as C from '../constants.js';
+import { IF97Error, OutOfRangeError } from '../types.js';
+export function normalizeRegion4Pressure(p) {
+    if (Math.abs(p - C.P_MIN) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.P_MIN;
+    if (Math.abs(p - C.Pt) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.Pt;
+    if (Math.abs(p - C.Pc) <= C.REGION4_PRESSURE_TOLERANCE)
+        return C.Pc;
+    return p;
+}
+export function normalizeRegion4Temperature(T) {
+    if (Math.abs(T - C.T_MIN) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.T_MIN;
+    if (Math.abs(T - C.Tt) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.Tt;
+    if (Math.abs(T - C.Tc) <= C.REGION4_TEMPERATURE_TOLERANCE)
+        return C.Tc;
+    return T;
+}
+export function clampRegion4TemperatureBelowCritical(T) {
+    return Math.max(C.Tt, Math.min(T, C.Tc - C.REGION4_SUBCRITICAL_TEMPERATURE_MARGIN));
+}
+export function isRegion4CriticalPressure(p) {
+    return normalizeRegion4Pressure(p) === C.Pc;
+}
+export function isRegion4CriticalTemperature(T) {
+    return normalizeRegion4Temperature(T) === C.Tc;
+}
+function criticalPointMessage(solverName) {
+    return `${solverName} does not support the critical point P=${C.Pc} MPa, T=${C.Tc} K as a Region 4 saturation state`;
+}
+export function assertRegion4PressureStateInput(p, solverName) {
+    const pressure = normalizeRegion4Pressure(p);
+    if (pressure < C.Pt || pressure > C.Pc) {
+        throw new OutOfRangeError('Pressure', p, C.Pt, C.Pc);
+    }
+    if (pressure === C.Pc) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+    return pressure;
+}
+export function assertRegion4TemperatureStateInput(T, solverName) {
+    const temperature = normalizeRegion4Temperature(T);
+    if (temperature < C.Tt || temperature > C.Tc) {
+        throw new OutOfRangeError('Temperature', T, C.Tt, C.Tc);
+    }
+    if (temperature === C.Tc) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+    return temperature;
+}
+export function assertRegion4StateAllowed(pressure, temperature, solverName) {
+    if (isRegion4CriticalPressure(pressure) || isRegion4CriticalTemperature(temperature)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function isCriticalRegion4Enthalpy(h) {
+    return Math.abs(h - C.R3_H_CRT) <= C.REGION4_CRITICAL_HS_TOLERANCE;
+}
+export function isCriticalRegion4Entropy(s) {
+    return Math.abs(s - C.R3_S_CRT) <= C.REGION4_CRITICAL_HS_TOLERANCE;
+}
+export function assertCriticalRegion4PHInput(p, h, solverName) {
+    if (isRegion4CriticalPressure(p) && isCriticalRegion4Enthalpy(h)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function assertCriticalRegion4PSInput(p, s, solverName) {
+    if (isRegion4CriticalPressure(p) && isCriticalRegion4Entropy(s)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+export function assertCriticalRegion4HSInput(h, s, solverName) {
+    if (isCriticalRegion4Enthalpy(h) && isCriticalRegion4Entropy(s)) {
+        throw new IF97Error(criticalPointMessage(solverName));
+    }
+}
+//# sourceMappingURL=region4-boundaries.js.map

package/dist/saturation/region4-boundaries.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-boundaries.js","sourceRoot":"","sources":["../../src/saturation/region4-boundaries.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,iBAAiB,CAAC;AACrC,OAAO,EAAE,SAAS,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAEzD,MAAM,UAAU,wBAAwB,CAAC,CAAS;IAChD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,KAAK,CAAC;IAC1E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACpE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,0BAA0B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACpE,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,2BAA2B,CAAC,CAAS;IACnD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,KAAK,CAAC;IAC7E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACvE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,6BAA6B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC;IACvE,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,oCAAoC,CAAC,CAAS;IAC5D,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,sCAAsC,CAAC,CAAC,CAAC;AACtF,CAAC;AAED,MAAM,UAAU,yBAAyB,CAAC,CAAS;IACjD,OAAO,wBAAwB,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;AAC9C,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS;IACpD,OAAO,2BAA2B,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;AACjD,CAAC;AAED,SAAS,oBAAoB,CAAC,UAAkB;IAC9C,OAAO,GAAG,UAAU,0CAA0C,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,EAAE,mCAAmC,CAAC;AACvH,CAAC;AAED,MAAM,UAAU,+BAA+B,CAAC,CAAS,EAAE,UAAkB;IAC3E,MAAM,QAAQ,GAAG,wBAAwB,CAAC,CAAC,CAAC,CAAC;IAC7C,IAAI,QAAQ,GAAG,CAAC,CAAC,EAAE,IAAI,QAAQ,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,IAAI,eAAe,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IACvD,CAAC;IACD,IAAI,QAAQ,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED,MAAM,UAAU,kCAAkC,CAAC,CAAS,EAAE,UAAkB;IAC9E,MAAM,WAAW,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IACnD,IAAI,WAAW,GAAG,CAAC,CAAC,EAAE,IAAI,WAAW,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC;QAC7C,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAC1D,CAAC;IACD,IAAI,WAAW,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;QACzB,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;IACD,OAAO,WAAW,CAAC;AACrB,CAAC;AAED,MAAM,UAAU,yBAAyB,CACvC,QAAgB,EAChB,WAAmB,EACnB,UAAkB;IAElB,IAAI,yBAAyB,CAAC,QAAQ,CAAC,IAAI,4BAA4B,CAAC,WAAW,CAAC,EAAE,CAAC;QACrF,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,yBAAyB,CAAC,CAAS;IACjD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,6BAA6B,CAAC;AACrE,CAAC;AAED,MAAM,UAAU,wBAAwB,CAAC,CAAS;IAChD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,6BAA6B,CAAC;AACrE,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,yBAAyB,CAAC,CAAC,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,MAAM,UAAU,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,UAAkB;IACnF,IAAI,yBAAyB,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,SAAS,CAAC,oBAAoB,CAAC,UAAU,CAAC,CAAC,CAAC;IACxD,CAAC;AACH,CAAC"}

package/dist/saturation/region4-hs.d.ts

@@ -0,0 +1,3 @@
+import type { BasicProperties } from '../types.js';
+export declare function tryRegion4HSState(h: number, s: number): BasicProperties | null;
+//# sourceMappingURL=region4-hs.d.ts.map

package/dist/saturation/region4-hs.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-hs.d.ts","sourceRoot":"","sources":["../../src/saturation/region4-hs.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAkFnD,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,GAAG,IAAI,CAkD9E"}

package/dist/saturation/region4-hs.js

@@ -0,0 +1,82 @@
+import * as C from '../constants.js';
+import { bracketedNewton } from '../solvers/bracketed-newton.js';
+import { mixSaturationState, saturationEndpointsAtTemperature } from './common.js';
+import { assertRegion4StateAllowed, clampRegion4TemperatureBelowCritical, } from './region4-boundaries.js';
+const REGION4_QUALITY_TOLERANCE = 1e-6;
+function isClose(a, b) {
+    return Math.abs(a - b) <= 1e-9 * Math.max(1, Math.abs(a), Math.abs(b));
+}
+function rawQualityFromSaturationProperty(liquidValue, vaporValue, mixtureValue) {
+    return (mixtureValue - liquidValue) / (vaporValue - liquidValue);
+}
+function isAdmissibleQuality(quality) {
+    return Number.isFinite(quality)
+        && quality >= -REGION4_QUALITY_TOLERANCE
+        && quality <= 1 + REGION4_QUALITY_TOLERANCE;
+}
+function entropyResidual(h, s, temperatureGuess) {
+    const endpoints = saturationEndpointsAtTemperature(clampRegion4TemperatureBelowCritical(temperatureGuess));
+    const quality = rawQualityFromSaturationProperty(endpoints.liquid.enthalpy, endpoints.vapor.enthalpy, h);
+    if (!Number.isFinite(quality)) {
+        return Number.NaN;
+    }
+    return endpoints.liquid.entropy * (1 - quality) + endpoints.vapor.entropy * quality - s;
+}
+function finalizeRegion4HSState(h, s, temperature) {
+    const endpoints = saturationEndpointsAtTemperature(temperature);
+    const qualityFromEnthalpy = rawQualityFromSaturationProperty(endpoints.liquid.enthalpy, endpoints.vapor.enthalpy, h);
+    const qualityFromEntropy = rawQualityFromSaturationProperty(endpoints.liquid.entropy, endpoints.vapor.entropy, s);
+    if (!isAdmissibleQuality(qualityFromEnthalpy) || !isAdmissibleQuality(qualityFromEntropy)) {
+        return null;
+    }
+    if (Math.abs(qualityFromEnthalpy - qualityFromEntropy) > REGION4_QUALITY_TOLERANCE) {
+        return null;
+    }
+    if (Math.abs(entropyResidual(h, s, temperature)) > C.REGION4_HS_RESIDUAL_TOLERANCE) {
+        return null;
+    }
+    try {
+        assertRegion4StateAllowed(endpoints.pressure, endpoints.temperature, 'Region4HS');
+    }
+    catch {
+        return null;
+    }
+    return mixSaturationState(endpoints, qualityFromEnthalpy);
+}
+export function tryRegion4HSState(h, s) {
+    const lower = C.Tt;
+    const upper = clampRegion4TemperatureBelowCritical(C.Tc);
+    const lowerState = finalizeRegion4HSState(h, s, lower);
+    if (lowerState !== null) {
+        return lowerState;
+    }
+    const upperState = finalizeRegion4HSState(h, s, upper);
+    if (upperState !== null) {
+        return upperState;
+    }
+    let previousTemperature = lower;
+    let previousResidual = entropyResidual(h, s, lower);
+    for (let i = 1; i <= C.REGION4_HS_BRACKET_SEGMENTS; i++) {
+        const temperature = lower + (upper - lower) * (i / C.REGION4_HS_BRACKET_SEGMENTS);
+        const residual = entropyResidual(h, s, temperature);
+        if (!Number.isFinite(previousResidual) || !Number.isFinite(residual)) {
+            previousTemperature = temperature;
+            previousResidual = residual;
+            continue;
+        }
+        if (isClose(previousResidual, 0)) {
+            return finalizeRegion4HSState(h, s, previousTemperature);
+        }
+        if (isClose(residual, 0)) {
+            return finalizeRegion4HSState(h, s, temperature);
+        }
+        if (previousResidual * residual < 0) {
+            const root = bracketedNewton((candidateTemperature) => entropyResidual(h, s, candidateTemperature), previousTemperature, temperature, (previousTemperature + temperature) / 2, { tolerance: C.REGION4_HS_RESIDUAL_TOLERANCE });
+            return finalizeRegion4HSState(h, s, root);
+        }
+        previousTemperature = temperature;
+        previousResidual = residual;
+    }
+    return null;
+}
+//# sourceMappingURL=region4-hs.js.map

package/dist/saturation/region4-hs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4-hs.js","sourceRoot":"","sources":["../../src/saturation/region4-hs.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,iBAAiB,CAAC;AACrC,OAAO,EAAE,eAAe,EAAE,MAAM,gCAAgC,CAAC;AAEjE,OAAO,EAAE,kBAAkB,EAAE,gCAAgC,EAAE,MAAM,aAAa,CAAC;AACnF,OAAO,EACL,yBAAyB,EACzB,oCAAoC,GACrC,MAAM,yBAAyB,CAAC;AAEjC,MAAM,yBAAyB,GAAG,IAAI,CAAC;AAEvC,SAAS,OAAO,CAAC,CAAS,EAAE,CAAS;IACnC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACzE,CAAC;AAED,SAAS,gCAAgC,CACvC,WAAmB,EACnB,UAAkB,EAClB,YAAoB;IAEpB,OAAO,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,CAAC;AACnE,CAAC;AAED,SAAS,mBAAmB,CAAC,OAAe;IAC1C,OAAO,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;WAC1B,OAAO,IAAI,CAAC,yBAAyB;WACrC,OAAO,IAAI,CAAC,GAAG,yBAAyB,CAAC;AAChD,CAAC;AAED,SAAS,eAAe,CAAC,CAAS,EAAE,CAAS,EAAE,gBAAwB;IACrE,MAAM,SAAS,GAAG,gCAAgC,CAChD,oCAAoC,CAAC,gBAAgB,CAAC,CACvD,CAAC;IACF,MAAM,OAAO,GAAG,gCAAgC,CAC9C,SAAS,CAAC,MAAM,CAAC,QAAQ,EACzB,SAAS,CAAC,KAAK,CAAC,QAAQ,EACxB,CAAC,CACF,CAAC;IAEF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC9B,OAAO,MAAM,CAAC,GAAG,CAAC;IACpB,CAAC;IAED,OAAO,SAAS,CAAC,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,OAAO,GAAG,OAAO,GAAG,CAAC,CAAC;AAC1F,CAAC;AAED,SAAS,sBAAsB,CAC7B,CAAS,EACT,CAAS,EACT,WAAmB;IAEnB,MAAM,SAAS,GAAG,gCAAgC,CAAC,WAAW,CAAC,CAAC;IAChE,MAAM,mBAAmB,GAAG,gCAAgC,CAC1D,SAAS,CAAC,MAAM,CAAC,QAAQ,EACzB,SAAS,CAAC,KAAK,CAAC,QAAQ,EACxB,CAAC,CACF,CAAC;IACF,MAAM,kBAAkB,GAAG,gCAAgC,CACzD,SAAS,CAAC,MAAM,CAAC,OAAO,EACxB,SAAS,CAAC,KAAK,CAAC,OAAO,EACvB,CAAC,CACF,CAAC;IAEF,IAAI,CAAC,mBAAmB,CAAC,mBAAmB,CAAC,IAAI,CAAC,mBAAmB,CAAC,kBAAkB,CAAC,EAAE,CAAC;QAC1F,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,IAAI,CAAC,GAAG,CAAC,mBAAmB,GAAG,kBAAkB,CAAC,GAAG,yBAAyB,EAAE,CAAC;QACnF,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,6BAA6B,EAAE,CAAC;QACnF,OAAO,IAAI,CAAC;IACd,CAAC;IAED,IAAI,CAAC;QACH,yBAAyB,CAAC,SAAS,CAAC,QAAQ,EAAE,SAAS,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;IACpF,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;IAED,OAAO,kBAAkB,CAAC,SAAS,EAAE,mBAAmB,CAAC,CAAC;AAC5D,CAAC;AAED,MAAM,UAAU,iBAAiB,CAAC,CAAS,EAAE,CAAS;IACpD,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;IACnB,MAAM,KAAK,GAAG,oCAAoC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACzD,MAAM,UAAU,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACvD,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;QACxB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,MAAM,UAAU,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACvD,IAAI,UAAU,KAAK,IAAI,EAAE,CAAC;QACxB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,IAAI,mBAAmB,GAAG,KAAK,CAAC;IAChC,IAAI,gBAAgB,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAEpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,2BAA2B,EAAE,CAAC,EAAE,EAAE,CAAC;QACxD,MAAM,WAAW,GAAG,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,2BAA2B,CAAC,CAAC;QAClF,MAAM,QAAQ,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QAEpD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC;YACrE,mBAAmB,GAAG,WAAW,CAAC;YAClC,gBAAgB,GAAG,QAAQ,CAAC;YAC5B,SAAS;QACX,CAAC;QAED,IAAI,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC,EAAE,CAAC;YACjC,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC,CAAC;QAC3D,CAAC;QAED,IAAI,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC;YACzB,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QACnD,CAAC;QAED,IAAI,gBAAgB,GAAG,QAAQ,GAAG,CAAC,EAAE,CAAC;YACpC,MAAM,IAAI,GAAG,eAAe,CAC1B,CAAC,oBAAoB,EAAE,EAAE,CAAC,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,oBAAoB,CAAC,EACrE,mBAAmB,EACnB,WAAW,EACX,CAAC,mBAAmB,GAAG,WAAW,CAAC,GAAG,CAAC,EACvC,EAAE,SAAS,EAAE,CAAC,CAAC,6BAA6B,EAAE,CAC/C,CAAC;YACF,OAAO,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QAC5C,CAAC;QAED,mBAAmB,GAAG,WAAW,CAAC;QAClC,gBAAgB,GAAG,QAAQ,CAAC;IAC9B,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/dist/saturation/two-phase.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"two-phase.d.ts","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAQnD;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAS7D;AAID;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAU7D"}
+{"version":3,"file":"two-phase.d.ts","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAYnD;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAS7D;AAID;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAQ7D"}

package/dist/saturation/two-phase.js

@@ -7,9 +7,9 @@
  * Reference: IAPWS-IF97, Region 4 (Saturation Line)
  */
 import { OutOfRangeError } from '../types.js';
-import { Pc, Tc, Pt, T_MIN } from '../constants.js';
 import { saturationPressure } from '../regions/region4.js';
 import { mixSaturationState, saturationEndpointsAtPressure } from './common.js';
+import { assertRegion4PressureStateInput, assertRegion4TemperatureStateInput, normalizeRegion4Pressure, } from './region4-boundaries.js';
 // ─── solvePx ────────────────────────────────────────────────────────────────
 /**
  * Solve for thermodynamic state given pressure and vapour quality.
@@ -20,10 +20,7 @@ export function solvePx(p, x) {
     if (x < 0 || x > 1) {
         throw new OutOfRangeError('Quality', x, 0, 1);
     }
-    if (p < Pt || p > Pc) {
-        throw new OutOfRangeError('Pressure', p, Pt, Pc);
-    }
-    return mixSaturationState(saturationEndpointsAtPressure(p), x);
+    return mixSaturationState(saturationEndpointsAtPressure(assertRegion4PressureStateInput(p, 'solvePx')), x);
 }
 // ─── solveTx ────────────────────────────────────────────────────────────────
 /**
@@ -35,10 +32,8 @@ export function solveTx(T, x) {
     if (x < 0 || x > 1) {
         throw new OutOfRangeError('Quality', x, 0, 1);
     }
-    if (T < T_MIN || T > Tc) {
-        throw new OutOfRangeError('Temperature', T, T_MIN, Tc);
-    }
-    const p = saturationPressure(T);
+    const temperature = assertRegion4TemperatureStateInput(T, 'solveTx');
+    const p = normalizeRegion4Pressure(saturationPressure(temperature));
     return solvePx(p, x);
 }
 //# sourceMappingURL=two-phase.js.map

package/dist/saturation/two-phase.js.map

@@ -1 +1 @@
-{"version":3,"file":"two-phase.js","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,MAAM,iBAAiB,CAAC;AACpD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,6BAA6B,EAAE,MAAM,aAAa,CAAC;AAEhF,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IACD,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACrB,MAAM,IAAI,eAAe,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACnD,CAAC;IAED,OAAO,kBAAkB,CAAC,6BAA6B,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACjE,CAAC;AAED,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IACD,IAAI,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACxB,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;IACzD,CAAC;IAED,MAAM,CAAC,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IAChC,OAAO,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACvB,CAAC"}
+{"version":3,"file":"two-phase.js","sourceRoot":"","sources":["../../src/saturation/two-phase.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,6BAA6B,EAAE,MAAM,aAAa,CAAC;AAChF,OAAO,EACL,+BAA+B,EAC/B,kCAAkC,EAClC,wBAAwB,GACzB,MAAM,yBAAyB,CAAC;AAEjC,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IAED,OAAO,kBAAkB,CACvB,6BAA6B,CAAC,+BAA+B,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,EAC5E,CAAC,CACF,CAAC;AACJ,CAAC;AAED,+EAA+E;AAE/E;;;;GAIG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC;IAED,MAAM,WAAW,GAAG,kCAAkC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACrE,MAAM,CAAC,GAAG,wBAAwB,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC;IACpE,OAAO,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACvB,CAAC"}

package/dist/saturation.d.ts

@@ -1,5 +1,10 @@
 /**
- * Public saturation helper exports.
+ * Public low-level saturation helper exports.
+ *
+ * These functions expose the mathematical Region 4 boundary relations.
+ * They normalize values near Pt/Pc/Tt/Tc, but they are not state constructors:
+ * high-level solvers such as solvePx/solveTx apply stricter triple-point and
+ * critical-point rules.
  */
 export { saturationPressure, saturationTemperature } from './regions/region4.js';
 //# sourceMappingURL=saturation.d.ts.map

package/dist/saturation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"saturation.d.ts","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"saturation.d.ts","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/saturation.js

@@ -1,5 +1,10 @@
 /**
- * Public saturation helper exports.
+ * Public low-level saturation helper exports.
+ *
+ * These functions expose the mathematical Region 4 boundary relations.
+ * They normalize values near Pt/Pc/Tt/Tc, but they are not state constructors:
+ * high-level solvers such as solvePx/solveTx apply stricter triple-point and
+ * critical-point rules.
  */
 export { saturationPressure, saturationTemperature } from './regions/region4.js';
 //# sourceMappingURL=saturation.js.map

package/dist/saturation.js.map

@@ -1 +1 @@
-{"version":3,"file":"saturation.js","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"saturation.js","sourceRoot":"","sources":["../src/saturation.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/types.d.ts

@@ -12,7 +12,7 @@ export declare enum Region {
     Region2 = 2,
     /** Supercritical fluid */
     Region3 = 3,
-    /** Two-phase (wet steam / saturation line) */
+    /** Two-phase (wet steam / saturation line, excluding the exact critical point) */
     Region4 = 4,
     /** High-temperature steam (T > 1073.15 K) */
     Region5 = 5
@@ -44,16 +44,16 @@ export interface SteamState {
     cv: number | null;
     /** Speed of sound [m/s] */
     speedOfSound: number | null;
-    /** Vapour quality (0–1), only defined in Region 4 */
+    /** Vapour quality (0–1), only defined on accepted Region 4 saturation states */
     quality: number | null;
-    /** Dynamic viscosity [Pa·s] */
-    viscosity: number;
-    /** Thermal conductivity [W/(m·K)] */
-    thermalConductivity: number;
+    /** Dynamic viscosity [Pa·s], null for two-phase mixtures */
+    viscosity: number | null;
+    /** Thermal conductivity [W/(m·K)], null for two-phase mixtures */
+    thermalConductivity: number | null;
     /** Surface tension [N/m], only meaningful at saturation */
     surfaceTension: number | null;
-    /** Static dielectric constant [dimensionless] */
-    dielectricConstant: number;
+    /** Static dielectric constant [dimensionless], null for two-phase mixtures */
+    dielectricConstant: number | null;
     /** Ionization constant pKw [dimensionless], null outside released validity range */
     ionizationConstant: number | null;
     /** Isobaric cubic expansion coefficient [1/K] */

package/dist/types.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"types.d.ts","sourceRoot":"","sources":["../src/types.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAIH,uCAAuC;AACvC,oBAAY,MAAM;IAChB,0CAA0C;IAC1C,OAAO,IAAI;IACX,wBAAwB;IACxB,OAAO,IAAI;IACX,0BAA0B;IAC1B,OAAO,IAAI;IACX,8CAA8C;IAC9C,OAAO,IAAI;IACX,6CAA6C;IAC7C,OAAO,IAAI;CACZ;AAID;;;GAGG;AACH,MAAM,WAAW,UAAU;IACzB,8BAA8B;IAC9B,MAAM,EAAE,MAAM,CAAC;IACf,qBAAqB;IACrB,QAAQ,EAAE,MAAM,CAAC;IACjB,sBAAsB;IACtB,WAAW,EAAE,MAAM,CAAC;IACpB,8BAA8B;IAC9B,cAAc,EAAE,MAAM,CAAC;IACvB,sBAAsB;IACtB,OAAO,EAAE,MAAM,CAAC;IAChB,uCAAuC;IACvC,cAAc,EAAE,MAAM,CAAC;IACvB,mCAAmC;IACnC,OAAO,EAAE,MAAM,CAAC;IAChB,gCAAgC;IAChC,QAAQ,EAAE,MAAM,CAAC;IACjB,yCAAyC;IACzC,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,0CAA0C;IAC1C,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,2BAA2B;IAC3B,YAAY,EAAE,MAAM,GAAG,IAAI,CAAC;IAC5B,qDAAqD;IACrD,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IAIvB,+BAA+B;IAC/B,SAAS,EAAE,MAAM,CAAC;IAClB,qCAAqC;IACrC,mBAAmB,EAAE,MAAM,CAAC;IAC5B,2DAA2D;IAC3D,cAAc,EAAE,MAAM,GAAG,IAAI,CAAC;IAC9B,iDAAiD;IACjD,kBAAkB,EAAE,MAAM,CAAC;IAC3B,oFAAoF;IACpF,kBAAkB,EAAE,MAAM,GAAG,IAAI,CAAC;IAIlC,iDAAiD;IACjD,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,yCAAyC;IACzC,yBAAyB,EAAE,MAAM,GAAG,IAAI,CAAC;CAC1C;AAID;;;GAGG;AACH,MAAM,WAAW,eAAe;IAC9B,MAAM,EAAE,MAAM,CAAC;IACf,QAAQ,EAAE,MAAM,CAAC;IACjB,WAAW,EAAE,MAAM,CAAC;IACpB,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IACvB,OAAO,EAAE,MAAM,CAAC;IAChB,QAAQ,EAAE,MAAM,CAAC;IACjB,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,YAAY,EAAE,MAAM,GAAG,IAAI,CAAC;IAC5B,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,yBAAyB,EAAE,MAAM,GAAG,IAAI,CAAC;CAC1C;AAID,iDAAiD;AACjD,MAAM,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;AAG9E,KAAK,gBAAgB,CAAC,KAAK,SAAS,MAAM,EAAE,IAAI,SAAS,MAAM,IAC3D,CAAC;KAAG,CAAC,IAAI,KAAK,GAAG,MAAM;CAAE,GAAG;KAAG,CAAC,IAAI,IAAI,CAAC,CAAC,EAAE,MAAM;CAAE,CAAC,GACrD,CAAC;KAAG,CAAC,IAAI,KAAK,CAAC,CAAC,EAAE,MAAM;CAAE,GAAG;KAAG,CAAC,IAAI,IAAI,GAAG,MAAM;CAAE,CAAC,CAAC;AAE1D,iDAAiD;AACjD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,CAAC;AAE7E,8CAA8C;AAC9C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,CAAC;AAE1E,6CAA6C;AAC7C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,6CAA6C;AAC7C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,0DAA0D;AAC1D,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,6DAA6D;AAC7D,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAE5E,iDAAiD;AACjD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,CAAC;AAE7E,gDAAgD;AAChD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAE5E,gCAAgC;AAChC,MAAM,MAAM,UAAU,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;AAIvG;;;;;;GAMG;AACH,MAAM,MAAM,gBAAgB,GAAG,SAAS,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,EAAE,CAAC;AAI9E,uDAAuD;AACvD,qBAAa,SAAU,SAAQ,KAAK;gBACtB,OAAO,EAAE,MAAM;CAI5B;AAED,mEAAmE;AACnE,qBAAa,eAAgB,SAAQ,SAAS;aAE1B,SAAS,EAAE,MAAM;aACjB,KAAK,EAAE,MAAM;aACb,GAAG,EAAE,MAAM;aACX,GAAG,EAAE,MAAM;gBAHX,SAAS,EAAE,MAAM,EACjB,KAAK,EAAE,MAAM,EACb,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM;CAO9B;AAED,uDAAuD;AACvD,qBAAa,gBAAiB,SAAQ,SAAS;aAE3B,MAAM,EAAE,MAAM;aACd,UAAU,EAAE,MAAM;gBADlB,MAAM,EAAE,MAAM,EACd,UAAU,EAAE,MAAM;CAKrC"}
+{"version":3,"file":"types.d.ts","sourceRoot":"","sources":["../src/types.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAIH,uCAAuC;AACvC,oBAAY,MAAM;IAChB,0CAA0C;IAC1C,OAAO,IAAI;IACX,wBAAwB;IACxB,OAAO,IAAI;IACX,0BAA0B;IAC1B,OAAO,IAAI;IACX,kFAAkF;IAClF,OAAO,IAAI;IACX,6CAA6C;IAC7C,OAAO,IAAI;CACZ;AAID;;;GAGG;AACH,MAAM,WAAW,UAAU;IACzB,8BAA8B;IAC9B,MAAM,EAAE,MAAM,CAAC;IACf,qBAAqB;IACrB,QAAQ,EAAE,MAAM,CAAC;IACjB,sBAAsB;IACtB,WAAW,EAAE,MAAM,CAAC;IACpB,8BAA8B;IAC9B,cAAc,EAAE,MAAM,CAAC;IACvB,sBAAsB;IACtB,OAAO,EAAE,MAAM,CAAC;IAChB,uCAAuC;IACvC,cAAc,EAAE,MAAM,CAAC;IACvB,mCAAmC;IACnC,OAAO,EAAE,MAAM,CAAC;IAChB,gCAAgC;IAChC,QAAQ,EAAE,MAAM,CAAC;IACjB,yCAAyC;IACzC,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,0CAA0C;IAC1C,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,2BAA2B;IAC3B,YAAY,EAAE,MAAM,GAAG,IAAI,CAAC;IAC5B,gFAAgF;IAChF,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IAIvB,4DAA4D;IAC5D,SAAS,EAAE,MAAM,GAAG,IAAI,CAAC;IACzB,kEAAkE;IAClE,mBAAmB,EAAE,MAAM,GAAG,IAAI,CAAC;IACnC,2DAA2D;IAC3D,cAAc,EAAE,MAAM,GAAG,IAAI,CAAC;IAC9B,8EAA8E;IAC9E,kBAAkB,EAAE,MAAM,GAAG,IAAI,CAAC;IAClC,oFAAoF;IACpF,kBAAkB,EAAE,MAAM,GAAG,IAAI,CAAC;IAIlC,iDAAiD;IACjD,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,yCAAyC;IACzC,yBAAyB,EAAE,MAAM,GAAG,IAAI,CAAC;CAC1C;AAID;;;GAGG;AACH,MAAM,WAAW,eAAe;IAC9B,MAAM,EAAE,MAAM,CAAC;IACf,QAAQ,EAAE,MAAM,CAAC;IACjB,WAAW,EAAE,MAAM,CAAC;IACpB,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IACvB,OAAO,EAAE,MAAM,CAAC;IAChB,QAAQ,EAAE,MAAM,CAAC;IACjB,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,EAAE,EAAE,MAAM,GAAG,IAAI,CAAC;IAClB,YAAY,EAAE,MAAM,GAAG,IAAI,CAAC;IAC5B,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,yBAAyB,EAAE,MAAM,GAAG,IAAI,CAAC;CAC1C;AAID,iDAAiD;AACjD,MAAM,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;AAG9E,KAAK,gBAAgB,CAAC,KAAK,SAAS,MAAM,EAAE,IAAI,SAAS,MAAM,IAC3D,CAAC;KAAG,CAAC,IAAI,KAAK,GAAG,MAAM;CAAE,GAAG;KAAG,CAAC,IAAI,IAAI,CAAC,CAAC,EAAE,MAAM;CAAE,CAAC,GACrD,CAAC;KAAG,CAAC,IAAI,KAAK,CAAC,CAAC,EAAE,MAAM;CAAE,GAAG;KAAG,CAAC,IAAI,IAAI,GAAG,MAAM;CAAE,CAAC,CAAC;AAE1D,iDAAiD;AACjD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,CAAC;AAE7E,8CAA8C;AAC9C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,CAAC;AAE1E,6CAA6C;AAC7C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,6CAA6C;AAC7C,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,0DAA0D;AAC1D,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAEzE,6DAA6D;AAC7D,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAE5E,iDAAiD;AACjD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,UAAU,CAAC,CAAC;AAE7E,gDAAgD;AAChD,MAAM,MAAM,OAAO,GAAG;IACpB,IAAI,EAAE,IAAI,CAAC;CACZ,GAAG,gBAAgB,CAAC,GAAG,EAAE,aAAa,CAAC,GAAG,gBAAgB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;AAE5E,gCAAgC;AAChC,MAAM,MAAM,UAAU,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;AAIvG;;;;;;GAMG;AACH,MAAM,MAAM,gBAAgB,GAAG,SAAS,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,EAAE,CAAC;AAI9E,uDAAuD;AACvD,qBAAa,SAAU,SAAQ,KAAK;gBACtB,OAAO,EAAE,MAAM;CAI5B;AAED,mEAAmE;AACnE,qBAAa,eAAgB,SAAQ,SAAS;aAE1B,SAAS,EAAE,MAAM;aACjB,KAAK,EAAE,MAAM;aACb,GAAG,EAAE,MAAM;aACX,GAAG,EAAE,MAAM;gBAHX,SAAS,EAAE,MAAM,EACjB,KAAK,EAAE,MAAM,EACb,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM;CAO9B;AAED,uDAAuD;AACvD,qBAAa,gBAAiB,SAAQ,SAAS;aAE3B,MAAM,EAAE,MAAM;aACd,UAAU,EAAE,MAAM;gBADlB,MAAM,EAAE,MAAM,EACd,UAAU,EAAE,MAAM;CAKrC"}

package/dist/types.js

@@ -14,7 +14,7 @@ export var Region;
     Region[Region["Region2"] = 2] = "Region2";
     /** Supercritical fluid */
     Region[Region["Region3"] = 3] = "Region3";
-    /** Two-phase (wet steam / saturation line) */
+    /** Two-phase (wet steam / saturation line, excluding the exact critical point) */
     Region[Region["Region4"] = 4] = "Region4";
     /** High-temperature steam (T > 1073.15 K) */
     Region[Region["Region5"] = 5] = "Region5";

package/dist/types.js.map

@@ -1 +1 @@
-{"version":3,"file":"types.js","sourceRoot":"","sources":["../src/types.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,+EAA+E;AAE/E,uCAAuC;AACvC,MAAM,CAAN,IAAY,MAWX;AAXD,WAAY,MAAM;IAChB,0CAA0C;IAC1C,yCAAW,CAAA;IACX,wBAAwB;IACxB,yCAAW,CAAA;IACX,0BAA0B;IAC1B,yCAAW,CAAA;IACX,8CAA8C;IAC9C,yCAAW,CAAA;IACX,6CAA6C;IAC7C,yCAAW,CAAA;AACb,CAAC,EAXW,MAAM,KAAN,MAAM,QAWjB;AA6ID,+EAA+E;AAE/E,uDAAuD;AACvD,MAAM,OAAO,SAAU,SAAQ,KAAK;IAClC,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,WAAW,CAAC;IAC1B,CAAC;CACF;AAED,mEAAmE;AACnE,MAAM,OAAO,eAAgB,SAAQ,SAAS;IAE1B;IACA;IACA;IACA;IAJlB,YACkB,SAAiB,EACjB,KAAa,EACb,GAAW,EACX,GAAW;QAE3B,KAAK,CACH,GAAG,SAAS,MAAM,KAAK,qBAAqB,GAAG,KAAK,GAAG,GAAG,CAC3D,CAAC;QAPc,cAAS,GAAT,SAAS,CAAQ;QACjB,UAAK,GAAL,KAAK,CAAQ;QACb,QAAG,GAAH,GAAG,CAAQ;QACX,QAAG,GAAH,GAAG,CAAQ;QAK3B,IAAI,CAAC,IAAI,GAAG,iBAAiB,CAAC;IAChC,CAAC;CACF;AAED,uDAAuD;AACvD,MAAM,OAAO,gBAAiB,SAAQ,SAAS;IAE3B;IACA;IAFlB,YACkB,MAAc,EACd,UAAkB;QAElC,KAAK,CAAC,GAAG,MAAM,6BAA6B,UAAU,aAAa,CAAC,CAAC;QAHrD,WAAM,GAAN,MAAM,CAAQ;QACd,eAAU,GAAV,UAAU,CAAQ;QAGlC,IAAI,CAAC,IAAI,GAAG,kBAAkB,CAAC;IACjC,CAAC;CACF"}
+{"version":3,"file":"types.js","sourceRoot":"","sources":["../src/types.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,+EAA+E;AAE/E,uCAAuC;AACvC,MAAM,CAAN,IAAY,MAWX;AAXD,WAAY,MAAM;IAChB,0CAA0C;IAC1C,yCAAW,CAAA;IACX,wBAAwB;IACxB,yCAAW,CAAA;IACX,0BAA0B;IAC1B,yCAAW,CAAA;IACX,kFAAkF;IAClF,yCAAW,CAAA;IACX,6CAA6C;IAC7C,yCAAW,CAAA;AACb,CAAC,EAXW,MAAM,KAAN,MAAM,QAWjB;AA6ID,+EAA+E;AAE/E,uDAAuD;AACvD,MAAM,OAAO,SAAU,SAAQ,KAAK;IAClC,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,WAAW,CAAC;IAC1B,CAAC;CACF;AAED,mEAAmE;AACnE,MAAM,OAAO,eAAgB,SAAQ,SAAS;IAE1B;IACA;IACA;IACA;IAJlB,YACkB,SAAiB,EACjB,KAAa,EACb,GAAW,EACX,GAAW;QAE3B,KAAK,CACH,GAAG,SAAS,MAAM,KAAK,qBAAqB,GAAG,KAAK,GAAG,GAAG,CAC3D,CAAC;QAPc,cAAS,GAAT,SAAS,CAAQ;QACjB,UAAK,GAAL,KAAK,CAAQ;QACb,QAAG,GAAH,GAAG,CAAQ;QACX,QAAG,GAAH,GAAG,CAAQ;QAK3B,IAAI,CAAC,IAAI,GAAG,iBAAiB,CAAC;IAChC,CAAC;CACF;AAED,uDAAuD;AACvD,MAAM,OAAO,gBAAiB,SAAQ,SAAS;IAE3B;IACA;IAFlB,YACkB,MAAc,EACd,UAAkB;QAElC,KAAK,CAAC,GAAG,MAAM,6BAA6B,UAAU,aAAa,CAAC,CAAC;QAHrD,WAAM,GAAN,MAAM,CAAQ;QACd,eAAU,GAAV,UAAU,CAAQ;QAGlC,IAAI,CAAC,IAAI,GAAG,kBAAkB,CAAC;IACjC,CAAC;CACF"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "iapws-if97",
-  "version": "2.0.2",
+  "version": "2.1.1",
   "description": "TypeScript implementation of IAPWS-IF97 for water and steam properties, including forward/backward solvers, saturation states, and transport properties.",
   "author": "jltonghui",
   "repository": {