iapws-if97 1.0.0

package/dist/core/solver.js

@@ -0,0 +1,92 @@
+import { Region, IF97Error, OutOfRangeError } from '../types.js';
+import * as C from '../constants.js';
+import { region1 } from '../regions/region1.js';
+import { region2 } from '../regions/region2.js';
+import { region3ByRhoT } from '../regions/region3.js';
+import { region5 } from '../regions/region5.js';
+import { region3Volume } from '../regions/region3-subregions.js';
+import { detectRegionPT } from './region-detector.js';
+import { newtonRaphson } from '../solvers/newton-raphson.js';
+import { viscosity, thermalConductivity, surfaceTension, dielectricConstant, ionizationConstant } from '../transport/properties.js';
+import { solvePH } from '../backward/ph.js';
+import { solvePS } from '../backward/ps.js';
+import { solveHS } from '../backward/hs.js';
+import { solvePx, solveTx } from '../saturation/two-phase.js';
+/**
+ * Compute full thermodynamic state from basic properties by adding transport properties.
+ */
+function enrichState(basic) {
+    const rho = 1 / basic.specificVolume;
+    const mu = viscosity(basic.temperature, rho);
+    // drhodP_T = ρ · κ_T where κ_T is isothermal compressibility [1/MPa]
+    const drhodP_T = rho * basic.isothermalCompressibility;
+    return {
+        ...basic,
+        density: rho,
+        viscosity: mu,
+        thermalConductivity: thermalConductivity(basic.temperature, rho, basic.cp, basic.cv, drhodP_T, mu),
+        surfaceTension: surfaceTension(basic.temperature),
+        dielectricConstant: dielectricConstant(basic.temperature, rho),
+        ionizationConstant: ionizationConstant(basic.temperature, rho),
+    };
+}
+/**
+ * Solve Region 3 from P and T by first computing density via Newton-Raphson.
+ */
+function solveRegion3PT(p, T) {
+    if (p === C.Pc && T === C.Tc) {
+        return region3ByRhoT(C.RHOc, T);
+    }
+    const v0 = region3Volume(p, T);
+    const rho = newtonRaphson((rho_x) => region3ByRhoT(rho_x, T).pressure - p, 1 / v0);
+    return region3ByRhoT(rho, T);
+}
+/**
+ * Solve for full thermodynamic state given pressure and temperature.
+ *
+ * @param p - Pressure [MPa]
+ * @param T - Temperature [K]
+ * @returns Complete steam state with transport properties
+ */
+export function solvePT(p, T) {
+    if (p < C.P_MIN || p > C.P_MAX) {
+        throw new OutOfRangeError('Pressure', p, C.P_MIN, C.P_MAX);
+    }
+    if (T < C.T_MIN || T > C.T_MAX) {
+        throw new OutOfRangeError('Temperature', T, C.T_MIN, C.T_MAX);
+    }
+    const region = detectRegionPT(p, T);
+    let basic;
+    switch (region) {
+        case Region.Region1:
+            basic = region1(p, T);
+            break;
+        case Region.Region2:
+            basic = region2(p, T);
+            break;
+        case Region.Region3:
+            basic = solveRegion3PT(p, T);
+            break;
+        case Region.Region5:
+            basic = region5(p, T);
+            break;
+        default:
+            throw new IF97Error(`Cannot determine region for P=${p} MPa, T=${T} K`);
+    }
+    return enrichState(basic);
+}
+/**
+ * Unified solver: compute full thermodynamic state from any supported input pair.
+ * Routes to the correct solver based on input mode, enriches with transport properties.
+ */
+export function solve(input) {
+    switch (input.mode) {
+        case 'PT': return solvePT(input.p, input.T);
+        case 'PH': return enrichState(solvePH(input.p, input.h));
+        case 'PS': return enrichState(solvePS(input.p, input.s));
+        case 'HS': return enrichState(solveHS(input.h, input.s));
+        case 'Px': return enrichState(solvePx(input.p, input.x));
+        case 'Tx': return enrichState(solveTx(input.T, input.x));
+    }
+}
+//# sourceMappingURL=solver.js.map

