pith. machine review for the scientific record. sign in

arxiv: 2604.19812 · v1 · submitted 2026-04-16 · ⚛️ physics.chem-ph

Recognition: unknown

An efficient method based on the evolutionary center algorithm for optimizing chemical-diffusive models for flame acceleration and DDT

Huahua Xiao , Xu Zhang , Mingbin Zhao , Congling Shi
Authors on Pith no claims yet

Pith reviewed 2026-05-10 10:23 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords chemical-diffusive modelsevolutionary center algorithmflame accelerationdeflagration-to-detonation transitionparameter optimizationhydrogen combustionNelder-Mead algorithm
0
0 comments X

The pith

A hybrid evolutionary and local search algorithm optimizes parameters for chemical-diffusive models that match detailed chemistry and experiments for flame acceleration and DDT.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces the ECA-NM method to optimize the reaction and diffusion parameters in chemical-diffusive models for use in simulations of flame acceleration and deflagration-to-detonation transitions. This approach combines the broad exploration of the evolutionary center algorithm with the precise tuning of the Nelder-Mead algorithm to achieve much lower errors at reduced computational expense compared to standard genetic algorithms. The resulting models for hydrogen mixtures reproduce key properties from detailed chemical mechanisms across various equivalence ratios and align with experimental observations of flame shapes and speeds. A sympathetic reader would care because such models enable practical, large-scale computations of complex combustion events that would otherwise require prohibitive resources.

Core claim

The central claim is that the ECA-NM hybrid optimization method efficiently determines the optimal parameters for chemical-diffusive models, enabling them to accurately reproduce the major properties of premixed flames and detonations over a wide range of equivalence ratios, while the simulations of flame acceleration and DDT in a channel match experimental results in terms of instabilities, speeds, and transition occurrence, and this method outperforms the traditional genetic algorithm by reducing global error by four orders of magnitude and computational cost by two orders of magnitude.

What carries the argument

The ECA-NM hybrid optimization algorithm that leverages global search from the evolutionary center algorithm and local optimization from the Nelder-Mead method to tune reaction and diffusion parameters in chemical-diffusive models against target combustion wave characteristics.

If this is right

  • The optimized CDMs accurately match flame and detonation properties from detailed mechanisms for hydrogen in air or oxygen.
  • Simulated flame acceleration and DDT show qualitative and quantitative agreement with experiments including tulip flame instabilities.
  • The method provides a way to develop CDMs for quantitative multi-scale simulations of transient flames and detonations.
  • Detailed comparisons confirm superior performance over genetic algorithms in error and cost.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This optimization technique could extend to other simplified combustion models that require parameter fitting for multi-scale problems.
  • Accurate CDMs might improve predictions of explosion hazards in confined spaces by allowing higher-resolution simulations.
  • Validation in additional mixture compositions or geometries would strengthen the models' reliability for engineering applications.

Load-bearing premise

The assumption that a small set of tunable reaction and diffusion parameters in the chemical-diffusive model, once fitted to canonical tests, can capture the full spectrum of flame instabilities and deflagration-to-detonation dynamics without overlooking key chemical or scale-dependent effects.

What would settle it

A direct comparison where the ECA-NM optimized CDM predicts incorrect flame speeds or fails to capture detonation transition in a new experimental setup that detailed chemistry models handle accurately.

read the original abstract

