Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution operators”

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

Quenching in the Right Place at the Right Time: Tracing the Shared History of Starbursts, Active Galactic Nuclei, and Poststarburst Galaxies Using Their Structures and Multiscale Environments

This work uses multiscale environments and the structures of galaxies in the Sloan Digital Sky Survey as consistency checks on the evolution from starburst to quiescence at redshift z < 0.2. The environmental indicators include fixed-aperture mass overdensities (δ xMpc , x ϵ {0.5, 1, 2, 4, 8} h –1 Mpc), k-nearest neighbor distances, the tidal parameter, halo mass (M h ), and satellite/central classification. The residuals of the specific star formation rates (Δ SSFRs) are used to select starbursts (ΔSSFR > 0.6 dex, N ≈ 8600). Quenched poststarbursts (QPSBs) are selected using Hα < 3 Å emission and Hδ A > 4 Å absorption (N ≈ 750). The environments of starbursts and QPSBs are compared with those of active galactic nuclei (AGNs) and inactive galaxies of varying ΔSSFR. The environments of starbursts, AGNs, and QPSBs are unlike the environments of most quiescent galaxies (QGs). About 70%–90% of starbursts, AGNs with Hδ A > 4 Å, and QPSBs are centrals; ~80%–90% have M h < 10 13 M ⊙ , and only ~2%–4% have M h > 10 14 M ⊙ or live in clusters. Their M h and satellite fractions are also different from those of QGs. All QPSBs are matched to some SFGs, starbursts, AGNs, and QGs of similar M $\star$ , environments, concentration indices, and velocity dispersions. A significant fraction (~20%–30%) of starbursts cannot be matched to QPSBs or QGs. The implications are that: (1) some starbursts do not quench rapidly; (2) the satellite-quenching mechanisms operating in high-density environments cannot account for most QPSBs; (3) the evolution from starbursts to QPSBs to QGs is not the dominant path at z < 0.2; and (4) starbursts are not mainly triggered by tidal interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2022

Revealing room temperature ferromagnetism in exfoliated Fe 5 GeTe 2 flakes with quantum magnetic imaging

Van der Waals (vdW) material Fe 5 GeTe 2 , with its long-range ferromagnetic ordering near room temperature, has significant potential to become an enabling platform for implementing novel spintronic and quantum devices. To pave the way for applications, it is crucial to determine the magnetic properties when the thickness of Fe 5 GeTe 2 reaches the few-layers regime. However, this is highly challenging due to the need for a characterization technique that is local, highly sensitive, artifact-free, and operational with minimal fabrication. Prior studies have indicated that Curie temperature T C can reach up to close to room temperature for exfoliated Fe 5 GeTe 2 flakes, as measured via electrical transport; there is a need to validate these results with a measurement that reveals magnetism more directly. Here, we investigate the magnetic properties of exfoliated thin flakes of vdW magnet Fe 5 GeTe 2 via quantum magnetic imaging technique based on nitrogen vacancy centers in diamond. Through imaging the stray fields, we confirm room-temperature magnetic order in Fe 5 GeTe 2 thin flakes with thickness down to 7 units cell. The stray field patterns and their response to magnetizing fields with different polarities is consistent with previously reported perpendicular easy-axis anisotropy. Furthermore, we perform imaging at different temperatures and determine the Curie temperature of the flakes at ≈300 K. These results provide the basis for realizing a room-temperature monolayer ferromagnet with Fe 5 GeTe 2 . This work also demonstrates that the imaging technique enables rapid screening of multiple flakes simultaneously as well as time-resolved imaging for monitoring time-dependent magnetic behaviors, thereby paving the way towards high throughput characterization of potential two-dimensional (2D) magnets near room temperature and providing critical insights into the evolution of domain behaviors in 2D magnets due to degradation.

2D magnet↗
DOE OSTI2019

Vibrational coherence transfer illuminates dark modes in models of the FeFe hydrogenase active site