package/dist/core/solver.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"solver.js","sourceRoot":"","sources":["../../src/core/solver.ts"],"names":[],"mappings":"AAKA,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AACjE,OAAO,KAAK,CAAC,MAAM,iBAAiB,CAAC;AACrC,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,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,kCAAkC,CAAC;AACjE,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAC;AACtD,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAC7D,OAAO,EAAE,SAAS,EAAE,mBAAmB,EAAE,cAAc,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,MAAM,4BAA4B,CAAC;AACpI,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AAE9D;;GAEG;AACH,SAAS,WAAW,CAAC,KAAsB;IACzC,MAAM,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC;IACrC,MAAM,EAAE,GAAG,SAAS,CAAC,KAAK,CAAC,WAAW,EAAE,GAAG,CAAC,CAAC;IAC7C,qEAAqE;IACrE,MAAM,QAAQ,GAAG,GAAG,GAAG,KAAK,CAAC,yBAAyB,CAAC;IACvD,OAAO;QACL,GAAG,KAAK;QACR,OAAO,EAAE,GAAG;QACZ,SAAS,EAAE,EAAE;QACb,mBAAmB,EAAE,mBAAmB,CACtC,KAAK,CAAC,WAAW,EAAE,GAAG,EACtB,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,EAAE,EAAE,QAAQ,EAAE,EAAE,CACjC;QACD,cAAc,EAAE,cAAc,CAAC,KAAK,CAAC,WAAW,CAAC;QACjD,kBAAkB,EAAE,kBAAkB,CAAC,KAAK,CAAC,WAAW,EAAE,GAAG,CAAC;QAC9D,kBAAkB,EAAE,kBAAkB,CAAC,KAAK,CAAC,WAAW,EAAE,GAAG,CAAC;KAC/D,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;QAC7B,OAAO,aAAa,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,MAAM,EAAE,GAAG,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/B,MAAM,GAAG,GAAG,aAAa,CACvB,CAAC,KAAa,EAAE,EAAE,CAAC,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,EACvD,CAAC,GAAG,EAAE,CACP,CAAC;IACF,OAAO,aAAa,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC;QAC/B,MAAM,IAAI,eAAe,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;IAC7D,CAAC;IACD,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC;QAC/B,MAAM,IAAI,eAAe,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;IAChE,CAAC;IAED,MAAM,MAAM,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,IAAI,KAAsB,CAAC;IAE3B,QAAQ,MAAM,EAAE,CAAC;QACf,KAAK,MAAM,CAAC,OAAO;YACjB,KAAK,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtB,MAAM;QACR,KAAK,MAAM,CAAC,OAAO;YACjB,KAAK,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtB,MAAM;QACR,KAAK,MAAM,CAAC,OAAO;YACjB,KAAK,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7B,MAAM;QACR,KAAK,MAAM,CAAC,OAAO;YACjB,KAAK,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtB,MAAM;QACR;YACE,MAAM,IAAI,SAAS,CAAC,iCAAiC,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;IAC5E,CAAC;IAED,OAAO,WAAW,CAAC,KAAK,CAAC,CAAC;AAC5B,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,KAAK,CAAC,KAAiB;IACrC,QAAQ,KAAK,CAAC,IAAI,EAAE,CAAC;QACnB,KAAK,IAAI,CAAC,CAAC,OAAO,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;QAC5C,KAAK,IAAI,CAAC,CAAC,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,KAAK,IAAI,CAAC,CAAC,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,KAAK,IAAI,CAAC,CAAC,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,KAAK,IAAI,CAAC,CAAC,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,KAAK,IAAI,CAAC,CAAC,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC;AACH,CAAC"}

package/dist/index.d.ts

@@ -0,0 +1,24 @@
+/**
+ * iapws-if97 — IAPWS-IF97 Steam/Water Properties Library
+ *
+ * Public API exports.
+ */
+export { solvePT, solve } from './core/solver.js';
+export { solvePH } from './backward/ph.js';
+export { solvePS } from './backward/ps.js';
+export { solveHS } from './backward/hs.js';
+export { solvePx, solveTx } from './saturation/two-phase.js';
+export type { SteamState, BasicProperties, SolveInput } from './types.js';
+export { Region, IF97Error, OutOfRangeError, ConvergenceError } from './types.js';
+export { region1 } from './regions/region1.js';
+export { region2 } from './regions/region2.js';
+export { region3ByRhoT } from './regions/region3.js';
+export { saturationPressure, saturationTemperature } from './regions/region4.js';
+export { region5 } from './regions/region5.js';
+export { boundary23_T_to_P, boundary23_P_to_T } from './regions/boundaries.js';
+export { region3Volume, region3SatVolume } from './regions/region3-subregions.js';
+export { viscosity, thermalConductivity, surfaceTension, dielectricConstant, ionizationConstant } from './transport/properties.js';
+export { detectRegionPT, detectRegionPH, detectRegionPS, detectRegionHS } from './core/region-detector.js';
+export { newtonRaphson } from './solvers/newton-raphson.js';
+export { nelderMead } from './solvers/nelder-mead.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,OAAO,EAAE,OAAO,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAG7D,YAAY,EAAE,UAAU,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AAC1E,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,eAAe,EAAE,gBAAgB,EAAE,MAAM,YAAY,CAAC;AAGlF,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,sBAAsB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC;AACjF,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAC/E,OAAO,EAAE,aAAa,EAAE,gBAAgB,EAAE,MAAM,iCAAiC,CAAC;AAGlF,OAAO,EAAE,SAAS,EAAE,mBAAmB,EAAE,cAAc,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AAGnI,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,2BAA2B,CAAC;AAG3G,OAAO,EAAE,aAAa,EAAE,MAAM,6BAA6B,CAAC;AAC5D,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC"}