This paper presents an efficient method based on Evolutionary Center Algorithm (ECA) for accurately and efficiently determining the optimal reaction and diffusion parameters for Chemical-Diffusive Models (CDM) to simulate flame acceleration (FA) and deflagration-to-detonation transition (DDT). The proposed method leverages the global search capability of the ECA and the local optimization strength of the Nelder-Mead (NM) algorithm. The hybrid approach (ECA-NM) can efficiently optimize CDM parameters that are capable of accurately reproducing the major properties of combustion waves. The CDMs for premixed flames and detonations of hydrogen in air or oxygen were developed using the present ECA-NM method and validated against canonical tests of combustion waves and previous experiments of FA and DDT. The results show that the major flame and detonation properties calculated using the developed CDMs match those obtained from detailed chemical reaction mechanisms over a wide range of equivalence ratio. The simulated FA and DDT in a channel also agree qualitatively and quantitatively with experiments in terms of complex flame instabilities (e.g., tulip and distorted tulip flames), flame displacement speed, and detonation occurrence. In addition, detailed comparisons to the traditional genetic algorithm demonstrate that the developed ECA-NM method diminishes the global error by four orders of magnitude while reducing the computational cost by two orders of magnitude. This work provides a significantly efficient method for developing chemical-diffusive models that allows quantitative multi-scale simulations of transient flames and detonations in complex scenarios.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents a hybrid Evolutionary Center Algorithm combined with Nelder-Mead (ECA-NM) optimizer for determining reaction and diffusion parameters in Chemical-Diffusive Models (CDMs) to simulate flame acceleration (FA) and deflagration-to-detonation transition (DDT) for hydrogen-air/oxygen premixed flames. It claims that the resulting CDMs reproduce major combustion wave properties from detailed chemical mechanisms across a wide range of equivalence ratios, match experimental FA/DDT observations including tulip flame instabilities and detonation onset, and that ECA-NM reduces global error by four orders of magnitude while cutting computational cost by two orders relative to a traditional genetic algorithm.

Significance. If the quantitative error reductions and experimental agreements hold under scrutiny, the work would provide a practical route to low-cost yet predictive simplified models for multi-scale transient combustion, enabling simulations of complex geometries and DDT scenarios that remain inaccessible to detailed chemistry. The hybrid optimizer's reported efficiency gains could accelerate model development for safety and propulsion applications.

major comments (2)
  1. [Abstract and results] Abstract and results sections: the central claim of a four-order global error reduction and two-order computational cost reduction versus genetic algorithm is load-bearing for the method's value, yet the manuscript provides no explicit definition of the global error metric, no description of the GA baseline implementation (population size, generations, fitness function), and no tabulated convergence data or timing benchmarks to allow verification of these specific factors.
  2. [Validation] Validation sections: the assertions of quantitative agreement with detailed mechanisms 'over a wide range of equivalence ratio' and with experiments on flame displacement speed and detonation occurrence lack accompanying error tables, specific equivalence ratio values tested, channel geometry/initial conditions, or direct comparison plots, which undermines assessment of whether the small CDM parameter set truly captures the dynamics without missing pathways.
minor comments (2)
  1. Define all acronyms (ECA, NM, CDM, FA, DDT) at first use and ensure consistent notation for reaction/diffusion parameters throughout.
  2. Add a summary table of the final optimized CDM parameters for each equivalence ratio to improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We have revised the manuscript to address the concerns about missing definitions, implementation details, and quantitative validation data, thereby strengthening the clarity and verifiability of our claims.

read point-by-point responses

  1. Referee: [Abstract and results] Abstract and results sections: the central claim of a four-order global error reduction and two-order computational cost reduction versus genetic algorithm is load-bearing for the method's value, yet the manuscript provides no explicit definition of the global error metric, no description of the GA baseline implementation (population size, generations, fitness function), and no tabulated convergence data or timing benchmarks to allow verification of these specific factors.

    Authors: We agree that explicit definitions and supporting data are required to substantiate the performance claims. The global error metric is defined as the sum of squared relative errors in key combustion properties (laminar flame speed, detonation speed, induction time, and heat release rate) across all tested conditions, normalized by the number of properties. In the revised manuscript we will insert this definition in the methods section, provide a full description of the GA baseline (population size 50, 100 generations, tournament selection, same fitness function as ECA-NM), and add tabulated convergence histories together with wall-clock timing benchmarks on identical hardware to document the reported four-order error reduction and two-order cost reduction. revision: yes

  2. Referee: [Validation] Validation sections: the assertions of quantitative agreement with detailed mechanisms 'over a wide range of equivalence ratio' and with experiments on flame displacement speed and detonation occurrence lack accompanying error tables, specific equivalence ratio values tested, channel geometry/initial conditions, or direct comparison plots, which undermines assessment of whether the small CDM parameter set truly captures the dynamics without missing pathways.

    Authors: We accept that additional quantitative detail is needed. The revised manuscript will include error tables listing maximum and average relative errors for each property at every equivalence ratio examined (0.5–2.0). We will state the precise channel geometry (length 1.0 m, height 0.1 m, closed at one end) and initial conditions (uniform mixture at 300 K, 1 atm). Direct comparison plots of flame displacement speed versus time and pressure histories at detonation onset will be added, together with quantitative agreement metrics against both detailed-chemistry simulations and the cited experiments. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper optimizes CDM parameters via ECA-NM to match external benchmarks (detailed chemical mechanisms on canonical waves plus direct FA/DDT experiments). Validation targets independent data on flame speeds, instabilities, and detonation onset rather than internal definitions or self-citations. No load-bearing step reduces a claimed prediction to a fitted input by construction, and the optimizer's performance claims are benchmarked against GA on the same external targets.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the premise that CDM parameters can be fitted to reproduce detailed chemistry and experiment; free parameters are the reaction and diffusion coefficients being optimized. No new physical entities are introduced. Standard assumptions about the convergence properties of ECA and NM are invoked without proof.

