iapws-if97 2.1.5 → 2.8.0

package/dist/backward/tolerances.js

@@ -1,21 +1,26 @@
 export const BACKWARD_REGION4_QUALITY_TOLERANCE = 1e-6;
+export const BACKWARD_HS_NEWTON_TOLERANCE = 1e-12;
+export const BACKWARD_PRESSURE_PROPERTY_TEMPERATURE_TOLERANCE = 1e-13;
+export function fixedTemperatureRootTolerance(property) {
+    return property === 'enthalpy' ? 1e-9 : 1e-12;
+}
 export function backwardConstraintTolerance(label, expected) {
     switch (label) {
         case 'pressure':
+            return 1e-9 * Math.abs(expected);
         case 'temperature':
             return 1e-9 * Math.max(1, Math.abs(expected));
         case 'enthalpy':
-            return 1e-6 * Math.max(1, Math.abs(expected));
+            return 1e-9 * Math.max(1, Math.abs(expected));
         case 'entropy':
-            return 1e-6 * Math.max(1, Math.abs(expected));
+            return 1e-9 * Math.max(1, Math.abs(expected));
         default:
             return 1e-9;
     }
 }
 export function region4SaturationPressureTolerance(pressure) {
-    return 1e-8 * Math.max(1, Math.abs(pressure));
+    return 1e-8 * Math.abs(pressure);
 }
 export function backwardSpecificVolumeTolerance(expected) {
     return 1e-6 * Math.max(1, Math.abs(expected));
 }
-//# sourceMappingURL=tolerances.js.map

package/dist/backward/ts.d.ts

@@ -1,16 +1,13 @@
 /**
  * Backward equation T,S → state.
  *
- * Given temperature and specific entropy, determines the IAPWS-IF97 region
- * and solves for pressure using a bracketed Newton method, then returns
- * the full set of basic thermodynamic properties.
+ * Returns a state only when temperature and entropy determine it uniquely.
  */
 import type { BasicProperties } from '../types.js';
 /**
  * Solve for full basic thermodynamic state from temperature and entropy.
  *
- * Detects the IAPWS-IF97 region, then inverts s(p, T) to recover pressure
- * within the appropriate bracket for that region.
+ * Searches every valid isotherm branch and rejects ambiguous inputs.
  *
  * Temperatures within ±CRITICAL_T_EXCLUSION_BAND of Tc are rejected because
  * the entropy surface is nearly flat there, making pressure recovery unreliable.
@@ -20,6 +17,6 @@ import type { BasicProperties } from '../types.js';
  * @returns Basic thermodynamic properties (without transport)
  * @throws {OutOfRangeError} if T is outside [T_MIN, T_MAX]
  * @throws {IF97Error} if T is within the critical exclusion band
+ * @throws {AmbiguousStateError} if multiple physical states match T and s
  */
 export declare function solveTS(T: number, s: number): BasicProperties;
-//# sourceMappingURL=ts.d.ts.map

package/dist/backward/ts.js

@@ -1,23 +1,11 @@
 import * as C from '../constants.js';
