iapws-if97 1.0.0

package/dist/solvers/nelder-mead.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"nelder-mead.js","sourceRoot":"","sources":["../../src/solvers/nelder-mead.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAmBH,MAAM,UAAU,UAAU,CAAC,CAA0B,EAAE,EAAY,EAAE,IAAwB;IAC3F,MAAM,OAAO,GAAG,IAAI,EAAE,aAAa,IAAI,EAAE,CAAC,MAAM,GAAG,GAAG,CAAC;IACvD,MAAM,YAAY,GAAG,IAAI,EAAE,YAAY,IAAI,IAAI,CAAC;IAChD,MAAM,SAAS,GAAG,IAAI,EAAE,SAAS,IAAI,KAAK,CAAC;IAC3C,MAAM,WAAW,GAAG,IAAI,EAAE,aAAa,IAAI,KAAK,CAAC;IACjD,MAAM,MAAM,GAAG,IAAI,EAAE,YAAY,IAAI,KAAK,CAAC;IAC3C,MAAM,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAI,EAAE,KAAK,IAAI,GAAG,CAAC;IAEjC,MAAM,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC;IAGpB,MAAM,OAAO,GAAa,IAAI,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3C,qBAAqB;IACrB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAW,CAAC;IAC7B,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACd,OAAO,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IAEhB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAW,CAAC;QAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,SAAS,CAAC;QACjD,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QACd,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IACtB,CAAC;IAED,MAAM,SAAS,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;IAExD,SAAS,WAAW,CAAC,GAAa,EAAE,EAAU,EAAE,EAAY,EAAE,EAAU,EAAE,EAAY;QACpF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACnC,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IAED,MAAM,QAAQ,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;IACzB,MAAM,SAAS,GAAG,CAAC,GAAG,EAAE,CAAW,CAAC;IACpC,MAAM,UAAU,GAAG,CAAC,GAAG,EAAE,CAAW,CAAC;IACrC,MAAM,QAAQ,GAAG,CAAC,GAAG,EAAE,CAAW,CAAC;IAEnC,SAAS,aAAa,CAAC,KAAa;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACrD,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAC3B,CAAC;IAED,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;QAC1C,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;QAExB,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvE,CAAC;QACD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,WAAW,IAAI,OAAO,GAAG,MAAM;YAAE,MAAM;QAErF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,QAAQ,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzD,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACnB,CAAC;QAED,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACzB,WAAW,CAAC,SAAS,EAAE,CAAC,GAAG,GAAG,EAAE,QAAQ,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;QACvD,SAAS,CAAC,EAAE,GAAG,CAAC,CAAC,SAAS,CAAC,CAAC;QAE5B,IAAI,SAAS,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC;YACjC,WAAW,CAAC,QAAQ,EAAE,CAAC,GAAG,GAAG,EAAE,QAAQ,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;YACtD,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;YAC1B,aAAa,CAAC,QAAQ,CAAC,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC;QACnE,CAAC;aAAM,IAAI,SAAS,CAAC,EAAE,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC;YAC7C,IAAI,YAAY,GAAG,KAAK,CAAC;YACzB,IAAI,SAAS,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,CAAC;gBAC5B,WAAW,CAAC,UAAU,EAAE,CAAC,GAAG,GAAG,EAAE,QAAQ,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;gBACxD,UAAU,CAAC,EAAE,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC;gBAC9B,IAAI,UAAU,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE;oBAAE,aAAa,CAAC,UAAU,CAAC,CAAC;;oBACnD,YAAY,GAAG,IAAI,CAAC;YAC3B,CAAC;iBAAM,CAAC;gBACN,WAAW,CAAC,UAAU,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,QAAQ,EAAE,GAAG,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC;gBACnE,UAAU,CAAC,EAAE,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC;gBAC9B,IAAI,UAAU,CAAC,EAAE,GAAG,SAAS,CAAC,EAAE;oBAAE,aAAa,CAAC,UAAU,CAAC,CAAC;;oBACvD,YAAY,GAAG,IAAI,CAAC;YAC3B,CAAC;YACD,IAAI,YAAY,EAAE,CAAC;gBACjB,IAAI,KAAK,IAAI,CAAC;oBAAE,MAAM;gBACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBACxC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBAClE,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChC,CAAC;YACH,CAAC;QACH,CAAC;aAAM,CAAC;YACN,aAAa,CAAC,SAAS,CAAC,CAAC;QAC3B,CAAC;IACH,CAAC;IAED,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IACxB,OAAO,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;AAC1D,CAAC"}

package/dist/solvers/newton-raphson.d.ts

@@ -0,0 +1,25 @@
+/**
+ * Newton-Raphson root-finding solver.
+ * Finds x such that f(x) = 0, with optional analytical derivative.
+ */
+export interface NewtonRaphsonOptions {
+    /** Convergence tolerance (default 1e-7) */
+    tolerance?: number;
+    /** Maximum iterations (default 50) */
+    maxIterations?: number;
+    /** Step size for numerical derivative (default 1e-4) */
+    h?: number;
+    /** Machine epsilon (default 2.22e-16) */
+    epsilon?: number;
+}
+/**
+ * Solve f(x) = 0 using Newton-Raphson method.
+ *
+ * @param f  - Function to find root of
+ * @param x0 - Initial guess
+ * @param fp - Optional analytical derivative f'(x)
+ * @param options - Solver options
+ * @returns x such that f(x) ≈ 0
+ */
+export declare function newtonRaphson(f: (x: number) => number, x0: number, fp?: ((x: number) => number) | null, options?: NewtonRaphsonOptions): number;
+//# sourceMappingURL=newton-raphson.d.ts.map