package/dist/index.js

@@ -0,0 +1,28 @@
+/**
+ * iapws-if97 — IAPWS-IF97 Steam/Water Properties Library
+ *
+ * Public API exports.
+ */
+// Core solvers
+export { solvePT, solve } from './core/solver.js';
+export { solvePH } from './backward/ph.js';
+export { solvePS } from './backward/ps.js';
+export { solveHS } from './backward/hs.js';
+export { solvePx, solveTx } from './saturation/two-phase.js';
+export { Region, IF97Error, OutOfRangeError, ConvergenceError } from './types.js';
+// Region functions (for advanced usage)
+export { region1 } from './regions/region1.js';
+export { region2 } from './regions/region2.js';
+export { region3ByRhoT } from './regions/region3.js';
+export { saturationPressure, saturationTemperature } from './regions/region4.js';
+export { region5 } from './regions/region5.js';
+export { boundary23_T_to_P, boundary23_P_to_T } from './regions/boundaries.js';
+export { region3Volume, region3SatVolume } from './regions/region3-subregions.js';
+// Transport properties (for advanced usage)
+export { viscosity, thermalConductivity, surfaceTension, dielectricConstant, ionizationConstant } from './transport/properties.js';
+// Region detection
+export { detectRegionPT, detectRegionPH, detectRegionPS, detectRegionHS } from './core/region-detector.js';
+// Solvers
+export { newtonRaphson } from './solvers/newton-raphson.js';
+export { nelderMead } from './solvers/nelder-mead.js';
+//# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,eAAe;AACf,OAAO,EAAE,OAAO,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAC3C,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAI7D,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,eAAe,EAAE,gBAAgB,EAAE,MAAM,YAAY,CAAC;AAElF,wCAAwC;AACxC,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,sBAAsB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC;AACjF,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAC/E,OAAO,EAAE,aAAa,EAAE,gBAAgB,EAAE,MAAM,iCAAiC,CAAC;AAElF,4CAA4C;AAC5C,OAAO,EAAE,SAAS,EAAE,mBAAmB,EAAE,cAAc,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AAEnI,mBAAmB;AACnB,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,2BAA2B,CAAC;AAE3G,UAAU;AACV,OAAO,EAAE,aAAa,EAAE,MAAM,6BAA6B,CAAC;AAC5D,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC"}

package/dist/regions/boundaries.d.ts

@@ -0,0 +1,50 @@
+/**
+ * IAPWS-IF97 Region Boundary Equations
+ *
+ * Boundary equations between regions, including:
+ * - B23: Boundary between Region 2 and Region 3
+ * - B3 subregion boundaries (ab, cd, ef, gh, ij, jk, mn, op, qu, rx, uv, wx)
+ *
+ * Reference: IAPWS-IF97 Section 4 & Supplementary Release on v(P,T) for Region 3
+ */
+/**
+ * B23 boundary: Temperature → Pressure
+ * Equation 5 (pp. 5), IAPWS-IF97
+ *
+ * @param T - Temperature [K]
+ * @returns Pressure [MPa] on the B23 boundary
+ */
+export declare function boundary23_T_to_P(T: number): number;
+/**
+ * B23 boundary: Pressure → Temperature
+ * Equation 6 (pp. 6), IAPWS-IF97
+ *
+ * @param p - Pressure [MPa]
+ * @returns Temperature [K] on the B23 boundary
+ */
+export declare function boundary23_P_to_T(p: number): number;
+/** B3ab boundary: P → T */
+export declare function b3ab(p: number): number;
+/** B3cd boundary: P → T */
+export declare function b3cd(p: number): number;
+/** B3ef boundary: P → T */
+export declare function b3ef(p: number): number;
+/** B3gh boundary: P → T */
+export declare function b3gh(p: number): number;
+/** B3ij boundary: P → T */
+export declare function b3ij(p: number): number;
+/** B3jk boundary: P → T */
+export declare function b3jk(p: number): number;
+/** B3mn boundary: P → T */
+export declare function b3mn(p: number): number;
+/** B3op boundary: P → T */
+export declare function b3op(p: number): number;
+/** B3qu boundary: P → T */
+export declare function b3qu(p: number): number;
+/** B3rx boundary: P → T */
+export declare function b3rx(p: number): number;
+/** B3uv boundary: P → T */
+export declare function b3uv(p: number): number;
+/** B3wx boundary: P → T */
+export declare function b3wx(p: number): number;
+//# sourceMappingURL=boundaries.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"boundaries.d.ts","sourceRoot":"","sources":["../../src/regions/boundaries.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAIH;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEnD;AAED;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEnD;AAwBD,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC;AAED,2BAA2B;AAC3B,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKtC"}