-import { IF97Error, OutOfRangeError, Region } from '../types.js';
-import { assertFiniteNumber } from '../core/input-validation.js';
-import { detectRegionTS } from '../core/region-detector.js';
-import { solveRegion3PTBasic } from '../core/region3-pt.js';
-import { region1 } from '../regions/region1.js';
-import { region2 } from '../regions/region2.js';
-import { region5 } from '../regions/region5.js';
-import { boundary23_T_to_P } from '../regions/boundaries.js';
-import { saturationPressure } from '../regions/region4.js';
-import { mixSaturationState, qualityFromSaturationProperty, saturationEndpointsAtTemperature, } from '../saturation/common.js';
-import { assertRegion4StateAllowed } from '../saturation/region4-boundaries.js';
-import { solveFixedTemperaturePressure } from './fixed-temperature-solver.js';
-import { validateBackwardState } from './solution-validation.js';
-const extractS = (s) => s.entropy;
+import { IF97Error, OutOfRangeError } from '../types.js';
+import { assertFiniteNumber, withinClosedLowerBound } from '../core/input-validation.js';
+import { uniqueFixedTemperatureState } from './fixed-temperature-candidates.js';
 /**
  * Solve for full basic thermodynamic state from temperature and entropy.
  *
- * Detects the IAPWS-IF97 region, then inverts s(p, T) to recover pressure
- * within the appropriate bracket for that region.
+ * Searches every valid isotherm branch and rejects ambiguous inputs.
  *
  * Temperatures within ±CRITICAL_T_EXCLUSION_BAND of Tc are rejected because
  * the entropy surface is nearly flat there, making pressure recovery unreliable.
@@ -27,6 +15,7 @@ const extractS = (s) => s.entropy;
  * @returns Basic thermodynamic properties (without transport)
  * @throws {OutOfRangeError} if T is outside [T_MIN, T_MAX]
  * @throws {IF97Error} if T is within the critical exclusion band
+ * @throws {AmbiguousStateError} if multiple physical states match T and s
  */
 export function solveTS(T, s) {
     assertFiniteNumber('Temperature', T);
@@ -34,67 +23,15 @@ export function solveTS(T, s) {
     if (T < C.T_MIN || T > C.T_MAX) {
         throw new OutOfRangeError('Temperature', T, C.T_MIN, C.T_MAX);
     }
-    if (s < C.S_MIN || s > C.S_MAX) {
+    if (!withinClosedLowerBound(s, C.S_MIN, 'entropy') || s > C.S_MAX) {
         throw new OutOfRangeError('Entropy', s, C.S_MIN, C.S_MAX);
     }
     if (Math.abs(T - C.Tc) < C.CRITICAL_T_EXCLUSION_BAND) {
         throw new IF97Error(`solveTS does not support temperatures within ${C.CRITICAL_T_EXCLUSION_BAND} K of the critical-point temperature Tc=${C.Tc} K; received T=${T} K`);
     }
-    const region = detectRegionTS(T, s);
-    switch (region) {
-        case Region.Region1:
-            return validateBackwardState(solveFixedTemperaturePressure(region1, T, s, saturationPressure(T), C.P_MAX, extractS, 'solveTS'), [
-                { label: 'temperature', expected: T },
-                { label: 'entropy', expected: s },
-            ], {
-                solverName: 'solveTS',
-                expectedRegion: Region.Region1,
-            });
-        case Region.Region2: {
-            const upper = T <= C.Tc
-                ? saturationPressure(T)
-                : T <= C.B23_T_MAX
-                    ? boundary23_T_to_P(T)
-                    : C.P_MAX;
-            return validateBackwardState(solveFixedTemperaturePressure(region2, T, s, C.P_MIN, upper, extractS, 'solveTS'), [
-                { label: 'temperature', expected: T },
-                { label: 'entropy', expected: s },
-            ], {
-                solverName: 'solveTS',
-                expectedRegion: Region.Region2,
-            });
-        }
-        case Region.Region3: {
-            const lower = T <= C.Tc
-                ? saturationPressure(T)
-                : boundary23_T_to_P(T);
-            return validateBackwardState(solveFixedTemperaturePressure(solveRegion3PTBasic, T, s, lower, C.P_MAX, extractS, 'solveTS'), [
-                { label: 'temperature', expected: T },
-                { label: 'entropy', expected: s },
-            ], {
-                solverName: 'solveTS',
-                expectedRegion: Region.Region3,
-            });
-        }
-        case Region.Region4: {
-            const endpoints = saturationEndpointsAtTemperature(T);
-            assertRegion4StateAllowed(endpoints.pressure, endpoints.temperature, 'solveTS');
-            const x = qualityFromSaturationProperty(endpoints.liquid.entropy, endpoints.vapor.entropy, s);
-            return validateBackwardState(mixSaturationState(endpoints, x), [
-                { label: 'temperature', expected: T },
-                { label: 'entropy', expected: s },
-            ], { solverName: 'solveTS', expectedRegion: Region.Region4 });
-        }
-        case Region.Region5:
-            return validateBackwardState(solveFixedTemperaturePressure(region5, T, s, C.P_MIN, C.R5_P_MAX, extractS, 'solveTS'), [
-                { label: 'temperature', expected: T },
-                { label: 'entropy', expected: s },
-            ], {
-                solverName: 'solveTS',
-                expectedRegion: Region.Region5,
-            });
-        default:
-            throw new IF97Error(`solveTS cannot determine a valid IF97 state for T=${T} K, s=${s} kJ/(kg·K)`);
+    const state = uniqueFixedTemperatureState(T, s, 'entropy');
+    if (state === null) {
+        throw new IF97Error(`solveTS cannot determine a valid IF97 state for T=${T} K, s=${s} kJ/(kg·K)`);
     }
+    return state;
 }
-//# sourceMappingURL=ts.js.map

package/dist/boundaries.d.ts

@@ -3,4 +3,3 @@
  */
 export { boundary23_T_to_P, boundary23_P_to_T } from './regions/boundaries.js';
 export { region3Volume, region3SatVolume } from './regions/region3-subregions.js';
-//# sourceMappingURL=boundaries.d.ts.map

package/dist/boundaries.js

@@ -3,4 +3,3 @@
  */
 export { boundary23_T_to_P, boundary23_P_to_T } from './regions/boundaries.js';
 export { region3Volume, region3SatVolume } from './regions/region3-subregions.js';
-//# sourceMappingURL=boundaries.js.map

package/dist/constants.d.ts

@@ -4,11 +4,7 @@
  * Physical constants, critical point properties, and range limits
  * for the IAPWS Industrial Formulation 1997.
  */
-/** Specific gas constant for water [kJ/(kg·K)].
- *  IAPWS-IF97, §2 (Table 1). Fixed value for IF97 internal consistency.
- *  Do not replace transport correlations with this value; thermal conductivity
- *  critical enhancement intentionally uses `R_IAPWS_2011_THERMAL`. */
-export declare const R = 0.461526;
+export { R } from './gas-constant.js';
 /** Specific gas constant used by the IAPWS 2011 thermal conductivity λ₂ correlation [kJ/(kg·K)].
  *  Intentionally differs from IF97's `R`; keep both constants distinct and named. */
 export declare const R_IAPWS_2011_THERMAL = 0.46151805;
@@ -16,20 +12,18 @@ export declare const R_IAPWS_2011_THERMAL = 0.46151805;
 export declare const Tc = 647.096;
 /** Temperature exclusion half-band around the critical point [K] */
 export declare const CRITICAL_T_EXCLUSION_BAND = 0.001;
-/** Region 4 endpoint snap tolerance for temperatures [K] */
+/** Region 4 tolerance for outward temperature-boundary roundoff [K] */
 export declare const REGION4_TEMPERATURE_TOLERANCE = 1e-9;
 /** Safety margin used when clamping Region 4 temperatures below Tc [K] */
 export declare const REGION4_SUBCRITICAL_TEMPERATURE_MARGIN = 1e-8;
 /** Critical pressure [MPa] */
 export declare const Pc = 22.064;
-/** Region 4 endpoint snap tolerance for pressures [MPa] */
+/** Relative tolerance for outward Region 4 pressure-boundary roundoff */
 export declare const REGION4_PRESSURE_TOLERANCE = 1e-9;
 /** Tolerance for identifying the critical h-s endpoint in Region 4 */
 export declare const REGION4_CRITICAL_HS_TOLERANCE = 0.000001;
 /** Residual tolerance for Region 4 h-s temperature inversion */
 export declare const REGION4_HS_RESIDUAL_TOLERANCE = 1e-10;
-/** Scan resolution used when bracketing Region 4 h-s temperature roots */
-export declare const REGION4_HS_BRACKET_SEGMENTS = 256;
 /** Critical density [kg/m³] */
 export declare const RHOc = 322;
 /** Triple point pressure [MPa] */
@@ -46,14 +40,14 @@ export declare const P_MAX = 100;
 export declare const T_MIN = 273.15;
 /** Maximum temperature [K] (2273.15 K = 2000 °C) */
 export declare const T_MAX = 2273.15;
-/** Minimum specific entropy [kJ/(kg·K)] */
-export declare const S_MIN = -0.00015454959194;
-/** Maximum specific entropy [kJ/(kg·K)] */
-export declare const S_MAX = 13.904956083429227;
-/** Minimum specific enthalpy [kJ/kg] */
-export declare const H_MIN = -0.041587825987;
-/** Maximum specific enthalpy [kJ/kg] */
-export declare const H_MAX = 7376.980263598506;
+/** Minimum specific entropy [kJ/(kg·K)], at T_MIN and P_MAX in Region 1. */
+export declare const S_MIN = -0.00858228709261889;
+/** Maximum specific entropy [kJ/(kg·K)], at T_MAX and P_MIN in Region 5. */
+export declare const S_MAX: number;
+/** Minimum specific enthalpy [kJ/kg], at T_MIN and P_MIN in Region 1. */
+export declare const H_MIN = -0.041587825987655716;
+/** Maximum specific enthalpy [kJ/kg], at T_MAX and P_MIN in Region 5. */
+export declare const H_MAX: number;
 /** Region 2 minimum temperature for high-P subregions [K] (623.15 K = 350 °C) */
 export declare const R2_T_MIN = 623.15;
 /** Region 2 maximum temperature for backward eqs [K] (1073.15 K = 800 °C) */
@@ -70,8 +64,6 @@ export declare const R2_P_CRT = 4;
 export declare const B23_P_MIN = 16.5291642526;
 /** B23 maximum temperature [K] */
 export declare const B23_T_MAX = 863.15;
-/** Region 3 minimum temperature [K] */
-export declare const R3_T_MIN = 623.15;
 /** Region 3 critical entropy [kJ/(kg·K)] */
 export declare const R3_S_CRT = 4.41202148223476;
 /** Region 3 critical enthalpy [kJ/kg] */
@@ -86,8 +78,3 @@ export declare const R5_T_MAX = 2273.15;
 export declare const B23_S_CURVE_MIN = 5.048096828;
 /** Maximum entropy along the B23 curve [kJ/(kg·K)], at T ≈ 644 K */
 export declare const B23_S_CURVE_MAX = 5.260578707;
-/** Minimum enthalpy along the B23 curve [kJ/kg], at T = 623.15 K (currently unused) */
-export declare const B23_H_CURVE_MIN = 2563.592004;
-/** Maximum enthalpy along the B23 curve [kJ/kg], at T = 863.15 K (currently unused) */
-export declare const B23_H_CURVE_MAX = 2812.942061;
-//# sourceMappingURL=constants.d.ts.map

package/dist/constants.js

@@ -4,12 +4,9 @@
  * Physical constants, critical point properties, and range limits
  * for the IAPWS Industrial Formulation 1997.
  */
+import { region5 } from './regions/region5.js';
 // ─── Fundamental Constants ──────────────────────────────────────────────────
-/** Specific gas constant for water [kJ/(kg·K)].
- *  IAPWS-IF97, §2 (Table 1). Fixed value for IF97 internal consistency.
- *  Do not replace transport correlations with this value; thermal conductivity
- *  critical enhancement intentionally uses `R_IAPWS_2011_THERMAL`. */
-export const R = 0.461526;
+export { R } from './gas-constant.js';
 /** Specific gas constant used by the IAPWS 2011 thermal conductivity λ₂ correlation [kJ/(kg·K)].
  *  Intentionally differs from IF97's `R`; keep both constants distinct and named. */
 export const R_IAPWS_2011_THERMAL = 0.46151805;
@@ -18,20 +15,18 @@ export const R_IAPWS_2011_THERMAL = 0.46151805;
 export const Tc = 647.096;
 /** Temperature exclusion half-band around the critical point [K] */
 export const CRITICAL_T_EXCLUSION_BAND = 1e-3;
-/** Region 4 endpoint snap tolerance for temperatures [K] */
+/** Region 4 tolerance for outward temperature-boundary roundoff [K] */
 export const REGION4_TEMPERATURE_TOLERANCE = 1e-9;
 /** Safety margin used when clamping Region 4 temperatures below Tc [K] */
 export const REGION4_SUBCRITICAL_TEMPERATURE_MARGIN = 1e-8;
 /** Critical pressure [MPa] */
 export const Pc = 22.064;
-/** Region 4 endpoint snap tolerance for pressures [MPa] */
+/** Relative tolerance for outward Region 4 pressure-boundary roundoff */
 export const REGION4_PRESSURE_TOLERANCE = 1e-9;
 /** Tolerance for identifying the critical h-s endpoint in Region 4 */
 export const REGION4_CRITICAL_HS_TOLERANCE = 1e-6;
 /** Residual tolerance for Region 4 h-s temperature inversion */
 export const REGION4_HS_RESIDUAL_TOLERANCE = 1e-10;
-/** Scan resolution used when bracketing Region 4 h-s temperature roots */
-export const REGION4_HS_BRACKET_SEGMENTS = 256;
 /** Critical density [kg/m³] */
 export const RHOc = 322;
 // ─── Triple Point ───────────────────────────────────────────────────────────
@@ -50,14 +45,16 @@ export const P_MAX = 100.0;
 export const T_MIN = 273.15;
 /** Maximum temperature [K] (2273.15 K = 2000 °C) */
 export const T_MAX = 2273.15;
-/** Minimum specific entropy [kJ/(kg·K)] */
-export const S_MIN = -0.00015454959194;
-/** Maximum specific entropy [kJ/(kg·K)] */
-export const S_MAX = 13.904956083429227;
-/** Minimum specific enthalpy [kJ/kg] */
-export const H_MIN = -0.041587825987;
-/** Maximum specific enthalpy [kJ/kg] */
-export const H_MAX = 7376.980263598506;
+// Match the Region 5 equation's rounding in the current runtime.
+const REGION5_UPPER = region5(P_MIN, T_MAX);
+/** Minimum specific entropy [kJ/(kg·K)], at T_MIN and P_MAX in Region 1. */
+export const S_MIN = -0.00858228709261889;
+/** Maximum specific entropy [kJ/(kg·K)], at T_MAX and P_MIN in Region 5. */
+export const S_MAX = REGION5_UPPER.entropy;
+/** Minimum specific enthalpy [kJ/kg], at T_MIN and P_MIN in Region 1. */
+export const H_MIN = -0.041587825987655716;
+/** Maximum specific enthalpy [kJ/kg], at T_MAX and P_MIN in Region 5. */
+export const H_MAX = REGION5_UPPER.enthalpy;
 // ─── Region 2 Boundaries ───────────────────────────────────────────────────
 /** Region 2 minimum temperature for high-P subregions [K] (623.15 K = 350 °C) */
 export const R2_T_MIN = 623.15;
@@ -77,8 +74,6 @@ export const B23_P_MIN = 16.5291642526;
 /** B23 maximum temperature [K] */
 export const B23_T_MAX = 863.15;
 // ─── Region 3 Boundaries ───────────────────────────────────────────────────
-/** Region 3 minimum temperature [K] */
-export const R3_T_MIN = 623.15;
 /** Region 3 critical entropy [kJ/(kg·K)] */
 export const R3_S_CRT = 4.41202148223476;
 /** Region 3 critical enthalpy [kJ/kg] */
@@ -90,12 +85,11 @@ export const R5_T_MIN = 1073.15;
 export const R5_P_MAX = 50;
 /** Region 5 maximum temperature [K] */
 export const R5_T_MAX = 2273.15;
-// ─── B23 Curve Extrema (Entropy / Enthalpy) ────────────────────────────────
+// ─── B23 Curve Extrema (Entropy) ───────────────────────────────────────────
 //
 // These are the GLOBAL EXTREMA along the B23 boundary curve, NOT values at
-// the curve endpoints. Entropy has a minimum at T ≈ 777 K and a maximum at
-// T ≈ 644 K. Enthalpy is monotonic: min at the low-T endpoint (623.15 K),
-// max at the high-T endpoint (863.15 K).
+// the curve endpoints: entropy has a minimum at T ≈ 777 K and a maximum at
+// T ≈ 644 K.
 //
 // Ref: Computed from Region 2 properties evaluated along B23(T) for
 //      T ∈ [623.15, 863.15] K.
@@ -103,8 +97,3 @@ export const R5_T_MAX = 2273.15;
 export const B23_S_CURVE_MIN = 5.048096828;
 /** Maximum entropy along the B23 curve [kJ/(kg·K)], at T ≈ 644 K */
 export const B23_S_CURVE_MAX = 5.260578707;
-/** Minimum enthalpy along the B23 curve [kJ/kg], at T = 623.15 K (currently unused) */
-export const B23_H_CURVE_MIN = 2563.592004;
-/** Maximum enthalpy along the B23 curve [kJ/kg], at T = 863.15 K (currently unused) */
-export const B23_H_CURVE_MAX = 2812.942061;
-//# sourceMappingURL=constants.js.map

package/dist/core/floating-point.d.ts

@@ -0,0 +1,2 @@
+/** Adjacent IEEE-754 double; value must be positive and finite. */
+export declare function nextPositiveFloat(value: number, direction: 1 | -1): number;

package/dist/core/floating-point.js

@@ -0,0 +1,7 @@
+const bits = new DataView(new ArrayBuffer(8));
+/** Adjacent IEEE-754 double; value must be positive and finite. */
+export function nextPositiveFloat(value, direction) {
+    bits.setFloat64(0, value);
+    bits.setBigUint64(0, bits.getBigUint64(0) + BigInt(direction));
+    return bits.getFloat64(0);
+}

package/dist/core/formulation-seams.d.ts

@@ -0,0 +1,27 @@
+import type { BasicProperties } from '../types.js';
+export type FormulationSeam = '13' | '23' | '25';
+export interface FormulationSeamPair {
+    seam: FormulationSeam;
+    pressure: number;
+    temperature: number;
+    /** Equations on the lower and higher temperature sides at fixed pressure. */
+    lower: BasicProperties;
+    upper: BasicProperties;
+    canonical: BasicProperties;
+}
+export interface FormulationSeamResolution {
+    pair: FormulationSeamPair;
+    kind: 'canonical' | 'gap' | 'overlap';
+    state: BasicProperties;
+    budgets: {
+        enthalpy: number;
+        entropy: number;
+    };
+}
+export declare function formulationSeamAtPressure(p: number, seam: FormulationSeam): FormulationSeamPair | null;
+/** B13/B25 cannot be crossed while retaining an independently supplied temperature. */
+export declare function formulationSeamAtTemperature(T: number): FormulationSeamPair | null;
+export declare function findPressureSeam(p: number, target: number, property: 'enthalpy' | 'entropy'): FormulationSeamResolution | null;
+export declare function matchFixedSeam(pair: FormulationSeamPair, property: 'enthalpy' | 'entropy', target: number, coordinate?: 'pressure' | 'temperature'): FormulationSeamResolution | null;
+/** Locate both constraints on one bounded ribbon between independently fitted EOS. */
+export declare function findHSSeam(h: number, s: number): FormulationSeamResolution | null;