Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2022

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure : Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.

diopside↗
DOE OSTI2024

Quantifying the Value of Technology and Policy Innovation in Water Resource Portfolios

Abstract Recent water infrastructure planning models demonstrate that explicitly accounting for hydrological changes in long‐term water planning reduces costs and increases system robustness. However, current models do not consider the effects of other changes occurring in the system over long time horizons, namely technology and policy innovation. By leveraging state‐of‐the‐art dynamic control techniques, we propose a suite of tools that explicitly relate hydrologic change, technology innovation, and policy interventions to system‐level water costs. We design robust and cost‐optimal strategies for water supply expansion under hydrological variability and hundreds of plausible innovation states, for example, treatment cost and diversification improvements, deployment times, and demand reduction campaigns. Our approach circumvents the need to pre‐assign probabilities to innovation scenarios and, instead, directly calculates the effects of technology, policy, and hydrology variations on system cost and planning strategies to identify avenues of impactful innovation. These tools can support technology development hubs, urban, state, and federal water planners in anticipating cost implications of technology and policy innovation at the local level, and prioritizing high‐impact innovation goals based on quantitative targets. Results for the case study of Santa Barbara, California identify high‐value innovation attributes for the city as well as combinations of innovation attributes across technologies and policies, and quantify their potential utility cost reductions.

13 HYDRO ENERGY↗
DOE OSTI2023

Upper critical solution temperature polymer assemblies via variable temperature liquid phase transmission electron microscopy and liquid resonant soft X-ray scattering

Abstract Here, we study the upper critical solution temperature triggered phase transition of thermally responsive poly(ethylene glycol)- block -poly(ethylene glycol) methyl ether acrylate- co -poly(ethylene glycol) phenyl ether acrylate- block -polystyrene nanoassemblies in isopropanol. To gain mechanistic insight into the organic solution-phase dynamics of the upper critical solution temperature polymer, we leverage variable temperature liquid-cell transmission electron microscopy correlated with variable temperature liquid resonant soft X-ray scattering. Heating above the upper critical solution temperature triggers a reduction in particle size and a morphological transition from a spherical core shell particle with a complex, multiphase core to a micelle with a uniform core and Gaussian polymer chains attached to the surface. These correlated solution phase methods, coupled with mass spectral validation and modeling, provide unique insight into these thermoresponsive materials. Moreover, we detail a generalizable workflow for studying complex, solution-phase nanomaterials via correlative methods.

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

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Aerosol Decline Accelerates the Increasing Extreme Precipitation in China

Extreme precipitation is becoming more intense and frequent. The increasing trends in extreme precipitation in China in warm season related to changes in aerosols and greenhouse gases (GHGs) are investigated using observations, reanalysis data and model simulations. A significant accelerating increase in extreme precipitation occurred around 2010, with the trend in accumulated extreme rainfall amount (R95pTOT) increasing from 2.88 mm per decade during 2000–2010 to 22.88 mm per decade during 2010–2023. The sudden acceleration of the increasing extreme precipitation is largely attributed to the reverse in aerosol trends associated with China’s clean air actions, which affects extreme precipitation through perturbing cloud microphysics and atmospheric dynamics, accounting for half of the change in R95pTOT trends. Future aerosol reduction to achieve carbon neutrality is shown to continue to intensify the extreme precipitation, which overweighs the effect induced by GHGs, highlighting the importance of aerosol changes in modulating future climate and weather extremes.

58 GEOSCIENCES↗
DOE OSTI2022

On-the-fly autonomous control of neutron diffraction via physics-informed Bayesian active learning

We demonstrate the first live, autonomous control over neutron diffraction experiments by developing and deploying ANDiE: the autonomous neutron diffraction explorer. Neutron scattering is a unique and versatile characterization technique for probing the magnetic structure and behavior of materials. However, instruments at neutron scattering facilities in the world is limited, and instruments at such facilities are perennially oversubscribed. We demonstrate a significant reduction in experimental time required for neutron diffraction experiments by implementation of autonomous navigation of measurement parameter space through machine learning. Prior scientific knowledge and Bayesian active learning are used to dynamically steer the sequence of measurements. We show that ANDiE can experimentally determine the magnetic ordering transition of both MnO and Fe 1.09 Te all while providing a fivefold enhancement in measurement efficiency. Furthermore, in a hypothesis testing post-processing step, ANDiE can determine transition behavior from a set of possible physical models. ANDiE's active learning approach is broadly applicable to a variety of neutron-based experiments and can open the door for neutron scattering as a tool of accelerated materials discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Warming increased bark beetle-induced tree mortality by 30% during an extreme drought in California

