Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “burn model”

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 325 records · Page 18

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗

Observations of the Bright Star in the Globular Cluster 47 Tucanae (NGC 104)

The Bright Star in the globular cluster 47 Tucanae (NGC 104) is a post-asymptotic giant branch (post-AGB) star of spectral type B8 III. The ultraviolet spectra of late-B stars exhibit myriad absorption features, many due to species unobservable from the ground. The Bright Star thus represents a unique window into the chemistry of 47 Tuc. We have analyzed observations obtained with the Far Ultraviolet Spectroscopic Explorer, the Cosmic Origins Spectrograph aboard the Hubble Space Telescope, and the Magellan Inamori Kyocera Echelle Spectrograph on the Magellan Telescope. By fitting these data with synthetic spectra, we determine various stellar parameters (T {sub eff} = 10,850 ± 250 K, logg=2.20±0.13) and the photospheric abundances of 26 elements, including Ne, P, Cl, Ga, Pd, In, Sn, Hg, and Pb, which have not previously been published for this cluster. Abundances of intermediate-mass elements (Mg through Ga) generally scale with Fe, while the heaviest elements (Pd through Pb) have roughly solar abundances. Its low C/O ratio indicates that the star did not undergo third dredge-up and suggests that its heavy elements were made by a previous generation of stars. If so, this pattern should be present throughout the cluster, not just in this star. Stellar-evolution models suggest that the Bright Star is powered by a He-burning shell, having left the AGB during or immediately after a thermal pulse. Its mass (0.54 ± 0.16M {sub ⊙}) implies that single stars in 47 Tuc lose 0.1–0.2 M {sub ⊙} on the AGB, only slightly less than they lose on the red giant branch.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Impact of Resolution on Double-detonation Models for Type Ia Supernovae

Abstract Thermonuclear supernovae are the result of the violent unbinding of a white dwarf (WD), but the precise nature of the explosion mechanism(s) is a matter of active debate. To this end, several specific scenarios have been proposed to explain the observable traits of Type Ia supernovae. A promising pathway is the double-detonation scenario, where a WD accretes a shell of helium-rich material from a companion and a detonation in the resulting helium shell is the primary cause of the explosion. Through a set of two-dimensional grid-based simulations of this scenario we clearly distinguish three phases of evolution: external helium-rich detonation, core compressive heating, and a final core carbon burn. Though final disruption of the whole system is achieved at all resolutions, only models with minimum resolutions of 4 km and better exhibit all three phases. Particularly, core compression detonation is only observed for higher resolutions, producing qualitatively different nucleosynthetic outcomes. We identify the effect of finer spatial resolution on the mixing of hot silicon at the interface between the detonating helium layer and the underlying C/O WD as a primary driver of these dynamic differences.

79 ASTRONOMY AND ASTROPHYSICS↗

Biomass Burning Emissions of Black Carbon over the Maritime Continent and ENSO Variability

Fire emissions from the Maritime Continent (MC) over the western tropical Pacific are strongly influenced by El Niño–Southern Oscillation (ENSO), posing various climate effects to the Earth system. Here, in this study, we show that the historical biomass burning emissions of black carbon (BC bb ) aerosol in the dry season from the MC are strengthened in El Niño years due to the dry conditions. The eastern Pacific type of El Niño exerts a stronger modulation in BC bb emissions over the MC region than the central Pacific type of El Niño. Based on simulations using the fully coupled Community Earth System Model (CESM), the impacts of increased BC bb emissions on ENSO variability and frequency are also investigated in this study. With BC bb emissions from the MC scaled up by a factor of 10, which enables the identification of climate response from the internal variability, the increased BC bb heats the local atmosphere and changes land–sea thermal contrast, which suppresses the westward transport of the eastern Pacific surface water. It leads to an increase in sea surface temperature in the eastern tropical Pacific, which further enhances ENSO variability and increases the frequency of extreme El Niño and La Niña events. This study highlights the potential role of BC bb emissions on extreme ENSO frequency, and this role may be increasingly important in the warming future with higher wildfire risks.

54 ENVIRONMENTAL SCIENCES↗

Augmentation of WRF-Hydro to simulate overland-flow- and streamflow-generated debris flow susceptibility in burn scars