Within the conceptual framework of Redfield theory, the optical response function arises from the dynamical evolution of the system’s density operator, where nonunitary relaxation is encoded in the Redfield relaxation superoperator. In the conventional approach, the so-called secular approximation neglects terms that induce transitions between distinct coherences and among coherences and populations. The rationale is that these nonsecular terms are small in comparison to the far more dominant population relaxation and coherence dephasing contributions. Since two-dimensional infrared (2D-IR) spectroscopy has significant contributions arising from population relaxation and transfer pathways, it can be challenging to isolate signatures of the nonsecular relaxation. We report om this paper that in three diiron dithiolate hexacarbonyl complexes that serve as small-molecule models of the [FeFe] hydrogenase H-cluster subsite, a fortuitous vibrational energy structure enables direct and clear signatures of vibrational coherence transfer in alkane solution. This finding holds promise towards developing a molecularly detailed understanding of the mechanism of vibrational coherence transfer processes, thanks to the ease of synthesizing derivatives based on the chemical modularity of these well studied diiron compounds. In addition to the fundamental need to characterize coherence transfer in molecular spectroscopy, we find in this set of molecules a practical utility for the nonsecular dynamics: the ability to determine the frequency of an IR-inactive mode. A coherence generated during the waiting time of the 2D-IR measurement transfers to a coherence involving the single dark CO stretching mode, which modulates some peak amplitudes in the 2D spectrum, revealing its transient excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Global polarization of hyperons and spin alignment of vector mesons in quark matters

Relativistic heavy ion collider (RHIC) as a dedicated nuclear facility has made a few major discoveries in physics. This year marks the 30th year STAR Collaboration formation and the 23th year of STAR detector operation and data collection at RHIC. In the last two decades, STAR has collected many datasets, exhibiting scientific versatility and flexibility of the RHIC facility. The total dataset in the first year is less than 1 million good events, and currently there are about 1 billion events per dataset. The Global Hyperon Polarization was proposed in 2004. This immediately prompted the STAR Collaboration to search for this phenomenon from the early datasets. The null results were presented at Quark Matter Conference in Shanghai in 2006 and subsequently published. Although there were peripheral and continuous efforts in the following decade, no positive result has been observed experimentally. This situation changed in the following decade with the upgrade of high data rate and time-of-flight (TOF) detector and the progress of the Beam Energy Scan Phase I (BES-I). The experimental discoveries of the global polarization of hyperons in 2017 and the spin alignment of vector mesons in 2023 at RHIC-STAR confirm the theory which was established nearly twenty years ago. The theory and these measurements open the way to studying the properties of the hot and dense nuclear matter created in high-energy heavy ion collisions from a new degree of freedom, spin. We briefly review these discoveries from the proposals of theory to the experimental measurements, and summarize the related measurements at the existing facilities and the theoretical explanations to the original proposal. The basic understanding and the original proposal are still valid and fundamental, that is, the angular momentum of system can transform into a spin effect observable in experiment. However, it appears that in each case a new model is needed to explain the new experimental observation. We need a more basic theory to help us unify all these spin related phenomena. Over the past five years, STAR has successfully installed 3 new detectors and we have begun to see the physical analysis results from datasets with those new functions. What makes the STAR detector viable after 20 years of operation is its continuous evolution through successful upgrades, with new scientific programs added year by year. The next big thing is to forward upgrade a tracking system (3 layers of silicon strips and 4 layers of sTGC chambers) and a calorimetry system (electromagnetic and hadronic calorimeters). In addition to studying the spin structure of protons by using the polarized proton beams at RHIC, the upgrades also provide a unique ability to investigate the origin of Λ Global Polarization as a function of rapidity and rapidity (de-)correlations in Au+Au collisions.

Physics↗
DOE OSTI2019

Reduced Model for Direct Induction Startup Scenario Development on MAST-U and NSTX-U

A reduced semi-empirical model using time-dependent axisymmetric vacuum field calculations is used to develop the prefill and feed-forward coil current targets required for reliable direct induction (DI) startup on the new MA-class spherical tokamaks, MAST-U and NSTX-U. The calculations are constrained by operational limits unique to each device, such as the geometry of the conductive elements and active coils, power supply specifications and coil heating and stress limits. The calculations are also constrained by semi-empirical models for sufficient breakdown, current drive, equilibrium and stability of the plasma developed from a shared database. A large database of DI startup on NSTX and NSTX-U is leveraged to quantify the requirements for achieving a reliable breakdown (Ip ~ 20 kA). It is observed that without pre-ionization, STs access the large E/P regime at modest loop voltage (Vloop) where the electrons in the weakly ionized plasma are continually accelerating along the open field lines. This ensures a rapid (order millisecond) breakdown of the neutral gas, even without pre-ionization or high-quality field nulls. The timescale of the initial increase in Ip on NSTX is reproduced in the reduced model provided a mechanism for impeding the applied electric field is included. Most discharges that fail in the startup phase are due to an inconsistency in the evolution of the plasma current (Ip) and equilibrium field or loss of vertical stability during the burn-through phase. The requirements for the self-consistent evolution of the fields in the weakly and full-ionized plasma states are derived from demonstrated DI startup on NSTX, NSTX-U and MAST. The predictive calculations completed for MAST-U and NSTX-U illustrate that the maximum Ip ramp rate (dIp/dt) in the early startup phase is limited by the voltage limits on the poloidal field coils on MAST-U and passive vertical stability on NSTX-U.

