Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2022

Cosmic rays in molecular clouds probed by H 2 rovibrational lines

Context. Low-energy cosmic rays (<1 TeV) play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H 2 . Their characterisation is therefore important both for the interpretation of observations and for the development of theoretical models. However, the methods used so far for estimating the cosmic-ray ionisation rate in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. Aims. Here, we refine and extend a previously proposed method to estimate the cosmic-ray ionisation rate in molecular clouds by observing rovibrational transitions of H 2 at near-infrared wavelengths, which are mainly excited by secondary cosmic-ray electrons. Methods. Combining models of interstellar cosmic-ray propagation and attenuation in molecular clouds with the rigorous calculation of the expected secondary electron spectrum and updated electron-H 2 excitation cross sections, we derive the intensity of the four H 2 rovibrational transitions observable in cold dense gas: (1—0)O(2), (1—0)Q(2), (1—0)S(0), and (1—0)O(4). Results. The proposed method allows the estimation of the cosmic-ray ionisation rate for a given observed line intensity and H 2 column density. We are also able to deduce the shape of the low-energy cosmic-ray proton spectrum impinging upon the molecular cloud. In addition, we present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar cosmic-ray proton spectrum. We finally comment on the capability of the James Webb Space Telescope to detect these near-infrared H 2 lines, which will make it possible to derive, for the first time, spatial variation in the cosmic-ray ionisation rate in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will finally be possible to test competing models of cosmic-ray propagation and attenuation in the interstellar medium, as well as compare cosmic-ray spectra in different Galactic regions

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Development and Test of High Temperature Surface Acoustic Wave Gas Sensors

The demand for sensors in hostile environments, such as power plant environments, aerospace environments, oil and gas extraction, and high-temperature metallurgy environments, has risen over the past decades in a continuous attempt to increase process control, improve energy and process efficiency in production, reduce operational and maintenance costs, increase safety, and perform condition-based maintenance in equipment and structures operating in high-temperature harsh-environment conditions. The increased reliability, improved performance, and development of new sensors and networks with a multitude of components, especially wireless networks, are the target for operation in harsh environments. Gas sensors, in particular H2 sensors, operating above 200 C are required in the instrumentation, process control and general safety of a number of industries including coal, natural gas, and nuclear power generation facilities, the aerospace and automotive industries, metallurgical production and defense-related applications. The surface acoustic wave (SAW) platform is a particularly promising option for high-temperature, harsh-environment gas sensing applications since the platform exhibits advantages, such as battery-free and wireless operation, small size, possibility for scale production using well-developed technologies from the semiconductor industry, and low cost of installation and operation. In this work, one-port SAW resonators (SAWRs) operating along five different orientations on a commercially available langasite (LGS) wafer employing Pt-Al2O3 electrodes and reflectors were designed, fabricated, and used as high-temperature H2 sensors. Two of the selected orientations were predicted and confirmed to have temperature-compensated operation above 150 C. A gas sensor test setup was developed, capable of gas cycling between N2, O2 and N2/H2 mixtures under extended high-temperature periods (up to 650 C for over 20 hours). Thin film Pt-Al2O3 was used as the electrode material for the transducers and reflectors capable of high-temperature operation, and also as H2 sensing film. In addition, yttria-stabilized zirconia (YSZ) thin films with Pt decoration were tested as sensing films aimed to enhance the SAWR sensor response to H2. The SAW devices were monitored in excess of 1700 hours in real-time during gas cycling sequences up to 600 C, leading to the following findings: i) the Pt-Al2O3 electrodes performed better for H2 sensing than the Pt-decorated YSZ sensing film, showing as much as 50% higher frequency variation response in the 200 C to 400 C range; ii) different crystallographic orientations operating on the same LGS wafer experienced different responses to H2 exposures up to 500 C; iii) the surface oxidation state of the SAWR sensors was shown to have an important impact on subsequent H2 exposure responses. In addition, a sensor system employing two LGS SAWRs, aligned along two different orientations, has been developed to simultaneously determine H2 presence and temperature. Finally, wireless interrogation of a SAWR sensor was successful within the gas cycling test fixture, and successful wireless H2 detection was achieved above 400 C.

Ayes Moncada, Armando↗
DOE OSTI2021

Data from: "Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine"