In steep wildfire-burned terrains, intense rainfall can produce large runoff that can trigger highly destructive debris flows. However, the ability to accurately characterize and forecast debris flow susceptibility in burned terrains using physics-based tools remains limited. Here, we augment the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) to simulate both overland and channelized flows and assess postfire debris flow susceptibility over a regional domain. We perform hindcast simulations using high-resolution weather-radar-derived precipitation and reanalysis data to drive non-burned baseline and burn scar sensitivity experiments. Our simulations focus on January 2021 when an atmospheric river triggered numerous debris flows within a wildfire burn scar in Big Sur – one of which destroyed California's famous Highway 1. Compared to the baseline, our burn scar simulation yields dramatic increases in total and peak discharge and shorter lags between rainfall onset and peak discharge, consistent with streamflow observations at nearby US Geological Survey (USGS) streamflow gage sites. For the 404 catchments located in the simulated burn scar area, median catchment-area-normalized peak discharge increases by ~ 450 % compared to the baseline. Catchments with anomalously high catchment-area-normalized peak discharge correspond well with post-event field-based and remotely sensed debris flow observations. We suggest that our regional postfire debris flow susceptibility analysis demonstrates WRF-Hydro as a compelling new physics-based tool whose utility could be further extended via coupling to sediment erosion and transport models and/or ensemble-based operational weather forecasts. Given the high-fidelity performance of our augmented version of WRF-Hydro, as well as its potential usage in probabilistic hazard forecasts, we argue for its continued development and application in postfire hydrologic and natural hazard assessments.

54 ENVIRONMENTAL SCIENCES↗

The Three-dimensional Collapse of a Rapidly Rotating 16 M ⊙ Star

I report on the three-dimensional (3D) hydrodynamic evolution of a rapidly rotating 16 M ⊙ star to iron core collapse. For the first time, I follow the 3D evolution of the angular momentum (AM) distribution in the iron core and convective shell burning regions for the final 10 minutes up to and including gravitational instability and core collapse. In 3D, convective regions show efficient AM transport that leads to an AM profile that differs in shape and magnitude from MESA within a few shell convective turnover timescales. For different progenitor models, such as those with tightly coupled Si/O convective shells, efficient AM transport in 3D simulations could lead to a significantly different AM distribution in the stellar interior affecting estimates of the natal neutron star or black hole spin. The results suggest that 3D AM transport in convective and rotating shell burning regions are critical components in models of massive stars and could qualitatively alter the explosion outcome and inferred compact remnant properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Soil carbon storage capacity of drylands under altered fire regimes

The determinants of fire-driven changes in soil organic carbon (SOC) across broad environmental gradients remains unclear, especially in global drylands. Here we combined datasets and field sampling of fire-manipulation experiments to evaluate where and why fire changes SOC and compared our statistical model to simulations from ecosystem models. Drier ecosystems experienced larger relative changes in SOC than humid ecosystems—in some cases exceeding losses from plant biomass pools—primarily explained by high fire-driven declines in tree biomass inputs in dry ecosystems. Many ecosystem models underestimated the SOC changes in drier ecosystems. Upscaling our statistical model predicted that soils in savannah–grassland regions may have gained 0.64 PgC due to net-declines in burned area over the past approximately two decades. Consequently, ongoing declines in fire frequencies have probably created an extensive carbon sink in the soils of global drylands that may have been underestimated by ecosystem models.

54 ENVIRONMENTAL SCIENCES↗

Post-fire time series of sensor and geochemistry sample data from surface water, groundwater, precipitation, soil, and vegetation across Oak Creek watershed, Washington

