Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLIP”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Stroboscopic Goggles as a Countermeasure for Dynamic Visual Acuity and Landing Sickness After Long-Duration Spaceflight

Gravity transitions cause changes in the vestibulo-occular reflex (VOR), which manifests as poor gaze control, a decrement in dynamic visual acuity (the ability to maintain gaze while in motion), both of which are caused by retinal slip. Retinal slip, the inability to keep an image focused on the retina, can drive or worsen sensory conflict, resulting in motion sickness (MS). Currently 100% of returning crewmembers report MS symptoms, which might affect their ability to perform mission critical tasks immediately after landing. Reschke et al. (2007) demonstrate that stroboscopic vision goggles improve motion sickness onset and symptom severity in motion sickness driven by retinal slip.

Rosenberg, M. J. F.↗
NASA NTRS2017

Mechanism of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This work presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand (AERTS) at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 50 and 80 m/sec. The temperature of the chamber was set at -20˚C. A monotonic droplet generator was used to produce droplets that fell from above, perpendicular to the path of the airfoil. The supercooled state of the droplets was determined by measurement of the temperature of the drops at various locations below the droplet generator exit. A temperature prediction code was also used to estimate the temperature of the droplets based on vertical velocity and the distance traveled by droplets from the droplet generator to the airfoil stagnation line. High speed imaging was employed to observe the interaction between the droplets and the airfoil. The high speed imaging provided droplet deformation information as the droplet approached the airfoil near the stagnation line. A tracking software program was used to measure the horizontal and vertical displacement of the droplet against time. It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip velocity to initial velocity ratios was selected for room temperature and supercooled droplet conditions. The airfoil velocity was 60 m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both environmental conditions indicates that water supercooling has no effect on particle deformation for the limited range of the weak supercooling of the droplets tested and the selected impact velocity. The assumption of a constant surface tension value was further supported by the equal trend of the Bond number obtained for supercooled and room temperature droplets.

in-flight icing↗
NASA NTRS2019

Combined Qualification Vibration Testing and Fixed Base Modal Testing Utilizing a Fixed Based Correction Method

Vibration testing spaceflight hardware is a vital, but time consuming and expensive endeavor. Traditionally modal tests are performed at the component, subassembly, or system level, preferably free-free with mass loaded interfaces or fixed base on a seismic mass to identify the fundamental structural dynamic (modal) characteristics. Vibration tests are then traditionally performed on single-axis slip tables at qualification levels that envelope the maximum predicted flight environment plus 3 dB and workmanship in order to verify the spaceflight hardware can survive its flight environment. These two tests currently require two significantly different test setups, facilities, and ultimately reconfiguration of the spaceflight hardware. The vision of this research is to show how traditional fixed-base modal testing can be accomplished using vibration qualification testing facilities, which not only streamlines testing and reduces test costs, but also opens up the possibility of performing modal testing to untraditionally high excitation levels that provide for test-correlated finite element models to be more representative of the spaceflight hardware's response in a flight environment. This paper documents the first steps towards this vision, which is the comparison of modal parameters identified from a traditional fixed-based modal test performed on a modal floor and those obtained by utilizing a fixed based correction method with a large single-axis electrodynamic shaker driving a slip table supplemented with additional small portable shakers driving on the slip table and test article. To show robustness of this approach, the test article chosen is a simple linear weldment, whose mass, size, and modal parameters couple well with the dynamics of the shaker/slip table. This paper will show that all dynamics due to the shaker/slip table were successfully removed resulting in true fixed-base modal parameters, including modal damping, being successfully extracted from a traditional style base-shake vibration test setup.

Winkel, James P.↗
NASA NTRS2019

Combined Qualification Vibration Testing and Fixed Base Modal Testing Utilizing a Fixed Based Correction Method

