Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methane”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Skeletal reaction models for methane combustion

A local-sensitivity-analysis technique is employed to generate new skeletal reaction models for methane combustion from the foundational fuel chemistry model (FFCM-1). Here, the sensitivities of the thermo-chemical variables with respect to the reaction rates are computed via the forced-optimally time dependent (f-OTD) methodology. In this methodology, the large sensitivity matrix containing all local sensitivities is modeled as a product of two low-rank time-dependent matrices. The evolution equations of these matrices are derived from the governing equations of the system. The modeled sensitivities are computed for the auto-ignition of methane at atmospheric and high pressures with different sets of initial temperatures, and equivalence ratios. These sensitivities are then analyzed to rank the most important (sensitive) species. A series of skeletal models with different number of species and levels of accuracy in reproducing the FFCM-1 results are suggested. The performances of the generated models are compared against FFCM-1 in predicting the ignition delay, the laminar flame speed, and the flame extinction. The results of this comparative assessment suggest the skeletal models with 24 and more species generate the FFCM-1 results with an excellent accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗
DOE OSTI2024

Gas-phase surface modification to control catalyst structure and yields in methane dehydroaromatization

Methane dehydroaromatization (MDA) is a promising approach for direct methane transformation to aromatics and hydrogen. The benchmark catalyst Mo/H-ZSM-5 struggles to find commercial adoption because of thermodynamically-limited yields and rapid coking on Brønsted acid and molybdenum carbide species, especially on zeolite external surfaces. Here, gas-phase atomic layer deposition (ALD) overcoats H-ZSM-5 external surfaces with SiO 2 or Al 2 O 3 . NH 3 -TPD, HRTEM, and textural properties show that these overcoats exclusively passivate zeolite external surfaces. Under MDA conditions, SiO 2 gives softer coke and increases cumulative benzene yields by 25%, while Al 2 O 3 strongly decreases yields. H 2 -TPR and UV-visible and Raman spectroscopy show how the overcoats redisperse the MoO x precatalysts, especially over multiple deactivation and isothermal oxidative regeneration cycles. Combined with 27 Al-MAS NMR, MoO x redistribution and dealumination are seen as the causes of long-term deactivation over multiple regeneration cycles, and this process continues to occur regardless of the overcoat. Altogether, the deposition of a small amount of silica on the outer surface of Mo/H-ZSM-5 reduces the formation of hard coke, which could be regenerated by milder methods such as hydrogen treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Evolution of Ni-Mo/MgO during catalytic methane pyrolysis to produce base-growth nanotubes

Catalytic pyrolysis of methane is a promising approach for affordable hydrogen production without CO 2 emissions. While this process is thermodynamically appealing compared to steam reforming, the high stability of methane requires severe conditions, making catalyst stability challenging. Here, we report the behavior of a highly promising Ni-Mo/MgO catalyst, which greatly outperforms its Ni/MgO, Mo/MgO, and Ni-Mo/SiO 2 counterparts at atmospheric pressure. At 800°C, nearly 229 g of carbon nanotubes per gram of Ni are produced. We propose that this superior performance results from the phase evolution of the catalyst, which exsolves stable nickel catalyst particles under reaction conditions. We further reveal that molybdenum carbide formation reduces sintering and adheres the active catalytic particles to the support throughout the reaction, enabling catalyst reuse. Several characterization techniques (same-spot TEM, XPS, XRD, and Raman) are employed to examine catalyst morphology at every step, fostering a deeper understanding of its catalytic activity and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deactivation of Mo/H-ZSM-5 in Microwave-Assisted and Thermal-Driven Methane Dehydroaromatization

A better understanding of catalyst deactivation is needed to improve catalyst design and performance in microwave-enhanced methane dehydroaromatization (MDA). Here, this study investigates the deactivation of a molybdenum supported H-ZSM-5 zeolite (Mo/H-ZSM-5) catalyst in MDA under microwave-heated conditions, comparing its performance to that of the same catalyst under conventional heating. While the microwave-assisted (MW) process achieved higher benzene yields, the catalyst experienced faster deactivation due to the selective and rapid deposition of coke within the pores of the zeolite, as confirmed through Brunauer–Emmett–Teller (BET), X-ray diffraction (XRD), ammonia-temperature programmed desorption (NH 3 -TPD), thermal gravimetric analysis (TGA), temperature programmed oxidation (TPO), and X-ray photoelectron spectroscopy (XPS) analyses. The quantification of total coke content via TGA/TPO and surface carbon (XPS) revealed that nearly twice as much coke was deposited on the catalyst under MW conditions compared to that on the conventionally heated material, and the coke exhibited a more conductive and graphitic nature. The accelerated deactivation rates were attributed to the formation of hot spots in the MW system, leading to enhanced coupling with coke formed in situ during the reaction and resulting in increased Mo reduction. Observations indicated that the CO activation used to carburize the catalyst prior to the reaction is not advantageous in the MW heating environment. The presence of large amounts of Mo oxides at elevated temperatures (through hot spots) exposed to methane leads to instability under the reaction conditions. Optimizing the activation environment and improvement of the Mo dispersion within the pores are potential strategies to improve catalyst stability.

