Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence optimization”

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 559 records · Page 31

An overview of 3D field optimization for control of transport and edge instabilities on KSTAR

An international team from several laboratories and universities has made key advances over the last few years in the control of plasma transport and edge instabilities with applied 3D fields in the KSTAR tokamak to optimize long pulse operation scenarios. This overview begins with the optimization of both core and edge resonant magnetic perturbations (RMP) to improve fast ion confinement to avoid excessive limiter heat loads due to fast ion losses and successful modeling of the experimental results. Integrated and advanced plasma control techniques with machine learning (ML) and adaptive control were then used to optimize the 3D field spectrum in real-time to control edge localized modes (ELMs) while avoiding core locked modes that could disrupt the plasma. Accelerating the offline model of 3D fields with a surrogate ML model can optimize ELM suppression in the edge while limiting the impact of the applied RMP fields deeper in the plasma core in real-time. In addition, the impact of the 3D fields on the divertor heat load has been modeled and compared with experimental measurements. An analysis of a multi-machine database including KSTAR has been performed to better understand the metrics for the observed RMP thresholds for ELM suppression and the resulting plasma performance. Predictive modeling of the operational space for ELM suppression and density pumpout due to RMP has shown the importance of magnetic islands in the plasma edge and their impact on plasma turbulence. This research has culminated in the development of successful long pulse operational scenarios on KSTAR while attempting to overcome challenges of the new tungsten divertor.

3D fields↗

Multi-Scale Simulation of Wind Farm Performance during a Frontal Passage

Predicting the response of wind farms to changing flow conditions is necessary for optimal design and operation. Here in this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine–flow interaction, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with 10 m grid spacing. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed background velocity, turbine wake structure, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD model framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

Stability of boundary layers with porous suction strips: Experiment and theory

Low turbulence tunnel experiments on the stability and transition of 2 D boundary layers on flat plates with and without suction are described. A number of general suction cases are discussed. Test results showed that the maximum stabilization occurred when the suction was moved toward the Branch I neutral point. An analytical study of the stability of two dimensional, incompressible boundary layer flows over plates with suction through porous strips was performed. The mean flow was calculated using linearized triple deck, closed form solutions. The stability results of the triple deck theory are shown to be in good agreement with those of the interacting boundary layers. An analytical optimization scheme for the suction configuration was developd. Numerical calculations were performed corresponding to the experimental configurations. In each case, the theory correctly predicts the experimental results.

Reynolds, G. A.↗

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

Validating and optimizing mismatch tolerance of Doppler backscattering measurements with the beam model (invited)

In this work, we use the beam model of Doppler backscattering (DBS), which was previously derived from beam tracing and the reciprocity theorem, to shed light on mismatch attenuation. This attenuation of the backscattered signal occurs when the wavevector of the probe beam’s electric field is not in the plane perpendicular to the magnetic field. Correcting for this effect is important for determining the amplitude of the actual density fluctuations. Previous preliminary comparisons between the model and Mega-Ampere Spherical Tokamak (MAST) plasmas were promising. In this work, we quantitatively account for this effect on DIII-D, a conventional tokamak. We compare the predicted and measured mismatch attenuation in various DIII-D, MAST, and MAST-U plasmas, showing that the beam model is applicable in a wide variety of situations. Finally, we performed a preliminary parameter sweep and found that the mismatch tolerance can be improved by optimizing the probe beam’s width and curvature at launch. This is potentially a design consideration for new DBS systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Numerical Simulations of High-Speed Turbulent Boundary Layers over Riblets

