P11.6 wiki regen + Phase 10 private-test rewrites where applicable

For all 11 nodes with auto-gen markers: wiki/Home.md sections 5 (topic
contract) and 9 (data model) regenerated via npm run wiki:all. New
Unit column shows '<measure> (default <unit>)' for declared topics,
'—' otherwise. Effect column now uses descriptor.description (P11.2
field) overriding the generic per-prefix fallback.

For rotatingMachine + reactor: Phase 10 test rewrites — 3 + 8 files
moved off private nodeClass internals (_attachInputHandler, _commands,
_pendingExtras, _registerChild, _tick, etc.) to the public
BaseNodeAdapter surface (node.handlers.input, node.source.*).
+6 / +7 net new tests.

Co-Authored-By: Claude Opus 4.7 (1M context) <noreply@anthropic.com>

.gitignore

@@ -1,136 +1,136 @@
-# Logs
+# Logs
-logs
+logs
-*.log
+*.log
-npm-debug.log*
+npm-debug.log*
-yarn-debug.log*
+yarn-debug.log*
-yarn-error.log*
+yarn-error.log*
-lerna-debug.log*
+lerna-debug.log*
-.pnpm-debug.log*
+.pnpm-debug.log*
-
+
-# Diagnostic reports (https://nodejs.org/api/report.html)
+# Diagnostic reports (https://nodejs.org/api/report.html)
-report.[0-9]*.[0-9]*.[0-9]*.[0-9]*.json
+report.[0-9]*.[0-9]*.[0-9]*.[0-9]*.json
-
+
-# Runtime data
+# Runtime data
-pids
+pids
-*.pid
+*.pid
-*.seed
+*.seed
-*.pid.lock
+*.pid.lock
-
+
-# Directory for instrumented libs generated by jscoverage/JSCover
+# Directory for instrumented libs generated by jscoverage/JSCover
-lib-cov
+lib-cov
-
+
-# Coverage directory used by tools like istanbul
+# Coverage directory used by tools like istanbul
-coverage
+coverage
-*.lcov
+*.lcov
-
+
-# nyc test coverage
+# nyc test coverage
-.nyc_output
+.nyc_output
-
+
-# Grunt intermediate storage (https://gruntjs.com/creating-plugins#storing-task-files)
+# Grunt intermediate storage (https://gruntjs.com/creating-plugins#storing-task-files)
-.grunt
+.grunt
-
+
-# Bower dependency directory (https://bower.io/)
+# Bower dependency directory (https://bower.io/)
-bower_components
+bower_components
-
+
-# node-waf configuration
+# node-waf configuration
-.lock-wscript
+.lock-wscript
-
+
-# Compiled binary addons (https://nodejs.org/api/addons.html)
+# Compiled binary addons (https://nodejs.org/api/addons.html)
-build/Release
+build/Release
-
+
-# Dependency directories
+# Dependency directories
-node_modules/
+node_modules/
-jspm_packages/
+jspm_packages/
-
+
-# Snowpack dependency directory (https://snowpack.dev/)
+# Snowpack dependency directory (https://snowpack.dev/)
-web_modules/
+web_modules/
-
+
-# TypeScript cache
+# TypeScript cache
-*.tsbuildinfo
+*.tsbuildinfo
-
+
-# Optional npm cache directory
+# Optional npm cache directory
-.npm
+.npm
-
+
-# Optional eslint cache
+# Optional eslint cache
-.eslintcache
+.eslintcache
-
+
-# Optional stylelint cache
+# Optional stylelint cache
-.stylelintcache
+.stylelintcache
-
+
-# Microbundle cache
+# Microbundle cache
-.rpt2_cache/
+.rpt2_cache/
-.rts2_cache_cjs/
+.rts2_cache_cjs/
-.rts2_cache_es/
+.rts2_cache_es/
-.rts2_cache_umd/
+.rts2_cache_umd/
-
+
-# Optional REPL history
+# Optional REPL history
-.node_repl_history
+.node_repl_history
-
+
-# Output of 'npm pack'
+# Output of 'npm pack'
-*.tgz
+*.tgz
-
+
-# Yarn Integrity file
+# Yarn Integrity file
-.yarn-integrity
+.yarn-integrity
-
+
-# dotenv environment variable files
+# dotenv environment variable files
-.env
+.env
-.env.development.local
+.env.development.local
-.env.test.local
+.env.test.local
-.env.production.local
+.env.production.local
-.env.local
+.env.local
-
+
-# parcel-bundler cache (https://parceljs.org/)
+# parcel-bundler cache (https://parceljs.org/)
-.cache
+.cache
-.parcel-cache
+.parcel-cache
-
+
-# Next.js build output
+# Next.js build output
-.next
+.next
-out
+out
-
+
-# Nuxt.js build / generate output
+# Nuxt.js build / generate output
-.nuxt
+.nuxt
-dist
+dist
-
+
-# Gatsby files
+# Gatsby files
-.cache/
+.cache/
-# Comment in the public line in if your project uses Gatsby and not Next.js
+# Comment in the public line in if your project uses Gatsby and not Next.js
-# https://nextjs.org/blog/next-9-1#public-directory-support
+# https://nextjs.org/blog/next-9-1#public-directory-support
-# public
+# public
-
+
-# vuepress build output
+# vuepress build output
-.vuepress/dist
+.vuepress/dist
-
+
-# vuepress v2.x temp and cache directory
+# vuepress v2.x temp and cache directory
-.temp
+.temp
-.cache
+.cache
-
+
-# vitepress build output
+# vitepress build output
-**/.vitepress/dist
+**/.vitepress/dist
-
+
-# vitepress cache directory
+# vitepress cache directory
-**/.vitepress/cache
+**/.vitepress/cache
-
+
-# Docusaurus cache and generated files
+# Docusaurus cache and generated files
-.docusaurus
+.docusaurus
-
+
-# Serverless directories
+# Serverless directories
-.serverless/
+.serverless/
-
+
-# FuseBox cache
+# FuseBox cache
-.fusebox/
+.fusebox/
-
+
-# DynamoDB Local files
+# DynamoDB Local files
-.dynamodb/
+.dynamodb/
-
+
-# TernJS port file
+# TernJS port file
-.tern-port
+.tern-port
-
+
-# Stores VSCode versions used for testing VSCode extensions
+# Stores VSCode versions used for testing VSCode extensions
-.vscode-test
+.vscode-test
-
+
-# yarn v2
+# yarn v2
-.yarn/cache
+.yarn/cache
-.yarn/unplugged
+.yarn/unplugged
-.yarn/build-state.yml
+.yarn/build-state.yml
-.yarn/install-state.gz
+.yarn/install-state.gz
-.pnp.*
+.pnp.*

CLAUDE.md

