iapws-if97 1.0.0

package/dist/regions/region3.js

@@ -0,0 +1,112 @@
+/**
+ * IAPWS-IF97 Region 3: Supercritical Fluid
+ *
+ * Basic equation for the specific Helmholtz free energy φ(δ,τ)
+ * Valid for: 623.15 K ≤ T ≤ T_B23(P), P_B23(T) ≤ P ≤ 100 MPa
+ *
+ * Reference: IAPWS-IF97, Section 7 (Equations for Region 3)
+ */
+import { R } from '../constants.js';
+import { Region } from '../types.js';
+// ─── Coefficient Table (Table 30, IAPWS-IF97) ──────────────────────────────
+const COEFFICIENTS = [
+    [0, 0, 1.0658070028513],
+    [0, 0, -0.15732845290239e2],
+    [0, 1, 0.20944396974307e2],
+    [0, 2, -0.76867707878716e1],
+    [0, 7, 0.26185947787954e1],
+    [0, 10, -0.28080781148620e1],
+    [0, 12, 0.12053369696517e1],
+    [0, 23, -0.84566812812502e-2],
+    [1, 2, -0.12654315477714e1],
+    [1, 6, -0.11524407806681e1],
+    [1, 15, 0.88521043984318],
+    [1, 17, -0.64207765181607],
+    [2, 0, 0.38493460186671],
+    [2, 2, -0.85214708824206],
+    [2, 6, 0.48972281541877e1],
+    [2, 7, -0.30502617256965e1],
+    [2, 22, 0.39420536879154e-1],
+    [2, 26, 0.12558408424308],
+    [3, 0, -0.27999329698710],
+    [3, 2, 0.13899799569460e1],
+    [3, 4, -0.20189915023570e1],
+    [3, 16, -0.82147637173963e-2],
+    [3, 26, -0.47596035734923],
+    [4, 0, 0.43984074473500e-1],
+    [4, 2, -0.44476435428739],
+    [4, 4, 0.90572070719733],
+    [4, 26, 0.70522450087967],
+    [5, 1, 0.10770512626332],
+    [5, 3, -0.32913623258954],
+    [5, 26, -0.50871062041158],
+    [6, 0, -0.22175400873096e-1],
+    [6, 2, 0.94260751665092e-1],
+    [6, 26, 0.16436278447961],
+    [7, 2, -0.13503372241348e-1],
+    [8, 26, -0.14834345352472e-1],
+    [9, 2, 0.57922953628084e-3],
+    [9, 26, 0.32308904703711e-2],
+    [10, 0, 0.80964802996215e-4],
+    [10, 1, -0.16557679795037e-3],
+    [11, 26, -0.44923899061815e-4],
+];
+// Reduction constants for Region 3
+const RHO_STAR = 322; // ρ* [kg/m³] (critical density)
+const T_STAR = 647.096; // T* [K] (critical temperature)
+/**
+ * Compute Region 3 thermodynamic properties from density and temperature.
+ *
+ * @param rho - Density [kg/m³]
+ * @param T   - Temperature [K]
+ * @returns Basic thermodynamic properties
+ */
+export function region3ByRhoT(rho, T) {
+    const delta = rho / RHO_STAR;
+    const tau = T_STAR / T;
+    const v = 1 / rho;
+    // First entry is the logarithmic term: n₁ · ln(δ)
+    const n1 = COEFFICIENTS[0][2];
+    let phi = n1 * Math.log(delta);
+    let phi_delta = n1 / delta;
+    let phi_deltadelta = -n1 / (delta * delta);
+    let phi_tau = 0;
+    let phi_tautau = 0;
+    let phi_deltatau = 0;
+    // Sum from i = 1 (skip index 0 which is the ln term)
+    for (let idx = 1; idx < COEFFICIENTS.length; idx++) {
+        const [I, J, N] = COEFFICIENTS[idx];
+        const deltaPow = Math.pow(delta, I);
+        const tauPow = Math.pow(tau, J);
+        phi += N * deltaPow * tauPow;
+        phi_delta += N * I * Math.pow(delta, I - 1) * tauPow;
+        phi_deltadelta += N * I * (I - 1) * Math.pow(delta, I - 2) * tauPow;
+        phi_tau += N * deltaPow * J * Math.pow(tau, J - 1);
+        phi_tautau += N * deltaPow * J * (J - 1) * Math.pow(tau, J - 2);
+        phi_deltatau += N * I * Math.pow(delta, I - 1) * J * Math.pow(tau, J - 1);
+    }
+    return {
+        region: Region.Region3,
+        pressure: rho * R * T * delta * phi_delta / 1000,
+        temperature: T,
+        specificVolume: v,
+        internalEnergy: R * T * tau * phi_tau,
+        entropy: R * (tau * phi_tau - phi),
+        enthalpy: R * T * (tau * phi_tau + delta * phi_delta),
+        cp: R * (-tau * tau * phi_tautau +
+            Math.pow(delta * phi_delta - delta * tau * phi_deltatau, 2) /
+                (2 * delta * phi_delta + delta * delta * phi_deltadelta)),
+        cv: R * (-tau * tau * phi_tautau),
+        speedOfSound: Math.sqrt(1000 * R * T * (2 * delta * phi_delta +
+            delta * delta * phi_deltadelta -
+            Math.pow(delta * phi_delta - delta * tau * phi_deltatau, 2) /
+                (tau * tau * phi_tautau))),
+        quality: null,
+        isobaricExpansion: (phi_delta - tau * phi_deltatau) /
+            (2 * phi_delta + delta * phi_deltadelta) / T,
+        isothermalCompressibility: 1 /
+            (2 * delta * phi_delta + delta * delta * phi_deltadelta) /
+            rho / R / T * 1000,
+    };
+}
+//# sourceMappingURL=region3.js.map