Direct induction↗
DOE OSTI2018

The Influence of Galvanizing Parameters on the Structural Development of Zn-Al-Based Coatings

The evolution of some layers’ characteristics, deposited by galvanization, was monitored according to the technological parameters of the process. Two types of alloys from the Zn-Al and Zn-Al-Ti-B systems were used for deposition. The structure of the deposited layers was analyzed by optical and electron microscopy. The chemical composition, phases and compounds in the structures were identified. The mechanical properties of the structural constituents of the layers were determined by nanoindentation. These techniques highlighted the typical non-uniformities in terms of composition and properties in the layers, mainly due to the diffusion and distribution of iron from the steel substrate in the deposited layers. The evaluation of iron diffusion in the deposited layer was done by differential thermal analysis. The values of the critical temperatures (melting temperature, eutectic temperature and solidification temperature) were correlated with the parameters used during deposition. It was observed that the structure of the coatings is layered, with a specific order from the steel substrate to the edge. In these layers, the constituents are mixtures of soft solid solutions of Al in Zn (η and β) and hard intermetallic compounds. Moreover, increasing the immersion time in the galvanizing bath and the bath temperature leads to an increase in iron diffusion in the coating layer, which contributes to lower coating thickness, thus reducing its durability in operation.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Simulations of droplet combustion under gas turbine conditions

In various applications with recirculation, liquid droplets can be immersed in gases that may have a wide range of possible compositions, from pure air to hot combustion products. In order to gain fundamental understanding of the behaviour of individual droplets in vitiated air, numerical simulations of kerosene single droplet evaporation, autoignition, and combustion in conditions relevant to gas turbines have been performed. The droplet autoignition behaviour has been analysed in both physical and mixture fraction space for a wide range of vitiated air compositions and initial droplet diameters. Results show that the autoignition time delay decreases with increasing level of dilution with hot combustion products and decreasing initial droplet diameter. Chemistry is confined up to a radius of almost 10 initial droplet diameters and the location of autoignition is influenced by both the initial droplet diameter and the level of dilution. The time evolution of species in the gaseous phase after autoignition shows similar trends for all the diameters and dilution levels investigated here with the peak of temperature and OH mass fraction moving towards the droplet surface as a consequence of the balance between fuel production and consumption. In mixture fraction space, the location of the peaks of temperature and OH mass fraction after autoignition do not change in time whereas other intermediate species such as CH 2 O and pyrolysis products still exhibit a quite variable behaviour. The long-time flame structure has been compared with gaseous laminar counterflow simulations and, although qualitatively similar, the flame structure in the two configurations has differences with implications for flamelet combustion models used in spray combustion. The droplet evaporation, autoignition, and combustion behaviour has been summarized through a regime diagram showing the evaporation and autoignition time delays as a function of both initial droplet diameter and vitiated air dilution. This allows the identification of different states in the droplet combustion scenario and the introduction of critical values of dilution and initial droplet diameter beyond which single droplet rather than cloud combustion can occur, which can be exploited in the design of lean burn gas turbine combustion systems. In conclusion, the approach presented here can be easily extended to other conditions and fuels allowing the generation of regime diagrams for various operating conditions.

42 ENGINEERING↗
DOE OSTI2016

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Reduced model for direct induction startup scenario development on MAST-U and NSTX-U

