1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
//! # Reactor Boundary Value Problem (BVP) Module
//!
//! This module provides solvers for steady-state gas-phase combustion and plug-flow reactor problems
//! using boundary value problem formulations with dimensionless variables.
//!
//! ## Mathematical Model
//!
//! ### Nomenclature
//!
//! | Symbol | Description | Units |
//! |--------|-------------|-------|
//! | `l` | Length scale | m |
//! | `dT` | Temperature displacement (dimensionless) | K |
//! | `Ts` | Temperature scale | K |
//! | `Tm` | Characteristic temperature at combustion front | K |
//! | `Cp` | Heat capacity | J/(kg·K) |
//! | `λ` | Thermal conductivity | W/(m·K) |
//! | `ρ` | Gas density | kg/m³ |
//!
//! ### Model Assumptions
//!
//! - Heat capacity, thermal conductivity, thermal effects, and gas density are constants
//! calculated at characteristic temperature `Tm`
//! - Only elementary kinetic expressions are considered (no fall-off, third-body reactions)
//! - reactor is one-dimensional with coordinate `x`
//! - Steady-state conditions (front velocity `m` is constant)
//! - reactions occur in the gas phase only
//! - reactor is adiabatic (no heat loss to surroundings)
//! - ideal gas behavior
//! - no radiation effects
//! - all reactions are irreversable
//!
//! ### Governing Equations
//!
//! **Dimensional form:**
//! ```text
//! m(dc_i)/dx = d/dx[ρD(dc_i)/dx] + Σ_j[ν_ij w_j]
//! mCp(dT)/dx = d/dx[λ(dT)/dx] + Σ_j[Q_j w_j]
//! ```
//!
//! **Dimensionless form:**
//!
//! With dimensionless variables `Θ = (T-dT)/Ts` and `ξ = x/l`:
//!
//! ```text
//! dΘ/dξ = q/λ
//! dq/dξ = Pe_T·q - a_T·Σ_j[Q_j w_j(Θ)]
//! dc_i/dξ = j_i/(ρD)
//! dj_i/dξ = Pe_c·j_i - a_c·Σ_j[ν_ij w_j]
//! ```
//!
//! Where:
//! - `Pe_T = lmCp/λ` (thermal Peclet number)
//! - `Pe_c = ml/(ρD)` (mass Peclet number)
//! - `a_T = l²/Ts`, `a_c = l²`
//!
//! ### Boundary Conditions
//!
//! - At `ξ = 0`: `c_i = c_i0`, `Θ = Θ0 = (T0-dT)/Ts`
//! - At `ξ = ∞`: `q = j_i = 0` (zero fluxes)
//!
//! ## Numerical Solution
//!
//! The system forms a boundary value problem with `2n+2` first-order ODEs for `n` substances.
//! Solutions are obtained using Newton-Raphson methods from the [RustedSciThe](https://crates.io/crates/RustedSciThe)
//! package.