package/dist/regions/region3.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region3.js","sourceRoot":"","sources":["../../src/regions/region3.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,CAAC,EAAE,MAAM,iBAAiB,CAAC;AAEpC,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAErC,8EAA8E;AAE9E,MAAM,YAAY,GAAqB;IACrC,CAAC,CAAC,EAAE,CAAC,EAAE,eAAe,CAAC;IACvB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,gBAAgB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC;IACxB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,gBAAgB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC;IACxB,CAAC,CAAC,EAAE,EAAE,EAAE,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC;IACxB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,gBAAgB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC5B,CAAC,EAAE,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC5B,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;CACtB,CAAC;AAEX,mCAAmC;AACnC,MAAM,QAAQ,GAAG,GAAG,CAAC,CAAC,gCAAgC;AACtD,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,gCAAgC;AAExD;;;;;;GAMG;AACH,MAAM,UAAU,aAAa,CAAC,GAAW,EAAE,CAAS;IAClD,MAAM,KAAK,GAAG,GAAG,GAAG,QAAQ,CAAC;IAC7B,MAAM,GAAG,GAAG,MAAM,GAAG,CAAC,CAAC;IACvB,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;IAElB,kDAAkD;IAClD,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9B,IAAI,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC/B,IAAI,SAAS,GAAG,EAAE,GAAG,KAAK,CAAC;IAC3B,IAAI,cAAc,GAAG,CAAC,EAAE,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;IAC3C,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,YAAY,GAAG,CAAC,CAAC;IAErB,qDAAqD;IACrD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,YAAY,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;QACnD,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC;QACpC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QAEhC,GAAG,IAAI,CAAC,GAAG,QAAQ,GAAG,MAAM,CAAC;QAC7B,SAAS,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QACrD,cAAc,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QACpE,OAAO,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACnD,UAAU,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAChE,YAAY,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5E,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,SAAS,GAAG,IAAI;QAChD,WAAW,EAAE,CAAC;QACd,cAAc,EAAE,CAAC;QACjB,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,OAAO;QACrC,OAAO,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;QAClC,QAAQ,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,KAAK,GAAG,SAAS,CAAC;QACrD,EAAE,EAAE,CAAC,GAAG,CACN,CAAC,GAAG,GAAG,GAAG,GAAG,UAAU;YACvB,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,SAAS,GAAG,KAAK,GAAG,GAAG,GAAG,YAAY,EAAE,CAAC,CAAC;gBAC3D,CAAC,CAAC,GAAG,KAAK,GAAG,SAAS,GAAG,KAAK,GAAG,KAAK,GAAG,cAAc,CAAC,CACzD;QACD,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,UAAU,CAAC;QACjC,YAAY,EAAE,IAAI,CAAC,IAAI,CACrB,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CACb,CAAC,GAAG,KAAK,GAAG,SAAS;YACrB,KAAK,GAAG,KAAK,GAAG,cAAc;YAC9B,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,SAAS,GAAG,KAAK,GAAG,GAAG,GAAG,YAAY,EAAE,CAAC,CAAC;gBAC3D,CAAC,GAAG,GAAG,GAAG,GAAG,UAAU,CAAC,CACzB,CACF;QACD,OAAO,EAAE,IAAI;QACb,iBAAiB,EAAE,CAAC,SAAS,GAAG,GAAG,GAAG,YAAY,CAAC;YACjD,CAAC,CAAC,GAAG,SAAS,GAAG,KAAK,GAAG,cAAc,CAAC,GAAG,CAAC;QAC9C,yBAAyB,EAAE,CAAC;YAC1B,CAAC,CAAC,GAAG,KAAK,GAAG,SAAS,GAAG,KAAK,GAAG,KAAK,GAAG,cAAc,CAAC;YACxD,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI;KACrB,CAAC;AACJ,CAAC"}