Vibration testing spaceflight hardware is a vital, but time consuming and expensive endeavor. Traditionally modal tests are performed at the component, subassembly, or system level, preferably free-free with mass loaded interfaces or fixed base on a seismic mass to identify the fundamental structural dynamic (modal) characteristics. Vibration tests are then traditionally performed on single-axis slip tables at qualification levels that envelope the maximum predicted flight environment and workmanship in order to verify the spaceflight hardware can survive its flight environment. These two tests currently require two significantly different test setups, facilities, and ultimately reconfiguration of the spaceflight hardware. The vision of this research is to show how traditional fixed-base modal testing can be accomplished using vibration qualification testing facilities, which not only streamlines testing and reduces test costs, but also opens up the possibility of performing modal testing to untraditionally high excitation levels that provide for test-correlated finite element models to be more representative of the spaceflight hardware's response in a flight environment. This paper documents the first steps towards this vision, which is the comparison of modal parameters identified from a traditional fixed-based modal test performed on a modal floor and those obtained by utilizing a fixed based correction method with a large single-axis electrodynamic shaker driving a slip table supplemented with additional small portable shakers driving on the slip table and test article. To show robustness of this approach, the test article chosen is a simple linear weldment, whose mass, size, and modal parameters couple well with the dynamics of the shaker/slip table. This paper will show that all dynamics due to the shaker/slip table were successfully removed resulting in true fixed-base modal parameters, including modal damping, being successfully extracted from a traditional style base-shake vibration test setup.

Base-Shake↗
DOE OSTI2012

Magnetic configuration effects on TAE-induced losses and a comparison with the orbit-following model in the Large Helical Device

Fast-ion losses from Large Helical Device (LHD) plasmas due to toroidal Alfven eigenmodes (TAEs) were measured by a scintillator-based lost fast-ion probe (SLIP) to understand the loss processes. TAE-induced losses measured by the SLIP appeared in energy E ranges of around 50-180 keV with pitch angles. between 35 degrees-45 degrees, and increased with the increase in TAE amplitudes. Position shifts of the magnetic axis due to a finite plasma pressure led not only to an increase in TAE-induced losses but also to a stronger scaling of fast-ion losses on TAE amplitudes. Characteristics of the observed fast-ion losses were compared with a numerical simulation based on orbit-following models in which the TAE fluctuations are taken into account. The calculation indicated that the number of lost fast ions reaching the SLIP increased with the increase in the TAE amplitude at the TAE gap. Moreover, the calculated dependence of fast-ion loss fluxes on the fluctuation amplitude became stronger in the case of large magnetic axis shifts, compared with the case of smaller shifts, as was observed in the experiments. The simulation results agreed qualitatively with the experimental observations in the LHD.

Ogawa, K.↗
DOE OSTI2013

Multi-Objective Optimization of a Wrought Magnesium Alloy for High Strength and Ductility

An optimization technique is coupled with crystal plasticity based finite element (CPFE) computations to aid the microstructural design of a wrought magnesium alloy for improved strength and ductility. The initial microstructure consists of a collection of sub-micron sized grains containing deformation twins. The variables used in the simulations are crystallographic texture, and twin spacing within the grains. It is assumed that plastic deformation occurs mainly by dislocation slip on two sets of slip systems classified as hard and soft modes. The hard modes are those slip systems that are inclined to the twin planes and the soft mode consists of dislocation glide along the twin plane. The CPFE code calculates the stress-strain response of the microstructure as a function of the microstructural parameters and the length-scale of the features. A failure criterion based on a critical shear strain and a critical hydrostatic stress is used to define ductility. The optimization is based on the sequential generation of an initial population defined by the texture and twin spacing variables. The CPFE code and the optimizer are coupled in parallel so that new generations are created and analyzed dynamically. In each successive generation, microstructures that satisfy at least 90% of the mean strength and mean ductility in the current generation are retained. Multiple generation runs based on the above procedure are carried out in order to obtain maximum strength-ductility combinations. The implications of the computations for the design of a wrought magnesium alloy are discussed. Research sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U. S. Department of Energy.

Radhakrishnan, Balasubramaniam↗
DOE OSTI2013

Investigating Short Circuit Models for Wind Turbine Generators

