SEARCH · Engineering Papers

Results for “constrained control”

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.

SourceResearch connectionFrom yearThrough year

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

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Particle phase-state variability in the North Atlantic free troposphere during summertime is determined by atmospheric transport patterns and sources

Abstract. Free tropospheric aerosol particles have important but poorly constrained climate effects due to transformations of their physicochemical properties during long-range transport. In this study, we investigate the chemical composition and provide an overview of the phase states of individual particles that have undergone long-range transport over the North Atlantic Ocean in June and July 2014, 2015, and 2017 to the Observatory of Mount Pico (OMP) in the Azores. The OMP is an ideal site for studying long-range-transported free tropospheric particles because local emissions have a negligible influence and contributions from the boundary layer are rare. We used the FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) to determine the origins and transport trajectories of sampled air masses and found that most of them originated from North America and recirculated over the North Atlantic Ocean. The FLEXPART analysis showed that the sampled air masses were highly aged (average plume age >10 d). Size-resolved chemical compositions of individual particles were probed using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX) and scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy (STXM-NEXAFS). CCSEM-EDX results showed that the most abundant particle types were carbonaceous (∼ 29.9 % to 82.0 %), sea salt (∼ 0.3 % to 31.6 %), and sea salt with sulfate (∼ 2.4 % to 31.5 %). We used a tilted stage interfaced within an environmental scanning electron microscope (ESEM) to determine the phase states of individual submicron particles. We found that most particles (∼ 47 % to 99 %) were in the liquid state at the time of collection due to inorganic inclusions. Moreover, we also observed substantial fractions of solid and semisolid particles (∼ 0 % to 30 % and ∼ 1 % to 42 %, respectively) during different transport patterns and events, reflecting the particles' phase-state variability for different atmospheric transport events and sources. Combining phase state measurements with FLEXPART CO tracer analysis, we found that wildfire-influenced plumes can result in particles with a wide range of viscosities after long-range transport in the free troposphere. We also used temperature and RH values extracted from the Global Forecast System (GFS) along the FLEXPART-simulated path to predict the phase state of the particles during transport and found that neglecting internal mixing with inorganics would lead to an overestimation of the viscosity of free tropospheric particles. Our findings warrant future investigation aiming at the quantitative assessment of the influence of internal mixing on the phase states of the individual particles. This study also provides insights into the chemical composition and phase state of free tropospheric particles, which can help models to reduce uncertainties about the effects of ambient aerosol particles on climate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Temperature dependence of uranium and thorium partitioning in igneous zircons

Zircon is a key mineral in geochronology because of its chemical and physical durability and tendency to incorporate radioactive trace elements such as U and Th. Quantifying the partitioning of the actinide elements is critical to constrain initial non-secular equilibrium amounts of 230 Th in zircon. An excess or deficit of 206 Pb will be produced from such an initial excess/deficit of 230 Th from the secular equilibrium condition, which influences the calculated 206 Pb/ 238 U age (Schärer, 1984, Mattinson, 1973). However, there is no standard way to calculate Th/U partitioning ratios when applying age corrections to young igneous zircon, making uncertainties hard to estimate. To better understand U and Th partitioning of zircon/melt that is the cornerstone of the secular disequilibrium correction, zircon was synthesized in one-atmosphere experiments using basaltic andesite, andesite, and rhyolite starting materials, doped with Zr, U, and Th. Different experimental temperatures and oxygen fugacity conditions (ΔQFM−4 to ΔQFM+4) were explored to examine their effects on U and Th partitioning. In addition, we specifically quantify the effects of sector zoning, fractional crystallization, and melt composition on U and Th partitioning. By combining experimental and natural zircon data, we find that temperature has the primary control on the partitioning of U and Th in the zircon and calibrate an inverse relationship between these partition coefficients and crystallization temperatures. The calibrated equation can be applied to the 230 Th correction for an improvement in the accuracy of Th-corrected 206 Pb/ 238 U dates when the zircon crystallization temperature is known.

206Pb/238U dating↗

DOE OSTI2020

Extended use of face masks during the COVID-19 pandemic - Thermal conditioning and spray-on surface disinfection