A reduced semi-empirical model using time-dependent axisymmetric vacuum field calculations is used to develop the prefill and feed-forward coil current targets required for reliable direct induction (DI) startup on the new MA-class spherical tokamaks, MAST-U and NSTX-U. The calculations are constrained by operational limits unique to each device, such as the geometry of the conductive elements and active coils, power supply specifications and coil heating and stress limits. The calculations are also constrained by semi-empirical models for sufficient breakdown, current drive, equilibrium and stability of the plasma developed from a shared database. A large database of DI startup on NSTX and NSTX-U is leveraged to quantify the requirements for achieving a reliable breakdown (I p ~ 20 kA). It is observed that without pre-ionization, STs access the large E/P regime at modest loop voltage (V loop ) where the electrons in the weakly ionized plasma are continually accelerating along the open field lines. This ensures a rapid (order millisecond) breakdown of the neutral gas, even without pre-ionization or high-quality field nulls. The timescale of the initial increase in I p on NSTX is reproduced in the reduced model provided a mechanism for impeding the applied electric field is included. Most discharges that fail in the startup phase are due to an inconsistency in the evolution of the plasma current (I p ) and equilibrium field or loss of vertical stability during the burn-through phase. Here, the requirements for the self-consistent evolution of the fields in the weakly and fully-ionized plasma states are derived from demonstrated DI startup on NSTX, NSTX-U and MAST. The predictive calculations completed for MAST-U and NSTX-U illustrate that the maximum I p ramp rate (dI p /dt) in the early startup phase is limited by the voltage limits on the poloidal field coils on MAST-U and passive vertical stability on NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Development of process parameters and post-build conditions for qualification of LPBF 316 SS

To harness the potential of laser powder bed fusion (LPBF) 316H stainless steel (SS) for use in advanced nuclear reactors, extensive research efforts are needed to develop optimal laser processing parameters and appropriate heat treatments, taking into account various manufacturing platforms and locations. Furthermore, it is crucial to assess the material properties in environments relevant to reactor operating conditions. In collaboration with Oak Ridge National Laboratory (ORNL) and Los Alamos National Laboratory (LANL), Argonne National Laboratory (ANL) is dedicated to gaining insights into how the manufacturing process and post-build treatments impact the performance of LPBF 316H SS. This report provides an overview of ANL's research findings for FY23, focusing on LPBF 316H SS produced using a Renishaw AM400 laser system. Specifically, our research has concentrated on three key areas: 1) Optimum laser processing parameters: Leveraging the high-throughput printing technology, we systematically varied the laser power, exposure time and point distance within the same build. By measuring the as-printed porosity, we identified an optimum printing parameter window that produced materials with near-complete density; 2) Post-built treatment development: Materials printed with the optimum laser parameters underwent treatments including stress relief, solution annealing, and hot isostatic pressing. Subsequently, we fabricated tension, creep, and fatigue specimens from materials subjected to these different conditions. Over the next year, a series of scoping tests will be conducted to facilitate the selection of the most suitable post-build treatment condition; 3) Thermal aging effects on microstructure and properties: thermal aging at 550°C, 650°C and 750°C was conducted on LPBF 316H SS specimens up to 2500 h. Microstructural characterization was performed with electron microscopy, and the evolution of the dislocation cell structures and secondary phases was studied. Microhardness tests and tension tests were conducted to quantify changes in mechanical properties. The results revealed that the aged LPBF 316H SS exhibited a high density of Cr-rich and Mo-rich fine precipitates within grains due to the presence of a high density of dislocations and the globally distributed dislocation cell structures, which served as heterogeneous nucleation sites for secondary phases. Along grain boundaries, a Mo-Cr-rich phase, observed in 750°C-aged conditions, displayed significant thermal coarsening. Throughout the aging process, various mechanisms, including stress relaxation, dislocation cell structure recovery, solution hardening, and precipitation hardening, competed to influence materials strength. These results contribute towards establishing the technical basis for qualifying LPBF 316H SS for advanced nuclear reactor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact of ethanol blending into gasoline on aromatic compound evaporation and particle emissions from a gasoline direct injection engine

