Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DRILLING”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

No 182 W excess in the Ontong Java Plateau source

We report that small-scale W isotope variations in ancient and modern terrestrial rocks provide insights into Earth's accretion and early differentiation history as well as the long-term evolution of the Earth's mantle. Tungsten isotope studies on such rocks have exploited advances in mass spectrometry, both NTIMS and MC-ICPMS, which now permit the determination of W isotope ratios at unprecedented precision. While W isotope studies performed in different labs by MC-ICPMS and NTIMS generally exhibit excellent agreement, obtaining accurate W isotope data at this level of precision remains analytically challenging. For example, a recent NTIMS study reported a relatively large, +24 ppm excess in 182 W/ 184 W for a Phanerozoic sample from the Ontong Java Plateau (OJP), but no such 182 W/ 184 W anomaly was found in another study by MC-ICPMS. The present study aims to resolve the discrepancy between these two previous studies, and more generally to evaluate the agreement between different recent W isotope studies by NTIMS and MC-ICPMS. To this end, we report new W isotope data for OJP drill core samples obtained by MC-ICPMS. Additionally, the OJP samples analyzed here exhibit no resolvable 182 W/ 184 W excess relative to the standards and most terrestrial rocks. Moreover, the OJP samples as well as the terrestrial rock standards analyzed here exhibit small but variable W isotope variations for ratios involving 183 W, producing coupled variations in both ‘radiogenic’ (i.e., 182 W/ 183 W) and ‘non-radiogenic’ (i.e., 183 W/ 184 W) ratios. These W isotope variations are analytical in origin, induced during sample preparation, and very likely caused by a nuclear field shift isotope fractionation affecting primarily the odd isotope ( 183 W). The recently reported 182 W excess for an OJP sample may result from this nuclear field shift effect, as the NTIMS analyses had to rely on a double normalization involving the 183 W/ 184 W ratio. More generally, these results demonstrate that using 183 W data from any MC-ICPMS or NTIMS study requires a careful quantification of any potential analytical 183 W effect. Nevertheless, once such effects are taken into account, then both 182 W/ 184 W and 183 W/ 184 W can accurately be determined to a very high level of precision.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2019

Regionally Correlated Oxygen and Carbon Isotope Zonation in Diagenetic Carbonates of the Bakken Formation

Diagenetic minerals preserve records of burial processes that overprint records of seawater chemistry and impact reservoir porosity. The Mississippian-Devonian aged Bakken Formation in the Williston Basin is a reservoir rock of economic importance whose productivity is affected by diagenetic carbonates, particularly dolomite-ankerite-series carbonates. To investigate how diagenetic carbonate alteration manifests in the Bakken and how that might change across the basin, a combined δ18O and δ13C isotope dataset for diagenetic carbonates was collected by in situ SIMS analysis of 10-μm spots from nine drill holes covering a ~250 km transect of the middle Bakken member. Observed core-to-rim isotopic variability in these small Fe-zoned dolomites and calcites frequently exceeds 10‰ in both δ 18 O and δ 13 C, indicating significant changes in thermal and chemical conditions during cementation. Individual ankeritic growth bands can be correlated across the basin by systematic similarities in minor-element compositions and isotope ratios. In the central part of the basin, δ 18 O and δ 13 C trends at sub-mm-scale decline consistently towards the rims of ankerite-series carbonates, which is interpreted to reflect mineral growth coincident with rising temperatures and an increasing organic contribution to inorganic carbon during burial. The most abrupt shifts in δ 18 O and δ 13 C (changes as large as 12‰ for δ 18 O and 6.5‰ for δ13C within distances < 15 μm) are observed along the basin margin, and are believed to preserve signals of hydrocarbon expulsion and fluid infiltration along basin margins. Based on the data, we conclude that ankerite-series carbonates preserve records of prolonged thermal and chemical processes that operated basinwide. The isotopic trends presented here show diagenetic rock evolution from shallow dolomite formation to deep burial ankerites, as well as the isotopic responses of carbonates to changing fluids and availability of organic δ 13 C. These results provide important constraints on how carbonate isotope records may be altered during burial in organic-rich sedimentary rocks and emphasize the need for caution in using bulk powder samples for geochemical analysis in carbonate systems.

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

Probing nanopore structure and confined fluid behavior in shale matrix: A review on small-angle neutron scattering studies