This dataset supports a broader study examining wildfire impacts on hydrologic connectivity across 5 sites within the Oak Creek watershed and the resulting biogeochemical impacts. Stream sites were selected using the Advanced Terrestrial Simulator (ATS) hydrologic model to identify locations with varying groundwater contributions and hydrologic responses across different burn severity scenarios. The Retreat Fire burned from July 23 to August 2 in 2024, affecting the five study sites at varying burn severities. Each site is equipped with YSI EXO2 sondes logging sub-hourly throughout the year, and grab samples are collected approximately every six weeks. YSI sondes are used to measure temporally resolved proxies for groundwater inputs (specific conductivity) and organic matter (fluorescent dissolved organic matter; fDOM) along with basic water quality and depth. Grab samples of surface water, groundwater, and precipitation are analyzed for water stable isotopes and conductivity to understand endmembers for hydrologic mixing Grab samples of surface water, groundwater, soil water, and litter/vegetation/soil leachates are analyzed for organic matter composition measured by Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) to understand organic matter dynamics. Game camera photos are provided in a separate data package available at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018598. Future versions of this dataset will include time series data from YSI EXO2 sondes (fDOM, dissolved oxygen, temperature, depth, specific conductance, turbidity, pH), BaroTROLL sensors (air temperature and barometric pressure), rain gauges (precipitation), and data from the soil and vegetation samples. Because this study is ongoing, this data package will be updated regularly to include newly collected data and the additional data types. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) a folder of field photos; (2) a folder of surface water sample data; (3) a folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data; (4) a data checks report; (5) file-level metadata; (6) data dictionary; (7) field metadata; (8) readme; (9) international generic sample number (IGSN) mapping file; and (10) field protocols. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) total dissolved nitrogen data and averages; (3) stable water isotopes and averages; (4) methods codes; (5) FTICR-MS methods; and (15) a subfolder of 9.4 Tesla (9.4T) FTICR-MS data. This folder contains the processed data and three subfolders, one containing the .xml files, one containing the water CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .d, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, .png, .mov, or .mp4.

Biogeochemistry↗

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗

Integrated operation scenarios: Chapter 6 of the special issue: on the path to tokamak burning plasma operation

Here we report the progress of the development and optimization of operational scenarios for ITER and beyond, focusing upon baseline, hybrid, and steady-state scenarios since 2007. This includes advancements made by the integrated operation scenarios (IOS) topical group of the international tokamak physical activity as well as contributions from the broader tokamak community. The key area of research involves developing IOSs that encompass tokamak physics, operation, and technology by utilizing integrated modeling and control strategies. This requires leveraging available actuators to simultaneously control plasma position and shape, MHD activities that could lead to disruptions, transport, plasma-wall interaction and power exhaust, fuel cycle, fusion burn, and tritium breeding. The control extends from the plasma initiation phase, through the current ramp-up, flattop, start and end of the fusion burn, and current ramp-down, to the plasma termination phase. A review of the currently developed scenarios and modeling is provided in terms of (i) optimizing plasma initiation in ITER, (ii) preparing for the low activation phase to fully commission all tokamak systems and establish and validate physics and scenario conditions in preparation for deuterim-tritium (DT) operation, (iii) developing and preparing baseline and hybrid scenarios to demonstrate the feasibility of achieving these regimes within device constraints, (iv) exploring steady-state scenarios to meet ITER’s steady-state goals, (v) evaluating and preparing actuators for ITER, (vi) developing integrated control solutions using shared actuators. The most notable achievements include; (i) the development of ITER demonstration discharges by matching various dimensionless parameters, (ii) the development of scenarios in an ITER-like tungsten environment and DT operation, and (iii) the development of scenarios in superconducting tokamaks, enabling long-pulse operations with similar coil constraints to ITER. Along with these significant achievements, outstanding issues and recommendations for further research and development are provided. Importantly, this study goes beyond simply updating the ITER Physics Basis; it carries profound implications for the broader field of burning plasma research, offering valuable insights and guidance for the next generation of fusion experiments and devices.

ITER↗

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

Data reduction considerations for spherical R-32(CH 2 F 2 )-air flame experiments

The burning velocities of mixtures of refrigerant R-32 (CH 2 F 2 ) with air over a range of equivalence ratios are studied via shadowgraph images of spherically expanding flames (SEFs) in a large, optically accessible spherical chamber at constant pressure. Numerical simulations of the 1-D, unsteady, spherical flames incorporating an optically thin radiation model and detailed kinetics accurately predict experimental results for a range of equivalence ratios. For these low burning velocity flames, the effects of stretch and radiation occur simultaneously and make extraction of the unstretched burning velocity from the experimental data difficult. Different data reduction approaches are shown to have large effects on the burning velocities inferred from the experiments and simulations. A new flame radius tracking approach for the experimental images is shown to provide improved agreement of the burned gas velocity variation with stretch predicted by the simulations and helps to compensate for mild flame distortion due to buoyancy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the impact of smoke from prescribed fire on road visibility