Mo/H-ZSM-5 zeolite↗
DOE OSTI2025

MgO Nanostructures on Au(111) as Catalysts for Low-Temperature Methane Activation and C-C Coupling

The selective conversion of methane (CH4) under mild conditions remains challenging due to strong C-H bonds and catalyst coking. We systematically investigated sub-monolayer MgO nanostructures on Au(111), where two-dimensional (2D) MgO islands with stable Mg-O-Au interfaces catalyze low-temperature CH4 activation and C-C coupling. Upon CH4 exposure at 300 K, surface-bound CHx and C2Hx intermediates formed and persisted post-evacuation, indicating robust CHx-O-Mg linkages. Temperature-programmed studies revealed that C-H activation and C-C coupling intensify with heat: the CHx signal grew continuously while the C2Hx signal reached a plateau at 400-500 K. O 1s and Mg 2p attenuation confirmed adsorption of the hydrocarbons on MgO. Catalytic tests at 500 K yielded C2H6 (70%) and C2H4 (30%) without coking, underscoring MgO's role as an active catalyst. These results offer new design principles for developing coke-resistant and low-temperature methane upgrading catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Stable Co-valorization of Carbon Dioxide and Methane via Dynamic Reconstruction of a Metal Oxide Solid Solution Catalyst

Dry reforming of methane (DRM) is a process that converts two greenhouse gases (methane and carbon dioxide) into syngas, a mixture of H 2 and CO, that can lead to a variety of value-added chemicals. Owing to its endothermic nature, high reaction temperatures up to 800 °C are typically required and the grand challenge lies in developing robust catalysts without sintering and coking-induced deactivation during the long-term on-stream operation. Towards this aim, herein, a robust complex oxide-supported NiCu alloy catalyst was generated in situ during DRM. By leveraging the configurational stability of a solid oxide solution precursor, tightly anchored NiCu bimetallic nanoparticles were in situ exsoluted and acted as the active sites in DRM. The as-afforded catalyst exhibited stable performance for DRM due to the ability to repel coke off the surface as the reaction proceeds. Kinetic experiments along with top surface characterization detail the reconstruction behavior of the solid oxide solution under DRM reaction conditions. In conclusion, the fundamental insights from this work provide guidance on generating resistant and flexible catalysts via in situ active sites formation from easily synthesized metal oxide solid solutions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗
DOE OSTI2025

Methane emission hotspots in a boreal forest-fen mosaic potentially linked to deep taliks

Permafrost thaw is transforming boreal forests into mosaics of wetlands and drier uplands. Topographic controls on hydrological and ecological conditions impact methane (CH 4 ) fluxes, contributing to uncertainty in local and regional CH 4 budgets and underlying drivers. The objective of this study was to explore CH 4 fluxes and their drivers in a transitioning boreal forest-fen ecosystem (Goldstream Valley, Alaska, USA). This landscape is characterized by thawing discontinuous permafrost and heterogeneous mosaics of fens, collapse-scar channels, and small mounds of permafrost soils. From a survey in July 2021, observed chamber CH4 fluxes included fen areas with intermediate to very high emissions (29.8–635.3 mg CH 4 m −2 d −1 ), clustered locations with CH 4 uptake (−2.11 to −0.7 mg CH 4 m −2 d −1 ), and three anomalous emission hotspots (342.4–772.4 mg CH 4 m −2 d −1 ) that were located near samples with lower emissions. Some surface and near-surface variables partially explained the spatial variation in CH 4 flux. Log-transformed CH 4 flux had a positive linear relationship with soil moisture at 20 cm depth ( R 2 = 0.31, p -value < 1e-5) and negative linear relationships with microtopography ( R 2 = 0.13, p -value < 0.006) and slope ( R 2 = 0.28, p -value < 2e-5). Methane emissions generally occurred in flat, wet, graminoid-dominated fens, whereas CH 4 uptake occurred on permafrost mounds dominated by feather mosses and woody vegetation. However, the CH 4 hotspots occurred on drier, slightly sloped locations with low or undetectable near-surface methanogen abundance, suggesting that CH 4 was produced in deeper soils. When the hotspot samples were omitted, log-transformed CH 4 flux had a positive linear relationship with near-surface methanogen abundance ( R 2 = 0.29, p -value = 0.0023), and stronger linear relationships with soil moisture, slope, and soil macronutrient concentrations. Our findings suggest that some CH 4 emission hotspots could arise from CH 4 in deep taliks. The inference that methanogenesis occurs in deep taliks was strengthened by the identification of intrapermafrost taliks across the study area using low-frequency geophysical induction. This study assesses surface spatial heterogeneity in the context of subsurface permafrost conditions and highlights the complexity of CH 4 flux patterns in transitioning forest-wetland ecosystems. To better inform regional CH 4 budgets, further research is needed to understand the spatial distribution of terrestrial CH 4 hotspots and to resolve their surface, near-surface, and subsurface drivers.

