Search NASA⌕ Search

SEARCH · Search NASA

Results for “RECIRCULATIVE FLUID FLOW”

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 19 records

The Complexity of Nonlinear Flow and non-Fickian Transport in Fractures Driven by Three-Dimensional Recirculation Zones

Subsurface fluid flow and solute transport are typically described by first-order or linear rate laws. However, deviations from the first-order rate laws (or anomalous behavior) are typical and lead to nonlinear flow and non-Fickian transport phenomena. The shared underlying mechanisms for these “anomalous” flow and transport phenomena, which are actually typical of fractures, are seldom analyzed. We investigated the co-occurrence of nonlinear flow and non-Fickian transport through numerical simulation experiments of pore-scale flow and solute transport processes through three-dimensional natural fractures. In these fractures, 3-D recirculation zones (RZs) developed and grew with increasing flow rates. Flow nonlinearity resulted from the RZ development and growth, which shrank the main flow channel. Moreover, the same RZs also captured solutes from the main flow channel and released them back later. This retardation resulted in non-Fickian transport. Thus, nonlinear flow and non-Fickian transport are intertwined via their shared dependence on RZs.

58 GEOSCIENCES↗

Understanding the role of flow dynamics in thermoacoustic combustion instability

Thermoacoustic combustion instability is one of the most challenging operational issues in several high-performance, low-emissions combustion technologies, including gas turbines, aircraft engines, rockets, and industrial boilers. Driven by the coupling between combustor acoustics and flame heat release rate fluctuations, thermoacoustic combustion instability can lead to reduced operability, increased emissions, and, in the most extreme cases, catastrophic failure of combustor components. The feedback loop between acoustics and combustion is often facilitated by fluid mechanic oscillations, referred to as “velocity coupling,” whereby acoustic oscillations drive flow fluctuations, which in turn create fluctuations in the flame. The character of these fluid mechanic oscillations is highly dependent on the structure of the flow field and the receptivity of the flow to external excitation. Combustor flow fields use features like fluid recirculation and shear to enhance flame holding and reduce emissions, but these are also the same features that can make the flow receptive to acoustic excitation or even drive self-excited oscillations. Here, in this paper, we discuss the basics of thermoacoustic instability with a focus on the role of hydrodynamic oscillations in typical combustor flows. To facilitate this discussion, we explore the hydrodynamic instability characteristics of several key combustor unit flows (wakes, swirling jets, etc.) and show how the hydrodynamic stability of a flow is an important consideration in determining a combustor’s propensity for thermoacoustic oscillations. Several examples of coupling between hydrodynamics and thermoacoustics are discussed to illustrate this important link. The paper concludes by discussing the potential for designing flow fields that are thermoacoustic instability resistant, either through a reduction in the receptivity of the flow or through nonlinear coupling mechanisms by which self-excited flow instabilities can suppress velocity-coupled combustion oscillations.

42 ENGINEERING↗

Geochemical evolution in Cacapon member: Fluid-rock interaction experiments and model insights for Appalachian Basin geothermal development

Here, this study combines recirculated flow-through experimental results conducted for 17 days at 90C and 200 PSI with reactive transport modeling to estimate fluid-rock interactions occurring in a sandy mudstone using an interbedded sandstone-shale sample from the Cacapon Member of the upper Tuscarora Sandstone/lower Rose Hill Formation for the purpose of geothermal exploration. Results suggest that the fluid and rock are likely to be in or near partial-equilibrium after approximately one year. In addition, after >400 h of continuous injection at 0.05 mL/min (5*10 −8 m 3 /min) the reactive front is restricted to the first ∼13% of the 4 cm*1.6 cm experimental rock length, whereas after >9100 h (∼1 year), the reactive front extends to 30% total length. The rate of changes in dissolution or precipitation are however, very minimal, with all major rock-forming minerals having rates <10 −11 mol/L porous media/s. Reservoir rock in the presence of dilute brine as may be the case during operation of an enhanced geothermal system would experience little alteration during the shut-in phase, and possibly up to one year. These results have utility in geothermal exploration for reservoirs at similar temperatures as well as general fluid-sandy mudstone rock interaction geochemistry.

Appalachian Basin↗

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗

Role of Three-Dimensional Swirl in Forced Convection Heat Transfer Enhancement in Wavy-Plate-Fin Channels

The influence of wall-corrugation-induced swirl flow on enhanced forced convection in wavy-plate-fin cores has been investigated. Three-dimensional computational simulations were carried out for steady-state periodically developed air flow (Pr ~ 0.71; 50 ≤ Re ≤ 4000) with channel walls subject to constant-uniform temperature conditions. The recirculation that develops in the wall troughs and grows to have an axially helical character is scaled by the Swirl number Sw. As Sw increases with higher flowrate and/or corrugation severity, tornado-shaped vortices appear in the wave trough region midway of the interfin channel height, then extend longitudinally to encompass majority of the flow channel. The local wall-shear and heat transfer coefficient variations indicate that boundary-layer thinning upstream of the wave peak aids in intensifying momentum and heat transfer. However, the flow recirculation in wall trough impedes heat transfer at low Sw due to flow stagnation but promotes it at high Sw because of the vortices-induced augmented fluid mixing. The effects of this secondary flow are quantified by Φf(or j), which is seen to increase log-linearly as fin corrugation aspect ratio γ and/or fin spacing ratio ζ increases; the influence of cross section aspect ratio α is marginal. Moreover, the pressure drag penalty due to swirl critically affects overall pressure loss, and its proportion remains nearly constant when α varies, but grows as Sw, γ, and/or ζ increases and can be as much as 80% of the total pressure drop.