package/dist/regions/region4.d.ts

@@ -0,0 +1,27 @@
+/**
+ * IAPWS-IF97 Region 4: Saturation Line
+ *
+ * Saturation-pressure equation (Eq. 30) and saturation-temperature equation (Eq. 31)
+ * Valid for: 273.15 K ≤ T ≤ 647.096 K
+ *
+ * Reference: IAPWS-IF97, Section 8 (Equations for Region 4)
+ */
+/**
+ * Saturation pressure for a given temperature.
+ * Equation 30 (pp. 33) of IAPWS-IF97.
+ *
+ * @param T - Temperature [K], 273.15 ≤ T ≤ 647.096
+ * @returns Saturation pressure [MPa]
+ * @throws {OutOfRangeError} if T is outside valid range
+ */
+export declare function saturationPressure(T: number): number;
+/**
+ * Saturation temperature for a given pressure.
+ * Equation 31 (pp. 34) of IAPWS-IF97.
+ *
+ * @param p - Pressure [MPa], Pt ≤ P ≤ Pc
+ * @returns Saturation temperature [K]
+ * @throws {OutOfRangeError} if P is outside valid range
+ */
+export declare function saturationTemperature(p: number): number;
+//# sourceMappingURL=region4.d.ts.map

package/dist/regions/region4.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4.d.ts","sourceRoot":"","sources":["../../src/regions/region4.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAqBH;;;;;;;GAOG;AACH,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAWpD;AAED;;;;;;;GAOG;AACH,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAcvD"}

package/dist/regions/region4.js

@@ -0,0 +1,62 @@
+/**
+ * IAPWS-IF97 Region 4: Saturation Line
+ *
+ * Saturation-pressure equation (Eq. 30) and saturation-temperature equation (Eq. 31)
+ * Valid for: 273.15 K ≤ T ≤ 647.096 K
+ *
+ * Reference: IAPWS-IF97, Section 8 (Equations for Region 4)
+ */
+import { Tc, Pt, Pc } from '../constants.js';
+import { OutOfRangeError } from '../types.js';
+// ─── Saturation Line Coefficients (Table 34, IAPWS-IF97) ───────────────────
+const n = [
+    0, // n[0] unused (1-indexed)
+    0.11670521452767e4,
+    -0.72421316703206e6,
+    -0.17073846940092e2,
+    0.12020824702470e5,
+    -0.32325550322333e7,
+    0.14915108613530e2,
+    -0.48232657361591e4,
+    0.40511340542057e6,
+    -0.23855557567849e0,
+    0.65017534844798e3,
+];
+/**
+ * Saturation pressure for a given temperature.
+ * Equation 30 (pp. 33) of IAPWS-IF97.
+ *
+ * @param T - Temperature [K], 273.15 ≤ T ≤ 647.096
+ * @returns Saturation pressure [MPa]
+ * @throws {OutOfRangeError} if T is outside valid range
+ */
+export function saturationPressure(T) {
+    if (T < 273.15 || T > Tc) {
+        throw new OutOfRangeError('Temperature', T, 273.15, Tc);
+    }
+    const theta = T + n[9] / (T - n[10]);
+    const A = theta * theta + n[1] * theta + n[2];
+    const B = n[3] * theta * theta + n[4] * theta + n[5];
+    const C = n[6] * theta * theta + n[7] * theta + n[8];
+    return Math.pow(2 * C / (-B + Math.sqrt(B * B - 4 * A * C)), 4);
+}
+/**
+ * Saturation temperature for a given pressure.
+ * Equation 31 (pp. 34) of IAPWS-IF97.
+ *
+ * @param p - Pressure [MPa], Pt ≤ P ≤ Pc
+ * @returns Saturation temperature [K]
+ * @throws {OutOfRangeError} if P is outside valid range
+ */
+export function saturationTemperature(p) {
+    if (p < Pt || p > Pc) {
+        throw new OutOfRangeError('Pressure', p, Pt, Pc);
+    }
+    const beta = Math.pow(p, 0.25);
+    const E = beta * beta + n[3] * beta + n[6];
+    const F = n[1] * beta * beta + n[4] * beta + n[7];
+    const G = n[2] * beta * beta + n[5] * beta + n[8];
+    const D = 2 * G / (-F - Math.sqrt(F * F - 4 * E * G));
+    return (n[10] + D - Math.sqrt(Math.pow(n[10] + D, 2) - 4 * (n[9] + n[10] * D))) / 2;
+}
+//# sourceMappingURL=region4.js.map