@@ -0,0 +1,23 @@
+# reactor — Claude Code context
+
+Biological reactor with ASM kinetics.
+Part of the [EVOLV](https://gitea.wbd-rd.nl/RnD/EVOLV) wastewater-automation platform.
+
+## S88 classification
+
+| Level | Colour | Placement lane |
+|---|---|---|
+| **Unit** | `#50a8d9` | L4 |
+
+## Flow layout rules
+
+When wiring this node into a multi-node demo or production flow, follow the
+placement rule set in the **EVOLV superproject**:
+
+> `.claude/rules/node-red-flow-layout.md` (in the EVOLV repo root)
+
+Key points for this node:
+- Place on lane **L4** (x-position per the lane table in the rule).
+- Stack same-level siblings vertically.
+- Parent/children sit on adjacent lanes (children one lane left, parent one lane right).
+- Wrap in a Node-RED group box coloured `#50a8d9` (Unit).

CONTRACT.md

@@ -0,0 +1,51 @@
+# reactor — Contract
+
+Hand-maintained for Phase 6; the `## Inputs` table is generated from
+`src/commands/index.js` (see Phase 9 generator). Keep ≤ 80 lines.
+
+## Inputs (msg.topic on Port 0)
+
+| Canonical | Aliases (deprecated) | Payload | Effect |
+|---|---|---|---|
+| `data.clock` | `clock` | `msg.timestamp` (ms since epoch) | Calls `source.updateState(timestamp)` — advances the ASM kinetics integrator by `n_iter` time steps that fit between `currentTime` and the supplied timestamp (scaled by `speedUpFactor`). |
+| `data.fluent` | `Fluent` | `{ inlet: number, F: number, C: number[13] }` | Writes the per-inlet flow rate (`F`, m³/d) and concentration vector (`C`) into `engine.Fs[inlet]` / `engine.Cs_in[inlet]`. |
+| `data.otr` | `OTR` | numeric | Sets the externally-supplied oxygen transfer rate (used when `kla` is NaN). |
+| `data.temperature` | `Temperature` | numeric or `{ value: number }` | Sets `engine.temperature` (°C). Non-numeric payloads are warned and ignored. |
+| `data.dispersion` | `Dispersion` | numeric | PFR only — sets the axial dispersion coefficient `D` (m²/d). |
+| `child.register` | `registerChild` | child node id (string) | Looks up the sibling node via `RED.nodes.getNode(id)` and delegates to `source.childRegistrationUtils.registerChild` with `msg.positionVsParent`. |
+
+Aliases log a one-time deprecation warning the first time they fire.
+
+## Outputs (msg.topic on Port 0/1/2)
+
+- **Port 0 (process):** every tick emits the engine's effluent:
+  `{ topic: 'Fluent', payload: { inlet: 0, F, C: number[13] }, timestamp }`.
+  For a PFR an additional `{ topic: 'GridProfile', payload: { grid, n_x, d_x, length, species, timestamp } }`
+  message goes out on the same port before the effluent.
+- **Port 1 (InfluxDB telemetry):** formatted via `outputUtils.formatMsg(..., 'influxdb')`
+  from `getOutput()` — carries `flow_total`, `temperature`, and one field per ASM3
+  species (`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`).
+- **Port 2 (registration):** at startup the node sends one
+  `{ topic: 'child.register', payload: <node.id>, positionVsParent, distance }`
+  to its parent.
+
+## Events emitted by `source.emitter`
+
+- `stateChange` — fires after every `updateState()` that advances the integrator.
+  Payload is the new `currentTime` (ms since epoch). Downstream reactors register
+  via `child.register` and subscribe to this event to pull the upstream
+  effluent on each advance.
+- `output-changed` — base notification fired by `updateState()` so the
+  BaseNodeAdapter pipeline pushes outputs (currently used only as a heartbeat;
+  effluent is emitted directly from the periodic tick).
+
+## Children accepted
+
+- `measurement` — subscribes to `<type>.measured.<position>` on the child's
+  `measurements.emitter`. Recognised reconciliations: `temperature.measured.atEquipment`
+  writes `engine.temperature`; PFR additionally honours
+  `quantity (oxygen).measured.<distance>` to reconcile dissolved-oxygen
+  concentration into the nearest grid cell.
+- `reactor` — registers as the upstream reactor; the downstream `updateState`
+  pulls the upstream effluent into `Fs[0]` / `Cs_in[0]` before integrating.

LICENSE

@@ -1,190 +1,190 @@
-EUROPEAN UNION PUBLIC LICENCE v. 1.2
+EUROPEAN UNION PUBLIC LICENCE v. 1.2
-EUPL © the European Union 2007, 2016
+EUPL © the European Union 2007, 2016
-
+
-This European Union Public Licence (the ‘EUPL’) applies to the Work (as defined below) which is provided under the
+This European Union Public Licence (the ‘EUPL’) applies to the Work (as defined below) which is provided under the
-terms of this Licence. Any use of the Work, other than as authorised under this Licence is prohibited (to the extent such
+terms of this Licence. Any use of the Work, other than as authorised under this Licence is prohibited (to the extent such
-use is covered by a right of the copyright holder of the Work).
+use is covered by a right of the copyright holder of the Work).
-The Work is provided under the terms of this Licence when the Licensor (as defined below) has placed the following
+The Work is provided under the terms of this Licence when the Licensor (as defined below) has placed the following
-notice immediately following the copyright notice for the Work:
+notice immediately following the copyright notice for the Work:
-                          Licensed under the EUPL
+                          Licensed under the EUPL
-or has expressed by any other means his willingness to license under the EUPL.
+or has expressed by any other means his willingness to license under the EUPL.
-
+
-1.Definitions
+1.Definitions
-In this Licence, the following terms have the following meaning:
+In this Licence, the following terms have the following meaning:
-— ‘The Licence’:this Licence.
+— ‘The Licence’:this Licence.
-— ‘The Original Work’:the work or software distributed or communicated by the Licensor under this Licence, available
+— ‘The Original Work’:the work or software distributed or communicated by the Licensor under this Licence, available
-as Source Code and also as Executable Code as the case may be.
+as Source Code and also as Executable Code as the case may be.
-— ‘Derivative Works’:the works or software that could be created by the Licensee, based upon the Original Work or
+— ‘Derivative Works’:the works or software that could be created by the Licensee, based upon the Original Work or
-modifications thereof. This Licence does not define the extent of modification or dependence on the Original Work
+modifications thereof. This Licence does not define the extent of modification or dependence on the Original Work
-required in order to classify a work as a Derivative Work; this extent is determined by copyright law applicable in
+required in order to classify a work as a Derivative Work; this extent is determined by copyright law applicable in
-the country mentioned in Article 15.
+the country mentioned in Article 15.
-— ‘The Work’:the Original Work or its Derivative Works.
+— ‘The Work’:the Original Work or its Derivative Works.
-— ‘The Source Code’:the human-readable form of the Work which is the most convenient for people to study and
+— ‘The Source Code’:the human-readable form of the Work which is the most convenient for people to study and
-modify.
+modify.
-— ‘The Executable Code’:any code which has generally been compiled and which is meant to be interpreted by
+— ‘The Executable Code’:any code which has generally been compiled and which is meant to be interpreted by
-a computer as a program.
+a computer as a program.
-— ‘The Licensor’:the natural or legal person that distributes or communicates the Work under the Licence.
+— ‘The Licensor’:the natural or legal person that distributes or communicates the Work under the Licence.
-— ‘Contributor(s)’:any natural or legal person who modifies the Work under the Licence, or otherwise contributes to
+— ‘Contributor(s)’:any natural or legal person who modifies the Work under the Licence, or otherwise contributes to
-the creation of a Derivative Work.
+the creation of a Derivative Work.
-— ‘The Licensee’ or ‘You’:any natural or legal person who makes any usage of the Work under the terms of the
+— ‘The Licensee’ or ‘You’:any natural or legal person who makes any usage of the Work under the terms of the
-Licence.
+Licence.
-— ‘Distribution’ or ‘Communication’:any act of selling, giving, lending, renting, distributing, communicating,
+— ‘Distribution’ or ‘Communication’:any act of selling, giving, lending, renting, distributing, communicating,
-transmitting, or otherwise making available, online or offline, copies of the Work or providing access to its essential
+transmitting, or otherwise making available, online or offline, copies of the Work or providing access to its essential
-functionalities at the disposal of any other natural or legal person.
+functionalities at the disposal of any other natural or legal person.
-
+
-2.Scope of the rights granted by the Licence
+2.Scope of the rights granted by the Licence
-The Licensor hereby grants You a worldwide, royalty-free, non-exclusive, sublicensable licence to do the following, for
+The Licensor hereby grants You a worldwide, royalty-free, non-exclusive, sublicensable licence to do the following, for
-the duration of copyright vested in the Original Work:
+the duration of copyright vested in the Original Work:
-— use the Work in any circumstance and for all usage,
+— use the Work in any circumstance and for all usage,
-— reproduce the Work,
+— reproduce the Work,
-— modify the Work, and make Derivative Works based upon the Work,
+— modify the Work, and make Derivative Works based upon the Work,
-— communicate to the public, including the right to make available or display the Work or copies thereof to the public
+— communicate to the public, including the right to make available or display the Work or copies thereof to the public
-and perform publicly, as the case may be, the Work,
+and perform publicly, as the case may be, the Work,
-— distribute the Work or copies thereof,
+— distribute the Work or copies thereof,
-— lend and rent the Work or copies thereof,
+— lend and rent the Work or copies thereof,
-— sublicense rights in the Work or copies thereof.
+— sublicense rights in the Work or copies thereof.
-Those rights can be exercised on any media, supports and formats, whether now known or later invented, as far as the
+Those rights can be exercised on any media, supports and formats, whether now known or later invented, as far as the
-applicable law permits so.
+applicable law permits so.
-In the countries where moral rights apply, the Licensor waives his right to exercise his moral right to the extent allowed
+In the countries where moral rights apply, the Licensor waives his right to exercise his moral right to the extent allowed
-by law in order to make effective the licence of the economic rights here above listed.
+by law in order to make effective the licence of the economic rights here above listed.
-The Licensor grants to the Licensee royalty-free, non-exclusive usage rights to any patents held by the Licensor, to the
+The Licensor grants to the Licensee royalty-free, non-exclusive usage rights to any patents held by the Licensor, to the
-extent necessary to make use of the rights granted on the Work under this Licence.
+extent necessary to make use of the rights granted on the Work under this Licence.
-
+
-3.Communication of the Source Code
+3.Communication of the Source Code
-The Licensor may provide the Work either in its Source Code form, or as Executable Code. If the Work is provided as
+The Licensor may provide the Work either in its Source Code form, or as Executable Code. If the Work is provided as
-Executable Code, the Licensor provides in addition a machine-readable copy of the Source Code of the Work along with
+Executable Code, the Licensor provides in addition a machine-readable copy of the Source Code of the Work along with
-each copy of the Work that the Licensor distributes or indicates, in a notice following the copyright notice attached to
+each copy of the Work that the Licensor distributes or indicates, in a notice following the copyright notice attached to
-the Work, a repository where the Source Code is easily and freely accessible for as long as the Licensor continues to
+the Work, a repository where the Source Code is easily and freely accessible for as long as the Licensor continues to
-distribute or communicate the Work.
+distribute or communicate the Work.
-
+
-4.Limitations on copyright
+4.Limitations on copyright
-Nothing in this Licence is intended to deprive the Licensee of the benefits from any exception or limitation to the
+Nothing in this Licence is intended to deprive the Licensee of the benefits from any exception or limitation to the
-exclusive rights of the rights owners in the Work, of the exhaustion of those rights or of other applicable limitations
+exclusive rights of the rights owners in the Work, of the exhaustion of those rights or of other applicable limitations
-thereto.
+thereto.
-
+
-5.Obligations of the Licensee
+5.Obligations of the Licensee
-The grant of the rights mentioned above is subject to some restrictions and obligations imposed on the Licensee. Those
+The grant of the rights mentioned above is subject to some restrictions and obligations imposed on the Licensee. Those
-obligations are the following:
+obligations are the following:
-
+
-Attribution right: The Licensee shall keep intact all copyright, patent or trademarks notices and all notices that refer to
+Attribution right: The Licensee shall keep intact all copyright, patent or trademarks notices and all notices that refer to
-the Licence and to the disclaimer of warranties. The Licensee must include a copy of such notices and a copy of the
+the Licence and to the disclaimer of warranties. The Licensee must include a copy of such notices and a copy of the
-Licence with every copy of the Work he/she distributes or communicates. The Licensee must cause any Derivative Work
+Licence with every copy of the Work he/she distributes or communicates. The Licensee must cause any Derivative Work
-to carry prominent notices stating that the Work has been modified and the date of modification.
+to carry prominent notices stating that the Work has been modified and the date of modification.
-
+
-Copyleft clause: If the Licensee distributes or communicates copies of the Original Works or Derivative Works, this
+Copyleft clause: If the Licensee distributes or communicates copies of the Original Works or Derivative Works, this
-Distribution or Communication will be done under the terms of this Licence or of a later version of this Licence unless
+Distribution or Communication will be done under the terms of this Licence or of a later version of this Licence unless
-the Original Work is expressly distributed only under this version of the Licence — for example by communicating
+the Original Work is expressly distributed only under this version of the Licence — for example by communicating
-‘EUPL v. 1.2 only’. The Licensee (becoming Licensor) cannot offer or impose any additional terms or conditions on the
+‘EUPL v. 1.2 only’. The Licensee (becoming Licensor) cannot offer or impose any additional terms or conditions on the
-Work or Derivative Work that alter or restrict the terms of the Licence.
+Work or Derivative Work that alter or restrict the terms of the Licence.
-
+
-Compatibility clause: If the Licensee Distributes or Communicates Derivative Works or copies thereof based upon both
+Compatibility clause: If the Licensee Distributes or Communicates Derivative Works or copies thereof based upon both
-the Work and another work licensed under a Compatible Licence, this Distribution or Communication can be done
+the Work and another work licensed under a Compatible Licence, this Distribution or Communication can be done
-under the terms of this Compatible Licence. For the sake of this clause, ‘Compatible Licence’ refers to the licences listed
+under the terms of this Compatible Licence. For the sake of this clause, ‘Compatible Licence’ refers to the licences listed
-in the appendix attached to this Licence. Should the Licensee's obligations under the Compatible Licence conflict with
+in the appendix attached to this Licence. Should the Licensee's obligations under the Compatible Licence conflict with
-his/her obligations under this Licence, the obligations of the Compatible Licence shall prevail.
+his/her obligations under this Licence, the obligations of the Compatible Licence shall prevail.
-
+
-Provision of Source Code: When distributing or communicating copies of the Work, the Licensee will provide
+Provision of Source Code: When distributing or communicating copies of the Work, the Licensee will provide
-a machine-readable copy of the Source Code or indicate a repository where this Source will be easily and freely available
+a machine-readable copy of the Source Code or indicate a repository where this Source will be easily and freely available
-for as long as the Licensee continues to distribute or communicate the Work.
+for as long as the Licensee continues to distribute or communicate the Work.
-Legal Protection: This Licence does not grant permission to use the trade names, trademarks, service marks, or names
+Legal Protection: This Licence does not grant permission to use the trade names, trademarks, service marks, or names
-of the Licensor, except as required for reasonable and customary use in describing the origin of the Work and
+of the Licensor, except as required for reasonable and customary use in describing the origin of the Work and
-reproducing the content of the copyright notice.
+reproducing the content of the copyright notice.
-
+
-6.Chain of Authorship
+6.Chain of Authorship
-The original Licensor warrants that the copyright in the Original Work granted hereunder is owned by him/her or
+The original Licensor warrants that the copyright in the Original Work granted hereunder is owned by him/her or
-licensed to him/her and that he/she has the power and authority to grant the Licence.
+licensed to him/her and that he/she has the power and authority to grant the Licence.
-Each Contributor warrants that the copyright in the modifications he/she brings to the Work are owned by him/her or
+Each Contributor warrants that the copyright in the modifications he/she brings to the Work are owned by him/her or
-licensed to him/her and that he/she has the power and authority to grant the Licence.
+licensed to him/her and that he/she has the power and authority to grant the Licence.
-Each time You accept the Licence, the original Licensor and subsequent Contributors grant You a licence to their contributions
+Each time You accept the Licence, the original Licensor and subsequent Contributors grant You a licence to their contributions
-to the Work, under the terms of this Licence.
+to the Work, under the terms of this Licence.
-
+
-7.Disclaimer of Warranty
+7.Disclaimer of Warranty
-The Work is a work in progress, which is continuously improved by numerous Contributors. It is not a finished work
+The Work is a work in progress, which is continuously improved by numerous Contributors. It is not a finished work
-and may therefore contain defects or ‘bugs’ inherent to this type of development.
+and may therefore contain defects or ‘bugs’ inherent to this type of development.
-For the above reason, the Work is provided under the Licence on an ‘as is’ basis and without warranties of any kind
+For the above reason, the Work is provided under the Licence on an ‘as is’ basis and without warranties of any kind
-concerning the Work, including without limitation merchantability, fitness for a particular purpose, absence of defects or
+concerning the Work, including without limitation merchantability, fitness for a particular purpose, absence of defects or
-errors, accuracy, non-infringement of intellectual property rights other than copyright as stated in Article 6 of this
+errors, accuracy, non-infringement of intellectual property rights other than copyright as stated in Article 6 of this
-Licence.
+Licence.
-This disclaimer of warranty is an essential part of the Licence and a condition for the grant of any rights to the Work.
+This disclaimer of warranty is an essential part of the Licence and a condition for the grant of any rights to the Work.
-
+
-8.Disclaimer of Liability
+8.Disclaimer of Liability
-Except in the cases of wilful misconduct or damages directly caused to natural persons, the Licensor will in no event be
+Except in the cases of wilful misconduct or damages directly caused to natural persons, the Licensor will in no event be
-liable for any direct or indirect, material or moral, damages of any kind, arising out of the Licence or of the use of the
+liable for any direct or indirect, material or moral, damages of any kind, arising out of the Licence or of the use of the
-Work, including without limitation, damages for loss of goodwill, work stoppage, computer failure or malfunction, loss
+Work, including without limitation, damages for loss of goodwill, work stoppage, computer failure or malfunction, loss
-of data or any commercial damage, even if the Licensor has been advised of the possibility of such damage. However,
+of data or any commercial damage, even if the Licensor has been advised of the possibility of such damage. However,
-the Licensor will be liable under statutory product liability laws as far such laws apply to the Work.
+the Licensor will be liable under statutory product liability laws as far such laws apply to the Work.
-
+
-9.Additional agreements
+9.Additional agreements
-While distributing the Work, You may choose to conclude an additional agreement, defining obligations or services
+While distributing the Work, You may choose to conclude an additional agreement, defining obligations or services
-consistent with this Licence. However, if accepting obligations, You may act only on your own behalf and on your sole
+consistent with this Licence. However, if accepting obligations, You may act only on your own behalf and on your sole
-responsibility, not on behalf of the original Licensor or any other Contributor, and only if You agree to indemnify,
+responsibility, not on behalf of the original Licensor or any other Contributor, and only if You agree to indemnify,
-defend, and hold each Contributor harmless for any liability incurred by, or claims asserted against such Contributor by
+defend, and hold each Contributor harmless for any liability incurred by, or claims asserted against such Contributor by
-the fact You have accepted any warranty or additional liability.
+the fact You have accepted any warranty or additional liability.
-
+
-10.Acceptance of the Licence
+10.Acceptance of the Licence
-The provisions of this Licence can be accepted by clicking on an icon ‘I agree’ placed under the bottom of a window
+The provisions of this Licence can be accepted by clicking on an icon ‘I agree’ placed under the bottom of a window
-displaying the text of this Licence or by affirming consent in any other similar way, in accordance with the rules of
+displaying the text of this Licence or by affirming consent in any other similar way, in accordance with the rules of
-applicable law. Clicking on that icon indicates your clear and irrevocable acceptance of this Licence and all of its terms
+applicable law. Clicking on that icon indicates your clear and irrevocable acceptance of this Licence and all of its terms
-and conditions.
+and conditions.
-Similarly, you irrevocably accept this Licence and all of its terms and conditions by exercising any rights granted to You
+Similarly, you irrevocably accept this Licence and all of its terms and conditions by exercising any rights granted to You
-by Article 2 of this Licence, such as the use of the Work, the creation by You of a Derivative Work or the Distribution
+by Article 2 of this Licence, such as the use of the Work, the creation by You of a Derivative Work or the Distribution
-or Communication by You of the Work or copies thereof.
+or Communication by You of the Work or copies thereof.
-
+
-11.Information to the public
+11.Information to the public
-In case of any Distribution or Communication of the Work by means of electronic communication by You (for example,
+In case of any Distribution or Communication of the Work by means of electronic communication by You (for example,
-by offering to download the Work from a remote location) the distribution channel or media (for example, a website)
+by offering to download the Work from a remote location) the distribution channel or media (for example, a website)
-must at least provide to the public the information requested by the applicable law regarding the Licensor, the Licence
+must at least provide to the public the information requested by the applicable law regarding the Licensor, the Licence
-and the way it may be accessible, concluded, stored and reproduced by the Licensee.
+and the way it may be accessible, concluded, stored and reproduced by the Licensee.
-
+
-12.Termination of the Licence
+12.Termination of the Licence
-The Licence and the rights granted hereunder will terminate automatically upon any breach by the Licensee of the terms
+The Licence and the rights granted hereunder will terminate automatically upon any breach by the Licensee of the terms
-of the Licence.
+of the Licence.
-Such a termination will not terminate the licences of any person who has received the Work from the Licensee under
+Such a termination will not terminate the licences of any person who has received the Work from the Licensee under
-the Licence, provided such persons remain in full compliance with the Licence.
+the Licence, provided such persons remain in full compliance with the Licence.
-
+
-13.Miscellaneous
+13.Miscellaneous
-Without prejudice of Article 9 above, the Licence represents the complete agreement between the Parties as to the
+Without prejudice of Article 9 above, the Licence represents the complete agreement between the Parties as to the
-Work.
+Work.
-If any provision of the Licence is invalid or unenforceable under applicable law, this will not affect the validity or
+If any provision of the Licence is invalid or unenforceable under applicable law, this will not affect the validity or
-enforceability of the Licence as a whole. Such provision will be construed or reformed so as necessary to make it valid
+enforceability of the Licence as a whole. Such provision will be construed or reformed so as necessary to make it valid
-and enforceable.
+and enforceable.
-The European Commission may publish other linguistic versions or new versions of this Licence or updated versions of
+The European Commission may publish other linguistic versions or new versions of this Licence or updated versions of
-the Appendix, so far this is required and reasonable, without reducing the scope of the rights granted by the Licence.
+the Appendix, so far this is required and reasonable, without reducing the scope of the rights granted by the Licence.
-New versions of the Licence will be published with a unique version number.
+New versions of the Licence will be published with a unique version number.
-All linguistic versions of this Licence, approved by the European Commission, have identical value. Parties can take
+All linguistic versions of this Licence, approved by the European Commission, have identical value. Parties can take
-advantage of the linguistic version of their choice.
+advantage of the linguistic version of their choice.
-
+
-14.Jurisdiction
+14.Jurisdiction
-Without prejudice to specific agreement between parties,
+Without prejudice to specific agreement between parties,
-— any litigation resulting from the interpretation of this License, arising between the European Union institutions,
+— any litigation resulting from the interpretation of this License, arising between the European Union institutions,
-bodies, offices or agencies, as a Licensor, and any Licensee, will be subject to the jurisdiction of the Court of Justice
+bodies, offices or agencies, as a Licensor, and any Licensee, will be subject to the jurisdiction of the Court of Justice
-of the European Union, as laid down in article 272 of the Treaty on the Functioning of the European Union,
+of the European Union, as laid down in article 272 of the Treaty on the Functioning of the European Union,
-— any litigation arising between other parties and resulting from the interpretation of this License, will be subject to
+— any litigation arising between other parties and resulting from the interpretation of this License, will be subject to
-the exclusive jurisdiction of the competent court where the Licensor resides or conducts its primary business.
+the exclusive jurisdiction of the competent court where the Licensor resides or conducts its primary business.
-
+
-15.Applicable Law
+15.Applicable Law
-Without prejudice to specific agreement between parties,
+Without prejudice to specific agreement between parties,
-— this Licence shall be governed by the law of the European Union Member State where the Licensor has his seat,
+— this Licence shall be governed by the law of the European Union Member State where the Licensor has his seat,
-resides or has his registered office,
+resides or has his registered office,
-— this licence shall be governed by Belgian law if the Licensor has no seat, residence or registered office inside
+— this licence shall be governed by Belgian law if the Licensor has no seat, residence or registered office inside
-a European Union Member State.
+a European Union Member State.
-
+
-
+
-                                                         Appendix
+                                                         Appendix
-
+
-‘Compatible Licences’ according to Article 5 EUPL are:
+‘Compatible Licences’ according to Article 5 EUPL are:
-— GNU General Public License (GPL) v. 2, v. 3
+— GNU General Public License (GPL) v. 2, v. 3
-— GNU Affero General Public License (AGPL) v. 3
+— GNU Affero General Public License (AGPL) v. 3
-— Open Software License (OSL) v. 2.1, v. 3.0
+— Open Software License (OSL) v. 2.1, v. 3.0
-— Eclipse Public License (EPL) v. 1.0
+— Eclipse Public License (EPL) v. 1.0
-— CeCILL v. 2.0, v. 2.1
+— CeCILL v. 2.0, v. 2.1
-— Mozilla Public Licence (MPL) v. 2
+— Mozilla Public Licence (MPL) v. 2
-— GNU Lesser General Public Licence (LGPL) v. 2.1, v. 3
+— GNU Lesser General Public Licence (LGPL) v. 2.1, v. 3
-— Creative Commons Attribution-ShareAlike v. 3.0 Unported (CC BY-SA 3.0) for works other than software
+— Creative Commons Attribution-ShareAlike v. 3.0 Unported (CC BY-SA 3.0) for works other than software
-— European Union Public Licence (EUPL) v. 1.1, v. 1.2
+— European Union Public Licence (EUPL) v. 1.1, v. 1.2
-— Québec Free and Open-Source Licence — Reciprocity (LiLiQ-R) or Strong Reciprocity (LiLiQ-R+).
+— Québec Free and Open-Source Licence — Reciprocity (LiLiQ-R) or Strong Reciprocity (LiLiQ-R+).
-
+
-The European Commission may update this Appendix to later versions of the above licences without producing
+The European Commission may update this Appendix to later versions of the above licences without producing
-a new version of the EUPL, as long as they provide the rights granted in Article 2 of this Licence and protect the
+a new version of the EUPL, as long as they provide the rights granted in Article 2 of this Licence and protect the
-covered Source Code from exclusive appropriation.
+covered Source Code from exclusive appropriation.
-All other changes or additions to this Appendix require the production of a new EUPL version.
+All other changes or additions to this Appendix require the production of a new EUPL version.

README.md

@@ -1,17 +1,17 @@
-# reactor
+# reactor
-
+
-Reactor: Advanced Hydraulic Tank & Biological Process Simulator
+Reactor: Advanced Hydraulic Tank & Biological Process Simulator
-
+
-A comprehensive reactor class for wastewater treatment simulation featuring plug flow hydraulics, ASM1-ASM3 biological modeling, and multi-sectional concentration tracking. Implements hydraulic retention time calculations, dispersion modeling, and real-time biological reaction kinetics for accurate process simulation.
+A comprehensive reactor class for wastewater treatment simulation featuring plug flow hydraulics, ASM1-ASM3 biological modeling, and multi-sectional concentration tracking. Implements hydraulic retention time calculations, dispersion modeling, and real-time biological reaction kinetics for accurate process simulation.
-
+
-Key Features:
+Key Features:
-
+
-Plug Flow Hydraulics: Multi-section reactor with configurable sectioning factor and dispersion modeling
+Plug Flow Hydraulics: Multi-section reactor with configurable sectioning factor and dispersion modeling
-ASM1 Integration: Complete biological nutrient removal modeling with 13 state variables (COD, nitrogen, phosphorus)
+ASM1 Integration: Complete biological nutrient removal modeling with 13 state variables (COD, nitrogen, phosphorus)
-Dynamic Volume Control: Automatic section management with overflow handling and retention time calculations
+Dynamic Volume Control: Automatic section management with overflow handling and retention time calculations
-Oxygen Transfer: Saturation-limited O2 transfer with Fick's law slowdown effects and solubility curves
+Oxygen Transfer: Saturation-limited O2 transfer with Fick's law slowdown effects and solubility curves
-Real-time Kinetics: Continuous biological reaction rate calculations with configurable time acceleration
+Real-time Kinetics: Continuous biological reaction rate calculations with configurable time acceleration
-Weighted Averaging: Volume-based concentration mixing for accurate mass balance calculations
+Weighted Averaging: Volume-based concentration mixing for accurate mass balance calculations
-Child Registration: Integration with diffuser systems and upstream/downstream reactor networks
+Child Registration: Integration with diffuser systems and upstream/downstream reactor networks
-Supports complex biological treatment train modeling with temperature compensation, sludge calculations, and comprehensive process monitoring for wastewater treatment plant optimization and regulatory compliance.
+Supports complex biological treatment train modeling with temperature compensation, sludge calculations, and comprehensive process monitoring for wastewater treatment plant optimization and regulatory compliance.
-
+

additional_nodes/recirculation-pump.html

@@ -1,57 +1,57 @@
-<script type="text/javascript">
+<script type="text/javascript">
-    RED.nodes.registerType("recirculation-pump", {
+    RED.nodes.registerType("recirculation-pump", {
-        category: "WWTP",
+        category: "WWTP",
-        color: "#e4a363",
+        color: "#e4a363",
-        defaults: {
+        defaults: {
-            name: { value: "" },
+            name: { value: "" },
-            F2: { value: 0, required: true },
+            F2: { value: 0, required: true },
-            inlet: { value: 1, required: true }
+            inlet: { value: 1, required: true }
-        },
+        },
-        inputs: 1,
+        inputs: 1,
-        outputs: 2,
+        outputs: 2,
-        outputLabels: ["Main effluent", "Recirculation effluent"],
+        outputLabels: ["Main effluent", "Recirculation effluent"],
-        icon: "font-awesome/fa-random",
+        icon: "font-awesome/fa-random",
-        label: function() {
+        label: function() {
-            return this.name || "Recirculation pump";
+            return this.name || "Recirculation pump";
-        },
+        },
-        oneditprepare: function() {
+        oneditprepare: function() {
-            $("#node-input-F2").typedInput({
+            $("#node-input-F2").typedInput({
-                type:"num",
+                type:"num",
-                types:["num"]
+                types:["num"]
-            });
+            });
-            $("#node-input-inlet").typedInput({
+            $("#node-input-inlet").typedInput({
-                type:"num",
+                type:"num",
-                types:["num"]
+                types:["num"]
-            });
+            });
-        },
+        },
-        oneditsave: function() {
+        oneditsave: function() {
-            let debit = parseFloat($("#node-input-F2").typedInput("value"));
+            let debit = parseFloat($("#node-input-F2").typedInput("value"));
-            if (isNaN(debit) || debit < 0) {
+            if (isNaN(debit) || debit < 0) {
-                RED.notify("Debit is not set correctly", {type: "error"});
+                RED.notify("Debit is not set correctly", {type: "error"});
-            }
+            }
-            let inlet = parseInt($("#node-input-n_inlets").typedInput("value"));
+            let inlet = parseInt($("#node-input-n_inlets").typedInput("value"));
-            if (inlet < 1) {
+            if (inlet < 1) {
-                RED.notify("Number of inlets not set correctly", {type: "error"});
+                RED.notify("Number of inlets not set correctly", {type: "error"});
-            }
+            }
-        }
+        }
-    });
+    });
-</script>
+</script>
-
+
-<script type="text/html" data-template-name="recirculation-pump">
+<script type="text/html" data-template-name="recirculation-pump">
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-name"><i class="fa fa-tag"></i> Name</label>
+        <label for="node-input-name"><i class="fa fa-tag"></i> Name</label>
-        <input type="text" id="node-input-name" placeholder="Name">
+        <input type="text" id="node-input-name" placeholder="Name">
-    </div>
+    </div>
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-F2"><i class="fa fa-tag"></i> Recirculation debit [m3 d-1]</label>
+        <label for="node-input-F2"><i class="fa fa-tag"></i> Recirculation debit [m3 d-1]</label>
-        <input type="text" id="node-input-F2" placeholder="m3 s-1">
+        <input type="text" id="node-input-F2" placeholder="m3 s-1">
-    </div>
+    </div>
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-inlet"><i class="fa fa-tag"></i> Assigned inlet recirculation</label>
+        <label for="node-input-inlet"><i class="fa fa-tag"></i> Assigned inlet recirculation</label>
-        <input type="text" id="node-input-inlet" placeholder="#">
+        <input type="text" id="node-input-inlet" placeholder="#">
-    </div>
+    </div>
-</script>
+</script>
-
+
-<script type="text/html" data-help-name="recirculation-pump">
+<script type="text/html" data-help-name="recirculation-pump">
-    <p>Recirculation-pump for splitting streams</p>
+    <p>Recirculation-pump for splitting streams</p>
-</script>
+</script>

additional_nodes/recirculation-pump.js

@@ -3,13 +3,12 @@ module.exports = function(RED) {
         RED.nodes.createNode(this, config);
         var node = this;
 
-        let name = config.name;
         let F2 = parseFloat(config.F2);
         const inlet_F2 = parseInt(config.inlet);
 
         node.on('input', function(msg, send, done) {
             switch (msg.topic) {
-                case "Fluent":
+                case "Fluent": {
                     // conserve volume flow debit
                     let F_in = msg.payload.F;
                     let F1 = Math.max(F_in - F2, 0);
@@ -24,6 +23,7 @@ module.exports = function(RED) {
 
                     send([msg_F1, msg_F2]);
                     break;
+                }
                 case "clock":
                     break;
                 default:

additional_nodes/settling-basin.html

@@ -1,57 +1,57 @@
-<script type="text/javascript">
+<script type="text/javascript">
-    RED.nodes.registerType("settling-basin", {
+    RED.nodes.registerType("settling-basin", {
-        category: "WWTP",
+        category: "WWTP",
-        color: "#e4a363",
+        color: "#e4a363",
-        defaults: {
+        defaults: {
-            name: { value: "" },
+            name: { value: "" },
-            TS_set: { value: 0.1, required: true },
+            TS_set: { value: 0.1, required: true },
-            inlet: { value: 1, required: true }
+            inlet: { value: 1, required: true }
-        },
+        },
-        inputs: 1,
+        inputs: 1,
-        outputs: 2,
+        outputs: 2,
-        outputLabels: ["Main effluent", "Sludge effluent"],
+        outputLabels: ["Main effluent", "Sludge effluent"],
-        icon: "font-awesome/fa-random",
+        icon: "font-awesome/fa-random",
-        label: function() {
+        label: function() {
-            return this.name || "Settling basin";
+            return this.name || "Settling basin";
-        },
+        },
-        oneditprepare: function() {
+        oneditprepare: function() {
-            $("#node-input-TS_set").typedInput({
+            $("#node-input-TS_set").typedInput({
-                type:"num",
+                type:"num",
-                types:["num"]
+                types:["num"]
-            });
+            });
-            $("#node-input-inlet").typedInput({
+            $("#node-input-inlet").typedInput({
-                type:"num",
+                type:"num",
-                types:["num"]
+                types:["num"]
-            });
+            });
-        },
+        },
-        oneditsave: function() {
+        oneditsave: function() {
-            let TS_set = parseFloat($("#node-input-TS_set").typedInput("value"));
+            let TS_set = parseFloat($("#node-input-TS_set").typedInput("value"));
-            if (isNaN(TS_set) || TS_set < 0) {
+            if (isNaN(TS_set) || TS_set < 0) {
-                RED.notify("TS is not set correctly", {type: "error"});
+                RED.notify("TS is not set correctly", {type: "error"});
-            }
+            }
-            let inlet = parseInt($("#node-input-n_inlets").typedInput("value"));
+            let inlet = parseInt($("#node-input-n_inlets").typedInput("value"));
-            if (inlet < 1) {
+            if (inlet < 1) {
-                RED.notify("Number of inlets not set correctly", {type: "error"});
+                RED.notify("Number of inlets not set correctly", {type: "error"});
-            }
+            }
-        }
+        }
-    });
+    });
-</script>
+</script>
-
+
-<script type="text/html" data-template-name="settling-basin">
+<script type="text/html" data-template-name="settling-basin">
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-name"><i class="fa fa-tag"></i> Name</label>
+        <label for="node-input-name"><i class="fa fa-tag"></i> Name</label>
-        <input type="text" id="node-input-name" placeholder="Name">
+        <input type="text" id="node-input-name" placeholder="Name">
-    </div>
+    </div>
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-TS_set"><i class="fa fa-tag"></i> Total Solids set point [g m-3]</label>
+        <label for="node-input-TS_set"><i class="fa fa-tag"></i> Total Solids set point [g m-3]</label>
-        <input type="text" id="node-input-TS_set" placeholder="">
+        <input type="text" id="node-input-TS_set" placeholder="">
-    </div>
+    </div>
-    <div class="form-row">
+    <div class="form-row">
-        <label for="node-input-inlet"><i class="fa fa-tag"></i> Assigned inlet return line</label>
+        <label for="node-input-inlet"><i class="fa fa-tag"></i> Assigned inlet return line</label>
-        <input type="text" id="node-input-inlet" placeholder="#">
+        <input type="text" id="node-input-inlet" placeholder="#">
-    </div>
+    </div>
-</script>
+</script>
-
+
-<script type="text/html" data-help-name="settling-basin">
+<script type="text/html" data-help-name="settling-basin">
-    <p>Settling tank</p>
+    <p>Settling tank</p>
-</script>
+</script>

additional_nodes/settling-basin.js

@@ -3,13 +3,12 @@ module.exports = function(RED) {
         RED.nodes.createNode(this, config);
         var node = this;
 
-        let name = config.name;
         let TS_set = parseFloat(config.TS_set);
         const inlet_sludge = parseInt(config.inlet);
 
         node.on('input', function(msg, send, done) {
             switch (msg.topic) {
-                case "Fluent":
+                case "Fluent": {
                     // conserve volume flow debit
                     let F_in = msg.payload.F;
                     let C_in = msg.payload.C;
@@ -41,6 +40,7 @@ module.exports = function(RED) {
 
                     send([msg_F1, msg_F2]);
                     break;
+                }
                 case "clock":
                     break;
                 default:

examples/README.md

@@ -0,0 +1,8 @@
+# reactor Example Flows
+
+Import-ready Node-RED examples for reactor.
+
+## Files
+- basic.flow.json
+- integration.flow.json
+- edge.flow.json

examples/basic.flow.json

@@ -0,0 +1,6 @@
+[
+  {"id":"reactor_basic_tab","type":"tab","label":"reactor basic","disabled":false,"info":"reactor basic example"},
+  {"id":"reactor_basic_node","type":"reactor","z":"reactor_basic_tab","name":"reactor basic","x":420,"y":180,"wires":[["reactor_basic_dbg"]]},
+  {"id":"reactor_basic_inj","type":"inject","z":"reactor_basic_tab","name":"basic trigger","props":[{"p":"topic","vt":"str"},{"p":"payload","vt":"str"}],"topic":"ping","payload":"1","payloadType":"str","x":160,"y":180,"wires":[["reactor_basic_node"]]},
+  {"id":"reactor_basic_dbg","type":"debug","z":"reactor_basic_tab","name":"reactor basic debug","active":true,"tosidebar":true,"console":false,"tostatus":false,"complete":"true","targetType":"full","x":660,"y":180,"wires":[]}
+]

examples/edge.flow.json

@@ -0,0 +1,6 @@
+[
+  {"id":"reactor_edge_tab","type":"tab","label":"reactor edge","disabled":false,"info":"reactor edge example"},
+  {"id":"reactor_edge_node","type":"reactor","z":"reactor_edge_tab","name":"reactor edge","x":420,"y":180,"wires":[["reactor_edge_dbg"]]},
+  {"id":"reactor_edge_inj","type":"inject","z":"reactor_edge_tab","name":"unknown topic","props":[{"p":"topic","vt":"str"},{"p":"payload","vt":"str"}],"topic":"doesNotExist","payload":"x","payloadType":"str","x":170,"y":180,"wires":[["reactor_edge_node"]]},
+  {"id":"reactor_edge_dbg","type":"debug","z":"reactor_edge_tab","name":"reactor edge debug","active":true,"tosidebar":true,"console":false,"tostatus":false,"complete":"true","targetType":"full","x":660,"y":180,"wires":[]}
+]

examples/integration.flow.json

@@ -0,0 +1,6 @@
+[
+  {"id":"reactor_int_tab","type":"tab","label":"reactor integration","disabled":false,"info":"reactor integration example"},
+  {"id":"reactor_int_node","type":"reactor","z":"reactor_int_tab","name":"reactor integration","x":420,"y":180,"wires":[["reactor_int_dbg"]]},
+  {"id":"reactor_int_inj","type":"inject","z":"reactor_int_tab","name":"registerChild","props":[{"p":"topic","vt":"str"},{"p":"payload","vt":"str"}],"topic":"registerChild","payload":"example-child-id","payloadType":"str","x":170,"y":180,"wires":[["reactor_int_node"]]},
+  {"id":"reactor_int_dbg","type":"debug","z":"reactor_int_tab","name":"reactor integration debug","active":true,"tosidebar":true,"console":false,"tostatus":false,"complete":"true","targetType":"full","x":680,"y":180,"wires":[]}
+]

package-lock.json

@@ -1,119 +1,119 @@
-{
+{
-  "name": "reactor",
+  "name": "reactor",
-  "version": "0.0.1",
+  "version": "0.0.1",
-  "lockfileVersion": 3,
+  "lockfileVersion": 3,
-  "requires": true,
+  "requires": true,
-  "packages": {
+  "packages": {
-    "": {
+    "": {
-      "name": "reactor",
+      "name": "reactor",
-      "version": "0.0.1",
+      "version": "0.0.1",
-      "license": "SEE LICENSE",
+      "license": "SEE LICENSE",
-      "dependencies": {
+      "dependencies": {
-        "generalFunctions": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git",
+        "generalFunctions": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git",
-        "mathjs": "^14.5.2"
+        "mathjs": "^14.5.2"
-      }
+      }
-    },
+    },
-    "node_modules/@babel/runtime": {
+    "node_modules/@babel/runtime": {
-      "version": "7.28.4",
+      "version": "7.28.4",
-      "resolved": "https://registry.npmjs.org/@babel/runtime/-/runtime-7.28.4.tgz",
+      "resolved": "https://registry.npmjs.org/@babel/runtime/-/runtime-7.28.4.tgz",
-      "integrity": "sha512-Q/N6JNWvIvPnLDvjlE1OUBLPQHH6l3CltCEsHIujp45zQUSSh8K+gHnaEX45yAT1nyngnINhvWtzN+Nb9D8RAQ==",
+      "integrity": "sha512-Q/N6JNWvIvPnLDvjlE1OUBLPQHH6l3CltCEsHIujp45zQUSSh8K+gHnaEX45yAT1nyngnINhvWtzN+Nb9D8RAQ==",
-      "license": "MIT",
+      "license": "MIT",
-      "engines": {
+      "engines": {
-        "node": ">=6.9.0"
+        "node": ">=6.9.0"
-      }
+      }
-    },
+    },
-    "node_modules/complex.js": {
+    "node_modules/complex.js": {
-      "version": "2.4.2",
+      "version": "2.4.2",
-      "resolved": "https://registry.npmjs.org/complex.js/-/complex.js-2.4.2.tgz",
+      "resolved": "https://registry.npmjs.org/complex.js/-/complex.js-2.4.2.tgz",
-      "integrity": "sha512-qtx7HRhPGSCBtGiST4/WGHuW+zeaND/6Ld+db6PbrulIB1i2Ev/2UPiqcmpQNPSyfBKraC0EOvOKCB5dGZKt3g==",
+      "integrity": "sha512-qtx7HRhPGSCBtGiST4/WGHuW+zeaND/6Ld+db6PbrulIB1i2Ev/2UPiqcmpQNPSyfBKraC0EOvOKCB5dGZKt3g==",
-      "license": "MIT",
+      "license": "MIT",
-      "engines": {
+      "engines": {
-        "node": "*"
+        "node": "*"
-      },
+      },
-      "funding": {
+      "funding": {
-        "type": "github",
+        "type": "github",
-        "url": "https://github.com/sponsors/rawify"
+        "url": "https://github.com/sponsors/rawify"
-      }
+      }
-    },
+    },
-    "node_modules/decimal.js": {
+    "node_modules/decimal.js": {
-      "version": "10.6.0",
+      "version": "10.6.0",
-      "resolved": "https://registry.npmjs.org/decimal.js/-/decimal.js-10.6.0.tgz",
+      "resolved": "https://registry.npmjs.org/decimal.js/-/decimal.js-10.6.0.tgz",
-      "integrity": "sha512-YpgQiITW3JXGntzdUmyUR1V812Hn8T1YVXhCu+wO3OpS4eU9l4YdD3qjyiKdV6mvV29zapkMeD390UVEf2lkUg==",
+      "integrity": "sha512-YpgQiITW3JXGntzdUmyUR1V812Hn8T1YVXhCu+wO3OpS4eU9l4YdD3qjyiKdV6mvV29zapkMeD390UVEf2lkUg==",
-      "license": "MIT"
+      "license": "MIT"
-    },
+    },
-    "node_modules/escape-latex": {
+    "node_modules/escape-latex": {
-      "version": "1.2.0",
+      "version": "1.2.0",
-      "resolved": "https://registry.npmjs.org/escape-latex/-/escape-latex-1.2.0.tgz",
+      "resolved": "https://registry.npmjs.org/escape-latex/-/escape-latex-1.2.0.tgz",
-      "integrity": "sha512-nV5aVWW1K0wEiUIEdZ4erkGGH8mDxGyxSeqPzRNtWP7ataw+/olFObw7hujFWlVjNsaDFw5VZ5NzVSIqRgfTiw==",
+      "integrity": "sha512-nV5aVWW1K0wEiUIEdZ4erkGGH8mDxGyxSeqPzRNtWP7ataw+/olFObw7hujFWlVjNsaDFw5VZ5NzVSIqRgfTiw==",
-      "license": "MIT"
+      "license": "MIT"
-    },
+    },
-    "node_modules/fraction.js": {
+    "node_modules/fraction.js": {
-      "version": "5.3.4",
+      "version": "5.3.4",
-      "resolved": "https://registry.npmjs.org/fraction.js/-/fraction.js-5.3.4.tgz",
+      "resolved": "https://registry.npmjs.org/fraction.js/-/fraction.js-5.3.4.tgz",
-      "integrity": "sha512-1X1NTtiJphryn/uLQz3whtY6jK3fTqoE3ohKs0tT+Ujr1W59oopxmoEh7Lu5p6vBaPbgoM0bzveAW4Qi5RyWDQ==",
+      "integrity": "sha512-1X1NTtiJphryn/uLQz3whtY6jK3fTqoE3ohKs0tT+Ujr1W59oopxmoEh7Lu5p6vBaPbgoM0bzveAW4Qi5RyWDQ==",
-      "license": "MIT",
+      "license": "MIT",
-      "engines": {
+      "engines": {
-        "node": "*"
+        "node": "*"
-      },
+      },
-      "funding": {
+      "funding": {
-        "type": "github",
+        "type": "github",
-        "url": "https://github.com/sponsors/rawify"
+        "url": "https://github.com/sponsors/rawify"
-      }
+      }
-    },
+    },
-    "node_modules/generalFunctions": {
+    "node_modules/generalFunctions": {
-      "version": "1.0.0",
+      "version": "1.0.0",
-      "resolved": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git#efc97d6cd17399391b011298e47e8c1b1599592d",
+      "resolved": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git#efc97d6cd17399391b011298e47e8c1b1599592d",
-      "license": "SEE LICENSE"
+      "license": "SEE LICENSE"
-    },
+    },
-    "node_modules/javascript-natural-sort": {
+    "node_modules/javascript-natural-sort": {
-      "version": "0.7.1",
+      "version": "0.7.1",
-      "resolved": "https://registry.npmjs.org/javascript-natural-sort/-/javascript-natural-sort-0.7.1.tgz",
+      "resolved": "https://registry.npmjs.org/javascript-natural-sort/-/javascript-natural-sort-0.7.1.tgz",
-      "integrity": "sha512-nO6jcEfZWQXDhOiBtG2KvKyEptz7RVbpGP4vTD2hLBdmNQSsCiicO2Ioinv6UI4y9ukqnBpy+XZ9H6uLNgJTlw==",
+      "integrity": "sha512-nO6jcEfZWQXDhOiBtG2KvKyEptz7RVbpGP4vTD2hLBdmNQSsCiicO2Ioinv6UI4y9ukqnBpy+XZ9H6uLNgJTlw==",
-      "license": "MIT"
+      "license": "MIT"
-    },
+    },
-    "node_modules/mathjs": {
+    "node_modules/mathjs": {
-      "version": "14.8.0",
+      "version": "14.8.0",
-      "resolved": "https://registry.npmjs.org/mathjs/-/mathjs-14.8.0.tgz",
+      "resolved": "https://registry.npmjs.org/mathjs/-/mathjs-14.8.0.tgz",
-      "integrity": "sha512-DN4wmAjNzFVJ9vHqpAJ3vX0UF306u/1DgGKh7iVPuAFH19JDRd9NAaQS764MsKbSwDB6uBSkQEmgVmKdgYaCoQ==",
+      "integrity": "sha512-DN4wmAjNzFVJ9vHqpAJ3vX0UF306u/1DgGKh7iVPuAFH19JDRd9NAaQS764MsKbSwDB6uBSkQEmgVmKdgYaCoQ==",
-      "license": "Apache-2.0",
+      "license": "Apache-2.0",
-      "dependencies": {
+      "dependencies": {
-        "@babel/runtime": "^7.26.10",
+        "@babel/runtime": "^7.26.10",
-        "complex.js": "^2.2.5",
+        "complex.js": "^2.2.5",
-        "decimal.js": "^10.4.3",
+        "decimal.js": "^10.4.3",
-        "escape-latex": "^1.2.0",
+        "escape-latex": "^1.2.0",
-        "fraction.js": "^5.2.1",
+        "fraction.js": "^5.2.1",
-        "javascript-natural-sort": "^0.7.1",
+        "javascript-natural-sort": "^0.7.1",
-        "seedrandom": "^3.0.5",
+        "seedrandom": "^3.0.5",
-        "tiny-emitter": "^2.1.0",
+        "tiny-emitter": "^2.1.0",
-        "typed-function": "^4.2.1"
+        "typed-function": "^4.2.1"
-      },
+      },
-      "bin": {
+      "bin": {
-        "mathjs": "bin/cli.js"
+        "mathjs": "bin/cli.js"
-      },
+      },
-      "engines": {
+      "engines": {
-        "node": ">= 18"
+        "node": ">= 18"
-      }
+      }
-    },
+    },
-    "node_modules/seedrandom": {
+    "node_modules/seedrandom": {
-      "version": "3.0.5",
+      "version": "3.0.5",
-      "resolved": "https://registry.npmjs.org/seedrandom/-/seedrandom-3.0.5.tgz",
+      "resolved": "https://registry.npmjs.org/seedrandom/-/seedrandom-3.0.5.tgz",
-      "integrity": "sha512-8OwmbklUNzwezjGInmZ+2clQmExQPvomqjL7LFqOYqtmuxRgQYqOD3mHaU+MvZn5FLUeVxVfQjwLZW/n/JFuqg==",
+      "integrity": "sha512-8OwmbklUNzwezjGInmZ+2clQmExQPvomqjL7LFqOYqtmuxRgQYqOD3mHaU+MvZn5FLUeVxVfQjwLZW/n/JFuqg==",
-      "license": "MIT"
+      "license": "MIT"
-    },
+    },
-    "node_modules/tiny-emitter": {
+    "node_modules/tiny-emitter": {
-      "version": "2.1.0",
+      "version": "2.1.0",
-      "resolved": "https://registry.npmjs.org/tiny-emitter/-/tiny-emitter-2.1.0.tgz",
+      "resolved": "https://registry.npmjs.org/tiny-emitter/-/tiny-emitter-2.1.0.tgz",
-      "integrity": "sha512-NB6Dk1A9xgQPMoGqC5CVXn123gWyte215ONT5Pp5a0yt4nlEoO1ZWeCwpncaekPHXO60i47ihFnZPiRPjRMq4Q==",
+      "integrity": "sha512-NB6Dk1A9xgQPMoGqC5CVXn123gWyte215ONT5Pp5a0yt4nlEoO1ZWeCwpncaekPHXO60i47ihFnZPiRPjRMq4Q==",
-      "license": "MIT"
+      "license": "MIT"
-    },
+    },
-    "node_modules/typed-function": {
+    "node_modules/typed-function": {
-      "version": "4.2.1",
+      "version": "4.2.1",
-      "resolved": "https://registry.npmjs.org/typed-function/-/typed-function-4.2.1.tgz",
+      "resolved": "https://registry.npmjs.org/typed-function/-/typed-function-4.2.1.tgz",
-      "integrity": "sha512-EGjWssW7Tsk4DGfE+5yluuljS1OGYWiI1J6e8puZz9nTMM51Oug8CD5Zo4gWMsOhq5BI+1bF+rWTm4Vbj3ivRA==",
+      "integrity": "sha512-EGjWssW7Tsk4DGfE+5yluuljS1OGYWiI1J6e8puZz9nTMM51Oug8CD5Zo4gWMsOhq5BI+1bF+rWTm4Vbj3ivRA==",
-      "license": "MIT",
+      "license": "MIT",
-      "engines": {
+      "engines": {
-        "node": ">= 18"
+        "node": ">= 18"
-      }
+      }
-    }
+    }
-  }
+  }
-}
+}

package.json

@@ -1,33 +1,36 @@
-{
+{
-  "name": "reactor",
+  "name": "reactor",
-  "version": "0.0.1",
+  "version": "0.0.1",
-  "description": "Implementation of the asm3 model for Node-Red",
+  "description": "Implementation of the asm3 model for Node-Red",
-  "repository": {
+  "repository": {
-    "type": "git",
+    "type": "git",
-    "url": "https://gitea.centraal.wbd-rd.nl/RnD/reactor.git"
+    "url": "https://gitea.centraal.wbd-rd.nl/RnD/reactor.git"
-  },
+  },
-  "keywords": [
+  "keywords": [
-    "asm3",
+    "asm3",
-    "activated sludge",
+    "activated sludge",
-    "wastewater",
+    "wastewater",
-    "biological model",
+    "biological model",
-    "node-red"
+    "node-red"
-  ],
+  ],
-  "license": "SEE LICENSE",
+  "license": "SEE LICENSE",
-  "author": "P.R. van der Wilt",
+  "author": "P.R. van der Wilt",
-  "main": "reactor.js",
+  "main": "reactor.js",
-  "scripts": {
+  "scripts": {
-    "test": "node reactor.js"
+    "test": "node --test test/basic/*.test.js test/integration/*.test.js test/edge/*.test.js",
-  },
+    "wiki:contract":  "node ../generalFunctions/scripts/wikiGen.js contract  ./src/commands/index.js   --write ./wiki/Home.md",
-  "node-red": {
+    "wiki:datamodel": "node ../generalFunctions/scripts/wikiGen.js datamodel ./src/specificClass.js     --write ./wiki/Home.md",
-    "nodes": {
+    "wiki:all":       "npm run wiki:contract && npm run wiki:datamodel"
-      "reactor": "reactor.js",
+  },
-      "recirculation-pump": "additional_nodes/recirculation-pump.js",
+  "node-red": {
-      "settling-basin": "additional_nodes/settling-basin.js"
+    "nodes": {
-    }
+      "reactor": "reactor.js",
-  },
+      "recirculation-pump": "additional_nodes/recirculation-pump.js",
-  "dependencies": {
+      "settling-basin": "additional_nodes/settling-basin.js"
-    "generalFunctions": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git",
+    }
-    "mathjs": "^14.5.2"
+  },
-  }
+  "dependencies": {
-}
+    "generalFunctions": "git+https://gitea.centraal.wbd-rd.nl/RnD/generalFunctions.git",
+    "mathjs": "^14.5.2"
+  }
+}

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 type="text/javascript">
     RED.nodes.registerType("reactor", {
-        category: "WWTP",
+        category: "EVOLV",
-        color: "#c4cce0",
+        color: "#50a8d9",
         defaults: {
             name: { value: "" },
             reactor_type: { value: "CSTR", required: true },
@@ -13,7 +23,7 @@
             alpha: {value: 0},
             n_inlets: { value: 1, required: true},
             kla: { value: null },
-            
+
             S_O_init: { value: 0., required: true },
             S_I_init: { value: 30., required: true },
             S_S_init: { value: 100., required: true },
@@ -29,6 +39,9 @@
             X_TS_init: { value: 125.0009, required: true },
 
             timeStep: { value: 1, required: true },
+            speedUpFactor: { value: 1 },
+            processOutputFormat: { value: "process" },
+            dbaseOutputFormat: { value: "influxdb" },
 
             enableLog: { value: false },
             logLevel: { value: "error" },
@@ -39,7 +52,7 @@
         outputs: 3,
         inputLabels: ["input"],
         outputLabels: ["process", "dbase", "parent"],
-        icon: "font-awesome/fa-recycle",
+        icon: "font-awesome/fa-flask",
         label: function() {
             return this.name || "Reactor";
         },
@@ -105,6 +118,10 @@
                 type:"num",
                 types:["num"]
             })
+            $("#node-input-speedUpFactor").typedInput({
+                type:"num",
+                types:["num"]
+            })
             // Set initial visibility on dialog open
             const initialType = $("#node-input-reactor_type").typedInput("value");
             if (initialType === "CSTR") {
@@ -120,8 +137,8 @@
             }
 
             // save position field
-            if (window.EVOLV?.nodes?.measurement?.positionMenu?.saveEditor) {
+            if (window.EVOLV?.nodes?.reactor?.positionMenu?.saveEditor) {
-                window.EVOLV.nodes.rotatingMachine.positionMenu.saveEditor(this);
+                window.EVOLV.nodes.reactor.positionMenu.saveEditor(this);
             }
 
             let volume = parseFloat($("#node-input-volume").typedInput("value"));
@@ -233,6 +250,27 @@
         <label for="node-input-timeStep"><i class="fa fa-tag"></i> Time step [s]</label>
         <input type="text" id="node-input-timeStep" placeholder="s">
     </div>
+    <div class="form-row">
+        <label for="node-input-speedUpFactor"><i class="fa fa-tag"></i> Speed-up factor</label>
+        <input type="text" id="node-input-speedUpFactor" placeholder="1 = real-time">
+    </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 -->
 	<div id="logger-fields-placeholder"></div>

reactor.js

@@ -1,26 +1,26 @@
-const nameOfNode = "reactor"; // name of the node, should match file name and node type in Node-RED
+const nameOfNode = "reactor"; // name of the node, should match file name and node type in Node-RED
-const nodeClass = require('./src/nodeClass.js'); // node class
+const nodeClass = require('./src/nodeClass.js'); // node class
-const { MenuManager } = require('generalFunctions');
+const { MenuManager } = require('generalFunctions');
-
+
-
+
-module.exports = function (RED) {
+module.exports = function (RED) {
-    // Register the node type
+    // Register the node type
-    RED.nodes.registerType(nameOfNode, function (config) {
+    RED.nodes.registerType(nameOfNode, function (config) {
-        // Initialize the Node-RED node first
+        // Initialize the Node-RED node first
-        RED.nodes.createNode(this, config);
+        RED.nodes.createNode(this, config);
-        // Then create your custom class and attach it
+        // Then create your custom class and attach it
-        this.nodeClass = new nodeClass(config, RED, this, nameOfNode);
+        this.nodeClass = new nodeClass(config, RED, this, nameOfNode);
-    });
+    });
-
+
-    const menuMgr = new MenuManager();
+    const menuMgr = new MenuManager();
-
+
-    // Serve /advancedReactor/menu.js
+    // Serve /advancedReactor/menu.js
-    RED.httpAdmin.get(`/${nameOfNode}/menu.js`, (req, res) => {
+    RED.httpAdmin.get(`/${nameOfNode}/menu.js`, (req, res) => {
-        try {
+        try {
-            const script = menuMgr.createEndpoint(nameOfNode, ['logger', 'position']);
+            const script = menuMgr.createEndpoint(nameOfNode, ['logger', 'position']);
-            res.type('application/javascript').send(script);
+            res.type('application/javascript').send(script);
-        } catch (err) {
+        } catch (err) {
-            res.status(500).send(`// Error generating menu: ${err.message}`);
+            res.status(500).send(`// Error generating menu: ${err.message}`);
-        }
+        }
-    });
+    });
-};
+};

src/commands/handlers.js

@@ -0,0 +1,25 @@
+'use strict';
+
+// Reactor input handlers. Each receives (source, msg, ctx) where source is
+// the Reactor domain and ctx is { node, RED, send, logger }. The handlers
+// either forward to engine setters or drive a synchronous state update.
+
+exports.dataClock = (source, msg) => {
+  source.updateState(msg.timestamp ?? Date.now());
+};
+
+exports.dataFluent     = (source, msg) => { source.setInfluent    = msg; };
+exports.dataOTR        = (source, msg) => { source.setOTR         = msg; };
+exports.dataTemperature = (source, msg) => { source.setTemperature = msg; };
+exports.dataDispersion = (source, msg) => { source.setDispersion  = msg; };
+
+exports.childRegister = (source, msg, ctx) => {
+  const childId = msg.payload;
+  const RED = ctx?.RED;
+  const childObj = RED?.nodes?.getNode?.(childId);
+  if (!childObj || !childObj.source) {
+    source?.logger?.warn?.(`registerChild skipped: missing child/source for id=${childId}`);
+    return;
+  }
+  source.childRegistrationUtils.registerChild(childObj.source, msg.positionVsParent);
+};

src/commands/index.js

@@ -0,0 +1,55 @@
+'use strict';
+
+// reactor command registry. Canonical names follow CONTRACTS.md §1.
+// Legacy names (clock, Fluent, OTR, Temperature, Dispersion, registerChild)
+// stay as aliases — they log a one-time deprecation warning on first use
+// and are removed in Phase 7.
+
+const handlers = require('./handlers');
+
+module.exports = [
+  {
+    topic: 'data.clock',
+    aliases: ['clock'],
+    payloadSchema: { type: 'any' },
+    description: 'Push the simulation clock tick (timestamp / dt) to the ASM solver.',
+    handler: handlers.dataClock,
+  },
+  {
+    topic: 'data.fluent',
+    aliases: ['Fluent'],
+    payloadSchema: { type: 'object' },
+    // Compound payload `{F, C: [...]}` — registry-level units normalisation is
+    // skipped (the handler converts per-field internally).
+    description: 'Push the influent stream (payload: {F: flow m3/h, C: [concentrations mg/L]}).',
+    handler: handlers.dataFluent,
+  },
+  {
+    topic: 'data.otr',
+    aliases: ['OTR'],
+    payloadSchema: { type: 'any' },
+    description: 'Push the current oxygen-transfer rate into the reactor.',
+    handler: handlers.dataOTR,
+  },
+  {
+    topic: 'data.temperature',
+    aliases: ['Temperature'],
+    payloadSchema: { type: 'any' },
+    description: 'Push the current reactor temperature.',
+    handler: handlers.dataTemperature,
+  },
+  {
+    topic: 'data.dispersion',
+    aliases: ['Dispersion'],
+    payloadSchema: { type: 'any' },
+    description: 'Push a dispersion/mixing parameter update.',
+    handler: handlers.dataDispersion,
+  },
+  {
+    topic: 'child.register',
+    aliases: ['registerChild'],
+    payloadSchema: { type: 'any' },
+    description: 'Register a child node (settler / measurement) with this reactor.',
+    handler: handlers.childRegister,
+  },
+];

src/kinetics/baseEngine.js

@@ -0,0 +1,139 @@
+'use strict';
+
+const EventEmitter = require('events');
+const ASM3 = require('../reaction_modules/asm3_class.js');
+const { create, all } = require('mathjs');
+const { childRegistrationUtils, logger, MeasurementContainer, POSITIONS } = require('generalFunctions');
+
+const math = create(all, { matrix: 'Array' });
+
+const S_O_INDEX = 0;
+const NUM_SPECIES = 13;
+
+// Abstract reactor engine. Holds the influent/OTR/temperature state plus
+// the parent-side child registration that the original Reactor class
+// exposed. Concrete CSTR / PFR subclasses provide tick().
+class BaseReactorEngine {
+  constructor(config) {
+    this.config = config;
+    this.logger = new logger(
+      this.config.general.logging.enabled,
+      this.config.general.logging.logLevel,
+      config.general.name,
+    );
+    this.emitter = new EventEmitter();
+    this.measurements = new MeasurementContainer();
+    this.upstreamReactor = null;
+    this.childRegistrationUtils = new childRegistrationUtils(this);
+
+    this.asm = new ASM3();
+    this.volume = config.volume;
+
+    this.Fs = Array(config.n_inlets).fill(0);
+    this.Cs_in = Array.from(Array(config.n_inlets), () => new Array(NUM_SPECIES).fill(0));
+    this.OTR = 0.0;
+    this.temperature = 20;
+
+    this.kla = config.kla;
+    this.currentTime = Date.now();
+    // timeStep stored in days (the integrator uses [d] internally).
+    this.timeStep = (1 / (24 * 60 * 60)) * this.config.timeStep;
+    this.speedUpFactor = config.speedUpFactor ?? 1;
+  }
+
+  set setInfluent(input) {
+    const index_in = input.payload.inlet;
+    this.Fs[index_in] = input.payload.F;
+    this.Cs_in[index_in] = input.payload.C;
+  }
+
+  set setOTR(input) { this.OTR = input.payload; }
+
+  set setTemperature(input) {
+    const p = input?.payload;
+    const raw = (p && typeof p === 'object' && p.value !== undefined) ? p.value : p;
+    const v = Number(raw);
+    if (!Number.isFinite(v)) { this.logger.warn(`Invalid temperature input: ${raw}`); return; }
+    this.temperature = v;
+  }
+
+  get getEffluent() {
+    const last = Array.isArray(this.state.at?.(-1)) ? this.state.at(-1) : this.state;
+    return { topic: 'Fluent', payload: { inlet: 0, F: math.sum(this.Fs), C: last }, timestamp: this.currentTime };
+  }
+
+  get getGridProfile() { return null; }
+
+  _calcOTR(S_O, T = 20.0) {
+    const sat = this._calcOxygenSaturation(T);
+    return this.kla * (sat - S_O);
+  }
+
+  _calcOxygenSaturation(T = 20.0) {
+    return 14.652 - 4.1022e-1 * T + 7.9910e-3 * T * T + 7.7774e-5 * T * T * T;
+  }
+
+  _capDissolvedOxygen(state) {
+    const sat = this._calcOxygenSaturation(this.temperature);
+    const capRow = (row) => {
+      if (!Array.isArray(row)) return row;
+      const next = row.slice();
+      if (Number.isFinite(next[S_O_INDEX])) next[S_O_INDEX] = Math.max(0, Math.min(next[S_O_INDEX], sat));
+      return next;
+    };
+    return (Array.isArray(state) && Array.isArray(state[0])) ? state.map(capRow) : capRow(state);
+  }
+
+  _arrayClip2Zero(arr) {
+    if (Array.isArray(arr)) return arr.map((x) => this._arrayClip2Zero(x));
+    return arr < 0 ? 0 : arr;
+  }
+
+  registerChild(child, softwareType) {
+    switch (softwareType) {
+      case 'measurement': this._connectMeasurement(child); break;
+      case 'reactor':     this._connectReactor(child);     break;
+      default:            this.logger.error(`Unrecognized softwareType: ${softwareType}`);
+    }
+  }
+
+  _connectMeasurement(measurement) {
+    if (!measurement) { this.logger.warn('Invalid measurement provided.'); return; }
+    const fn = measurement.config.functionality;
+    const position = fn.distance !== 'undefined' ? fn.distance : fn.positionVsParent;
+    const measurementType = measurement.config.asset.type;
+    const eventName = `${measurementType}.measured.${position}`;
+    measurement.measurements.emitter.on(eventName, (eventData) => {
+      this.measurements
+        .type(measurementType).variant('measured').position(position)
+        .value(eventData.value, eventData.timestamp, eventData.unit);
+      this._updateMeasurement(measurementType, eventData.value, position, eventData);
+    });
+  }
+
+  _connectReactor(reactor) {
+    if (!reactor) { this.logger.warn('Invalid reactor provided.'); return; }
+    this.upstreamReactor = reactor;
+    reactor.emitter.on('stateChange', (data) => this.updateState(data));
+  }
+
+  _updateMeasurement(measurementType, value, position) {
+    if (measurementType === 'temperature' && position === POSITIONS.AT_EQUIPMENT) {
+      this.temperature = value;
+      return;
+    }
+    this.logger.error(`Type '${measurementType}' not recognized for measured update.`);
+  }
+
+  updateState(newTime = Date.now()) {
+    const day2ms = 1000 * 60 * 60 * 24;
+    if (this.upstreamReactor) this.setInfluent = this.upstreamReactor.getEffluent;
+    const n_iter = Math.floor(this.speedUpFactor * (newTime - this.currentTime) / (this.timeStep * day2ms));
+    if (!n_iter) return;
+    for (let n = 0; n < n_iter; n += 1) this.tick(this.timeStep);
+    this.currentTime += (n_iter * this.timeStep * day2ms) / this.speedUpFactor;
+    this.emitter.emit('stateChange', this.currentTime);
+  }
+}
+
+module.exports = { BaseReactorEngine, math, S_O_INDEX, NUM_SPECIES };

src/kinetics/cstr.js

@@ -0,0 +1,27 @@
+'use strict';
+
+const { BaseReactorEngine, math, S_O_INDEX, NUM_SPECIES } = require('./baseEngine.js');
+
+class Reactor_CSTR extends BaseReactorEngine {
+  constructor(config) {
+    super(config);
+    this.state = config.initialState;
+  }
+
+  // Forward Euler step over `time_step` days.
+  tick(time_step) {
+    const inflow  = math.multiply(math.divide([this.Fs], this.volume), this.Cs_in)[0];
+    const outflow = math.multiply(-1 * math.sum(this.Fs) / this.volume, this.state);
+    const reaction = this.asm.compute_dC(this.state, this.temperature);
+    const transfer = Array(NUM_SPECIES).fill(0.0);
+    transfer[S_O_INDEX] = isNaN(this.kla)
+      ? this.OTR
+      : this._calcOTR(this.state[S_O_INDEX], this.temperature);
+
+    const dC_total = math.multiply(math.add(inflow, outflow, reaction, transfer), time_step);
+    this.state = this._capDissolvedOxygen(this._arrayClip2Zero(math.add(this.state, dC_total)));
+    return this.state;
+  }
+}
+
+module.exports = Reactor_CSTR;

src/kinetics/pfr.js

@@ -0,0 +1,132 @@
+'use strict';
+
+const { assertNoNaN } = require('../utils.js');
+const { BaseReactorEngine, math, S_O_INDEX, NUM_SPECIES } = require('./baseEngine.js');
+
+class Reactor_PFR extends BaseReactorEngine {
+  constructor(config) {
+    super(config);
+
+    this.length = config.length;
+    this.n_x = config.resolution_L;
+    this.d_x = this.length / this.n_x;
+    this.A = this.volume / this.length;
+    this.alpha = config.alpha;
+
+    this.state = Array.from(Array(this.n_x), () => config.initialState.slice());
+    this.D = 0.0;
+
+    this.D_op  = this._makeDoperator(true, true);
+    this.D2_op = this._makeD2operator();
+    assertNoNaN(this.D_op,  'Derivative operator');
+    assertNoNaN(this.D2_op, 'Second derivative operator');
+  }
+
+  get getGridProfile() {
+    return {
+      grid: this.state.map((row) => row.slice()),
+      n_x: this.n_x,
+      d_x: this.d_x,
+      length: this.length,
+      species: ['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'],
+      timestamp: this.currentTime,
+    };
+  }
+
+  set setDispersion(input) { this.D = input.payload; }
+
+  updateState(newTime) {
+    super.updateState(newTime);
+    const Pe_local = (this.d_x * math.sum(this.Fs)) / (this.D * this.A);
+    const Co_D = (this.D * this.timeStep) / (this.d_x * this.d_x);
+    if (Pe_local >= 2)  this.logger.warn(`Local Peclet number (${Pe_local}) is too high! Increase reactor resolution.`);
+    if (Co_D >= 0.5)    this.logger.warn(`Courant number (${Co_D}) is too high! Reduce time step size.`);
+  }
+
+  // Explicit finite-difference step over `time_step` days.
+  tick(time_step) {
+    const dispersion = math.multiply(this.D / (this.d_x * this.d_x), this.D2_op, this.state);
+    const advection  = math.multiply(-1 * math.sum(this.Fs) / (this.A * this.d_x), this.D_op, this.state);
+    const reaction   = this.state.map((slice) => this.asm.compute_dC(slice, this.temperature));
+    const transfer   = Array.from(Array(this.n_x), () => new Array(NUM_SPECIES).fill(0));
+
+    const klaIsNaN = isNaN(this.kla);
+    for (let i = 1; i < this.n_x - 1; i += 1) {
+      const otr = klaIsNaN ? this.OTR : this._calcOTR(this.state[i][S_O_INDEX], this.temperature);
+      transfer[i][S_O_INDEX] = otr * this.n_x / (this.n_x - 2);
+    }
+
+    const dC_total = math.multiply(math.add(dispersion, advection, reaction, transfer), time_step);
+    const stateNew = math.add(this.state, dC_total);
+    this._applyBoundaryConditions(stateNew);
+    this.state = this._capDissolvedOxygen(this._arrayClip2Zero(stateNew));
+    return stateNew;
+  }
+
+  _updateMeasurement(measurementType, value, position, context) {
+    if (measurementType === 'quantity (oxygen)') {
+      if (!Number.isFinite(position) || !Number.isFinite(value) || this.config.length <= 0) {
+        this.logger.warn(`Ignoring oxygen measurement update with invalid data (position=${position}, value=${value}).`);
+        return;
+      }
+      const rawIndex = Math.round((position / this.config.length) * this.n_x);
+      const grid_pos = Math.max(0, Math.min(this.n_x - 1, rawIndex));
+      this.state[grid_pos][S_O_INDEX] = value;
+      return;
+    }
+    super._updateMeasurement(measurementType, value, position, context);
+  }
+
+  // Generalised Danckwerts at inlet when flow > 0; Neumann (no-flux) at outlet
+  // and at inlet when there is no flow.
+  _applyBoundaryConditions(state) {
+    if (math.sum(this.Fs) > 0) {
+      const BC_C_in = math.multiply(1 / math.sum(this.Fs), [this.Fs], this.Cs_in)[0];
+      const BC_disp = ((1 - this.alpha) * this.D * this.A) / (math.sum(this.Fs) * this.d_x);
+      state[0] = math.multiply(1 / (1 + BC_disp), math.add(BC_C_in, math.multiply(BC_disp, state[1])));
+    } else {
+      state[0] = state[1];
+    }
+    state[this.n_x - 1] = state[this.n_x - 2];
+  }
+
+  _makeDoperator(central = false, higher_order = false) {
+    if (higher_order) {
+      if (!central) throw new Error('Upwind higher order method not implemented! Use central scheme instead.');
+      const I = math.resize(math.diag(Array(this.n_x).fill(1 / 12), -2),  [this.n_x, this.n_x]);
+      const A = math.resize(math.diag(Array(this.n_x).fill(-2 / 3), -1), [this.n_x, this.n_x]);
+      const B = math.resize(math.diag(Array(this.n_x).fill(2 / 3),  1),  [this.n_x, this.n_x]);
+      const C = math.resize(math.diag(Array(this.n_x).fill(-1 / 12), 2), [this.n_x, this.n_x]);
+      const D = math.add(I, A, B, C);
+      // Preserve the pre-refactor aliasing: D[1] = NearBoundary; NearBoundary.reverse()
+      // mutates D[1] in place; then D[n_x-2] = -1 * NearBoundary uses the reversed view.
+      const nb = Array(this.n_x).fill(0.0);
+      nb[0] = -1 / 4; nb[1] = -5 / 6; nb[2] = 3 / 2; nb[3] = -1 / 2; nb[4] = 1 / 12;
+      D[1] = nb;
+      nb.reverse();
+      D[this.n_x - 2] = math.multiply(-1, nb);
+      D[0] = Array(this.n_x).fill(0);
+      D[this.n_x - 1] = Array(this.n_x).fill(0);
+      return D;
+    }
+    const I = math.resize(math.diag(Array(this.n_x).fill(1 / (1 + central)), central),  [this.n_x, this.n_x]);
+    const A = math.resize(math.diag(Array(this.n_x).fill(-1 / (1 + central)), -1),       [this.n_x, this.n_x]);
+    const D = math.add(I, A);
+    D[0] = Array(this.n_x).fill(0);
+    D[this.n_x - 1] = Array(this.n_x).fill(0);
+    return D;
+  }
+
+  _makeD2operator() {
+    const I = math.diag(Array(this.n_x).fill(-2), 0);
+    const A = math.resize(math.diag(Array(this.n_x).fill(1),  1),  [this.n_x, this.n_x]);
+    const B = math.resize(math.diag(Array(this.n_x).fill(1), -1),  [this.n_x, this.n_x]);
+    const D2 = math.add(I, A, B);
+    D2[0] = Array(this.n_x).fill(0);
+    D2[this.n_x - 1] = Array(this.n_x).fill(0);
+    return D2;
+  }
+}
+
+module.exports = Reactor_PFR;

src/nodeClass.js

@@ -1,165 +1,53 @@
-const { Reactor_CSTR, Reactor_PFR } = require('./specificClass.js');
+'use strict';
 
+const { BaseNodeAdapter } = require('generalFunctions');
+const Reactor = require('./specificClass.js');
+const commands = require('./commands');
 
-class nodeClass {
+const SPECIES = ['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'];
-    * Node-RED node class for advanced-reactor.
-    * @param {object} uiConfig - Node-RED node configuration
-    * @param {object} RED - Node-RED runtime API
-    * @param {object} nodeInstance - Node-RED node instance
-    * @param {string} nameOfNode - Name of the node
-    */
-    constructor(uiConfig, RED, nodeInstance, nameOfNode) {
-        // Preserve RED reference for HTTP endpoints if needed
-        this.node = nodeInstance;
-        this.RED = RED;
-        this.name = nameOfNode;
-        this.source = null;
 
-        this._loadConfig(uiConfig)
+class nodeClass extends BaseNodeAdapter {
-        this._setupClass();
+  static DomainClass = Reactor;
+  static commands = commands;
+  // Tick-driven: ASM kinetics integrate over wall-clock time. The engine's
+  // updateState computes how many internal Euler/FD steps fit in the elapsed
+  // ms; without a periodic tick the integrator never advances.
+  static tickInterval = 1000;
+  static statusInterval = 1000;
 
-        this._attachInputHandler();
+  buildDomainConfig(uiConfig) {
-        this._registerChild();
+    const initialState = {};
-        this._startTickLoop();
+    for (const k of SPECIES) initialState[k] = parseFloat(uiConfig[`${k}_init`]);
-        this._attachCloseHandler();
+    return {
-    }
+      reactor: {
+        reactor_type: uiConfig.reactor_type,
+        volume:       parseFloat(uiConfig.volume),
+        length:       parseFloat(uiConfig.length),
+        resolution_L: parseInt(uiConfig.resolution_L, 10),
+        alpha:        parseFloat(uiConfig.alpha),
+        n_inlets:     parseInt(uiConfig.n_inlets, 10),
+        kla:          parseFloat(uiConfig.kla),
+        timeStep:     parseFloat(uiConfig.timeStep),
+        speedUpFactor: Number(uiConfig.speedUpFactor) || 1,
+      },
+      initialState,
+    };
+  }
 
-    /**
+  // The kinetics engine drives Port-0 effluent + grid-profile shapes that
-     * Handle node-red input messages
+  // don't fit BaseNodeAdapter's delta-compressed payload. Override the
-     */
+  // periodic emission so the Fluent / GridProfile contract is preserved.
-    _attachInputHandler() {
+  _emitOutputs() {
-        this.node.on('input', (msg, send, done) => {
+    const src = this.source;
-
+    if (!src?.engine) return;
-            switch (msg.topic) {
+    src.updateState(Date.now());
-                case "clock":
+    const grid = src.getGridProfile;
-                    this.source.updateState(msg.timestamp);
+    if (grid) this.node.send([{ topic: 'GridProfile', payload: grid }, null, null]);
-                    send([msg, null, null]);
+    const raw = src.getOutput();
-                    break;
+    const influx = this._output.formatMsg(raw, src.config || this.config, 'influxdb');
-                case "Fluent":
+    this.node.send([src.getEffluent, influx, null]);
-                    this.source.setInfluent = msg;
+  }
-                    break;
-                case "OTR":
-                    this.source.setOTR = msg;
-                    break;
-                case "Temperature":
-                    this.source.setTemperature = msg;
-                    break;
-                case "Dispersion":
-                    this.source.setDispersion = msg;
-                    break;
-                case 'registerChild':
-                    // Register this node as a parent of the child node
-                    const childId = msg.payload;
-                    const childObj = this.RED.nodes.getNode(childId); 
-                    this.source.childRegistrationUtils.registerChild(childObj.source, msg.positionVsParent);
-                    break;
-                default:
-                    console.log("Unknown topic: " + msg.topic);
-            }
-
-            if (done) {
-                done();
-            }
-        });
-    }
-
-    /**
-     * Parse node configuration
-     * @param {object} uiConfig Config set in UI in node-red
-     */
-    _loadConfig(uiConfig) {
-        this.config = {
-            general: {
-                name: uiConfig.name || this.name,
-                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,
-            volume: parseFloat(uiConfig.volume),
-            length: parseFloat(uiConfig.length),
-            resolution_L: parseInt(uiConfig.resolution_L),
-            alpha: parseFloat(uiConfig.alpha),
-            n_inlets: parseInt(uiConfig.n_inlets),
-            kla: parseFloat(uiConfig.kla),
-            initialState: [
-                parseFloat(uiConfig.S_O_init),
-                parseFloat(uiConfig.S_I_init),
-                parseFloat(uiConfig.S_S_init),
-                parseFloat(uiConfig.S_NH_init),
-                parseFloat(uiConfig.S_N2_init),
-                parseFloat(uiConfig.S_NO_init),
-                parseFloat(uiConfig.S_HCO_init),
-                parseFloat(uiConfig.X_I_init),
-                parseFloat(uiConfig.X_S_init),
-                parseFloat(uiConfig.X_H_init),
-                parseFloat(uiConfig.X_STO_init),
-                parseFloat(uiConfig.X_A_init),
-                parseFloat(uiConfig.X_TS_init)
-            ],
-            timeStep: parseFloat(uiConfig.timeStep)
-        }
-    }
-
-    /**
-     * Register this node as a child upstream and downstream.
-     * Delayed to avoid Node-RED startup race conditions.
-     */
-    _registerChild() {
-        setTimeout(() => {
-        this.node.send([
-            null,
-            null,
-            { topic: 'registerChild', payload: this.node.id, positionVsParent: this.config?.functionality?.positionVsParent || 'atEquipment' }
-        ]);
-        }, 100);
-    }
-
-    /**
-     * Setup reactor class based on config
-     */
-    _setupClass() {
-        let new_reactor;
-
-        switch (this.config.reactor_type) {
-            case "CSTR":
-                new_reactor = new Reactor_CSTR(this.config);
-                break;
-            case "PFR":
-                new_reactor = new Reactor_PFR(this.config);
-                break;
-            default:
-                console.warn("Unknown reactor type: " + uiConfig.reactor_type);
-        }
-
-        this.source = new_reactor; // protect from reassignment
-        this.node.source = this.source;
-    }
-
-    _startTickLoop() {
-        setTimeout(() => {
-            this._tickInterval = setInterval(() => this._tick(), 1000);
-        }, 1000);
-    }
-
-    _tick(){
-        this.node.send([this.source.getEffluent, null, null]);
-    }
-
-    _attachCloseHandler() {
-        this.node.on('close', (done) => {
-            clearInterval(this._tickInterval);
-            done();
-        });
-    }
 }
 
-module.exports = nodeClass;
+module.exports = nodeClass;

src/reaction_modules/asm3_class Koch.js

@@ -13,7 +13,7 @@ class ASM3 {
     this.kin_params = {
       // Hydrolysis
       k_H:      9.,   // hydrolysis rate constant       [g X_S g-1 X_H d-1]
-      K_X:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H] 
+      K_X:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H]
       // Heterotrophs
       k_STO:    12.,  // storage rate constant          [g S_S g-1 X_H d-1]
       nu_NO:    0.5,  // anoxic reduction factor        [-]
@@ -30,7 +30,7 @@ class ASM3 {
       b_STO_NO: 0.15, // anoxic respitation rate X_STO  [d-1]
       // Autotrophs
       mu_A_max: 1.3,  // maximum specific growth rate   [d-1]
-      K_A_NH:   1.4,  // saturation constant S_NH3      [g NH3-N m-3] 
+      K_A_NH:   1.4,  // saturation constant S_NH3      [g NH3-N m-3]
       K_A_O:    0.5,  // saturation constant S_0        [g O2 m-3]
       K_A_HCO:  0.5,  // saturation constant S_HCO      [mole HCO3 m-3]
       b_A_O:    0.20, // aerobic respiration rate       [d-1]
@@ -132,7 +132,7 @@ class ASM3 {
   /**
    * Computes the inverse Monod equation rate value for a given concentration and half-saturation constant. Used for inhibition.
    * @param {number} c - Concentration of reaction species.
-   * @param {number} K - Half-saturation constant for the reaction species. 
+   * @param {number} K - Half-saturation constant for the reaction species.
    * @returns {number} - Inverse Monod equation rate value for the given concentration and half-saturation constant.
    */
   _inv_monod(c, K) {
@@ -171,7 +171,7 @@ class ASM3 {
   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
     const rates = Array(12);
-    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 [S_O, , S_S, S_NH, , S_NO, S_HCO, , X_S, X_H, X_STO, X_A] = 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 { 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;
 
@@ -187,12 +187,12 @@ class ASM3 {
     rates[6] = this._arrhenius(b_H_NO, theta_b_H_NO, T) * this._inv_monod(S_O, K_O) * this._monod(S_NO, K_NO) * X_H;
     rates[7] = this._arrhenius(b_STO_O, theta_b_STO_O, T) * this._monod(S_O, K_O) * X_H;
     rates[8] = this._arrhenius(b_STO_NO, theta_b_STO_NO, T) * this._inv_monod(S_O, K_O) * this._monod(S_NO, K_NO) * X_STO;
-    
+
     // Autotrophs
     rates[9] = this._arrhenius(mu_A_max, theta_mu_A, T) * this._monod(S_O, K_A_O) * this._monod(S_NH, K_A_NH) * this._monod(S_HCO, K_A_HCO) * X_A;
     rates[10] = this._arrhenius(b_A_O, theta_b_A_O, T) * this._monod(S_O, K_O) * X_A;
     rates[11] = this._arrhenius(b_A_NO, theta_b_A_NO, T) * this._inv_monod(S_O, K_A_O) * this._monod(S_NO, K_NO) * X_A;
-    
+
     return rates;
   }
 
@@ -208,4 +208,4 @@ class ASM3 {
   }
 }
 
-module.exports = ASM3;
+module.exports = ASM3;

src/reaction_modules/asm3_class.js

@@ -13,7 +13,7 @@ class ASM3 {
     this.kin_params = {
       // Hydrolysis
       k_H:      3.,   // hydrolysis rate constant       [g X_S g-1 X_H d-1]
-      K_X:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H] 
+      K_X:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H]
       // Heterotrophs
       k_STO:    5.,   // storage rate constant          [g S_S g-1 X_H d-1]
       nu_NO:    0.6,  // anoxic reduction factor        [-]
@@ -30,7 +30,7 @@ class ASM3 {
       b_STO_NO: 0.1,  // anoxic respitation rate X_STO  [d-1]
       // Autotrophs
       mu_A_max: 1.0,  // maximum specific growth rate   [d-1]
-      K_A_NH:   1.,   // saturation constant S_NH3      [g NH3-N m-3] 
+      K_A_NH:   1.,   // saturation constant S_NH3      [g NH3-N m-3]
       K_A_O:    0.5,  // saturation constant S_0        [g O2 m-3]
       K_A_HCO:  0.5,  // saturation constant S_HCO      [mole HCO3 m-3]
       b_A_O:    0.15, // aerobic respiration rate       [d-1]
@@ -132,7 +132,7 @@ class ASM3 {
   /**
    * Computes the inverse Monod equation rate value for a given concentration and half-saturation constant. Used for inhibition.
    * @param {number} c - Concentration of reaction species.
-   * @param {number} K - Half-saturation constant for the reaction species. 
+   * @param {number} K - Half-saturation constant for the reaction species.
    * @returns {number} - Inverse Monod equation rate value for the given concentration and half-saturation constant.
    */
   _inv_monod(c, K) {
@@ -171,7 +171,7 @@ class ASM3 {
   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
     const rates = Array(12);
-    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 [S_O, , S_S, S_NH, , S_NO, S_HCO, , X_S, X_H, X_STO, X_A] = 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 { 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;
 
@@ -187,12 +187,12 @@ class ASM3 {
     rates[6] = this._arrhenius(b_H_NO, theta_b_H_NO, T) * this._inv_monod(S_O, K_O) * this._monod(S_NO, K_NO) * X_H;
     rates[7] = this._arrhenius(b_STO_O, theta_b_STO_O, T) * this._monod(S_O, K_O) * X_H;
     rates[8] = this._arrhenius(b_STO_NO, theta_b_STO_NO, T) * this._inv_monod(S_O, K_O) * this._monod(S_NO, K_NO) * X_STO;
-    
+
     // Autotrophs
     rates[9] = this._arrhenius(mu_A_max, theta_mu_A, T) * this._monod(S_O, K_A_O) * this._monod(S_NH, K_A_NH) * this._monod(S_HCO, K_A_HCO) * X_A;
     rates[10] = this._arrhenius(b_A_O, theta_b_A_O, T) * this._monod(S_O, K_O) * X_A;
     rates[11] = this._arrhenius(b_A_NO, theta_b_A_NO, T) * this._inv_monod(S_O, K_A_O) * this._monod(S_NO, K_NO) * X_A;
-    
+
     return rates;
   }
 
@@ -208,4 +208,4 @@ class ASM3 {
   }
 }
 
-module.exports = ASM3;
+module.exports = ASM3;

src/specificClass.js

@@ -1,419 +1,134 @@
-const ASM3 = require('./reaction_modules/asm3_class.js');
+'use strict';
-const { create, all, isArray } = require('mathjs');
-const { assertNoNaN } = require('./utils.js');
-const { childRegistrationUtils, logger, MeasurementContainer } = require('generalFunctions');
-const EventEmitter = require('events');
 
-const mathConfig = {
+const { BaseDomain, statusBadge, POSITIONS } = require('generalFunctions');
-  matrix: 'Array' // use Array as the matrix type
+const Reactor_CSTR = require('./kinetics/cstr.js');
-};
+const Reactor_PFR  = require('./kinetics/pfr.js');
 
-const math = create(all, mathConfig);
+const SPECIES_KEYS = ['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 S_O_INDEX = 0;
+// Reactor — biological reactor orchestrator (Unit-level). Wraps a CSTR or
-const NUM_SPECIES = 13;
+// PFR kinetics engine and exposes the BaseDomain surface to BaseNodeAdapter.
-const DEBUG = false;
+// The engines own the ASM3 integration; this class wires child registration
+// through ChildRouter, holds the validated config, and presents getOutput /
+// getStatusBadge.
+class Reactor extends BaseDomain {
+  static name = 'reactor';
 
-class Reactor {
+  configure() {
-  /**
+    const flat = this._flattenEngineConfig(this.config);
-   * Reactor base class.
+    this.engine = this._buildEngine(flat);
-   * @param {object} config - Configuration object containing reactor parameters.
-   */
-  constructor(config) {
-    this.config = config;
-    // EVOLV stuff
-    this.logger = new logger(this.config.general.logging.enabled, this.config.general.logging.logLevel, config.general.name);
-    this.emitter = new EventEmitter();
-    this.measurements = new MeasurementContainer();
-    this.upstreamReactor = null;
-    this.childRegistrationUtils = new childRegistrationUtils(this); // Child registration utility
 
-    this.asm = new ASM3();
+    // Re-emit upstream-reactor stateChange and engine stateChange events on
+    // the BaseDomain emitter so adapter listeners pick them up uniformly.
+    this.engine.emitter.on('stateChange', (t) => this.emitter.emit('stateChange', t));
 
-    this.volume = config.volume; // fluid volume reactor [m3]
+    // ChildRouter dispatches to engine handlers — keeps the existing
+    // _connectMeasurement / _connectReactor wiring intact, just centralised.
+    this.router.onRegister('measurement', (child) => this.engine._connectMeasurement(child));
+    this.router.onRegister('reactor',     (child) => this.engine._connectReactor(child));
 
-    this.Fs = Array(config.n_inlets).fill(0); // fluid debits per inlet [m3 d-1]
+    // Bridge engine.measurements into the BaseDomain measurements container
-    this.Cs_in = Array.from(Array(config.n_inlets), () => new Array(NUM_SPECIES).fill(0)); // composition influents
+    // so getFlattenedOutput surfaces temperature / oxygen series.
-    this.OTR = 0.0; // oxygen transfer rate [g O2 d-1 m-3]
+    this.measurements = this.engine.measurements;
-    this.temperature = 20; // temperature [C]
-
-    this.kla = config.kla; // if NaN, use externaly provided OTR [d-1]
-
-    this.currentTime = Date.now(); // milliseconds since epoch [ms]
-    this.timeStep = 1 / (24*60*60) * this.config.timeStep; // time step in seconds, converted to days.
-    this.speedUpFactor = 60; // speed up factor for simulation, 60 means 1 minute per simulated second
   }
 
-  /**
+  // Translate the nested schema config (reactor.*, initialState.*) into the
-   * Setter for influent data.
+  // flat shape the kinetics engines accept.
-   * @param {object} input - Input object (msg) containing payload with inlet index, flow rate, and concentrations.
+  _flattenEngineConfig(config) {
-  */
+    const reactor = config.reactor || {};
-  set setInfluent(input) {
+    const init = config.initialState || {};
-    let index_in = input.payload.inlet;
+    const initialState = SPECIES_KEYS.map((k) => Number(init[k] ?? 0));
-    this.Fs[index_in] = input.payload.F;
+    return {
-    this.Cs_in[index_in] = input.payload.C;
+      general: config.general,
+      functionality: config.functionality,
+      reactor_type: reactor.reactor_type ?? 'CSTR',
+      volume: Number(reactor.volume),
+      length: Number(reactor.length),
+      resolution_L: Number(reactor.resolution_L),
+      alpha: Number(reactor.alpha),
+      n_inlets: Number(reactor.n_inlets),
+      kla: Number(reactor.kla),
+      timeStep: Number(reactor.timeStep),
+      speedUpFactor: Number(reactor.speedUpFactor) || 1,
+      initialState,
+    };
   }
 
-  /**
+  _buildEngine(flat) {
-   * Setter for OTR (Oxygen Transfer Rate).
+    // The schema enum validator lowercases the configured value, so accept
-   * @param {object} input - Input object (msg) containing payload with OTR value [g O2 d-1 m-3].
+    // either case.
-   */
+    switch (String(flat.reactor_type || '').toUpperCase()) {
-  set setOTR(input) {
+      case 'CSTR': return new Reactor_CSTR(flat);
-    this.OTR = input.payload;
+      case 'PFR':  return new Reactor_PFR(flat);
-  }
-
-  /**
-   * Getter for effluent data.
-   * @returns {object} Effluent data object (msg), defaults to inlet 0.
-   */
-  get getEffluent() { // getter for Effluent, defaults to inlet 0
-    if (isArray(this.state.at(-1))) {
-      return { topic: "Fluent", payload: { inlet: 0, F: math.sum(this.Fs), C: this.state.at(-1) }, timestamp: this.currentTime };
-    }
-    return { topic: "Fluent", payload: { inlet: 0, F: math.sum(this.Fs), C: this.state }, timestamp: this.currentTime };
-  }
-
-  /**
-   * Calculate the oxygen transfer rate (OTR) based on the dissolved oxygen concentration and temperature.
-   * @param {number} S_O - Dissolved oxygen concentration [g O2 m-3].
-   * @param {number} T - Temperature in Celsius, default to 20 C.
-   * @returns {number} - Calculated OTR [g O2 d-1 m-3].
-   */
-  _calcOTR(S_O, T = 20.0) { // caculate the OTR using basic correlation, default to temperature: 20 C
-    let S_O_sat = 14.652 - 4.1022e-1 * T + 7.9910e-3 * T*T + 7.7774e-5 * T*T*T;
-    return this.kla * (S_O_sat - S_O);
-  }
-
-  /**
-   * Clip values in an array to zero.
-   * @param {Array} arr - Array of values to clip.
-   * @returns {Array} - New array with values clipped to zero.
-   */
-  _arrayClip2Zero(arr) {
-    if (Array.isArray(arr)) {
-      return arr.map(x => this._arrayClip2Zero(x));
-    } else {
-      return arr < 0 ? 0 : arr;
-    }
-  }
-
-  registerChild(child, softwareType) {
-    switch (softwareType) {
-      case "measurement":
-        this.logger.debug(`Registering measurement child.`);
-        this._connectMeasurement(child);
-        break;
-      case "reactor":
-        this.logger.debug(`Registering reactor child.`);
-        this._connectReactor(child);
-        break;
-
       default:
-        this.logger.error(`Unrecognized softwareType: ${softwareType}`);
+        this.logger.warn(`Unknown reactor type: ${flat.reactor_type}. Falling back to CSTR.`);
+        return new Reactor_CSTR(flat);
     }
   }
 
-  _connectMeasurement(measurement) {
+  // Adapter input setters — forwarded straight to the engine.
-    if (!measurement) {
+  set setInfluent(msg)    { this.engine.setInfluent = msg; }
-      this.logger.warn("Invalid measurement provided.");
+  set setOTR(msg)         { this.engine.setOTR = msg; }
-      return;
+  set setTemperature(msg) { this.engine.setTemperature = msg; }
-    }
+  set setDispersion(msg)  { if (this.engine instanceof Reactor_PFR) this.engine.setDispersion = msg; }
 
-    let position;
+  updateState(t) { this.engine.updateState(t); this.notifyOutputChanged(); }
-    if (measurement.config.functionality.distance !== 'undefined') {
-      position = measurement.config.functionality.distance;
-    } else {
-      position = measurement.config.functionality.positionVsParent;
-    }
-    const measurementType = measurement.config.asset.type;
-    const key = `${measurementType}_${position}`;
-    const eventName = `${measurementType}.measured.${position}`;
 
-    // Register event listener for measurement updates
+  // Engine pass-through — needed so the BaseNodeAdapter tick loop (and
-    measurement.measurements.emitter.on(eventName, (eventData) => {
+  // tests calling reactor.tick(dt) directly) drive the ASM integration.
-      this.logger.debug(`${position} ${measurementType} from ${eventData.childName}: ${eventData.value} ${eventData.unit}`);
+  // Without this the Node-RED tick fires `source.tick?.()`, gets undefined,
-
+  // and the kinetics state never advances.
-      // Store directly in parent's measurement container
+  tick(timeStep) {
-      this.measurements
+    const result = this.engine.tick(timeStep);
-        .type(measurementType)
+    this.notifyOutputChanged();
-        .variant("measured")
+    return result;
-        .position(position)
-        .value(eventData.value, eventData.timestamp, eventData.unit);
-      
-      this._updateMeasurement(measurementType, eventData.value, position, eventData);
-    });
   }
 
+  get getEffluent()    { return this.engine.getEffluent; }
+  get getGridProfile() { return this.engine.getGridProfile; }
+  get temperature()    { return this.engine.temperature; }
 
-  _connectReactor(reactor) {
+  // Per-tick output for Port 0 / Port 1. Carries the effluent vector plus
-    if (!reactor) {
+  // a flat per-species block keyed by SPECIES_KEYS for InfluxDB telemetry.
-      this.logger.warn("Invalid reactor provided.");
+  getOutput() {
-      return;
+    const eff = this.engine.getEffluent;
+    const C = Array.isArray(eff?.payload?.C) ? eff.payload.C : [];
+    const out = {
+      flow_total: Number(eff?.payload?.F),
+      temperature: Number(this.engine.temperature),
+    };
+    for (let i = 0; i < Math.min(SPECIES_KEYS.length, C.length); i += 1) {
+      const v = Number(C[i]);
+      if (Number.isFinite(v)) out[SPECIES_KEYS[i]] = v;
     }
-
+    return out;
-    this.upstreamReactor = reactor;
-
-    reactor.emitter.on("stateChange", (data) => {
-      this.logger.debug(`State change of upstream reactor detected.`);
-      this.updateState(data);
-    });
   }
 
-  
+  getStatusBadge() {
-  _updateMeasurement(measurementType, value, position, context) {
+    const eff = this.engine.getEffluent;
-    this.logger.debug(`---------------------- updating ${measurementType} ------------------ `);
+    const F = Number(eff?.payload?.F) || 0;
-    switch (measurementType) {
+    const SO = Array.isArray(eff?.payload?.C) ? Number(eff.payload.C[0]) : NaN;
-      case "temperature":
+    const so = Number.isFinite(SO) ? SO.toFixed(2) : '—';
-        if (position == "atEquipment") {
+    return statusBadge.compose(
-          this.temperature = value;
+      [`${this.engine.constructor.name.replace('Reactor_', '')}`,
-        }
+       `T=${Number(this.engine.temperature).toFixed(1)} C`,
-        break;
+       `F=${F.toFixed(2)} m³/d`,
-      default:
+       `S_O=${so} mg/L`],
-        this.logger.error(`Type '${measurementType}' not recognized for measured update.`);
+      { fill: 'green', shape: 'dot' },
-        return;
+    );
-    }
   }
 
-  /**
+  close() {
-   * Update the reactor state based on the new time.
+    this.engine?.emitter?.removeAllListeners?.();
-   * @param {number} newTime - New time to update reactor state to, in milliseconds since epoch.
+    super.close();
-   */
-  updateState(newTime = Date.now()) { // expect update with timestamp
-    const day2ms = 1000 * 60 * 60 * 24;
-
-    if (this.upstreamReactor) {
-      this.setInfluent = this.upstreamReactor.getEffluent;
-    }
-
-    let n_iter = Math.floor(this.speedUpFactor * (newTime-this.currentTime) / (this.timeStep*day2ms));
-    if (n_iter) {
-      let n = 0;
-      while (n < n_iter) {
-        this.tick(this.timeStep);
-        n += 1;
-      }
-      this.currentTime += n_iter * this.timeStep * day2ms / this.speedUpFactor;
-      this.emitter.emit("stateChange", this.currentTime);
-    }
   }
 }
 
-class Reactor_CSTR extends Reactor {
+module.exports = Reactor;
-  /**
+module.exports.Reactor = Reactor;
-   * Reactor_CSTR class for Continuous Stirred Tank Reactor.
+module.exports.Reactor_CSTR = Reactor_CSTR;
-   * @param {object} config - Configuration object containing reactor parameters.
+module.exports.Reactor_PFR  = Reactor_PFR;
-   */
+// POSITIONS is consumed by older test setups; surface it here so they don't
-  constructor(config) {
+// need to chase down generalFunctions internals.
-      super(config);
+module.exports.POSITIONS    = POSITIONS;
-      this.state = config.initialState;
-  }
-
-  /**
-   * Tick the reactor state using the forward Euler method.
-   * @param {number} time_step - Time step for the simulation [d].
-   * @returns {Array} - New reactor state.
-   */
-  tick(time_step) { // tick reactor state using forward Euler method
-    const inflow = math.multiply(math.divide([this.Fs], this.volume), this.Cs_in)[0];
-    const outflow = math.multiply(-1 * math.sum(this.Fs) / this.volume, this.state);
-    const reaction = this.asm.compute_dC(this.state, this.temperature);
-    const transfer = Array(NUM_SPECIES).fill(0.0);
-    transfer[S_O_INDEX] = isNaN(this.kla) ? this.OTR : this._calcOTR(this.state[S_O_INDEX], this.temperature); // calculate OTR if kla is not NaN, otherwise use externaly calculated OTR
-
-    const dC_total = math.multiply(math.add(inflow, outflow, reaction, transfer), time_step)
-    this.state = this._arrayClip2Zero(math.add(this.state, dC_total)); // clip value element-wise to avoid negative concentrations
-    if(DEBUG){
-      assertNoNaN(dC_total, "change in state");
-      assertNoNaN(this.state, "new state");
-    }
-    return this.state;
-  }
-}
-
-class Reactor_PFR extends Reactor {
-  /**
-   * Reactor_PFR class for Plug Flow Reactor.
-   * @param {object} config - Configuration object containing reactor parameters.
-   */
-  constructor(config) {
-    super(config);
-
-    this.length = config.length; // reactor length [m]
-    this.n_x = config.resolution_L; // number of slices
-
-    this.d_x = this.length / this.n_x;
-    this.A = this.volume / this.length; // crosssectional area [m2]
-
-    this.alpha = config.alpha;
-
-    this.state = Array.from(Array(this.n_x), () => config.initialState.slice())
-
-    this.D = 0.0; // axial dispersion [m2 d-1]
-
-    this.D_op = this._makeDoperator(true, true);
-    assertNoNaN(this.D_op, "Derivative operator");
-
-    this.D2_op = this._makeD2operator();
-    assertNoNaN(this.D2_op, "Second derivative operator");
-  }
-
-  /**
-   * Setter for axial dispersion.
-   * @param {object} input - Input object (msg) containing payload with dispersion value [m2 d-1].
-   */
-  set setDispersion(input) {
-    this.D = input.payload;
-  }
-
-  updateState(newTime) {
-    super.updateState(newTime);
-    let Pe_local = this.d_x*math.sum(this.Fs)/(this.D*this.A)
-    let Co_D = this.D*this.timeStep/(this.d_x*this.d_x);
-
-    (Pe_local >= 2) && this.logger.warn(`Local Péclet number (${Pe_local}) is too high! Increase reactor resolution.`);
-    (Co_D >= 0.5) && this.logger.warn(`Courant number (${Co_D}) is too high! Reduce time step size.`);
-
-    if(DEBUG) {
-      console.log("Inlet state max " + math.max(this.state[0]))
-      console.log("Pe total " + this.length*math.sum(this.Fs)/(this.D*this.A));
-      console.log("Pe local " + Pe_local);
-      console.log("Co ad " + math.sum(this.Fs)*this.timeStep/(this.A*this.d_x));
-      console.log("Co D " + Co_D);
-    }
-  }
-
-  /**
-   * Tick the reactor state using explicit finite difference method.
-   * @param {number} time_step - Time step for the simulation [d].
-   * @returns {Array} - New reactor state.
-   */
-  tick(time_step) {
-    const dispersion = math.multiply(this.D / (this.d_x*this.d_x), this.D2_op, this.state);
-    const advection = math.multiply(-1 * math.sum(this.Fs) / (this.A*this.d_x), this.D_op, this.state);
-    const reaction = this.state.map((state_slice) => this.asm.compute_dC(state_slice, this.temperature));
-    const transfer = Array.from(Array(this.n_x), () => new Array(NUM_SPECIES).fill(0));
-
-    if (isNaN(this.kla)) { // calculate OTR if kla is not NaN, otherwise use externally calculated OTR
-      for (let i = 1; i < this.n_x - 1; i++) {
-        transfer[i][S_O_INDEX] = this.OTR * this.n_x/(this.n_x-2);
-      }
-    } else {
-    for (let i = 1; i < this.n_x - 1; i++) {
-        transfer[i][S_O_INDEX] = this._calcOTR(this.state[i][S_O_INDEX], this.temperature)  * this.n_x/(this.n_x-2);
-      }
-    }
-
-    const dC_total = math.multiply(math.add(dispersion, advection, reaction, transfer), time_step);
-
-    const stateNew = math.add(this.state, dC_total);
-    this._applyBoundaryConditions(stateNew);
-
-    if (DEBUG) {
-      assertNoNaN(dispersion, "dispersion");
-      assertNoNaN(advection, "advection");
-      assertNoNaN(reaction, "reaction");
-      assertNoNaN(dC_total, "change in state");
-      assertNoNaN(stateNew, "new state post BC");
-    }
-
-    this.state = this._arrayClip2Zero(stateNew);
-    return stateNew;
-  }
-
-  _updateMeasurement(measurementType, value, position, context) {
-    switch(measurementType) {
-      case "quantity (oxygen)":
-        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
-        break;
-      default:
-        super._updateMeasurement(measurementType, value, position, context);
-    }
-  }
-
-  /**
-   * Apply boundary conditions to the reactor state.
-   *  for inlet, apply generalised Danckwerts BC, if there is not flow, apply Neumann BC with no flux
-   *  for outlet, apply regular Danckwerts BC (Neumann BC with no flux)
-   * @param {Array} state - Current reactor state without enforced BCs.
-   */
-  _applyBoundaryConditions(state) {
-    if (math.sum(this.Fs) > 0) { // Danckwerts BC
-      const BC_C_in = math.multiply(1 / math.sum(this.Fs), [this.Fs], this.Cs_in)[0];
-      const BC_dispersion_term = (1-this.alpha)*this.D*this.A/(math.sum(this.Fs)*this.d_x);
-      state[0] = math.multiply(1/(1+BC_dispersion_term), math.add(BC_C_in, math.multiply(BC_dispersion_term, state[1])));
-    } else {
-      state[0] = state[1];
-    }
-    // Neumann BC (no flux)
-    state[this.n_x-1] = state[this.n_x-2];
-  }
-
-  /**
-   * Create finite difference first derivative operator.
-   * @param {boolean} central - Use central difference scheme if true, otherwise use upwind scheme.
-   * @param {boolean} higher_order - Use higher order scheme if true, otherwise use first order scheme.
-   * @returns {Array} - First derivative operator matrix.
-   */
-  _makeDoperator(central = false, higher_order = false) { // create gradient operator
-    if (higher_order) {
-      if (central) {
-        const I = math.resize(math.diag(Array(this.n_x).fill(1/12), -2), [this.n_x, this.n_x]);
-        const A = math.resize(math.diag(Array(this.n_x).fill(-2/3), -1), [this.n_x, this.n_x]);
-        const B = math.resize(math.diag(Array(this.n_x).fill(2/3), 1), [this.n_x, this.n_x]);
-        const C = math.resize(math.diag(Array(this.n_x).fill(-1/12), 2), [this.n_x, this.n_x]);
-        const D = math.add(I, A, B, C);
-        const NearBoundary = Array(this.n_x).fill(0.0);
-        NearBoundary[0] = -1/4;
-        NearBoundary[1] = -5/6;
-        NearBoundary[2] = 3/2;
-        NearBoundary[3] = -1/2;
-        NearBoundary[4] = 1/12;
-        D[1] = NearBoundary;
-        NearBoundary.reverse();
-        D[this.n_x-2] = math.multiply(-1, NearBoundary);
-        D[0] = Array(this.n_x).fill(0); // set by BCs elsewhere
-        D[this.n_x-1] = Array(this.n_x).fill(0);
-        return D;
-      } else {
-        throw new Error("Upwind higher order method not implemented! Use central scheme instead.");
-      }
-    } else {
-      const I = math.resize(math.diag(Array(this.n_x).fill(1 / (1+central)), central), [this.n_x, this.n_x]);
-      const A = math.resize(math.diag(Array(this.n_x).fill(-1 / (1+central)), -1), [this.n_x, this.n_x]);
-      const D = math.add(I, A);
-      D[0] = Array(this.n_x).fill(0); // set by BCs elsewhere
-      D[this.n_x-1] = Array(this.n_x).fill(0);
-      return D;
-    }
-  }
-
-  /**
-   * Create central finite difference second derivative operator.
-   * @returns {Array} - Second derivative operator matrix.
-   */
-  _makeD2operator() { // create the central second derivative operator
-    const I = math.diag(Array(this.n_x).fill(-2), 0);
-    const A = math.resize(math.diag(Array(this.n_x).fill(1), 1), [this.n_x, this.n_x]);
-    const B = math.resize(math.diag(Array(this.n_x).fill(1), -1), [this.n_x, this.n_x]);
-    const D2 = math.add(I, A, B);
-    D2[0] = Array(this.n_x).fill(0); // set by BCs elsewhere
-    D2[this.n_x - 1] = Array(this.n_x).fill(0);
-    return D2;
-  }
-}
-
-module.exports = { Reactor_CSTR, Reactor_PFR };
-
-// DEBUG
-// 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
-// let initial_state = [0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1];
-// const Reactor = new Reactor_PFR(200, 10, 10, 1, 100, initial_state);
-// Reactor.Cs_in[0] = [0.0, 30., 100., 16., 0., 0., 5., 25., 75., 30., 0., 0., 125.];
-// Reactor.Fs[0] = 10;
-// Reactor.D = 0.01;
-// let N = 0;
-// while (N < 5000) {
-//     console.log(Reactor.tick(0.001));
-//     N += 1;
-// }

src/utils.js

@@ -1,18 +1,18 @@
-/**
+/**
- * Assert that no NaN values are present in an array.
+ * Assert that no NaN values are present in an array.
- * @param {Array} arr 
+ * @param {Array} arr 
- * @param {string} label 
+ * @param {string} label 
- */
+ */
-function assertNoNaN(arr, label = "array") {
+function assertNoNaN(arr, label = "array") {
-    if (Array.isArray(arr)) {
+    if (Array.isArray(arr)) {
-        for (const el of arr) {
+        for (const el of arr) {
-            assertNoNaN(el, label);
+            assertNoNaN(el, label);
-        }
+        }
-    } else {
+    } else {
-        if (Number.isNaN(arr)) {
+        if (Number.isNaN(arr)) {
-            throw new Error(`NaN detected in ${label}!`);
+            throw new Error(`NaN detected in ${label}!`);
-        }
+        }
-    }
+    }
-}
+}
-
+
 module.exports = { assertNoNaN };

test/README.md

@@ -0,0 +1,12 @@
+# reactor Test Suite Layout
+
+Required EVOLV layout:
+- basic/
+- integration/
+- edge/
+- helpers/
+
+Baseline structure tests:
+- basic/structure-module-load.basic.test.js
+- integration/structure-examples.integration.test.js
+- edge/structure-examples-node-type.edge.test.js

test/basic/constructor.basic.test.js

@@ -0,0 +1,83 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// The pre-refactor _loadConfig / _setupClass private methods are gone —
+// config build is exposed via buildDomainConfig (override hook in
+// CONTRACTS.md §2), and engine selection is observable via
+// `inst.source.engine instanceof Reactor_CSTR | Reactor_PFR` after a
+// full `new nodeClass(...)` construction.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { Reactor_CSTR, Reactor_PFR } = require('../../src/specificClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('buildDomainConfig coerces numeric fields and builds initial state vector', () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+  try {
+    const dc = inst.buildDomainConfig(
+      makeUiConfig({
+        volume: '12.5',
+        length: '9',
+        resolution_L: '7',
+        alpha: '0.5',
+        n_inlets: '3',
+        timeStep: '2',
+        S_O_init: '1.1',
+      }),
+    );
+
+    assert.equal(dc.reactor.volume, 12.5);
+    assert.equal(dc.reactor.length, 9);
+    assert.equal(dc.reactor.resolution_L, 7);
+    assert.equal(dc.reactor.alpha, 0.5);
+    assert.equal(dc.reactor.n_inlets, 3);
+    assert.equal(dc.reactor.timeStep, 2);
+    assert.equal(Object.keys(dc.initialState).length, 13);
+    assert.equal(dc.initialState.S_O, 1.1);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('Reactor wrapper instantiates CSTR engine when configured as CSTR', () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig({ reactor_type: 'CSTR' }), makeRED(), node, 'reactor');
+  try {
+    assert.ok(inst.source.engine instanceof Reactor_CSTR);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('Reactor wrapper instantiates PFR engine when configured as PFR', () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig({ reactor_type: 'PFR' }), makeRED(), node, 'reactor');
+  try {
+    assert.ok(inst.source.engine instanceof Reactor_PFR);
+  } finally {
+    closeNode(node);
+  }
+});

test/basic/cstr-tick.basic.test.js

@@ -0,0 +1,42 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+const NUM_SPECIES = 13;
+
+test('Reactor_CSTR tick clips negative concentrations to zero', () => {
+  const reactor = new Reactor_CSTR(
+    makeReactorConfig({
+      reactor_type: 'CSTR',
+      volume: 1,
+      n_inlets: 1,
+      kla: NaN,
+      S_O_init: 0.1,
+      S_I_init: 0.1,
+      S_S_init: 0.1,
+      S_NH_init: 0.1,
+      S_N2_init: 0.1,
+      S_NO_init: 0.1,
+      S_HCO_init: 0.1,
+      X_I_init: 0.1,
+      X_S_init: 0.1,
+      X_H_init: 0.1,
+      X_STO_init: 0.1,
+      X_A_init: 0.1,
+      X_TS_init: 0.1,
+    }),
+  );
+
+  reactor.asm = {
+    compute_dC: () => Array(NUM_SPECIES).fill(0),
+  };
+  reactor.Fs[0] = 1;
+  reactor.Cs_in[0] = Array(NUM_SPECIES).fill(0);
+
+  reactor.tick(1);
+
+  assert.equal(reactor.state.every((v) => Number.isFinite(v) && v >= 0), true);
+  assert.equal(reactor.state.every((v) => v === 0), true);
+});

test/basic/effluent-shape.basic.test.js

@@ -0,0 +1,38 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR, Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+test('CSTR getEffluent returns flat concentration vector', () => {
+  const reactor = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR', n_inlets: 1 }));
+  reactor.state = Array.from({ length: 13 }, (_, i) => i + 1);
+  reactor.Fs[0] = 5;
+
+  const effluent = reactor.getEffluent;
+
+  assert.equal(effluent.topic, 'Fluent');
+  assert.equal(effluent.payload.inlet, 0);
+  assert.equal(effluent.payload.F, 5);
+  assert.deepEqual(effluent.payload.C, reactor.state);
+});
+
+test('PFR getEffluent returns last slice concentration vector', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', n_inlets: 1, length: 10, resolution_L: 4 }),
+  );
+
+  reactor.state = [
+    Array(13).fill(10),
+    Array(13).fill(20),
+    Array(13).fill(30),
+    Array(13).fill(40),
+  ];
+  reactor.Fs[0] = 7;
+
+  const effluent = reactor.getEffluent;
+
+  assert.equal(effluent.topic, 'Fluent');
+  assert.equal(effluent.payload.F, 7);
+  assert.deepEqual(effluent.payload.C, Array(13).fill(40));
+});

test/basic/grid-profile.basic.test.js

@@ -0,0 +1,45 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR, Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+test('CSTR getGridProfile returns null', () => {
+  const reactor = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR' }));
+  assert.equal(reactor.getGridProfile, null);
+});
+
+test('PFR getGridProfile returns state matrix with correct dimensions', () => {
+  const n_x = 8;
+  const length = 40;
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', resolution_L: n_x, length }),
+  );
+
+  const profile = reactor.getGridProfile;
+  assert.notEqual(profile, null);
+  assert.equal(profile.n_x, n_x);
+  assert.equal(profile.d_x, length / n_x);
+  assert.equal(profile.length, length);
+  assert.equal(profile.grid.length, n_x, 'grid should have n_x rows');
+  assert.equal(profile.grid[0].length, 13, 'each row should have 13 species');
+  assert.ok(Array.isArray(profile.species), 'species list should be an array');
+  assert.equal(profile.species.length, 13);
+  assert.equal(profile.species[3], 'S_NH');
+  assert.equal(typeof profile.timestamp, 'number');
+});
+
+test('PFR getGridProfile is mutation-safe', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', resolution_L: 5, length: 10 }),
+  );
+
+  const profile = reactor.getGridProfile;
+  const originalValue = reactor.state[0][3]; // S_NH at cell 0
+
+  // Mutate the returned grid
+  profile.grid[0][3] = 999;
+
+  // Reactor internal state should be unchanged
+  assert.equal(reactor.state[0][3], originalValue, 'mutating grid copy must not affect reactor state');
+});

test/basic/input-routing.basic.test.js

@@ -0,0 +1,111 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// The pre-refactor _attachInputHandler private switch is gone — input
+// dispatch goes through the commands registry that BaseNodeAdapter builds
+// at construction. Tests fire msgs through `node.handlers.input` and
+// observe via `node.sends`, `inst.source.engine.*`, and per-fire calls
+// captured on a child stub registered through `RED.nodes.getNode(id)`.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeNode(id = 'reactor-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function makeRED(nodeMap = {}) {
+  return { nodes: { getNode: (id) => nodeMap[id] || null } };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('legacy alias topics drive engine setters and updateState', async () => {
+  const childSource = {
+    id: 'child-source-A',
+    config: { general: { id: 'child-source-A' }, functionality: { softwareType: 'measurement', positionVsParent: 'upstream' }, asset: { type: 'temperature' } },
+  };
+  const node = makeNode();
+  const RED = makeRED({ childA: { source: childSource } });
+  const inst = new nodeClass(makeUiConfig(), RED, node, 'reactor');
+
+  try {
+    let doneCount = 0;
+    const done = () => { doneCount += 1; };
+
+    // data.clock alias → updateState(timestamp). Capture currentTime
+    // before/after to verify the engine advanced.
+    const t0 = inst.source.engine.currentTime;
+    await node.handlers.input({ topic: 'clock', timestamp: t0 + 1 }, () => {}, done);
+
+    // Fluent alias → engine setInfluent setter.
+    await node.handlers.input(
+      { topic: 'Fluent', payload: { inlet: 0, F: 7, C: [1,2,3,4,5,6,7,8,9,10,11,12,13] } },
+      () => {}, done,
+    );
+    assert.equal(inst.source.engine.Fs[0], 7);
+    assert.deepEqual(inst.source.engine.Cs_in[0], [1,2,3,4,5,6,7,8,9,10,11,12,13]);
+
+    // OTR alias → engine setOTR setter.
+    await node.handlers.input({ topic: 'OTR', payload: 3.5 }, () => {}, done);
+    assert.equal(inst.source.engine.OTR, 3.5);
+
+    // Temperature alias → engine setTemperature setter.
+    await node.handlers.input({ topic: 'Temperature', payload: 18.2 }, () => {}, done);
+    assert.equal(inst.source.engine.temperature, 18.2);
+
+    // Dispersion alias — CSTR engine does not own a setDispersion setter
+    // (only PFR does); the Reactor wrapper guards on engine type and the
+    // dispatch should silently return without throwing.
+    await node.handlers.input({ topic: 'Dispersion', payload: 0.2 }, () => {}, done);
+
+    // registerChild alias → registers via childRegistrationUtils.
+    // The handler resolves the child via RED.nodes.getNode(payload).source.
+    await node.handlers.input(
+      { topic: 'registerChild', payload: 'childA', positionVsParent: 'upstream' },
+      () => {}, done,
+    );
+
+    assert.equal(doneCount, 6);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('canonical topics are accepted (data.fluent, data.otr, data.temperature)', async () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let done = 0;
+    await node.handlers.input(
+      { topic: 'data.fluent', payload: { inlet: 0, F: 11, C: [0,0,0,0,0,0,0,0,0,0,0,0,0] } },
+      () => {}, () => { done += 1; },
+    );
+    assert.equal(inst.source.engine.Fs[0], 11);
+
+    await node.handlers.input({ topic: 'data.otr', payload: 4.2 }, () => {}, () => { done += 1; });
+    assert.equal(inst.source.engine.OTR, 4.2);
+
+    await node.handlers.input({ topic: 'data.temperature', payload: 19.7 }, () => {}, () => { done += 1; });
+    assert.equal(inst.source.engine.temperature, 19.7);
+
+    assert.equal(done, 3);
+  } finally {
+    closeNode(node);
+  }
+});

test/basic/pfr-operators.basic.test.js

@@ -0,0 +1,27 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+test('Reactor_PFR derivative operators have expected dimensions and boundary rows', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({
+      reactor_type: 'PFR',
+      length: 12,
+      resolution_L: 6,
+      volume: 60,
+      n_inlets: 1,
+    }),
+  );
+
+  assert.equal(reactor.D_op.length, reactor.n_x);
+  assert.equal(reactor.D2_op.length, reactor.n_x);
+  assert.equal(reactor.D_op.every((row) => row.length === reactor.n_x), true);
+  assert.equal(reactor.D2_op.every((row) => row.length === reactor.n_x), true);
+
+  assert.deepEqual(reactor.D_op[0], Array(reactor.n_x).fill(0));
+  assert.deepEqual(reactor.D_op[reactor.n_x - 1], Array(reactor.n_x).fill(0));
+  assert.deepEqual(reactor.D2_op[0], Array(reactor.n_x).fill(0));
+  assert.deepEqual(reactor.D2_op[reactor.n_x - 1], Array(reactor.n_x).fill(0));
+});

test/basic/register-child.basic.test.js

@@ -0,0 +1,84 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// The pre-refactor _registerChild method was renamed to
+// _scheduleRegistration inside BaseNodeAdapter and now fires automatically
+// 100ms after construction. We verify the emission by capturing the Port-2
+// message on `node.sends` after the registration delay elapses.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('scheduled child.register message lands on Port 2 after construction', async () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ positionVsParent: 'downstream' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+
+  try {
+    // BaseNodeAdapter._scheduleRegistration uses a 100ms setTimeout; wait
+    // slightly longer to let it fire.
+    await new Promise((r) => setTimeout(r, 130));
+
+    // The registration send is the [null, null, {child.register}] triple.
+    const regSends = node.sends.filter(
+      (s) => Array.isArray(s) && s[0] === null && s[1] === null && s[2] && s[2].topic === 'child.register',
+    );
+    assert.equal(regSends.length, 1, 'exactly one child.register message expected');
+    const msg = regSends[0][2];
+    assert.equal(msg.topic, 'child.register');
+    assert.equal(msg.payload, node.id);
+    assert.equal(msg.positionVsParent, 'downstream');
+    // After construction the source is exposed on the node for sibling lookup.
+    assert.strictEqual(node.source, inst.source);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('child.register handler ignores unknown child ids without throwing', async () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let done = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { topic: 'child.register', payload: 'missing-child', positionVsParent: 'upstream' },
+        () => {},
+        () => { done += 1; },
+      );
+    });
+    assert.equal(done, 1);
+    // No child should have been registered into the engine's registry.
+    const registered = inst.source.engine.childRegistrationUtils;
+    assert.ok(registered, 'childRegistrationUtils exists on engine');
+  } finally {
+    closeNode(node);
+  }
+});

test/basic/speedup-factor.basic.test.js

@@ -0,0 +1,91 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface for
+// the nodeClass-level checks, and the public Reactor_CSTR engine surface
+// for the domain-level checks. The pre-refactor private nodeClass methods
+// are gone — `buildDomainConfig` is the documented override hook
+// (CONTRACTS.md §2) and is fair game to call on a real constructed
+// instance.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { Reactor_CSTR } = require('../../src/specificClass');
+const { makeReactorConfig, makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('Reactor_CSTR engine defaults speedUpFactor to 1 when not in config', () => {
+  const config = makeReactorConfig();
+  const reactor = new Reactor_CSTR(config);
+  assert.equal(reactor.speedUpFactor, 1, 'speedUpFactor should default to 1');
+});
+
+test('Reactor_CSTR engine accepts speedUpFactor from config', () => {
+  const config = makeReactorConfig();
+  config.speedUpFactor = 10;
+  const reactor = new Reactor_CSTR(config);
+  assert.equal(reactor.speedUpFactor, 10, 'speedUpFactor should be read from config');
+});
+
+test('Reactor_CSTR engine accepts speedUpFactor = 60 for accelerated simulation', () => {
+  const config = makeReactorConfig();
+  config.speedUpFactor = 60;
+  const reactor = new Reactor_CSTR(config);
+  assert.equal(reactor.speedUpFactor, 60, 'speedUpFactor=60 should be accepted');
+});
+
+test('buildDomainConfig propagates speedUpFactor from uiConfig', () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+  try {
+    const dc = inst.buildDomainConfig(makeUiConfig({ speedUpFactor: 5 }));
+    assert.equal(dc.reactor.speedUpFactor, 5);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('buildDomainConfig defaults speedUpFactor to 1 when missing from uiConfig', () => {
+  const node = makeNode();
+  const inst = new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+  try {
+    const ui = makeUiConfig();
+    delete ui.speedUpFactor;
+    const dc = inst.buildDomainConfig(ui);
+    assert.equal(dc.reactor.speedUpFactor, 1);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('updateState with speedUpFactor=1 advances roughly real-time', () => {
+  const config = makeReactorConfig();
+  config.speedUpFactor = 1;
+  config.n_inlets = 1;
+  const reactor = new Reactor_CSTR(config);
+
+  const t0 = reactor.currentTime;
+  reactor.updateState(t0 + 2000);
+
+  const elapsed = reactor.currentTime - t0;
+  assert.ok(elapsed < 5000, `Elapsed ${elapsed}ms should be close to 2000ms, not 120000ms (old 60x factor)`);
+});

test/basic/structure-module-load.basic.test.js

@@ -0,0 +1,8 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+test('reactor module load smoke', () => {
+  assert.doesNotThrow(() => {
+    require('../../reactor.js');
+  });
+});

test/edge/invalid-reactor-type.edge.test.js

@@ -0,0 +1,65 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// The schema validator coerces `reactor_type` through the enum — values
+// outside `CSTR` / `PFR` are remapped to the default `CSTR` at validation
+// time. The Reactor wrapper additionally falls back to CSTR if anything
+// unrecognised slips through (defensive guard). Either way, the observable
+// effect after `new nodeClass(...)` is `inst.source.engine instanceof
+// Reactor_CSTR`.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { Reactor_CSTR } = require('../../src/specificClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('Reactor wrapper falls back to CSTR when reactor_type is unknown', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: 'UNKNOWN_TYPE' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+  try {
+    assert.ok(inst.source.engine instanceof Reactor_CSTR);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('Reactor wrapper falls back to CSTR when reactor_type is empty string', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: '' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+  try {
+    assert.ok(inst.source.engine instanceof Reactor_CSTR);
+  } finally {
+    closeNode(node);
+  }
+});

test/edge/invalid-topic.edge.test.js

@@ -0,0 +1,68 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface. The
+// commands registry built by BaseNodeAdapter logs a warn on unknown topics
+// and still calls done — no throw.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('unknown input topic does not throw and still calls done', async () => {
+  const node = makeNode();
+  new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let doneCalled = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { topic: 'somethingUnknown', payload: 1 },
+        () => {},
+        () => { doneCalled += 1; },
+      );
+    });
+    assert.equal(doneCalled, 1);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('missing topic field is handled gracefully', async () => {
+  const node = makeNode();
+  new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let doneCalled = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { payload: 'no-topic-here' },
+        () => {},
+        () => { doneCalled += 1; },
+      );
+    });
+    assert.equal(doneCalled, 1);
+  } finally {
+    closeNode(node);
+  }
+});

test/edge/missing-child.edge.test.js

@@ -0,0 +1,91 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// A child.register / registerChild msg with an unknown id should resolve
+// to no-op (the handler logs warn, no throw) and still call done.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED(nodeMap = {}) {
+  return { nodes: { getNode: (id) => nodeMap[id] || null } };
+}
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+test('registerChild alias with unknown id is ignored without throwing', async () => {
+  const node = makeNode();
+  new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let done = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { topic: 'registerChild', payload: 'missing-child', positionVsParent: 'upstream' },
+        () => {},
+        () => { done += 1; },
+      );
+    });
+    assert.equal(done, 1);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('child.register canonical topic with unknown id is ignored without throwing', async () => {
+  const node = makeNode();
+  new nodeClass(makeUiConfig(), makeRED(), node, 'reactor');
+
+  try {
+    let done = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { topic: 'child.register', payload: 'missing-child', positionVsParent: 'upstream' },
+        () => {},
+        () => { done += 1; },
+      );
+    });
+    assert.equal(done, 1);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('child.register with a child that has no .source is ignored without throwing', async () => {
+  const node = makeNode();
+  // The looked-up RED node exists but lacks a `.source` — the handler
+  // guards against this and logs warn.
+  new nodeClass(makeUiConfig(), makeRED({ orphan: {} }), node, 'reactor');
+
+  try {
+    let done = 0;
+    await assert.doesNotReject(async () => {
+      await node.handlers.input(
+        { topic: 'child.register', payload: 'orphan', positionVsParent: 'upstream' },
+        () => {},
+        () => { done += 1; },
+      );
+    });
+    assert.equal(done, 1);
+  } finally {
+    closeNode(node);
+  }
+});

test/edge/pfr-measurement-grid.edge.test.js

@@ -0,0 +1,16 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+test('oxygen measurement at exact reactor length is clamped to the last PFR grid index', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', length: 10, resolution_L: 5, n_inlets: 1 }),
+  );
+
+  assert.doesNotThrow(() => {
+    reactor._updateMeasurement('quantity (oxygen)', 2.5, 10, {});
+  });
+  assert.equal(reactor.state[reactor.n_x - 1][0], 2.5);
+});

test/edge/structure-examples-node-type.edge.test.js

@@ -0,0 +1,11 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+const fs = require('node:fs');
+const path = require('node:path');
+
+const flow = JSON.parse(fs.readFileSync(path.resolve(__dirname, '../../examples/basic.flow.json'), 'utf8'));
+
+test('basic example includes node type reactor', () => {
+  const count = flow.filter((n) => n && n.type === 'reactor').length;
+  assert.equal(count >= 1, true);
+});

test/edge/zero-dispersion.edge.test.js

@@ -0,0 +1,27 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+const DAY_MS = 1000 * 60 * 60 * 24;
+
+test('updateState warns when local Peclet number is too high at zero dispersion', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', length: 10, resolution_L: 5, volume: 50, n_inlets: 1 }),
+  );
+
+  const warnings = [];
+  reactor.logger.warn = (msg) => warnings.push(String(msg));
+
+  reactor.currentTime = 0;
+  reactor.timeStep = 1;
+  reactor.speedUpFactor = 1;
+  reactor.Fs[0] = 2;
+  reactor.D = 0;
+  reactor.tick = () => reactor.state;
+
+  reactor.updateState(DAY_MS);
+
+  assert.equal(warnings.some((w) => w.includes('Péclet number') || w.includes('Peclet number')), true);
+});

test/helpers/factories.js

@@ -0,0 +1,149 @@
+const EventEmitter = require('node:events');
+
+function makeUiConfig(overrides = {}) {
+  return {
+    name: 'reactor-test',
+    reactor_type: 'CSTR',
+    volume: 100,
+    length: 10,
+    resolution_L: 5,
+    alpha: 0,
+    n_inlets: 1,
+    kla: NaN,
+    S_O_init: 0,
+    S_I_init: 30,
+    S_S_init: 100,
+    S_NH_init: 16,
+    S_N2_init: 0,
+    S_NO_init: 0,
+    S_HCO_init: 5,
+    X_I_init: 25,
+    X_S_init: 75,
+    X_H_init: 30,
+    X_STO_init: 0,
+    X_A_init: 0.001,
+    X_TS_init: 125,
+    timeStep: 1,
+    enableLog: false,
+    logLevel: 'error',
+    positionVsParent: 'atEquipment',
+    ...overrides,
+  };
+}
+
+function makeReactorConfig(overrides = {}) {
+  const ui = makeUiConfig(overrides);
+  return {
+    general: {
+      id: 'reactor-node-1',
+      name: ui.name,
+      unit: null,
+      logging: {
+        enabled: ui.enableLog,
+        logLevel: ui.logLevel,
+      },
+    },
+    functionality: {
+      positionVsParent: ui.positionVsParent || 'atEquipment',
+      softwareType: 'reactor',
+    },
+    reactor_type: ui.reactor_type,
+    volume: Number(ui.volume),
+    length: Number(ui.length),
+    resolution_L: Number(ui.resolution_L),
+    alpha: Number(ui.alpha),
+    n_inlets: Number(ui.n_inlets),
+    kla: Number(ui.kla),
+    initialState: [
+      Number(ui.S_O_init),
+      Number(ui.S_I_init),
+      Number(ui.S_S_init),
+      Number(ui.S_NH_init),
+      Number(ui.S_N2_init),
+      Number(ui.S_NO_init),
+      Number(ui.S_HCO_init),
+      Number(ui.X_I_init),
+      Number(ui.X_S_init),
+      Number(ui.X_H_init),
+      Number(ui.X_STO_init),
+      Number(ui.X_A_init),
+      Number(ui.X_TS_init),
+    ],
+    timeStep: Number(ui.timeStep),
+  };
+}
+
+function makeNodeStub() {
+  const handlers = {};
+  const sent = [];
+  const warns = [];
+  const errors = [];
+  const statuses = [];
+
+  return {
+    id: 'reactor-node-1',
+    source: null,
+    on(event, cb) {
+      handlers[event] = cb;
+    },
+    send(msg) {
+      sent.push(msg);
+    },
+    warn(msg) {
+      warns.push(msg);
+    },
+    error(msg) {
+      errors.push(msg);
+    },
+    status(msg) {
+      statuses.push(msg);
+    },
+    _handlers: handlers,
+    _sent: sent,
+    _warns: warns,
+    _errors: errors,
+    _statuses: statuses,
+  };
+}
+
+function makeREDStub(nodeMap = {}) {
+  return {
+    nodes: {
+      getNode(id) {
+        return nodeMap[id] || null;
+      },
+      createNode() {},
+      registerType() {},
+    },
+    httpAdmin: {
+      get() {},
+    },
+  };
+}
+
+function makeMeasurementChild({
+  id = 'measurement-1',
+  name = 'temp-sensor-1',
+  distance = 'atEquipment',
+  positionVsParent = 'atEquipment',
+  type = 'temperature',
+} = {}) {
+  return {
+    config: {
+      general: { id, name },
+      functionality: { distance, positionVsParent, softwareType: 'measurement' },
+      asset: { type },
+    },
+    measurements: {
+      emitter: new EventEmitter(),
+    },
+  };
+}
+
+module.exports = {
+  makeUiConfig,
+  makeReactorConfig,
+  makeNodeStub,
+  makeREDStub,
+  makeMeasurementChild,
+};

test/integration/measurement-temperature.integration.test.js

@@ -0,0 +1,26 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR } = require('../../src/specificClass');
+const { makeReactorConfig, makeMeasurementChild } = require('../helpers/factories');
+
+test('measurement child temperature event updates reactor temperature', () => {
+  const reactor = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR' }));
+
+  const measurement = makeMeasurementChild({
+    type: 'temperature',
+    distance: 'atEquipment',
+    positionVsParent: 'upstream',
+  });
+
+  reactor.registerChild(measurement, 'measurement');
+
+  measurement.measurements.emitter.emit('temperature.measured.atEquipment', {
+    childName: 'T-1',
+    value: 27.5,
+    unit: 'C',
+    timestamp: Date.now(),
+  });
+
+  assert.equal(reactor.temperature, 27.5);
+});

test/integration/otr-kla.integration.test.js

@@ -0,0 +1,91 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR, Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+const NUM_SPECIES = 13;
+
+test('CSTR uses external OTR when kla is NaN', () => {
+  const reactor = new Reactor_CSTR(
+    makeReactorConfig({ reactor_type: 'CSTR', kla: NaN, n_inlets: 1 }),
+  );
+
+  reactor.asm = {
+    compute_dC: () => Array(NUM_SPECIES).fill(0),
+  };
+  reactor.Fs[0] = 0;
+  reactor.OTR = 4;
+  reactor.state = Array(NUM_SPECIES).fill(0);
+
+  reactor.tick(1);
+
+  assert.equal(reactor.state[0], 4);
+});
+
+test('CSTR uses kla-based oxygen transfer when kla is finite', () => {
+  const reactor = new Reactor_CSTR(
+    makeReactorConfig({ reactor_type: 'CSTR', kla: 2, n_inlets: 1 }),
+  );
+
+  reactor.asm = {
+    compute_dC: () => Array(NUM_SPECIES).fill(0),
+  };
+  reactor.Fs[0] = 0;
+  reactor.OTR = 1;
+  reactor.state = Array(NUM_SPECIES).fill(0);
+
+  const expected = Math.min(
+    reactor._calcOTR(0, reactor.temperature),
+    reactor._calcOxygenSaturation(reactor.temperature),
+  );
+  reactor.tick(1);
+
+  assert.ok(Math.abs(reactor.state[0] - expected) < 1e-9);
+});
+
+test('PFR uses external OTR branch when kla is NaN', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', kla: NaN, n_inlets: 1, length: 8, resolution_L: 6, volume: 40 }),
+  );
+
+  reactor.asm = {
+    compute_dC: () => Array(NUM_SPECIES).fill(0),
+  };
+  reactor.Fs[0] = 0;
+  reactor.D = 0;
+  reactor.OTR = 3;
+  reactor.state = Array.from({ length: reactor.n_x }, () => Array(NUM_SPECIES).fill(0));
+
+  reactor.tick(1);
+
+  assert.equal(reactor.state[1][0], 4.5);
+  assert.equal(reactor.state[2][0], 4.5);
+  assert.equal(reactor.state[3][0], 4.5);
+  assert.equal(reactor.state[4][0], 4.5);
+});
+
+test('PFR uses kla-based transfer branch when kla is finite', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', kla: 1, n_inlets: 1, length: 8, resolution_L: 6, volume: 40 }),
+  );
+
+  reactor.asm = {
+    compute_dC: () => Array(NUM_SPECIES).fill(0),
+  };
+  reactor.Fs[0] = 0;
+  reactor.D = 0;
+  reactor.OTR = 0;
+  reactor.state = Array.from({ length: reactor.n_x }, () => Array(NUM_SPECIES).fill(0));
+
+  const expected = Math.min(
+    reactor._calcOTR(0, reactor.temperature) * (reactor.n_x / (reactor.n_x - 2)),
+    reactor._calcOxygenSaturation(reactor.temperature),
+  );
+  reactor.tick(1);
+
+  assert.ok(Math.abs(reactor.state[1][0] - expected) < 1e-9);
+  assert.ok(Math.abs(reactor.state[2][0] - expected) < 1e-9);
+  assert.ok(Math.abs(reactor.state[3][0] - expected) < 1e-9);
+  assert.ok(Math.abs(reactor.state[4][0] - expected) < 1e-9);
+});

test/integration/pfr-boundary.integration.test.js

@@ -0,0 +1,35 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_PFR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+test('_applyBoundaryConditions enforces Danckwerts inlet and Neumann outlet for flowing case', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', n_inlets: 1, length: 10, resolution_L: 5, volume: 50, alpha: 0.2 }),
+  );
+
+  reactor.Fs[0] = 2;
+  reactor.Cs_in[0] = Array(13).fill(9);
+  reactor.D = 1;
+
+  const state = Array.from({ length: reactor.n_x }, (_, i) => Array(13).fill(i));
+  reactor._applyBoundaryConditions(state);
+
+  assert.deepEqual(state[reactor.n_x - 1], state[reactor.n_x - 2]);
+  assert.equal(state[0].every((v) => Number.isFinite(v)), true);
+});
+
+test('_applyBoundaryConditions copies first interior slice when no flow is present', () => {
+  const reactor = new Reactor_PFR(
+    makeReactorConfig({ reactor_type: 'PFR', n_inlets: 1, length: 10, resolution_L: 5, volume: 50 }),
+  );
+
+  reactor.Fs[0] = 0;
+  const state = Array.from({ length: reactor.n_x }, (_, i) => Array(13).fill(i + 10));
+
+  reactor._applyBoundaryConditions(state);
+
+  assert.deepEqual(state[0], state[1]);
+  assert.deepEqual(state[reactor.n_x - 1], state[reactor.n_x - 2]);
+});

test/integration/structure-examples.integration.test.js

@@ -0,0 +1,23 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+const fs = require('node:fs');
+const path = require('node:path');
+
+const dir = path.resolve(__dirname, '../../examples');
+
+function loadJson(file) {
+  return JSON.parse(fs.readFileSync(path.join(dir, file), 'utf8'));
+}
+
+test('examples package exists for reactor', () => {
+  for (const file of ['README.md', 'basic.flow.json', 'integration.flow.json', 'edge.flow.json']) {
+    assert.equal(fs.existsSync(path.join(dir, file)), true, file + ' missing');
+  }
+});
+
+test('example flows are parseable arrays for reactor', () => {
+  for (const file of ['basic.flow.json', 'integration.flow.json', 'edge.flow.json']) {
+    const parsed = loadJson(file);
+    assert.equal(Array.isArray(parsed), true);
+  }
+});

test/integration/tick-loop.integration.test.js

@@ -0,0 +1,156 @@
+'use strict';
+
+// Phase 10 rewrite: drives only the public BaseNodeAdapter surface.
+// The pre-refactor _tick / _startTickLoop methods are gone — periodic
+// emission lives in `_emitOutputs()` (overridden in the reactor nodeClass
+// to preserve the Fluent / GridProfile Port-0 contract; delta-compressed
+// payloads can't carry the C-vector). The override is part of the
+// documented BaseNodeAdapter override surface, so we exercise it
+// directly. The fully-constructed adapter wires `inst.source.engine`,
+// `inst._output`, etc. so we don't have to assemble stub bags.
+
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const nodeClass = require('../../src/nodeClass');
+const { makeUiConfig } = require('../helpers/factories');
+
+function makeRED() { return { nodes: { getNode: () => null } }; }
+
+function makeNode(id = 'reactor-node-1') {
+  const sends = [];
+  const statuses = [];
+  const handlers = {};
+  return {
+    id, sends, statuses, handlers,
+    send(arr)  { sends.push(arr); },
+    status(b)  { statuses.push(b); },
+    on(ev, fn) { handlers[ev] = fn; },
+    warn() {}, error() {},
+  };
+}
+
+function closeNode(node) {
+  if (node.handlers.close) node.handlers.close(() => {});
+}
+
+function pickEffluentSends(node) {
+  return node.sends.filter((s) => Array.isArray(s) && s[0] && s[0].topic === 'Fluent');
+}
+
+function pickGridSends(node) {
+  return node.sends.filter((s) => Array.isArray(s) && s[0] && s[0].topic === 'GridProfile');
+}
+
+test('_emitOutputs sends the effluent message on process output (CSTR)', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: 'CSTR' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+
+  try {
+    // Reset sends so any construction-time emissions don't pollute the
+    // assertion (the registration triple lands on the same buffer).
+    node.sends.length = 0;
+    inst._emitOutputs();
+
+    const fluentSends = pickEffluentSends(node);
+    assert.equal(fluentSends.length, 1, 'exactly one Fluent message');
+    const triple = fluentSends[0];
+    assert.equal(triple[0].topic, 'Fluent');
+    assert.ok(triple[0].payload && Array.isArray(triple[0].payload.C));
+    // CSTR has no grid profile.
+    assert.equal(pickGridSends(node).length, 0);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('_emitOutputs emits a GridProfile message when engine exposes one (PFR)', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: 'PFR' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+
+  try {
+    node.sends.length = 0;
+    inst._emitOutputs();
+
+    assert.equal(pickGridSends(node).length, 1, 'exactly one GridProfile message');
+    assert.equal(pickEffluentSends(node).length, 1, 'exactly one Fluent message');
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('_emitOutputs formats per-species influx telemetry via outputUtils', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: 'CSTR' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+
+  try {
+    // Stub updateState so the engine integration does not overwrite the
+    // engineered state we want the telemetry formatter to see.
+    inst.source.updateState = () => {};
+    inst.source.engine.setInfluent = {
+      payload: { inlet: 0, F: 42, C: [2.1, 30, 100, 16, 0, 1, 8, 25, 75, 1500, 0, 15, 2500] },
+    };
+    inst.source.engine.state = [2.1, 30, 100, 16, 0, 1, 8, 25, 75, 1500, 0, 15, 2500];
+    inst.source.engine.temperature = 19.5;
+
+    let captured = null;
+    const realFormat = inst._output.formatMsg.bind(inst._output);
+    inst._output.formatMsg = (output, cfg, format) => {
+      if (format === 'influxdb') captured = { output, format };
+      return realFormat(output, cfg, format);
+    };
+
+    node.sends.length = 0;
+    inst._emitOutputs();
+
+    assert.ok(captured, 'formatMsg was called with influxdb format');
+    assert.equal(captured.format, 'influxdb');
+    assert.equal(captured.output.flow_total, 42);
+    assert.equal(captured.output.temperature, 19.5);
+    assert.equal(captured.output.S_O, 2.1);
+    assert.equal(captured.output.S_NH, 16);
+    assert.equal(captured.output.X_TS, 2500);
+  } finally {
+    closeNode(node);
+  }
+});
+
+test('Reactor.tick(dt) drives the kinetics engine and advances state', () => {
+  const node = makeNode();
+  const inst = new nodeClass(
+    makeUiConfig({ reactor_type: 'CSTR' }),
+    makeRED(),
+    node,
+    'reactor',
+  );
+
+  try {
+    // Feed an influent so the integrator has something to chew on.
+    inst.source.engine.setInfluent = {
+      payload: { inlet: 0, F: 5, C: [0,30,100,16,0,0,5,25,75,30,0,0.001,125] },
+    };
+
+    const stateBefore = JSON.stringify(inst.source.engine.state);
+    inst.source.tick(0.001);
+    const stateAfter = JSON.stringify(inst.source.engine.state);
+
+    assert.notEqual(stateBefore, stateAfter, 'engine state should advance after tick(dt)');
+  } finally {
+    closeNode(node);
+  }
+});

test/integration/upstream-reactor.integration.test.js

@@ -0,0 +1,48 @@
+const test = require('node:test');
+const assert = require('node:assert/strict');
+
+const { Reactor_CSTR } = require('../../src/specificClass');
+const { makeReactorConfig } = require('../helpers/factories');
+
+const DAY_MS = 1000 * 60 * 60 * 24;
+
+test('registering upstream reactor subscribes to upstream stateChange events', () => {
+  const downstream = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR' }));
+  const upstream = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR' }));
+
+  let calledWith = null;
+  downstream.updateState = (timestamp) => {
+    calledWith = timestamp;
+  };
+
+  downstream.registerChild(upstream, 'reactor');
+  upstream.emitter.emit('stateChange', 12345);
+
+  assert.equal(downstream.upstreamReactor, upstream);
+  assert.equal(calledWith, 12345);
+});
+
+test('updateState pulls influent from upstream reactor effluent when linked', () => {
+  const downstream = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR', n_inlets: 1, timeStep: 1 }));
+  const upstream = new Reactor_CSTR(makeReactorConfig({ reactor_type: 'CSTR', n_inlets: 1 }));
+
+  upstream.Fs[0] = 3;
+  upstream.state = Array(13).fill(11);
+
+  downstream.upstreamReactor = upstream;
+  downstream.currentTime = 0;
+  downstream.timeStep = 1;
+  downstream.speedUpFactor = 1;
+
+  let ticks = 0;
+  downstream.tick = () => {
+    ticks += 1;
+    return downstream.state;
+  };
+
+  downstream.updateState(DAY_MS);
+
+  assert.equal(ticks, 1);
+  assert.equal(downstream.Fs[0], 3);
+  assert.deepEqual(downstream.Cs_in[0], Array(13).fill(11));
+});

test/specificClass.test.js

@@ -0,0 +1,346 @@
+/**
+ * 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);
+    });
+  });
+});

wiki/Home.md

@@ -0,0 +1,285 @@
+# reactor
+
+> **Reflects code as of `b8247fc` · regenerated `2026-05-11` via `npm run wiki:all`**
+> If this banner is stale, the page may be out of date. Treat as informative, not authoritative.
+
+## 1. What this node is
+
+**reactor** is an S88 Unit that wraps an ASM3 biological-process engine — either a CSTR (fully mixed tank) or a PFR (plug-flow with axial dispersion). It integrates 13 species (S_O, S_NH, X_H, X_TS, …) and emits the effluent vector each tick. Drives a settler downstream and accepts a recirculation pump child.
+
+## 2. Position in the platform
+
+```mermaid
+flowchart LR
+    upstream[reactor<br/>upstream<br/>Unit]:::unit
+    reactor[reactor<br/>Unit]:::unit
+    settler[settler<br/>downstream<br/>Unit]:::unit
+    pump[rotatingMachine<br/>downstream<br/>Equipment]:::equip
+    tsens[measurement<br/>temperature<br/>atequipment]:::ctrl
+    osens[measurement<br/>oxygen<br/>position]:::ctrl
+
+    upstream -.stateChange.-> reactor
+    reactor -->|Fluent inlet=0| settler
+    pump -->|child.register downstream| reactor
+    tsens -->|temperature.measured.atequipment| reactor
+    osens -->|quantity (oxygen).measured.&lt;position&gt;| reactor
+    classDef unit fill:#50a8d9,color:#000
+    classDef equip fill:#86bbdd,color:#000
+    classDef ctrl fill:#a9daee,color:#000
+```
+
+S88 colours: Unit `#50a8d9`, Equipment `#86bbdd`, Control Module `#a9daee`. Source of truth: `.claude/rules/node-red-flow-layout.md`.
+
+## 3. Capability matrix
+
+| Capability | Status | Notes |
+|---|---|---|
+| ASM3 13-species ODE integration | ✅ | CSTR + PFR engines under `kinetics/`. |
+| CSTR (fully mixed) | ✅ | Single concentration vector per tick. |
+| PFR (axial discretization) | ✅ | `resolution_L` grid cells; emits `GridProfile` alongside `Fluent`. |
+| Multi-inlet mixing | ✅ | `n_inlets`; each inlet receives its own `data.fluent` with `inlet` index. |
+| Temperature reconcile from measurement | ✅ | `temperature.measured.atEquipment` writes `engine.temperature`. |
+| Oxygen reconcile (PFR) | ✅ | `quantity (oxygen).measured.<distance>` maps to nearest grid cell. |
+| KLa-driven aeration | ✅ | `reactor.kla` > 0 enables internal mass transfer; falls back to `data.otr`. |
+| Speed-up factor (sim time) | ✅ | `reactor.speedUpFactor` accelerates wall-clock → process time. |
+| Dispersion override (PFR) | ✅ | `data.dispersion` updates axial `D`. |
+| Hot-swap engine type | ❌ | `reactor_type` is read once in `configure()`. |
+
+## 4. Code map
+
+```mermaid
+flowchart TB
+    subgraph nodeRED["nodeClass.js — adapter (BaseNodeAdapter)"]
+        nc["buildDomainConfig()<br/>static DomainClass = Reactor<br/>static commands"]
+    end
+    subgraph domain["specificClass.js — orchestrator (BaseDomain)"]
+        sc["Reactor.configure()<br/>flatten config → build engine<br/>ChildRouter rules"]
+    end
+    subgraph kinetics["src/kinetics/"]
+        be["baseEngine.js<br/>shared ASM3 rate vector"]
+        cstr["cstr.js<br/>0-D integrator"]
+        pfr["pfr.js<br/>spatial discretization + dispersion"]
+    end
+    subgraph commands["src/commands/"]
+        cmds["index.js + handlers.js<br/>6 input topics"]
+    end
+    sc --> be
+    sc --> cstr
+    sc --> pfr
+    nc --> sc
+    nc --> cmds
+```
+
+| Module | Owns | Read first if you're changing… |
+|---|---|---|
+| `kinetics/baseEngine.js` | ASM3 stoichiometry + rate vector + species list. | Stoichiometric matrix, kinetic constants. |
+| `kinetics/cstr.js` | 0-D CSTR integrator + `_connectMeasurement` + `_connectReactor`. | Mixed-tank behaviour, child wiring. |
+| `kinetics/pfr.js` | Axial discretization, dispersion, grid profile emission. | PFR-specific behaviour, grid math. |
+| `commands/` | 6 input descriptors + handlers (clock, fluent, OTR, temperature, dispersion, child). | Inbound topic API, alias deprecation. |
+| `reaction_modules/` | Optional plug-in reaction modules (legacy — not yet refactored). | Adding new bio-process modules. |
+| `additional_nodes/` | Sibling Node-RED nodes (`recirculation-pump`, `settling-basin`) shipped from this repo. | Cross-node deploy in same package. |
+
+## 5. Topic contract
+
+> **Auto-generated** from `src/commands/index.js`. Do NOT hand-edit between the markers. Re-run `npm run wiki:contract`.
+
+<!-- BEGIN AUTOGEN: topic-contract -->
+
+| Canonical topic | Aliases | Payload | Unit | Effect |
+|---|---|---|---|---|
+| `data.clock` | `clock` | `any` | — | Push the simulation clock tick (timestamp / dt) to the ASM solver. |
+| `data.fluent` | `Fluent` | `object` | — | Push the influent stream (payload: {F: flow m3/h, C: [concentrations mg/L]}). |
+| `data.otr` | `OTR` | `any` | — | Push the current oxygen-transfer rate into the reactor. |
+| `data.temperature` | `Temperature` | `any` | — | Push the current reactor temperature. |
+| `data.dispersion` | `Dispersion` | `any` | — | Push a dispersion/mixing parameter update. |
+| `child.register` | `registerChild` | `any` | — | Register a child node (settler / measurement) with this reactor. |
+
+<!-- END AUTOGEN: topic-contract -->
+
+## 6. Child registration
+
+```mermaid
+flowchart LR
+    subgraph kids["accepted children (softwareType)"]
+        m_t["measurement<br/>temperature"]:::ctrl
+        m_o["measurement<br/>quantity (oxygen)"]:::ctrl
+        r_up["reactor<br/>upstream"]:::unit
+    end
+    m_t -->|temperature.measured.atEquipment| h_meas[engine._connectMeasurement]
+    m_o -->|quantity (oxygen).measured.&lt;pos&gt;| h_meas
+    r_up -.stateChange.-> h_react[engine._connectReactor]
+    h_meas --> reconcile[reconcile T / O2 into engine state]
+    h_react --> pull[pull upstream effluent → Fs/Cs_in]
+    classDef ctrl fill:#a9daee,color:#000
+    classDef unit fill:#50a8d9,color:#000
+```
+
+| softwareType | filter | wired to | side-effect |
+|---|---|---|---|
+| `measurement` | any | `engine._connectMeasurement` | `temperature.measured.atEquipment` → `engine.temperature`. PFR additionally honours `quantity (oxygen).measured.<distance>` → nearest grid cell DO. |
+| `reactor` | upstream | `engine._connectReactor` | Subscribes to upstream reactor's `stateChange`; pulls effluent into `Fs[0]` / `Cs_in[0]` before next integration step. |
+
+## 7. Lifecycle — what one `data.clock` advance does
+
+```mermaid
+sequenceDiagram
+    participant clock as clock injector
+    participant reactor as reactor
+    participant engine as kinetics engine
+    participant downstream as settler / next reactor
+    participant out as Port-0 output
+
+    clock->>reactor: data.clock { timestamp }
+    reactor->>engine: updateState(timestamp)
+    Note over engine: n_iter steps,<br/>each timeStep × speedUpFactor
+    engine->>engine: integrate ASM3 rates
+    engine->>engine: emit 'stateChange'
+    reactor->>reactor: notifyOutputChanged
+    reactor->>out: Fluent { inlet=0, F, C[13] }
+    alt PFR
+        reactor->>out: GridProfile { grid, n_x, d_x, … }
+    end
+    out->>downstream: Fluent envelope
+```
+
+`stateChange` re-emits on `reactor.emitter` (BaseDomain emitter) so downstream reactors / settlers can listen. The effluent emission goes through the BaseNodeAdapter tick pipeline.
+
+## 8. Data model — `getOutput()`
+
+Port-0 process payload is the `Fluent` envelope (+ optional `GridProfile` for PFR). Port-1 telemetry is the scalar snapshot below.
+
+<!-- BEGIN AUTOGEN: data-model -->
+
+| Key | Type | Unit | Sample |
+|---|---|---|---|
+| `S_HCO` | number | — | `5` |
+| `S_I` | number | — | `30` |
+| `S_N2` | number | — | `0` |
+| `S_NH` | number | — | `25` |
+| `S_NO` | number | — | `0` |
+| `S_O` | number | — | `0` |
+| `S_S` | number | — | `70` |
+| `X_A` | number | — | `200` |
+| `X_H` | number | — | `2000` |
+| `X_I` | number | — | `1000` |
+| `X_S` | number | — | `100` |
+| `X_STO` | number | — | `0` |
+| `X_TS` | number | — | `3500` |
+| `flow_total` | number | — | `0` |
+| `temperature` | number | — | `20` |
+
+<!-- END AUTOGEN: data-model -->
+
+**Concrete sample** (CSTR mid-integration, nitrifying):
+
+```json
+{
+  "flow_total": 1000,
+  "temperature": 15.2,
+  "S_O": 2.1,
+  "S_I": 30,
+  "S_S": 12.4,
+  "S_NH": 0.8,
+  "S_N2": 4.3,
+  "S_NO": 18.6,
+  "S_HCO": 4.2,
+  "X_I": 1050,
+  "X_S": 65,
+  "X_H": 2150,
+  "X_STO": 4.5,
+  "X_A": 215,
+  "X_TS": 3680
+}
+```
+
+Species ordering follows ASM3: indices 0–6 are soluble, 7–12 are particulate. `flow_total` is the effluent flow (m³/d); the reactor uses days as the time unit internally.
+
+## 9. Configuration — editor form ↔ config keys
+
+```mermaid
+flowchart TB
+    subgraph editor["Node-RED editor form"]
+        f1[Reactor type CSTR / PFR]
+        f2[Volume m3]
+        f3[Length m + resolution]
+        f4[Alpha dispersion]
+        f5[KLa 1/h]
+        f6[Time step + speed-up]
+        f7[Initial state 13 species]
+    end
+    subgraph config["Domain config slice"]
+        c1[reactor.reactor_type]
+        c2[reactor.volume]
+        c3[reactor.length<br/>reactor.resolution_L]
+        c4[reactor.alpha]
+        c5[reactor.kla]
+        c6[reactor.timeStep<br/>reactor.speedUpFactor]
+        c7[initialState.* ASM3 keys]
+    end
+    f1 --> c1
+    f2 --> c2
+    f3 --> c3
+    f4 --> c4
+    f5 --> c5
+    f6 --> c6
+    f7 --> c7
+```
+
+| Form field | Config key | Default | Range | Where used |
+|---|---|---|---|---|
+| Reactor type | `reactor.reactor_type` | `CSTR` | enum: `CSTR` / `PFR` | engine selection in `_buildEngine` |
+| Volume (m³) | `reactor.volume` | `1000` | > 0 | residence time, mass balance |
+| Length (m) | `reactor.length` | `10` | > 0 | PFR only — axial extent |
+| Resolution L | `reactor.resolution_L` | `10` | ≥ 1 | PFR grid cell count |
+| Alpha | `reactor.alpha` | `0.5` | 0–1 | dispersion vs plug-flow blend |
+| Inlets | `reactor.n_inlets` | `1` | ≥ 1 | `Fs[]` / `Cs_in[]` array sizes |
+| KLa (1/h) | `reactor.kla` | `0` | ≥ 0 | aeration mass transfer (NaN → use `data.otr`) |
+| Time step (h) | `reactor.timeStep` | `0.001` | ≥ 0.0001 | integrator inner step |
+| Speed-up factor | `reactor.speedUpFactor` | `1` | ≥ 1 | wall-clock → process-time multiplier |
+| Initial S_NH | `initialState.S_NH` | `25` | ≥ 0 (mg/L) | starting ammonium |
+| Initial X_H | `initialState.X_H` | `2000` | ≥ 0 (mg/L) | starting heterotroph biomass |
+| Initial X_A | `initialState.X_A` | `200` | ≥ 0 (mg/L) | starting autotroph biomass — must be ≥ ~50 for nitrification |
+| Initial X_TS | `initialState.X_TS` | `3500` | ≥ 0 (mg/L) | starting TSS — drives settler split |
+
+## 10. State chart
+
+Skipped — reactor has no FSM. It runs continuous-state ODE integration; the engine's only stateful event is `stateChange`, fired after every successful integration advance. See section 7 for the integration sequence.
+
+## 11. Examples
+
+| Tier | File | What it shows | Status |
+|---|---|---|---|
+| Basic | `examples/basic.flow.json` | CSTR with one inlet, watch `Fluent` effluent | ✅ in repo |
+| Integration | `examples/integration.flow.json` | upstream reactor → reactor → settler chain | ✅ in repo |
+| Edge | `examples/edge.flow.json` | PFR with dispersion + multi-inlet | ✅ in repo |
+| Companions | `additional_nodes/*` | recirculation-pump + settling-basin Node-RED nodes shipped from this repo | ✅ in repo |
+
+One screenshot per tier where helpful. PNG ≤ 200 KB under `wiki/_partial-screenshots/reactor/`.
+
+## 12. Debug recipes
+
+| Symptom | First thing to check | Where to look |
+|---|---|---|
+| Nitrification doesn't proceed (S_NH stays high) | `initialState.X_A` must be ≥ ~50 mg/L. Defaulting to `0.001` (a known footgun) means no autotrophs. | `generalFunctions/src/configs/reactor.json` |
+| `Fluent` effluent flow zero | No `data.clock` ticks arriving, or `data.fluent` never set `Fs[0] > 0`. | `commands/handlers.js`, engine `setInfluent` |
+| PFR `GridProfile` not emitted | `reactor_type` set to `CSTR` — only PFR emits grid. | `_buildEngine` switch |
+| Settler downstream not updating | `stateChange` event listener path: settler must subscribe to `reactor.emitter`, NOT `reactor.measurements.emitter`. | settler `_connectReactor` |
+| Temperature reconcile silently ignored | Child measurement's `asset.type` not `temperature` exactly, or `positionVsParent` not `atEquipment`. | `engine._connectMeasurement` |
+| Integrator slow / stalls | `reactor.timeStep` too small for `speedUpFactor`. Internal `n_iter` count blows up. | `engine.updateState` |
+| `wiki:datamodel` script slow / hangs | `mathjs` cold-start ~13 s; instantiation depends on it transitively. See known-limitations row 1. | `kinetics/baseEngine.js` |
+
+## 13. When you would NOT use this node
+
+- Use reactor for **ASM3 biological treatment** modelling (activated sludge, nitrification, denitrification). For aerobic-only or simpler kinetics, the ASM3 species vector is overkill.
+- Don't use reactor for a passive equalisation tank — the kinetics engines assume reactions are happening.
+- Skip reactor when you only need a residence-time delay; a simple buffer node is lighter and doesn't require `mathjs`.
+
+## 14. Known limitations / current issues
+
+| # | Issue | Tracked in |
+|---|---|---|
+| 1 | `mathjs` cold-start adds ~13 s to first `require()` — `wiki:datamodel` auto-gen may time out on the 60 s wrapper. Falls back to the hand-curated `concrete sample` block. | `.claude/refactor/OPEN_QUESTIONS.md` — "mathjs slow load" |
+| 2 | `initialState.X_A` default of `200` mg/L is correct; older config snapshots used `0.001` which silently disabled nitrification. Verify on every new deploy. | `generalFunctions/src/configs/reactor.json` |
+| 3 | `getEffluent` shape historically varied (array vs single envelope) — settler's `_connectReactor` tolerates both. Don't break the contract without updating settler. | `nodes/settler/src/specificClass.js → _connectReactor` |
+| 4 | `additional_nodes/recirculation-pump` and `settling-basin` are legacy companions — not yet refactored to BaseDomain. | P6.5 follow-up |
+| 5 | `reaction_modules/` is a legacy plug-in directory not consumed by the current engines. Removal pending. | P6.5 follow-up |