This archive contains data used to support conclusions drawn in “Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine”, by Reinhardt et al., 2011. Data were collected over one summer season in plots within the Alpine Treeline Warming Experiment (ATWE), before climate manipulations began. The experiment was located on Niwot Ridge, in the Front Range of the Colorado Rocky Mountains. This data package includes five comma-separated-values (.csv) files, five Microsoft Excel (.xlsx) files, one .pdf file, and two types of geospatial files: keyhole markup language (.kml), and ESRI shapefiles (.shp). .csv files can be opened using any simple text-editing software (such as Notepad and TextEdit), R, and Microsoft Excel. .xlsx files can only be opened using Microsoft Excel. The .pdf file can be opened using Adobe Acrobat Reader or any other compatible file viewing software. The .kml file can be opened using Google Earth and Google Maps, and shapefiles can be opened using any software compatible with the file type, such as ESRI’s ArcGIS suite and QGIS.Data archived contain gas exchange and plant physiology measurements, non-structural carbohydrate data, among others. Geospatial files are also provided for additional locational context. The files and their contents in this data package are summarized under "Data Summary" in the included Data User's Guide. All files (excluding geospatial) are available in both Microsoft Excel and in .csv format, and are indicated in the Data Summary list as well.-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Climate change is predicted to cause upward shifts in forest tree distributions, which will require seedling recruitment beyond current forest boundaries. However, predicting the likelihood of successful plant establishment beyond current species’ ranges under changing climate is complicated by the interaction of genetic and environmental controls on seedling establishment. To determine how genetics and climate may interact to affect seedling establishment, we transplanted recently germinated seedlings from high- and low-elevation provenances (HI and LO, respectively) of Pinus flexilis in common gardens arrayed along an elevation and canopy gradient from subalpine forest into the alpine zone and examined differences in physiology and morphology between provenances and among sites. Plant dry mass, projected leaf area and shoot:root ratios were 12–40% greater in LO compared with HI seedlings at each elevation. There were no significant changes in these variables among sites except for decreased dry mass of LO seedlings in the alpine site. Photosynthesis, carbon balance (photosynthesis/respiration) and conductance increased >2× with elevation for both provenances, and were 35–77% greater in LO seedlings compared with HI seedlings. There were no differences in dark-adapted chlorophyll fluorescence (Fv/Fm) among sites or between provenances. Our results suggest that for P. flexilis seedlings, provenances selected for above-ground growth may outperform those selected for stress resistance in the absence of harsh climatic conditions, even well above the species’ range limits in the alpine zone. This indicates that forest genetics may be important to understanding and managing species’ range adjustments due to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Predicting Plant-Soil Feedback in the Field: Meta-Analysis Reveals That Competition and Environmental Stress Differentially Influence PSF

Past research on plant-soil feedbacks (PSF), largely undertaken in highly controlled greenhouse conditions, has established that plant species differentially alter abiotic and biotic soil conditions that in turn affect growth of other conspecific and heterospecific individuals in that soil. Yet, whether feedbacks under controlled greenhouse conditions reflect feedbacks in natural environments where plants are exposed to a range of abiotic and biotic pressures is still unresolved. To address how environmental context affects PSF, we conducted a meta-analysis of previously published studies that examined plant growth responses to multiple forms of competition, stress, and disturbance across various PSF methodology. We asked the following questions: (1) Can competition, stress, and disturbance alter the direction and/or strength of PSF? (2) Do particular types of competition, stress, or disturbance affect the direction and/or strength of PSF more than others? and (3) Do methods of conducting PSF research (i.e., greenhouse vs. field experiments and whether the source of soil inoculum conditioning is from the field vs. greenhouse) affect plant growth responses to PSF or competition, stress, and disturbance, or their interactions? We discovered four patterns that may be predictive of what future PSF studies conducted under more realistic conditions might reveal. First, relatively little is known about how PSF responds to environmental stress and disturbance compared to plant-plant competition. Second, specific types of competition enhanced negative effects of soil microbes on plant growth, and specific environmental stressors enhanced positive effects of soil microbes on plant growth. Third, whether PSF experiments are conducted in the field or greenhouse can change plant growth responses. And, fourth, how the soil conditioning phase is conducted can change plant growth responses. With more detail than previously shown, these results confirm that environmental context writ large can change plant growth responses in PSF experiments. These data should aid theory and predictions for conservation and restoration applications by showing the relative importance of competition, stress, and disturbance in PSF studies over time. Lastly, these data demonstrate how variation in experimental methods can alter interpretation and conclusions of PSF studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2023

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Hybridization History and Repetitive Element Content in the Genome of a Homoploid Hybrid, Yucca gloriosa (Asparagaceae)