Although continued growth in unconventional oil and gas production is generally projected, its long-term growth potential and sustainability have significant uncertainties. A critical problem is the low hydrocarbon recovery rates from shale and other tight formations using the horizontal well drilling and hydraulic fracturing techniques: < 10% for tight oil and ~ 20% for shale gas. Moreover, the production rate for a given well typically declines rapidly within one year. The low recoveries and declining production of shale oil and gas reservoirs are apparently related to the small porosity (a few to a few hundred nm) and low permeability (10 –16 –10 –20 m 2 ) of shale matrix, which make the enclosed hydrocarbon fluids difficult to access. Furthermore, to enhance the hydrocarbon recovery from shale matrix, it is essential to study its nanopore structure and confined fluid behavior.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2025

Macro-level mechanical interlocking: A rapid joining approach for additively manufactured compression molded composite panels

Composite joining typically involves multiple steps, such as drilling and surface treatment, as part of the manufacturing process, which leads to low throughput and long cycle times. In the present study, we demonstrated a macro-level mechanical interlocking (MI) based, rapid joining technique to assemble additively manufactured compression molded (AMCM) panels, enabling the production of parts larger than the mold dimensions. Composite panels made of 20 wt% short carbon fiber reinforced acrylonitrile butadiene styrene (CF/ABS) were joined using MI features of various geometries, namely tree (TR), dovetail (Dov), rectangle 2 (Rect2), and rectangle 1 (Rect1), and their in-plane strength was evaluated. The resultant strength of the tested MI joints reached up to 74 % of the baseline tensile strength (i.e., the ‘no joint’ case). Observations from optical and scanning electron microscopy revealed inadequate polymer diffusion between the adherends, indicating that the joint strength was primarily derived from mechanical interlocking. Additionally, the fracture surfaces exhibited stress-whitening marks, which were characterized using differential scanning calorimetry (DSC). The increase in melting enthalpy suggested local stretching of polymer chains due to MI. Finite element analysis (FEA) indicated that the Rect1 MI feature, which generated the lowest stress concentration, outperformed the others in terms of joint strength, achieving 42 MPa. As a demonstration of the MI joining method, a battery box tray measuring 108 cm × 34 cm using a mold with an effective dimension of 36 cm × 34 cm successfully manufactured, resulting in a part with an area three times larger than the mold. In conclusion, this study presents a promising approach to improving composite joining techniques while minimizing production complexities.

In-plane joining↗
DOE OSTI2015

Energy requirements of the switchable polarity solvent forward osmosis (SPS-FO) water purification process

A model was developed to estimate the process energy requirements of a switchable polarity solvent forward osmosis (SPS FO) system for water purification from aqueous NaCl feed solution concentrations ranging from 0.5 to 4.0 molal at an operational scale of 480 m3/day (feed stream). The model indicates recovering approximately 90% of the water from a feed solution with NaCl concentration similar to seawater using SPS FO would have total equivalent energy requirements between 2.4 and 4.3 kWh per m 3 of purified water product. The process is predicted to be competitive with current costs for disposal/treatment of produced water from oil and gas drilling operations. As a result, once scaled up the SPS FO process may be a thermally driven desalination process that can compete with the cost of seawater reverse osmosis.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Co-treating flue gas desulfurized effluent and produced water enables novel waste management and recovery of critical minerals

Herein this study reports a novel approach of resource recovery from co-managing two geographically co-located and chemically complementary wastewaters using a pilot-scale treatment process. Designed to treat flue gas desulfurized (FGD) effluent from combustion powerplants and produced water (PW) from energy industries, the process consists of soda-ash softening, nanofiltration (NF), and reverse osmosis (RO). Recovered products are barite, calcite, and low-salinity water. Using field-collected waters, the results show that softening at pH 8.5 produces calcite (yield: 30 kg/m 3 treated water), a chemical used as SO ₂(g) scrubbers. NF treatment under an applied pressure of 3.5 MPa yields a permeate stream laden with monovalent ions (water recovery 60%) and a concentrate stream with a sulfate concentration 1.8 times of the feedwater concentration. Mixing the NF concentrate and PW at a volumetric ratio of 1.0 precipitates a high-density barite material (4.1 g/cm 3 , yield: ~7.5 kg/m 3 mixture) – a critical mineral commonly used as a weighting agent in drilling. The RO treatment recovers >64% water as the permeate, which can be readily used as cooling make-up water at the powerplants. The RO concentrate stream can be further processed in a thermal evaporative system for additional water recovery and brine production.

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

Injection and Monitoring at the Wallula Basalt Pilot Project