This study examined the interactions of ethanol with aromatic compounds on aromatic species evolution during gasoline evaporation and the consequent impacts on particle emissions from a single-cylinder, gasoline direct injection engine. From a chemical kinetic standpoint, the dilution of the aromatic species in gasoline by ethanol blending should reduce particle emissions because aromatics are the primary species that form particles during gasoline combustion. However, increased evaporative cooling from the ethanol and ethanol's effect on when aromatics in the fuel droplet (or surface film or pool) evaporate can drive an increase in particle emissions. The objective of this study was to quantify these competing effects. The results show that a combination of ethanol's increased evaporative cooling and impact on aromatic compound vapor-liquid equilibrium can extend droplet lifetime and cause aromatics to evaporate later in the droplet evaporation process than would be the case without ethanol. When fuels were burned in a direct-injection single-cylinder engine at high speed (2500 rpm) and relatively high load, ethanol blending caused an increase in particle emissions for fuels containing a low vapor pressure aromatic. At a lower speed (1500 rpm), no statistically significant increase in particle emissions was observed, likely because more time was available for evaporation and mixing to occur. In contrast, particle emissions from a fuel blended with a high vapor pressure aromatic were either insensitive to or reduced by ethanol blending, likewise to engine operating conditions. Current models for predicting particle emissions from gasoline engines, such as particulate matter index (PMI), generally lack variables for non-linear interactions between ethanol and gasoline hydrocarbons. To identify nonlinear interactions between predictor variables that might better explain the observed particle emissions, the least absolute shrinkage and selection operator (LASSO) regularized regression approach was employed. Among the predictor variables evaluated was the Yield Sooting Index (YSI) which is as an improved measure to quantify the aromatic's chemical tendency to generate particle emissions. Following this analysis, two variables were selected based on their frequent appearance. The first was ethanol mole percent multiplied by aromatic mole percent divided by aromatic vapor pressure [(EtOH * Aro%)/Aro VP@443 K]. The second was aromatic mole percent multiplied by aromatic YSI divided by aromatic molecular weight [(Aro% * Aro YSI)/Aro MW]. The results of linear regression including only these two combined explanatory variables achieved an r 2 = 0.959, substantially improved over the PMI model (r 2 = 0.688) and provides insight into how the PMI formalism might be modified to account for non-linear oxygenate-hydrocarbon coupling.

09 BIOMASS FUELS↗
DOE OSTI2023

Depletion-driven thermochemistry of molten salt reactors: review, method, and analysis

Molten salt reactors (MSRs) are innovative advanced nuclear reactors that utilize nuclear fuel by dissolving it in a high-temperature liquid salt. This unique feature differentiates MSRs from other types of reactors and allows for enhanced safety and economic performance. The liquid fuel also entails several multiphysics effects that can complicate reactor design and operation. One primary effect termed here as depletion-driven thermochemistry is a driving force in altering the multiphysics behavior of the reactor. Essentially, depletion-driven thermochemistry is the effect that fuel depletion has on changing the chemical redox potential of the fuel salt over time. As the fuel is consumed, the redox potential shifts toward a more oxidizing state. Without active control, the changing chemistry due to depletion increases corrosion thereby limiting reactor component lifetimes. Additionally, the changing redox potential of the fuel salt alters the vapor pressures of chemical species dissolved in the fuel salt. Changing vapor pressures of species in the fuel salt is an important parameter to understand when off-gassing volatile species during normal reactor operation, and for source term characterization during accident scenario transients. The present work represents a fundamental step toward modeling and coupling the driving physics (i.e., neutronics and chemistry) involved in altering the redox potential in an MSR. Here, the neutronic code Griffin models the depletion of the fuel-salt system, while the chemical equilibrium code Thermochimica calculates the thermochemical state of the isotopic inventory, using the Molten Salt Thermodynamic Database - Thermochemical (MSTDB-TC). These two codes are tightly coupled to predict the impact of fuel depletion in altering the chemistry in MSR systems. Redox potential control methods are discussed and can be modeled using this multiphysics approach. The vapor pressures of chemical species that could be extracted to an off-gas system, as determined by the reactor’s thermochemical state, are examined. The neutronics-chemistry coupling developed in this work is expected to have potential application for analyzing corrosion, source term evolution, and material safeguards in MSR systems. Lastly, suggestions for areas of further improvements of the models to expand these capabilities by incorporating other coupled physics effects is provided.

Walker, Samuel A.↗
DOE OSTI2023

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗
DOE OSTI2021

An experimental and one-dimensional modeling analysis of turbulent gas ejection in pre-chamber engines

In this work, experimental results from a study on the evolution of gas jets ejected through the orifices of a pre-chamber in a heavy-duty optical engine are presented. The work examines conditions without fuel inside the main-chamber, which helps to describe the dynamics of the ejected gas jets without the interference of subsequent combustion in the main-chamber. Experimental diagnostics consist of high-speed visible intensified imaging and low-speed infrared imaging. Additionally a one-dimensional gas jet model is used to characterize the spatial distribution of the ejected flow, including parameters such as tip penetration, which are then validated based on experimental results. Different stages in the ejection of pre-chamber jets are identified, with chemical activity restricted to a maximum distance of 5 to 10 orifice diameters downstream of the orifice as indicated by the recorded visible radiation. Sensitivity of cycle-to-cycle variations in pre-chamber jet development to the air-to-fuel ratio in the pre-chamber observed in the experiments is in most part attributed to the variations in the timing of combustion initiation in the pre-chamber. The influence of the ejection flow on the penetration of the gas jet on a cycle-to-cycle basis is presented using the one-dimensional model. The one-dimensional model also indicates that the local flow exhibits highest sensitivity to operating conditions during the start of ejection until the timing when maximum flow is attained. Differences that exist during the decreasing mass-flow ejection time-period tend to smear out in part due to the transient slowdown of the ejection process.