Hybridization in plants results in phenotypic and genotypic perturbations that can have dramatic effects on hybrid physiology, ecology, and overall fitness. Hybridization can also perturb epigenetic control of transposable elements, resulting in their proliferation. Understanding the mechanisms that maintain genomic integrity after hybridization is often confounded by changes in ploidy that occur in hybrid plant species. Homoploid hybrid species, which have no change in chromosome number relative to their parents, offer an opportunity to study the genomic consequences of hybridization in the absence of change in ploidy. Yucca gloriosa (Asparagaceae) is a young homoploid hybrid species, resulting from a cross between Yucca aloifolia and Yucca filamentosa. Previous analyses of ~11 kb of the chloroplast genome and nuclear-encoded microsatellites implicated a single Y. aloifolia genotype as the maternal parent of Y. gloriosa. Using whole genome resequencing, we assembled chloroplast genomes from 41 accessions of all three species to re-assess the hybrid origins of Y. gloriosa. We further used re-sequencing data to annotate transposon abundance in the three species and mRNA-seq to analyze transcription of transposons. The chloroplast phylogeny and haplotype analysis suggest multiple hybridization events contributing to the origin of Y. gloriosa, with both parental species acting as the maternal donor. Transposon abundance at the superfamily level was significantly different between the three species; the hybrid was frequently intermediate to the parental species in TE superfamily abundance or appeared more similar to one or the other parent. In only one case—Copia LTR transposons—did Y. gloriosa have a significantly higher abundance relative to either parent. Expression patterns across the three species showed little increased transcriptional activity of transposons, suggesting that either no transposon release occurred in Y. gloriosa upon hybridization, or that any transposons that were activated via hybridization were rapidly silenced. The identification and quantification of transposon families paired with expression evidence paves the way for additional work seeking to link epigenetics with the important trait variation seen in this homoploid hybrid system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Barajas-Olalde, César↗
DOE OSTI2021

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Impact of crowding on the diversity of expanding populations

Crowding effects critically impact the self-organization of densely packed cellular assemblies, such as biofilms, solid tumors, and developing tissues. When cells grow and divide, they push each other apart, remodeling the structure and extent of the population’s range. Recent work has shown that crowding has a strong impact on the strength of natural selection. However, the impact of crowding on neutral processes, which controls the fate of new variants as long as they are rare, remains unclear. Here, we quantify the genetic diversity of expanding microbial colonies and uncover signatures of crowding in the site frequency spectrum. By combining Luria–Delbrück fluctuation tests, lineage tracing in a novel microfluidic incubator, cell-based simulations, and theoretical modeling, we find that the majority of mutations arise behind the expanding frontier, giving rise to clones that are mechanically “pushed out” of the growing region by the proliferating cells in front. These excluded-volume interactions result in a clone-size distribution that solely depends on where the mutation first arose relative to the front and is characterized by a simple power law for low-frequency clones. Our model predicts that the distribution depends on a single parameter—the characteristic growth layer thickness—and hence allows estimation of the mutation rate in a variety of crowded cellular populations. Combined with previous studies on high-frequency mutations, our finding provides a unified picture of the genetic diversity in expanding populations over the whole frequency range and suggests a practical method to assess growth dynamics by sequencing populations across spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

State preparation of lattice field theories using quantum optimal control

Here, we explore the application of quantum optimal control (QOC) techniques to state preparation of lattice field theories on quantum computers. As a first example, we focus on the Schwinger model, quantum electrodynamics in 1+1 dimensions. We demonstrate that QOC can significantly speed up the ground state preparation compared to gate-based methods, even for models with long-range interactions. Using classical simulations, we explore the dependence on the interqubit coupling strength and the device connectivity, and we study the optimization in the presence of noise. While our simulations indicate potential speedups, the results strongly depend on the device specifications. In addition, we perform exploratory studies on the preparation of thermal states. Our results motivate further studies of QOC techniques in the context of quantum simulations for fundamental physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Predicting defects in laser powder bed fusion using in-situ thermal imaging data and machine learning