boreal↗
DOE OSTI2025

Catalyst Deactivation Modes of PdO/γ-Al 2 O 3 Catalysts for Lean Methane Oxidation

PdO/γ-Al 2 O 3 catalysts are one of the most active catalytic components for the complete oxidation of methane. Under reaction conditions, especially in a wet feed, the catalysts suffer severe performance degradation. This study establishes a series of testing protocols to systematically investigate the causes of catalyst deactivation under methane oxidation reaction conditions. Four distinct catalyst deactivation modes are identified. Two of the deactivation modes are directly related to H 2 O, either from the feed gas or as a part of the reaction products, with one (Mode 2) being attributed to the formation of surface hydroxyl groups and the other (Mode 3) to the competitive adsorption of H 2 O on the catalysts. The impact of the two deactivation modes is acute and severe but reversible. In contrast, the other two deactivation modes are gradual and persistent but irreversible. Both modes are induced by CH 4 oxidation reaction, with the impact of a wet feed (Mode 4) being substantially more severe than that of a dry feed (Mode 1). The major cause of the irreversible catalyst deactivation is attributed to surface reconstruction of PdO nanoparticles, which behaves as a passivation layer lowering the number of coordinately unsaturated Pd sites for CH 4 activation. Although the passivation layer is relatively stable against thermal or hydrothermal treatment, it is not completely inert. Formation and partial regeneration of the passivation layer is a highly dynamic process and heavily depends on the reaction temperature: a lower reaction temperature (≤ 450 ℃) can lead to quicker catalyst deactivation; but a higher reaction temperature (between 500 – 550 ℃) can result in a greater extent of catalyst deactivation.

PdO/γ-Al2O3↗
DOE OSTI2025

The Effect of Metal Promoters in an Mo-Supported HZSM-5 Catalyst for Microwave-Assisted Methane Dehydroaromatization to Aromatics

Microwave (MW)-assisted methane dehydroaromatization (MDHA) using an Mo-supported HZSM-5 catalyst (Mo/HZ5) can convert methane into value-added aromatic products in modular microwave reactor systems, enabling producers to generate revenue from an otherwise wasted resource. Modifying the local environments of active Mo species with metal promoters potentially regulates the reaction/deactivation pathways and improves the heating properties of the Mo/HZ5 under microwaves. Herein, metal promoters (M), including monovalent K+ and bivalent Co2+ and Ni2+, were incorporated to form M-Mo/HZ5 and their MDHA performance was investigated.

metal promoters↗
DOE OSTI2026

Developing Methane Emissions Inventories by Fusing Airborne, Satellite, and Modeled Assessments: Comprehensive Surveys of the Anadarko and Haynesville Basin

The objective of this program is to create basin-specific methane emissions inventories of the Haynesville and Anadarko (Woodford shale) basins, based on truly basin-scale data collection and analysis. Emissions of natural gas, predominantly composed of methane, present a significant financial loss to the oil and gas sector and reduce the competitiveness of American energy sources. While large emissions, above 10 kg/h, occur at only a small percentage of oil and gas sites, they account for a disproportionately large volume of emissions. To address this, we executed a concurrent campaign combining cutting-edge bottom-up modeling with top-down aerial and satellite surveys to capture the heavy-tailed distribution of these large emissions to accurately characterize total emissions throughout both regions.

03 NATURAL GAS↗
DOE OSTI2025

Crankcase Gas Rerouting/Filtration System to Reduce Crankcase Methane Emissions from Leanburn NG Engines

The U.S. has the opportunity to transform the energy sector by aggressively expanding the use of abundant, domestic natural gas (NG) as a fuel. With the intrinsic composition of NG having a higher H:C ratio (~3.6) than that of diesel (H:C ~ 1.8), lower levels of CO 2 , the preeminent greenhouse gas (GHG), are generated during combustion. A major challenge is the emissions abatement of the primary NG constituent methane (CH4). Its inherent stability makes it a potent GHG. Lean-burn NG engines used to drive compressors in the NG supply chain, generate electricity, and repower ships are substantial sources of methane emissions.