1-D gas jet model↗
DOE OSTI2025

Regulatory helix plays a key role in genetic ON-OFF switching for the 2’-deoxyguanosine sensing mRNA element

Transcriptional riboswitches, noncoding mRNA elements that operate in cis to regulate gene expression, have a promising potential in medicine, synthetic biology and directed evolution. They bind to cellular metabolites or metal ions with high specificity, leading to conformational rearrangements that facilitate the activation or premature termination of transcription for downstream genes. This elegant mechanism for feedback regulation of metabolic pathways has been identified in prokaryotes and a few in eukaryotes. Our chemical probing of the 2’-deoxyguanosine (2’-dG)-sensing riboswitch demonstrates that the overall conformational state of the full-length riboswitch (dGsw-fl) is unresponsive to the 2’-dG. Although binding proceeds as expected, dGsw-fl exclusively populates an OFF state of transcriptional inhibition. We chemically probed the structure of a known dGsw transcriptional intermediate (dGsw-int) to evaluate the possibility of a cotranscriptional regulatory role. Interestingly, apo dGsw-int adopts an alternative conformation in which a stable anti-terminator helix is formed, leading to an ON state where transcription can proceed. In the presence of 2’-dG, this anti-terminator helix is destabilized to produce a conformation reminiscent of the full-length, OFF-state dGsw. Using a fluorescence quenching assay, we demonstrate that binding 2’-dG to early transcriptional intermediates can inhibit the formation of the anti-terminator helix, locking dGsw in an OFF state. These data suggest that metabolite sensing occurs during a brief window of time between the synthesis of two transcriptional intermediates. Our studies indicate that dGsw does not function as a binary ON−OFF switch, but instead fine-tunes the transcription of downstream genes during RNA synthesis using key intermediates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

WE-AB-BRA-07: Operating Room Quality Assurance (ORQA) for Spine Surgery Using Known-Component 3D-2D Image Registration

Purpose: Intraoperative x-ray radiography/fluoroscopy is commonly used to qualitatively assess delivery of surgical devices (e.g., spine pedicle screws) but can fail to reliably detect suboptimal placement (e.g., breach of adjacent critical structures). We present a method wherein prior knowledge of the patient and surgical components is leveraged to match preoperative CT and intraoperative radiographs for quantitative assessment of 3D pose. The method presents a new means of operating room quantitative quality assurance (ORQA) that could improve quality and safety, and reduce the frequency of revision surgeries. Methods: The algorithm (known-component registration, KC-Reg) uses patient-specific preoperative CT and parametrically defined surgical component models within a robust 3D-2D registration method to iteratively optimize gradient similarity using the covariance matrix adaptation evolution strategy. Advances from previous work address key challenges to clinical translation: i) absolving the need for offline geometric calibration of the C-arm; and ii) solving multiple component bodies simultaneously, thereby allowing QA in a single step (e.g., spinal construct with 4–20 screws), rather than sequential QA of each component. Performance was tested in a spine phantom with 10 pedicle screws, and first results from clinical studies are reported. Results: Phantom experiments demonstrated median target registration error (TRE) of (1.0±0.3) mm at the screw tip and (0.7°±0.4°) in angulation. The simultaneous multi-body registration approach improved TRE from the previous (sequential) method by 42%, reduced outliers, and fits into the natural workflow. Initial application of KC-Reg in clinical data shows TRE of (2.5±4.5) mm and (4.7°±0.5°). Conclusion: The KC-Reg algorithm offers a potentially valuable method for quantitative QA of the surgical product, using radiographic systems that are already within the surgical arsenal. For spine surgery, the method offers a near-real-time independent check on the quality of surgical product, facilitating immediate revision if necessary and potentially avoiding postoperative morbidity and/or revision surgery. Gerhard Kleinszig and Sebastian Vogt are employees of Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