Variation in the local thermal history during the Laser Powder Bed Fusion (LPBF) process in Additive Manufacturing (AM) can cause micropore defects, which add to the uncertainty of the mechanical properties (e.g., fatigue life, tensile strength) of the built materials. In-situ sensing has been proposed for monitoring the AM process to minimize defects, but successful minimization requires establishing a quantitative relationship between the sensing data and the porosity, which is particularly challenging with a large number of variables (e.g., laser speed, power, scan path, powder property). Physics-based modeling can simulate such an in-situ sensing-porosity relationship, but it is computationally costly. In this work, we develop Machine Learning (ML) models that can use in-situ thermographic data to predict the micropore of LPBF stainless steel materials. This work considers two identified key features from the thermal histories: the time above the apparent melting threshold ($\tau$) and the maximum radiance ( T max ). These features are computed, stored for each voxel in the built material, and then used as inputs. The binary state of each voxel, either defective or normal, is the output. Different ML models are trained and tested for the binary classification task. In addition to using the thermal features of each voxel to predict its own state, the thermal features of neighboring voxels are also included as inputs. This is shown to improve the prediction accuracy, which is consistent with thermal transport physics around each voxel contributing to its final state. Among the models trained, the F1 scores on test sets reach above 0.96 for Random Forests. Feature importance analysis based on the ML models shows that T max is more important to the voxel state than $\tau$ . The analysis also finds that the thermal history of the voxels above the present voxel is more influential than those beneath it. Our study significantly extends the capability of using in-situ thermographic data to predict porosity in LPBF materials. Finally, since ML models are fast, they may play integral roles in the optimization and control of such AM technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deterministic and Statistical Modeling of a New Thermal Breakout Technology for Measuring the Maximum Horizontal In-Situ Stress

The current state-of-the-art technology for in-situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., “mini-fracs”). This methodology has achieved wide application in the oil and gas industry but has several limitations that often prevent successful in-situ stress measurements. One major limitation is that breakouts do not appear in all boreholes and are generally only a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to directly measure the maximum horizontal in-situ stress magnitude, the absence of borehole breakouts presents a major data gap for in-situ stress measurements. In response to this data gap, a new US Department of Energy (US DOE)-sponsored thermal breakout technology that will provide a method for thermally inducing borehole breakouts and allow the consistent measurement of the maximum horizontal stress magnitude is currently in development. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which can be directly correlated to the maximum horizontal stress magnitude. The first step in this project was an analytical modeling study of the thermal breakout process. Based on the Kirsch solution (Kirsch 1898), a deterministic and statistical analysis was performed on the pertinent parameters that influence the maximum horizontal stress calculation. As a result of the analysis, the findings indicate that the thermal breakout technology is feasible and provides improved accuracy and/or an enhanced ability to measure the maximum horizontal stress. Future work as part of this US DOE-sponsored project includes additional validation through more detailed numerical modeling, laboratory testing, and field testing of the thermal breakout technology.

58 GEOSCIENCES↗
DOE OSTI2023

Quantification of Energy Savings and Demand Reduction for a Heat Pump Integrated with Thermal Energy Storage (Final Report)

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries or installing new peaking power plants. Stakeholders and policy makers across the United States have expressed interests in promoting TES, as demonstrated by the US Department of Energy’s Grid-Interactive Efficient Buildings program and the efforts of various state legislatures. However, the cost value provided by TES are unclear. If reliable cost benefits were determined, stakeholders would have a clearer picture of the financial returns that can be gained from their investment in TES. The study in this report is conducted by ORNL with collaboration with Emerson the Helix Innovation Center. In the first part of this report, EnergyPlus was used to perform whole-building simulations for two residential buildings in Indianapolis and Atlanta. The HVAC system in both buildings were equipped with phase change material TES. The TES tank was charged in off-peak hours and discharged in peak hours to perform load shifting. First, the economic value implied by existing time-of-use (TOU) rates offered by utility companies was analyzed via whole-building simulation. Second, existing demand reduction (DR) incentives sourced from 3 different electrical grid administrators (i.e., California, Texas, and New England region) were surveyed to determine their implied value. The study suggests that the traditional value analysis that focuses on ROI for the building owner significantly undervalues TES technology making economic viability difficult. A more comprehensive value analysis that includes peak demand management and deferred capital for peaking power plants shows that TES should be economically viable but here the value is greater for the utility and requires large market penetration and aggregation to fully realize the benefits. Therefore, to facilitate commercialization, new business models are needed that include a broader range of stakeholders and distribute the value of TES proportionally. In the 2 nd part of this project, the benefits of a novel phase change material (PCM) integrated heat pump configuration were evaluated via detailed component based simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation. To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Using CESM-RESFire to understand climate–fire–ecosystem interactions and the implications for decadal climate variability

