Search NASA⌕ Search

SEARCH · Search NASA

Results for “varying data rates”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

SAS4A/SASSYS-1 Validation with the FFTF LOFWOS Tests #10-12

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. The results of modeling and simulation of LOFWOS Tests #10, #11, and #12 are presented in this report to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. These three tests did not use the primary loop pony motors and thus transitioned to natural circulation flow rates from varying initial power levels. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project (CRP) for benchmark analysis of the FFTF LOFWOS Test #13. The SAS4A/SASSYS-1 model originally developed for the CRP has been extended to include LOFWOS Tests #10-12 and modernized to utilize new code features and modeling practices. Predicted results from SAS4A/SASSYS-1 were compared to measured test data and differences were investigated. Overall, it was concluded that, for LOFWOS Tests #10-12, there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Thermodynamic control on the decomposition of organic matter across different electron acceptors

The increasing availability of high-resolution characterization of natural organic matter (OM) data has shifted the paradigm of lumped descriptions of OM components and potential microbial activities. Our recent development of a substrate-explicit thermodynamic model uniquely enables incorporating complex OM pools to formulate biogeochemical reaction models based on their elemental compositions. While this previous work facilitates prediction of aerobic respiration of complex OM, it is equally imperative to consider the role of non-oxygenic electron acceptors in regulating OM turnover and the fate of carbon. In this study, we significantly expand our previous model by flexibly incorporating both detailed OM chemistry and electron acceptors other than oxygen. Here, our modeling analysis has revealed substantial variations in the energy status of OM molecules across different soils, which drive the co-occurrence of different electron-accepting processes. We demonstrated the effectiveness of the proposed model using a consistency check with experimental data. Through systematic evaluation of the impact of diverse chemical inputs (both electron donors and acceptors) on OM decomposition, the new model also revealed how key microbial growth parameters such as carbon use efficiency (CUE) and reaction rates vary across different electron-accepting processes. Our model provides a unified framework integrating thermodynamic and kinetic constraints on microbial metabolic activities. It complements traditional kinetic models, which are often designed solely to capture mass fluxes. We conclude that thermodynamic modeling emerges as a powerful tool for describing the mechanisms underlying the interplay between microbial growth and OM chemistry and cycling across different electron acceptors, enhancing our ability to project complex ecosystem behaviors in dynamic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Frequency–Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress

Understanding the nucleation phase of earthquakes is key for advancing earthquake hazard assessment and improving earthquake early warning systems. However, little progress has been made in this area due to a poor understanding of nucleation processes and incomplete seismic and fault zone measurements. The ability to integrate measured fault zone properties with seismic data could significantly improve our understanding of how earthquakes begin and whether there are systematic variations in seismic properties preceding failure. In this work, we use high-resolution laboratory measurements of fault zone properties along with acoustic emission data to document temporal variations of foreshock properties. Our data show that foreshock size increases with shear stress, loading rate, and fault slip rate. We propose that the preseismic fault slip rate and fault zone thickness (i.e., porosity) work in concert to modulate foreshock properties.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2020

Degradation of Commercial Lithium-Ion Cells as a Function of Chemistry and Cycling Conditions

Energy storage systems with Li-ion batteries are increasingly deployed to maintain a robust and resilient grid and facilitate the integration of renewable energy resources. However, appropriate selection of cells for different applications is difficult due to limited public data comparing the most commonly used off-the-shelf Li-ion chemistries under the same operating conditions. This article details a multi-year cycling study of commercial LiFePO 4 (LFP), LiNi x Co y Al 1–x–y O 2 (NCA), and LiNi x Mn y Co 1–x–y O 2 (NMC) cells, varying the discharge rate, depth of discharge (DOD), and environment temperature. The capacity and discharge energy retention, as well as the round-trip efficiency, were compared. Even when operated within manufacturer specifications, the range of cycling conditions had a profound effect on cell degradation, with time to reach 80% capacity varying by thousands of hours and cycle counts among cells of each chemistry. The degradation of cells in this study was compared to that of similar cells in previous studies to identify universal trends and to provide a standard deviation for performance. All cycling files have been made publicly available at batteryarchive.org, a recently developed repository for visualization and comparison of battery data, to facilitate future experimental and modeling efforts.