package/dist/solvers/newton-raphson.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"newton-raphson.d.ts","sourceRoot":"","sources":["../../src/solvers/newton-raphson.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,MAAM,WAAW,oBAAoB;IACnC,2CAA2C;IAC3C,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,sCAAsC;IACtC,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB,wDAAwD;IACxD,CAAC,CAAC,EAAE,MAAM,CAAC;IACX,yCAAyC;IACzC,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED;;;;;;;;GAQG;AACH,wBAAgB,aAAa,CAC3B,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EACxB,EAAE,EAAE,MAAM,EACV,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC,GAAG,IAAI,EACnC,OAAO,CAAC,EAAE,oBAAoB,GAC7B,MAAM,CA8BR"}

package/dist/solvers/newton-raphson.js

@@ -0,0 +1,44 @@
+/**
+ * Newton-Raphson root-finding solver.
+ * Finds x such that f(x) = 0, with optional analytical derivative.
+ */
+/**
+ * Solve f(x) = 0 using Newton-Raphson method.
+ *
+ * @param f  - Function to find root of
+ * @param x0 - Initial guess
+ * @param fp - Optional analytical derivative f'(x)
+ * @param options - Solver options
+ * @returns x such that f(x) ≈ 0
+ */
+export function newtonRaphson(f, x0, fp, options) {
+    const tol = options?.tolerance ?? 1e-7;
+    const maxIter = options?.maxIterations ?? 50;
+    const h = options?.h ?? 1e-4;
+    const eps = options?.epsilon ?? 2.220446049250313e-16;
+    const hr = 1 / h;
+    let x = x0;
+    for (let iter = 0; iter < maxIter; iter++) {
+        const y = f(x);
+        let yp;
+        if (fp) {
+            yp = fp(x);
+        }
+        else {
+            // 5-point numerical derivative
+            const yph = f(x + h);
+            const ymh = f(x - h);
+            const yp2h = f(x + 2 * h);
+            const ym2h = f(x - 2 * h);
+            yp = ((ym2h - yp2h) + 8 * (yph - ymh)) * hr / 12;
+        }
+        if (Math.abs(yp) <= eps * Math.abs(y))
+            return x;
+        const x1 = x - y / yp;
+        if (Math.abs(x1 - x) <= tol * Math.abs(x1))
+            return x1;
+        x = x1;
+    }
+    return x;
+}
+//# sourceMappingURL=newton-raphson.js.map