Large wildfires exert strong disturbance on regional and global climate systems and ecosystems by perturbing radiative forcing as well as the carbon and water balance between the atmosphere and land surface, while short- and long-term variations in fire weather, terrestrial ecosystems, and human activity modulate fire intensity and reshape fire regimes. The complex climate–fire–ecosystem interactions were not fully integrated in previous climate model studies, and the resulting effects on the projections of future climate change are not well understood. Here we use the fully interactive REgion-Specific ecosystem feedback Fire model (RESFire) that was developed in the Community Earth System Model (CESM) to investigate these interactions and their impacts on climate systems and fire activity. We designed two sets of decadal simulations using CESM-RESFire for present-day (2001–2010) and future (2051–2060) scenarios, respectively, and conducted a series of sensitivity experiments to assess the effects of individual feedback pathways among climate, fire, and ecosystems. Our implementation of RESFire, which includes online land–atmosphere coupling of fire emissions and fire-induced land cover change (LCC), reproduces the observed aerosol optical depth (AOD) from space-based Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products and ground-based AErosol RObotic NETwork (AERONET) data; it agrees well with carbon budget benchmarks from previous studies. We estimate the global averaged net radiative effect of both fire aerosols and fire-induced LCC at –0.59±0.52 W m –2 , which is dominated by fire aerosol–cloud interactions (–0.82±0.19 W m –2 ), in the present-day scenario under climatological conditions of the 2000s. The fire-related net cooling effect increases by ~170 % to –1.60±0.27 W m –2 in the 2050s under the conditions of the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Such considerably enhanced radiative effect is attributed to the largely increased global burned area (+19 %) and fire carbon emissions (+100 %) from the 2000s to the 2050s driven by climate change. The net ecosystem exchange (NEE) of carbon between the land and atmosphere components in the simulations increases by 33 % accordingly, implying that biomass burning is an increasing carbon source at short-term timescales in the future. High-latitude regions with prevalent peatlands would be more vulnerable to increased fire threats due to climate change, and the increase in fire aerosols could counter the projected decrease in anthropogenic aerosols due to air pollution control policies in many regions. We also evaluate two distinct feedback mechanisms that are associated with fire aerosols and fire-induced LCC, respectively. On a global scale, the first mechanism imposes positive feedbacks to fire activity through enhanced droughts with suppressed precipitation by fire aerosol–cloud interactions, while the second one manifests as negative feedbacks due to reduced fuel loads by fire consumption and post-fire tree mortality and recovery processes. These two feedback pathways with opposite effects compete at regional to global scales and increase the complexity of climate–fire–ecosystem interactions and their climatic impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Technical Evaluation of Accelerated Basin De-Inventory Material Addition to Sludge Batch 11 (Rev. 1)

