{"title":"Chemical timescale effects on detonation convergence","authors":"Shivam Barwey , Michael Ullman , Ral Bielawski , Venkat Raman","doi":"10.1016/j.compfluid.2025.106550","DOIUrl":null,"url":null,"abstract":"<div><div>Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next-generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver (Sharma et al., 2024)—which enables resolutions up to <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mn>1000</mn><mo>)</mo></mrow></mrow></math></span> cells-per-induction length—the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen-air mechanism (Mével et al. (2009), Shepherd (2018)), where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. However, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"289 ","pages":"Article 106550"},"PeriodicalIF":2.5000,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers & Fluids","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0045793025000118","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next-generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver (Sharma et al., 2024)—which enables resolutions up to cells-per-induction length—the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen-air mechanism (Mével et al. (2009), Shepherd (2018)), where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. However, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.
含爆流的数值模拟已经成为设计下一代动力和推进装置的重要工具。随着这些工具的成熟,燃烧界在模拟爆轰时正确理解和隔离网格分辨率影响是很重要的。为此,本工作的目的是对非定常爆轰模拟的数值收敛性进行全面分析,重点是在算子分裂的情况下分离化学时间尺度修改对收敛特性的影响。借助基于amrex的自适应网格细化流求解器(Sharma等人,2024)-其分辨率可达O(1000)个细胞-每个诱导长度-使用两种动力学配置进行收敛分析:(1) Short and Quirk(1997)简化的基于arrhenius的三步模型机制,通过调整引发反应和分支反应中的活化能来修改爆炸中的化学时间尺度;(2)详细的氢-空气机制(m等人(2009),Shepherd(2018)),通过改变环境压力来调整化学时间尺度。参考稳态理论基线解,利用这些机制分析了一维通道中非定常自持续爆轰的收敛性。分析的目的是提供网格分辨率对宏观(峰值压力和波速)和微观(波结构)量的影响的详细比较,并在偏离稳态基线和最小化学时间尺度之间建立联系。特别是,发现化学时间尺度的缩减对宏观性质收敛的影响最小。然而,对微观收敛趋势的分析,特别是在反应前沿位置,揭示了一个关键的见解:通过机制简化和燃烧建模,在保持诱导时间的同时消除令人望而却步的化学时间尺度,可以显著提高爆轰收敛性能。最后,本研究揭示了与分辨率相关的非定常爆轰收敛机制,并强调了化学时间尺度以及化学时间尺度与感应时间之间的比值对爆轰波结构数值收敛的重要作用。
期刊介绍:
Computers & Fluids is multidisciplinary. The term ''fluid'' is interpreted in the broadest sense. Hydro- and aerodynamics, high-speed and physical gas dynamics, turbulence and flow stability, multiphase flow, rheology, tribology and fluid-structure interaction are all of interest, provided that computer technique plays a significant role in the associated studies or design methodology.