package/dist/regions/region4.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region4.js","sourceRoot":"","sources":["../../src/regions/region4.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAE9C,8EAA8E;AAE9E,MAAM,CAAC,GAAG;IACR,CAAC,EAAE,0BAA0B;IAC7B,kBAAkB;IAClB,CAAC,kBAAkB;IACnB,CAAC,kBAAkB;IACnB,kBAAkB;IAClB,CAAC,kBAAkB;IACnB,kBAAkB;IAClB,CAAC,kBAAkB;IACnB,kBAAkB;IAClB,CAAC,kBAAkB;IACnB,kBAAkB;CACV,CAAC;AAEX;;;;;;;GAOG;AACH,MAAM,UAAU,kBAAkB,CAAC,CAAS;IAC1C,IAAI,CAAC,GAAG,MAAM,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACzB,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,CAAC;IAC1D,CAAC;IAED,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,MAAM,CAAC,GAAG,KAAK,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAErD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAClE,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CAAC,CAAS;IAC7C,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,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IAC/B,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAEtD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAC3B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAChD,CAAC,GAAG,CAAC,CAAC;AACT,CAAC"}

package/dist/regions/region5.d.ts

@@ -0,0 +1,18 @@
+/**
+ * IAPWS-IF97 Region 5: High-Temperature Steam
+ *
+ * Basic equation for the specific Gibbs free energy γ(π,τ) = γ°(π,τ) + γʳ(π,τ)
+ * Valid for: 1073.15 K < T ≤ 2273.15 K, P ≤ 50 MPa
+ *
+ * Reference: IAPWS-IF97, Section 9 (Equations for Region 5)
+ */
+import type { BasicProperties } from '../types.js';
+/**
+ * Compute Region 5 thermodynamic properties from pressure and temperature.
+ *
+ * @param p - Pressure [MPa]
+ * @param T - Temperature [K]
+ * @returns Basic thermodynamic properties
+ */
+export declare function region5(p: number, T: number): BasicProperties;
+//# sourceMappingURL=region5.d.ts.map

package/dist/regions/region5.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region5.d.ts","sourceRoot":"","sources":["../../src/regions/region5.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,KAAK,EAAE,eAAe,EAAoB,MAAM,aAAa,CAAC;AAyBrE;;;;;;GAMG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,eAAe,CAoE7D"}

package/dist/regions/region5.js