Engineering↗

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗

Falling particle solar receivers

Falling particle solar receivers, systems, and methods are disclosed that include one non-linear falling particle curtain or two or more falling particle curtains within a solar receiver that receives incident solar radiation. The particles heated in the solar receiver may be used to heat a secondary fluid. In an embodiment, the particles may be recirculated to improve energy capture and thermal efficiency. In other embodiments, an air curtain may be used across the aperture of the receiver, and flow-control devices may be used to evenly spread particles across the width of the receiver inlet. Finally, feed particles may be preheated using heat from the solar receiver.

Ho, Clifford K.↗

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beam-water Interaction System for PFAS Treatment

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants that resist conventional water treatment due to the strength of their carbon fluorine bonds. Electron-beam irradiation has shown potential for PFAS degradation; however, its effectiveness is limited by shallow energy penetration, restricting treatment to the upper region of the water surface. This project addresses this limitation through the design of a compact mixing chamber and recirculation system to enhance surface renewal and improve treatment efficiency. The system, constrained to approximately 10 cm 20 cm 20 cm, was developed using computer-aided design, computational fluid dynamics (CFD), and MATLAB-based hydraulic analysis to evaluate flow behavior, head loss, and pump compatibility. A perforated-baffle geometry and controlled flow range of 1 3 GPM were implemented to promote uniform mixing while maintaining a stable free surface for electron-beam exposure. Modeling results indicate improved flow distribution and reduced dead zones, supporting more consistent interaction between contaminated water and the beam. The final design enhances hydrodynamic performance within system constraints, contributing to scalable electron-beam water treatment technologies for improved PFAS degradation.

Ramos, Gerardo [Northern Illinois U.]↗

Numerical study on aerosol sampling in a nuclear facility duct with a 90-degree elbow

Due to the challenging design requirements, elbows are often unavoidable in duct configuration, and these 90-degree bends introduce swirl, velocity variations, recirculation, and secondary flow. These disturbances make it difficult for nuclear facilities to meet particle sampling standards. A series of numerical analyses are conducted to track aerosols in a nuclear facility duct having a 90-degree elbow with the assistance of computational fluid dynamics (CFD). A turbulence model, a continuity, and a momentum, a discrete phase model, and species transport equations are solved simultaneously to track aerosols in the duct. The effect of turbulence models, turbulent dispersion models, droplet drag model, aerosol amount, aerosol spray configuration, guide vanes, and mixers are investigated. Simulation results are analyzed per relevant testing codes such as DOE-HDBK-1169, ASME AG1, ISO 14644-3, ACGIH, and ANSI/HPS N13.1.

Han, Kai [Savannah River Nuclear Solutions (SRNS),↗

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗

Scale-up of microwave-assisted, continuous flow, liquid phase reactors: Application to 5-Hydroxymethylfurfural production

Microwave (MW) technology can be powerful for the electrification and intensification of chemical manufacturing. However, very few examples of scaled MW processes exist. In this work, we build a scaled, MW-assisted, continuous flow reactor for the processing of liquid phase chemistries with specific application to 5-hydroxymethylfurfural (HMF) production. Here, we demonstrate a corresponding computational fluid dynamics model to simulate the reactor’s temperature profile and performance, both of which are in good agreement with experiments. We construct a surrogate model to relate operating parameters to performance for guiding experiments, without demanding simulations, via active learning. Heat recirculation is demonstrated in this scaled reactor, further extending its scale and reducing the energy demand. An HMF yield of ~55 % and a productivity of 0.1 kg/hr, 8x higher than any other reactor, at a flowrate >20x than prior work, are demonstrated while maintaining energy efficiency of >98 %. A basic economic analysis estimates a $1.85/kg cost of HMF from fructose and >60 % reduction in CO 2 emissions than a conventional system without considering the impact of green electricity. We demonstrate that MW technology is well-applied to systems larger than the laboratory scale, and strategies to further reduction in cost and CO 2 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Closed impeller with self-recirculation casing treatment

A closed impeller includes: a housing provided with an opening; an insert inserted into one side of the opening; and an impeller body inserted into the other side of the opening, and including a plurality of blades and a shroud covering the plurality of blades, wherein the insert is provided with an inflow channel guiding a fluid into an inlet of the impeller body, and wherein a clearance between an outer surface of the insert and an inner surface of the housing forms a circulation flow channel connecting the other end of the inflow channel with a portion where one end of the inflow channel meets the inlet of the impeller body.

Pelton, Robert Young↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