package/dist/regions/boundaries.js

@@ -0,0 +1,98 @@
+/**
+ * IAPWS-IF97 Region Boundary Equations
+ *
+ * Boundary equations between regions, including:
+ * - B23: Boundary between Region 2 and Region 3
+ * - B3 subregion boundaries (ab, cd, ef, gh, ij, jk, mn, op, qu, rx, uv, wx)
+ *
+ * Reference: IAPWS-IF97 Section 4 & Supplementary Release on v(P,T) for Region 3
+ */
+// ─── B23 Boundary (Region 2 ↔ Region 3) ────────────────────────────────────
+/**
+ * B23 boundary: Temperature → Pressure
+ * Equation 5 (pp. 5), IAPWS-IF97
+ *
+ * @param T - Temperature [K]
+ * @returns Pressure [MPa] on the B23 boundary
+ */
+export function boundary23_T_to_P(T) {
+    return 348.05185628969 - 1.1671859879975 * T + 0.0010192970039326 * T * T;
+}
+/**
+ * B23 boundary: Pressure → Temperature
+ * Equation 6 (pp. 6), IAPWS-IF97
+ *
+ * @param p - Pressure [MPa]
+ * @returns Temperature [K] on the B23 boundary
+ */
+export function boundary23_P_to_T(p) {
+    return 572.54459862746 + Math.sqrt((p - 13.91883977887) / 0.0010192970039326);
+}
+// ─── B3 Subregion Boundaries ────────────────────────────────────────────────
+// Supplementary Release on Backward Equations for v(P,T) in Region 3
+/** Helper: evaluate polynomial Σ N[i] * P^I[i] */
+function polyEval(p, I, N) {
+    let result = 0;
+    for (let i = 0; i < I.length; i++) {
+        result += N[i] * Math.pow(p, I[i]);
+    }
+    return result;
+}
+/** Helper: evaluate polynomial with ln(P) as argument */
+function polyEvalLn(p, I, N) {
+    const lnP = Math.log(p);
+    let result = 0;
+    for (let i = 0; i < I.length; i++) {
+        result += N[i] * Math.pow(lnP, I[i]);
+    }
+    return result;
+}
+/** B3ab boundary: P → T */
+export function b3ab(p) {
+    return polyEvalLn(p, [0, 1, 2, -1, -2], [1.54793642129415e3, -1.87661219490113e2, 2.13144632222113e1, -1.91887498864292e3, 9.18419702359447e2]);
+}
+/** B3cd boundary: P → T */
+export function b3cd(p) {
+    return polyEval(p, [0, 1, 2, 3], [5.85276966696349e2, 2.78233532206915, -1.27283549295878e-2, 1.59090746562729e-4]);
+}
+/** B3ef boundary: P → T */
+export function b3ef(p) {
+    return 3.727888004 * (p - 22.064) + 647.096;
+}
+/** B3gh boundary: P → T */
+export function b3gh(p) {
+    return polyEval(p, [0, 1, 2, 3, 4], [-2.49284240900418e4, 4.28143584791546e3, -2.6902917314013e2, 7.51608051114157, -7.87105249910383e-2]);
+}
+/** B3ij boundary: P → T */
+export function b3ij(p) {
+    return polyEval(p, [0, 1, 2, 3, 4], [5.84814781649163e2, -6.16179320924617e-1, 2.60763050899562e-1, -5.87071076864459e-3, 5.15308185433082e-5]);
+}
+/** B3jk boundary: P → T */
+export function b3jk(p) {
+    return polyEval(p, [0, 1, 2, 3, 4], [6.17229772068439e2, -7.70600270141675, 6.97072596851896e-1, -1.57391839848015e-2, 1.37897492684194e-4]);
+}
+/** B3mn boundary: P → T */
+export function b3mn(p) {
+    return polyEval(p, [0, 1, 2, 3], [5.35339483742384e2, 7.61978122720128, -1.58365725441648e-1, 1.92871054508108e-3]);
+}
+/** B3op boundary: P → T */
+export function b3op(p) {
+    return polyEvalLn(p, [0, 1, 2, -1, -2], [0.969461372400213e3, -0.332500170441278e3, 0.642859598466067e2, 0.773845935768222e3, -0.152313732937084e4]);
+}
+/** B3qu boundary: P → T */
+export function b3qu(p) {
+    return polyEval(p, [0, 1, 2, 3], [5.65603648239126e2, 5.29062258221222, -1.02020639611016e-1, 1.22240301070145e-3]);
+}
+/** B3rx boundary: P → T */
+export function b3rx(p) {
+    return polyEval(p, [0, 1, 2, 3], [5.84561202520006e2, -1.02961025163669, 2.43293362700452e-1, -2.94905044740799e-3]);
+}
+/** B3uv boundary: P → T */
+export function b3uv(p) {
+    return polyEval(p, [0, 1, 2, 3], [5.28199646263062e2, 8.90579602135307, -0.222814134903755, 0.286791682263697e-2]);
+}
+/** B3wx boundary: P → T */
+export function b3wx(p) {
+    return polyEvalLn(p, [0, 1, 2, -1, -2], [7.2805260914538, 9.73505869861952e1, 1.47370491183191e1, 0.329196213998375e3, 8.73371668682417e2]);
+}
+//# sourceMappingURL=boundaries.js.map