The current COVID-19 pandemic has resulted in globally constrained supplies for face masks and personal protective equipment (PPE). Production capacity is limited in many countries and the future course of the pandemic will likely continue with shortages for high quality masks and PPE in the foreseeable future. Therefore, expectations are that mask reuse, extended wear and similar approaches will enhance the availability of personal protective measures. Repeated thermal disinfection could be an important option and likely easier implemented in some situations, at least on the small scale, than UV illumination, irradiation or hydrogen peroxide vapor exposure. In this paper, an overview on thermal responses and ongoing filtration performance of multiple face mask types is provided. Most masks have adequate material properties to survive a few cycles (i.e. 30 min disinfection steps) of thermal exposure in the 75°C regime. Some are more easily affected, as seen by the fusing of plastic liner or warping, given that preferred conditioning temperatures are near the softening point for some of the plastics and fibers used in these masks. Hence adequate temperature control is equally important. As guidance, disinfectants sprayed via dilute solutions maintain a surface presence over extended time at 25 and 37°C. Some spray-on alcohol-based solutions containing disinfectants were gently applied to the top surface of masks. Neither moderate thermal aging (less than 24 h at 80 and 95°C) nor gentle application of surface disinfectant sprays resulted in measurable loss of mask filter performance. Subject to bio-medical concurrence (additional checks for virus kill efficiency) and the use of low risk non-toxic disinfectants, such strategies, either individually or combined, by offering additional anti-viral properties or short term refreshing, may complement reuse options of professional masks or the now ubiquitous custom-made face masks with their often unknown filtration effectiveness.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Latest Three-Flavor Neutrino Oscillation Results from NOvA

NOvA, is a two-detector, long-baseline neutrino oscillation experiment located at Fermilab, Batavia, IL, USA. It aims to constrain neutrino oscillation parameters by analyzing $\nu_\mu (\bar{\nu}_\mu)$ disappearance and $\nu_e (\bar{\nu}_e)$ appearance data. The experiment uses the Neutrinos at Main Injector (NuMI) beamline at Fermilab, which delivers a high-purity 900 KW beam of neutrinos and anti-neutrinos. The detectors are functionally identical finely granulated liquid tracking calorimeters, both situated 14.6 mrad off-axis to the beam direction. The NOvA Near Detector (ND), situated 100 meters underground and 1 kilometer from the beam source, detects the un-oscillated $\nu_\mu (\bar{\nu}_\mu)$ and beam $\nu_e (\bar{\nu}_e)$ events. The Far Detector (FD), located in Ash River, MN, USA, 809 kilometers from the ND, records the oscillated $\nu_e (\bar{\nu}_e)$ and the un-oscillated $\nu_\mu (\bar{\nu}_\mu)$ events. NOvA employs an extrapolation technique to predict the expected events at the Far Detector based on the Near Detector data, thereby providing a significant constraint on systematic uncertainties in the oscillation analyses. As NOvA accumulates more data, controlling these systematic uncertainties becomes increasingly important. This talk will detail the NOvA neutrino oscillation analysis framework and its approach to minimizing dominant systematic uncertainties using Near Detector data. The latest three flavor neutrino oscillation results based on a neutrino-beam exposure of $26.60 \times 10^{20}$ POT and an anti-neutrino beam exposure of $12.50\times 10^{20}$ POT and a novel low energy $\nu_e$ sample, will also be presented.

Choudhary, Brajesh↗

DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

DOE OSTI2021

Probing Quantum Speed Limits with Ultracold Gases

Here, quantum speed limits (QSLs) rule the minimum time for a quantum state to evolve into a distinguishable state in an arbitrary physical process. These fundamental results constrain a notion of distance traveled by the quantum state, known as the Bures angle, in terms of the speed of evolution set by nonadiabatic energy fluctuations. I theoretically propose how to measure QSLs in an ultracold quantum gas confined in a time-dependent harmonic trap. In this highly-dimensional system of continuous variables, quantum tomography is prohibited. Yet, QSLs can be probed whenever the dynamics is self-similar by measuring as a function of time the cloud size of the ultracold gas. This makes it possible to determine the Bures angle and energy fluctuations, as I discuss for various ultracold atomic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗

DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗

DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part II: Modeling Study

Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.

Meteorology & Atmospheric Sciences↗

DOE OSTI2023

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗

DOE OSTI2022

Inferring and evaluating satellite-based constraints on NO x emissions estimates in air quality simulations

Satellite observations of tropospheric NO 2 columns can provide top-down observational constraints on emissions estimates of nitrogen oxides (NO x ). Mass-balance-based methods are often applied for this purpose but do not isolate near-surface emissions from those aloft, such as lightning emissions. Here, we introduce an inverse modeling framework that couples satellite chemical data assimilation to a chemical transport model. In the framework, satellite-constrained emissions totals are inferred using model simulations with and without data assimilation in the iterative finite-difference mass-balance method. The approach improves the finite-difference mass-balance inversion by isolating the near-surface emissions increment. We apply the framework to separately estimate lightning and anthropogenic NO x emissions over the Northern Hemisphere for 2019. Using overlapping observations from the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI), we compare separate NO x emissions inferences from these satellite instruments, as well as the impacts of emissions changes on modeled NO 2 and O 3 . OMI inferences of anthropogenic emissions consistently lead to larger emissions than TROPOMI inferences, attributed to a low bias in TROPOMI NO 2 retrievals. Updated lightning NO x emissions from either satellite improve the chemical transport model's low tropospheric O 3 bias. The combined lighting and anthropogenic emissions updates improve the model's ability to reproduce measured ozone by adjusting natural, long-range, and local pollution contributions. Thus, the framework informs and supports the design of domestic and international control strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Assessing the Behavior of Moderately Volatile Elements on the Moon in Order to Constrain Processes of Magmatic Evolution and Planetary Accretion

