formulab 0.1.11 → 0.5.0

package/dist/battery/types.d.ts

@@ -0,0 +1,175 @@
+/**
+ * Battery Domain Types
+ */
+export type BatteryChemistry = 'LFP' | 'NMC' | 'NCA' | 'LTO' | 'LCO' | 'LeadAcid' | 'NiMH';
+/**
+ * Energy Density Types
+ */
+export interface EnergyDensityInput {
+    capacityAh: number;
+    nominalVoltage: number;
+    massKg?: number;
+    volumeL?: number;
+}
+export interface EnergyDensityResult {
+    energyWh: number;
+    gravimetricWhPerKg: number | null;
+    volumetricWhPerL: number | null;
+}
+/**
+ * C-Rate Types
+ */
+export type CRateMode = 'currentToRate' | 'rateToCurrent';
+export type CRateInput = {
+    mode: 'currentToRate';
+    capacityAh: number;
+    currentA: number;
+} | {
+    mode: 'rateToCurrent';
+    capacityAh: number;
+    cRate: number;
+};
+export interface CRateResult {
+    cRate: number;
+    currentA: number;
+    theoreticalTimeH: number;
+    theoreticalTimeMin: number;
+}
+/**
+ * State of Health Types
+ */
+export interface StateOfHealthInput {
+    measuredCapacityAh: number;
+    ratedCapacityAh: number;
+}
+export type SohStatus = 'excellent' | 'good' | 'degraded' | 'poor' | 'endOfLife';
+export interface StateOfHealthResult {
+    sohPercent: number;
+    capacityLoss: number;
+    capacityLossPercent: number;
+    status: SohStatus;
+}
+/**
+ * Battery Pack Config Types
+ */
+export interface BatteryPackConfigInput {
+    cellVoltage: number;
+    cellCapacityAh: number;
+    targetVoltage: number;
+    targetCapacityAh: number;
+}
+export interface BatteryPackConfigResult {
+    seriesCells: number;
+    parallelCells: number;
+    totalCells: number;
+    actualVoltage: number;
+    actualCapacityAh: number;
+    totalEnergyWh: number;
+    totalEnergyKWh: number;
+}
+/**
+ * Cycle Life Types
+ */
+export interface CycleLifeInput {
+    chemistry: BatteryChemistry;
+    depthOfDischarge: number;
+    temperatureC: number;
+}
+export interface CycleLifeResult {
+    estimatedCycles: number;
+    baseCycles: number;
+    dodFactor: number;
+    temperatureFactor: number;
+    chemistry: BatteryChemistry;
+}
+/**
+ * Internal Resistance Types
+ */
+export interface InternalResistanceInput {
+    openCircuitVoltage: number;
+    loadVoltage: number;
+    loadCurrentA: number;
+}
+export interface InternalResistanceResult {
+    resistanceOhm: number;
+    resistanceMilliOhm: number;
+    voltageDrop: number;
+    powerLossW: number;
+}
+/**
+ * Self-Discharge Types
+ */
+export interface SelfDischargeInput {
+    initialVoltage: number;
+    finalVoltage: number;
+    days: number;
+    nominalVoltage: number;
+}
+export interface SelfDischargeResult {
+    voltageDropPerDay: number;
+    monthlyRatePercent: number;
+    dailyRatePercent: number;
+}
+/**
+ * Thermal Runaway Types
+ */
+export interface ThermalRunawayInput {
+    ambientTempC: number;
+    currentA: number;
+    internalResistanceOhm: number;
+    heatTransferCoeff: number;
+    surfaceAreaM2: number;
+    runawayTempC: number;
+}
+export interface ThermalRunawayResult {
+    steadyStateTempC: number;
+    temperatureRiseC: number;
+    heatGenerationW: number;
+    heatDissipationW: number;
+    safetyMarginC: number;
+    isSafe: boolean;
+}
+/**
+ * BMS Balancing Types
+ */
+export interface BmsBalancingInput {
+    cellVoltages: number[];
+    balancingCurrentMA: number;
+    cellCapacityAh: number;
+}
+export interface BmsBalancingResult {
+    maxVoltage: number;
+    minVoltage: number;
+    voltageDelta: number;
+    averageVoltage: number;
+    maxBalancingTimeMin: number;
+    maxBalancingTimeH: number;
+    isBalanced: boolean;
+    cellDetails: {
+        cellIndex: number;
+        voltage: number;
+        deltaFromAvg: number;
+        balancingTimeMin: number;
+    }[];
+}
+/**
+ * Charging Profile Types
+ */
+export interface ChargingProfileInput {
+    capacityAh: number;
+    chargingCurrentA: number;
+    cutoffCurrentA: number;
+    ccEndSocPercent?: number;
+}
+export interface ChargingProfileResult {
+    ccPhaseTimeH: number;
+    ccPhaseTimeMin: number;
+    cvPhaseTimeH: number;
+    cvPhaseTimeMin: number;
+    totalTimeH: number;
+    totalTimeMin: number;
+    ccPhaseAh: number;
+    cvPhaseAh: number;
+    averageCRate: number;
+}
+//# sourceMappingURL=types.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"types.d.ts","sourceRoot":"","sources":["../../src/battery/types.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,MAAM,MAAM,gBAAgB,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,UAAU,GAAG,MAAM,CAAC;AAE3F;;GAEG;AACH,MAAM,WAAW,kBAAkB;IACjC,UAAU,EAAE,MAAM,CAAC;IACnB,cAAc,EAAE,MAAM,CAAC;IACvB,MAAM,CAAC,EAAE,MAAM,CAAC;IAChB,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED,MAAM,WAAW,mBAAmB;IAClC,QAAQ,EAAE,MAAM,CAAC;IACjB,kBAAkB,EAAE,MAAM,GAAG,IAAI,CAAC;IAClC,gBAAgB,EAAE,MAAM,GAAG,IAAI,CAAC;CACjC;AAED;;GAEG;AACH,MAAM,MAAM,SAAS,GAAG,eAAe,GAAG,eAAe,CAAC;AAE1D,MAAM,MAAM,UAAU,GAClB;IAAE,IAAI,EAAE,eAAe,CAAC;IAAC,UAAU,EAAE,MAAM,CAAC;IAAC,QAAQ,EAAE,MAAM,CAAA;CAAE,GAC/D;IAAE,IAAI,EAAE,eAAe,CAAC;IAAC,UAAU,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC;AAEjE,MAAM,WAAW,WAAW;IAC1B,KAAK,EAAE,MAAM,CAAC;IACd,QAAQ,EAAE,MAAM,CAAC;IACjB,gBAAgB,EAAE,MAAM,CAAC;IACzB,kBAAkB,EAAE,MAAM,CAAC;CAC5B;AAED;;GAEG;AACH,MAAM,WAAW,kBAAkB;IACjC,kBAAkB,EAAE,MAAM,CAAC;IAC3B,eAAe,EAAE,MAAM,CAAC;CACzB;AAED,MAAM,MAAM,SAAS,GAAG,WAAW,GAAG,MAAM,GAAG,UAAU,GAAG,MAAM,GAAG,WAAW,CAAC;AAEjF,MAAM,WAAW,mBAAmB;IAClC,UAAU,EAAE,MAAM,CAAC;IACnB,YAAY,EAAE,MAAM,CAAC;IACrB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,MAAM,EAAE,SAAS,CAAC;CACnB;AAED;;GAEG;AACH,MAAM,WAAW,sBAAsB;IACrC,WAAW,EAAE,MAAM,CAAC;IACpB,cAAc,EAAE,MAAM,CAAC;IACvB,aAAa,EAAE,MAAM,CAAC;IACtB,gBAAgB,EAAE,MAAM,CAAC;CAC1B;AAED,MAAM,WAAW,uBAAuB;IACtC,WAAW,EAAE,MAAM,CAAC;IACpB,aAAa,EAAE,MAAM,CAAC;IACtB,UAAU,EAAE,MAAM,CAAC;IACnB,aAAa,EAAE,MAAM,CAAC;IACtB,gBAAgB,EAAE,MAAM,CAAC;IACzB,aAAa,EAAE,MAAM,CAAC;IACtB,cAAc,EAAE,MAAM,CAAC;CACxB;AAED;;GAEG;AACH,MAAM,WAAW,cAAc;IAC7B,SAAS,EAAE,gBAAgB,CAAC;IAC5B,gBAAgB,EAAE,MAAM,CAAC;IACzB,YAAY,EAAE,MAAM,CAAC;CACtB;AAED,MAAM,WAAW,eAAe;IAC9B,eAAe,EAAE,MAAM,CAAC;IACxB,UAAU,EAAE,MAAM,CAAC;IACnB,SAAS,EAAE,MAAM,CAAC;IAClB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,SAAS,EAAE,gBAAgB,CAAC;CAC7B;AAED;;GAEG;AACH,MAAM,WAAW,uBAAuB;IACtC,kBAAkB,EAAE,MAAM,CAAC;IAC3B,WAAW,EAAE,MAAM,CAAC;IACpB,YAAY,EAAE,MAAM,CAAC;CACtB;AAED,MAAM,WAAW,wBAAwB;IACvC,aAAa,EAAE,MAAM,CAAC;IACtB,kBAAkB,EAAE,MAAM,CAAC;IAC3B,WAAW,EAAE,MAAM,CAAC;IACpB,UAAU,EAAE,MAAM,CAAC;CACpB;AAED;;GAEG;AACH,MAAM,WAAW,kBAAkB;IACjC,cAAc,EAAE,MAAM,CAAC;IACvB,YAAY,EAAE,MAAM,CAAC;IACrB,IAAI,EAAE,MAAM,CAAC;IACb,cAAc,EAAE,MAAM,CAAC;CACxB;AAED,MAAM,WAAW,mBAAmB;IAClC,iBAAiB,EAAE,MAAM,CAAC;IAC1B,kBAAkB,EAAE,MAAM,CAAC;IAC3B,gBAAgB,EAAE,MAAM,CAAC;CAC1B;AAED;;GAEG;AACH,MAAM,WAAW,mBAAmB;IAClC,YAAY,EAAE,MAAM,CAAC;IACrB,QAAQ,EAAE,MAAM,CAAC;IACjB,qBAAqB,EAAE,MAAM,CAAC;IAC9B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,aAAa,EAAE,MAAM,CAAC;IACtB,YAAY,EAAE,MAAM,CAAC;CACtB;AAED,MAAM,WAAW,oBAAoB;IACnC,gBAAgB,EAAE,MAAM,CAAC;IACzB,gBAAgB,EAAE,MAAM,CAAC;IACzB,eAAe,EAAE,MAAM,CAAC;IACxB,gBAAgB,EAAE,MAAM,CAAC;IACzB,aAAa,EAAE,MAAM,CAAC;IACtB,MAAM,EAAE,OAAO,CAAC;CACjB;AAED;;GAEG;AACH,MAAM,WAAW,iBAAiB;IAChC,YAAY,EAAE,MAAM,EAAE,CAAC;IACvB,kBAAkB,EAAE,MAAM,CAAC;IAC3B,cAAc,EAAE,MAAM,CAAC;CACxB;AAED,MAAM,WAAW,kBAAkB;IACjC,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,YAAY,EAAE,MAAM,CAAC;IACrB,cAAc,EAAE,MAAM,CAAC;IACvB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,UAAU,EAAE,OAAO,CAAC;IACpB,WAAW,EAAE;QACX,SAAS,EAAE,MAAM,CAAC;QAClB,OAAO,EAAE,MAAM,CAAC;QAChB,YAAY,EAAE,MAAM,CAAC;QACrB,gBAAgB,EAAE,MAAM,CAAC;KAC1B,EAAE,CAAC;CACL;AAED;;GAEG;AACH,MAAM,WAAW,oBAAoB;IACnC,UAAU,EAAE,MAAM,CAAC;IACnB,gBAAgB,EAAE,MAAM,CAAC;IACzB,cAAc,EAAE,MAAM,CAAC;IACvB,eAAe,CAAC,EAAE,MAAM,CAAC;CAC1B;AAED,MAAM,WAAW,qBAAqB;IACpC,YAAY,EAAE,MAAM,CAAC;IACrB,cAAc,EAAE,MAAM,CAAC;IACvB,YAAY,EAAE,MAAM,CAAC;IACrB,cAAc,EAAE,MAAM,CAAC;IACvB,UAAU,EAAE,MAAM,CAAC;IACnB,YAAY,EAAE,MAAM,CAAC;IACrB,SAAS,EAAE,MAAM,CAAC;IAClB,SAAS,EAAE,MAAM,CAAC;IAClB,YAAY,EAAE,MAAM,CAAC;CACtB"}

