2. Nomenclature

This is the single, authoritative symbol table for the Pyrolysis.jl Technical Reference Guide and User Guide. Every writer MUST conform to these symbols, meanings, and SI units verbatim. Where a symbol is used differently in different source subsystems, the conflict has been resolved here and the resolution is binding. Conflicts and overloads are called out explicitly in the "Overloaded / reserved symbols" section at the end.

Conventions:

  • The spatial coordinate is z (axial, through-thickness). z = 0 is the bottom / substrate (boundary id 2, tag :bottom); z = L is the top / exposed surface (boundary id 1, tag :top). Heat enters from the top; the material shrinks downward. The transverse/lateral plane carries the cross-sectional area A.
  • Sign convention for fluxes: positive flux = transport in the +z direction (bottom → top). Divergence (F_R − F_L)·A/V is positive when the quantity flows out of a cell.
  • Heat-of-reaction sign convention: h > 0 endothermic (cools), h < 0 exothermic (heats); the volumetric heat source is Q_rxn = −Σ h_r r_r, so endothermic reactions give Q_rxn < 0. (This is the storage convention used internally; note ThermaKin/Gpyro publish h > 0 = exothermic — see overload note H1.)
  • "Bulk" density ρ_j is the pure-phase density used by mixing rules and mixture density; "intrinsic/skeletal" density ρ_{i,j} (cell-wall density) is used ONLY in the porosity formula.
  • Mass concentration ξ_j (kg/m³) is the primary species state variable, not mass fraction. ξ_j = Y_j ρ.

1. Latin letters