25 ENERGY STORAGE↗
DOE OSTI2024

Optical-pump–terahertz-probe spectroscopy in high magnetic fields with kHz single-shot detection

We demonstrate optical pump–THz probe (OPTP) spectroscopy with a variable external magnetic field (0–9 T), in which the time-dependent THz signal is measured by echelon-based single-shot detection at a repetition rate of 1 kHz. The method reduces data acquisition times by more than an order of magnitude compared to conventional electro-optic sampling using a scanning delay stage. The approach illustrates the wide applicability of the single-shot measurement approach to non-equilibrium systems that are studied through OPTP spectroscopy, especially in cases where parameters such as magnetic field strength (B) or other experimental parameters are varied. We demonstrate the capabilities of our measurement by performing cyclotron resonance experiments in bulk silicon, where we observe B-field-dependent carrier relaxation and distinct relaxation rates for different carrier types. We use a pair of economical linear array detectors to measure 500 time points on each shot, offering an equivalent performance to camera-based detection with possibilities for higher repetition rates.

Instruments & Instrumentation↗
DOE OSTI2024

Ensemble modeling of global lake evaporation under climate change

Global projections of lake evaporation are typically based on simulations using single mechanistic models. However, because of its complex interactions with various lake physical properties, environmental and anthropogenic drivers, lake evaporation is highly variable and sensitive to the choice of model used. In this study, we present a multi-model analysis to investigate differences across global simulations of lake evaporation during the warm-season using three different lake models driven by outputs from four general circulation models (GCM) (i.e. 12 model combinations in total) for historic and future scenarios. Our results suggest substantial differences among lake-climate model simulations of lake evaporation. These differences varied throughout the 20th and 21st century, with model driver data explaining 74% of the variance in future projections of warm- season lake evaporation. Our projections indicate that, by the end of the 21st century (2070–2099), global annual lake evaporation rates will increase by 10–27% under Representative Concentration Pathways (RCPs) 2.6–8.5. Finally, we highlight the importance of using a multi-model approach for the prediction of future global lake evaporation responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Role of self-generated magnetic fields in the inertial fusion ignition threshold

Magnetic fields spontaneously grow at unstable interfaces around hot-spot asymmetries during inertial confinement fusion implosions. Although difficult to measure, theoretical considerations and numerical simulations predict field strengths exceeding 5 kT in current National Ignition Facility experiments. Magnetic confinement of electrons then reduces the rate of hot-spot heat loss by >5%. We demonstrate this via magnetic post-processing of two-dimensional xRAGE hydrodynamic simulation data at bang time. We then derive a model for the self-magnetization, finding that it varies with the square of the hot-spot temperature and inversely with the areal density. The self-magnetized Lawson analysis then gives a slightly reduced ignition threshold. Time-dependent hot-spot energy balance models corroborate this finding, with the magnetic field quadrupling the fusion yield for near-threshold parameters. The inclusion of magnetized multi-dimensional fluid instabilities could further alter the ignition threshold and will be the subject of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental Characterization of High-Surface Area Thermal Energy Storage