package/dist/battery/types.js

@@ -0,0 +1,5 @@
+/**
+ * Battery Domain Types
+ */
+export {};
+//# sourceMappingURL=types.js.map

package/dist/battery/types.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"types.js","sourceRoot":"","sources":["../../src/battery/types.ts"],"names":[],"mappings":"AAAA;;GAEG"}

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

@@ -1 +1 @@
-{"version":3,"file":"dilution.d.ts","sourceRoot":"","sources":["../../src/chemical/dilution.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,aAAa,EAAE,cAAc,EAAE,MAAM,YAAY,CAAC;AAEhE;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,aAAa,GAAG,cAAc,CAgE7D"}
+{"version":3,"file":"dilution.d.ts","sourceRoot":"","sources":["../../src/chemical/dilution.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,aAAa,EAAE,cAAc,EAAE,MAAM,YAAY,CAAC;AAEhE;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,aAAa,GAAG,cAAc,CA2D7D"}

package/dist/chemical/dilution.js

@@ -8,54 +8,51 @@
  * @throws Error if division by zero would occur
  */
 export function dilution(input) {
-    const { solveFor, c1, v1, c2, v2 } = input;
     let resultC1;
     let resultV1;
     let resultC2;
     let resultV2;
-    switch (solveFor) {
+    switch (input.solveFor) {
         case 'c2':
             // C2 = C1 * V1 / V2
-            if (v2 === undefined || v2 === 0) {
+            if (input.v2 === 0) {
                 throw new Error('Cannot solve for c2: v2 must be non-zero');
             }
-            resultC1 = c1;
-            resultV1 = v1;
-            resultV2 = v2;
-            resultC2 = (c1 * v1) / v2;
+            resultC1 = input.c1;
+            resultV1 = input.v1;
+            resultV2 = input.v2;
+            resultC2 = (input.c1 * input.v1) / input.v2;
             break;
         case 'v2':
             // V2 = C1 * V1 / C2
-            if (c2 === undefined || c2 === 0) {
+            if (input.c2 === 0) {
                 throw new Error('Cannot solve for v2: c2 must be non-zero');
             }
-            resultC1 = c1;
-            resultV1 = v1;
-            resultC2 = c2;
-            resultV2 = (c1 * v1) / c2;
+            resultC1 = input.c1;
+            resultV1 = input.v1;
+            resultC2 = input.c2;
+            resultV2 = (input.c1 * input.v1) / input.c2;
             break;
         case 'c1':
             // C1 = C2 * V2 / V1
-            if (v1 === undefined || v1 === 0) {
+            if (input.v1 === 0) {
                 throw new Error('Cannot solve for c1: v1 must be non-zero');
             }
-            resultV1 = v1;
-            resultC2 = c2;
-            resultV2 = v2;
-            resultC1 = (c2 * v2) / v1;
+            resultV1 = input.v1;
+            resultC2 = input.c2;
+            resultV2 = input.v2;
+            resultC1 = (input.c2 * input.v2) / input.v1;
             break;
         case 'v1':
             // V1 = C2 * V2 / C1
-            if (c1 === undefined || c1 === 0) {
+            if (input.c1 === 0) {
                 throw new Error('Cannot solve for v1: c1 must be non-zero');
             }
-            resultC1 = c1;
-            resultC2 = c2;
-            resultV2 = v2;
-            resultV1 = (c2 * v2) / c1;
+            resultC1 = input.c1;
+            resultC2 = input.c2;
+            resultV2 = input.v2;
+            resultV1 = (input.c2 * input.v2) / input.c1;
             break;
-        default:
-            throw new Error(`Unknown solveFor value: ${solveFor}`);
     }
     return {
         c1: resultC1,

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,MAAM,EAAE,QAAQ,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC;IAE3C,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IACrB,IAAI,QAAgB,CAAC;IAErB,QAAQ,QAAQ,EAAE,CAAC;QACjB,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACjC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAE,CAAC;YACd,QAAQ,GAAG,CAAC,EAAG,GAAG,EAAG,CAAC,GAAG,EAAE,CAAC;YAC5B,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACjC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAE,CAAC;YACd,QAAQ,GAAG,CAAC,EAAG,GAAG,EAAG,CAAC,GAAG,EAAE,CAAC;YAC5B,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACjC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,EAAE,CAAC;YACd,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,CAAC,EAAG,GAAG,EAAG,CAAC,GAAG,EAAE,CAAC;YAC5B,MAAM;QAER,KAAK,IAAI;YACP,oBAAoB;YACpB,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACjC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,CAAC;YACD,QAAQ,GAAG,EAAE,CAAC;YACd,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,EAAG,CAAC;YACf,QAAQ,GAAG,CAAC,EAAG,GAAG,EAAG,CAAC,GAAG,EAAE,CAAC;YAC5B,MAAM;QAER;YACE,MAAM,IAAI,KAAK,CAAC,2BAA2B,QAAQ,EAAE,CAAC,CAAC;IAC3D,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,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"}

package/dist/chemical/heatTransfer.d.ts

@@ -0,0 +1,16 @@
+import type { HeatTransferInput, HeatTransferResult } from './types.js';
+/**
+ * Heat Transfer Calculator
+ *
+ * Supports three modes of heat transfer:
+ * - **Conduction** (Fourier's Law): Q = k × A × ΔT / L
+ * - **Convection** (Newton's Law of Cooling): Q = h × A × ΔT
+ * - **Radiation** (Stefan-Boltzmann Law): Q = ε × σ × A × (T_hot⁴ - T_cold⁴)
+ *
+ * All temperatures in °C internally converted to K for radiation.
+ *
+ * @param input - discriminated union by `mode`
+ * @returns heat transfer rate and related quantities
+ */
+export declare function heatTransfer(input: HeatTransferInput): HeatTransferResult;
+//# sourceMappingURL=heatTransfer.d.ts.map

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

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

package/dist/chemical/heatTransfer.js

@@ -0,0 +1,85 @@
+import { roundTo } from '../utils.js';
+/**
+ * Stefan-Boltzmann constant (W/(m²·K⁴))
+ */
+const STEFAN_BOLTZMANN = 5.670374419e-8;
+/**
+ * Heat Transfer Calculator
+ *
+ * Supports three modes of heat transfer:
+ * - **Conduction** (Fourier's Law): Q = k × A × ΔT / L
+ * - **Convection** (Newton's Law of Cooling): Q = h × A × ΔT
+ * - **Radiation** (Stefan-Boltzmann Law): Q = ε × σ × A × (T_hot⁴ - T_cold⁴)
+ *
+ * All temperatures in °C internally converted to K for radiation.
+ *
+ * @param input - discriminated union by `mode`
+ * @returns heat transfer rate and related quantities
+ */
+export function heatTransfer(input) {
+    let heatRate;
+    let area;
+    let tempDifference;
+    let thermalResistance;
+    switch (input.mode) {
+        case 'conduction': {
+            const { conductivity, thickness, tempHot, tempCold } = input;
+            area = input.area;
+            if (conductivity <= 0)
+                throw new Error('Thermal conductivity must be positive');
+            if (area <= 0)
+                throw new Error('Area must be positive');
+            if (thickness <= 0)
+                throw new Error('Thickness must be positive');
+            tempDifference = tempHot - tempCold;
+            // Fourier's Law: Q = k × A × ΔT / L
+            heatRate = (conductivity * area * tempDifference) / thickness;
+            // Thermal resistance: R = L / (k × A)
+            thermalResistance = thickness / (conductivity * area);
+            break;
+        }
+        case 'convection': {
+            const { coefficient, tempSurface, tempFluid } = input;
+            area = input.area;
+            if (coefficient <= 0)
+                throw new Error('Convection coefficient must be positive');
+            if (area <= 0)
+                throw new Error('Area must be positive');
+            tempDifference = tempSurface - tempFluid;
+            // Newton's Law of Cooling: Q = h × A × ΔT
+            heatRate = coefficient * area * tempDifference;
+            // Thermal resistance: R = 1 / (h × A)
+            thermalResistance = 1 / (coefficient * area);
+            break;
+        }
+        case 'radiation': {
+            const { emissivity, tempHot, tempCold } = input;
+            area = input.area;
+            if (emissivity <= 0 || emissivity > 1)
+                throw new Error('Emissivity must be between 0 and 1');
+            if (area <= 0)
+                throw new Error('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');
+            // 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;
+            break;
+        }
+    }
+    const heatFluxDensity = area > 0 ? heatRate / area : 0;
+    return {
+        heatRate: roundTo(heatRate, 4),
+        heatFluxDensity: roundTo(heatFluxDensity, 4),
+        tempDifference: roundTo(tempDifference, 2),
+        thermalResistance: roundTo(thermalResistance, 6),
+    };
+}
+//# sourceMappingURL=heatTransfer.js.map

package/dist/chemical/heatTransfer.js.map

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

package/dist/chemical/index.d.ts

@@ -1,4 +1,6 @@
-export type { 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, } from './types.js';
+export { heatTransfer } from './heatTransfer.js';
+export { pipeFlow } from './pipeFlow.js';
 export { batch } from './batch.js';
 export { concentration } from './concentration.js';
 export { dilution } from './dilution.js';

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,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,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,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"}

package/dist/chemical/index.js

@@ -1,4 +1,6 @@
 // Functions
+export { heatTransfer } from './heatTransfer.js';
+export { pipeFlow } from './pipeFlow.js';
 export { batch } from './batch.js';
 export { concentration } from './concentration.js';
 export { dilution } from './dilution.js';

package/dist/chemical/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/chemical/index.ts"],"names":[],"mappings":"AAiCA,YAAY;AACZ,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":"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"}

package/dist/chemical/pipeFlow.d.ts

@@ -0,0 +1,17 @@
+import type { PipeFlowInput, PipeFlowResult } from './types.js';
+/**
+ * Pipe Flow Calculator (Darcy-Weisbach + Swamee-Jain)
+ *
+ * Calculates pressure drop and head loss for fluid flow through pipes.
+ * Uses:
+ * - Laminar: f = 64/Re (Re < 2300)
+ * - Turbulent: Swamee-Jain explicit formula (Re > 4000)
+ * - Transitional: Linear interpolation (2300 ≤ Re ≤ 4000)
+ *
+ * Darcy-Weisbach: ΔP = f × (L/D) × (ρv²/2)
+ *
+ * @param input - flow rate, pipe geometry, fluid properties
+ * @returns pressure drop, velocity, Reynolds number, flow regime
+ */
+export declare function pipeFlow(input: PipeFlowInput): PipeFlowResult;
+//# sourceMappingURL=pipeFlow.d.ts.map

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

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

package/dist/chemical/pipeFlow.js

@@ -0,0 +1,112 @@
+import { roundTo } from '../utils.js';
+/**
+ * Pipe roughness values in mm
+ * Source: Engineering Toolbox, Pipe Flow Software, ASTM
+ */
+const PIPE_ROUGHNESS = {
+    commercialSteel: 0.045,
+    stainlessSteel: 0.015,
+    castIron: 0.26,
+    copper: 0.0015,
+    pvc: 0.0015,
+    concrete: 1.0,
+    galvanizedSteel: 0.15,
+};
+/**
+ * Swamee-Jain explicit friction factor for turbulent flow
+ *
+ * f = 0.25 / [log10(ε/(3.7D) + 5.74/Re^0.9)]²
+ *
+ * Valid for: 5000 ≤ Re ≤ 10⁸, 10⁻⁶ ≤ ε/D ≤ 10⁻²
+ * Accuracy: ±1% vs Colebrook-White
+ *
+ * Reference: Swamee & Jain (1976), J. Hydraulics Division, ASCE
+ */
+function swameeJainFriction(Re, epsilonM, diameterM) {
+    const relativeRoughness = epsilonM / diameterM;
+    const term = relativeRoughness / 3.7 + 5.74 / Math.pow(Re, 0.9);
+    const logTerm = Math.log10(term);
+    return 0.25 / (logTerm * logTerm);
+}
+/**
+ * Pipe Flow Calculator (Darcy-Weisbach + Swamee-Jain)
+ *
+ * Calculates pressure drop and head loss for fluid flow through pipes.
+ * Uses:
+ * - Laminar: f = 64/Re (Re < 2300)
+ * - Turbulent: Swamee-Jain explicit formula (Re > 4000)
+ * - Transitional: Linear interpolation (2300 ≤ Re ≤ 4000)
+ *
+ * Darcy-Weisbach: ΔP = f × (L/D) × (ρv²/2)
+ *
+ * @param input - flow rate, pipe geometry, fluid properties
+ * @returns pressure drop, velocity, Reynolds number, flow regime
+ */
+export function pipeFlow(input) {
+    const { flowRate, pipeDiameter, pipeLength, pipeMaterial, fluidDensity, fluidViscosity, customRoughness } = input;
+    if (flowRate <= 0)
+        throw new Error('Flow rate must be positive');
+    if (pipeDiameter <= 0)
+        throw new Error('Pipe diameter must be positive');
+    if (pipeLength <= 0)
+        throw new Error('Pipe length must be positive');
+    if (fluidDensity <= 0)
+        throw new Error('Fluid density must be positive');
+    if (fluidViscosity <= 0)
+        throw new Error('Fluid viscosity must be positive');
+    // Get pipe roughness
+    let roughnessMm;
+    if (pipeMaterial === 'custom') {
+        if (customRoughness == null || customRoughness < 0) {
+            throw new Error('Custom roughness must be provided and non-negative');
+        }
+        roughnessMm = customRoughness;
+    }
+    else {
+        roughnessMm = PIPE_ROUGHNESS[pipeMaterial];
+    }
+    // Unit conversions
+    const diameterM = pipeDiameter / 1000; // mm → m
+    const epsilonM = roughnessMm / 1000; // mm → m
+    const flowRateM3s = flowRate / 60000; // L/min → m³/s
+    // Velocity: v = Q / A = 4Q / (πD²)
+    const area = (Math.PI * diameterM * diameterM) / 4;
+    const velocity = flowRateM3s / area;
+    // Reynolds number: Re = ρvD / μ
+    const Re = (fluidDensity * velocity * diameterM) / fluidViscosity;
+    // Determine flow regime and friction factor
+    let frictionFactor;
+    let flowRegime;
+    if (Re < 2300) {
+        flowRegime = 'laminar';
+        frictionFactor = 64 / Re;
+    }
+    else if (Re > 4000) {
+        flowRegime = 'turbulent';
+        frictionFactor = swameeJainFriction(Re, epsilonM, diameterM);
+    }
+    else {
+        // Transitional: linear interpolation
+        flowRegime = 'transitional';
+        const fLam = 64 / 2300;
+        const fTurb = swameeJainFriction(4000, epsilonM, diameterM);
+        const t = (Re - 2300) / (4000 - 2300);
+        frictionFactor = fLam + t * (fTurb - fLam);
+    }
+    // Darcy-Weisbach: ΔP = f × (L/D) × (ρv²/2)
+    const pressureDrop = frictionFactor * (pipeLength / diameterM) * (fluidDensity * velocity * velocity / 2);
+    // Head loss: h_f = ΔP / (ρg)
+    const g = 9.81;
+    const headLoss = pressureDrop / (fluidDensity * g);
+    return {
+        velocity: roundTo(velocity, 4),
+        reynoldsNumber: roundTo(Re, 0),
+        flowRegime,
+        frictionFactor: roundTo(frictionFactor, 6),
+        pressureDrop: roundTo(pressureDrop, 2),
+        pressureDropKpa: roundTo(pressureDrop / 1000, 4),
+        pressureDropBar: roundTo(pressureDrop / 100000, 6),
+        headLoss: roundTo(headLoss, 4),
+    };
+}
+//# sourceMappingURL=pipeFlow.js.map

package/dist/chemical/pipeFlow.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pipeFlow.js","sourceRoot":"","sources":["../../src/chemical/pipeFlow.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAGtC;;;GAGG;AACH,MAAM,cAAc,GAAoD;IACtE,eAAe,EAAE,KAAK;IACtB,cAAc,EAAE,KAAK;IACrB,QAAQ,EAAE,IAAI;IACd,MAAM,EAAE,MAAM;IACd,GAAG,EAAE,MAAM;IACX,QAAQ,EAAE,GAAG;IACb,eAAe,EAAE,IAAI;CACtB,CAAC;AAEF;;;;;;;;;GASG;AACH,SAAS,kBAAkB,CAAC,EAAU,EAAE,QAAgB,EAAE,SAAiB;IACzE,MAAM,iBAAiB,GAAG,QAAQ,GAAG,SAAS,CAAC;IAC/C,MAAM,IAAI,GAAG,iBAAiB,GAAG,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAChE,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACjC,OAAO,IAAI,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC;AACpC,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,QAAQ,CAAC,KAAoB;IAC3C,MAAM,EAAE,QAAQ,EAAE,YAAY,EAAE,UAAU,EAAE,YAAY,EAAE,YAAY,EAAE,cAAc,EAAE,eAAe,EAAE,GAAG,KAAK,CAAC;IAElH,IAAI,QAAQ,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;IACjE,IAAI,YAAY,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;IACzE,IAAI,UAAU,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IACrE,IAAI,YAAY,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;IACzE,IAAI,cAAc,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;IAE7E,qBAAqB;IACrB,IAAI,WAAmB,CAAC;IACxB,IAAI,YAAY,KAAK,QAAQ,EAAE,CAAC;QAC9B,IAAI,eAAe,IAAI,IAAI,IAAI,eAAe,GAAG,CAAC,EAAE,CAAC;YACnD,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;QACxE,CAAC;QACD,WAAW,GAAG,eAAe,CAAC;IAChC,CAAC;SAAM,CAAC;QACN,WAAW,GAAG,cAAc,CAAC,YAAY,CAAC,CAAC;IAC7C,CAAC;IAED,mBAAmB;IACnB,MAAM,SAAS,GAAG,YAAY,GAAG,IAAI,CAAC,CAAQ,SAAS;IACvD,MAAM,QAAQ,GAAG,WAAW,GAAG,IAAI,CAAC,CAAW,SAAS;IACxD,MAAM,WAAW,GAAG,QAAQ,GAAG,KAAK,CAAC,CAAW,eAAe;IAE/D,mCAAmC;IACnC,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,WAAW,GAAG,IAAI,CAAC;IAEpC,gCAAgC;IAChC,MAAM,EAAE,GAAG,CAAC,YAAY,GAAG,QAAQ,GAAG,SAAS,CAAC,GAAG,cAAc,CAAC;IAElE,4CAA4C;IAC5C,IAAI,cAAsB,CAAC;IAC3B,IAAI,UAAoD,CAAC;IAEzD,IAAI,EAAE,GAAG,IAAI,EAAE,CAAC;QACd,UAAU,GAAG,SAAS,CAAC;QACvB,cAAc,GAAG,EAAE,GAAG,EAAE,CAAC;IAC3B,CAAC;SAAM,IAAI,EAAE,GAAG,IAAI,EAAE,CAAC;QACrB,UAAU,GAAG,WAAW,CAAC;QACzB,cAAc,GAAG,kBAAkB,CAAC,EAAE,EAAE,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC/D,CAAC;SAAM,CAAC;QACN,qCAAqC;QACrC,UAAU,GAAG,cAAc,CAAC;QAC5B,MAAM,IAAI,GAAG,EAAE,GAAG,IAAI,CAAC;QACvB,MAAM,KAAK,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,EAAE,SAAS,CAAC,CAAC;QAC5D,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;QACtC,cAAc,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,CAAC;IAC7C,CAAC;IAED,2CAA2C;IAC3C,MAAM,YAAY,GAAG,cAAc,GAAG,CAAC,UAAU,GAAG,SAAS,CAAC,GAAG,CAAC,YAAY,GAAG,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC;IAE1G,6BAA6B;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC;IACf,MAAM,QAAQ,GAAG,YAAY,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;IAEnD,OAAO;QACL,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;QAC9B,cAAc,EAAE,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC;QAC9B,UAAU;QACV,cAAc,EAAE,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,YAAY,EAAE,OAAO,CAAC,YAAY,EAAE,CAAC,CAAC;QACtC,eAAe,EAAE,OAAO,CAAC,YAAY,GAAG,IAAI,EAAE,CAAC,CAAC;QAChD,eAAe,EAAE,OAAO,CAAC,YAAY,GAAG,MAAM,EAAE,CAAC,CAAC;QAClD,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC;KAC/B,CAAC;AACJ,CAAC"}

package/dist/chemical/reactor.js

@@ -16,13 +16,13 @@ import { roundTo } from '../utils.js';
  * @returns Reactor volumes and surface area
  */
 export function reactor(input) {
-    const { shape, diameter, height, fillRatio } = input;
+    const { shape, diameter, fillRatio } = input;
     const radius = diameter / 2;
     let totalVolume;
     let surfaceArea;
     if (shape === 'cylindrical') {
         // Cylindrical reactor
-        const h = height;
+        const h = input.height;
         // V = PI * r^2 * h
         totalVolume = Math.PI * radius * radius * h;
         // Surface Area = 2 * PI * r^2 (ends) + 2 * PI * r * h (lateral)

package/dist/chemical/reactor.js.map

@@ -1 +1 @@
-{"version":3,"file":"reactor.js","sourceRoot":"","sources":["../../src/chemical/reactor.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAEtC;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,OAAO,CAAC,KAAmB;IACzC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,GAAG,KAAK,CAAC;IAErD,MAAM,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IAC5B,IAAI,WAAmB,CAAC;IACxB,IAAI,WAAmB,CAAC;IAExB,IAAI,KAAK,KAAK,aAAa,EAAE,CAAC;QAC5B,sBAAsB;QACtB,MAAM,CAAC,GAAG,MAAO,CAAC;QAClB,mBAAmB;QACnB,WAAW,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,CAAC,CAAC;QAC5C,gEAAgE;QAChE,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,CAAC,CAAC;IACzE,CAAC;SAAM,CAAC;QACN,oBAAoB;QACpB,uBAAuB;QACvB,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3D,8BAA8B;QAC9B,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;IAC9C,CAAC;IAED,MAAM,aAAa,GAAG,WAAW,GAAG,SAAS,CAAC;IAC9C,MAAM,iBAAiB,GAAG,WAAW,GAAG,IAAI,CAAC,CAAC,gBAAgB;IAC9D,MAAM,mBAAmB,GAAG,aAAa,GAAG,IAAI,CAAC;IACjD,MAAM,oBAAoB,GAAG,WAAW,GAAG,WAAW,CAAC;IAEvD,OAAO;QACL,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC;QACpC,aAAa,EAAE,OAAO,CAAC,aAAa,EAAE,CAAC,CAAC;QACxC,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,mBAAmB,EAAE,OAAO,CAAC,mBAAmB,EAAE,CAAC,CAAC;QACpD,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC;QACpC,oBAAoB,EAAE,OAAO,CAAC,oBAAoB,EAAE,CAAC,CAAC;KACvD,CAAC;AACJ,CAAC;AAED;;GAEG"}
+{"version":3,"file":"reactor.js","sourceRoot":"","sources":["../../src/chemical/reactor.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAEtC;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,OAAO,CAAC,KAAmB;IACzC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,GAAG,KAAK,CAAC;IAE7C,MAAM,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IAC5B,IAAI,WAAmB,CAAC;IACxB,IAAI,WAAmB,CAAC;IAExB,IAAI,KAAK,KAAK,aAAa,EAAE,CAAC;QAC5B,sBAAsB;QACtB,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QACvB,mBAAmB;QACnB,WAAW,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,CAAC,CAAC;QAC5C,gEAAgE;QAChE,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,CAAC,CAAC;IACzE,CAAC;SAAM,CAAC;QACN,oBAAoB;QACpB,uBAAuB;QACvB,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3D,8BAA8B;QAC9B,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;IAC9C,CAAC;IAED,MAAM,aAAa,GAAG,WAAW,GAAG,SAAS,CAAC;IAC9C,MAAM,iBAAiB,GAAG,WAAW,GAAG,IAAI,CAAC,CAAC,gBAAgB;IAC9D,MAAM,mBAAmB,GAAG,aAAa,GAAG,IAAI,CAAC;IACjD,MAAM,oBAAoB,GAAG,WAAW,GAAG,WAAW,CAAC;IAEvD,OAAO;QACL,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC;QACpC,aAAa,EAAE,OAAO,CAAC,aAAa,EAAE,CAAC,CAAC;QACxC,iBAAiB,EAAE,OAAO,CAAC,iBAAiB,EAAE,CAAC,CAAC;QAChD,mBAAmB,EAAE,OAAO,CAAC,mBAAmB,EAAE,CAAC,CAAC;QACpD,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC;QACpC,oBAAoB,EAAE,OAAO,CAAC,oBAAoB,EAAE,CAAC,CAAC;KACvD,CAAC;AACJ,CAAC;AAED;;GAEG"}