Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stabilizing a strongly nonlinear structure through shaker dynamics in fixed frequency voltage control tests

Bifurcations are commonly encountered during force controlled swept and stepped sine testing of nonlinear structures, which generally leads to the so-called jump-down or jump-up phenomena between stable solutions. There are various experimental closed-loop control algorithms, such as control-based continuation and phase-locked loop, to stabilize dynamical systems through these bifurcations, but they generally rely on specialized control algorithms that are not readily available with many commercial data acquisition software packages. A recent method was developed to experimentally apply sequential continuation using the shaker voltage that can be readily deployed using commercially available software. By utilizing the stabilizing effects of electrodynamic shakers and the force dropout phenomena in fixed frequency voltage control sine tests, this approach has been demonstrated to stabilize the unstable branch of a nonlinear system with three branches, allowing for three multivalued solutions to be identified within a specific frequency bandwidth near resonance. Recent testing on a strongly nonlinear system with vibro-impact nonlinearity has revealed jumping behavior when performing sequential continuation along the voltage parameter, like the jump phenomena seen during more traditional force controlled swept and stepped sine testing. Here, this paper investigates the stabilizing effects of an electrodynamic shaker on strongly nonlinear structures in fixed frequency voltage control tests using both numerical and experimental methods. The harmonic balance method is applied to the coupled shaker-structure system with an electromechanical model to simulate the fixed voltage control tests and predict the stabilization for different parameters of the model. The simulated results are leveraged to inform the design of a set of experiments to demonstrate the stabilization characteristics on a fixture-pylon assembly with a vibro-impact nonlinearity. Through numerical simulation and experimental testing on two different strongly nonlinear systems, the various parameters that influence the stability of the coupled shaker-structure are revealed to better understand the performance of fixed frequency voltage control tests.

42 ENGINEERING↗
DOE OSTI2023

Aerosol Indirect Effects on Water Vapor in the UTLS of Typhoon Saomai (2006)