package/dist/regions/boundaries.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"boundaries.js","sourceRoot":"","sources":["../../src/regions/boundaries.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,CAAS;IACzC,OAAO,eAAe,GAAG,eAAe,GAAG,CAAC,GAAG,kBAAkB,GAAG,CAAC,GAAG,CAAC,CAAC;AAC5E,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,CAAS;IACzC,OAAO,eAAe,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,GAAG,kBAAkB,CAAC,CAAC;AAChF,CAAC;AAED,+EAA+E;AAC/E,qEAAqE;AAErE,kDAAkD;AAClD,SAAS,QAAQ,CAAC,CAAS,EAAE,CAAoB,EAAE,CAAoB;IACrE,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,yDAAyD;AACzD,SAAS,UAAU,CAAC,CAAS,EAAE,CAAoB,EAAE,CAAoB;IACvE,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxB,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,UAAU,CAAC,CAAC,EACjB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EACjB,CAAC,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CACvG,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACZ,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,CAAC,CAClF,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,WAAW,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,OAAO,CAAC;AAC9C,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACf,CAAC,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,CAAC,iBAAiB,EAAE,gBAAgB,EAAE,CAAC,mBAAmB,CAAC,CACtG,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACf,CAAC,kBAAkB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,CAAC,CAC3G,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACf,CAAC,kBAAkB,EAAE,CAAC,gBAAgB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,CAAC,CACxG,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACZ,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,CAAC,CAClF,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,UAAU,CAAC,CAAC,EACjB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EACjB,CAAC,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,CAAC,CAC5G,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACZ,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,CAAC,CAClF,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACZ,CAAC,kBAAkB,EAAE,CAAC,gBAAgB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,CAAC,CACnF,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,QAAQ,CAAC,CAAC,EACf,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACZ,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,iBAAiB,EAAE,oBAAoB,CAAC,CACjF,CAAC;AACJ,CAAC;AAED,2BAA2B;AAC3B,MAAM,UAAU,IAAI,CAAC,CAAS;IAC5B,OAAO,UAAU,CAAC,CAAC,EACjB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EACjB,CAAC,eAAe,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,mBAAmB,EAAE,kBAAkB,CAAC,CACnG,CAAC;AACJ,CAAC"}

package/dist/regions/region1.d.ts

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

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

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

package/dist/regions/region1.js