@@ -0,0 +1,92 @@
+/**
+ * IAPWS-IF97 Region 5: High-Temperature Steam
+ *
+ * Basic equation for the specific Gibbs free energy γ(π,τ) = γ°(π,τ) + γʳ(π,τ)
+ * Valid for: 1073.15 K < T ≤ 2273.15 K, P ≤ 50 MPa
+ *
+ * Reference: IAPWS-IF97, Section 9 (Equations for Region 5)
+ */
+import { R } from '../constants.js';
+import { Region } from '../types.js';
+// ─── Ideal Part Coefficients (Table 37) ─────────────────────────────────────
+const IDEAL_COEFFICIENTS = [
+    [0, -13.179983674201],
+    [1, 6.8540841634434],
+    [-3, -0.024805148933466],
+    [-2, 0.36901534980333],
+    [-1, -3.1161318213925],
+    [2, -0.32961626538917],
+];
+// ─── Residual Part Coefficients (Table 38) ──────────────────────────────────
+const RESIDUAL_COEFFICIENTS = [
+    [1, 1, 1.5736404855259e-3],
+    [1, 2, 9.0153761673944e-4],
+    [1, 3, -5.0270077677648e-3],
+    [2, 3, 2.2440037409485e-6],
+    [2, 9, -4.1163275453471e-6],
+    [3, 7, 3.7919454822955e-8],
+];
+/**
+ * Compute Region 5 thermodynamic properties from pressure and temperature.
+ *
+ * @param p - Pressure [MPa]
+ * @param T - Temperature [K]
+ * @returns Basic thermodynamic properties
+ */
+export function region5(p, T) {
+    const pi = p;
+    const tau = 1000 / T;
+    // ── Ideal part ──
+    let gi = Math.log(pi);
+    let gi_pi = 1 / pi;
+    // gi_pipi = -1/(pi*pi) — not needed for Region 5 properties
+    let gi_tau = 0;
+    let gi_tautau = 0;
+    for (const [J, N] of IDEAL_COEFFICIENTS) {
+        gi += N * Math.pow(tau, J);
+        gi_tau += N * J * Math.pow(tau, J - 1);
+        gi_tautau += N * J * (J - 1) * Math.pow(tau, J - 2);
+    }
+    // ── Residual part ──
+    let gr = 0;
+    let gr_pi = 0;
+    let gr_pipi = 0;
+    let gr_tau = 0;
+    let gr_tautau = 0;
+    let gr_pitau = 0;
+    for (const [I, J, N] of RESIDUAL_COEFFICIENTS) {
+        const piPow = Math.pow(pi, I);
+        const tauPow = Math.pow(tau, J);
+        gr += N * piPow * tauPow;
+        gr_pi += N * I * Math.pow(pi, I - 1) * tauPow;
+        gr_pipi += N * I * (I - 1) * Math.pow(pi, I - 2) * tauPow;
+        gr_tau += N * piPow * J * Math.pow(tau, J - 1);
+        gr_tautau += N * piPow * J * (J - 1) * Math.pow(tau, J - 2);
+        gr_pitau += N * I * Math.pow(pi, I - 1) * J * Math.pow(tau, J - 1);
+    }
+    const gp = gi_pi + gr_pi;
+    const gt = gi_tau + gr_tau;
+    const gtt = gi_tautau + gr_tautau;
+    return {
+        region: Region.Region5,
+        pressure: p,
+        temperature: T,
+        specificVolume: R * T * pi * gp / (1000 * p),
+        internalEnergy: R * T * (tau * gt - pi * gp),
+        entropy: R * (tau * gt - (gi + gr)),
+        enthalpy: R * T * tau * gt,
+        cp: -R * tau * tau * gtt,
+        cv: R * (-tau * tau * gtt -
+            Math.pow(1 + pi * gr_pi - tau * pi * gr_pitau, 2) /
+                (1 - pi * pi * gr_pipi)),
+        speedOfSound: Math.sqrt(1000 * R * T *
+            (1 + 2 * pi * gr_pi + pi * pi * gr_pi * gr_pi) /
+            ((1 - pi * pi * gr_pipi) +
+                Math.pow(1 + pi * gr_pi - tau * pi * gr_pitau, 2) /
+                    (tau * tau * gtt))),
+        quality: null,
+        isobaricExpansion: (1 + pi * gr_pi - tau * pi * gr_pitau) / (1 + pi * gr_pi) / T,
+        isothermalCompressibility: (1 - pi * pi * gr_pipi) / (1 + pi * gr_pi) / p,
+    };
+}
+//# sourceMappingURL=region5.js.map

package/dist/regions/region5.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region5.js","sourceRoot":"","sources":["../../src/regions/region5.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,OAAO,EAAE,CAAC,EAAE,MAAM,iBAAiB,CAAC;AAEpC,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAErC,+EAA+E;AAE/E,MAAM,kBAAkB,GAA2C;IACjE,CAAC,CAAC,EAAE,CAAC,eAAe,CAAC;IACrB,CAAC,CAAC,EAAE,eAAe,CAAC;IACpB,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC;IACxB,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,eAAe,CAAC;IACtB,CAAC,CAAC,EAAE,CAAC,gBAAgB,CAAC;CACd,CAAC;AAEX,+EAA+E;AAE/E,MAAM,qBAAqB,GAAqB;IAC9C,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;CAClB,CAAC;AAEX;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,MAAM,EAAE,GAAG,CAAC,CAAC;IACb,MAAM,GAAG,GAAG,IAAI,GAAG,CAAC,CAAC;IAErB,mBAAmB;IACnB,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,GAAG,EAAE,CAAC;IACnB,4DAA4D;IAC5D,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,SAAS,GAAG,CAAC,CAAC;IAElB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,kBAAkB,EAAE,CAAC;QACxC,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QAC3B,MAAM,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACvC,SAAS,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACtD,CAAC;IAED,sBAAsB;IACtB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,IAAI,QAAQ,GAAG,CAAC,CAAC;IAEjB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,qBAAqB,EAAE,CAAC;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC9B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QAEhC,EAAE,IAAI,CAAC,GAAG,KAAK,GAAG,MAAM,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAC9C,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAC1D,MAAM,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/C,SAAS,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5D,QAAQ,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACrE,CAAC;IAED,MAAM,EAAE,GAAG,KAAK,GAAG,KAAK,CAAC;IACzB,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;IAC3B,MAAM,GAAG,GAAG,SAAS,GAAG,SAAS,CAAC;IAElC,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ,EAAE,CAAC;QACX,WAAW,EAAE,CAAC;QACd,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC;QAC5C,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC5C,OAAO,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACnC,QAAQ,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,EAAE;QAC1B,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG;QACxB,EAAE,EAAE,CAAC,GAAG,CACN,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG;YAChB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,GAAG,GAAG,EAAE,GAAG,QAAQ,EAAE,CAAC,CAAC;gBACjD,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,OAAO,CAAC,CACxB;QACD,YAAY,EAAE,IAAI,CAAC,IAAI,CACrB,IAAI,GAAG,CAAC,GAAG,CAAC;YACZ,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,KAAK,CAAC;YAC9C,CACE,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,OAAO,CAAC;gBACvB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,GAAG,GAAG,EAAE,GAAG,QAAQ,EAAE,CAAC,CAAC;oBACjD,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAClB,CACF;QACD,OAAO,EAAE,IAAI;QACb,iBAAiB,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,GAAG,GAAG,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC;QAChF,yBAAyB,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC;KAC1E,CAAC;AACJ,CAAC"}

