formulab 0.5.1 → 0.10.1

package/dist/battery/cRate.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"cRate.d.ts","sourceRoot":"","sources":["../../src/battery/cRate.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAE1D;;;;;;GAMG;AACH,wBAAgB,KAAK,CAAC,KAAK,EAAE,UAAU,GAAG,WAAW,CAqBpD"}
+{"version":3,"file":"cRate.d.ts","sourceRoot":"","sources":["../../src/battery/cRate.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAE1D;;;;;;GAMG;AACH,wBAAgB,KAAK,CAAC,KAAK,EAAE,UAAU,GAAG,WAAW,CAyBpD"}

package/dist/battery/cRate.js

@@ -7,13 +7,19 @@ import { roundTo } from '../utils.js';
  * @returns C-rate, current, and theoretical charge/discharge time
  */
 export function cRate(input) {
+    if (input.capacityAh <= 0)
+        throw new RangeError('capacityAh must be positive');
     let cRateValue;
     let currentA;
     if (input.mode === 'currentToRate') {
+        if (input.currentA <= 0)
+            throw new RangeError('currentA must be positive');
         currentA = input.currentA;
         cRateValue = currentA / input.capacityAh;
     }
     else {
+        if (input.cRate <= 0)
+            throw new RangeError('cRate must be positive');
         cRateValue = input.cRate;
         currentA = cRateValue * input.capacityAh;
     }

package/dist/battery/cRate.js.map

@@ -1 +1 @@
-{"version":3,"file":"cRate.js","sourceRoot":"","sources":["../../src/battery/cRate.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;GAMG;AACH,MAAM,UAAU,KAAK,CAAC,KAAiB;IACrC,IAAI,UAAkB,CAAC;IACvB,IAAI,QAAgB,CAAC;IAErB,IAAI,KAAK,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;QACnC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,UAAU,GAAG,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;IAC3C,CAAC;SAAM,CAAC;QACN,UAAU,GAAG,KAAK,CAAC,KAAK,CAAC;QACzB,QAAQ,GAAG,UAAU,GAAG,KAAK,CAAC,UAAU,CAAC;IAC3C,CAAC;IAED,MAAM,gBAAgB,GAAG,CAAC,GAAG,UAAU,CAAC;IACxC,MAAM,kBAAkB,GAAG,gBAAgB,GAAG,EAAE,CAAC;IAEjD,OAAO;QACL,KAAK,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;QAC7B,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,gBAAgB,EAAE,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC;QAC9C,kBAAkB,EAAE,OAAO,CAAC,kBAAkB,EAAE,CAAC,CAAC;KACnD,CAAC;AACJ,CAAC"}
+{"version":3,"file":"cRate.js","sourceRoot":"","sources":["../../src/battery/cRate.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;GAMG;AACH,MAAM,UAAU,KAAK,CAAC,KAAiB;IACrC,IAAI,KAAK,CAAC,UAAU,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,CAAC,CAAC;IAE/E,IAAI,UAAkB,CAAC;IACvB,IAAI,QAAgB,CAAC;IAErB,IAAI,KAAK,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;QACnC,IAAI,KAAK,CAAC,QAAQ,IAAI,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,2BAA2B,CAAC,CAAC;QAC3E,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,UAAU,GAAG,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;IAC3C,CAAC;SAAM,CAAC;QACN,IAAI,KAAK,CAAC,KAAK,IAAI,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,wBAAwB,CAAC,CAAC;QACrE,UAAU,GAAG,KAAK,CAAC,KAAK,CAAC;QACzB,QAAQ,GAAG,UAAU,GAAG,KAAK,CAAC,UAAU,CAAC;IAC3C,CAAC;IAED,MAAM,gBAAgB,GAAG,CAAC,GAAG,UAAU,CAAC;IACxC,MAAM,kBAAkB,GAAG,gBAAgB,GAAG,EAAE,CAAC;IAEjD,OAAO;QACL,KAAK,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;QAC7B,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,gBAAgB,EAAE,OAAO,CAAC,gBAAgB,EAAE,CAAC,CAAC;QAC9C,kBAAkB,EAAE,OAAO,CAAC,kBAAkB,EAAE,CAAC,CAAC;KACnD,CAAC;AACJ,CAAC"}

package/dist/battery/guards.d.ts

@@ -0,0 +1,7 @@
+import type { CRateInput } from './types.js';
+/**
+ * Runtime type guard for {@link CRateInput} discriminated union.
+ * Validates the `mode` discriminant and required fields for each variant.
+ */
+export declare function isCRateInput(input: unknown): input is CRateInput;
+//# sourceMappingURL=guards.d.ts.map

package/dist/battery/guards.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"guards.d.ts","sourceRoot":"","sources":["../../src/battery/guards.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AAE7C;;;GAGG;AACH,wBAAgB,YAAY,CAAC,KAAK,EAAE,OAAO,GAAG,KAAK,IAAI,UAAU,CAShE"}

package/dist/battery/guards.js

@@ -0,0 +1,17 @@
+/**
+ * Runtime type guard for {@link CRateInput} discriminated union.
+ * Validates the `mode` discriminant and required fields for each variant.
+ */
+export function isCRateInput(input) {
+    if (typeof input !== 'object' || input === null)
+        return false;
+    const o = input;
+    if (typeof o.capacityAh !== 'number')
+        return false;
+    switch (o.mode) {
+        case 'currentToRate': return typeof o.currentA === 'number';
+        case 'rateToCurrent': return typeof o.cRate === 'number';
+        default: return false;
+    }
+}
+//# sourceMappingURL=guards.js.map

package/dist/battery/guards.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"guards.js","sourceRoot":"","sources":["../../src/battery/guards.ts"],"names":[],"mappings":"AAEA;;;GAGG;AACH,MAAM,UAAU,YAAY,CAAC,KAAc;IACzC,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC;IAC9D,MAAM,CAAC,GAAG,KAAgC,CAAC;IAC3C,IAAI,OAAO,CAAC,CAAC,UAAU,KAAK,QAAQ;QAAE,OAAO,KAAK,CAAC;IACnD,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC;QACf,KAAK,eAAe,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,QAAQ,KAAK,QAAQ,CAAC;QAC5D,KAAK,eAAe,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,KAAK,KAAK,QAAQ,CAAC;QACzD,OAAO,CAAC,CAAC,OAAO,KAAK,CAAC;IACxB,CAAC;AACH,CAAC"}

package/dist/battery/index.d.ts

@@ -1,4 +1,5 @@
 export type { BatteryChemistry, EnergyDensityInput, EnergyDensityResult, CRateMode, CRateInput, CRateResult, StateOfHealthInput, SohStatus, StateOfHealthResult, BatteryPackConfigInput, BatteryPackConfigResult, CycleLifeInput, CycleLifeResult, InternalResistanceInput, InternalResistanceResult, SelfDischargeInput, SelfDischargeResult, ThermalRunawayInput, ThermalRunawayResult, BmsBalancingInput, BmsBalancingResult, ChargingProfileInput, ChargingProfileResult, } from './types.js';
+export { isCRateInput } from './guards.js';
 export { energyDensity } from './energyDensity.js';
 export { cRate } from './cRate.js';
 export { stateOfHealth } from './stateOfHealth.js';

package/dist/battery/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/battery/index.ts"],"names":[],"mappings":"AACA,YAAY,EACV,gBAAgB,EAChB,kBAAkB,EAClB,mBAAmB,EACnB,SAAS,EACT,UAAU,EACV,WAAW,EACX,kBAAkB,EAClB,SAAS,EACT,mBAAmB,EACnB,sBAAsB,EACtB,uBAAuB,EACvB,cAAc,EACd,eAAe,EACf,uBAAuB,EACvB,wBAAwB,EACxB,kBAAkB,EAClB,mBAAmB,EACnB,mBAAmB,EACnB,oBAAoB,EACpB,iBAAiB,EACjB,kBAAkB,EAClB,oBAAoB,EACpB,qBAAqB,GACtB,MAAM,YAAY,CAAC;AAGpB,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/battery/index.ts"],"names":[],"mappings":"AACA,YAAY,EACV,gBAAgB,EAChB,kBAAkB,EAClB,mBAAmB,EACnB,SAAS,EACT,UAAU,EACV,WAAW,EACX,kBAAkB,EAClB,SAAS,EACT,mBAAmB,EACnB,sBAAsB,EACtB,uBAAuB,EACvB,cAAc,EACd,eAAe,EACf,uBAAuB,EACvB,wBAAwB,EACxB,kBAAkB,EAClB,mBAAmB,EACnB,mBAAmB,EACnB,oBAAoB,EACpB,iBAAiB,EACjB,kBAAkB,EAClB,oBAAoB,EACpB,qBAAqB,GACtB,MAAM,YAAY,CAAC;AAGpB,OAAO,EAAE,YAAY,EAAE,MAAM,aAAa,CAAC;AAG3C,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/battery/index.js

@@ -1,3 +1,5 @@
+// Guards
+export { isCRateInput } from './guards.js';
 // Functions
 export { energyDensity } from './energyDensity.js';
 export { cRate } from './cRate.js';

package/dist/battery/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/battery/index.ts"],"names":[],"mappings":"AA2BA,YAAY;AACZ,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/battery/index.ts"],"names":[],"mappings":"AA2BA,SAAS;AACT,OAAO,EAAE,YAAY,EAAE,MAAM,aAAa,CAAC;AAE3C,YAAY;AACZ,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/chemical/dilution.js

@@ -16,7 +16,7 @@ export function dilution(input) {
         case 'c2':
             // C2 = C1 * V1 / V2
             if (input.v2 === 0) {
-                throw new Error('Cannot solve for c2: v2 must be non-zero');
+                throw new RangeError('Cannot solve for c2: v2 must be non-zero');
             }
             resultC1 = input.c1;
             resultV1 = input.v1;
@@ -26,7 +26,7 @@ export function dilution(input) {
         case 'v2':
             // V2 = C1 * V1 / C2
             if (input.c2 === 0) {
-                throw new Error('Cannot solve for v2: c2 must be non-zero');
+                throw new RangeError('Cannot solve for v2: c2 must be non-zero');
             }
             resultC1 = input.c1;
             resultV1 = input.v1;
@@ -36,7 +36,7 @@ export function dilution(input) {
         case 'c1':
             // C1 = C2 * V2 / V1
             if (input.v1 === 0) {
-                throw new Error('Cannot solve for c1: v1 must be non-zero');
+                throw new RangeError('Cannot solve for c1: v1 must be non-zero');
             }
             resultV1 = input.v1;
             resultC2 = input.c2;
@@ -46,7 +46,7 @@ export function dilution(input) {
         case 'v1':
             // V1 = C2 * V2 / C1
             if (input.c1 === 0) {
-                throw new Error('Cannot solve for v1: c1 must be non-zero');
+                throw new RangeError('Cannot solve for v1: c1 must be non-zero');
             }
             resultC1 = input.c1;
             resultC2 = input.c2;

package/dist/chemical/dilution.js.map

@@ -1 +1 @@
-{"version":3,"file":"dilution.js","sourceRoot":"","sources":["../../src/chemical/dilution.ts"],"names":[],"mappings":"AAEA;;;;;;;;GAQG;AACH,MAAM,UAAU,QAAQ,CAAC,KAAoB;IAC3C,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IAErB,QAAQ,KAAK,CAAC,QAAQ,EAAE,CAAC;QACvB,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;IACV,CAAC;IAED,OAAO;QACL,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,YAAY,EAAE,QAAQ,GAAG,QAAQ;KAClC,CAAC;AACJ,CAAC"}
+{"version":3,"file":"dilution.js","sourceRoot":"","sources":["../../src/chemical/dilution.ts"],"names":[],"mappings":"AAEA;;;;;;;;GAQG;AACH,MAAM,UAAU,QAAQ,CAAC,KAAoB;IAC3C,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IAErB,QAAQ,KAAK,CAAC,QAAQ,EAAE,CAAC;QACvB,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;YACnE,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;YACnE,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;YACnE,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;gBACnB,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;YACnE,CAAC;YACD,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,QAAQ,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YAC5C,MAAM;IACV,CAAC;IAED,OAAO;QACL,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;QACZ,YAAY,EAAE,QAAQ,GAAG,QAAQ;KAClC,CAAC;AACJ,CAAC"}

package/dist/chemical/flowControl.d.ts

@@ -0,0 +1,14 @@
+import type { FlowControlInput, FlowControlResult } from './types.js';
+/**
+ * Control Valve Cv Calculator
+ *
+ * @formula
+ *   - Liquid: Cv = Q × √(SG / ΔP)  (Q in GPM, ΔP in psi)
+ *   - Gas: Cv = Q / (N × P1) × √(SG × T / ΔP)
+ *   - Kv = 0.865 × Cv
+ *
+ * @reference ISA-75.01.01-2012 — Flow equations for sizing control valves
+ * @reference IEC 60534-2-1 — Industrial-process control valves
+ */
+export declare function flowControl(input: FlowControlInput): FlowControlResult;
+//# sourceMappingURL=flowControl.d.ts.map

package/dist/chemical/flowControl.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"flowControl.d.ts","sourceRoot":"","sources":["../../src/chemical/flowControl.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,gBAAgB,EAAE,iBAAiB,EAAE,MAAM,YAAY,CAAC;AAEtE;;;;;;;;;;GAUG;AACH,wBAAgB,WAAW,CAAC,KAAK,EAAE,gBAAgB,GAAG,iBAAiB,CAkEtE"}

package/dist/chemical/flowControl.js

@@ -0,0 +1,69 @@
+import { roundTo } from '../utils.js';
+/**
+ * Control Valve Cv Calculator
+ *
+ * @formula
+ *   - Liquid: Cv = Q × √(SG / ΔP)  (Q in GPM, ΔP in psi)
+ *   - Gas: Cv = Q / (N × P1) × √(SG × T / ΔP)
+ *   - Kv = 0.865 × Cv
+ *
+ * @reference ISA-75.01.01-2012 — Flow equations for sizing control valves
+ * @reference IEC 60534-2-1 — Industrial-process control valves
+ */
+export function flowControl(input) {
+    const { flowRate, inletPressure, outletPressure, fluidDensity, fluidType, temperature = 20, molecularWeight, specificHeatRatio = 1.4, } = input;
+    const pressureDrop = inletPressure - outletPressure; // kPa
+    // Reference density: water = 999 kg/m³
+    const sg = fluidDensity / 999;
+    const pressureRatio = pressureDrop / inletPressure;
+    let cv;
+    let isChoked = false;
+    if (fluidType === 'liquid') {
+        // ISA liquid sizing: Cv = Q(m³/h) × √(SG / ΔP(bar))
+        // Convert kPa to bar: 1 bar = 100 kPa
+        const dpBar = pressureDrop / 100;
+        if (dpBar <= 0) {
+            cv = 0;
+        }
+        else {
+            cv = flowRate * Math.sqrt(sg / dpBar);
+        }
+    }
+    else {
+        // Gas/Steam sizing
+        // Critical pressure ratio for choked flow
+        const xT = 0.7; // typical for globe valve
+        const Fk = specificHeatRatio / 1.4;
+        const xCritical = Fk * xT;
+        const x = pressureDrop / inletPressure;
+        if (x >= xCritical) {
+            isChoked = true;
+        }
+        const xEffective = Math.min(x, xCritical);
+        const p1Bar = inletPressure / 100;
+        const T = temperature + 273.15;
+        const MW = molecularWeight ?? 29; // default air
+        // N8 = 94.8 for metric (m³/h, bar, K)
+        const N8 = 94.8;
+        const Y = 1 - xEffective / (3 * xCritical); // Expansion factor
+        if (p1Bar <= 0 || xEffective <= 0) {
+            cv = 0;
+        }
+        else {
+            cv = flowRate / (N8 * p1Bar * Y) * Math.sqrt(MW * T / xEffective);
+        }
+    }
+    const kv = cv * 0.865;
+    // Estimate pipe velocity (assuming DN50 pipe)
+    const pipeArea = Math.PI * (0.05 / 2) ** 2; // m²
+    const velocity = (flowRate / 3600) / pipeArea; // m/s
+    return {
+        cv: roundTo(cv, 2),
+        kv: roundTo(kv, 2),
+        pressureDrop: roundTo(pressureDrop, 2),
+        pressureRatio: roundTo(pressureRatio, 4),
+        isChoked,
+        velocity: roundTo(velocity, 2),
+    };
+}
+//# sourceMappingURL=flowControl.js.map

package/dist/chemical/flowControl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"flowControl.js","sourceRoot":"","sources":["../../src/chemical/flowControl.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,WAAW,CAAC,KAAuB;IACjD,MAAM,EACJ,QAAQ,EAAE,aAAa,EAAE,cAAc,EAAE,YAAY,EAAE,SAAS,EAChE,WAAW,GAAG,EAAE,EAAE,eAAe,EAAE,iBAAiB,GAAG,GAAG,GAC3D,GAAG,KAAK,CAAC;IAEV,MAAM,YAAY,GAAG,aAAa,GAAG,cAAc,CAAC,CAAC,MAAM;IAE3D,uCAAuC;IACvC,MAAM,EAAE,GAAG,YAAY,GAAG,GAAG,CAAC;IAC9B,MAAM,aAAa,GAAG,YAAY,GAAG,aAAa,CAAC;IAEnD,IAAI,EAAU,CAAC;IACf,IAAI,QAAQ,GAAG,KAAK,CAAC;IAErB,IAAI,SAAS,KAAK,QAAQ,EAAE,CAAC;QAC3B,oDAAoD;QACpD,sCAAsC;QACtC,MAAM,KAAK,GAAG,YAAY,GAAG,GAAG,CAAC;QACjC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;YACf,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,CAAC;YACN,EAAE,GAAG,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC;QACxC,CAAC;IACH,CAAC;SAAM,CAAC;QACN,mBAAmB;QACnB,0CAA0C;QAC1C,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,0BAA0B;QAC1C,MAAM,EAAE,GAAG,iBAAiB,GAAG,GAAG,CAAC;QACnC,MAAM,SAAS,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1B,MAAM,CAAC,GAAG,YAAY,GAAG,aAAa,CAAC;QAEvC,IAAI,CAAC,IAAI,SAAS,EAAE,CAAC;YACnB,QAAQ,GAAG,IAAI,CAAC;QAClB,CAAC;QAED,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,aAAa,GAAG,GAAG,CAAC;QAClC,MAAM,CAAC,GAAG,WAAW,GAAG,MAAM,CAAC;QAC/B,MAAM,EAAE,GAAG,eAAe,IAAI,EAAE,CAAC,CAAC,cAAc;QAEhD,sCAAsC;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC;QAChB,MAAM,CAAC,GAAG,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,mBAAmB;QAE/D,IAAI,KAAK,IAAI,CAAC,IAAI,UAAU,IAAI,CAAC,EAAE,CAAC;YAClC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,CAAC;YACN,EAAE,GAAG,QAAQ,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC;QACpE,CAAC;IACH,CAAC;IAED,MAAM,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC;IAEtB,8CAA8C;IAC9C,MAAM,QAAQ,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK;IACjD,MAAM,QAAQ,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,GAAG,QAAQ,CAAC,CAAC,MAAM;IAErD,OAAO;QACL,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,YAAY,EAAE,OAAO,CAAC,YAAY,EAAE,CAAC,CAAC;QACtC,aAAa,EAAE,OAAO,CAAC,aAAa,EAAE,CAAC,CAAC;QACxC,QAAQ;QACR,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;KAC/B,CAAC;AACJ,CAAC"}

package/dist/chemical/guards.d.ts

@@ -0,0 +1,17 @@
+import type { DilutionInput, ReactorInput, HeatTransferInput } from './types.js';
+/**
+ * Runtime type guard for {@link DilutionInput} discriminated union.
+ * Validates the `solveFor` discriminant and required fields for each variant.
+ */
+export declare function isDilutionInput(input: unknown): input is DilutionInput;
+/**
+ * Runtime type guard for {@link ReactorInput} discriminated union.
+ * Validates the `shape` discriminant and required fields for each variant.
+ */
+export declare function isReactorInput(input: unknown): input is ReactorInput;
+/**
+ * Runtime type guard for {@link HeatTransferInput} discriminated union.
+ * Validates the `mode` discriminant and required fields for each variant.
+ */
+export declare function isHeatTransferInput(input: unknown): input is HeatTransferInput;
+//# sourceMappingURL=guards.d.ts.map

package/dist/chemical/guards.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"guards.d.ts","sourceRoot":"","sources":["../../src/chemical/guards.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,aAAa,EAAE,YAAY,EAAE,iBAAiB,EAAE,MAAM,YAAY,CAAC;AAEjF;;;GAGG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,OAAO,GAAG,KAAK,IAAI,aAAa,CAUtE;AAED;;;GAGG;AACH,wBAAgB,cAAc,CAAC,KAAK,EAAE,OAAO,GAAG,KAAK,IAAI,YAAY,CASpE;AAED;;;GAGG;AACH,wBAAgB,mBAAmB,CAAC,KAAK,EAAE,OAAO,GAAG,KAAK,IAAI,iBAAiB,CAe9E"}

package/dist/chemical/guards.js

@@ -0,0 +1,54 @@
+/**
+ * Runtime type guard for {@link DilutionInput} discriminated union.
+ * Validates the `solveFor` discriminant and required fields for each variant.
+ */
+export function isDilutionInput(input) {
+    if (typeof input !== 'object' || input === null)
+        return false;
+    const o = input;
+    switch (o.solveFor) {
+        case 'c1': return typeof o.v1 === 'number' && typeof o.c2 === 'number' && typeof o.v2 === 'number';
+        case 'v1': return typeof o.c1 === 'number' && typeof o.c2 === 'number' && typeof o.v2 === 'number';
+        case 'c2': return typeof o.c1 === 'number' && typeof o.v1 === 'number' && typeof o.v2 === 'number';
+        case 'v2': return typeof o.c1 === 'number' && typeof o.v1 === 'number' && typeof o.c2 === 'number';
+        default: return false;
+    }
+}
+/**
+ * Runtime type guard for {@link ReactorInput} discriminated union.
+ * Validates the `shape` discriminant and required fields for each variant.
+ */
+export function isReactorInput(input) {
+    if (typeof input !== 'object' || input === null)
+        return false;
+    const o = input;
+    if (typeof o.diameter !== 'number' || typeof o.fillRatio !== 'number')
+        return false;
+    switch (o.shape) {
+        case 'cylindrical': return typeof o.height === 'number';
+        case 'spherical': return true;
+        default: return false;
+    }
+}
+/**
+ * Runtime type guard for {@link HeatTransferInput} discriminated union.
+ * Validates the `mode` discriminant and required fields for each variant.
+ */
+export function isHeatTransferInput(input) {
+    if (typeof input !== 'object' || input === null)
+        return false;
+    const o = input;
+    switch (o.mode) {
+        case 'conduction':
+            return typeof o.conductivity === 'number' && typeof o.area === 'number' &&
+                typeof o.thickness === 'number' && typeof o.tempHot === 'number' && typeof o.tempCold === 'number';
+        case 'convection':
+            return typeof o.coefficient === 'number' && typeof o.area === 'number' &&
+                typeof o.tempSurface === 'number' && typeof o.tempFluid === 'number';
+        case 'radiation':
+            return typeof o.emissivity === 'number' && typeof o.area === 'number' &&
+                typeof o.tempHot === 'number' && typeof o.tempCold === 'number';
+        default: return false;
+    }
+}
+//# sourceMappingURL=guards.js.map

package/dist/chemical/guards.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"guards.js","sourceRoot":"","sources":["../../src/chemical/guards.ts"],"names":[],"mappings":"AAEA;;;GAGG;AACH,MAAM,UAAU,eAAe,CAAC,KAAc;IAC5C,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC;IAC9D,MAAM,CAAC,GAAG,KAAgC,CAAC;IAC3C,QAAQ,CAAC,CAAC,QAAQ,EAAE,CAAC;QACnB,KAAK,IAAI,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,CAAC;QACnG,KAAK,IAAI,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,CAAC;QACnG,KAAK,IAAI,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,CAAC;QACnG,KAAK,IAAI,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,CAAC;QACnG,OAAO,CAAC,CAAC,OAAO,KAAK,CAAC;IACxB,CAAC;AACH,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,cAAc,CAAC,KAAc;IAC3C,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC;IAC9D,MAAM,CAAC,GAAG,KAAgC,CAAC;IAC3C,IAAI,OAAO,CAAC,CAAC,QAAQ,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,SAAS,KAAK,QAAQ;QAAE,OAAO,KAAK,CAAC;IACpF,QAAQ,CAAC,CAAC,KAAK,EAAE,CAAC;QAChB,KAAK,aAAa,CAAC,CAAC,OAAO,OAAO,CAAC,CAAC,MAAM,KAAK,QAAQ,CAAC;QACxD,KAAK,WAAW,CAAC,CAAC,OAAO,IAAI,CAAC;QAC9B,OAAO,CAAC,CAAC,OAAO,KAAK,CAAC;IACxB,CAAC;AACH,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,mBAAmB,CAAC,KAAc;IAChD,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC;IAC9D,MAAM,CAAC,GAAG,KAAgC,CAAC;IAC3C,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC;QACf,KAAK,YAAY;YACf,OAAO,OAAO,CAAC,CAAC,YAAY,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,IAAI,KAAK,QAAQ;gBACrE,OAAO,CAAC,CAAC,SAAS,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,OAAO,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,QAAQ,KAAK,QAAQ,CAAC;QACvG,KAAK,YAAY;YACf,OAAO,OAAO,CAAC,CAAC,WAAW,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,IAAI,KAAK,QAAQ;gBACpE,OAAO,CAAC,CAAC,WAAW,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,SAAS,KAAK,QAAQ,CAAC;QACzE,KAAK,WAAW;YACd,OAAO,OAAO,CAAC,CAAC,UAAU,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,IAAI,KAAK,QAAQ;gBACnE,OAAO,CAAC,CAAC,OAAO,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,QAAQ,KAAK,QAAQ,CAAC;QACpE,OAAO,CAAC,CAAC,OAAO,KAAK,CAAC;IACxB,CAAC;AACH,CAAC"}

package/dist/chemical/heatTransfer.d.ts

@@ -9,6 +9,12 @@ import type { HeatTransferInput, HeatTransferResult } from './types.js';
  *
  * All temperatures in °C internally converted to K for radiation.
  *
+ * @throws {RangeError} Thermal conductivity must be positive
+ * @throws {RangeError} Area must be positive
+ * @throws {RangeError} Thickness must be positive
+ * @throws {RangeError} Convection coefficient must be positive
+ * @throws {RangeError} Emissivity must be between 0 and 1
+ * @throws {RangeError} Absolute temperature must be positive
  * @param input - discriminated union by `mode`
  * @returns heat transfer rate and related quantities
  */

package/dist/chemical/heatTransfer.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"heatTransfer.d.ts","sourceRoot":"","sources":["../../src/chemical/heatTransfer.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,YAAY,CAAC;AAOxE;;;;;;;;;;;;GAYG;AACH,wBAAgB,YAAY,CAAC,KAAK,EAAE,iBAAiB,GAAG,kBAAkB,CA4DzE"}
+{"version":3,"file":"heatTransfer.d.ts","sourceRoot":"","sources":["../../src/chemical/heatTransfer.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,YAAY,CAAC;AAOxE;;;;;;;;;;;;;;;;;;GAkBG;AACH,wBAAgB,YAAY,CAAC,KAAK,EAAE,iBAAiB,GAAG,kBAAkB,CA4DzE"}

package/dist/chemical/heatTransfer.js

@@ -13,6 +13,12 @@ const STEFAN_BOLTZMANN = 5.670374419e-8;
  *
  * All temperatures in °C internally converted to K for radiation.
  *
+ * @throws {RangeError} Thermal conductivity must be positive
+ * @throws {RangeError} Area must be positive
+ * @throws {RangeError} Thickness must be positive
+ * @throws {RangeError} Convection coefficient must be positive
+ * @throws {RangeError} Emissivity must be between 0 and 1
+ * @throws {RangeError} Absolute temperature must be positive
  * @param input - discriminated union by `mode`
  * @returns heat transfer rate and related quantities
  */
@@ -26,11 +32,11 @@ export function heatTransfer(input) {
             const { conductivity, thickness, tempHot, tempCold } = input;
             area = input.area;
             if (conductivity <= 0)
-                throw new Error('Thermal conductivity must be positive');
+                throw new RangeError('Thermal conductivity must be positive');
             if (area <= 0)
-                throw new Error('Area must be positive');
+                throw new RangeError('Area must be positive');
             if (thickness <= 0)
-                throw new Error('Thickness must be positive');
+                throw new RangeError('Thickness must be positive');
             tempDifference = tempHot - tempCold;
             // Fourier's Law: Q = k × A × ΔT / L
             heatRate = (conductivity * area * tempDifference) / thickness;
@@ -42,9 +48,9 @@ export function heatTransfer(input) {
             const { coefficient, tempSurface, tempFluid } = input;
             area = input.area;
             if (coefficient <= 0)
-                throw new Error('Convection coefficient must be positive');
+                throw new RangeError('Convection coefficient must be positive');
             if (area <= 0)
-                throw new Error('Area must be positive');
+                throw new RangeError('Area must be positive');
             tempDifference = tempSurface - tempFluid;
             // Newton's Law of Cooling: Q = h × A × ΔT
             heatRate = coefficient * area * tempDifference;
@@ -56,21 +62,21 @@ export function heatTransfer(input) {
             const { emissivity, tempHot, tempCold } = input;
             area = input.area;
             if (emissivity <= 0 || emissivity > 1)
-                throw new Error('Emissivity must be between 0 and 1');
+                throw new RangeError('Emissivity must be between 0 and 1');
             if (area <= 0)
-                throw new Error('Area must be positive');
+                throw new RangeError('Area must be positive');
             tempDifference = tempHot - tempCold;
             // Convert to Kelvin for Stefan-Boltzmann
             const tHotK = tempHot + 273.15;
             const tColdK = tempCold + 273.15;
             if (tHotK <= 0 || tColdK <= 0)
-                throw new Error('Absolute temperature must be positive');
+                throw new RangeError('Absolute temperature must be positive');
             // Stefan-Boltzmann Law: Q = ε × σ × A × (T_h⁴ - T_c⁴)
             heatRate = emissivity * STEFAN_BOLTZMANN * area * (Math.pow(tHotK, 4) - Math.pow(tColdK, 4));
             // Linearized thermal resistance (approximate): R ≈ 1 / (h_rad × A)
             // where h_rad = ε × σ × (T_h² + T_c²) × (T_h + T_c)
             const hRad = emissivity * STEFAN_BOLTZMANN * (tHotK * tHotK + tColdK * tColdK) * (tHotK + tColdK);
-            thermalResistance = hRad > 0 ? 1 / (hRad * area) : Infinity;
+            thermalResistance = 1 / (hRad * area);
             break;
         }
     }

package/dist/chemical/heatTransfer.js.map

@@ -1 +1 @@
-{"version":3,"file":"heatTransfer.js","sourceRoot":"","sources":["../../src/chemical/heatTransfer.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;GAEG;AACH,MAAM,gBAAgB,GAAG,cAAc,CAAC;AAExC;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,YAAY,CAAC,KAAwB;IACnD,IAAI,QAAgB,CAAC;IACrB,IAAI,IAAY,CAAC;IACjB,IAAI,cAAsB,CAAC;IAC3B,IAAI,iBAAyB,CAAC;IAE9B,QAAQ,KAAK,CAAC,IAAI,EAAE,CAAC;QACnB,KAAK,YAAY,CAAC,CAAC,CAAC;YAClB,MAAM,EAAE,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,KAAK,CAAC;YAC7D,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,YAAY,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;YAChF,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uBAAuB,CAAC,CAAC;YACxD,IAAI,SAAS,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;YAClE,cAAc,GAAG,OAAO,GAAG,QAAQ,CAAC;YACpC,oCAAoC;YACpC,QAAQ,GAAG,CAAC,YAAY,GAAG,IAAI,GAAG,cAAc,CAAC,GAAG,SAAS,CAAC;YAC9D,sCAAsC;YACtC,iBAAiB,GAAG,SAAS,GAAG,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC;YACtD,MAAM;QACR,CAAC;QACD,KAAK,YAAY,CAAC,CAAC,CAAC;YAClB,MAAM,EAAE,WAAW,EAAE,WAAW,EAAE,SAAS,EAAE,GAAG,KAAK,CAAC;YACtD,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,WAAW,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;YACjF,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uBAAuB,CAAC,CAAC;YACxD,cAAc,GAAG,WAAW,GAAG,SAAS,CAAC;YACzC,0CAA0C;YAC1C,QAAQ,GAAG,WAAW,GAAG,IAAI,GAAG,cAAc,CAAC;YAC/C,sCAAsC;YACtC,iBAAiB,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,IAAI,CAAC,CAAC;YAC7C,MAAM;QACR,CAAC;QACD,KAAK,WAAW,CAAC,CAAC,CAAC;YACjB,MAAM,EAAE,UAAU,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,KAAK,CAAC;YAChD,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,UAAU,IAAI,CAAC,IAAI,UAAU,GAAG,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;YAC7F,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uBAAuB,CAAC,CAAC;YACxD,cAAc,GAAG,OAAO,GAAG,QAAQ,CAAC;YACpC,yCAAyC;YACzC,MAAM,KAAK,GAAG,OAAO,GAAG,MAAM,CAAC;YAC/B,MAAM,MAAM,GAAG,QAAQ,GAAG,MAAM,CAAC;YACjC,IAAI,KAAK,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;YACxF,sDAAsD;YACtD,QAAQ,GAAG,UAAU,GAAG,gBAAgB,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7F,mEAAmE;YACnE,oDAAoD;YACpD,MAAM,IAAI,GAAG,UAAU,GAAG,gBAAgB,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;YAClG,iBAAiB,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;YAC5D,MAAM;QACR,CAAC;IACH,CAAC;IAED,MAAM,eAAe,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAEvD,OAAO;QACL,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,cAAc,EAAE,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;KACjD,CAAC;AACJ,CAAC"}
+{"version":3,"file":"heatTransfer.js","sourceRoot":"","sources":["../../src/chemical/heatTransfer.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;GAEG;AACH,MAAM,gBAAgB,GAAG,cAAc,CAAC;AAExC;;;;;;;;;;;;;;;;;;GAkBG;AACH,MAAM,UAAU,YAAY,CAAC,KAAwB;IACnD,IAAI,QAAgB,CAAC;IACrB,IAAI,IAAY,CAAC;IACjB,IAAI,cAAsB,CAAC;IAC3B,IAAI,iBAAyB,CAAC;IAE9B,QAAQ,KAAK,CAAC,IAAI,EAAE,CAAC;QACnB,KAAK,YAAY,CAAC,CAAC,CAAC;YAClB,MAAM,EAAE,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,KAAK,CAAC;YAC7D,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,YAAY,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,uCAAuC,CAAC,CAAC;YACrF,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,uBAAuB,CAAC,CAAC;YAC7D,IAAI,SAAS,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;YACvE,cAAc,GAAG,OAAO,GAAG,QAAQ,CAAC;YACpC,oCAAoC;YACpC,QAAQ,GAAG,CAAC,YAAY,GAAG,IAAI,GAAG,cAAc,CAAC,GAAG,SAAS,CAAC;YAC9D,sCAAsC;YACtC,iBAAiB,GAAG,SAAS,GAAG,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC;YACtD,MAAM;QACR,CAAC;QACD,KAAK,YAAY,CAAC,CAAC,CAAC;YAClB,MAAM,EAAE,WAAW,EAAE,WAAW,EAAE,SAAS,EAAE,GAAG,KAAK,CAAC;YACtD,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,WAAW,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,yCAAyC,CAAC,CAAC;YACtF,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,uBAAuB,CAAC,CAAC;YAC7D,cAAc,GAAG,WAAW,GAAG,SAAS,CAAC;YACzC,0CAA0C;YAC1C,QAAQ,GAAG,WAAW,GAAG,IAAI,GAAG,cAAc,CAAC;YAC/C,sCAAsC;YACtC,iBAAiB,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,IAAI,CAAC,CAAC;YAC7C,MAAM;QACR,CAAC;QACD,KAAK,WAAW,CAAC,CAAC,CAAC;YACjB,MAAM,EAAE,UAAU,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,KAAK,CAAC;YAChD,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,IAAI,UAAU,IAAI,CAAC,IAAI,UAAU,GAAG,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,oCAAoC,CAAC,CAAC;YAClG,IAAI,IAAI,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,uBAAuB,CAAC,CAAC;YAC7D,cAAc,GAAG,OAAO,GAAG,QAAQ,CAAC;YACpC,yCAAyC;YACzC,MAAM,KAAK,GAAG,OAAO,GAAG,MAAM,CAAC;YAC/B,MAAM,MAAM,GAAG,QAAQ,GAAG,MAAM,CAAC;YACjC,IAAI,KAAK,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,uCAAuC,CAAC,CAAC;YAC7F,sDAAsD;YACtD,QAAQ,GAAG,UAAU,GAAG,gBAAgB,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7F,mEAAmE;YACnE,oDAAoD;YACpD,MAAM,IAAI,GAAG,UAAU,GAAG,gBAAgB,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;YAClG,iBAAiB,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;YACtC,MAAM;QACR,CAAC;IACH,CAAC;IAED,MAAM,eAAe,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAEvD,OAAO;QACL,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,eAAe,EAAE,OAAO,CAAC,eAAe,EAAE,CAAC,CAAC;QAC5C,cAAc,EAAE,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;KACjD,CAAC;AACJ,CAAC"}

package/dist/chemical/index.d.ts

@@ -1,4 +1,5 @@
-export type { HeatTransferMode, HeatTransferInput, HeatTransferResult, PipeMaterial, PipeFlowInput, PipeFlowResult, BatchIngredient, BatchInput, ScaledIngredient, BatchResult, ConcentrationUnit, ConcentrationInput, ConcentrationResult, DilutionSolveFor, DilutionInput, DilutionResult, BufferSystem, PhInput, PhResult, ReactorShape, ReactorInput, ReactorResult, ShelfLifeInput, ShelfLifeResult, ResinType, InjectionCycleInput, InjectionCyclePhase, InjectionCycleResult, } from './types.js';
+export type { HeatTransferMode, HeatTransferInput, HeatTransferResult, PipeMaterial, PipeFlowInput, PipeFlowResult, BatchIngredient, BatchInput, ScaledIngredient, BatchResult, ConcentrationUnit, ConcentrationInput, ConcentrationResult, DilutionSolveFor, DilutionInput, DilutionResult, BufferSystem, PhInput, PhResult, ReactorShape, ReactorInput, ReactorResult, ShelfLifeInput, ShelfLifeResult, ResinType, InjectionCycleInput, InjectionCyclePhase, InjectionCycleResult, FlowControlInput, FlowControlResult, ReliefValveInput, ReliefValveResult, PidInput, PidResult, } from './types.js';
+export { isDilutionInput, isReactorInput, isHeatTransferInput } from './guards.js';
 export { heatTransfer } from './heatTransfer.js';
 export { pipeFlow } from './pipeFlow.js';
 export { batch } from './batch.js';
@@ -8,4 +9,7 @@ export { ph } from './ph.js';
 export { reactor } from './reactor.js';
 export { shelfLife } from './shelfLife.js';
 export { injectionCycle } from './injectionCycle.js';
+export { flowControl } from './flowControl.js';
+export { reliefValve } from './reliefValve.js';
+export { pid } from './pid.js';
 //# sourceMappingURL=index.d.ts.map

package/dist/chemical/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/chemical/index.ts"],"names":[],"mappings":"AACA,YAAY,EAEV,gBAAgB,EAChB,iBAAiB,EACjB,kBAAkB,EAElB,YAAY,EACZ,aAAa,EACb,cAAc,EAEd,eAAe,EACf,UAAU,EACV,gBAAgB,EAChB,WAAW,EAEX,iBAAiB,EACjB,kBAAkB,EAClB,mBAAmB,EAEnB,gBAAgB,EAChB,aAAa,EACb,cAAc,EAEd,YAAY,EACZ,OAAO,EACP,QAAQ,EAER,YAAY,EACZ,YAAY,EACZ,aAAa,EAEb,cAAc,EACd,eAAe,EAEf,SAAS,EACT,mBAAmB,EACnB,mBAAmB,EACnB,oBAAoB,GACrB,MAAM,YAAY,CAAC;AAGpB,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,EAAE,EAAE,MAAM,SAAS,CAAC;AAC7B,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/chemical/index.ts"],"names":[],"mappings":"AACA,YAAY,EAEV,gBAAgB,EAChB,iBAAiB,EACjB,kBAAkB,EAElB,YAAY,EACZ,aAAa,EACb,cAAc,EAEd,eAAe,EACf,UAAU,EACV,gBAAgB,EAChB,WAAW,EAEX,iBAAiB,EACjB,kBAAkB,EAClB,mBAAmB,EAEnB,gBAAgB,EAChB,aAAa,EACb,cAAc,EAEd,YAAY,EACZ,OAAO,EACP,QAAQ,EAER,YAAY,EACZ,YAAY,EACZ,aAAa,EAEb,cAAc,EACd,eAAe,EAEf,SAAS,EACT,mBAAmB,EACnB,mBAAmB,EACnB,oBAAoB,EAEpB,gBAAgB,EAChB,iBAAiB,EAEjB,gBAAgB,EAChB,iBAAiB,EAEjB,QAAQ,EACR,SAAS,GACV,MAAM,YAAY,CAAC;AAGpB,OAAO,EAAE,eAAe,EAAE,cAAc,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAGnF,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,EAAE,EAAE,MAAM,SAAS,CAAC;AAC7B,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC"}

package/dist/chemical/index.js

@@ -1,3 +1,5 @@
+// Guards
+export { isDilutionInput, isReactorInput, isHeatTransferInput } from './guards.js';
 // Functions
 export { heatTransfer } from './heatTransfer.js';
 export { pipeFlow } from './pipeFlow.js';
@@ -8,4 +10,7 @@ export { ph } from './ph.js';
 export { reactor } from './reactor.js';
 export { shelfLife } from './shelfLife.js';
 export { injectionCycle } from './injectionCycle.js';
+export { flowControl } from './flowControl.js';
+export { reliefValve } from './reliefValve.js';
+export { pid } from './pid.js';
 //# sourceMappingURL=index.js.map

package/dist/chemical/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/chemical/index.ts"],"names":[],"mappings":"AAyCA,YAAY;AACZ,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,EAAE,EAAE,MAAM,SAAS,CAAC;AAC7B,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/chemical/index.ts"],"names":[],"mappings":"AAkDA,SAAS;AACT,OAAO,EAAE,eAAe,EAAE,cAAc,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAEnF,YAAY;AACZ,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,EAAE,EAAE,MAAM,SAAS,CAAC;AAC7B,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAC;AACvC,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,cAAc,EAAE,MAAM,qBAAqB,CAAC;AACrD,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC"}

package/dist/chemical/pid.d.ts

@@ -0,0 +1,14 @@
+import type { PidInput, PidResult } from './types.js';
+/**
+ * PID Tuning Calculator — Ziegler-Nichols / Cohen-Coon
+ *
+ * @formula
+ *   Z-N Step Response: Kp = T/(K×L), Ti = 2L, Td = 0.5L (for PID)
+ *   Z-N Ultimate:      Kp = 0.6×Ku, Ti = 0.5×Pu, Td = 0.125×Pu (for PID)
+ *   Cohen-Coon:        Kp/K × (T/L) × f(L/T) with dead-time ratio corrections
+ *
+ * @reference Ziegler, J.G. & Nichols, N.B. (1942). Optimum Settings for Automatic Controllers.
+ * @reference Cohen, G.H. & Coon, G.A. (1953). Theoretical Consideration of Retarded Control.
+ */
+export declare function pid(input: PidInput): PidResult;
+//# sourceMappingURL=pid.d.ts.map

package/dist/chemical/pid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"pid.d.ts","sourceRoot":"","sources":["../../src/chemical/pid.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM,YAAY,CAAC;AAEtD;;;;;;;;;;GAUG;AACH,wBAAgB,GAAG,CAAC,KAAK,EAAE,QAAQ,GAAG,SAAS,CA8F9C"}

package/dist/chemical/pid.js

@@ -0,0 +1,103 @@
+import { roundTo } from '../utils.js';
+/**
+ * PID Tuning Calculator — Ziegler-Nichols / Cohen-Coon
+ *
+ * @formula
+ *   Z-N Step Response: Kp = T/(K×L), Ti = 2L, Td = 0.5L (for PID)
+ *   Z-N Ultimate:      Kp = 0.6×Ku, Ti = 0.5×Pu, Td = 0.125×Pu (for PID)
+ *   Cohen-Coon:        Kp/K × (T/L) × f(L/T) with dead-time ratio corrections
+ *
+ * @reference Ziegler, J.G. & Nichols, N.B. (1942). Optimum Settings for Automatic Controllers.
+ * @reference Cohen, G.H. & Coon, G.A. (1953). Theoretical Consideration of Retarded Control.
+ */
+export function pid(input) {
+    const { method, controllerType } = input;
+    let kp = 0;
+    let ti = Infinity;
+    let td = 0;
+    let methodName = method;
+    if (method === 'ziegler-nichols-step') {
+        const K = input.processGain ?? 1;
+        const L = input.deadTime ?? 1;
+        const T = input.timeConstant ?? 1;
+        if (K <= 0 || L <= 0 || T <= 0) {
+            return { kp: 0, ki: 0, kd: 0, ti: 0, td: 0, method: methodName };
+        }
+        if (controllerType === 'P') {
+            kp = T / (K * L);
+            ti = Infinity;
+            td = 0;
+        }
+        else if (controllerType === 'PI') {
+            kp = 0.9 * T / (K * L);
+            ti = L / 0.3;
+            td = 0;
+        }
+        else {
+            // PID
+            kp = 1.2 * T / (K * L);
+            ti = 2 * L;
+            td = 0.5 * L;
+        }
+    }
+    else if (method === 'ziegler-nichols-ultimate') {
+        const Ku = input.ultimateGain ?? 1;
+        const Pu = input.ultimatePeriod ?? 1;
+        if (Ku <= 0 || Pu <= 0) {
+            return { kp: 0, ki: 0, kd: 0, ti: 0, td: 0, method: methodName };
+        }
+        if (controllerType === 'P') {
+            kp = 0.5 * Ku;
+            ti = Infinity;
+            td = 0;
+        }
+        else if (controllerType === 'PI') {
+            kp = 0.45 * Ku;
+            ti = Pu / 1.2;
+            td = 0;
+        }
+        else {
+            // PID
+            kp = 0.6 * Ku;
+            ti = 0.5 * Pu;
+            td = 0.125 * Pu;
+        }
+    }
+    else if (method === 'cohen-coon') {
+        const K = input.processGain ?? 1;
+        const L = input.deadTime ?? 1;
+        const T = input.timeConstant ?? 1;
+        if (K <= 0 || L <= 0 || T <= 0) {
+            return { kp: 0, ki: 0, kd: 0, ti: 0, td: 0, method: methodName };
+        }
+        const r = L / T; // dead time ratio
+        if (controllerType === 'P') {
+            kp = (1 / K) * (T / L) * (1 + r / 3);
+            ti = Infinity;
+            td = 0;
+        }
+        else if (controllerType === 'PI') {
+            kp = (1 / K) * (T / L) * (0.9 + r / 12);
+            ti = L * (30 + 3 * r) / (9 + 20 * r);
+            td = 0;
+        }
+        else {
+            // PID
+            kp = (1 / K) * (T / L) * (4 / 3 + r / 4);
+            ti = L * (32 + 6 * r) / (13 + 8 * r);
+            td = L * 4 / (11 + 2 * r);
+        }
+    }
+    // Ki = Kp / Ti, Kd = Kp × Td
+    const ki = ti > 0 && Number.isFinite(ti) ? kp / ti : 0;
+    const kd = kp * td;
+    return {
+        kp: roundTo(kp, 4),
+        ki: roundTo(ki, 6),
+        kd: roundTo(kd, 4),
+        ti: Number.isFinite(ti) ? roundTo(ti, 4) : 0,
+        td: roundTo(td, 4),
+        method: methodName,
+    };
+}
+//# sourceMappingURL=pid.js.map

package/dist/chemical/pid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pid.js","sourceRoot":"","sources":["../../src/chemical/pid.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,GAAG,CAAC,KAAe;IACjC,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,GAAG,KAAK,CAAC;IAEzC,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,QAAQ,CAAC;IAClB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,UAAU,GAAG,MAAM,CAAC;IAExB,IAAI,MAAM,KAAK,sBAAsB,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,WAAW,IAAI,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,KAAK,CAAC,QAAQ,IAAI,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,IAAI,CAAC,CAAC;QAElC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YAC/B,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC;QACnE,CAAC;QAED,IAAI,cAAc,KAAK,GAAG,EAAE,CAAC;YAC3B,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACjB,EAAE,GAAG,QAAQ,CAAC;YACd,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;YACnC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvB,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;YACb,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,CAAC;YACN,MAAM;YACN,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACX,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC;QACf,CAAC;IAEH,CAAC;SAAM,IAAI,MAAM,KAAK,0BAA0B,EAAE,CAAC;QACjD,MAAM,EAAE,GAAG,KAAK,CAAC,YAAY,IAAI,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,KAAK,CAAC,cAAc,IAAI,CAAC,CAAC;QAErC,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;YACvB,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC;QACnE,CAAC;QAED,IAAI,cAAc,KAAK,GAAG,EAAE,CAAC;YAC3B,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;YACd,EAAE,GAAG,QAAQ,CAAC;YACd,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;YACnC,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;YACf,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC;YACd,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,CAAC;YACN,MAAM;YACN,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;YACd,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;YACd,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;QAClB,CAAC;IAEH,CAAC;SAAM,IAAI,MAAM,KAAK,YAAY,EAAE,CAAC;QACnC,MAAM,CAAC,GAAG,KAAK,CAAC,WAAW,IAAI,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,KAAK,CAAC,QAAQ,IAAI,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,IAAI,CAAC,CAAC;QAElC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YAC/B,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC;QACnE,CAAC;QAED,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,kBAAkB;QAEnC,IAAI,cAAc,KAAK,GAAG,EAAE,CAAC;YAC3B,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACrC,EAAE,GAAG,QAAQ,CAAC;YACd,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;YACnC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;YACxC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;YACrC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC;aAAM,CAAC;YACN,MAAM;YACN,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACzC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACrC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5B,CAAC;IACH,CAAC;IAED,6BAA6B;IAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAEnB,OAAO;QACL,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,EAAE,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5C,EAAE,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAClB,MAAM,EAAE,UAAU;KACnB,CAAC;AACJ,CAAC"}

package/dist/chemical/pipeFlow.d.ts

@@ -10,6 +10,12 @@ import type { PipeFlowInput, PipeFlowResult } from './types.js';
  *
  * Darcy-Weisbach: ΔP = f × (L/D) × (ρv²/2)
  *
+ * @throws {RangeError} Flow rate must be positive
+ * @throws {RangeError} Pipe diameter must be positive
+ * @throws {RangeError} Pipe length must be positive
+ * @throws {RangeError} Fluid density must be positive
+ * @throws {RangeError} Fluid viscosity must be positive
+ * @throws {RangeError} Custom roughness must be provided and non-negative
  * @param input - flow rate, pipe geometry, fluid properties
  * @returns pressure drop, velocity, Reynolds number, flow regime
  */

package/dist/chemical/pipeFlow.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"pipeFlow.d.ts","sourceRoot":"","sources":["../../src/chemical/pipeFlow.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,aAAa,EAAE,cAAc,EAAgB,MAAM,YAAY,CAAC;AAiC9E;;;;;;;;;;;;;GAaG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,aAAa,GAAG,cAAc,CAoE7D"}
+{"version":3,"file":"pipeFlow.d.ts","sourceRoot":"","sources":["../../src/chemical/pipeFlow.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,aAAa,EAAE,cAAc,EAAgB,MAAM,YAAY,CAAC;AAiC9E;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,aAAa,GAAG,cAAc,CAoE7D"}