Prescribed fires are planned to achieve conservation and fuel reduction objectives while minimizing smoke ground concentration to limit health impacts and road visibility impairment. Prescribed burns cannot indeed be conducted if those hazards are not within predefined limits. This paper proposes a new framework to evaluate road visibility that overcomes the limitation of the state of the art model, VSMOKE. The framework leverages the fast-running framework QUIC-Fire/QUIC-SMOKE to capture fire and smoke dynamics, and the timing and duration of hazardous conditions on the road network close to the burn unit (within 50–100 km). The paper presents parametric study using a real burn plot at Fort Stewart (GA, USA), under hypothetical wind conditions to understand the interplay between buoyancy and smoke dilution. Here, results showed that faster winds caused fire escape while slower winds did not achieve a complete burn. Furthermore, faster winds featured brief road visibility reduction below braking distance.

54 ENVIRONMENTAL SCIENCES↗

A statistical model to predict ignition probability

Ignition capability is a critical design constraint for aeronautical gas turbines. However the current trend toward overall lean burn is detrimental to the engine ignition and relight and the ignition system must be adapted to ensure a fast and reliable light-round in all circumstances. As ignition is a stochastic phenomenon, the optimization of an ignition system requires to build ignition probability maps, which is difficult and costly with either experiment or numerical simulation as both require many tests. Here, this work proposes a model to predict the ignition probability map, knowing only flow statistics in non-reacting conditions, i.e., with only one test. The originality of the model is to construct statistics of the flame kernel trajectory, which are then combined with local flow indicators to evaluate the ignition probability at the considered sparking location. Application to a swirled burner operated in premixed, non-premixed and spray combustion modes illustrates the model concepts and demonstrates its ability to recover the experimental ignition map with good accuracy.

42 ENGINEERING↗

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time resolved ablator areal density during peak fusion burn on inertial confinement fusion implosions

Near peak compression, inertial confinement fusion implosions release both deuterium–tritium (DT) fusion gamma rays and neutron induced gamma rays from carbon from the areal density of the remaining ablator shell. The gamma reaction history diagnostic makes a time resolved measurement of both. Across many recent implosions, the carbon gamma ray peak arrives systematically 11 ± 10 ps later compared to DT fusion burn. The timing shift is consistent with the carbon areal density increasing throughout the peak of the fusion burn, implying that the carbon portion of the capsule continues to converge. A model finds that the observed timing shift is consistent with a 4π averaged carbon ablator inward velocity of 80 μm/ns for the contemporary National Ignition Facility implosions. In this work, the timing shift is possibly related to the energy balance of the implosion, with the expectation that a high performing, igniting capsule would see the carbon gamma rays arrive before the DT fusion peak.

Nuclear fusion↗

Investigating the Origins of Cyclic Variability in Internal Combustion Engines Using Wall-Resolved Large Eddy Simulations

Modern internal combustion engines (ICE) operate at the ragged edge of stable operation characterized by high cycle-to-cycle variations (CCV). A key scientific challenge for ICE is the understanding, modeling, and control of CCV in engine performance, which can contribute to partial burns, misfire, and knock. The objective of this study is to use high-fidelity numerical simulations to improve the understanding of the causes of CCV. Nek5000, a leading high-order spectral element, open source code, is used to simulate the turbulent flow in the engine combustion chamber. Multicycle, wall-resolved large-eddy simulations (LESs) are performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. The mean and root-mean-square (rms) of the in-cylinder flow fields at various piston positions are validated using particle image velocimetry (PIV) measurements during the intake and compression strokes. The large-scale flow structures, including the swirl and tumble flow patterns, are analyzed in detail and the causes for cyclic variabilities in these flow features are explained. The energy distribution across the different scales of the flow are quantified using one-dimensional (1D) energy spectra, and the effect of the tumble breakdown process on the energy distribution is examined. Finally, the insights from this study can help us develop improved engine designs with reduced cyclic variabilities in the in-cylinder flow leading to enhanced engine performance.

33 ADVANCED PROPULSION SYSTEMS↗