@@ -0,0 +1,95 @@
+/**
+ * IAPWS-IF97 Region 1: Subcooled Liquid
+ *
+ * Basic equation for the specific Gibbs free energy γ(π,τ)
+ * Valid for: 273.15 K ≤ T ≤ 623.15 K, Psat(T) ≤ P ≤ 100 MPa
+ *
+ * Reference: IAPWS-IF97, Section 5 (Equations for Region 1)
+ */
+import { R } from '../constants.js';
+import { Region } from '../types.js';
+// ─── Coefficient Table (Table 2, IAPWS-IF97) ───────────────────────────────
+const COEFFICIENTS = [
+    [0, -2, 0.14632971213167],
+    [0, -1, -0.84548187169114],
+    [0, 0, -3.756360367204],
+    [0, 1, 3.3855169168385],
+    [0, 2, -0.95791963387872],
+    [0, 3, 0.15772038513228],
+    [0, 4, -0.016616417199501],
+    [0, 5, 8.1214629983568e-4],
+    [1, -9, 2.8319080123804e-4],
+    [1, -7, -6.0706301565874e-4],
+    [1, -1, -0.018990068218419],
+    [1, 0, -0.032529748770505],
+    [1, 1, -0.021841717175414],
+    [1, 3, -5.283835796993e-5],
+    [2, -3, -4.7184321073267e-4],
+    [2, 0, -3.0001780793026e-4],
+    [2, 1, 4.7661393906987e-5],
+    [2, 3, -4.4141845330846e-6],
+    [2, 17, -7.2694996297594e-16],
+    [3, -4, -3.1679644845054e-5],
+    [3, 0, -2.8270797985312e-6],
+    [3, 6, -8.5205128120103e-10],
+    [4, -5, -2.2425281908e-6],
+    [4, -2, -6.5171222895601e-7],
+    [4, 10, -1.4341729937924e-13],
+    [5, -8, -4.0516996860117e-7],
+    [8, -11, -1.2734301741641e-9],
+    [8, -6, -1.7424871230634e-10],
+    [21, -29, -6.8762131295531e-19],
+    [23, -31, 1.4478307828521e-20],
+    [29, -38, 2.6335781662795e-23],
+    [30, -39, -1.1947622640071e-23],
+    [31, -40, 1.8228094581404e-24],
+    [32, -41, -9.3537087292458e-26],
+];
+// Reduction constants for Region 1
+const PI_STAR = 16.53; // P* [MPa]
+const TAU_STAR = 1386.0; // T* [K]
+/**
+ * Compute Region 1 thermodynamic properties from pressure and temperature.
+ *
+ * @param p - Pressure [MPa]
+ * @param T - Temperature [K]
+ * @returns Basic thermodynamic properties
+ */
+export function region1(p, T) {
+    const pi = p / PI_STAR;
+    const tau = TAU_STAR / T;
+    // Reduced Gibbs free energy and its derivatives
+    let g = 0;
+    let g_pi = 0;
+    let g_pipi = 0;
+    let g_tau = 0;
+    let g_tautau = 0;
+    let g_pitau = 0;
+    for (const [I, J, N] of COEFFICIENTS) {
+        const piTerm = Math.pow(7.1 - pi, I);
+        const tauTerm = Math.pow(tau - 1.222, J);
+        g += N * piTerm * tauTerm;
+        g_pi += -N * I * Math.pow(7.1 - pi, I - 1) * tauTerm;
+        g_pipi += N * I * (I - 1) * Math.pow(7.1 - pi, I - 2) * tauTerm;
+        g_tau += N * piTerm * J * Math.pow(tau - 1.222, J - 1);
+        g_tautau += N * piTerm * J * (J - 1) * Math.pow(tau - 1.222, J - 2);
+        g_pitau += -N * I * Math.pow(7.1 - pi, I - 1) * J * Math.pow(tau - 1.222, J - 1);
+    }
+    return {
+        region: Region.Region1,
+        pressure: p,
+        temperature: T,
+        specificVolume: R * T * pi * g_pi / (1000 * p),
+        internalEnergy: R * T * (tau * g_tau - pi * g_pi),
+        entropy: R * (tau * g_tau - g),
+        enthalpy: R * T * tau * g_tau,
+        cp: R * (-tau * tau * g_tautau),
+        cv: R * (-tau * tau * g_tautau + Math.pow(g_pi - tau * g_pitau, 2) / g_pipi),
+        speedOfSound: Math.sqrt(1000 * R * T * g_pi * g_pi /
+            (Math.pow(g_pi - tau * g_pitau, 2) / (tau * tau * g_tautau) - g_pipi)),
+        quality: null,
+        isobaricExpansion: (1 - tau * g_pitau / g_pi) / T,
+        isothermalCompressibility: -pi * g_pipi / g_pi / p,
+    };
+}
+//# sourceMappingURL=region1.js.map