There is growing interest in energy storage technologies due to the expansion of renewable energy sources that are inherently intermittent and the increasing frequency of extreme weather events that disturb the power grid. Power consumption in buildings makes up approximately 76% of all electricity usage on the grid and is primarily used for thermal applications such as space conditioning, hot water, and cooking. This makes thermal energy storage (TES) an ideal solution for many of these applications. Many TES technologies rely on latent energy storage, which utilizes the melting/solidification of phase change materials (PCM) to store energy. Typically, TES designs suffer from low power density due to their low inherent thermal conductivity. This limitation makes the deployment of TES in active applications difficult as the ease of access to energy is essential for effective use. Common routes for improving power density include high thermal conductivity additives or extended features such as fins that increase cost. This study presents an alternative approach to improving performance through increasing the overall surface area to volume ratio of the device, to increase the available area for convection to occur between the working fluid and PCM. In the study a commercial PCM was selected with a transition temperature ideal for space heating applications. The heat exchanger design utilizes a unique application of triply periodic minimal surfaces (TPMS) for macro-encapsulation of the PCM. The use of TPMS for heat exchangers has been growing in interest due to their high-surface area to volume ratios, which were previously unmanufacturable until the development of additive manufacturing. A modular system was designed and manufactured with a stereolithography resin printing system that is then backfilled with PCM. An experiment test set-up is designed to test the charge and discharge performance of the thermal storage using a conditioned air stream. The pressure drop of the design is tested across a variety of flow rates. When compared to existing experimental data within literature, there is excellent agreement based on the Reynolds number at similar hydraulic diameters. Several inlet temperatures are tested at consistent temperature differences from the phase change temperature for both charging and discharging. Additionally, the volumetric flow rate is varied for each temperature set point. It was found that increasing flow rate had diminishing returns in reducing the overall charge time of the TES. The temperature delta from the melting temperature was the primary contributor to the change in average heat flux with limited variation in average heat transfer rate between charging and discharging at similar inlet temperatures and flow rates. The TPMS heat exchanger design has a high air-side pressure drop but it provides high heat transfer rates. This helps maintain a high outlet temperature during discharge, which is important to thermal comfort applications. The design, manufacturing, and experimental characterization of the TES device will be presented as part of this study.

25 ENERGY STORAGE↗
DOE OSTI2026

Rasterization with Data-Parallel Primitives

Parallel rasterization can suffer from race conditions during fragment generation, which is traditionally addressed by using specialized hardware accessible via vendor graphics APIs. Unfortunately, graphics APIs are increasingly problematic on high-performance computers, either because they are not provided or because of concerns about dependencies with in situ visualization. In response, we present a hardware-agnostic rasterization algorithm that handles race conditions using only data-parallel primitives (DPPs), enabling efficient rendering on HPC systems without graphics API dependencies and aligning with recent efforts to deliver visualization software with DPPs. Our evaluation consists of three phases: (1) evaluating portability across different CPU and GPU architectures, (2) evaluating competitiveness with a community standard, and (3) evaluating performance across varying workloads and available parallelism. The supporting experiments run on both AMD and NVIDIA GPUs, considering data sets as large as 460 million triangles and 160 million pixels. While performance generally falls short of graphics API baselines, it achieves interactive frame rates on most workloads. As a result, we conclude our approach is a viable solution for rasterization on high-performance computers since our approach is portably performant across different architectures without the need for specialized vendor support.