Short circuit study needs to be performed at the planning stage before connecting large gen erators with the electric grid. This study will reveal the impact of new generation on the fault levels of the grid. The fault levels will determine the rating of nearby equipment and setting of relays. However, a short circuit model of the new generation is required for this study. This report investigates and quantifies the short circuit models of different wind turbine generators (WTGs) and of wind farms, and using these models, explores the effects of the wind farms on the nearby switchgear. WTGs, commonly classified as Type 1, Type 2, Type 3, and Type 4, exhibit different short circuit behavior depending upon their architecture and controls. The short circuit model (synonymous to Thevenin model) of Type 1 WTG is well known; it comprises of a voltage source behind a transient reactance, X'. The transient reactance X' is derived using the reactance of stator, rotor and the magnetizing reactance. In case of Type 2 WTG, an external rotor resistance is used to increase the range of operating slips. The value of the external rotor resistance increases at higher operating slips. The effect of the external rotor resistance on the fault current response is analyzed by simulating a validated model of Type 2 WTG. The validation is performed using the field data from an actual 67-machine Type 2 wind farm located in Wyoming, USA. The results shows that the external rotor resistance should be considered in its Thevenin impedance for better accuracy. A wind farm contains a large number of WTGs, pad mounted transformers, and cables typically connected in a ‘daisy chain” configuration. The effect of these cable impedances on the short circuit model of Type 2 wind farm is examined using the validated models. It is found that the collector circuit cables inside the farm do not have a significant effect on the short circuit model of the farm. Based on this finding, a simple equivalencing procedure for the wind farm is proposed and validated. Faults inside and outside the wind farm are analyzed and recommendations to select main circuit breakers and generator circuit breakers are made. Type 3 WTG, also known as Doubly Fed Asynchronous Generator (DFAG), uses back to back converters in the rotor circuit to inject slip frequency voltage into the rotor circuit. Until recently, these machines used crowbar to short the rotor circuit during faults to protect the converter com ponents form high rotor currents resulting from faults. This type of WTG is simulated and the model is validated using field data from a 66-machine Type 3 wind farm in Wyoming, USA. The simulation study reveals that the crowbar feature is activated almost immediately after the fault takes place, and remains activated throughout the fault-period. This confirms that the short circuit model for Type 3 WTG with crowbar is similar to the model for Type 2 WTG. In the past, wind farms were allowed to disconnect during faults. However, recent grid codes around the world and in the United States demand that the wind farms be kept connected to the grid during faults, and provide certain reactive power quantified based on the voltage dip at the point of interconnection. This requirement, also called low voltage ride through (LVRT), can be accomplished by using Type 4 WTGs that employ full back to back converter interface between generators (permanent magnet synchronous generators are typically used) and the grid. However, the controls are proprietary, and hence the short circuit models are not well defined. This report proposes controls that enable the WTG to comply with the strictest grid codes, and based on this model, proposes a short circuit model of the Type 4 WTG. This model is not validated due to lack of validation data, but the use of grid code requirements as control objectives lends good credibility to the model. PSCAD/EMTDC ® , a powerful and time-tested time-domain simulation software is used to implement all models. This software is capable of modeling in detail all components of wind farms, and all components and controls of converters in time domain.

17 WIND ENERGY↗
DOE OSTI2014

Anode protection system for shutdown of solid oxide fuel cell system

An Anode Protection Systems for a SOFC system, having a Reductant Supply and safety subsystem, a SOFC anode protection subsystem, and a Post Combustion and slip stream control subsystem. The Reductant Supply and safety subsystem includes means for generating a reducing gas or vapor to prevent re-oxidation of the Ni in the anode layer during the course of shut down of the SOFC stack. The underlying ammonia or hydrogen based material used to generate a reducing gas or vapor to prevent the re-oxidation of the Ni can be in either a solid or liquid stored inside a portable container. The SOFC anode protection subsystem provides an internal pressure of 0.2 to 10 kPa to prevent air from entering into the SOFC system. The Post Combustion and slip stream control subsystem provides a catalyst converter configured to treat any residual reducing gas in the slip stream gas exiting from SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2016

Fermilab Recycler Collimation System Design

To provide 700 kW proton beams for neutrino production in the NuMI facility, we employ slip stacking in the Recycler with transfer to the Main Injector for recapture and acceleration. Slip stacking with 12 Booster batches per 1.33 sec cycle of the Main Injector has been implemented and briefly tested while extensive operation with 8 batches and 10 batches per MI cycle has been demonstrated. Operation in this mode since 2013 shows that loss localization is an essential component for long term operation. Beam loss in the Recycler will be localized in a collimation region with design capability for absorbing up to 2 kW of lost protons in a pair of 20-Ton collimators (absorbers). This system will employ a two stage collimation with a thin molybdenum scattering foil to define the bottom edge of both the injected and decelerated-for-slipping beams. Optimization and engineering design of the collimator components and radiation shielding are based on comprehensive MARS15 simulations predicting high collimation efficiency as well as tolerable levels of prompt and residual radiation. The system installation during the Fermilab 2016 facility shutdown will permit commissioning in the subsequent operating period.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Building the Next Generation of Regional Geothermal Potential Maps: Examples from the Great Basin Region, Western USA