Heterogeneous isotopic compositions of Zn and Ga in lunar basalts and highlands rocks indicate that secondary processing affected the stable isotopic composition of major lunar reservoirs. It follows that these fractionating processes must be identified and characterized in order to make robust estimates for the isotopic composition of the mare basalt source regions or the bulk silicate Moon (BSM). Contrary to previous assertions (e.g., Kato & Moynier, 2017), MVEs such as Zn and Ga do not behave similarly on the Moon, in line with differences in their chemical behavior and volatility. There is clear evidence that impact reworking fractionates zinc isotopes but, in addition, zinc isotopes are also highly sensitive to the formation of Ti-bearing oxides such as ilmenite. In contrast, gallium is far less susceptible to volatile reworking, but its isotopic composition is fractionated during incorporation into plagioclase (Wimpenny et al., 2022; Render et al., 2023). Based on these findings we conclude that current isotopic compositions of Zn and Ga in lunar rocks, such as those measured for the mare basalt suite, are unlikely to be representative of the bulk Moon when it accreted. As such, estimates for the MVE isotopic composition of the BSM derived from mare basalts may not be accurate, which may have implications for the incorporation of these data into dynamical models of the Moon’s formation. In the following, we discuss the main findings from this project in more detail, focusing on the behaviour of zinc and gallium isotopic systematics in lunar rocks and the processes that control their isotopic fractionation on the Moon.

58 GEOSCIENCES↗

DOE OSTI2023

Scalable Bayesian optimization with randomized prior networks

Several fundamental problems in science and engineering consist of global optimization tasks involving unknown high-dimensional (black-box) functions that map a set of controllable variables to the outcomes of an expensive experiment. Bayesian Optimization (BO) techniques are known to be effective in tackling global optimization problems using a relatively small number objective function evaluations, but their performance suffers when dealing with high-dimensional outputs. To overcome the major challenge of dimensionality, here we propose a deep learning framework for BO and sequential decision making based on bootstrapped ensembles of neural architectures with randomized priors. Using appropriate architecture choices, we show that the proposed framework can approximate functional relationships between design variables and quantities of interest, even in cases where the latter take values in high-dimensional vector spaces or even infinite-dimensional function spaces. In the context of BO, we augmented the proposed probabilistic surrogates with re-parameterized Monte Carlo approximations of multiple-point (parallel) acquisition functions, as well as methodological extensions for accommodating black-box constraints and multi-fidelity information sources. We test the proposed framework against state-of-the-art methods for BO and demonstrate superior performance across several challenging tasks with high-dimensional outputs, including a constrained multi-fidelity optimization task involving shape optimization of rotor blades in turbo-machinery.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗

DOE OSTI2020

Landscape responses to intraplate deformation in the Kalahari constrained by sediment provenance and chronology in the Okavango Basin

The structural depression that occupies the Okavango Basin in southern Africa comprises a depo-centre within the intracratonic Kalahari Basin where sediments of the Cenozoic Kalahari Group have accumulated. The Okavango Basin has been formed due to stretching and subsidence at an area of diffused deformation, southwestwards to the main East African Rift System (EARS). Sediments from two full Kalahari Group sequences, located on opposite sides of the Gumare Fault that forms a major fault within the Okavango Basin, were studied to determine their provenance and chronology. Terrestrial Cosmogenic Nuclide (TCN) 26 Al/ 10 Be burial dating was used to constrain a chronostratigraphical framework, and Pb, Sr, and Nd isotopic ratios combined with geochemical and sedimentological analyses were applied to track the source areas of the sediments.Results indicate the following sequence of basin filling: (a) Accumulation between ca. 4–3 Ma during which the currently downthrown (southern) block received a mixture of sediments mostly from the Choma-Kalomo, Ghanzi-Chobe, and Damara terranes, and possibly from the Lufilian Belt and/or Karoo basalts during earlier stages of deposition. Simultaneously, the upthrown (northern) block received sediments from more distant Archean sources in the Zimbabwe and/or Kasai cratons, (b) Hiatus in sedimentation occurred at both sites between ca. 3–2 Ma, (c) Sediments on both sides of the Gumare Fault share a similar source (Angolan Shield) with minor distinct contributions to the downthrown block from the Kasai Craton and local sources input to the upthrown block, and (d) Regional distribution of aeolian sand since at least 1 Ma. The change in source areas is attributed to rearrangements of the drainage systems that were probably linked to vertical crustal movements on the margins of the Okavango Basin. The tectonically induced morphodynamics controlled the landscape evolution of the endorheic basin where vast lakes, wetlands and salt pans have developed through time.

58 GEOSCIENCES↗

DOE OSTI2026

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