Quantifying the responses of forest disturbances to climate warming is critical to our understanding of carbon cycles and energy balances of the Earth system. The impact of warming on bark beetle outbreaks is complex as multiple drivers of these events may respond differently to warming. Using a novel model of bark beetle biology and host tree interactions, we assessed how contemporary warming affected western pine beetle (Dendroctonus brevicomis) populations and mortality of its host, ponderosa pine (Pinus ponderosa), during an extreme drought in the Sierra Nevada, California, United States. When compared with the field data, our model captured the western pine beetle flight timing and rates of ponderosa pine mortality observed during the drought. In assessing the influence of temperature on western pine beetles, we found that contemporary warming increased the development rate of the western pine beetle and decreased the overwinter mortality rate of western pine beetle larvae leading to increased population growth during periods of lowered tree defense. We attribute a 29.9% (95% CI: 29.4%–30.2%) increase in ponderosa pine mortality during drought directly to increases in western pine beetle voltinism (i.e., associated with increased development rates of western pine beetle) and, to a much lesser extent, reductions in overwintering mortality. These findings, along with other studies, suggest each degree (°C) increase in temperature may have increased the number of ponderosa pine killed by upwards of 35%–40% °C-1 if the effects of compromised tree defenses (15%–20%) and increased western pine beetle populations (20%) are additive. Due to the warming ability to considerably increase mortality through the mechanism of bark beetle populations, models need to consider climate's influence on both host tree stress and the bark beetle population dynamics when determining future levels of tree mortality.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Coupling of electroelastic dynamics and direct contact ultrasonic drying formulation for annular piezoelectric bimorph transducers

A newly developed technique for drying clothes without thermal energy has been developed through the utilization of ultrasonic vibrations from piezoelectric transducers. The novel technique incorporates the actuation of a thin stainless steel disk in contact with wet fabric via annular piezoelectric rings, where water in the liquid form is atomized, transported through microchannels in the disk, and ejected as a mist. In such a system, resonance matching between the actuation portion of the transducer and the portion contacting fabric must be realized, with theoretical results from the developed electromechanical model showing a reduction in energy consumption by 50% when resonance matching is achieved. The electrically coupled distributed parameter model for an annular bimorph piezoelectric transducer is developed for optimization of ultrasonic drying technology. The thickness mode vibrations are shown to dominate the behavior of the system, where the analytically developed model can be optimized to increase the output acceleration of the transducer, thus increasing drying performance. Here, the electromechanical equation developed will be connected to the drying rates of fabrics in contact with said vibrations, where the novelty of the coupled equations and its description of the physics of ultrasonic drying will be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Distributed deep reinforcement learning for simulation control

Abstract Several applications in the scientific simulation of physical systems can be formulated as control/optimization problems. The computational models for such systems generally contain hyperparameters, which control solution fidelity and computational expense. The tuning of these parameters is non-trivial and the general approach is to manually ‘spot-check’ for good combinations. This is because optimal hyperparameter configuration search becomes intractable when the parameter space is large and when they may vary dynamically. To address this issue, we present a framework based on deep reinforcement learning (RL) to train a deep neural network agent that controls a model solve by varying parameters dynamically. First, we validate our RL framework for the problem of controlling chaos in chaotic systems by dynamically changing the parameters of the system. Subsequently, we illustrate the capabilities of our framework for accelerating the convergence of a steady-state computational fluid dynamics solver by automatically adjusting the relaxation factors of the discretized Navier–Stokes equations during run-time. The results indicate that the run-time control of the relaxation factors by the learned policy leads to a significant reduction in the number of iterations for convergence compared to the random selection of the relaxation factors. Our results point to potential benefits for learning adaptive hyperparameter learning strategies across different geometries and boundary conditions with implications for reduced computational campaign expenses 4 4 Data and codes available at https://github.com/Romit-Maulik/PAR-RL . .

42 ENGINEERING↗
DOE OSTI2026

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

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

Beam dynamics corrections in the measurement of the anomalous precession frequency at the Muon $g-2$ experiment at Fermilab