03 NATURAL GAS↗
DOE OSTI2024

Mechanistic and Kinetic Analysis of Complete Methane Oxidation on a Practical PtPd/Al 2 O 3 Catalyst

A PtPd/Al 2 O 3 catalyst developed for the complete oxidation of methane from the ventilation air of underground coal mines is compared against a model PdO/Al 2 O 3 catalyst. Although the PtPd/Al 2 O 3 catalyst is substantially more active and stable than the model catalyst, the nature of active sites between the two catalysts is deemed to be fundamentally the same based on their response to different feed gas compositions and the evolution of surface CO adsorption complexes during time-resolved CO adsorption DRIFTS experiment. For both catalysts, coordinatively unsaturated Pd sites are considered the active centers for methane activation and the subsequent oxidation reaction. H 2 O competes with CH 4 for the same active sites, resulting in severe inhibition. Additionally, the CH 4 oxidation reaction also causes self-inhibition. Taking both inhibition effects into consideration, a relatively simple kinetic model is developed. The model provides a good fit of the 72 sets of kinetic data collected on the PtPd/Al 2 O 3 catalyst under practically relevant reaction conditions with CH 4 concentration in the range of 0.05–0.4%, H 2 O concentration of 1.0–5.0%, and reaction temperatures of 450–700 °C. Kinetic parameters based on the model suggest that the CH 4 activation energy on the PtPd/Al 2 O 3 catalyst is 96.7 kJ/mol, and the H 2 O adsorption energy is –31.0 kJ/mol. Both values are consistent with the parameters reported in the literature. The model can be used to develop catalyst sizing guidelines and be incorporated into the control algorithm of the catalytic system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Structurally Regenerable High Entropy Aluminate Spinel Catalysts for Dry Reforming of Methane

The dry reforming of methane reaction is a promising means to convert two potent greenhouse gases, methane and carbon dioxide, into industrially valuable synthesis gas. However, the presence of reducing gases and high operating temperatures degrade conventional nickel catalysts via excessive coke formation and particle sintering. These catalysts are not readily regenerated because the oxidative heat treatments employed to remove coke further promote active particle sintering. In this work, we designed high entropy aluminate spinel oxides (MAl 2 O 4 where M = Co, Mg, Ni, and divalent site vacancies in nominal equimolar concentration) as selective and regenerable reforming catalysts. Under reaction conditions, reducible nickel and cobalt cations exsolved from the spinel lattice to form highly selective bimetallic particles on the oxide surface. Instead of sintering, these particles uniquely redissolved back into the aluminate lattice upon reoxidation and regained the original spinel structure. This phenomenon is ascribed to entropic stabilization, wherein an increase in configurational entropy creates a thermodynamic driving force for redispersing supported metal particles back into the multi‐cationic oxide structure. During the dry reforming reaction, nickel atoms similarly exsolved from a NiAl 2 O 4 sample and reduced to form metallic nickel particles. However, subsequent oxidation of this sample promoted sintering and oxidation of the nickel particles to an inactive state. High entropy materials thus provide a unique mechanism of regeneration, which is inaccessible in conventional catalysts.

coke deposition↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.

03 NATURAL GAS↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 2: Experimental determination of hydrogen isotopic equilibrium for ethane-H2 from 30 to 200 °C and propane-H2 from 75 to 200 °C

The stable isotopic compositions of light n-alkanes, including methane, ethane, and propane, are often used to identify the sources and thermal maturity of natural gas samples. Though stable isotopic compositions of these molecules are commonly assumed to be controlled by kinetic isotope effects, recent studies have proposed both carbon and hydrogen isotopic equilibrium may also occur in some samples. Assessing whether samples are in isotopic equilibrium requires knowledge of light alkane equilibrium fractionation factors over geologically relevant temperatures for formation and storage (up to ∼300 °C). In this study, we report experimental results of hydrogen isotopic equilibrium between ethane and H2 from 30 to 200 °C and propane and H2 from 75 to 200 °C. We compare these results with high-level theoretical calculations and provide a preferred polynomial fit to describe equilibrium fractionation factors. Comparison of these fractionation factors with a compilation of ∼500 compiled environmental gas samples supports the proposal that many (∼50%) of these natural gas samples exhibit hydrogen isotopic compositions consistent with having formed in or attained methane-ethane-propane hydrogen isotopic equilibrium over geologically relevant temperatures for formation and storage (50–300 °C).

Turner, Andrew C↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