Over the past decade research has greatly advanced our understanding of structurally controlled geothermal systems and the potential for “deep stratigraphic reservoirs” in the Great Basin region, western USA. For example, studies have shown that structurally controlled geothermal systems are associated with specific structural settings, such as terminations of major normal faults, accommodation zones, pull-aparts in strike-slip faults, displacement transfer zones, or step-overs in range-front faults. However, we know that not every structure (e.g., termination of a major normal fault or step-over in a range-front fault) hosts an active geothermal system. Recent research has also shown that regional strain rate model as determined from GPS geodesy, the composition of that strain (extensional versus shear or compressive strain), recency of Quaternary faulting, and slip and dilation tendency of faults correlate positively with the distribution of both high and low-temperature, structurally controlled geothermal systems within the Great Basin. Combining these strain data and fault slip history data sets with existing temperature gradient data, heat flow data, spring temperatures, fluid geothermometry, and geophysical data provides an effective platform for evaluating the geothermal potential of respective structural plays. Recent studies focused on deep stratigraphic reservoirs indicate the best conditions for economically viable development are associated with thick sedimentary basins with low-thermal-conductivity sedimentary cap rocks overlying permeable host rocks (e.g. carbonates) and co-located with areas of high heat flow. Some of the factors that contribute to the potential for deep stratigraphic reservoirs and conventional structurally controlled geothermal systems apply to both types of resources, and in some places the potential for these two types of resources overlap in space, suggesting the potential for hybrid reservoirs (structurally controlled deep stratigraphic reservoirs). Effective and complete evaluation of geothermal potential in the Great Basin should include parameters that cover both types of resources. In this study, two areas with a combined area of 26,000 km 2 were evaluated for geothermal resource potential, equating to <4% of the total 680,000 km 2 size of the Great Basin region. Within these two study areas, 73 potential structural targets and 8 potential deep sedimentary target areas were identified. Four of the 73 potential structural target areas are associated with known thermal anomalies. All 73 structural target areas were ranked according to multiple data sets.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Geothermal Exploration North of Mount St. Helens: Washington State Play-Fairway Project

Active seismicity within a broad zone in proximity to volcanism in the St. Helens Seismic Zone north of Washington State’s Mount Saint Helens provide key ingredients for hydrothermal circulation at depth. This seismicity extends well north of the volcanic edifice below where several faults and associated fractures in outcrop record repeated slip, dilation and alteration indicative of localized fluid flow. Candidate reservoir rocks include marine metasediments overlain by extrusive volcanics. This study characterizes the faults and fractures in the extrusives which are available in outcrop and sampled along well TGH17-24 and Plate Boundary Observatory strainmeter borehole, B201. Outcrop mapping and borehole image log analysis down to 250 m document abundant volcaniclastic deposits and basalt flows that are highly fractured, although intervening ash layers truncate the vertical extent of many of these structures. However, strike slip faults with well-developed gouge zones and adjacent localized high fracture density provide potential for vertically extensive fracture flow paths. Where these rocks are altered by intrusion, they develop high fracture density that can support heat transfer. Borehole image logs further record breakout that indicate the azimuth of SHmax in the shear zone is approximately 054/234°±23.6°, consistent with seismically active strike slip faults in the seismic zone, faults mapped in outcrop and evident in the borehole image logs, and with of plate convergence data. Furthermore, temperature profiles in two wells indicate isothermal conditions at average groundwater temperatures, consistent with rapidly flowing water localized within fractures.

drill↗
DOE OSTI2012

A fourth polymorph in the family of BEDT-TTF salts with thiocyanatocuprate(I) anions: (BEDT-TTF){<_4}Cu{<_3}(NCS){<_5}.

A new crystallographic modification was found in the family of BEDT-TTF salts containing thiocyanatocuprate(I) anions. The cation radical salt crystallizes in the P{sup {bar 1}} space group with a = 8.2720(9) {angstrom}, b = 19.099(2) {angstrom}, c = 43.758(5) {angstrom}, {alpha} = 81.088(4){sup o}, {beta} = 87.253(4){sup o}, {gamma} = 83.790(4){sup o}, V = 6786.2(13) {angstrom}{sup 3} at 150 K. The BEDT-TTF radical cations pack in an {alpha}-type motif, but every four columns, there is a slip in the layer. Where this slip occurs, the stacks are coupled in a {beta}'' manner. The material exhibits semiconductive behavior with an activation energy of 0.13 eV above 250 K and 0.22 eV below. Packing motif of the (BEDT-TTF){sub 4}Cu{sub 3}(NCS){sub 5} structure.