The Accelerated Basin De-inventory (ABD) program involves discarding spent nuclear fuel that is currently stored in L-Basin to the Defense Waste Processing Facility (DWPF) for vitrification. The first ABD discards will occur during the preparation of Sludge Batch (SB) 11. Savannah River Mission Completion has requested that the Savannah River National Laboratory assess the technical gaps related to the increased gadolinium poisoning requirement and the impacts of performing the Low Temperature Aluminum Dissolution (LTAD) process in Tank 51 with H-Canyon discards present. The following summarizes the evaluation of the impacts of increasing the quantity of gadolinium (and related topics) from what was previously evaluated in the SRNL studies of gadolinium-poisoned ABD material solubility, the overall ABD flowsheet review, and increasing the fissile mass loading in glass: 1) Based on literature surveys, there is no indication that organic interactions with gadolinium will be significant at the high pH (typically >13) conditions of the Concentration, Storage, and Transfer Facilities. Any interactions of gadolinium with organics in DWPF are not expected to adversely impact DWPF or downstream facilities. Thus, there is little-to-no residual risk from organic interactions with gadolinium [Gap closed]; 2) Adding depleted uranium to ABD material, targeting 235 U enrichment of 4.90% within each transfer window, will mitigate potential impacts from an increase in soluble 235 U enrichment during sludge washing and LTAD. The plan to take advantage of previous transfers and allow 235 U enrichment of >5% during the final transfer window carries a risk that Tank 51 supernate will have a 235 U enrichment of >5%, which should be evaluated for acceptance; 3) Increasing the gadolinium mass ratio to 3.0:1 Gd: 235 U(eq SLU ) should lead to the same or higher partitioning of gadolinium into the solid phase within the DWPF Chemical Process Cell, resulting in both liquid and solid phases with expected partitioning of Gd consistent with the prior solubility study [Gap closed for SB11]; 4) There are no expected impacts on DWPF melt temperature and melter operations due to the minimal ~0.2 weight percent (wt%) increase in Gd concentration relative to previous sludge batches [Gap closed for SB11]; 5) As observed previously, Gd is expected to enter the off-gas system via physical entrainment, but at a slightly higher concentration than what was observed for SB9 melter off-gas pluggage deposits (0.07 wt%) [Gap closed for SB11] ; 6) There are no expected impacts on DWPF recycle or the Recycle Collection Tank glycolate destruction process. [Gap closed for SB11]; 7) Gd is projected to be a trace component in the SB11 glass (<0.5 wt%) and can be ignored for process control. Trace components do not significantly impact glass durability, thus the conclusions of the previous Product Consistency Test evaluation at a fissile mass loading of 2,500 g fissile/m3 glass still applies to SB11. The ~0.1 wt% increase in Gd2O3 concentration relative to the previous study will not impact the predictability of SB11 glass with the DWPF Product Composition Control System (PCCS) models for durability or the acceptability of glass according to the Waste Acceptance Product Specifications (WAPS) criterion for product consistency [Gap closed for SB11]; 8) No additional Toxicity Characteristic Leaching Procedure testing is necessary for SB11 and the hazardous waste specification of the SB11 DWPF waste form is unchanged after the addition of the ABD stream [Gap closed for SB11]. The following summarizes the evaluation of the impacts of adding two-thirds of the ABD material to Tank 51 prior to LTAD: 1) The addition of two-thirds of the ABD increases overall aluminum mass from 1.39×10 4 kg to 1.64×10 4 kg (15.5% ABD Al). The form of the insoluble portion of the Al resulting from ABD addition should be the more readily dissolved Al(OH) 3 and amorphous forms. The portion of the ABD aluminum that is processed by LTAD is expected to be completely soluble, thus requiring that less of the boehmite in the sludge be dissolved to reach the same Al target in the SB. [Gap closed for SB11]; The expected LTAD impact on other components, as related primarily to the components in ABD, are discussed. Gd is expected to remain insoluble during LTAD and not impact the solubility of other components. [Gap closed for SB11]; The addition of two-thirds of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 13,100 kg (63% ABD U) and the projected plutonium mass from 86.0 kg to 89.5 kg (3.9% ABD Pu). The addition of all of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 16,100 kg (70% ABD U) and the projected plutonium mass from 86.0 kg to 90.4 kg (5.3% ABD Pu). The 235 U enrichment will be ≤5%. The fissile uranium will be adequately poisoned by Gd and the fissile Pu will be adequately poisoned by Fe from the sludge. [Gap closed for SB11]; There is a low risk that ABD addition will impact the rheology or pumpability of the slurry. There is a low but higher risk of ABD addition prior to LTAD impacting the settling rate; Based on the evaluation of adding two-thirds of the ABD material and all of the ABD material prior to the LTAD process, there is no volume or mass limit that would need to be imposed on ABD additions prior to LTAD. [Gap closed for SB11]. Revision 1 of this report addresses a variation on the ABD additions and LTAD strategy where sodium hydroxide additions for LTAD may be performed intermittently or concurrently with an ABD addition window. The proposed change does not alter the conclusions of this evaluation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