free parameters (1)
  • reaction and diffusion parameters in CDM
    Multiple scalar parameters per equivalence ratio are adjusted via optimization to match flame speed, detonation speed, and instability behavior from detailed mechanisms.
axioms (1)
  • domain assumption ECA combined with NM can locate parameter sets that make the CDM reproduce major combustion-wave properties across equivalence ratios.
    Invoked when claiming the hybrid optimizer succeeds where genetic algorithms fail.

pith-pipeline@v0.9.0 · 5576 in / 1409 out tokens · 70394 ms · 2026-05-10T10:23:26.486920+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 33 canonical work pages

  1. [1]

    Origins of the deflagration -to-detonation transition in gas -phase combustion

    Oran ES, Gamezo VN. Origins of the deflagration -to-detonation transition in gas -phase combustion. Combustion and Flame 2007;148:4 –47. https://doi.org/10.1016/j.combustflame.2006.07.010

    work page doi:10.1016/j.combustflame.2006.07.010 2007

  2. [2]

    Detonative propulsion.Proceedings of the Combustion Institute, 34(1):125–158, 2013

    Wolański P. Detonative propulsion. Proceedings of the Combustion Institute 2013;34:125–58. https://doi.org/10.1016/j.proci.2012.10.005

    work page doi:10.1016/j.proci.2012.10.005 2013

  3. [3]

    https://www.science.org/doi/10.1126/science.1078129 (accessed February 9, 2026)

    Thermonuclear Supernovae: Simulations of the Deflagration Stage and Their Implications | Science n.d. https://www.science.org/doi/10.1126/science.1078129 (accessed February 9, 2026)

    work page doi:10.1126/science.1078129 2026

  4. [4]

    Astrophysical combustion

    Oran ES. Astrophysical combustion. Proceedings of the Combustion Institute 2005;30:1823–

    2005

  5. [5]

    https://doi.org/10.1016/j.proci.2004.08.278

    work page doi:10.1016/j.proci.2004.08.278 2004

  6. [6]

    Mechanisms and occurrence of detonations in vapor cloud explosions

    Oran ES, Chamberlain G, Pekals ki A. Mechanisms and occurrence of detonations in vapor cloud explosions. Progress in Energy and Combustion Science 2020;77:100804. https://doi.org/10.1016/j.pecs.2019.100804

    work page doi:10.1016/j.pecs.2019.100804 2020

  7. [7]

    Shock focusing and detonation initiation at a flame front

    Xiao H, Oran ES. Shock focusing and detonation initiation at a flame front. Combustion and Flame 2019;203:397–406. https://doi.org/10.1016/j.combustflame.2019.02.012

    work page doi:10.1016/j.combustflame.2019.02.012 2019

  8. [8]

    A Numerical Study of a Two- Dimensional H2-O2-Ar Detonation Using a Detailed Chemical Reaction Model

    Oran ES, Weber JW, Stefaniw EI, Lefebvre MH, Anderson JD. A Numerical Study of a Two- Dimensional H2-O2-Ar Detonation Using a Detailed Chemical Reaction Model. Combustion and Flame 1998;113:147–63. https://doi.org/10.1016/S0010-2180(97)00218-6

    work page doi:10.1016/s0010-2180(97)00218-6 1998

  9. [9]

    Premixed flame stability and transition to detonation in a supersonic combustor

    Goodwin GB, Oran ES. Premixed flame stability and transition to detonation in a supersonic combustor. Combustion and Flame 2018;197:145 –60. https://doi.org/10.1016/j.combustflame.2018.07.008. 33

    work page doi:10.1016/j.combustflame.2018.07.008 2018

  10. [10]

    Comparison of the performance of several recent hydrogen combustion mechanisms

    Olm C, Zsély IG, Pálvölgyi R. Comparison of the performance of several recent hydrogen combustion mechanisms. Combustion and Flame 2014;161:2219 –34. https://doi.org/10.1016/j.combustflame.2014.03.006

    work page doi:10.1016/j.combustflame.2014.03.006 2014

  11. [11]

    Effect of H2 addition on laminar burning velocity of NH3/DME blends by experimental and numerical method using a reduced mechanism

    Li H, Xiao H. Effect of H2 addition on laminar burning velocity of NH3/DME blends by experimental and numerical method using a reduced mechanism. Combustion and Flame 2023;257:113000. https://doi.org/10.1016/j.combustflame.2023.113000

    work page doi:10.1016/j.combustflame.2023.113000 2023

  12. [12]

    Turbulent combustion modeling for internal combustion engine CFD: A review

    Posch S, Gößnitzer C, Lang M, Novella R, Steiner H, Wimmer A. Turbulent combustion modeling for internal combustion engine CFD: A review. Progress in Energy and Combustion Science 2025;106:101200. https://doi.org/10.1016/j.pecs.2024.101200

    work page doi:10.1016/j.pecs.2024.101200 2025

  13. [13]

    Simulations of f lame acceleration and deflagration -to- detonation transitions in methane –air systems

    Kessler DA, Gamezo VN, Oran ES. Simulations of f lame acceleration and deflagration -to- detonation transitions in methane –air systems. Combustion and Flame 2010;157:2063 –77. https://doi.org/10.1016/j.combustflame.2010.04.011

    work page doi:10.1016/j.combustflame.2010.04.011 2010

  14. [14]

    Formation and evolution of distorted tulip flames

    Xiao H, Houim RW, Oran ES. Formation and evolution of distorted tulip flames. Combustion and Flame 2015;162:4084–101. https://doi.org/10.1016/j.combustflame.2015.08.020

    work page doi:10.1016/j.combustflame.2015.08.020 2015

  15. [15]

    The interaction of high -speed turbulence with flames: Global properties and internal flame structure

    Poludnenko AY , Oran ES. The interaction of high -speed turbulence with flames: Global properties and internal flame structure. Combustion and Flame 2010;157: 995–1011. https://doi.org/10.1016/j.combustflame.2009.11.018

    work page doi:10.1016/j.combustflame.2009.11.018 2010

  16. [16]

    The interaction of high-speed turbulence with flames: Turbulent flame speed

    Poludnenko AY , Oran ES. The interaction of high-speed turbulence with flames: Turbulent flame speed. Combustion and Flame 2011;158:301 –26. https://doi.org/10.1016/j.combustflame.2010.09.002

    work page doi:10.1016/j.combustflame.2010.09.002 2011

  17. [17]

    Three -dimensional shock–flame interactions: effect of wall friction

    Xiao H, Shen T, Si T, others. Three -dimensional shock–flame interactions: effect of wall friction. Combustion Theory and Modelling 2023:1–19. 34

    2023

  18. [18]

    Numerical study of the stability of premixed flames propagating in half-open tubes

    Shen T, Xiao H. Numerical study of the stability of premixed flames propagating in half-open tubes. Co mbustion Theory and Modelling 2022;26:774 –95. https://doi.org/10.1080/13647830.2022.2069601

    work page doi:10.1080/13647830.2022.2069601 2022

  19. [19]

    Propagation of premixed hydrogen -air flame initiated by a planar ignition in a closed tube

    Shen T, Li M, Xiao H. Propagation of premixed hydrogen -air flame initiated by a planar ignition in a closed tube. International Journal of Hydrogen Energy 202 2;47:4903–15. https://doi.org/10.1016/j.ijhydene.2021.11.123

    work page doi:10.1016/j.ijhydene.2021.11.123 2021

  20. [20]

    Flame acceleration and deflagration -to-detonation transition in hydrogen- air mixture in a channel with an array of obstacles of different shapes

    Xiao H, Oran ES. Flame acceleration and deflagration -to-detonation transition in hydrogen- air mixture in a channel with an array of obstacles of different shapes. Combustion and Flame 2020;220:378–93. https://doi.org/10.1016/j.combustflame.2020.07.013

    work page doi:10.1016/j.combustflame.2020.07.013 2020

  21. [21]

    DDT run -up distance in uniform and non -uniform hydrogen–air mixtures in sparsely-obstructed channels

    Xiao H, Fan J. DDT run -up distance in uniform and non -uniform hydrogen–air mixtures in sparsely-obstructed channels. Combustion and Flame 2025;273:113919. https://doi.org/10.1016/j.combustflame.2024.113919

    work page doi:10.1016/j.combustflame.2024.113919 2025

  22. [22]

    Choked flame and its transition to detonation in an obstructed channel

    Fan J, Xiao H. Choked flame and its transition to detonation in an obstructed channel. Combustion and Flame 2026;285:114765. https://doi.org/10.1016/j.combustflame.2026.114765

    work page doi:10.1016/j.combustflame.2026.114765 2026

  23. [23]

    Chemical -diffusive m odels for flame acceleration and transition-to-detonation: genetic algorithm and optimisation procedure

    Kaplan CR, Ozgen A, Oran ES. Chemical -diffusive m odels for flame acceleration and transition-to-detonation: genetic algorithm and optimisation procedure. Combustion Theory and Modelling 2019;23:67–86. https://doi.org/10.1080/13647830.2018.1481228

    work page doi:10.1080/13647830.2018.1481228 2019

  24. [24]

    Available: https://doi.org/10.1093/comjnl/7.4.308

    Nelder JA, Mead R. A Simplex Method for Function Min imization. The Computer Journal 1965;7:308–13. https://doi.org/10.1093/comjnl/7.4.308

    work page doi:10.1093/comjnl/7.4.308 1965

  25. [25]

    An efficient and accurate optimization method for the chemical- diffusive model

    Lu X, Kaplan CR, Oran ES. An efficient and accurate optimization method for the chemical- diffusive model. Combustion and Flame 2021;232. 35 https://doi.org/10.1016/j.combustflame.2021.111517

    work page doi:10.1016/j.combustflame.2021.111517 2021

  26. [26]

    A new Evolutionary Structural Optimization method and application for aided design to reinforced concrete components

    Wang L, Zhang H, Zhu M, Chen YF. A new Evolutionary Structural Optimization method and application for aided design to reinforced concrete components. Struct Multidisc Optim 2020;62:2599–613. https://doi.org/10.1007/s00158-020-02626-z

    work page doi:10.1007/s00158-020-02626-z 2020

  27. [27]

    A chemical-diffusive model for simulating detonative combustion with constrained detonation cell sizes

    Lu X, Kaplan CR, Oran ES. A chemical-diffusive model for simulating detonative combustion with constrained detonation cell sizes. Combustion and Flame 2021;230:111417. https://doi.org/10.1016/j.combustflame.2021.111417

    work page doi:10.1016/j.combustflame.2021.111417 2021

  28. [28]

    Cantera: An Object -oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes 2023

    Goodwin DG, Moffat HK, Schoegl I, Speth RL, Weber BW. Cantera: An Object -oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes 2023

    2023

  29. [29]

    Numerical Solution Methods for Shock and Detonation Jump Conditions, 2004

    Browne S, Ziegler J. Numerical Solution Methods for Shock and Detonation Jump Conditions, 2004

    2004

  30. [30]

    Implementing the Nelder-Mead simplex algorithm with adaptive parameters , url =

    Gao F, Han L. Implementing the Nelder-Mead simplex algorithm with adaptive parameters. Computational Optimization and Applications 2012;51:259 –77. https://doi.org/10.1007/s10589-010-9329-3

    work page doi:10.1007/s10589-010-9329-3 2012

  31. [31]

    Comprehensive H2/O2 kinetic model for highpressure combustion

    Burke MP, Chaos M, Ju Y , Dryer FL, Klippenstein SJ. Comprehensive H2/O2 kinetic model for highpressure combustion. International Journal of Chemical Kinetics 2012;44:444–74

    2012

  32. [32]

    An experimental study of distorted tulip flame formation in a closed duct

    Xiao H, wang Q, Shen X, Guo S, Sun J. An experimental study of distorted tulip flame formation in a closed duct. Combustion and Flam e 2013;160:1725 –8. https://doi.org/10.1016/j.combustflame.2013.03.011

    work page doi:10.1016/j.combustflame.2013.03.011 2013

  33. [33]

    Experimental and numerical investigation of premixed flame propagation with distorted tulip shape in a closed duct

    Xiao H, Makarov D, Sun J, Molkov V . Experimental and numerical investigation of premixed flame propagation with distorted tulip shape in a closed duct. Combustion and Flame 2012;159:1523–38. https://doi.org/10.1016/j.combustflame.2011.12.003. 36

    work page doi:10.1016/j.combustflame.2011.12.003 2012

  34. [34]

    2019, Journal of Open Source Software, 4, 1370, doi: 10.21105/joss.01370

    Zhang W, Almgren A, Beckner V , Bell J, Blaschke J, Chan C, et al. AMReX: a framework for block-structured adaptive mesh refinement. Journal of Open Source Software 2019;4:1370. https://doi.org/10.21105/joss.01370

    work page doi:10.21105/joss.01370 2019

  35. [35]

    Navier –Stokes characteristic boundary conditions for real fluids with kinetic-energy- and pressure-equilibrium-preserving schemes

    Li Z, Abdellatif A, Yang R, Jofre L, Capuano F. Navier –Stokes characteristic boundary conditions for real fluids with kinetic-energy- and pressure-equilibrium-preserving schemes. Journal of Computational Physics 2025;535:1140 35. https://doi.org/10.1016/j.jcp.2025.114035

    work page doi:10.1016/j.jcp.2025.114035 2025

  36. [36]

    On the use of immersed boundary methods for shock/obstacle interactions

    Chaudhuri A, Hadjadj A, Chinnayya A. On the use of immersed boundary methods for shock/obstacle interactions. Journal of Computational Physics 2011;230:1731 –48. https://doi.org/10.1016/j.jcp.2010.11.016

    work page doi:10.1016/j.jcp.2010.11.016 2011

  37. [37]

    Flame acceleration and DDT in a channel with fence-type obstacles: Effect of obstacle shape and arrangement

    Liu Z, Li X, Li M, Xiao H. Flame acceleration and DDT in a channel with fence-type obstacles: Effect of obstacle shape and arrangement. International Journal of Hydrogen Energy 2022;39:2787–96. https://doi.org/10.1016/j.proci.2022.08.046

    work page doi:10.1016/j.proci.2022.08.046 2022

  38. [38]

    A New Evolutionary Optimization Method Based on Center of Mass

    Mejía-de-Dios J-A, Mezura-Montes E. A New Evolutionary Optimization Method Based on Center of Mass. In: Deep K, Jain M, Salhi S, editors. Decision Science in Action: Theory and Applications of Modern Decision Analytic Optimisation, Singapore: Springer; 2019, p. 6 5–

    2019

  39. [39]

    https://doi.org/10.1007/978-981-13-0860-4_6

    work page doi:10.1007/978-981-13-0860-4_6