SymbolMeaningSI units
ACross-sectional area of the 1D column (may vary in time via lateral-shrinkage law); also face area in the divergence operator
A_0Initial / calibration-reference cross-sectional area (fixed at setup)
A_{αβ}, A_iArrhenius pre-exponential factor for reaction i (or αβ)s⁻¹ (1st-order); m³·kg⁻¹·s⁻¹ (bimolecular)
BSpalding mass-transfer number (liquid evaporation, comparison only)
c_p, c_{p,j}Specific heat capacity (of component j)J·kg⁻¹·K⁻¹
cSpeed of light (radiation energy density diagnostic only)m·s⁻¹
D_{AB}, D_jBinary (Chapman–Enskog) diffusion coefficientm²·s⁻¹
D_radP1 radiation diffusion coefficient, 1/(3(α+σ_s))m
d, d_L, d_R, d_{LR}Distance from cell center to face / between adjacent cell centersm
E, E_{αβ}, E_iActivation energy for reaction iJ·mol⁻¹
E_matrixMatrix-only sensible energy of the domain (gas excluded)J
E_totalTotal sensible energy of the domain (matrix + gas)J
eEuler's number / base of natural log (literal)
FRadiative view factor
F_L, F_RGeneric face flux (left/right) in the divergence operator(flux units)·m⁻²
fGeneric field, or Newton-residual function f(T_s)varies
GIncident radiation intensity (P1 primary variable)W·m⁻²
G_extAmbient radiation field at boundary, 4σT_∞⁴W·m⁻²
gGravitational acceleration (Darcy buoyancy, comparison only)m·s⁻²
H(·)Heaviside step (implemented as smooth tanh ramp)
h (kinetics)Heat of reaction, per kg of first reactantJ·kg⁻¹
h_convConvective heat-transfer coefficientW·m⁻²·K⁻¹
h_mExternal mass-transfer (film) coefficientm·s⁻¹
h_PPressure-transfer coefficient (convective pressure BC)m·Pa⁻¹·s⁻¹
Δh, Δh_gSensible enthalpy difference (midpoint-rule integral)J·kg⁻¹
I(z), I_inRadiation intensity along path / at cell entrance (Beer–Lambert)W·m⁻²
I_surface, I_extIncident radiation intensity at domain surface / external sourceW·m⁻²
J_jMass flux of (gas) component j at a facekg·m⁻²·s⁻¹
J_j^diff, J_j^advDiffusive / advective part of the gas fluxkg·m⁻²·s⁻¹
K, κPermeability — use κ; K reserved for Gpyro-comparison text only
k, k_jThermal conductivity (of component j)W·m⁻¹·K⁻¹
k_eff, k_fEffective (mixture) / face-averaged thermal conductivityW·m⁻¹·K⁻¹
k_BBoltzmann constant (reduced temperature T* = k_B T/ε)J·K⁻¹
LDomain thickness (total material length); also characteristic length in convection correlationsm
M, M_jMolar mass (of component j)kg·mol⁻¹
M_refReference molar mass in the Chapman–Enskog combining rules (= 0.029, air; §5.10)kg·mol⁻¹
M_dryTotal initial dry-solid mass in the columnkg
m_jTotal mass of component j in domainkg
m_{dry,i}Initial dry-solid mass in cell i (χ denominator)kg
ṁ_gOutward gas mass flux at surface (transpiration)kg·m⁻²·s⁻¹
NPartial-pressure sum (numerator of ideal-gas pressure)Pa
N_C (NC)Number of chemical components (type parameter)
N_R (NR)Number of reactions (type parameter)
nNumber of finite-volume cells
n_nodes, n_facesNumber of mesh nodes / faces (= n+1 in 1D)
n_{i,j}, n_sReaction order of reaction i w.r.t. reactant j
NuNusselt number (convection BC)
PGas pressure in the porous mediumPa
P_refReference pressure (= 101325)Pa
Q_rxnVolumetric reaction heat source, −Σ h_r r_rW·m⁻³
Q_radVolumetric radiation absorption sourceW·m⁻³
Q_em, Q_absP1 volumetric emission 4σαT⁴ / absorption αGW·m⁻³
q, q_cond(Conductive) heat fluxW·m⁻²
q_BCBoundary-condition heat flux (positive into domain)W·m⁻²
q''', q_rad'''Volumetric radiation source (Beer–Lambert)W·m⁻³
R_g (R, R_gas)Universal gas constant (= 8.314462618)J·mol⁻¹·K⁻¹
r, r_iReaction rate (per unit volume) of reaction ikg·m⁻³·s⁻¹
r (limiters)Slope ratio argument of a flux limiter ψ(r)
SSutherland constant (= 110.4 K, air) in the viscosity lawK
S_j, S_{j,i}Species source term for component j (in cell i)kg·m⁻³·s⁻¹
S_convVolumetric gas-advection (convection) energy sourceW·m⁻³
S_genVolumetric gas-generation enthalpy sink (optional)W·m⁻³
S_radNet volumetric P1 radiation source, α(G − 4σT⁴)W·m⁻³
S_dry,iDry-pyrolysis-gas source term in cell ikg·m⁻³·s⁻¹
TAbsolute temperatureK
T_s, T_wSurface (boundary face) temperatureK
T_cellTemperature at the cell center adjacent to a boundaryK
T_∞Ambient / surroundings temperatureK
T_refReference temperature for enthalpy datum (= 298.15)K
T_face, T_fArithmetic-mean temperature at a faceK
T_min,i, T_max,iLower / upper temperature gate for reaction iK
T*Reduced temperature k_B T / ε_AB (collision integral)
tTimes
uFlat ODE/DAE state vectormixed
duTime derivative of state vectormixed
V, V_iCell volume (of cell i)
v_gBulk (Darcy) gas velocity at a facem·s⁻¹
v_jVolume fraction of component j
v_j^solidSolid volume fraction (excludes non-swelling gases)
WImplicit-RK "W" matrix, I/(γΔt) − J
w, w_iALE mesh node velocity (at node i)m·s⁻¹
w_cellCell-centered ALE mesh velocitym·s⁻¹
X_αVolume fraction (FDS/Gpyro comparison notation; internally v_j)
x, yTransverse/lateral coordinates (multi-D comparison only)m
Y_jMass fraction of component j, ξ_j/ρ
ZPre-exponential factor (Gpyro/FDS comparison notation; internally A)varies
zAxial / through-thickness coordinate (primary spatial variable)m
z_iPosition of mesh node i (ALE state)m
z_i^cCell-center position of cell im
Δz, Δz_iCell thickness (of cell i)m

2. Greek letters