package/dist/regions/region1.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region1.js","sourceRoot":"","sources":["../../src/regions/region1.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,CAAC,EAAE,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,gBAAgB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,cAAc,CAAC;IACvB,CAAC,CAAC,EAAE,CAAC,EAAE,eAAe,CAAC;IACvB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC;IACxB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC5B,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,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,eAAe,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC/B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC/B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC;CACvB,CAAC;AAEX,mCAAmC;AACnC,MAAM,OAAO,GAAG,KAAK,CAAC,CAAE,WAAW;AACnC,MAAM,QAAQ,GAAG,MAAM,CAAC,CAAC,SAAS;AAElC;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,MAAM,EAAE,GAAG,CAAC,GAAG,OAAO,CAAC;IACvB,MAAM,GAAG,GAAG,QAAQ,GAAG,CAAC,CAAC;IAEzB,gDAAgD;IAChD,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,IAAI,OAAO,GAAG,CAAC,CAAC;IAEhB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,YAAY,EAAE,CAAC;QACrC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;QACrC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,EAAE,CAAC,CAAC,CAAC;QAEzC,CAAC,IAAI,CAAC,GAAG,MAAM,GAAG,OAAO,CAAC;QAC1B,IAAI,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;QACrD,MAAM,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;QAChE,KAAK,IAAI,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACvD,QAAQ,IAAI,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACpE,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACnF,CAAC;IAED,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,OAAO;QACtB,QAAQ,EAAE,CAAC;QACX,WAAW,EAAE,CAAC;QACd,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC;QAC9C,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,GAAG,EAAE,GAAG,IAAI,CAAC;QACjD,OAAO,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,GAAG,CAAC,CAAC;QAC9B,QAAQ,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,KAAK;QAC7B,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,QAAQ,CAAC;QAC/B,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,OAAO,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC;QAC5E,YAAY,EAAE,IAAI,CAAC,IAAI,CACrB,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,IAAI;YAC1B,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,OAAO,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,QAAQ,CAAC,GAAG,MAAM,CAAC,CACtE;QACD,OAAO,EAAE,IAAI;QACb,iBAAiB,EAAE,CAAC,CAAC,GAAG,GAAG,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC;QACjD,yBAAyB,EAAE,CAAC,EAAE,GAAG,MAAM,GAAG,IAAI,GAAG,CAAC;KACnD,CAAC;AACJ,CAAC"}

package/dist/regions/region2.d.ts

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

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

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

package/dist/regions/region2.js