The evolution of landfalling Typhoon Saomai (2006) is simulated using the Weather Research and Forecasting model with spectral bin microphysics scheme, showing tropical cyclone (TC) could vertically transport water vapor from troposphere to upper troposphere and lower stratosphere (UTLS). Three aerosol sensitivity experiments are conducted to explore the indirect effects of aerosol served as cloud condensation nuclei (CCN) on overshooting convection (OTC) and water vapor content in UTLS, namely MAR (maritime), MIX (semi‐continental) and CON (continental). Results revealed that CON simulation with high initial aerosol concentration showing invigorated convection at the periphery of the TC, which interfered with the reformation of TC eyewall, leading to TC weakening finally. In addition, CON simulated the least water content in UTLS through suppressing the number of convection overshooting the tropopause level. Most overshooting events take place during TC developing stage in the eyewall region before landfall. Higher levels of CCN concentration could lead to decreased number of OTC in TC, consistently reduced the water vapor content in UTLS. CON and MIX experiments simulated respectively 21.9% (16.2%) and 9.96% (1.5%) less OTC (water vapor content in UTLS) than MAR simulation. Further, water vapor content in UTLS is affected by the interaction of OTC and ice‐phase microphysics. The interesting point is that the CCN sensitivity experiments show that water vapor content in UTLS of a TC is sensitive to the CCN concentration in the boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterizing the GD-1 Stream with DESI DR2 Data: Thin Stream and Hot Cocoon

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($σ_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($σ_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.

Jarvis, Emma [Toronto U.] (ORCID:0009000656127336)↗
DOE OSTI2025

The new ITER baseline, research plan and open R&D issues

A new baseline (NB) has been proposed by the ITER Project to ensure a robust achievement of the Projects’ goals, in view of past challenges including delays incurred due to the Covid-19 pandemic, technical challenges in completing first-of-a-kind components and in nuclear licensing. The NB includes modifications to the configuration of the ITER device and its ancillaries (e.g. change from beryllium to tungsten as first wall material, modification of the heating and current drive mix, etc.) as well as additional testing of components (e.g. toroidal field coils) or phased installation (start with inertially cooled first wall before later installation of the final actively water-cooled components) to minimise operational risks. In the NB, the ITER research plan (IRP) will be divided into three main phases: (a) start of research operation, with 40 MW of ECH and 10 MW of ICH, which will focus on the demonstration of 15 MA operation in L-mode, commissioning of all required systems, including disruption mitigation, and the demonstration of H-mode plasma operation in deuterium; (b) DT-1, with 60–67 MW of ECH, 33 MW of neutral beam injection (NBI) and 10–20 MW of ICH, which will demonstrate robust operation in high confinement H-mode plasmas in DT up to Q ⩾ 10 and for burn durations of 300–500 s within an accumulated neutron fluence of ∼1% of the ITER machine’s lifetime total, and; (c) DT-2, with up to 67 MW of ECH, up to 49.5 MW of NBI and up to 20 MW of ICH, with the ITER tokamak and ancillaries in their final configuration to demonstrate routine operation in DT plasmas at high Q and the Q ⩾ 5 long-pulse and steady-state scenarios to the final neutron fluence and to perform R&D on nuclear fusion reactor issues. The logic, physics basis, modelling and experimental evaluations carried out to support the NB and the associated IRP are described. These include the impact of the tungsten wall on plasma scenarios and associated risk mitigation measures, as well as the optimisation of the tokamak components and ancillaries to minimise Project risks. Open R&D issues related to these evaluations and mitigation measures are also described together with experimental, modelling and validation activities required to address them.

ITER↗
DOE OSTI2025

Searches for New Long-Lived Particles and Upgrade to the ATLAS Inner Detector (Final Technical Report)

The search for new fundamental particles is one of the defining goals of the Large Hadron Collider (LHC). The discovery of the Higgs Boson by the ATLAS and CMS collaborations provided the capstone of the Standard Model of particle physics, but outstanding questions remain. Why does the Higgs boson have a mass of 125 GeV when its natural mass would be many orders of magnitude larger? Is there a universal symmetry which unites all three forces described by the Standard Model? Can that symmetry be extended to include gravity? Is dark matter, evidenced by astronomical observations, made of a particle that interacts via Standard Model forces with the rest of matter? Together, these motivations provide compelling arguments that new physical processes await discovery. This project addressed some outstanding questions about the fundamental particles and their interactions with the ATLAS experiment at the Large Hadron Collider. In particular, the project improved the discovery potential for new, long- lived particles produced via electroweak processes in proton-proton collisions and set world-leading limits on their existence for certain values of their potential mass and lifetime. To achieve this, the project developed new data analysis methods, developed new triggers to select events with new long-lived particles during data-taking of the ATLAS experiment, and analyzed the largest proton–proton collision dataset ever produced. The project also supported significant development of the data acquisition software for the upgrade to the ATLAS inner detector, the Inner TracKer (ITk). The upgrade of the ATLAS inner detector is essential to the success of the entire Phase II physics program on ATLAS. Personnel supported by the project provided support for integration, assembly, and testing of the inner two layers of the ITk pixel system during its prototype and pre-production phase. Four PhD students and two post-doctoral scholars were supported by the grant and received invaluable scientific training as part of the research endeavor. The students and postdocs gained essential professional skills in the areas of advanced data analysis techniques, statistical analysis of data and simulation, programming in C++ and Python, hardware and instrumentation development, and presentation and collaboration skills. Additionally, approximately ten undergraduate students supported through other funding sources participated in research activities synergistic with the goals of this project, receiving essential mentorship from the personnel supported by this project.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗
DOE OSTI2025

Measurement of three-dimensional inclusive muon-neutrino charged-current cross sections on argon with the MicroBooNE detector

We report the measurement of the triple-differential cross section d 3 σ/dE vis d cos(θ μ )dP μ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4 x 10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. The mapping from reconstructed kinematics to truth quantities is validated within uncertainties by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after applying a conditional constraint from the muon momentum distribution in data. The success of this validation provides confidence that the energy transfer in the MicroBooNE detector is well-modeled within simulation uncertainties, enabling a reliable unfolding to a triple-differential cross section defined at the nominal neutrino flux over muon momentum, muon scattering angle, and visible neutrino energy. This validation not only supports accurate cross-section extraction, but also establishes a critical foundation for tuning interaction models used in future neutrino oscillation measurements. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Formation, cooling history and age of impact events on the IIE iron parent body: Evidence from the Miles meteorite

Most iron meteorites formed in planetary cores during differentiation, but the IIE iron meteorites have chemical and physical features that are inconsistent with this origin. By combining mineral chemistry, mineral modes and three-dimensional petrography, we reconstruct the bulk chemistry of the felsic silicate-bearing Miles IIE iron meteorite and demonstrate that the silicate inclusion compositions are similar to partial melts produced experimentally from an H chondrite composition. We use the reconstructed bulk composition, mineralogy and thermodynamic modelling to show that melting above ~ 1200°C under reducing conditions formed metal (Fe-Ni alloy) and felsic silicate partial melts. Upon cooling, the melts crystallized Mg-rich pyroxenes, Na- and K-rich feldspars, and tridymite. Importantly, this mechanism enriches cosmochemically volatile elements (i.e., those with a 50% condensation temperature of ~430–830°C, like Na and K) to the level found in the felsic silicate inclusions. The presence of crystallographically disordered srilankite (only stable above 1160°C) and an absence of Widmanstätten texture require both high peak temperatures and rapid cooling, which cannot be explained by core formation. Instead, they point to small melt volumes, a transient heat pulse, and small thermal mass, and imply efficient physical segregation of silicate and metallic melts through buoyancy separation followed by rapid cooling that arrested the separation of metal and silicate liquid phases. In situ 207 Pb/ 206 Pb age of 4542.3 ± 4.0 Ma in Zr-oxide and phosphate minerals dates the melting event that formed the silicate inclusions. This age aligns with the earliest ages found in other IIE iron meteorite silicates and requires a heating event ~25 million years after the solar system formed. We found 39 Ar/ 40 Ar ages of 3495 ± 52 Ma (low-T) and 4303 ± 7 Ma (high-T) in a K-feldspar grain, with the 3495 Ma age aligning with later thermal events recorded in other IIE iron meteorites. Dating reveals the complex petrogenetic and thermal history of Miles and the IIE iron meteorites. This is the first IIE iron meteorite found to record evidence of heating at 4.5 and 3.5 Ga likely from impact events. We propose that high-velocity impact(s) into an iron-rich, porous chondritic parent body at ~4.54 Ga produced immiscible metal and silicate melts that cooled rapidly and trapped low density silicate inclusions within high density metal. Other IIE irons that formed at lower peak temperatures (900–1000°C) contain chondritic silicate inclusions and relict chondrules, supporting this conceptual model. In conclusion, this petrogenesis is consistent with thermodynamic modelling, experimental data and the wide range of peak temperatures and cooling rates observed in the IIE iron meteorites.

58 GEOSCIENCES↗
DOE OSTI2025

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-throughput search for magnetic topological materials using spin-orbit spillage, machine learning, and experiments

Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation. Here, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from the ≈ 40000 three-dimensional materials in the JARVIS-DFT database. First, we screen materials with net magnetic moment > 0.5 μB and spin-orbit spillage > 0.25, resulting in 25 insulating and 564 metallic candidates. The spillage acts as a signature of spin-orbit induced band-inversion. Then, we carry out calculations of Wannier charge centers, Chern numbers, anomalous Hall conductivities, surface bandstructures, and Fermi-surfaces to determine interesting topological characteristics of the screened compounds. We also train machine learning models for predicting the spillage, bandgaps, and magnetic moments of new compounds, to further accelerate the screening process. We experimentally synthesize and characterize a few candidate materials to support our theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

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

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2024

Nuclear Fuel and Pu Redox Studies from The Glenn T. Seaborg Institute at Idaho National Laboratory

The Glenn T. Seaborg Institute at Idaho National Laboratory (INL-GTSI) focuses on advancing fundamental research in the actinide sciences by providing unique opportunities to early career scientists and engineers to gain experience studying the actinide elements and their associated systems. The INL-GTSI is built from three focus areas that are based on the expertise and supporting infrastructure at INL and include solid state chemistry and physics, solution phase chemistry and physics, and forensic and isotope science. INL is the lead Laboratory for nuclear energy research and development in the U. S. and the research on nuclear fuels performed under the INL-GTSI gives good examples of solid state studies. Uranium-Molybdenum (U-Mo) alloys are leading fuel candidates for conversion of high performance research and test reactors to low-enriched fuels. During irradiation, generated fission gas accumulates into bubbles and self-organizes into a gas bubble superlattice (GBS) that effectively stores fission gases and inhibits fuel swelling. A study on the early self- organizing behavior of the GBS shows that not only grain boundaries but the interfaces between the U-Mo matrix and uranium carbide (UC) impurities are important to GBS formation.[1] In solution, understanding the complex redox behavior of plutonium in aqueous environments is critical for establishing optimized nuclear waste reprocessing solvent systems and storage tank environments. INL-GTSI researchers have produced an experimentally validated multi-scale model of the gamma radiation induced behavior of plutonium ions in concentrated aqueous HNO3 solutions.[2] Here, gamma radiation effected only minimal steady state changes in the redox distribution of the plutonium oxidation states. The redox cycling between Pu(IV) and Pu(III) is demonstrated to be mediated by the •OH/NO3• radical oxidation of Pu(III) and the H2O2/HNO3 driven reduction of Pu(IV). The INL-GTSI offers young researchers the unique chance to work directly with actinide bearing materials in a U. S. National Laboratory environment. Further topical areas of interest to the INL-GTSI include, but are not limited to, fundamental actinide properties, structure/property (electronic, magnetic, thermal) relations, actinide quantum criticality, f- electron interactions, electron correlations, computational studies, new phases, defect effects, interface interactions, isotope production and separation, forensic analytical chemistry, structure and dynamic properties of actinides in non-aqueous media, separations chemistry and kinetics for advanced nuclear fuel cycles, radiation effects, and innovative and advanced ligand design for complexation of the actinides.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Design of Phase-Change Thermal Storage Device in a Heat Pump for Building Electric Peak Load Shaving: Preprint

Replacing carbon-intensive fossil fuel heating systems with electric heat pumps powered by renewables is a promising approach to decarbonize the building sector. However, one of the technical barriers of this approach is the large scale of heat demand, which will put excessive stress on the electricity grid. Integrating thermal energy storage (TES) into the heating systems can help alleviate this problem, by shifting thermal load and thus shaving peaks in the building electric load. Therefore, it is critical to understand how to design a thermal storage device in a heat pump for peak load shaving. In this study, we developed a numerical model for a cascaded vapor compression heat pump system integrating a phase change thermal storage device. This novel system can control the net thermal charging and discharging rate of the TES independently from the building's thermal load, which allows for precise control of electric power use. In the current study, we controlled the system to shave building electric peak load during a cold winter morning, and to charge TES during the relatively warm afternoon while still providing space heating. We used the model to evaluate the system performance with and without peak shaving. We also investigated the effect of PCM transition Tt on the peak reduction and electric energy saving potentials. The results show that peak shaving scheme effectively reduces the peak electric power consumption during a predefined discharge time window. When comparing to the no shaving case, for PCM Tt = 10 degrees C, the peak electric load reduction is 23.5%. When comparing to an air-source heat pump with back up electric heater, as Tt increases from 0 to 20 degrees C, the peak reduction increases from 46.1% to 50.9%. Integrating PCM Tt = 10 degrees C with peak shaving leads to the 45.5% of electric energy saving, which is the highest among the three transition temperatures.

electric energy saving↗
DOE OSTI2021

Biophysical analysis of the plant-specific GIPC sphingolipids reveals multiple modes of membrane regulation

The plant plasma membrane (PM) is an essential barrier between the cell and the external environment, controlling signal perception and transmission. It consists of an asymmetrical lipid bilayer made up of three different lipid classes: sphingolipids, sterols, and phospholipids. The glycosyl inositol phosphoryl ceramides (GIPCs), representing up to 40% of total sphingolipids, are assumed to be almost exclusively in the outer leaflet of the PM. However, their biological role and properties are poorly defined. In this study, we investigated the role of GIPCs in membrane organization. Because GIPCs are not commercially available, we developed a protocol to extract and isolate GIPC-enriched fractions from eudicots (cauliflower and tobacco) and monocots (leek and rice). Lipidomic analysis confirmed the presence of trihydroxylated long chain bases and 2-hydroxylated very long-chain fatty acids up to 26 carbon atoms. The glycan head groups of the GIPCs from monocots and dicots were analyzed by gas chromatograph-mass spectrometry, revealing different sugar moieties. Multiple biophysics tools, namely Langmuir monolayer, ζ-Potential, light scattering, neutron reflectivity, solid state 2H-NMR, and molecular modeling, were used to investigate the physical properties of the GIPCs, as well as their interaction with free and conjugated phytosterols. We showed that GIPCs increase the thickness and electronegativity of model membranes, interact differentially with the different phytosterols species, and regulate the gel-to-fluid phase transition during temperature variations. These results unveil the multiple roles played by GIPCs in the plant PM.

59 BASIC BIOLOGICAL SCIENCES↗