package/dist/saturation/common.d.ts

@@ -0,0 +1,12 @@
+import type { BasicProperties } from '../types.js';
+export interface SaturationEndpoints {
+    pressure: number;
+    temperature: number;
+    liquid: BasicProperties;
+    vapor: BasicProperties;
+}
+export declare function saturationEndpointsAtPressure(p: number): SaturationEndpoints;
+export declare function saturationEndpointsAtTemperature(T: number): SaturationEndpoints;
+export declare function qualityFromSaturationProperty(liquidValue: number, vaporValue: number, mixtureValue: number): number;
+export declare function mixSaturationState(endpoints: SaturationEndpoints, qualityInput: number): BasicProperties;
+//# sourceMappingURL=common.d.ts.map

package/dist/saturation/common.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"common.d.ts","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAUnD,MAAM,WAAW,mBAAmB;IAClC,QAAQ,EAAE,MAAM,CAAC;IACjB,WAAW,EAAE,MAAM,CAAC;IACpB,MAAM,EAAE,eAAe,CAAC;IACxB,KAAK,EAAE,eAAe,CAAC;CACxB;AAoBD,wBAAgB,6BAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,mBAAmB,CAuB5E;AAED,wBAAgB,gCAAgC,CAAC,CAAC,EAAE,MAAM,GAAG,mBAAmB,CAE/E;AAED,wBAAgB,6BAA6B,CAC3C,WAAW,EAAE,MAAM,EACnB,UAAU,EAAE,MAAM,EAClB,YAAY,EAAE,MAAM,GACnB,MAAM,CAER;AAED,wBAAgB,kBAAkB,CAChC,SAAS,EAAE,mBAAmB,EAC9B,YAAY,EAAE,MAAM,GACnB,eAAe,CA2BjB"}

package/dist/saturation/common.js

@@ -0,0 +1,75 @@
+import { Region } from '../types.js';
+import { Pc, R2_T_MIN, RHOc } 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';
+function withRegion4Metadata(state, quality) {
+    return { ...state, region: Region.Region4, quality };
+}
+function solveR3Density(p, T, x) {
+    const v0 = region3SatVolume(p, T, x);
+    return newtonRaphson((rho) => region3ByRhoT(rho, T).pressure - p, 1 / v0);
+}
+function clampQuality(x) {
+    if (x <= 0 && x > -1e-12)
+        return 0;
+    if (x >= 1 && x < 1 + 1e-12)
+        return 1;
+    return Math.max(0, Math.min(1, x));
+}
+export function saturationEndpointsAtPressure(p) {
+    const temperature = saturationTemperature(p);
+    if (Math.abs(p - Pc) < 1e-10) {
+        const state = region3ByRhoT(RHOc, temperature);
+        return { pressure: p, temperature, liquid: state, vapor: state };
+    }
+    if (temperature > R2_T_MIN) {
+        return {
+            pressure: p,
+            temperature,
+            liquid: region3ByRhoT(solveR3Density(p, temperature, 0), temperature),
+            vapor: region3ByRhoT(solveR3Density(p, temperature, 1), temperature),
+        };
+    }
+    return {
+        pressure: p,
+        temperature,
+        liquid: region1(p, temperature),
+        vapor: region2(p, temperature),
+    };
+}
+export function saturationEndpointsAtTemperature(T) {
+    return saturationEndpointsAtPressure(saturationPressure(T));
+}
+export function qualityFromSaturationProperty(liquidValue, vaporValue, mixtureValue) {
+    return clampQuality((mixtureValue - liquidValue) / (vaporValue - liquidValue));
+}
+export function mixSaturationState(endpoints, qualityInput) {
+    const quality = clampQuality(qualityInput);
+    const { liquid, vapor, pressure, temperature } = endpoints;
+    if (quality === 0) {
+        return withRegion4Metadata(liquid, 0);
+    }
+    if (quality === 1) {
+        return withRegion4Metadata(vapor, 1);
+    }
+    return {
+        region: Region.Region4,
+        pressure,
+        temperature,
+        specificVolume: liquid.specificVolume * (1 - quality) + vapor.specificVolume * quality,
+        internalEnergy: liquid.internalEnergy * (1 - quality) + vapor.internalEnergy * quality,
+        entropy: liquid.entropy * (1 - quality) + vapor.entropy * quality,
+        enthalpy: liquid.enthalpy * (1 - quality) + vapor.enthalpy * quality,
+        cp: NaN,
+        cv: NaN,
+        speedOfSound: NaN,
+        quality,
+        isobaricExpansion: NaN,
+        isothermalCompressibility: NaN,
+    };
+}
+//# sourceMappingURL=common.js.map