SymbolMeaningSI units
α (α_eff)Volumetric absorption (extinction) coefficientm⁻¹
α_jMass-basis absorption coefficient of component jm²·kg⁻¹
α_convConversion / extent variable (Gpyro kinetics comparison)
βWeighting parameter of the WEIGHTED mixing rule, k = βk_∥ + (1−β)k_series
γ (γ_j)Swelling factor of component j (1 = solids, 0 = non-swelling gas)
γ_RKImplicit-RK stage coefficient (in W = I/(γΔt) − J)
ΔForward-difference / change operator (prefix)varies
δ_{ij}Kronecker delta
δT, δξ_jFinite-difference perturbation stepK; kg·m⁻³
ε, ε_jEmissivity (of component j)
ε_effEffective (composition-weighted) emissivity
ε_ABLennard-Jones well depth (collision integral)J
ζDummy integration variable along the optical pathm
η_iAMR error-indicator value for cell ivaries
θ, θ_iVolumetric strain rate (dilation) of cell i, (1/V) dV/dts⁻¹
θ_AColumn-global radial volumetric strain rate, (1/A) dA/dts⁻¹
θ_LAxial volumetric strain rate (θ − θ_A)s⁻¹
θ_i^jStoichiometric coefficient (ThermaKin comparison notation; internally ν_{i,j})
κ (κ_j, κ_face)Permeability (Darcy) of component j / at a face
κ_sSolid absorption coefficient (FDS comparison notation; internally α)m⁻¹
λ (λ_j, λ_face, λ_eff)Gas transfer (diffusion) coefficientm²·s⁻¹
Characteristic particle/pore diameter (Carman–Kozeny)m
μ, μ_faceGas dynamic (Sutherland) viscosityPa·s
μ_refReference viscosity (= 1.716×10⁻⁵ at 273.15 K, air)Pa·s
νKinematic viscosity (Gpyro Darcy comparison)m²·s⁻¹
ν_{i,j}Stoichiometric coefficient (yield) of component j in reaction i (negative reactant, positive product)
ξ_j, ξ_{j,i}Mass concentration of component j (in cell i) — primary species statekg·m⁻³
tanh(ξ/ξ_th)Depletion rate factor (per reactant)
ρMixture / bulk density (Σ ξ_j for total; mixture-density formula for effective)kg·m⁻³
ρ_jBulk (pure-phase) density of component jkg·m⁻³
ρ_{i,j}Intrinsic / skeletal density of component j (porosity formula only)kg·m⁻³
ρ_gGas-phase density (ideal-gas law)kg·m⁻³
ρc_p^effEffective volumetric heat capacity (matrix only)J·m⁻³·K⁻¹
ρc_p^totalTotal volumetric heat capacity (matrix + gas)J·m⁻³·K⁻¹
σStefan–Boltzmann constant (= 5.670374419×10⁻⁸)W·m⁻²·K⁻⁴
σ_sScattering coefficient (≈ 0 for pyrolysis)m⁻¹
σ_ABLennard-Jones collision diameterÅ
σ_monGaussian width of the interface monitor (r-refinement)m
τ, τ_iOptical thickness (total / cell-local), α Δz
τ_ILUILU drop tolerance
φPorosity (gas void fraction)
χ, χ_iPyrolysis progress in cell i (cumulative dry-gas mass released per unit initial dry-solid mass)
χ̄Mass-weighted column-average pyrolysis progress
ψ(r)Flux-limiter function
Ω_DDiffusion collision integral (Neufeld)
ω (ω(z))Monitor function (r-refinement equidistribution)varies
ω̇'''Volumetric reaction rate (Gpyro comparison notation; internally r)kg·m⁻³·s⁻¹

3. Subscripts

SubscriptMeaning
jComponent / species index (1…N_C)
iCell index (1…n); also reaction index in kinetics context (see overload I1)
r, αβReaction index
gGas phase / gaseous component
sSolid / condensed phase
, liqLiquid phase
L, RLeft / right of a face
fFace-centered quantity
faceFace-averaged property
effEffective (mixture) property
cellCell-centered quantity
s (boundary)Surface (boundary face) — e.g. T_s, ξ_{j,s}
Ambient / surroundings
refReference value
inCell-entrance (radiation)
extExternal source
min, maxLower / upper bound
dryDry-pyrolysis-gas (excludes water vapor)
init, 0Initial value
top, botTop (z=L) / bottom (z=0) boundary
rxnReaction-source contribution
cond, conv, rad, genConduction / convection / radiation / gas-generation contribution
, seriesParallel / series mixing limit
matrix, totalMatrix-only / matrix+gas (energy, heat capacity)

4. Superscripts and accents

MarkMeaning
°, ^o, _initValue at previous time step / initial value
+, (radiation)Forward / backward hemispheric intensity (two-flux)
solidRestricted to solid/liquid (non-swelling) components
diff, advDiffusive / advective flux component
^cCell-center (e.g. z_i^c)
(overbar)Column / spatial average (e.g. χ̄, J̄_g) or Gpyro-style mixture average (ρ̄)
· (overdot)Time rate / flux (e.g. , ω̇)
', '', '''Per-length / per-area / per-volume (FDS/Gpyro comparison notation)

5. Operators

OperatorMeaning
∂/∂tPartial time derivative (Eulerian, fixed point)
∂/∂t|_χALE material-point (moving-mesh) time derivative
d/dtTotal time derivative
, ∇·Gradient / divergence (1D: ∂/∂z, (F_R−F_L)A/V)
∇²Laplacian / second derivative (curvature indicator)
Σ, Sum / integral
clamp(x,a,b)Clamp x to [a,b]
tanhHyperbolic tangent (smooth gates and depletion limiter)
‖·‖_∞Infinity norm
diag(v)Diagonal matrix from vector v
J, A, C_jJacobian / its sparse block / structured coupling j