Buckley, Makani [University of Oregon] (ORCID:0009↗
DOE OSTI2020

Iterative and Self-Consistent Estimation of Degradation and Soiling Loss in PV Systems - A Case Study

In many parts of the world, soiling is a significant loss mechanism for PV energy production. Methods for accurately quantifying historical and current soiling losses are needed. Furthermore, state-of-the-art methods to estimate degradation rates based on production data do not handle soiling in a satisfactory manner. In order to accurately estimate either of these quantities, it is in many cases helpful to consider their combined effect on PV system performance. In this work we make use of the Combined Degradation and Soiling (CODS) algorithm, which iteratively decomposes time series of daily normalized energy yield into four components: a soiling ratio, a degradation trend, seasonal variations, and noise. We demonstrate the use of the model on field performance data from a utility-scale power plant in a high-soiling environment in the Middle East. We apply the algorithm to the temperature-corrected performance index of the DC power on inverter level, quantifying how soiling varies over time and through the plant. The irradiation-weighted average soiling loss in this plant is 2.3 %, with a median soiling rate of -0.086 %/day. Between the inverters in the site in study, estimated degradation rates range between -0.86 %/year and -0.15 %/year, with a mean rate of -0.48 %/year.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Population Dynamics and Temporal Trends of Bull Trout in the East Fork Salmon River, Idaho

Abstract Because of their long-term listing under the Endangered Species Act, much interest has been placed on estimating population vital rates for Bull Trout Salvelinus confluentus, but the biotic and abiotic factors that influence the interannual variability in those vital rates have rarely been evaluated. We used mark–recapture data to estimate fish growth, survival, and trends in abundance for fluvial adult Bull Trout in the East Fork Salmon River, Idaho. Over an 8-year period, a total of 1,205 individual Bull Trout were collected at a weir on the East Fork Salmon River (29 km upstream of its confluence with the Salmon River) during June–September, of which 420 were recaptures from prior years. Bull Trout varied in length from 215 to 756 mm and achieved a slightly larger asymptotic length and a slightly lower rate of growth relative to other fluvial and adfluvial Bull Trout populations. Apparent survival averaged 0.42 across all years, which was similar to previous studies estimating apparent survival for Bull Trout. The number of emigrating anadromous salmonid smolts in the upper Salmon River basin positively influenced East Fork Salmon River Bull Trout growth and survival, and survival was higher in years with lower annual discharge. Assessment of population growth via linear regression analysis indicated that the Bull Trout population was increasing during the study period (λ = 1.08; 95% confidence interval = 1.03–1.14). Our findings highlight the ecological link between the abundance of wild and hatchery Chinook Salmon Oncorhynchus tshawytscha, Sockeye Salmon O. nerka, and steelhead O. mykiss smolts and growth and survival of adult Bull Trout in systems where these species occur in sympatry.

Roth, Curtis J.↗
DOE OSTI2025

Temporal Variability in Reservoir Surface Area Is an Important Source of Uncertainty in GHG Emission Estimates

Ebullitive methane (CH 4 ) emissions in lentic ecosystems tend to concentrate at river-lake interfaces and within shallow littoral zones. However, inconsistent definitions of the littoral zone and static representations of the lake or reservoir surface area contribute to major uncertainties in greenhouse gas (GHG) emissions estimates, particularly in reservoirs with large water-level fluctuations. This study examines temporal variation in littoral and total surface areas of US reservoirs and demonstrates how different methods and data sources lead to discrepencies in reservoir GHG emissions at large scales and over time. We also explore variability in remotely sensed water occurrence according to maximum surface area, reservoir purposes, and hydrologic regions. Notably, the largest relative variability in surface area is exhibited by small reservoirs with a maximum surface area <1 km 2 and non-hydroelectric reservoirs. Additionally, we use a case study of measured CH 4 emissions from the southeastern United States (Douglas Reservoir) to illustrate the effects of varying surface area on reservoir-wide GHG estimates. Upscaled CH 4 emissions in Douglas Reservoir differed by nearly two-fold depending on the source of total surface area data and whether estimates accounted for seasonal fluctuations in surface area. During seasonal drawdown in Douglas Reservoir, relative littoral area varies non-linearly; periods of lower pool elevation (and thus larger relative littoral area) likely contribute disproportionately high CH 4 emission rates compared to the commonly sampled summer season when water levels are at full-pool elevation. Improved GHG monitoring and upscaling techniques require accounting for temporal variability in reservoir surface extent and littoral area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Long-term impact of light- and elevated temperature-induced degradation on photovoltaic arrays

Low levelized cost of electricity (LCOE) has been identified as critical for widespread adoption of photovoltaics (PV) without subsidies. Maintaining decades-long high-energy production is often an under-recognized opportunity in meeting cost goals because component lifetimes are not fully quantified at the time of manufacture. Whereas certain standardized tests minimize risk of early failure, there is little guidance to quantitatively predict degradation (or lack thereof) over decades, based on accelerated tests. In this article, we move toward bridging the understanding between indoor accelerated tests and outdoor performance data, with the goal of predicting energy yield with enough accuracy to inform financial decisions. Light- and elevated temperature-induced degradation (LETID) in p-type Si modules is analyzed in terms of impact on long-term module performance and thus LCOE. A method to predict the progression of LETID, using fixed kinetic constants and a numerical solution to the basic reaction rate equations, is detailed. Predictions are compared against both published data and that new to this study. These data include both indoor accelerated tests and fielded modules. We use the results in financial models to derive LCOE of modules in different climates with varying amounts of LETID, including uncertainty. Cost models based on the predictions indicate that LETID has a significant and climate-dependent impact on LCOE. Finally, we show that - even given the uncertainties identified in the study - these financial calculations can provide useful guidance to quantify risk based on accelerated test results. The analysis serves as an example of developing a predictive approach to PV reliability using physics of failure.

14 SOLAR ENERGY↗
DOE OSTI2022

Stand and environmental conditions drive functional shifts associated with mesophication in eastern US forests

There is a growing body of evidence that mesic tree species are increasing in importance across much of the eastern US. This increase is often observed in tandem with a decrease in the abundance and importance of species considered to be better adapted to disturbance and drier conditions (e.g., Quercus species). Concern over this transition is related to several factors, including the potential that this transition is self-reinforcing (termed “mesophication”), will result in decreased resiliency of forests to a variety of disturbances, and may negatively impact ecosystem functioning, timber value, and wildlife habitat. Evidence for shifts in composition provide broad-scale support for mesophication, but we lack information on the fine-scale factors that drive the associated functional changes. Understanding this variability is particularly important as managers work to develop site-and condition-specific management practices to target stands or portions of the landscape where this transition is occurring or is likely to occur in the future. To address this knowledge gap and identify forests that are most susceptible to mesophication (which we evaluate as a functional shift to less drought or fire tolerant, or more shade tolerant, forests), we used data from the USDA Forest Service Forest Inventory and Analysis program to determine what fine-scale factors impact the rate (change through time) and degree (difference between the overstory and midstory) of change in eastern US forests. We found that mesophication varies along stand and environmental gradients, but this relationship depended on the functional trait examined. For example, shade and drought tolerance suggest mesophication is greatest at sites with more acidic soils, while fire tolerance suggests mesophication increases with soil pH. Mesophication was also generally more pronounced in older stands, stands with more variable diameters, and in wetter sites, but plots categorized as “hydric” were often highly variable. Our results provide evidence that stand-scale conditions impact current and potential future changes in trait conditions and composition across eastern US forests. We provide a starting point for managers looking to prioritize portions of the landscape most at risk and developing treatments to address the compositional and functional changes associated with mesophication.

Woodbridge, Margaret↗
DOE OSTI2021

Microbial Biomass in Soils Receiving Varying Levels of Organic Inputs - A Data Compilation

This dataset contains soil microbial biomass data compiled from published laboratory studies and field manipulations of organic inputs to soils. Each study included at least three levels of a given type of organic input applied to a given type of soil. Studies were selected from the literature that had experimentally manipulated organic input rates to soils using unburned, non-synthetic material (e.g. plant litters, manures, and composts). Fifty-five (55) studies, conducted from 1997 to October 2019, were chosen and compiled into a dataset with 358 observations from 96 combinations of a given organic material applied to a given soil in agricultural systems, laboratory incubations, forests, or grasslands. • Studies were selected when biomass measurements used the chloroform fumigation-extraction procedure followed by carbon analysis, or, in two instances, studies estimated microbial biomass using phospholipid fatty acid analysis. • Studies that used fewer than three input rates or characterized microbial biomass via substrate induced respiration were excluded. • The temporal and spatial aspects of the laboratory and field studies were not essential to the analyses (Craig et al., 2021) and were not included with the measurement data. These data were used to assess the relationship between organic input amounts and microbial biomass. Complete reference citations are provided if additional information is needed. The compiled data are provided in two .csv files: (1) soil treatments and microbial biomass data and (2) literature citations. The “Study” field is the cross reference between files. Also provided is the R code for reproducing the results in the related publication.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