package/dist/solvers/newton-raphson.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"newton-raphson.js","sourceRoot":"","sources":["../../src/solvers/newton-raphson.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAaH;;;;;;;;GAQG;AACH,MAAM,UAAU,aAAa,CAC3B,CAAwB,EACxB,EAAU,EACV,EAAmC,EACnC,OAA8B;IAE9B,MAAM,GAAG,GAAG,OAAO,EAAE,SAAS,IAAI,IAAI,CAAC;IACvC,MAAM,OAAO,GAAG,OAAO,EAAE,aAAa,IAAI,EAAE,CAAC;IAC7C,MAAM,CAAC,GAAG,OAAO,EAAE,CAAC,IAAI,IAAI,CAAC;IAC7B,MAAM,GAAG,GAAG,OAAO,EAAE,OAAO,IAAI,qBAAqB,CAAC;IACtD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IAEjB,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,EAAU,CAAC;QAEf,IAAI,EAAE,EAAE,CAAC;YACP,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACb,CAAC;aAAM,CAAC;YACN,+BAA+B;YAC/B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACrB,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACrB,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1B,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1B,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;QACnD,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC;QAEhD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QACtB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAAE,OAAO,EAAE,CAAC;QACtD,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/dist/transport/properties.d.ts

@@ -0,0 +1,47 @@
+/**
+ * IAPWS Transport Properties
+ *
+ * - Viscosity (IAPWS 2008)
+ * - Thermal Conductivity (IAPWS 2011)
+ * - Surface Tension (IAPWS 2014)
+ * - Dielectric Constant (IAPWS 1997)
+ * - Ionization Constant (IAPWS 2007)
+ */
+/**
+ * Dynamic viscosity of water/steam [Pa·s]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export declare function viscosity(T: number, rho: number): number;
+/**
+ * Thermal conductivity of water/steam [W/(m·K)]
+ *
+ * When optional thermodynamic state parameters are provided, the full IAPWS 2011
+ * critical enhancement (λ₂) is computed. Otherwise λ₂ = 0 (backward compatible).
+ *
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ * @param cp  Isobaric heat capacity [kJ/(kg·K)] (optional)
+ * @param cv  Isochoric heat capacity [kJ/(kg·K)] (optional)
+ * @param drhodP_T  Isothermal density derivative ∂ρ/∂P|_T [kg/m³/MPa] (optional)
+ * @param mu  Dynamic viscosity [Pa·s] (optional)
+ */
+export declare function thermalConductivity(T: number, rho: number, cp?: number, cv?: number, drhodP_T?: number, mu?: number): number;
+/**
+ * Surface tension of water [N/m]
+ * @param T Temperature [K], valid for 273.15 ≤ T ≤ 647.096
+ */
+export declare function surfaceTension(T: number): number;
+/**
+ * Static dielectric constant [dimensionless]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export declare function dielectricConstant(T: number, rho: number): number;
+/**
+ * Ionization constant pKw [dimensionless]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export declare function ionizationConstant(T: number, rho: number): number;
+//# sourceMappingURL=properties.d.ts.map

package/dist/transport/properties.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"properties.d.ts","sourceRoot":"","sources":["../../src/transport/properties.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAiBH;;;;GAIG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,CA4BxD;AAyBD;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CACjC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACtB,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,QAAQ,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,GACvD,MAAM,CA0BR;AAgDD;;;GAGG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAGhD;AAUD;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,CAoBjE;AAMD;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,CAgBjE"}

package/dist/transport/properties.js

@@ -0,0 +1,205 @@
+/**
+ * IAPWS Transport Properties
+ *
+ * - Viscosity (IAPWS 2008)
+ * - Thermal Conductivity (IAPWS 2011)
+ * - Surface Tension (IAPWS 2014)
+ * - Dielectric Constant (IAPWS 1997)
+ * - Ionization Constant (IAPWS 2007)
+ */
+import { Tc, Pc, RHOc } from '../constants.js';
+// ─── Viscosity (IAPWS 2008) ─────────────────────────────────────────────────
+const MU0_H = [1.67752, 2.20462, 0.6366564, -0.241605];
+const MU1_H = [
+    5.20094e-1, 8.50895e-2, -1.08374, -2.89555e-1, 0, 0,
+    2.22531e-1, 9.99115e-1, 1.88797, 1.26613, 0, 1.20573e-1,
+    -2.81378e-1, -9.06851e-1, -7.72479e-1, -4.89837e-1, -2.5704e-1, 0,
+    1.61913e-1, 2.57399e-1, 0, 0, 0, 0,
+    -3.25372e-2, 0, 0, 6.98452e-2, 0, 0,
+    0, 0, 0, 0, 8.72102e-3, 0,
+    0, 0, 0, -4.35673e-3, 0, -5.93264e-4,
+];
+/**
+ * Dynamic viscosity of water/steam [Pa·s]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export function viscosity(T, rho) {
+    const That = T / Tc;
+    const rhat = rho / RHOc;
+    // μ₀ — dilute gas contribution
+    let sum0 = 0;
+    for (let i = 0; i < 4; i++)
+        sum0 += MU0_H[i] / Math.pow(That, i);
+    const mu0 = 100 * Math.sqrt(That) / sum0;
+    // μ₁ — finite-density contribution
+    let mu1sum = 0;
+    for (let j = 0; j < 6; j++) {
+        const x = Math.pow(1 / That - 1, j);
+        let y = 0;
+        for (let z = 0; z < 7; z++) {
+            y += MU1_H[z * 6 + j] * Math.pow(rhat - 1, z);
+        }
+        mu1sum += x * y;
+    }
+    const mu1 = Math.exp(rhat * mu1sum);
+    // μ₂ — critical enhancement
+    // Per IAPWS 2008, μ₂ is negligible except extremely close to the critical point
+    // (within ~1 K and ~1 MPa). The original reference library also omits it.
+    // Keeping μ₂ = 1 (no enhancement) for practical use.
+    const mu2 = 1;
+    return mu0 * mu1 * mu2 / 1e6; // Convert μPa·s to Pa·s
+}
+// ─── Thermal Conductivity (IAPWS 2011) ──────────────────────────────────────
+const K0_L = [2.443221e-3, 1.323095e-2, 6.770357e-3, -3.454586e-3, 4.096266e-4];
+const K1_L = [
+    1.60397357, -0.646013523, 0.111443906, 0.102997357, -0.0504123634, 0.00609859258,
+    2.33771842, -2.78843778, 1.53616167, -0.463045512, 0.0832827019, -0.00719201245,
+    2.19650529, -4.54580785, 3.55777244, -1.40944978, 0.275418278, -0.0205938816,
+    -1.21051378, 1.60812989, -0.621178141, 0.0716373224, 0, 0,
+    -2.720337, 4.57586331, -3.18369245, 1.1168348, -0.19268305, 0.012913842,
+];
+// Reference compressibility piecewise polynomial coefficients for λ₂ (IAPWS 2011)
+const DRHO_REF_COEFFS = [
+    [0.310559006, [6.53786807199516, -5.61149954923348, 3.39624167361325, -2.27492629730878, 10.2631854662709, 1.97815050331519]],
+    [0.776397516, [6.52717759281799, -6.30816983387575, 8.08379285492595, -9.82240510197603, 12.1358413791395, -5.54349664571295]],
+    [1.242236025, [5.35500529896124, -3.96415689925446, 8.91990208918795, -12.0338729505790, 9.19494865194302, -2.16866274479712]],
+    [1.863354037, [1.55225959906681, 0.464621290821181, 8.93237374861479, -11.0321960061126, 6.16780999933360, -0.965458722086812]],
+    [Infinity, [1.11999926419994, 0.595748562571649, 9.88952565078920, -10.3255051147040, 4.66861294457414, -0.503243546373828]],
+];
+/** Specific gas constant for water [kJ/(kg·K)] */
+const R_WATER = 0.46151805;
+/**
+ * Thermal conductivity of water/steam [W/(m·K)]
+ *
+ * When optional thermodynamic state parameters are provided, the full IAPWS 2011
+ * critical enhancement (λ₂) is computed. Otherwise λ₂ = 0 (backward compatible).
+ *
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ * @param cp  Isobaric heat capacity [kJ/(kg·K)] (optional)
+ * @param cv  Isochoric heat capacity [kJ/(kg·K)] (optional)
+ * @param drhodP_T  Isothermal density derivative ∂ρ/∂P|_T [kg/m³/MPa] (optional)
+ * @param mu  Dynamic viscosity [Pa·s] (optional)
+ */
+export function thermalConductivity(T, rho, cp, cv, drhodP_T, mu) {
+    const That = T / Tc;
+    const rhat = rho / RHOc;
+    // λ₀
+    let sum0 = 0;
+    for (let i = 0; i < 5; i++)
+        sum0 += K0_L[i] / Math.pow(That, i);
+    const k0 = Math.sqrt(That) / sum0;
+    // λ₁
+    let k1sum = 0;
+    for (let j = 0; j < 5; j++) {
+        const x = Math.pow(1 / That - 1, j);
+        let y = 0;
+        for (let z = 0; z < 6; z++)
+            y += K1_L[j * 6 + z] * Math.pow(rhat - 1, z);
+        k1sum += x * y;
+    }
+    const k1 = Math.exp(rhat * k1sum);
+    // λ₂ — critical enhancement (IAPWS 2011)
+    let k2 = 0;
+    if (cp !== undefined && cv !== undefined && drhodP_T !== undefined && mu !== undefined) {
+        k2 = computeK2(That, rhat, cp, cv, drhodP_T, mu);
+    }
+    return (k0 * k1 + k2) / 1000; // Convert mW/(m·K) to W/(m·K)
+}
+/**
+ * Compute the critical enhancement term λ₂ [mW/(m·K)] per IAPWS 2011.
+ */
+function computeK2(That, rhat, cp, cv, drhodP_T, mu) {
+    // Reference compressibility from piecewise polynomial
+    let ai = DRHO_REF_COEFFS[DRHO_REF_COEFFS.length - 1][1];
+    for (const [threshold, coeffs] of DRHO_REF_COEFFS) {
+        if (rhat <= threshold) {
+            ai = coeffs;
+            break;
+        }
+    }
+    let tempD = 0;
+    for (let i = 0; i < 6; i++)
+        tempD += ai[i] * Math.pow(rhat, i);
+    const drho_ref = (1 / tempD) * (RHOc / Pc); // [kg/m³/MPa]
+    // Compressibility difference
+    let DeltaX = rhat * (Pc / RHOc * drhodP_T - Pc / RHOc * drho_ref * 1.5 / That);
+    if (DeltaX < 0)
+        DeltaX = 0;
+    // Correlation length ξ [nm]
+    const X = 0.13 * Math.pow(DeltaX / 0.06, 0.63 / 1.239);
+    // qD·ξ
+    const Y = X / 0.4;
+    let Z;
+    if (Y < 1.2e-7) {
+        Z = 0;
+    }
+    else {
+        const kappa = cp / cv;
+        Z = 2 / (Math.PI * Y) * (((1 - 1 / kappa) * Math.atan(Y) + Y / kappa)
+            - (1 - Math.exp(-1 / (1 / Y + Y * Y / (3 * rhat * rhat)))));
+    }
+    // Viscosity in μPa·s (mu is in Pa·s)
+    const mu_uPas = mu * 1e6;
+    // λ₂ in mW/(m·K)
+    return 177.8514 * rhat * (cp / R_WATER) * That / mu_uPas * Z;
+}
+// ─── Surface Tension (IAPWS 2014) ───────────────────────────────────────────
+/**
+ * Surface tension of water [N/m]
+ * @param T Temperature [K], valid for 273.15 ≤ T ≤ 647.096
+ */
+export function surfaceTension(T) {
+    const tau = 1 - T / Tc;
+    return 235.8e-3 * Math.pow(tau, 1.256) * (1 - 0.625 * tau);
+}
+// ─── Dielectric Constant (IAPWS 1997) ───────────────────────────────────────
+const DC_Nh = [0.978224486826, -0.957771379375, 0.237511794148, 0.714692244396,
+    -0.298217036956, -0.108863472196, 0.0949327488264, -0.00980469816509,
+    0.16516763497e-4, 0.937359795772e-4, -0.12317921872e-9];
+const DC_Ih = [1, 1, 1, 2, 3, 3, 4, 5, 6, 7, 10];
+const DC_Jh = [0.25, 1, 2.5, 1.5, 1.5, 2.5, 2, 2, 5, 0.5, 10];
+/**
+ * Static dielectric constant [dimensionless]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export function dielectricConstant(T, rho) {
+    const MW = 0.018015268;
+    const rhoM = rho / MW;
+    const rhoRc = rhoM / (RHOc / MW);
+    const alpha = 1.636e-40;
+    const Na = 6.0221367e23;
+    const mu = 6.138e-30;
+    const eps0 = 1 / (4e-7 * Math.PI * 299792458 ** 2);
+    const k = 1.380658e-23;
+    let g = 0;
+    for (let i = 0; i < 11; i++) {
+        g += DC_Nh[i] * Math.pow(rhoRc, DC_Ih[i]) * Math.pow(Tc / T, DC_Jh[i]);
+    }
+    g = 1 + g + 0.00196096504426 * rhoRc * Math.pow(T / 228 - 1, -1.2);
+    const A = Na * mu * mu * rhoM * g / (eps0 * k * T);
+    const B = Na * alpha * rhoM / (3 * eps0);
+    return (1 + A + 5 * B + Math.sqrt(9 + 2 * A + 18 * B + A * A + 10 * A * B + 9 * B * B)) / (4 - 4 * B);
+}
+// ─── Ionization Constant (IAPWS 2024) ───────────────────────────────────────
+function log10(x) { return Math.log(x) / Math.LN10; }
+/**
+ * Ionization constant pKw [dimensionless]
+ * @param T   Temperature [K]
+ * @param rho Density [kg/m³]
+ */
+export function ionizationConstant(T, rho) {
+    const density = rho / 1000;
+    const n = 6;
+    const molarMass = 18.015268;
+    const pKwG = 0.61415 + 48251.33 / T - 67707.93 / (T * T) + 10102100 / (T * T * T);
+    const Z = density * Math.exp(-0.702132
+        + 8681.05 / T
+        - 24145.1 * Math.pow(density, 2 / 3) / (T * T));
+    const beta = 0.813876 - 51.4471 / T - 0.469920 * density;
+    return -2 * n * (log10(1 + Z)
+        - Z / (Z + 1) * density * beta) + pKwG + 2 * log10(molarMass / 1000);
+}
+//# sourceMappingURL=properties.js.map

package/dist/transport/properties.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"properties.js","sourceRoot":"","sources":["../../src/transport/properties.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,MAAM,iBAAiB,CAAC;AAE/C,+EAA+E;AAE/E,MAAM,KAAK,GAAG,CAAC,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,CAAC,QAAQ,CAAC,CAAC;AACvD,MAAM,KAAK,GAAG;IACZ,UAAU,EAAE,UAAU,EAAE,CAAC,OAAO,EAAE,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC;IACnD,UAAU,EAAE,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,EAAE,UAAU;IACvD,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,SAAS,EAAE,CAAC;IACjE,UAAU,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;IAClC,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC;IACnC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC;IACzB,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,UAAU;CACrC,CAAC;AAEF;;;;GAIG;AACH,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,GAAW;IAC9C,MAAM,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC;IACpB,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;IAExB,+BAA+B;IAC/B,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACjE,MAAM,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;IAEzC,mCAAmC;IACnC,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,IAAI,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,CAAC;QACD,MAAM,IAAI,CAAC,GAAG,CAAC,CAAC;IAClB,CAAC;IACD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,MAAM,CAAC,CAAC;IAEpC,4BAA4B;IAC5B,gFAAgF;IAChF,0EAA0E;IAC1E,qDAAqD;IACrD,MAAM,GAAG,GAAG,CAAC,CAAC;IAEd,OAAO,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,wBAAwB;AACxD,CAAC;AAED,+EAA+E;AAE/E,MAAM,IAAI,GAAG,CAAC,WAAW,EAAE,WAAW,EAAE,WAAW,EAAE,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;AAChF,MAAM,IAAI,GAAG;IACX,UAAU,EAAE,CAAC,WAAW,EAAE,WAAW,EAAE,WAAW,EAAE,CAAC,YAAY,EAAE,aAAa;IAChF,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,EAAE,CAAC,WAAW,EAAE,YAAY,EAAE,CAAC,aAAa;IAC/E,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,EAAE,CAAC,UAAU,EAAE,WAAW,EAAE,CAAC,YAAY;IAC5E,CAAC,UAAU,EAAE,UAAU,EAAE,CAAC,WAAW,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC;IACzD,CAAC,QAAQ,EAAE,UAAU,EAAE,CAAC,UAAU,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,WAAW;CACxE,CAAC;AAEF,kFAAkF;AAClF,MAAM,eAAe,GAAyB;IAC5C,CAAC,WAAW,EAAE,CAAC,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,gBAAgB,CAAC,CAAC;IAC7H,CAAC,WAAW,EAAE,CAAC,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,CAAC,CAAC;IAC9H,CAAC,WAAW,EAAE,CAAC,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,CAAC,CAAC;IAC9H,CAAC,WAAW,EAAE,CAAC,gBAAgB,EAAE,iBAAiB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,iBAAiB,CAAC,CAAC;IAC/H,CAAC,QAAQ,EAAE,CAAC,gBAAgB,EAAE,iBAAiB,EAAE,gBAAgB,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,CAAC,iBAAiB,CAAC,CAAC;CAC7H,CAAC;AAEF,kDAAkD;AAClD,MAAM,OAAO,GAAG,UAAU,CAAC;AAE3B;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CACjC,CAAS,EAAE,GAAW,EACtB,EAAW,EAAE,EAAW,EAAE,QAAiB,EAAE,EAAW;IAExD,MAAM,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC;IACpB,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;IAExB,KAAK;IACL,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,IAAI,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAChE,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;IAElC,KAAK;IACL,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACzE,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC;IACjB,CAAC;IACD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;IAElC,yCAAyC;IACzC,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,IAAI,QAAQ,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE,CAAC;QACvF,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,EAAE,CAAC,CAAC;IACnD,CAAC;IAED,OAAO,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,8BAA8B;AAC9D,CAAC;AAED;;GAEG;AACH,SAAS,SAAS,CAChB,IAAY,EAAE,IAAY,EAC1B,EAAU,EAAE,EAAU,EAAE,QAAgB,EAAE,EAAU;IAEpD,sDAAsD;IACtD,IAAI,EAAE,GAAa,eAAe,CAAC,eAAe,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,KAAK,MAAM,CAAC,SAAS,EAAE,MAAM,CAAC,IAAI,eAAe,EAAE,CAAC;QAClD,IAAI,IAAI,IAAI,SAAS,EAAE,CAAC;YAAC,EAAE,GAAG,MAAM,CAAC;YAAC,MAAM;QAAC,CAAC;IAChD,CAAC;IACD,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAC/D,MAAM,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,cAAc;IAE1D,6BAA6B;IAC7B,IAAI,MAAM,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,QAAQ,GAAG,EAAE,GAAG,IAAI,GAAG,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC;IAC/E,IAAI,MAAM,GAAG,CAAC;QAAE,MAAM,GAAG,CAAC,CAAC;IAE3B,4BAA4B;IAC5B,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,EAAE,IAAI,GAAG,KAAK,CAAC,CAAC;IAEvD,OAAO;IACP,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;IAElB,IAAI,CAAS,CAAC;IACd,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;QACf,CAAC,GAAG,CAAC,CAAC;IACR,CAAC;SAAM,CAAC;QACN,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;QACtB,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CACtB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;cAC1C,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAC3D,CAAC;IACJ,CAAC;IAED,qCAAqC;IACrC,MAAM,OAAO,GAAG,EAAE,GAAG,GAAG,CAAC;IAEzB,iBAAiB;IACjB,OAAO,QAAQ,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,IAAI,GAAG,OAAO,GAAG,CAAC,CAAC;AAC/D,CAAC;AAED,+EAA+E;AAE/E;;;GAGG;AACH,MAAM,UAAU,cAAc,CAAC,CAAS;IACtC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IACvB,OAAO,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,CAAC;AAC7D,CAAC;AAED,+EAA+E;AAE/E,MAAM,KAAK,GAAG,CAAC,cAAc,EAAE,CAAC,cAAc,EAAE,cAAc,EAAE,cAAc;IAC5E,CAAC,cAAc,EAAE,CAAC,cAAc,EAAE,eAAe,EAAE,CAAC,gBAAgB;IACpE,gBAAgB,EAAE,iBAAiB,EAAE,CAAC,gBAAgB,CAAC,CAAC;AAC1D,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;AACjD,MAAM,KAAK,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC;AAE9D;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,CAAS,EAAE,GAAW;IACvD,MAAM,EAAE,GAAG,WAAW,CAAC;IACvB,MAAM,IAAI,GAAG,GAAG,GAAG,EAAE,CAAC;IACtB,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,SAAS,CAAC;IACxB,MAAM,EAAE,GAAG,YAAY,CAAC;IACxB,MAAM,EAAE,GAAG,SAAS,CAAC;IACrB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,GAAG,SAAS,IAAI,CAAC,CAAC,CAAC;IACnD,MAAM,CAAC,GAAG,YAAY,CAAC;IAEvB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACzE,CAAC;IACD,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,gBAAgB,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IAEnE,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACnD,MAAM,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IAEzC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAC,CAAC,GAAG,EAAE,GAAC,CAAC,GAAG,CAAC,GAAC,CAAC,GAAG,EAAE,GAAC,CAAC,GAAC,CAAC,GAAG,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAC,CAAC,CAAC,CAAC;AACxF,CAAC;AAED,+EAA+E;AAE/E,SAAS,KAAK,CAAC,CAAS,IAAY,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;AAErE;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,CAAS,EAAE,GAAW;IACvD,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC;IAC3B,MAAM,CAAC,GAAG,CAAC,CAAC;IACZ,MAAM,SAAS,GAAG,SAAS,CAAC;IAC5B,MAAM,IAAI,GAAG,OAAO,GAAG,QAAQ,GAAG,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAClF,MAAM,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAC1B,CAAC,QAAQ;UACP,OAAO,GAAG,CAAC;UACX,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAC/C,CAAC;IACF,MAAM,IAAI,GAAG,QAAQ,GAAG,OAAO,GAAG,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC;IAEzD,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CACd,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;UACV,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,IAAI,CAC/B,GAAG,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC;AACzC,CAAC"}

package/dist/types.d.ts

@@ -0,0 +1,161 @@
+/**
+ * IAPWS-IF97 Core Types & Interfaces
+ *
+ * Type definitions for the IAPWS Industrial Formulation 1997
+ * for the Thermodynamic Properties of Water and Steam.
+ */
+/** IAPWS-IF97 thermodynamic regions */
+export declare enum Region {
+    /** Subcooled liquid (compressed water) */
+    Region1 = 1,
+    /** Superheated steam */
+    Region2 = 2,
+    /** Supercritical fluid */
+    Region3 = 3,
+    /** Two-phase (wet steam / saturation line) */
+    Region4 = 4,
+    /** High-temperature steam (T > 1073.15 K) */
+    Region5 = 5
+}
+/**
+ * Complete thermodynamic state of water/steam.
+ * All properties use SI-based units consistent with IAPWS-IF97.
+ */
+export interface SteamState {
+    /** IAPWS-IF97 region (1–5) */
+    region: Region;
+    /** Pressure [MPa] */
+    pressure: number;
+    /** Temperature [K] */
+    temperature: number;
+    /** Specific volume [m³/kg] */
+    specificVolume: number;
+    /** Density [kg/m³] */
+    density: number;
+    /** Specific internal energy [kJ/kg] */
+    internalEnergy: number;
+    /** Specific entropy [kJ/(kg·K)] */
+    entropy: number;
+    /** Specific enthalpy [kJ/kg] */
+    enthalpy: number;
+    /** Isobaric heat capacity [kJ/(kg·K)] */
+    cp: number;
+    /** Isochoric heat capacity [kJ/(kg·K)] */
+    cv: number;
+    /** Speed of sound [m/s] */
+    speedOfSound: number;
+    /** Vapour quality (0–1), only defined in Region 4 */
+    quality: number | null;
+    /** Dynamic viscosity [Pa·s] */
+    viscosity: number;
+    /** Thermal conductivity [W/(m·K)] */
+    thermalConductivity: number;
+    /** Surface tension [N/m], only meaningful at saturation */
+    surfaceTension: number;
+    /** Static dielectric constant [dimensionless] */
+    dielectricConstant: number;
+    /** Ionization constant pKw [dimensionless] */
+    ionizationConstant: number;
+    /** Isobaric cubic expansion coefficient [1/K] */
+    isobaricExpansion: number;
+    /** Isothermal compressibility [1/MPa] */
+    isothermalCompressibility: number;
+}
+/**
+ * Core thermodynamic properties computed from a region's fundamental equation.
+ * Does not include transport properties (viscosity, conductivity, etc.)
+ */
+export interface BasicProperties {
+    region: Region;
+    pressure: number;
+    temperature: number;
+    specificVolume: number;
+    internalEnergy: number;
+    entropy: number;
+    enthalpy: number;
+    cp: number;
+    cv: number;
+    speedOfSound: number;
+    quality: number | null;
+    isobaricExpansion: number;
+    isothermalCompressibility: number;
+}
+/** Supported thermodynamic input combinations */
+export type InputMode = 'PT' | 'PH' | 'PS' | 'HS' | 'Px' | 'Tx';
+/** Input for pressure-temperature calculation */
+export interface PTInput {
+    mode: 'PT';
+    /** Pressure [MPa] */
+    p: number;
+    /** Temperature [K] */
+    T: number;
+}
+/** Input for pressure-enthalpy calculation */
+export interface PHInput {
+    mode: 'PH';
+    /** Pressure [MPa] */
+    p: number;
+    /** Specific enthalpy [kJ/kg] */
+    h: number;
+}
+/** Input for pressure-entropy calculation */
+export interface PSInput {
+    mode: 'PS';
+    /** Pressure [MPa] */
+    p: number;
+    /** Specific entropy [kJ/(kg·K)] */
+    s: number;
+}
+/** Input for enthalpy-entropy calculation */
+export interface HSInput {
+    mode: 'HS';
+    /** Specific enthalpy [kJ/kg] */
+    h: number;
+    /** Specific entropy [kJ/(kg·K)] */
+    s: number;
+}
+/** Input for pressure-quality (saturation) calculation */
+export interface PxInput {
+    mode: 'Px';
+    /** Pressure [MPa] */
+    p: number;
+    /** Vapour quality (0–1) */
+    x: number;
+}
+/** Input for temperature-quality (saturation) calculation */
+export interface TxInput {
+    mode: 'Tx';
+    /** Temperature [K] */
+    T: number;
+    /** Vapour quality (0–1) */
+    x: number;
+}
+/** Union of all valid inputs */
+export type SolveInput = PTInput | PHInput | PSInput | HSInput | PxInput | TxInput;
+/**
+ * Typed coefficient table for IAPWS polynomial evaluations.
+ * Each entry is [I, J, N] where:
+ *   - I = pressure exponent
+ *   - J = temperature exponent
+ *   - N = coefficient value
+ */
+export type CoefficientTable = readonly (readonly [number, number, number])[];
+/** Base error class for all IF97 calculation errors */
+export declare class IF97Error extends Error {
+    constructor(message: string);
+}
+/** Thrown when input values are outside IAPWS-IF97 valid ranges */
+export declare class OutOfRangeError extends IF97Error {
+    readonly parameter: string;
+    readonly value: number;
+    readonly min: number;
+    readonly max: number;
+    constructor(parameter: string, value: number, min: number, max: number);
+}
+/** Thrown when a numerical solver fails to converge */
+export declare class ConvergenceError extends IF97Error {
+    readonly solver: string;
+    readonly iterations: number;
+    constructor(solver: string, iterations: number);
+}
+//# sourceMappingURL=types.d.ts.map

package/dist/types.d.ts.map

@@ -0,0 +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,CAAC;IACX,0CAA0C;IAC1C,EAAE,EAAE,MAAM,CAAC;IACX,2BAA2B;IAC3B,YAAY,EAAE,MAAM,CAAC;IACrB,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,CAAC;IACvB,iDAAiD;IACjD,kBAAkB,EAAE,MAAM,CAAC;IAC3B,8CAA8C;IAC9C,kBAAkB,EAAE,MAAM,CAAC;IAI3B,iDAAiD;IACjD,iBAAiB,EAAE,MAAM,CAAC;IAC1B,yCAAyC;IACzC,yBAAyB,EAAE,MAAM,CAAC;CACnC;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,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,YAAY,EAAE,MAAM,CAAC;IACrB,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,yBAAyB,EAAE,MAAM,CAAC;CACnC;AAID,iDAAiD;AACjD,MAAM,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;AAEhE,iDAAiD;AACjD,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,qBAAqB;IACrB,CAAC,EAAE,MAAM,CAAC;IACV,sBAAsB;IACtB,CAAC,EAAE,MAAM,CAAC;CACX;AAED,8CAA8C;AAC9C,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,qBAAqB;IACrB,CAAC,EAAE,MAAM,CAAC;IACV,gCAAgC;IAChC,CAAC,EAAE,MAAM,CAAC;CACX;AAED,6CAA6C;AAC7C,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,qBAAqB;IACrB,CAAC,EAAE,MAAM,CAAC;IACV,mCAAmC;IACnC,CAAC,EAAE,MAAM,CAAC;CACX;AAED,6CAA6C;AAC7C,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,gCAAgC;IAChC,CAAC,EAAE,MAAM,CAAC;IACV,mCAAmC;IACnC,CAAC,EAAE,MAAM,CAAC;CACX;AAED,0DAA0D;AAC1D,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,qBAAqB;IACrB,CAAC,EAAE,MAAM,CAAC;IACV,2BAA2B;IAC3B,CAAC,EAAE,MAAM,CAAC;CACX;AAED,6DAA6D;AAC7D,MAAM,WAAW,OAAO;IACtB,IAAI,EAAE,IAAI,CAAC;IACX,sBAAsB;IACtB,CAAC,EAAE,MAAM,CAAC;IACV,2BAA2B;IAC3B,CAAC,EAAE,MAAM,CAAC;CACX;AAED,gCAAgC;AAChC,MAAM,MAAM,UAAU,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;AAInF;;;;;;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

@@ -0,0 +1,56 @@
+/**
+ * IAPWS-IF97 Core Types & Interfaces
+ *
+ * Type definitions for the IAPWS Industrial Formulation 1997
+ * for the Thermodynamic Properties of Water and Steam.
+ */
+// ─── Region Enum ────────────────────────────────────────────────────────────
+/** IAPWS-IF97 thermodynamic regions */
+export var Region;
+(function (Region) {
+    /** Subcooled liquid (compressed water) */
+    Region[Region["Region1"] = 1] = "Region1";
+    /** Superheated steam */
+    Region[Region["Region2"] = 2] = "Region2";
+    /** Supercritical fluid */
+    Region[Region["Region3"] = 3] = "Region3";
+    /** Two-phase (wet steam / saturation line) */
+    Region[Region["Region4"] = 4] = "Region4";
+    /** High-temperature steam (T > 1073.15 K) */
+    Region[Region["Region5"] = 5] = "Region5";
+})(Region || (Region = {}));
+// ─── Errors ─────────────────────────────────────────────────────────────────
+/** Base error class for all IF97 calculation errors */
+export class IF97Error extends Error {
+    constructor(message) {
+        super(message);
+        this.name = 'IF97Error';
+    }
+}
+/** Thrown when input values are outside IAPWS-IF97 valid ranges */
+export class OutOfRangeError extends IF97Error {
+    parameter;
+    value;
+    min;
+    max;
+    constructor(parameter, value, min, max) {
+        super(`${parameter} = ${value} is out of range [${min}, ${max}]`);
+        this.parameter = parameter;
+        this.value = value;
+        this.min = min;
+        this.max = max;
+        this.name = 'OutOfRangeError';
+    }
+}
+/** Thrown when a numerical solver fails to converge */
+export class ConvergenceError extends IF97Error {
+    solver;
+    iterations;
+    constructor(solver, iterations) {
+        super(`${solver} failed to converge after ${iterations} iterations`);
+        this.solver = solver;
+        this.iterations = iterations;
+        this.name = 'ConvergenceError';
+    }
+}
+//# sourceMappingURL=types.js.map

package/dist/types.js.map

@@ -0,0 +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;AAsJD,+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

@@ -0,0 +1,62 @@
+{
+  "name": "iapws-if97",
+  "version": "1.0.0",
+  "description": "TypeScript implementation of IAPWS-IF97 for water and steam properties, including forward/backward solvers, saturation states, and transport properties.",
+  "author": "jltonghui",
+  "repository": {
+    "type": "git",
+    "url": "git+https://github.com/jltonghui/iapws-if97.git"
+  },
+  "bugs": {
+    "url": "https://github.com/jltonghui/iapws-if97/issues"
+  },
+  "homepage": "https://github.com/jltonghui/iapws-if97#readme",
+  "type": "module",
+  "main": "./dist/index.js",
+  "types": "./dist/index.d.ts",
+  "exports": {
+    ".": {
+      "import": "./dist/index.js",
+      "types": "./dist/index.d.ts"
+    },
+    "./package.json": "./package.json"
+  },
+  "files": [
+    "dist",
+    "LICENSE",
+    "README.md"
+  ],
+  "engines": {
+    "node": ">=18"
+  },
+  "publishConfig": {
+    "access": "public"
+  },
+  "sideEffects": false,
+  "scripts": {
+    "build": "tsc",
+    "test": "vitest run",
+    "test:watch": "vitest",
+    "test:coverage": "vitest run --coverage",
+    "test:standards": "vitest run tests/standards-asme tests/standards-cn",
+    "test:standards:asme": "vitest run tests/standards-asme",
+    "test:standards:cn": "vitest run tests/standards-cn",
+    "prepack": "npm run build && npm test",
+    "prepublishOnly": "npm run build && npm test"
+  },
+  "devDependencies": {
+    "@vitest/coverage-v8": "^4.1.0",
+    "typescript": "^5.4.0",
+    "vitest": "^4.1.0"
+  },
+  "keywords": [
+    "iapws",
+    "if97",
+    "steam",
+    "water",
+    "thermodynamics",
+    "properties",
+    "engineering"
+  ],
+  "license": "MIT"
+}