@@ -0,0 +1,134 @@
+/**
+ * IAPWS-IF97 Region 2: Superheated Steam
+ *
+ * Basic equation for the specific Gibbs free energy γ(π,τ) = γ°(π,τ) + γʳ(π,τ)
+ * Valid for: 273.15 K ≤ T ≤ 623.15 K at Psat(T), and T ≤ 1073.15 K at P ≤ 100 MPa
+ *
+ * Reference: IAPWS-IF97, Section 6 (Equations for Region 2)
+ */
+import { R } from '../constants.js';
+import { Region } from '../types.js';
+// ─── Ideal Part Coefficients (Table 10) ─────────────────────────────────────
+/** γ° coefficients: [_, J, N] (I is implicit for ideal part) */
+const IDEAL_COEFFICIENTS = [
+    [0, -9.6927686500217],
+    [1, 10.086655968018],
+    [-5, -0.005608791128302],
+    [-4, 0.071452738081455],
+    [-3, -0.40710498223928],
+    [-2, 1.4240819171444],
+    [-1, -4.383951131945],
+    [2, -0.28408632460772],
+    [3, 0.021268463753307],
+];
+// ─── Residual Part Coefficients (Table 11) ──────────────────────────────────
+const RESIDUAL_COEFFICIENTS = [
+    [1, 0, -1.7731742473213e-3],
+    [1, 1, -0.017834862292358],
+    [1, 2, -0.045996013696365],
+    [1, 3, -0.057581259083432],
+    [1, 6, -0.05032527872793],
+    [2, 1, -3.3032641670203e-5],
+    [2, 2, -1.8948987516315e-4],
+    [2, 4, -3.9392777243355e-3],
+    [2, 7, -0.043797295650573],
+    [2, 36, -2.6674547914087e-5],
+    [3, 0, 2.0481737692309e-8],
+    [3, 1, 4.3870667284435e-7],
+    [3, 3, -3.227767723857e-5],
+    [3, 6, -1.5033924542148e-3],
+    [3, 35, -0.040668253562649],
+    [4, 1, -7.8847309559367e-10],
+    [4, 2, 1.2790717852285e-8],
+    [4, 3, 4.8225372718507e-7],
+    [5, 7, 2.2922076337661e-6],
+    [6, 3, -1.6714766451061e-11],
+    [6, 16, -2.1171472321355e-3],
+    [6, 35, -23.895741934104],
+    [7, 0, -5.905956432427e-18],
+    [7, 11, -1.2621808899101e-6],
+    [7, 25, -0.038946842435739],
+    [8, 8, 1.1256211360459e-11],
+    [8, 36, -8.2311340897998],
+    [9, 13, 1.9809712802088e-8],
+    [10, 4, 1.0406965210174e-19],
+    [10, 10, -1.0234747095929e-13],
+    [10, 14, -1.0018179379511e-9],
+    [16, 29, -8.0882908646985e-11],
+    [16, 50, 0.10693031879409],
+    [18, 57, -0.33662250574171],
+    [20, 20, 8.9185845355421e-25],
+    [20, 35, 3.0629316876232e-13],
+    [20, 48, -4.2002467698208e-6],
+    [21, 21, -5.9056029685639e-26],
+    [22, 53, 3.7826947613457e-6],
+    [23, 39, -1.2768608934681e-15],
+    [24, 26, 7.3087610595061e-29],
+    [24, 40, 5.5414715350778e-17],
+    [24, 58, -9.436970724121e-7],
+];
+/**
+ * Compute Region 2 thermodynamic properties from pressure and temperature.
+ *
+ * @param p - Pressure [MPa]
+ * @param T - Temperature [K]
+ * @returns Basic thermodynamic properties
+ */
+export function region2(p, T) {
+    const pi = p;
+    const tau = 540.0 / T;
+    // ── Ideal part derivatives ──
+    let gi = Math.log(pi);
+    let gi_pi = 1 / pi;
+    // gi_pipi = -1/(pi*pi) — not needed for Region 2 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 derivatives ──
+    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 tauShift = tau - 0.5;
+        const tauPow = Math.pow(tauShift, 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(tauShift, J - 1);
+        gr_tautau += N * piPow * J * (J - 1) * Math.pow(tauShift, J - 2);
+        gr_pitau += N * I * Math.pow(pi, I - 1) * J * Math.pow(tauShift, J - 1);
+    }
+    const gp = gi_pi + gr_pi;
+    const gt = gi_tau + gr_tau;
+    const gtt = gi_tautau + gr_tautau;
+    return {
+        region: Region.Region2,
+        pressure: p,
+        temperature: T,
+        specificVolume: (R * T / (1000 * p)) * pi * gp,
+        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=region2.js.map

package/dist/regions/region2.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"region2.js","sourceRoot":"","sources":["../../src/regions/region2.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,gEAAgE;AAChE,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,iBAAiB,CAAC;IACvB,CAAC,CAAC,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACvB,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC;IACrB,CAAC,CAAC,CAAC,EAAE,CAAC,cAAc,CAAC;IACrB,CAAC,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACtB,CAAC,CAAC,EAAE,iBAAiB,CAAC;CACd,CAAC;AAEX,+EAA+E;AAE/E,MAAM,qBAAqB,GAAqB;IAC9C,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,gBAAgB,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,iBAAiB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,iBAAiB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,CAAC;IAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,mBAAmB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,eAAe,CAAC;IACzB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,kBAAkB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC5B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,iBAAiB,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC3B,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,eAAe,CAAC;IACzB,CAAC,CAAC,EAAE,EAAE,EAAE,kBAAkB,CAAC;IAC3B,CAAC,EAAE,EAAE,CAAC,EAAE,mBAAmB,CAAC;IAC5B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,EAAE,EAAE,gBAAgB,CAAC;IAC1B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,gBAAgB,CAAC;IAC3B,CAAC,EAAE,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,kBAAkB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,EAAE,EAAE,kBAAkB,CAAC;IAC5B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,mBAAmB,CAAC;IAC9B,CAAC,EAAE,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,mBAAmB,CAAC;IAC7B,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,iBAAiB,CAAC;CACpB,CAAC;AAEX;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,MAAM,EAAE,GAAG,CAAC,CAAC;IACb,MAAM,GAAG,GAAG,KAAK,GAAG,CAAC,CAAC;IAEtB,+BAA+B;IAC/B,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,kCAAkC;IAClC,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,QAAQ,GAAG,GAAG,GAAG,GAAG,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC;QAErC,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,QAAQ,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACpD,SAAS,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACjE,QAAQ,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC1E,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,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE;QAC9C,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/regions/region3-data-am.d.ts

@@ -0,0 +1,19 @@
+/**
+ * Region 3 subregions a–m coefficient data.
+ * Reference: Supplementary Release on Backward Equations for v(P,T) in Region 3
+ */
+import type { SubregionConfig } from './region3-eval.js';
+export declare const R3A: SubregionConfig;
+export declare const R3B: SubregionConfig;
+export declare const R3C: SubregionConfig;
+export declare const R3D: SubregionConfig;
+export declare const R3E: SubregionConfig;
+export declare const R3F: SubregionConfig;
+export declare const R3G: SubregionConfig;
+export declare const R3H: SubregionConfig;
+export declare const R3I: SubregionConfig;
+export declare const R3J: SubregionConfig;
+export declare const R3K: SubregionConfig;
+export declare const R3L: SubregionConfig;
+export declare const R3M: SubregionConfig;
+//# sourceMappingURL=region3-data-am.d.ts.map

package/dist/regions/region3-data-am.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"region3-data-am.d.ts","sourceRoot":"","sources":["../../src/regions/region3-data-am.ts"],"names":[],"mappings":"AAAA;;;GAGG;AACH,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEzD,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC;AACF,eAAO,MAAM,GAAG,EAAE,eAGjB,CAAC"}