Compare commits
1 Commits
1da55fc3f5
...
7c8722b324
| Author | SHA1 | Date | |
|---|---|---|---|
|
|
7c8722b324 |
@@ -3,12 +3,13 @@ module.exports = function(RED) {
|
|||||||
RED.nodes.createNode(this, config);
|
RED.nodes.createNode(this, config);
|
||||||
var node = this;
|
var node = this;
|
||||||
|
|
||||||
|
let name = config.name;
|
||||||
let F2 = parseFloat(config.F2);
|
let F2 = parseFloat(config.F2);
|
||||||
const inlet_F2 = parseInt(config.inlet);
|
const inlet_F2 = parseInt(config.inlet);
|
||||||
|
|
||||||
node.on('input', function(msg, send, done) {
|
node.on('input', function(msg, send, done) {
|
||||||
switch (msg.topic) {
|
switch (msg.topic) {
|
||||||
case "Fluent": {
|
case "Fluent":
|
||||||
// conserve volume flow debit
|
// conserve volume flow debit
|
||||||
let F_in = msg.payload.F;
|
let F_in = msg.payload.F;
|
||||||
let F1 = Math.max(F_in - F2, 0);
|
let F1 = Math.max(F_in - F2, 0);
|
||||||
@@ -23,7 +24,6 @@ module.exports = function(RED) {
|
|||||||
|
|
||||||
send([msg_F1, msg_F2]);
|
send([msg_F1, msg_F2]);
|
||||||
break;
|
break;
|
||||||
}
|
|
||||||
case "clock":
|
case "clock":
|
||||||
break;
|
break;
|
||||||
default:
|
default:
|
||||||
|
|||||||
@@ -3,12 +3,13 @@ module.exports = function(RED) {
|
|||||||
RED.nodes.createNode(this, config);
|
RED.nodes.createNode(this, config);
|
||||||
var node = this;
|
var node = this;
|
||||||
|
|
||||||
|
let name = config.name;
|
||||||
let TS_set = parseFloat(config.TS_set);
|
let TS_set = parseFloat(config.TS_set);
|
||||||
const inlet_sludge = parseInt(config.inlet);
|
const inlet_sludge = parseInt(config.inlet);
|
||||||
|
|
||||||
node.on('input', function(msg, send, done) {
|
node.on('input', function(msg, send, done) {
|
||||||
switch (msg.topic) {
|
switch (msg.topic) {
|
||||||
case "Fluent": {
|
case "Fluent":
|
||||||
// conserve volume flow debit
|
// conserve volume flow debit
|
||||||
let F_in = msg.payload.F;
|
let F_in = msg.payload.F;
|
||||||
let C_in = msg.payload.C;
|
let C_in = msg.payload.C;
|
||||||
@@ -40,7 +41,6 @@ module.exports = function(RED) {
|
|||||||
|
|
||||||
send([msg_F1, msg_F2]);
|
send([msg_F1, msg_F2]);
|
||||||
break;
|
break;
|
||||||
}
|
|
||||||
case "clock":
|
case "clock":
|
||||||
break;
|
break;
|
||||||
default:
|
default:
|
||||||
|
|||||||
35
reactor.html
35
reactor.html
@@ -1,9 +1,19 @@
|
|||||||
|
<!--
|
||||||
|
| S88-niveau | Primair (blokkleur) | Tekstkleur |
|
||||||
|
| ---------------------- | ------------------- | ---------- |
|
||||||
|
| **Area** | `#0f52a5` | wit |
|
||||||
|
| **Process Cell** | `#0c99d9` | wit |
|
||||||
|
| **Unit** | `#50a8d9` | zwart |
|
||||||
|
| **Equipment (Module)** | `#86bbdd` | zwart |
|
||||||
|
| **Control Module** | `#a9daee` | zwart |
|
||||||
|
|
||||||
|
-->
|
||||||
<script src="/reactor/menu.js"></script>
|
<script src="/reactor/menu.js"></script>
|
||||||
|
|
||||||
<script type="text/javascript">
|
<script type="text/javascript">
|
||||||
RED.nodes.registerType("reactor", {
|
RED.nodes.registerType("reactor", {
|
||||||
category: "WWTP",
|
category: "EVOLV",
|
||||||
color: "#c4cce0",
|
color: "#50a8d9",
|
||||||
defaults: {
|
defaults: {
|
||||||
name: { value: "" },
|
name: { value: "" },
|
||||||
reactor_type: { value: "CSTR", required: true },
|
reactor_type: { value: "CSTR", required: true },
|
||||||
@@ -29,8 +39,6 @@
|
|||||||
X_TS_init: { value: 125.0009, required: true },
|
X_TS_init: { value: 125.0009, required: true },
|
||||||
|
|
||||||
timeStep: { value: 1, required: true },
|
timeStep: { value: 1, required: true },
|
||||||
processOutputFormat: { value: "process" },
|
|
||||||
dbaseOutputFormat: { value: "influxdb" },
|
|
||||||
|
|
||||||
enableLog: { value: false },
|
enableLog: { value: false },
|
||||||
logLevel: { value: "error" },
|
logLevel: { value: "error" },
|
||||||
@@ -41,7 +49,7 @@
|
|||||||
outputs: 3,
|
outputs: 3,
|
||||||
inputLabels: ["input"],
|
inputLabels: ["input"],
|
||||||
outputLabels: ["process", "dbase", "parent"],
|
outputLabels: ["process", "dbase", "parent"],
|
||||||
icon: "font-awesome/fa-recycle",
|
icon: "font-awesome/fa-flask",
|
||||||
label: function() {
|
label: function() {
|
||||||
return this.name || "Reactor";
|
return this.name || "Reactor";
|
||||||
},
|
},
|
||||||
@@ -235,23 +243,6 @@
|
|||||||
<label for="node-input-timeStep"><i class="fa fa-tag"></i> Time step [s]</label>
|
<label for="node-input-timeStep"><i class="fa fa-tag"></i> Time step [s]</label>
|
||||||
<input type="text" id="node-input-timeStep" placeholder="s">
|
<input type="text" id="node-input-timeStep" placeholder="s">
|
||||||
</div>
|
</div>
|
||||||
<h3>Output Formats</h3>
|
|
||||||
<div class="form-row">
|
|
||||||
<label for="node-input-processOutputFormat"><i class="fa fa-random"></i> Process Output</label>
|
|
||||||
<select id="node-input-processOutputFormat" style="width:60%;">
|
|
||||||
<option value="process">process</option>
|
|
||||||
<option value="json">json</option>
|
|
||||||
<option value="csv">csv</option>
|
|
||||||
</select>
|
|
||||||
</div>
|
|
||||||
<div class="form-row">
|
|
||||||
<label for="node-input-dbaseOutputFormat"><i class="fa fa-database"></i> Database Output</label>
|
|
||||||
<select id="node-input-dbaseOutputFormat" style="width:60%;">
|
|
||||||
<option value="influxdb">influxdb</option>
|
|
||||||
<option value="json">json</option>
|
|
||||||
<option value="csv">csv</option>
|
|
||||||
</select>
|
|
||||||
</div>
|
|
||||||
|
|
||||||
<!-- Logger fields injected here -->
|
<!-- Logger fields injected here -->
|
||||||
<div id="logger-fields-placeholder"></div>
|
<div id="logger-fields-placeholder"></div>
|
||||||
|
|||||||
@@ -1,5 +1,4 @@
|
|||||||
const { Reactor_CSTR, Reactor_PFR } = require('./specificClass.js');
|
const { Reactor_CSTR, Reactor_PFR } = require('./specificClass.js');
|
||||||
const { configManager } = require('generalFunctions');
|
|
||||||
|
|
||||||
|
|
||||||
class nodeClass {
|
class nodeClass {
|
||||||
@@ -49,15 +48,14 @@ class nodeClass {
|
|||||||
case "Dispersion":
|
case "Dispersion":
|
||||||
this.source.setDispersion = msg;
|
this.source.setDispersion = msg;
|
||||||
break;
|
break;
|
||||||
case 'registerChild': {
|
case 'registerChild':
|
||||||
// Register this node as a parent of the child node
|
// Register this node as a parent of the child node
|
||||||
const childId = msg.payload;
|
const childId = msg.payload;
|
||||||
const childObj = this.RED.nodes.getNode(childId);
|
const childObj = this.RED.nodes.getNode(childId);
|
||||||
this.source.childRegistrationUtils.registerChild(childObj.source, msg.positionVsParent);
|
this.source.childRegistrationUtils.registerChild(childObj.source, msg.positionVsParent);
|
||||||
break;
|
break;
|
||||||
}
|
|
||||||
default:
|
default:
|
||||||
this.source.logger.warn(`Unknown topic: ${msg.topic}`);
|
console.log("Unknown topic: " + msg.topic);
|
||||||
}
|
}
|
||||||
|
|
||||||
if (done) {
|
if (done) {
|
||||||
@@ -71,10 +69,20 @@ class nodeClass {
|
|||||||
* @param {object} uiConfig Config set in UI in node-red
|
* @param {object} uiConfig Config set in UI in node-red
|
||||||
*/
|
*/
|
||||||
_loadConfig(uiConfig) {
|
_loadConfig(uiConfig) {
|
||||||
const cfgMgr = new configManager();
|
this.config = {
|
||||||
|
general: {
|
||||||
// Build config: base sections + reactor-specific domain config
|
name: uiConfig.name || this.name,
|
||||||
this.config = cfgMgr.buildConfig('reactor', uiConfig, this.node.id, {
|
id: this.node.id,
|
||||||
|
unit: null,
|
||||||
|
logging: {
|
||||||
|
enabled: uiConfig.enableLog,
|
||||||
|
logLevel: uiConfig.logLevel
|
||||||
|
}
|
||||||
|
},
|
||||||
|
functionality: {
|
||||||
|
positionVsParent: uiConfig.positionVsParent || 'atEquipment', // Default to 'atEquipment' if not specified
|
||||||
|
softwareType: "reactor" // should be set in config manager
|
||||||
|
},
|
||||||
reactor_type: uiConfig.reactor_type,
|
reactor_type: uiConfig.reactor_type,
|
||||||
volume: parseFloat(uiConfig.volume),
|
volume: parseFloat(uiConfig.volume),
|
||||||
length: parseFloat(uiConfig.length),
|
length: parseFloat(uiConfig.length),
|
||||||
@@ -98,7 +106,7 @@ class nodeClass {
|
|||||||
parseFloat(uiConfig.X_TS_init)
|
parseFloat(uiConfig.X_TS_init)
|
||||||
],
|
],
|
||||||
timeStep: parseFloat(uiConfig.timeStep)
|
timeStep: parseFloat(uiConfig.timeStep)
|
||||||
});
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
/**
|
/**
|
||||||
@@ -129,7 +137,7 @@ class nodeClass {
|
|||||||
new_reactor = new Reactor_PFR(this.config);
|
new_reactor = new Reactor_PFR(this.config);
|
||||||
break;
|
break;
|
||||||
default:
|
default:
|
||||||
throw new Error(`Unknown reactor type: ${this.config.reactor_type}`);
|
console.warn("Unknown reactor type: " + uiConfig.reactor_type);
|
||||||
}
|
}
|
||||||
|
|
||||||
this.source = new_reactor; // protect from reassignment
|
this.source = new_reactor; // protect from reassignment
|
||||||
|
|||||||
@@ -19,7 +19,7 @@ class ASM3 {
|
|||||||
nu_NO: 0.5, // anoxic reduction factor [-]
|
nu_NO: 0.5, // anoxic reduction factor [-]
|
||||||
K_O: 0.2, // saturation constant S_0 [g O2 m-3]
|
K_O: 0.2, // saturation constant S_0 [g O2 m-3]
|
||||||
K_NO: 0.5, // saturation constant S_NO [g NO3-N m-3]
|
K_NO: 0.5, // saturation constant S_NO [g NO3-N m-3]
|
||||||
K_S: 10.0, // saturation constant S_s [g COD m-3]
|
K_S: 10., // saturation constant S_s [g COD m-3]
|
||||||
K_STO: 0.1, // saturation constant X_STO [g X_STO g-1 X_H]
|
K_STO: 0.1, // saturation constant X_STO [g X_STO g-1 X_H]
|
||||||
mu_H_max: 3., // maximum specific growth rate [d-1]
|
mu_H_max: 3., // maximum specific growth rate [d-1]
|
||||||
K_NH: 0.01, // saturation constant S_NH3 [g NH3-N m-3]
|
K_NH: 0.01, // saturation constant S_NH3 [g NH3-N m-3]
|
||||||
@@ -171,7 +171,7 @@ class ASM3 {
|
|||||||
compute_rates(state, T = 20) {
|
compute_rates(state, T = 20) {
|
||||||
// state: S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
|
// state: S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
|
||||||
const rates = Array(12);
|
const rates = Array(12);
|
||||||
const [S_O, , S_S, S_NH, , S_NO, S_HCO, , X_S, X_H, X_STO, X_A] = state;
|
const [S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS] = state;
|
||||||
const { k_H, K_X, k_STO, nu_NO, K_O, K_NO, K_S, K_STO, mu_H_max, K_NH, K_HCO, b_H_O, b_H_NO, b_STO_O, b_STO_NO, mu_A_max, K_A_NH, K_A_O, K_A_HCO, b_A_O, b_A_NO } = this.kin_params;
|
const { k_H, K_X, k_STO, nu_NO, K_O, K_NO, K_S, K_STO, mu_H_max, K_NH, K_HCO, b_H_O, b_H_NO, b_STO_O, b_STO_NO, mu_A_max, K_A_NH, K_A_O, K_A_HCO, b_A_O, b_A_NO } = this.kin_params;
|
||||||
const { theta_H, theta_STO, theta_mu_H, theta_b_H_O, theta_b_H_NO, theta_b_STO_O, theta_b_STO_NO, theta_mu_A, theta_b_A_O, theta_b_A_NO } = this.temp_params;
|
const { theta_H, theta_STO, theta_mu_H, theta_b_H_O, theta_b_H_NO, theta_b_STO_O, theta_b_STO_NO, theta_mu_A, theta_b_A_O, theta_b_A_NO } = this.temp_params;
|
||||||
|
|
||||||
|
|||||||
@@ -171,7 +171,7 @@ class ASM3 {
|
|||||||
compute_rates(state, T = 20) {
|
compute_rates(state, T = 20) {
|
||||||
// state: S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
|
// state: S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
|
||||||
const rates = Array(12);
|
const rates = Array(12);
|
||||||
const [S_O, , S_S, S_NH, , S_NO, S_HCO, , X_S, X_H, X_STO, X_A] = state;
|
const [S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS] = state;
|
||||||
const { k_H, K_X, k_STO, nu_NO, K_O, K_NO, K_S, K_STO, mu_H_max, K_NH, K_HCO, b_H_O, b_H_NO, b_STO_O, b_STO_NO, mu_A_max, K_A_NH, K_A_O, K_A_HCO, b_A_O, b_A_NO } = this.kin_params;
|
const { k_H, K_X, k_STO, nu_NO, K_O, K_NO, K_S, K_STO, mu_H_max, K_NH, K_HCO, b_H_O, b_H_NO, b_STO_O, b_STO_NO, mu_A_max, K_A_NH, K_A_O, K_A_HCO, b_A_O, b_A_NO } = this.kin_params;
|
||||||
const { theta_H, theta_STO, theta_mu_H, theta_b_H_O, theta_b_H_NO, theta_b_STO_O, theta_b_STO_NO, theta_mu_A, theta_b_A_O, theta_b_A_NO } = this.temp_params;
|
const { theta_H, theta_STO, theta_mu_H, theta_b_H_O, theta_b_H_NO, theta_b_STO_O, theta_b_STO_NO, theta_mu_A, theta_b_A_O, theta_b_A_NO } = this.temp_params;
|
||||||
|
|
||||||
|
|||||||
@@ -1,7 +1,7 @@
|
|||||||
const ASM3 = require('./reaction_modules/asm3_class.js');
|
const ASM3 = require('./reaction_modules/asm3_class.js');
|
||||||
const { create, all, isArray } = require('mathjs');
|
const { create, all, isArray } = require('mathjs');
|
||||||
const { assertNoNaN } = require('./utils.js');
|
const { assertNoNaN } = require('./utils.js');
|
||||||
const { childRegistrationUtils, logger, MeasurementContainer, POSITIONS } = require('generalFunctions');
|
const { childRegistrationUtils, logger, MeasurementContainer } = require('generalFunctions');
|
||||||
const EventEmitter = require('events');
|
const EventEmitter = require('events');
|
||||||
|
|
||||||
const mathConfig = {
|
const mathConfig = {
|
||||||
@@ -126,6 +126,7 @@ class Reactor {
|
|||||||
position = measurement.config.functionality.positionVsParent;
|
position = measurement.config.functionality.positionVsParent;
|
||||||
}
|
}
|
||||||
const measurementType = measurement.config.asset.type;
|
const measurementType = measurement.config.asset.type;
|
||||||
|
const key = `${measurementType}_${position}`;
|
||||||
const eventName = `${measurementType}.measured.${position}`;
|
const eventName = `${measurementType}.measured.${position}`;
|
||||||
|
|
||||||
// Register event listener for measurement updates
|
// Register event listener for measurement updates
|
||||||
@@ -159,11 +160,11 @@ class Reactor {
|
|||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
_updateMeasurement(measurementType, value, position, _context) {
|
_updateMeasurement(measurementType, value, position, context) {
|
||||||
this.logger.debug(`---------------------- updating ${measurementType} ------------------ `);
|
this.logger.debug(`---------------------- updating ${measurementType} ------------------ `);
|
||||||
switch (measurementType) {
|
switch (measurementType) {
|
||||||
case "temperature":
|
case "temperature":
|
||||||
if (position == POSITIONS.AT_EQUIPMENT) {
|
if (position == "atEquipment") {
|
||||||
this.temperature = value;
|
this.temperature = value;
|
||||||
}
|
}
|
||||||
break;
|
break;
|
||||||
@@ -319,15 +320,14 @@ class Reactor_PFR extends Reactor {
|
|||||||
return stateNew;
|
return stateNew;
|
||||||
}
|
}
|
||||||
|
|
||||||
_updateMeasurement(measurementType, value, position, _context) {
|
_updateMeasurement(measurementType, value, position, context) {
|
||||||
switch(measurementType) {
|
switch(measurementType) {
|
||||||
case "quantity (oxygen)": {
|
case "quantity (oxygen)":
|
||||||
let grid_pos = Math.round(position / this.config.length * this.n_x);
|
let grid_pos = Math.round(position / this.config.length * this.n_x);
|
||||||
this.state[grid_pos][S_O_INDEX] = value; // naive approach for reconciling measurements and simulation
|
this.state[grid_pos][S_O_INDEX] = value; // naive approach for reconciling measurements and simulation
|
||||||
break;
|
break;
|
||||||
}
|
|
||||||
default:
|
default:
|
||||||
super._updateMeasurement(measurementType, value, position, _context);
|
super._updateMeasurement(measurementType, value, position, context);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|||||||
@@ -1,346 +0,0 @@
|
|||||||
/**
|
|
||||||
* Tests for reactor specificClass (domain logic).
|
|
||||||
*
|
|
||||||
* Two reactor classes are exported: Reactor_CSTR and Reactor_PFR.
|
|
||||||
* Both extend a base Reactor class.
|
|
||||||
*
|
|
||||||
* Key methods tested:
|
|
||||||
* - _calcOTR: oxygen transfer rate calculation
|
|
||||||
* - _arrayClip2Zero: clip negative values to zero
|
|
||||||
* - setInfluent / getEffluent: influent/effluent data flow
|
|
||||||
* - setOTR: external OTR override
|
|
||||||
* - tick (CSTR): forward Euler state update
|
|
||||||
* - tick (PFR): finite difference state update
|
|
||||||
* - registerChild: dispatches to measurement / reactor handlers
|
|
||||||
*/
|
|
||||||
|
|
||||||
const { Reactor_CSTR, Reactor_PFR } = require('../src/specificClass');
|
|
||||||
|
|
||||||
// --------------- helpers ---------------
|
|
||||||
|
|
||||||
const NUM_SPECIES = 13;
|
|
||||||
|
|
||||||
function makeCSTRConfig(overrides = {}) {
|
|
||||||
return {
|
|
||||||
general: {
|
|
||||||
name: 'TestCSTR',
|
|
||||||
id: 'cstr-test-1',
|
|
||||||
logging: { enabled: false, logLevel: 'error' },
|
|
||||||
},
|
|
||||||
functionality: {
|
|
||||||
softwareType: 'reactor',
|
|
||||||
positionVsParent: 'atEquipment',
|
|
||||||
},
|
|
||||||
volume: 1000,
|
|
||||||
n_inlets: 1,
|
|
||||||
kla: 240,
|
|
||||||
timeStep: 1, // 1 second
|
|
||||||
initialState: new Array(NUM_SPECIES).fill(1.0),
|
|
||||||
...overrides,
|
|
||||||
};
|
|
||||||
}
|
|
||||||
|
|
||||||
function makePFRConfig(overrides = {}) {
|
|
||||||
return {
|
|
||||||
general: {
|
|
||||||
name: 'TestPFR',
|
|
||||||
id: 'pfr-test-1',
|
|
||||||
logging: { enabled: false, logLevel: 'error' },
|
|
||||||
},
|
|
||||||
functionality: {
|
|
||||||
softwareType: 'reactor',
|
|
||||||
positionVsParent: 'atEquipment',
|
|
||||||
},
|
|
||||||
volume: 200,
|
|
||||||
length: 10,
|
|
||||||
resolution_L: 10,
|
|
||||||
n_inlets: 1,
|
|
||||||
kla: 240,
|
|
||||||
alpha: 0.5,
|
|
||||||
timeStep: 1,
|
|
||||||
initialState: new Array(NUM_SPECIES).fill(0.1),
|
|
||||||
...overrides,
|
|
||||||
};
|
|
||||||
}
|
|
||||||
|
|
||||||
// --------------- CSTR tests ---------------
|
|
||||||
|
|
||||||
describe('Reactor_CSTR', () => {
|
|
||||||
|
|
||||||
describe('constructor / initialization', () => {
|
|
||||||
it('should create an instance and set state from initialState', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
expect(r).toBeDefined();
|
|
||||||
expect(r.state).toEqual(new Array(NUM_SPECIES).fill(1.0));
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should initialize Fs and Cs_in arrays based on n_inlets', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig({ n_inlets: 3 }));
|
|
||||||
expect(r.Fs).toHaveLength(3);
|
|
||||||
expect(r.Cs_in).toHaveLength(3);
|
|
||||||
expect(r.Fs.every(v => v === 0)).toBe(true);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should store volume from config', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig({ volume: 500 }));
|
|
||||||
expect(r.volume).toBe(500);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should initialize temperature to 20', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
expect(r.temperature).toBe(20);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('_calcOTR()', () => {
|
|
||||||
let r;
|
|
||||||
beforeAll(() => { r = new Reactor_CSTR(makeCSTRConfig({ kla: 240 })); });
|
|
||||||
|
|
||||||
it('should return a positive value when S_O < saturation', () => {
|
|
||||||
const otr = r._calcOTR(0, 20);
|
|
||||||
expect(otr).toBeGreaterThan(0);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should return approximately zero when S_O equals saturation', () => {
|
|
||||||
// S_O_sat at T=20: 14.652 - 4.1022e-1*20 + 7.9910e-3*400 + 7.7774e-5*8000
|
|
||||||
const T = 20;
|
|
||||||
const S_O_sat = 14.652 - 4.1022e-1 * T + 7.9910e-3 * T * T + 7.7774e-5 * T * T * T;
|
|
||||||
const otr = r._calcOTR(S_O_sat, T);
|
|
||||||
expect(otr).toBeCloseTo(0, 5);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should return a negative value when S_O > saturation (supersaturated)', () => {
|
|
||||||
const otr = r._calcOTR(100, 20);
|
|
||||||
expect(otr).toBeLessThan(0);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should use T=20 as default temperature', () => {
|
|
||||||
const otr1 = r._calcOTR(0);
|
|
||||||
const otr2 = r._calcOTR(0, 20);
|
|
||||||
expect(otr1).toBe(otr2);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('_arrayClip2Zero()', () => {
|
|
||||||
let r;
|
|
||||||
beforeAll(() => { r = new Reactor_CSTR(makeCSTRConfig()); });
|
|
||||||
|
|
||||||
it('should clip negative values to zero', () => {
|
|
||||||
expect(r._arrayClip2Zero([-5, 3, -1, 0, 7])).toEqual([0, 3, 0, 0, 7]);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should leave all-positive arrays unchanged', () => {
|
|
||||||
expect(r._arrayClip2Zero([1, 2, 3])).toEqual([1, 2, 3]);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should handle nested arrays (2D)', () => {
|
|
||||||
const result = r._arrayClip2Zero([[-1, 2], [3, -4]]);
|
|
||||||
expect(result).toEqual([[0, 2], [3, 0]]);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should handle a single scalar', () => {
|
|
||||||
expect(r._arrayClip2Zero(-5)).toBe(0);
|
|
||||||
expect(r._arrayClip2Zero(5)).toBe(5);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('setInfluent / getEffluent', () => {
|
|
||||||
it('should store influent data via setter', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig({ n_inlets: 2 }));
|
|
||||||
const input = {
|
|
||||||
payload: {
|
|
||||||
inlet: 0,
|
|
||||||
F: 100,
|
|
||||||
C: new Array(NUM_SPECIES).fill(5),
|
|
||||||
},
|
|
||||||
};
|
|
||||||
r.setInfluent = input;
|
|
||||||
expect(r.Fs[0]).toBe(100);
|
|
||||||
expect(r.Cs_in[0]).toEqual(new Array(NUM_SPECIES).fill(5));
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should return effluent with the sum of Fs and the current state', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
r.Fs[0] = 50;
|
|
||||||
const eff = r.getEffluent;
|
|
||||||
expect(eff.topic).toBe('Fluent');
|
|
||||||
expect(eff.payload.F).toBe(50);
|
|
||||||
expect(eff.payload.C).toEqual(r.state);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('setOTR', () => {
|
|
||||||
it('should set the OTR value', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig({ kla: NaN }));
|
|
||||||
r.setOTR = { payload: 42 };
|
|
||||||
expect(r.OTR).toBe(42);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('tick()', () => {
|
|
||||||
it('should return a new state array of correct length', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
const result = r.tick(0.001);
|
|
||||||
expect(result).toHaveLength(NUM_SPECIES);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not produce NaN values', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
r.Fs[0] = 10;
|
|
||||||
r.Cs_in[0] = new Array(NUM_SPECIES).fill(5);
|
|
||||||
const result = r.tick(0.001);
|
|
||||||
result.forEach(v => expect(Number.isNaN(v)).toBe(false));
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not produce negative concentrations', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
// Run multiple ticks
|
|
||||||
for (let i = 0; i < 100; i++) {
|
|
||||||
r.tick(0.001);
|
|
||||||
}
|
|
||||||
r.state.forEach(v => expect(v).toBeGreaterThanOrEqual(0));
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should reach steady state with zero flow (concentrations change only via reaction)', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
// No inflow
|
|
||||||
const initial = [...r.state];
|
|
||||||
r.tick(0.0001);
|
|
||||||
// State should have changed due to reaction/OTR
|
|
||||||
const changed = r.state.some((v, i) => v !== initial[i]);
|
|
||||||
expect(changed).toBe(true);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('registerChild()', () => {
|
|
||||||
it('should not throw for "measurement" software type', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
// Passing null child will trigger warn but not crash
|
|
||||||
expect(() => r.registerChild(null, 'measurement')).not.toThrow();
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not throw for "reactor" software type', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
expect(() => r.registerChild(null, 'reactor')).not.toThrow();
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not throw for unknown software type', () => {
|
|
||||||
const r = new Reactor_CSTR(makeCSTRConfig());
|
|
||||||
expect(() => r.registerChild(null, 'unknown')).not.toThrow();
|
|
||||||
});
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
// --------------- PFR tests ---------------
|
|
||||||
|
|
||||||
describe('Reactor_PFR', () => {
|
|
||||||
|
|
||||||
describe('constructor / initialization', () => {
|
|
||||||
it('should create an instance with 2D state grid', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
expect(r).toBeDefined();
|
|
||||||
expect(r.state).toHaveLength(10); // resolution_L = 10
|
|
||||||
expect(r.state[0]).toHaveLength(NUM_SPECIES);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should compute d_x = length / n_x', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ length: 10, resolution_L: 5 }));
|
|
||||||
expect(r.d_x).toBe(2);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should compute cross-sectional area A = volume / length', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ volume: 200, length: 10 }));
|
|
||||||
expect(r.A).toBe(20);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should initialize D (dispersion) to 0', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
expect(r.D).toBe(0);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should create derivative operators of correct size', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ resolution_L: 8 }));
|
|
||||||
expect(r.D_op).toHaveLength(8);
|
|
||||||
expect(r.D_op[0]).toHaveLength(8);
|
|
||||||
expect(r.D2_op).toHaveLength(8);
|
|
||||||
expect(r.D2_op[0]).toHaveLength(8);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('setDispersion', () => {
|
|
||||||
it('should set the axial dispersion value', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
r.setDispersion = { payload: 0.5 };
|
|
||||||
expect(r.D).toBe(0.5);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('tick()', () => {
|
|
||||||
it('should return a 2D state grid of correct dimensions', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
r.D = 0.01;
|
|
||||||
const result = r.tick(0.0001);
|
|
||||||
expect(result).toHaveLength(10);
|
|
||||||
expect(result[0]).toHaveLength(NUM_SPECIES);
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not produce NaN values with small time step and dispersion', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
r.D = 0.01;
|
|
||||||
r.Fs[0] = 10;
|
|
||||||
r.Cs_in[0] = new Array(NUM_SPECIES).fill(5);
|
|
||||||
const result = r.tick(0.0001);
|
|
||||||
result.forEach(row => {
|
|
||||||
row.forEach(v => expect(Number.isNaN(v)).toBe(false));
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should not produce negative concentrations', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
r.D = 0.01;
|
|
||||||
for (let i = 0; i < 10; i++) {
|
|
||||||
r.tick(0.0001);
|
|
||||||
}
|
|
||||||
r.state.forEach(row => {
|
|
||||||
row.forEach(v => expect(v).toBeGreaterThanOrEqual(0));
|
|
||||||
});
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('_applyBoundaryConditions()', () => {
|
|
||||||
it('should apply Neumann BC at outlet (last = second to last)', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ resolution_L: 5 }));
|
|
||||||
const state = Array.from({ length: 5 }, () => new Array(NUM_SPECIES).fill(1));
|
|
||||||
state[3] = new Array(NUM_SPECIES).fill(7);
|
|
||||||
r._applyBoundaryConditions(state);
|
|
||||||
// outlet BC: state[4] = state[3]
|
|
||||||
expect(state[4]).toEqual(new Array(NUM_SPECIES).fill(7));
|
|
||||||
});
|
|
||||||
|
|
||||||
it('should apply Neumann BC at inlet when no flow', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ resolution_L: 5 }));
|
|
||||||
r.Fs[0] = 0;
|
|
||||||
const state = Array.from({ length: 5 }, () => new Array(NUM_SPECIES).fill(1));
|
|
||||||
state[1] = new Array(NUM_SPECIES).fill(3);
|
|
||||||
r._applyBoundaryConditions(state);
|
|
||||||
// No flow: state[0] = state[1]
|
|
||||||
expect(state[0]).toEqual(new Array(NUM_SPECIES).fill(3));
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('_arrayClip2Zero() (inherited)', () => {
|
|
||||||
it('should clip 2D arrays correctly', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig());
|
|
||||||
const result = r._arrayClip2Zero([[-1, 2], [3, -4]]);
|
|
||||||
expect(result).toEqual([[0, 2], [3, 0]]);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
|
|
||||||
describe('_calcOTR() (inherited)', () => {
|
|
||||||
it('should work the same as in CSTR', () => {
|
|
||||||
const r = new Reactor_PFR(makePFRConfig({ kla: 240 }));
|
|
||||||
const otr = r._calcOTR(0, 20);
|
|
||||||
expect(otr).toBeGreaterThan(0);
|
|
||||||
});
|
|
||||||
});
|
|
||||||
});
|
|
||||||
Reference in New Issue
Block a user