Continental flood basalts represent one of the largest geologic structures on earth but have received comparatively little attention for geologic storage of CO2. Flood basalt lava flows have flow tops that are porous, permeable, and have large potential capacity for storage of CO2. In appropriate geologic settings, interbedded sediment layers and dense low-permeability basalt rock flow interior sections may act as effective seals allowing time for mineralization reactions to occur. Previous laboratory experiments showed the relatively rapid chemical reaction of CO2-saturated pore water with basalts to form stable carbonate minerals. However, recent laboratory tests with water-saturated supercritical CO2 show that mineralization reactions occur in this phase as well, providing a second and potentially more important mineralization pathway than was previously understood. Field testing of these concepts is proceeding with drilling of the world’s first supercritical CO2 injection well in flood basalt being completed in May 2009 near the township of Wallula in Washington State and corresponding CO2 injection permit granted by the State of Washington in March 2011. Injection of a nominal 1000 MT of CO2 was completed in August 2013 and site monitoring is in progress. Well logging conducted immediately after injection termination confirmed the presence of CO2 predominantly within the upper flow top region, and showed no evidence of vertical CO2 migration outside the well casing. Shallow soil gas samples collected around the injection well show no evidence of leakage and fluid and gas samples collected from the injection zone show strongly elevated concentrations of Ca, Mg, Mn, and Fe and 13C/18O isotopic shifts that are consistent with basalt-water chemical reactions. If proven viable by this field test and others that are in progress or being planned, major flood basalts in the U.S., India, and perhaps Australia would provide significant additional CO2 storage capacity and additional geologic sequestration options in regions of these countries where conventional storage options are limited.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

The Modular Borehole Monitoring Program. A research program to optimize well-based monitoring for geologic carbon sequestration

Understanding the impacts caused by injection of large volumes of CO 2 in the deep subsurface necessitates a comprehensive monitoring strategy. While surface-based and other remote geophysical methods can provide information on the general morphology of a CO 2 plume, verification of the geochemical conditions and validation of the remote sensing data requires measurements from boreholes that penetrate the storage formation. Unfortunately, the high cost of drilling deep wellbores and deploying instrumentation systems constrains the number of dedicated monitoring borings as well as limits the technologies that can be incorporated in a borehole completion. The objective of the Modular Borehole Monitoring (MBM) Program was to develop a robust suite of well-based tools optimized for subsurface monitoring of CO 2 that could meet the needs of a comprehensive well-based monitoring program. It should have enough flexibility to be easily reconfigured for various reservoir geometries and geologies. The MBM Program sought to provide storage operators with a turn-key fully engineered design that incorporated key technologies, function over the decades long time-span necessary for post-closure reservoir monitoring, and meet industry acceptable risk profiles for deep-well installations. While still within the conceptual design phase of the MBM program, the SECARB Anthropogenic Test in Citronelle, Alabama, USA was identified as a deployment site for our engineered monitoring systems. The initial step in designing the Citronelle MBM system was to down-select from the various monitoring tools available to include technologies that we considered essential to any program. Monitoring methods selected included U-tube geochemical sampling, discrete quartz pressure and temperature gauges, an integrated fibre-optic bundle consisting of distributed temperature and heat-pulse sensing, and a sparse string of conventional 3C-geophones. While not originally planned within the initial MBM work scope, the fibre-optic cable was able to also be used for the emergent technology of distributed acoustic sensing. The MBM monitoring string was installed in March, 2012. To date, the Citronelle MBM instruments continue to operate reliably. Results and lessons learned from the Citronelle MBM deployment are addressed along with examples of data being collected.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

FutureGen 2.0 Monitoring Program: An Overview of the Monitoring Approach and Technologies Selected for Implementation

The FutureGen 2.0 Project will design and build a first-of-its-kind, near-zero emissions coal-fueled power plant with carbon capture and storage (CCS). To assess storage site performance and meet the regulatory requirements of the Class VI Underground Injection Control (UIC) Program for CO2 Geologic Sequestration, the FutureGen 2.0 project will implement a suite of monitoring technologies designed to 1) evaluate CO2 mass balance and 2) detect any unforeseen loss in CO2 containment. The monitoring program will include direct monitoring of the injection stream and reservoir, and early-leak-detection monitoring directly above the primary confining zone. It will also implement an adaptive monitoring strategy whereby monitoring results are continually evaluated and the monitoring network is modified as required, including the option to drill additional wells in out-years. Wells will be monitored for changes in CO2 concentration and formation pressure, and other geochemical/isotopic signatures that provide indication of CO2 or brine leakage. Indirect geophysical monitoring technologies that were selected for implementation include passive seismic, integrated surface deformation, time-lapse gravity, and pulsed neutron capture logging. Near-surface monitoring approaches that have been initiated include surficial aquifer and surface- water monitoring, soil-gas monitoring, atmospheric monitoring, and hyperspectral data acquisition for assessment of vegetation conditions. Initially, only the collection of baseline data sets is planned; the need for additional near- surface monitoring will be continually evaluated throughout the design and operational phases of the project, and selected approaches may be reinstituted if conditions warrant. Given the current conceptual understanding of the subsurface environment, early and appreciable impacts to near-surface environments are not expected.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Co-optimization of CO 2 -EOR and Storage Processes under Geological Uncertainty