package/dist/saturation/common.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"common.js","sourceRoot":"","sources":["../../src/saturation/common.ts"],"names":[],"mappings":"AACA,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;AAS7D,SAAS,mBAAmB,CAAC,KAAsB,EAAE,OAAe;IAClE,OAAO,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC;AACvD,CAAC;AAED,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,SAAS,YAAY,CAAC,CAAS;IAC7B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK;QAAE,OAAO,CAAC,CAAC;IACnC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK;QAAE,OAAO,CAAC,CAAC;IACtC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED,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,MAAM,UAAU,gCAAgC,CAAC,CAAS;IACxD,OAAO,6BAA6B,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED,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,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,GAAG;QACP,EAAE,EAAE,GAAG;QACP,YAAY,EAAE,GAAG;QACjB,OAAO;QACP,iBAAiB,EAAE,GAAG;QACtB,yBAAyB,EAAE,GAAG;KAC/B,CAAC;AACJ,CAAC"}

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

@@ -0,0 +1,22 @@
+/**
+ * IAPWS-IF97 Two-Phase Saturation Solvers
+ *
+ * Compute thermodynamic state given pressure+quality (solvePx)
+ * or temperature+quality (solveTx).
+ *
+ * Reference: IAPWS-IF97, Region 4 (Saturation Line)
+ */
+import type { BasicProperties } from '../types.js';
+/**
+ * Solve for thermodynamic state given pressure and vapour quality.
+ * @param p - Pressure [MPa], Pt <= p <= Pc
+ * @param x - Vapour quality [0, 1]
+ */
+export declare function solvePx(p: number, x: number): BasicProperties;
+/**
+ * Solve for thermodynamic state given temperature and vapour quality.
+ * @param T - Temperature [K], T_MIN <= T <= Tc
+ * @param x - Vapour quality [0, 1]
+ */
+export declare function solveTx(T: number, x: number): BasicProperties;
+//# sourceMappingURL=two-phase.d.ts.map

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

@@ -0,0 +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"}

package/dist/saturation/two-phase.js

@@ -0,0 +1,44 @@
+/**
+ * IAPWS-IF97 Two-Phase Saturation Solvers
+ *
+ * Compute thermodynamic state given pressure+quality (solvePx)
+ * or temperature+quality (solveTx).
+ *
+ * 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';
+// ─── solvePx ────────────────────────────────────────────────────────────────
+/**
+ * Solve for thermodynamic state given pressure and vapour quality.
+ * @param p - Pressure [MPa], Pt <= p <= Pc
+ * @param x - Vapour quality [0, 1]
+ */
+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);
+}
+// ─── solveTx ────────────────────────────────────────────────────────────────
+/**
+ * Solve for thermodynamic state given temperature and vapour quality.
+ * @param T - Temperature [K], T_MIN <= T <= Tc
+ * @param x - Vapour quality [0, 1]
+ */
+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);
+    return solvePx(p, x);
+}
+//# sourceMappingURL=two-phase.js.map

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

@@ -0,0 +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"}

package/dist/solvers/nelder-mead.d.ts

@@ -0,0 +1,21 @@
+/**
+ * Nelder-Mead (downhill simplex) minimizer.
+ * Minimizes f(x) where x is an array of numbers.
+ */
+export interface NelderMeadOptions {
+    maxIterations?: number;
+    nonZeroDelta?: number;
+    zeroDelta?: number;
+    minErrorDelta?: number;
+    minTolerance?: number;
+    rho?: number;
+    chi?: number;
+    psi?: number;
+    sigma?: number;
+}
+export interface NelderMeadResult {
+    x: number[];
+    fx: number;
+}
+export declare function nelderMead(f: (x: number[]) => number, x0: number[], opts?: NelderMeadOptions): NelderMeadResult;
+//# sourceMappingURL=nelder-mead.d.ts.map