The Muon $g-2$ experiment at Fermilab (E989) aims to measure the anomalous magnetic moment of the muon with an accuracy of 140 ppb (parts per billions). This accuracy, obtained by adding in quadrature a statistical and a systematic contribution of comparable value (100 ppb), will allow to reduce the experimental uncertainty (from the previous E821 experiment at BNL \cite{bnl}) of a factor of 4, and represents one of the most precise tests of the Standard Model (SM) theory of elementary particles. The first result on the Run-1 dataset \cite{prl} was released on April 7, 2021, showing a very good agreement with the previous result from BNL experiment, with a slightly better uncertainty. The corresponding experimental average increases the significance of the discrepancy between the measured and Standard Model prediction of 4.2$\sigma$ \cite{white_paper}.\\ \noindent The measured quantity is the muon magnetic anomaly a$_\mu$=$\frac{g_\mu-2}{2}$ where g$_\mu$ is the gyromagnetic factor of the muon. Dirac's equation predicts g$_\mu$ = 2, while radiative corrections, dominated by the QED contribution due to an exchange of a virtual photon, causes a per-mille correction on this quantity. By including all the SM contributions, a$_\mu$ is known at 370 ppb. The E989 experiment measures $a_\mu$ injecting positive muons with momentum of 3.1 GeV/$c$ polarized longitudinally in bunches (called $fills$) with an average rate of 12 Hz, in a storage ring of 14 meters diameter. Due to the parity violation in the weak muon decay, high energy positrons produced are emitted preferably in the muon's spin direction. By counting the number of positrons with energy greater then 1.7 GeV in function of the time, the frequency precession of the muon spin is measured, that together with the measurement of the magnetic field, allows to extract $a_\mu$. The positrons are detected with 24 electromagnetic calorimeters, that measure the energy and the arrival time of the positrons, each made of 54 crystals of lead fluoride (PbF$_2$) read by silicon photomultipliers (SiPM). Together with calorimeters, two tracking detectors are used to make non-destructive measurements of the muon beam characteristic by reconstructing the muon decay position extrapolating backward the decay positrons. The knowledge of the beam motion inside the ring plays a fundamental role in the analysis of $a_{\mu}$, where the measured anomalous precession frequency $\omega_a$ must be corrected for four main beam dynamics effects. A first correction is associated with the presence of an electric field responsible of the vertical focusing on the storage ring, where vertical direction is orthogonal to the orbit and horizontal direction is along the storage ring radius. Due to the oscillations in the vertical direction, the so called vertical betatron oscillations, a second correction is necessary to account for an average angle associated with the muon motion off the ideal orbit. A third correction is caused by lost muons in the ring which have a different spin phase at the injection respect to the decay ones. Finally, due to the correlation between the vertical and horizontal beam motion of the muons and the acceptance of the calorimeter, a correction (called ``phase-acceptance") arises. In Run1, due to the presence of two damaged resistors in one electrostatic quadrupole, this correction had a prominent role. The replacement of the damaged resistors before Run2 decreased this effect by one order of magnitude, and a further reduction in Run-3 was provided by an improved orbit.\\ \noindent The work of this Thesis focuses on the beam dynamics corrections on $\omega_a$. Due to the prominent role in Run-1 a special attention has been put to the phase-acceptance correction\footnote{The effect of this correction on the E821 BNL measurement of the $g-2$ has been evaluated to be within the quoted systematic error.}. Chapter \ref{ch:anomalous_magnetic_moment} introduces the anomalous magnetic moment of the muon. Chapter \ref{ch:early_experiment} describes the history of the Muon $g-2$ experiments. Chapter \ref{ch:theory} discusses the Standard Model prediction and possible new physics scenario. Chapter \ref{ch:muon_g2_experiment} describes the E989 experiment, whose experimental technique and the apparatus are discussed focusing on the improvements needed to reach the final goal on $a_\mu$ measurement. The original work of the Thesis is discussed in the last 5 chapters. Chapter \ref{ch:experiment_technique} presents the analysis technique to extract $\omega_a$, Chapter \ref{ch:beam_dynamics} describes the beam dynamics in the E989 experiment, Chapter \ref{ch:beam_dynamics_correction} discusses the beam dynamics correction to $\omega_a$, and Chapter \ref{sec:pa_corr} focuses on the phase acceptance correction. Finally Chapter \ref{ch:run23_analysis_improvements} contains the Run-2/3 improvements which are expected to allow for an increased precision ons the beam dynamics co...

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Disentangling climate and policy uncertainties for the Colorado River post-2026 operations

Abstract Lakes Mead and Powell in the Colorado River Basin underpin water and hydropower supply for the western United States. While the policies currently regulating the basin will expire by 2026, planning remains challenging due to intertwined climate variability and policy uncertainties. Based on streamflow projections from 10 dynamically downscaled CMIP6 global climate models and unique methods that add and remove internal variability, we evaluate future conditions at Powell and Mead under existing and alternative policies. Due to projected streamflow declines, under existing policy, both reservoirs will face substantial risks (>80% likelihood) of reaching dead pool before 2060. Adopting recently proposed alternative policies reduces but doesn’t eliminate such risks. All policies also exhibit tipping points where reservoir levels can change rapidly with a slight change in streamflow. A sustainable policy may require larger reductions to further reduce the reservoirs’ dead pool risks and provide better buffers from sudden changes.

Science & Technology - Other Topics↗
DOE OSTI2022

liionpack: A Python package for simulating packs of batteries with PyBaMM

Electrification of transport and other energy intensive activities is of growing importance as it provides an underpinning method to reduce carbon emissions. With an increase in reliance on renewable sources of energy and a reduction in the use of more predictable fossil fuels in both stationary and mobile applications, energy storage will play a pivotal role and batteries are currently the most widely adopted and versatile form. Therefore, understanding how batteries work, how they degrade, and how to optimize and manage their operation at large scales is critical to achieving emission reduction targets. The electric vehicle (EV) industry requires a considerable number of batteries even for a single vehicle, sometimes numbering in the thousands if smaller cells are used, and the dynamics and degradation of these systems, as well as large stationary power systems, is not that well understood. As increases in the efficiency of a single battery become diminishing for standard commercially available chemistries, gains made at the system level become more important and can potentially be realised more quickly compared with developing new chemistries. Mathematical models and simulations provide a way to address these challenging questions and can aid the engineer and designers of batteries and battery management systems to provide longer lasting and more efficient energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Evaluating the Effectiveness of Soil Profile Rehabilitation for Pluvial Flood Mitigation Through Two-Dimensional Hydrodynamic Modeling

Pluvial flooding, driven by increasingly impervious surfaces and intense storm events, presents a growing challenge for urban areas worldwide. In Baltimore City, MD, USA, climate change, rapid urbanization, and aging stormwater infrastructure are exacerbating flooding impacts, resulting in significant socio-economic consequences. This study evaluated the effectiveness of a soil profile rehabilitation scenario using a 2D hydrodynamic modeling approach for the Tiffany Run watershed, Baltimore City. This study utilized different extreme storm events, a high-resolution (1 m) LiDAR Digital Terrain Model (DTM), building footprints, and hydrological soil data. These datasets were integrated into a fully coupled 2D hydrodynamic model, the City Catchment Analysis Tool (CityCAT), to simulate urban flood dynamics. The pre-soil rehabilitation simulation revealed a maximum water depth of 3.00 m in most areas, with hydrologic soil groups C and D, especially downstream of the study area. The post-soil rehabilitation simulation was targeted at vacant lots and public parcels, accounting for 33.20% of the total area of the watershed. This resulted in a reduced water depth of 2.50 m. Additionally, the baseline runoff coefficient of 0.49 decreased to 0.47 following the rehabilitation, and the model consistently recorded a peak runoff reduction rate of 4.10 across varying rainfall intensities. The validation using a contingency matrix demonstrated true-positive rates of 0.75, 0.50, 0.64, and 0 for the selected events, confirming the model’s capability at capturing real-world flood occurrences.

Baltimore City↗
DOE OSTI2025

Rate expressions and kinetic parameters for metal ferrites in relation to applications of fossil fuel conversion to hydrogen: Part 1 of 2

Here, the goal of the present work was to provide the necessary reaction emulation information to enable detailed process simulation of a chemical looping H 2 production system from fossil fuels using CaFe 2 O 4 . This specifically pertained to the necessary kinetic data, reaction model development, and model rate parameters required for reaction emulation in both reducing and oxidizing environments. A logical methodology was defined, which included discretization of the reaction network, establishing a core model for reaction emulation that could be adapted based on the system phenomena, and development of a rate parameter regression tool designed around the core model. An extensive array of data sets was acquired by which parametric regressions were performed. The work presented and tabulated a comprehensive set of rate parameters for the reduction and oxidation reactions of CaFe 2 O 4 and descendent phases of Ca 2 Fe 2 O 5 , FeO, Fe 3 O 4 , Fe, and CaO to emulate reaction behavior in a looping-based process environment. This included direct reduction using CH 4 , H 2 , and CO, and direct oxidation reactions with steam, CO 2 and O 2 . Dynamic equilibrium was quantified for reactions that could utilize H 2 O and CO 2 as soft oxidants to re-saturate lattice oxygen in the depleted structure/phases. The kinetics associated with the oxidative mechanisms with the soft oxidants were quantified and compared to those of the reducing counterparts. The analysis provided critical insight to emulate reactions for a process that seeks to use natural gas (NG) or other fossil fuels as a direct reductant for the end goal of H 2 production.

calcium ferrite oxygen carriers↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