V-BO-04 — Uniform zero-order burnout: exact thickness through park cascades

Tags: ale, burnout, depletion, kinetics, exact

References:

  • closed form L(t) = L₀ exp(−r t/ρ_s)

Problem statement

An adiabatic slab at uniform T₀ with one zero-order reaction solid → gas, γgas = 0 and no heat of reaction. The temperature stays at T₀ exactly (no heat source, no gradient), so the volumetric rate r = A·exp(−E/RT₀) is a constant and every cell loses the same fraction r/ρs of its current volume per unit time. Under ALE the solid concentration stays at ρ_s (dilation removes volume, not concentration) and the exact thickness is

L(t) = L₀ · exp(−r t / ρ_s).

Every cell reaches the depletion threshold at the same instant, so the parking machinery meets its worst case: a cascade of park events, one per step, surface-first, halving the live cell count at tk = ln(2^k·20) ρs/r. L(t) and the solid mass must follow the exact law through every cascade, and ξs must stay at ρs to roundoff.

Quantities of interest (n = 40)

QoIvalueexacterrortolerancewithin tolprovenance
thickness follows L₀ exp(−r t/ρ) through the park cascades8.23e-058.23e-051.932e-070.0001yesmax relative error over all output times
solid concentration stays at ρ_s (dilation removes volume only)100010002.274e-161e-08yesmax |ξs − ρs|/ρ_s over live cells and times
isothermal: T stays at T₀6006001.137e-131e-06yesK; no heat of reaction, no gradient
park cascade: 40 → 5 live cells5500.5yes35 park events; every cell hits the threshold at once

V-BO-04 thickness_history

V-BO-04 thickness_history_error

V-BO-04 live_cells_history

Solver configuration

settingvalue
integratorKenCarp4 (default)
use_aletrue
handle_depletiontrue
depletion_threshold0.05
min_cells2
min_thickness1.0e-6
abstol1.0e-9
reltol1.0e-7