package/dist/solvers/nelder-mead.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"nelder-mead.d.ts","sourceRoot":"","sources":["../../src/solvers/nelder-mead.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,MAAM,WAAW,iBAAiB;IAChC,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,GAAG,CAAC,EAAE,MAAM,CAAC;IACb,KAAK,CAAC,EAAE,MAAM,CAAC;CAChB;AAED,MAAM,WAAW,gBAAgB;IAC/B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,wBAAgB,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,iBAAiB,GAAG,gBAAgB,CAiG/G"}

package/dist/solvers/nelder-mead.js

@@ -0,0 +1,99 @@
+/**
+ * Nelder-Mead (downhill simplex) minimizer.
+ * Minimizes f(x) where x is an array of numbers.
+ */
+export function nelderMead(f, x0, opts) {
+    const maxIter = opts?.maxIterations ?? x0.length * 200;
+    const nonZeroDelta = opts?.nonZeroDelta ?? 1.05;
+    const zeroDelta = opts?.zeroDelta ?? 0.001;
+    const minErrDelta = opts?.minErrorDelta ?? 1e-11;
+    const minTol = opts?.minTolerance ?? 1e-11;
+    const rho = opts?.rho ?? 1;
+    const chi = opts?.chi ?? 2;
+    const psi = opts?.psi ?? -0.5;
+    const sigma = opts?.sigma ?? 0.5;
+    const N = x0.length;
+    const simplex = new Array(N + 1);
+    // Initialize simplex
+    const v0 = [...x0];
+    v0.fx = f(x0);
+    simplex[0] = v0;
+    for (let i = 0; i < N; i++) {
+        const pt = [...x0];
+        pt[i] = pt[i] ? pt[i] * nonZeroDelta : zeroDelta;
+        pt.fx = f(pt);
+        simplex[i + 1] = pt;
+    }
+    const sortOrder = (a, b) => a.fx - b.fx;
+    function weightedSum(ret, w1, v1, w2, v2) {
+        for (let j = 0; j < ret.length; j++) {
+            ret[j] = w1 * v1[j] + w2 * v2[j];
+        }
+        return ret;
+    }
+    const centroid = [...x0];
+    const reflected = [...x0];
+    const contracted = [...x0];
+    const expanded = [...x0];
+    function updateSimplex(value) {
+        for (let i = 0; i < N; i++)
+            simplex[N][i] = value[i];
+        simplex[N].fx = value.fx;
+    }
+    for (let iter = 0; iter < maxIter; iter++) {
+        simplex.sort(sortOrder);
+        let maxDiff = 0;
+        for (let i = 0; i < N; i++) {
+            maxDiff = Math.max(maxDiff, Math.abs(simplex[0][i] - simplex[1][i]));
+        }
+        if (Math.abs(simplex[0].fx - simplex[N].fx) < minErrDelta && maxDiff < minTol)
+            break;
+        for (let i = 0; i < N; i++) {
+            centroid[i] = 0;
+            for (let j = 0; j < N; j++)
+                centroid[i] += simplex[j][i];
+            centroid[i] /= N;
+        }
+        const worst = simplex[N];
+        weightedSum(reflected, 1 + rho, centroid, -rho, worst);
+        reflected.fx = f(reflected);
+        if (reflected.fx < simplex[0].fx) {
+            weightedSum(expanded, 1 + chi, centroid, -chi, worst);
+            expanded.fx = f(expanded);
+            updateSimplex(expanded.fx < reflected.fx ? expanded : reflected);
+        }
+        else if (reflected.fx >= simplex[N - 1].fx) {
+            let shouldReduce = false;
+            if (reflected.fx > worst.fx) {
+                weightedSum(contracted, 1 + psi, centroid, -psi, worst);
+                contracted.fx = f(contracted);
+                if (contracted.fx < worst.fx)
+                    updateSimplex(contracted);
+                else
+                    shouldReduce = true;
+            }
+            else {
+                weightedSum(contracted, 1 - psi * rho, centroid, psi * rho, worst);
+                contracted.fx = f(contracted);
+                if (contracted.fx < reflected.fx)
+                    updateSimplex(contracted);
+                else
+                    shouldReduce = true;
+            }
+            if (shouldReduce) {
+                if (sigma >= 1)
+                    break;
+                for (let i = 1; i < simplex.length; i++) {
+                    weightedSum(simplex[i], 1 - sigma, simplex[0], sigma, simplex[i]);
+                    simplex[i].fx = f(simplex[i]);
+                }
+            }
+        }
+        else {
+            updateSimplex(reflected);
+        }
+    }
+    simplex.sort(sortOrder);
+    return { fx: simplex[0].fx, x: simplex[0].slice(0, N) };
+}
+//# sourceMappingURL=nelder-mead.js.map