Direct numerical simulations (DNS) of spatially developing turbulent boundary layers over riblets with a broad range of riblet spacings are conducted to investigate the effects of riblets on skin friction at high speeds. Zero-pressure gradient boundary layers under two flow conditions (Mach 2:5 with T(sub w)/T(sub r) = 1 and Mach 7:2 with T(sub w)/T(sub r) = 0:5) are considered. The DNS results show that the drag-reduction curve (delta C(sub f)/C(sub f) vs l(sup +)(sub g )) at both supersonic speeds follows the trend of low-speed data and consists of a `viscous' regime for small riblet size, a `breakdown' regime with optimal drag reduction, and a `drag-increasing' regime for larger riblet sizes. At l l(sup +)(sub g) approx. 10 (corresponding to s+ approx 20 for the current triangular riblets), drag reduction of approximately 7% is achieved at both Mach numbers, and con rms the observations of the few existing experiments under supersonic conditions. The Mach- number dependence of the drag-reduction curve occurs for riblet sizes that are larger than the optimal size, with smaller slopes of (delta C(sub f)/C(sub f) for larger freestream Mach numbers. The Reynolds analogy holds with 2(C(sub h)=C(sub f) approximately equal to that of at plates for both drag-reducing and drag-increasing configurations.

Duan, Lian↗

A review of some approximate methods used in aerodynamic heating analyses

It is pointed out that preliminary design and optimization studies for new aerospace vehicles require techniques which can calculate aerodynamic heating rates accurately and efficiently. The method employed to calculate the flow field depends to a large extent on the shape of the vehicle, Mach number, Reynolds number, and Knudsen number. In the case of the aero-assisted orbital transfer vehicle (AOTV), a substantial portion of the flight will be in the transitional regime between continuum and free molecule flow. The present paper discusses some approximate methods which have been used to calculate heating rates on high-speed vehicles. Attention is given to the stagnation point and leading edges, the downstream region, the axisymmetric analog, laminar and turbulent heating rates, transition heating rates, gas models, and three-dimensional applications.

Dejarnette, F. R.↗

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

Magnetohydrodynamics with physics informed neural operators

Abstract The modeling of multi-scale and multi-physics complex systems typically involves the use of scientific software that can optimally leverage extreme scale computing. Despite major developments in recent years, these simulations continue to be computationally intensive and time consuming. Here we explore the use of AI to accelerate the modeling of complex systems at a fraction of the computational cost of classical methods, and present the first application of physics informed neural operators (NOs) (PINOs) to model 2D incompressible magnetohydrodynamics (MHD) simulations. Our AI models incorporate tensor Fourier NOs as their backbone, which we implemented with the TensorLY package. Our results indicate that PINOs can accurately capture the physics of MHD simulations that describe laminar flows with Reynolds numbers Re ⩽ 250 . We also explore the applicability of our AI surrogates for turbulent flows, and discuss a variety of methodologies that may be incorporated in future work to create AI models that provide a computationally efficient and high fidelity description of MHD simulations for a broad range of Reynolds numbers. The scientific software developed in this project is released with this manuscript.

97 MATHEMATICS AND COMPUTING↗

Multiple Mode Actuation of a Turbulent Jet

The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Pack, LaTunia G.↗

Multiple Mode Actuation of a Turbulent Jet

The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Pack, LaTunia G.↗

Truncated ARQ Statistical Link Analysis for Dynamic Links

The future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non Gaussian and non linear effects such as atmospheric turbulence, scintillation, antenna mis-pointing, jitter, etc. These dynamic links thus will experience various degrees of fading loss, and some of these link disruptions cannot be effectively mitigated by forward error correction coding and/or interleaving. One effective way to ensure reliable communication is by using Automatic Repeat Request (ARQ) protocol, where the receiver acknowledges to the transmitter whether or not a data unit is successfully received. If a data unit is not successfully received (such as after a pre-set time-out), the transmitter would then re-transmit the lost data unit to the receiver. In a previous paper, we derived a statistical link analysis method of finding the optimal operating Signal-to-Noise Ratio (SNR) and estimating the latency of an ARQ scheme. In a more recent paper, we demonstrated the above method using the SNR distribution constructed from the Ka-band (32 GHz) flight data. To simplify the discussion, we considered the academic approach that the ARQ scheme allows for an infinite number of retransmissions. In this paper, we consider the more practical case of a truncated ARQ scheme, where there is a limit on the number of retransmissions. We derive the error probability, the optimal SNR setting, and the latency statistics of the correctly received frames of the truncated ARQ schemes. We first discuss the truncated ARQ link analysis principles using the Gaussian assumption for SNR distribution with a large variance. Next, we demonstrate the statistical truncated ARQ link analysis using the SNR distribution constructed from the Ka-band flight data. The results in this paper can be applied in the design of reliable communication systems such as the Consultative Committee for Space Data System (CCSDS) File Transfer Protocol (CFTP) and the Delay Tolerant Network (DTN).

Morabito, David↗

Global Estimates of PBL Depth from Space-Borne LIDAR

The planetary boundary layer (PBL) is responsible for communicating the exchange of energy, moisture, momentum, pollutants, and aerosols between the surface and the free atmosphere and is therefore crucial to many studies of the atmosphere. Unfortunately, there have historically been few observations of this important layer due to the complexity involved in its measurement. However, with the advent of more advanced satellites, global measurements of the PBL are now becoming possible. The PBL is often characterized by a high concentration of aerosols within the layer and low level clouds capping it and these are observable from space. The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite was launched in 2006 with the intention of observing aerosols and clouds and was the first space-based LIDAR optimized for this purpose. CALIPSO observations are therefore well suited to observing the depth of the PBL. Since it was launched, CALIPSO has been making nearly continuous measurements enabling a global picture of PBL depth. We plan to present a global PBL depth product and how it evolves throughout the year. The product is able to identify deeper PBL depths in the summer hemisphere over land and deeper depths along the northern hemisphere oceanic storm tracks in winter associated with cold air traveling over warm water. Large seasonal cycles are also evident in the subtropical desert locations among other features. In addition, comparisons will be made between several estimates of PBL depth based on turbulent intensity, meteorology profiles, and aerosol profiles from the GEOS5 model.

McGrath-Spangler, Erica lynn↗

The Results of a Simulator Study to Determine the Effects on Pilot Performance of Two Different Motion Cueing Algorithms and Various Delays, Compensated and Uncompensated

A study was conducted employing the Visual Motion Simulator (VMS) at the NASA Langley Research Center, Hampton, Virginia. This study compared two motion cueing algorithms, the NASA adaptive algorithm and a new optimal control based algorithm. Also, the study included the effects of transport delays and the compensation thereof. The delay compensation algorithm employed is one developed by Richard McFarland at NASA Ames Research Center. This paper reports on the analyses of the results of analyzing the experimental data collected from preliminary simulation tests. This series of tests was conducted to evaluate the protocols and the methodology of data analysis in preparation for more comprehensive tests which will be conducted during the spring of 2003. Therefore only three pilots were used. Nevertheless some useful results were obtained. The experimental conditions involved three maneuvers; a straight-in approach with a rotating wind vector, an offset approach with turbulence and gust, and a takeoff with and without an engine failure shortly after liftoff. For each of the maneuvers the two motion conditions were combined with four delay conditions (0, 50, 100 & 200ms), with and without compensation.

Guo, Li-Wen↗

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗

Computational Study of Flow Interactions in Coaxial Rotors

Although the first idea of coaxial rotors appeared more than 150 years ago, most helicopters have used single main-rotor/tail-rotor combination. Since reactive moments of coaxial rotors are canceled by contra-rotation, no tail rotor is required to counter the torque generated by the main rotor. Unlike the single main rotor design that distributes power to both main and tail rotors, all of the power for coaxial rotors is used for vertical thrust. Thus, no power is wasted for anti-torque or directional control. The saved power helps coaxial rotors reach a higher hover ceiling than single rotor helicopters. Another advantage of coaxial rotors is that the overall rotor diameter can be reduced for a given vehicle gross weight because each rotor provides a maximum contribution to vertical thrust to overcome vehicle weight. However, increased mechanical complexity of the hub has been one of the challenges for manufacturing coaxial rotorcraft. Only the Kamov Design Bureau of Russia had been notably successful in production of coaxial helicopters until Sikorsky built X2, an experimental compound helicopter. Recent developments in unmanned aircraft systems and high-speed rotorcraft have renewed interest in the coaxial configuration. Multi-rotors are frequently used for small electric unmanned rotorcraft partly due to mechanical simplicity. The use of multiple motors provides redundancy as well as cost-efficiency. The multi-rotor concept has rarely been used until recently because of its inherent stability and control problems. However, advances in inexpensive electronic flight control systems have opened the floodgates for small drones using multirotors. Coaxial rotors have started to appear in some multi-rotor configurations. Small coaxial rotors have often been designed using a hundred year old approach that is "sketch, build, fly, and iterate." In that approach, there is no systematic way to explore trade-offs or determine logical next steps. It is neither possible to account for multiple real-world constraints up front in design nor possible to know what performance is possible with a given design. Since unmanned vehicles are sized and optimized for the particular mission, a modern low-fidelity conceptual design and sizing tool that has been used for the design of large helicopters can be used for design of small coaxial rotorcraft. However, unlike most helicopters with single main rotor, the interactions between the upper and lower rotors emerge as an important factor to consider in design because an increase in performance of a multi-rotor system is not proportional to the number of rotors. Interference losses and differences in thrusts between the upper and lower rotors were investigated by theoretical methods as well as a computational fluid dynamics (CFD) method using the Reynolds-Averaged Navier-Stokes (RANS) equations. In this work, hybrid turbulence models are used to investigate the physics of interactions between coaxial rotors and a fuselage that are not well understood. Present study covers not only small-scale drones but also large-scale coaxial rotors for heavy-lifting missions. Considering the recently proposed FAA drone rules that require the flight only in visual line-of-sight, a large multirotor might be used as an airborne carrier for launch and recovery of unmanned aircraft systems with a human operator onboard. For applications to civil operations, their aerodynamic performance and noise levels need to be assessed. Noise is one of the largest limiting factors to rotorcraft operations in urban area. Since the high-frequency noise of multi-rotors may increase the annoyance, noise may turn out to be a key issue that must be addressed for market acceptability. One of the objectives of the present work is to study the effects of inter-rotor spacing and collectives on the performance, efficiency, and acoustics of coaxial rotor systems.

Study↗

Performance analysis of a noncoherently combined large aperture optical heterodyne receiver

The performance of a noncoherently combined (NC) multiple-mirror heterodyne receiver is analyzed. With no atmospheric turbulence (AT), the performance of the NC receiver is shown to the inferior to that of a monolithic diffraction-limited receiver with equivalent aperture area. When AT is taken into consideration, however, the efficiency of a monolithic receiver is limited by the phase coherence length of the atmosphere, and generally does not improve with increasing aperture size. In contrast, the performance of an NC system improves with an increasing number of receivers. Consequently, given a fixed collecting area, the NC system can offer superior performance. The performance of the NC heterodyne receiver is studied by analyzing the combining loss of the receiver SNR. It is shown that, given a constant collecting area, the performance of the NC receiver is optimized when the diameter of each of the individual receivers is on the order of the phase coherence length of the AT.

Chen, Chien-Chung↗

Plasma-Assisted Pre-Chamber Ignition System for Highly Dilute Stoichiometric Heavy-Duty Natural Gas Engines (Final Technical Report)

This project explored advanced ignition technologies to significantly enhance efficiency and reduce operating costs for heavy-duty natural gas engines operating at stoichiometric conditions, while meeting ultra-low NOx emission standards. The main goal was to develop and validate a plasma-assisted pre-chamber ignition system that could deliver at least a 2% increase in brake thermal efficiency (BTE) and a 4% decrease in total cost of ownership (TCO) compared to a typical multi-cylinder engine with three-way catalyst aftertreatment, ensuring compatibility with the expected 2027 EPA/CARB regulations. In the first half of the project, the research team concentrated on developing and testing plasma-assisted pre-chamber ignition using nanosecond pulsed discharges. Extensive experiments were conducted in an optically accessible rapid-compression and expansion machine, a constant-volume chamber, and an optical single-cylinder engine. Experiments were coupled with CFD simulations. The work produced unique insights into pre-chamber flame formation, jet ignition, dilution effects, and flame quenching at pressures, temperatures, and dilution levels relevant to engines. Although plasma-assisted ignition showed promise in controlled lab settings, the research also identified fundamental and practical challenges when applying this technology to real engine conditions. Midway through the project, a crucial pivot was made, guided by three key findings. First, the power electronics required for nanosecond plasma discharges were found to be too costly for commercial use, undermining the project’s cost-of-ownership goals. Second, nanosecond plasma ignition was highly sensitive to turbulent flow in the pre-chamber, resulting in lower ignition reliability than traditional spark under engine-like conditions. Third, achieving a truly diffuse low-temperature plasma at high pressures near top dead center was not possible, reducing the anticipated chemical enhancement benefits. These results collectively suggested that continuing with plasma-assisted ignition was unlikely to meet both efficiency and cost objectives. In response, the project shifted focus to a more realistic approach: enhancing traditional spark-based pre-chamber ignition with significantly less spark energy. Using insights gained earlier in the project, the team redesigned the pre-chamber to maintain high dilution tolerance and quick combustion, even with lower ignition energy. Testing confirmed that with optimized pre-chamber design and combustion timing, a lower-energy spark could reliably ignite highly diluted stoichiometric mixtures, reduce burn time, and boost thermal efficiency. Final engine testing and techno-economic analysis verified that this revised approach successfully achieved the project goals. The optimized pre-chamber ignition system provided over a 2% increase in calculated brake thermal efficiency compared to the baseline engine. Notably, the lower ignition energy and simplified hardware reduced component stress, extended maintenance intervals, and lowered the total cost of ownership. When used with stoichiometric operation and traditional three-way aftertreatment, the system remained compatible with near-zero NOx emissions targets without increasing cost or complexity in the emissions control system. In summary, although the project deviated from its initial plasma-assisted ignition idea, the work produced a more practical and commercially viable solution. The results show that precisely optimized, low-energy pre-chamber spark ignition can significantly improve efficiency and reduce overall ownership costs for heavy-duty natural gas engines. This directly aligns with DOE goals for cleaner, more efficient, and cost-effective transportation technologies.

03 NATURAL GAS↗