Efficient two-dimensional simulation of primary reference fuel ignition under engine-relevant thermal stratification

Despite vast research on engine knock, there remains a limited understanding of the interaction between reaction front propagation, pressure oscillations, and fuel chemistry. To explore this through computational fluid dynamics, the adoption of advanced numerical methods is necessary. In this context, the current study introduces ARCFoam, a computational framework that combines dynamic mesh balancing, chemistry balancing, and adaptive mesh refinement with an explicit, density-based solver designed for simulating high-speed flows in OpenFOAM. First, the validity and performance of the solver are assessed by simulating directly initiated detonation in a hydrogen/air mixture. Second, the study explores the one/two-dimensional (1D/2D) hotspot ignition for the primary reference fuel and illuminates the impact of transitioning to 2D simulations on the predicted combustion modes. The 2D hotspot simulations reveal a variety of 2D physical phenomena, including the appearance of converging shock/detonation fronts as a result of negative temperature coefficient (NTC) behavior and shock wave reflection-induced detonation. The main results of the paper are as follows: (1) NTC chemistry is capable of drastically changing the anticipated reaction front propagation mode by manipulating the local/global reactivity distribution inside and outside the hotspot, (2) subsonic hotspot ignition can induce detonation (superknock) through the generation of shock waves and subsequent wall reflections, and (3) while the 1D framework predicts the initial combustion mode within the hotspot, significant differences between 1D and 2D results may emerge in scenarios involving ignition-to-detonation transitions and curvature effect on shock/detonation front propagation.