6. Physical constants (fixed values)

ConstantSymbolValueUnits
Universal gas constantR_g8.314462618J·mol⁻¹·K⁻¹
Stefan–Boltzmann constantσ5.670374419×10⁻⁸W·m⁻²·K⁻⁴
Reference pressureP_ref101325Pa
Reference temperatureT_ref298.15K
Reference molar mass (air)M_ref0.029kg·mol⁻¹
Reference viscosity (air, 273.15 K)μ_ref1.716×10⁻⁵Pa·s
Sutherland constant (air)S110.4K
Kozeny constant (spheres)180

7. Overloaded / reserved symbols — resolutions (BINDING)

These symbols are overloaded somewhere in the source subsystems. The resolution below is authoritative; writers must use the disambiguating subscript/context.

  • χ (chi) — A1 [pyrolysis progress vs. cross-section ratio]. The volume_change, problem_residual, and geometry subsystems define χ as the pyrolysis progress (dimensionless, cumulative dry-gas mass released per unit initial dry-solid mass). The jacobian KB once glosses it as "cross-section area ratio." Resolution: χ is always pyrolysis progress; the cross-section ratio is A/A_0 = law(χ̄) and is never written as χ.

  • θ (theta) — A2 [strain rate vs. stoichiometry]. Internally θ is the volumetric strain rate (s⁻¹), with θ_A radial and θ_L axial. ThermaKin uses θ_i^j for stoichiometric coefficients. Resolution: use ν_{i,j} for stoichiometry everywhere; reserve θ for strain rate. θ_i^j appears only in the ThermaKin-comparison passage and is immediately mapped to ν_{i,j}.

  • A — A3 [area vs. pre-exponential factor]. A is the cross-sectional area; the Arrhenius pre-exponential is A_i / A_{αβ} (always subscripted by reaction). In Jacobian/linear-algebra context A also denotes the sparse Jacobian block. Resolution: bare A = area in physics chapters; A_i = pre-exponential (kinetics); A (linear algebra) only inside the Jacobian chapter where the meaning is unambiguous and stated. The FDS/Gpyro pre-exponential Z is mapped to A_i on first use.

  • κ vs. K — A4 [permeability]. Use κ for permeability throughout. K (Gpyro permeability) appears only in comparison text and is mapped to κ.

  • α — A5 [absorption vs. conversion vs. heat-transfer]. α (or α_eff) = volumetric absorption coefficient (m⁻¹); α_j = mass-basis absorption coefficient (m²·kg⁻¹). Conversion (Gpyro) is α_conv; FDS κ_s maps to α. The mass-transfer "absorptivity" of a surface is also α but always paired with a BC context and Kirchhoff α = ε_eff.

  • λ — A6 [gas transfer vs. anything else]. λ is reserved for the gas transfer (diffusion) coefficient (m²·s⁻¹). Do not use λ for thermal conductivity (that is k) or eigenvalues.

  • ρ subscripts — A7 [bulk vs. intrinsic]. ρ_j = bulk density (mixing, mixture density); ρ_{i,j} = intrinsic/skeletal density (porosity only). Never write a bare ρ_i for skeletal density without the explicit "intrinsic" qualifier; prefer ρ_{i,j}.

  • h — H1 [heat of reaction; sign]. h (kinetics) = heat of reaction per kg of first reactant, with h > 0 endothermic internally; ThermaKin/Gpyro publish h > 0 = exothermic. The conversion to the volumetric source carries the minus sign: Q_rxn = −Σ h_r r_r. The convective coefficient is h_conv, the mass coefficient h_m, the pressure coefficient h_P, all subscripted.

  • i — I1 [cell index vs. reaction index]. i is the cell index in discretization/geometry/ALE; in the kinetics chapter i is the reaction index (with j the component index). Each chapter states its convention in its first paragraph. Where both appear together, use i = cell, r = reaction.

  • φ — A8 [porosity only]. φ is porosity; do not reuse for view factor (that is F) or any potential. Two historical porosity definitions exist (volume-fraction Σ v_j^gas vs. intrinsic-density 1 − Σ ξ_j/ρ_{i,j}); the intrinsic-density form is canonical and is what φ denotes.

  • S — S1 [source vs. Sutherland constant]. S with a subscript (S_conv, S_gen, S_rad, S_j) is a source term; bare S is the Sutherland constant in the viscosity law only.