This paper presents an integrated numerical framework to co-optimize EOR and CO 2 storage performance in the Farnsworth field unit (FWU), Ochiltree County, Texas. The framework includes a field-scale compositional reservoir flow model, an uncertainty quantification model and a neural network optimization process. The reservoir flow model has been constructed based on the field geophysical, geological, and engineering data. A laboratory fluid analysis was tuned to an equation of state and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). A history match of primary and secondary recovery processes was conducted to estimate the reservoir and multiphase flow parameters as the baseline case for analyzing the effect of recycling produced gas, infill drilling and water alternating gas (WAG) cycles on oil recovery and CO 2 storage. A multi-objective optimization model was defined for maximizing both oil recovery and CO 2 storage. The uncertainty quantification model comprising the Latin Hypercube sampling, Monte Carlo simulation, and sensitivity analysis, was used to study the effects of uncertain variables on the defined objective functions. Uncertain variables such as bottom hole injection pressure, WAG cycle, injection and production group rates, and gas-oil ratio among others were selected. The most significant variables were selected as control variables to be used for the optimization process. A neural network optimization algorithm was utilized to optimize the objective function both with and without geological uncertainty. The vertical permeability anisotropy (Kv/Kh) was selected as one of the uncertain parameters in the optimization process. The simulation results were compared to a scenario baseline case that predicted CO 2 storage of 74%. The results showed an improved approach for optimizing oil recovery and CO 2 storage in the FWU. The optimization process predicted more than 94% of CO 2 storage and most importantly about 28% of incremental oil recovery. The sensitivity analysis reduced the number of control variables to decrease computational time. A risk aversion factor was used to represent results at various confidence levels to assist management in the decision-making process. The defined objective functions were proved to be a robust approach to co-optimize oil recovery and CO 2 storage. The Farnsworth CO 2 project will serve as a benchmark for future CO 2 –EOR or CCUS projects in the Anadarko basin or geologically similar basins throughout the world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The clumped-isotope geochemistry of exhumed marbles from Naxos, Greece

Exhumation and accompanying retrograde metamorphism alter the compositions and textures of metamorphic rocks through deformation, mineral–mineral reactions, water–rock reactions, and diffusion-controlled intra- and inter-mineral atomic mobility. Here, we demonstrate that these processes are recorded in the clumped- and single-isotope (δ 13 C and δ 18 O) compositions of marbles, which can be used to constrain retrograde metamorphic histories. We collected 27 calcite and dolomite marbles along a transect from the rim to the center of the metamorphic core-complex of Naxos (Greece), and analyzed their carbonate single- and clumped-isotope compositions. The majority of Δ 47 values of whole-rock samples are consistent with exhumation- controlled cooling of the metamorphic complex. However, the data also reveal that water–rock interaction, deformation driven recrystallization and thermal shock associated with hydrothermal alteration may considerably impact the overall distribution of Δ 47 values. We analyzed specific carbonate fabrics influenced by deformation and fluid–rock reaction to study how these processes register in the carbonate clumped-isotope system. Δ 47 values of domains drilled from a calcite marble show a bimodal distribution. Low Δ 47 values correspond to an apparent temperature of 260 °C and are common in static fabrics; high Δ 47 values correspond to an apparent temperature of 200 °C and are common in dynamically recrystallized fabrics. We suggest that the low Δ 47 values reflect diffusion-controlled isotopic reordering during cooling, whereas high Δ 47 values reflect isotopic reordering driven by dynamic recrystallization. We further studied the mechanism by which dynamic recrystallization may alter Δ 47 values by controlled heating experiments. Results show no significant difference between laboratory reactions rates in the static and dynamic fabrics, consistent with a mineral-extrinsic mechanism, in which slip along crystal planes was associated with atomic-scale isotopic reordering in the calcite lattice. An intrinsic mechanism (enhanced isotopic reordering rate in deformed minerals) is contraindicated by these experiments. We suggest that Δ 47 values of dynamically recrystallized fabrics that form below the diffusion-controlled blocking-temperature for calcite constrain the temperature of deformation. We find that Δ 47 -based temperatures of static fabrics from Naxos marbles are ~60–80 °C higher than commonly observed in slowly cooled metamorphic rocks, and would suggest cooling rates of ~10 5 °CMyr -1 . A similar thermal history is inferred for dolomite marbles from the core vicinity, which preserve apparent temperatures up to 200 °C higher than a typical blocking temperature (~300 °C). This finding could be explained by a hydrothermal event driving a brief thermal pulse and locally resetting Δ 47 values. Rapid cooling of the core-complex region is consistent with a compilation of published cooling ages and a new apatite U–Th/He age, associating the thermal event with the emplacement of a granodiorite pluton at ~12 Ma.

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