36 MATERIALS SCIENCE↗
DOE OSTI2012

The enigma of post-perovskite anisotropy: deformation versus transformation textures

The D'' region that lies just above the core mantle boundary exhibits complex anisotropy that this is likely due to preferred orientation (texturing) of the constituent minerals. (Mg,Fe)SiO{sub 3} post-perovskite is widely thought to be the major mineral phase of the D''. Texture development has been studied in various post-perovskite phases (MgSiO{sub 3}, MgGeO{sub 3}, and CaIrO{sub 3}), and different results were obtained. To clarify this controversy, we report on transformation and deformation textures in MgGeO{sub 3} post-perovskite synthesized and deformed at room temperature in the diamond anvil cell. Transformed from the enstatite phase, MgGeO{sub 3} post-perovskite exhibits a transformation texture characterized by (100) planes at high angles to the direction of compression. Upon subsequent deformation, this texture changes and (001) lattice planes become oriented nearly perpendicular to compression, consistent with dominant (001)[100] slip. When MgGeO{sub 3} post-perovskite is synthesized from the perovskite phase, a different transformation texture is observed. This texture has (001) planes at high angle to compression and becomes slightly stronger upon compression. We also find that the yield strength of MgGeO{sub 3} post-perovskite is dependent on grain size and texture. Finer-grained samples exhibit higher yield strength and are harder to induce plastic deformation. Strong textures also affect the yield strength and can result in higher differential stresses. The inferred dominant (001) slip for pPv is significant for geophysics, because, when applied to geodynamic convection models, it predicts the observed anisotropies of S-waves as well as an anti-correlation between P- and S-waves.

36 MATERIALS SCIENCE↗
DOE OSTI2017

In-Process Measurements of Flow Characteristics of Wood Plastic Composites

The flow characteristics of polypropylene filled with high percentages of wood flour (50 and 70 wt%) are investigated at a feasible processing temperature (195 °C). Rotational rheometry in oscillatory mode in conjunction with the Cox–Merz rule is ineffective at this temperature due to insufficient linear viscoelastic region (LVR) size and ease of fibres degradation. An in-line rheometer, directly attached to a single screw extruder, is used to avoid these problems: this method allows measurements in real processing conditions with no LVR issues and with reduced degradation risks. The 70 wt% displays plug flow, thus it has been impossible to obtain the flow curve but only the slip velocity as a function of shear stress. The 50 wt% could be corrected for slip using the Mooney procedure. Comparing the flow curve of this material with the ones of the 30 wt% and the unfilled polypropylene matrix, a single master-curve can be obtained and modelled with a Carreau–Yasuda equation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plastic deformation mechanism of calcium-silicate hydrates determined by deviatoric-stress Raman spectroscopy

Highlights: • C-(A-)S-H samples under stress were characterized at the bond and grain level. • Small angle scattering showed that 720 MPa reduced grain thickness by 30%. • In-situ Raman showed that cross-linked C-A-S-H slips along its intralayer. • However, in-situ Raman showed that C-S-H slips along the interlayer. • In-situ Raman also showed a growth in CaCO{sub 3}, indicating fracture in both samples. Creep of the cement matrix affects the structural stability of concrete. In Portland cements, the creep is largely controlled by the binding phase calcium-(aluminum-)silicate-hydrate, or C-(A-)S-H. This phase has a lamellar structure and under deviatoric stress aligns its c-axis with the principal stress. However, the limiting resistance to this reorientation is unknown at the nanocrystalline level. Small-angle X-ray scattering shows that the lamellae thickness decreases under 100's MPa deviatoric stress. Deviatoric stress Raman spectroscopy shows that there are two ways that this break-up can occur. If the material's silicate chains are cross-linked, then strain in SiO bonds does not increase above certain stresses, indicating a relaxation adjacent to the SiO bond. If the chains are not cross-linked, then the silicate chains are broken up by rastering against each other, introducing defects. These results show that the plastic deformation of C-(A-)S-H is relevant for Portland cement creep.

36 MATERIALS SCIENCE↗
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↗