Clustered BSRs: Evidence for gas hydrate-bearing turbidite complexes in folded regions, example from the Perdido Fold Belt, northern Gulf of Mexico

We detail previously undocumented but extensive gas hydrate accumulations in the mouth of Perdido Canyon in the northern Gulf of Mexico. The accumulations are located within central parts of structural domes (four-way closures) and are characterized by stacked, high-amplitude bottom simulating reflections (BSRs) that we call clustered BSRs. Seismic data from Perdido Canyon show two clustered BSRs associated with turbidite sequences within two dome folds formed from tectonic folding and salt diapir rise. The northwestern (NW) and southeastern (SE) clustered BSRs have aerial extents of ~25 km2 and 50 km2, respectively. Well log data confirm gas hydrate occurs above the NW clustered BSR, within a 225 m-thick consistently high-resistivity interval that we interpret as gas hydrate in near-vertical fractures and turbidite sands. The SE dome is only drilled at the edge of the BSR; nevertheless, the well log data indicate that a 30 m-thick gas hydrate accumulation is present. Gas chromatographic logs in both domes suggest a gradual transition from predominantly microbial gas below the BSR (500–1000 m below seafloor (mbsf)) to thermogenic gas at 1000–2000 mbsf. Based on the well log data and seismic stratigraphic analysis, we find gas hydrate is concentrated in fractures in marine mud, as well as in the pores of submarine fan turbidities, where saturations reach as high as 75%. An estimate of the total gas hydrate-bound gas volume at standard temperature and pressure is between 0.04 and 0.17 trillion cubic meters (TCM) assuming average hydrate saturation of 5-20% in a ~45 m thick turbidite sand unit above the Perdido Canyon BSR area. Measured BSR extent and gas volume estimates indicate that the NW and SE reservoirs are among the largest gas hydrate occurrences known in the Gulf of Mexico.

58 GEOSCIENCES↗
DOE OSTI2017

Internal and near nozzle measurements of Engine Combustion Network “Spray G” gasoline direct injectors

Gasoline direct injection (GDI) sprays are complex multiphase flows. When compared to multi-hole diesel sprays, the plumes are closely spaced, and the sprays are more likely to interact. The effects of multi-jet interaction on entrainment and spray targeting can be influenced by small variations in the mass fluxes from the holes, which in turn depend on transients in the needle movement and small-scale details of the internal geometry. In this paper, we present a comprehensive overview of a multi-institutional effort to experimentally characterize the internal geometry and near-nozzle flow of the Engine Combustion Network (ECN) Spray G gasoline injector. In order to develop a complete picture of the near-nozzle flow, a standardized setup was shared between facilities. A wide range of techniques were employed, including both X-ray and visible-light diagnostics. The novel aspects of this work include both new experimental measurements, and a comparison of the results across different techniques and facilities. The breadth and depth of the data reveal phenomena which were not apparent from analysis of the individual data sets. We show that plume-to-plume variations in the mass fluxes from the holes can cause large-scale asymmetries in the entrainment field and spray structure. Both internal flow transients and small-scale geometric features can have an effect on the external flow. The sharp turning angle of the flow into the holes also causes an inward vectoring of the plumes relative to the hole drill angle, which increases with time due to entrainment of gas into a low-pressure region between the plumes. In conclusion, these factors increase the likelihood of spray collapse with longer injection durations.

33 ADVANCED PROPULSION SYSTEMS↗
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↗