Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional distributions”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Summary of DASSH-CFD Inter-Assembly Heat Transfer Comparison

Fast reactor core design requires integration of detailed power distributions obtained from neutronics models and thermal fluids calculations to determine temperature distributions and flow conditions. During the design phase, calculations may be repeated as specifications change and computational expense must be balanced with solution accuracy, creating an opportunity for the use of intermediate-fidelity codes. To that end, the Ducted Assembly Steady-State Heat Transfer Software (DASSH) has been developed at Argonne National Laboratory to perform full-core, steady state coolant and fuel pin temperature calculations in reactors consisting of hexagonal, ducted assemblies. DASSH is a first-order, finite volume subchannel code based on methodologies initially introduced in the legacy thermal hydraulics codes SUPERENERGY-2 and SE2-ANL. It can work coupled with the Argonne Reactor Computational (ARC) code suite to obtain neutron and gamma power distributions based on neutron transport solutions and features many improvements and new capabilities compared to its predecessors. DASSH calculates temperatures by balancing energy between subchannels and relies on correlations to account for flow distribution and coolant mixing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Which Combinations of Environmental Conditions and Microphysical Parameter Values Produce a Given Orographic Precipitation Distribution?

This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.

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

Development of high-performance roll-to-roll-coated gas-diffusion-electrode-based fuel cells

This study focuses on determining fabrication conditions to create high-performance roll-to-roll-coated (R2R-coated) gas-diffusion electrodes (GDEs) for proton-exchange-membrane fuel cells (PEMFCs). Here, we examine how process conditions influence the distribution of ionomer in the electrode, which is shown to be critical for high performance. Using a combination of Kelvin probe, X-ray photoelectron spectroscopy, and nano-scale X-ray computed tomography we show that formation of an ionomer-rich surface is promoted by using a higher drying rate. We show that R2R-coated GDEs have higher surface ionomer concentration than spray-coated GDEs, which enables these R2R-coated GDEs to not need an additional ionomer overlayer, as is typically the case for spray-coated GDEs. This will reduce the number of processing steps and lower material costs in a manufacturing setting. This work shows that with the appropriate selection of materials, ink formulation, and processing conditions, direct-coated GDEs are a viable pathway for fuel cell manufacturing.

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

A comparative analysis of residual stresses from friction stir processing of aluminum cast 380 and wrought 7075 alloy sheets: experimental characterization and modeling

Residual stresses are often overlooked in friction stir processing (FSP), but their significant impact on fatigue performance necessitates their consideration in optimizing processing parameters. The first step in this effort is understanding how process conditions influence residual stress distributions, especially across different alloys. This study focuses on determining and explaining the through-thickness residual stress variations and the effect of process temperature on the residual stress magnitude in wrought AA7075 and cast AA380.0 alloys. Additionally, for AA380.0, the impact of a second FSP pass was investigated. To achieve this, hole-drilling electronic speckle pattern interferometry (ESPI) and the thermal pseudo-mechanical (TPM) model within finite element analysis were employed to study the 3D distributions of in-plane residual stresses in processed samples under various conditions. A key finding was the varying impact of process temperatures on residual stress magnitudes. Higher process temperatures reduced stresses in AA380.0 but increased them in AA7075. Additionally, the through-thickness stress distributions differed between the two alloys. Further analysis revealed that yield stresses are crucial in explaining these phenomena and the effects of additional FSP passes. Further, this fundamental understanding will be vital in guiding the efforts to mitigate residual stresses and assess their impact on the performance of FSP aluminum alloys.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Electron Beam Welding of Pure Tungsten Hex Cans for Nuclear Thermal Propulsion Engines

Nuclear thermal propulsion (NTP) is an in-space propulsion method currently being developed at the NASA Marshall Space Flight Center (MSFC). NTP systems are a high specific impulse (750–1,100 s), high thrust (15,000–250,000 lbf ) method of propulsion which have the potential to allow for faster transit times when optimizing for high ΔV. In the nuclear rocket engine, the heat from the nuclear fission reaction is transferred to a low molecular mass propellant (such as hydrogen). Hot propellant is expanded through a nozzle to generate thrust. Development of ceramic metal (cermet) fuel systems for NTP applications is currently ongoing at MSFC. In cermet fuel systems, ceramic fissile fuel particles such as uranium nitride or uranium dioxide are dispersed within a net-shaped, high-density structural matrix. The composite material is cladded by a protective metal structure to make up an NTP fuel element. Cladding materials must be able to withstand the demanding operating conditions required of the engine as well as retain a hermetic seal to allow for retention of fuel element structural integrity, prevent hydrogen attack or migration of the ceramic fuel, and limit release of fission products during operation. For NTP applications, tungsten is a prime material for both the metal matrix and cladding in cermet fuel systems because of its high melting point, high temperature strength, and compatibility with hot hydrogen. If a weld in tungsten with the capability of holding a hermetic seal is achievable, tungsten becomes a strong candidate for NTP applications. This Technical Memorandum focuses on determining the weldability of pure tungsten using electron beam welding (EBW). Tungsten appears well suited for NTP applications, but it has a high ductile to brittle transition temperature (DBTT) dependent upon chemical composition, structure/stress distribution, and mechanical conditions. Therefore, it is highly subject to brittle fracture. Because of its high susceptibility to brittle fracture, it is very difficult to weld. EBW was chosen for joining pure tungsten because of its low heat input compared to gas tungsten arc welding. Reduced heat input can be directly correlated with an increase in ductility of a tungsten weld. EBW is a high energy density welding process in which a stream of electrons penetrates a weld joint in a deep, narrow spike in contrast to a broad gas tungsten arc weld pool. The investigation initially focused on EBW of tungsten plates of both 0.01 in and 0.03 in thickness to determine if EBW could weld pure tungsten without the presence of visual defects—particularly cracking—in the welds. Variation in the weld procedure and post-weld heat treatment (PWHT) was used to improve the surface appearance of flat EBWs on a pure tungsten sheet. The investigation moved on to weld 0.05-in-thick hexagonal tungsten cans with a weld joint thickness of 0.025 in. The goal for welding the pure tungsten hex cans was to avoid any visual surface defects and generate a weld capable of a hermetic seal. This proved difficult. Cold welds commonly exhibited porosity that leaked air. Hot welds exhibited cracks, typically observed immediately after welding. Later welds were preheated to increase ductility and decrease the likelihood of through-thickness cracking. PWHT was used to arrest microcrack growth both in the flat weld samples and hexagonal weld samples.

Courtright, Z. S.↗
DOE OSTI2017

Warm Ionized Medium throughout the Sagittarius–Carina Arm

Wisconsin H-Alpha Mapper observations of Hα and [S ii] λ6716 emission are used to trace the vertical distribution and physical conditions of the warm ionized medium along the Sagittarius–Carina arm. CO emission, tracing cold molecular gas in the plane of the Galaxy, is used as a guide to isolate Hα and [S ii] λ6716 emission along individual spiral arms. Exponential scale heights of electron density squared (or emission measure) are determined using Hα emission above (below) the midplane to be 330±80 pc (550±230 pc) along the near Sagittarius arm, 300±100 pc (250±30 pc) along the near Carina arm, and >1000 pc along the far Carina arm. The emission measure scale height tends to increase as a function of Galactocentric radius along the Sagittarius–Carina arm for R{sub G}>8 kpc. Physical conditions of the ionized gas are analyzed using the [S ii]/Hα line ratio, which more closely traces I {sub Hα} than height above the plane, z, suggesting a stronger relationship with the in situ electron density. We interpret this result as further evidence for the majority of the observed diffuse emission originating from in situ ionized gas as opposed to scattered light from classical H ii regions in the plane.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Power Hardware-in-the-Loop-Based Anti-Islanding Evaluation and Demonstration

The National Renewable Energy Laboratory (NREL) teamed with Southern California Edison (SCE), Clean Power Research (CPR), Quanta Technology (QT), and Electrical Distribution Design (EDD) to conduct a U.S. Department of Energy (DOE) and California Public Utility Commission (CPUC) California Solar Initiative (CSI)-funded research project investigating the impacts of integrating high-penetration levels of photovoltaics (PV) onto the California distribution grid. One topic researched in the context of high-penetration PV integration onto the distribution system is the ability of PV inverters to (1) detect islanding conditions (i.e., when the distribution system to which the PV inverter is connected becomes disconnected from the utility power connection) and (2) disconnect from the islanded system within the time specified in the performance specifications outlined in IEEE Standard 1547. This condition may cause damage to other connected equipment due to insufficient power quality (e.g., over-and under-voltages) and may also be a safety hazard to personnel that may be working on feeder sections to restore service. NREL teamed with the Florida State University (FSU) Center for Advanced Power Systems (CAPS) to investigate a new way of testing PV inverters for IEEE Standard 1547 unintentional islanding performance specifications using power hardware-in-loop (PHIL) laboratory testing techniques.

14 SOLAR ENERGY↗
DOE OSTI2016

A flux-splitting method for hyperbolic-equation system of magnetized electron fluids in quasi-neutral plasmas

A flux-splitting method is proposed for the hyperbolic-equation system (HES) of magnetized electron fluids in quasi-neutral plasmas. The numerical fluxes are split into four categories, which are computed by using an upwind method which incorporates a flux-vector splitting (FVS) and advection upstream splitting method (AUSM). The method is applied to a test calculation condition of uniformly distributed and angled magnetic lines of force. All of the pseudo-time advancement terms converge monotonically and the conservation laws are strictly satisfied in the steady state. The calculation results are compared with those computed by using the elliptic–parabolic-equation system (EPES) approach using a magnetic-field-aligned mesh (MFAM). Both qualitative and quantitative comparisons yield good agreements of results, indicating that the HES approach with the flux-splitting method attains a high computational accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Behavior of helium cavities in ion-irradiated W-Ni-Fe ductile-phase toughened tungsten

This study reports on the distribution of helium (He) cavities in a hot-rolled W-Ni-Fe ductile-phase toughened tungsten (DPT W) composite irradiated to a dose and a helium concentration that are comparable to those in the material after 5-year irradiation in a conceptual fusion power plant. Here, the DPT W sample consists of W particles embedded in a ductile-phase NiFeW matrix with a nominal composition of 90W-7Ni-3Fe by weight. It was hot-rolled to a thickness reduction by 87% (87R DPT W). Sequential irradiation of the material with 1.2 MeV Ni + ions to a fluence of 2.15×10 16 Ni + /cm 2 and 90 keV He + ions to 6.5×10 15 He + /cm 2 was performed at 973 K. Larger He cavities with a lower number density are observed in NiFeW than W. Helium cavities are aggregated preferentially along the NiFeW/W interphase boundary. This behavior is not observed along the W/W grain boundary under the same irradiation conditions. A similar distribution of He cavities with comparable sizes is observed in 87R DPT W irradiated with He + ions only, suggesting that the He cavity distribution is not significantly affected by pre-irradiation of Ni + ions at 973 K. The data from this study does not provide any evidence for formation of visible voids or Ni precipitates in W phase irradiated with Ni + ions only. Diffusion and clustering of He atoms and mobile vacancies during ion irradiation at 973 K are responsible for the formation and growth of the He cavities.

36 MATERIALS SCIENCE↗
DOE OSTI2016

CFD simulation of simultaneous monotonic cooling and surface heat transfer coefficient

The monotonic heating regime method for determination of thermal diffusivity is based on the analysis of an unsteady-state (stabilised) thermal process characterised by an independence of the space-time temperature distribution on initial conditions. At the first kind of the monotonic regime a sample of simple geometry is heated / cooled at constant ambient temperature. The determination of thermal diffusivity requires the determination rate of a temperature change and simultaneous determination of the first eigenvalue. According to a characteristic equation the first eigenvalue is a function of the Biot number defined by a surface heat transfer coefficient and thermal conductivity of an analysed material. Knowing the surface heat transfer coefficient and the first eigenvalue the thermal conductivity can be determined. The surface heat transport coefficient during the monotonic regime can be determined by the continuous measurement of long-wave radiation heat flow and the photoelectric measurement of the air refractive index gradient in a boundary layer. CFD simulation of the cooling process was carried out to analyse local convective and radiative heat transfer coefficients more in detail. Influence of ambient air flow was analysed. The obtained eigenvalues and corresponding surface heat transfer coefficient values enable to determine thermal conductivity of the analysed specimen together with its thermal diffusivity during a monotonic heating regime.

42 ENGINEERING↗
NASA NTRSDate not supplied

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗
DOE OSTI2025

Probing the Effect of Electrode Thermodynamics on Reaction Heterogeneity in Thick Battery Electrodes

Thick electrodes present a viable strategy for enhancing energy density and reducing manufacturing costs of lithium-ion batteries. However, reaction heterogeneity during cycling compromises their rate capability and cycle life. While this nonuniformity is commonly attributed to sluggish charge transport, it is demonstrated here that the thermodynamic properties of the electrode material play an equally critical role. Through combined X-ray fluorescence microscopy and absorption near-edge structure spectroscopy, reaction distributions in LiFePO 4 (LFP) and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC) thick electrodes with matched porosity and tortuosity are compared. LFP electrodes develop pronounced depth-oriented state-of-charge (SOC) gradients that worsen with increasing discharge rates, whereas NMC maintains much more uniform SOC distributions under such conditions. This difference originates from their distinct SOC dependence of equilibrium potentials and is quantifiable through a dimensionless “reaction uniformity” number. Intriguingly, LFP thick electrodes also exhibit lateral SOC variations that strengthen during slow discharge. In conclusion, the enhanced reaction uniformity in NMC correlates with better active material utilization and slower capacity fade than LFP, highlighting electrode thermodynamics as a key design consideration for thick electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2015

On-demand in situ generation of oxygen in a nanofluidic embedded planar microband electrochemical reactor

In situ generation of reagents and their subsequent use downstream presents new opportunities to amplify the utility of nanofluidic devices by exploiting the confined geometry to address mass transport limitations on reaction kinetics and efficiency. Oxygen, an inherently valuable reactant, can be produced from electrolysis of water, a process that can be conveniently integrated within a nanofluidic system. Here, we construct and characterize a nanofluidic device consisting of a planar microband electrode embedded within a nanochannel for in situ electrochemical generation and optical monitoring of O 2 . Fluorescein, a dye with a pH-sensitive emission intensity, was used to monitor the spatiotemporal characteristics of the oxidation of H 2 O, using the co-produced H + . Application of anodic potentials at the nanochannel-embedded electrode results in a decrease in fluorescence intensity, which reflects the decreasing solution pH. A combination of fluorescence intensity and chronoamperometric response was used to quantitatively determine proton generation, and the H + /O 2 stoichiometry was then used to determine the concentration of the O 2 in the channel. Comparison of the experimental results to finite element simulations validates the use of fluorescein emission intensity to spectroscopically determine the local oxygen concentration in the nanochannel. By varying the applied potential, spatially averaged O 2 concentrations ranging from 0.13 to 0.41 mM were generated. In conclusion, the results demonstrate a convenient route to in situ modulation of the dissolved O 2 level in a nanofluidic device and the use of an optical probe to monitor its spatial and temporal distribution under flow conditions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

A stochastic scan strategy for grain structure control in complex geometries using electron beam powder bed fusion

Spatial control of microstructure within a three-dimensional component has been a dream of materials scientists for centuries. However, limitations in traditional manufacturing processes prevent detailed control over the distribution of microstructures in a single part. Here, we demonstrate the ability to control grain structure and crystallographic texture during metal additive manufacturing for arbitrary cross-sections of a practical size, with profound implications for the design and optimization of next-generation products. The key to this advance is a new geometry agnostic scan path algorithm that manipulates the spatial distribution of solidification conditions. Utilizing a fundamental understanding of solidification dynamics and a model of the heat transfer during processing, we have designed this algorithm to manipulate the natural competition between epitaxial dendrite growth and grain nucleation. With this algorithm, we successfully controlled the grain structure of Ni-based superalloy IN718 in the shape of the Mona Lisa.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sustainable ammonia synthesis from nitrogen wet with sea water by single-step plasma catalysis

Ammonia synthesis at ambient conditions employing intermittent distributed green sources of energy and feedstocks is globally sought to replace the centralized Haber-Bosch (H-B) process operating at high temperature and pressure. We report herein for the first time an effective and sustainable ammonia synthesis pathway from N 2 wet with seawater vapor over spherical SiO 2 and M/SiO 2 (M: Ag, Cu, and Co) catalysts driven by non-thermal plasma (NTP). Experimental results indicate that the presence of a catalyst is required for ammonia production from seawater vapor and N 2 . The Co/SiO 2 catalyst delivered the highest ammonia synthesis rate (r NH 3 ) of 3.7 mmol g cat -1 h -1 and energy yield of 3.2 g NH 3 ∙kW -1 ∙h -1 at a relatively low input power of 2 W. The extraction of H atoms from H 2 O molecules plays an important role in the ammonia synthesis from seawater vapor. Finally, this work unfolds a novel platform for the subsequent optimization of sustainable ammonia production from endless resources such as seawater and N 2 through catalytic non-thermal plasma potentially powered by renewable sources.

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

Mixing behavior of a model cellulosic biomass slurry during settling and resuspension

Thorough mixing during biochemical deconstruction of biomass is crucial for achieving maximum process yields and economic success. However, due to the complex morphology and surface chemistry of biomass particles, biomass mixing is challenging and currently it is not well understood. This study investigates the bulk rheology of negatively buoyant, non-Brownian α-cellulose particles during settling and resuspension. The torque signal of a vane mixer across two distinct experimental setups (vane-in-cup and vane-in-beaker) was used to understand how mixing conditions affect the distribution of biomass particles. During experimentation, a bifurcated torque response as a function of vane speed was observed, indicating that the slurry transitions from a “settling-dominant” regime to a “suspension-dominant” regime. The torque response of well-characterized fluids (i.e., DI water) were then used to empirically identify when sufficient mixing turbulence was established in each experimental setup. The predicted critical mixing speeds were in agreement with measured values, suggesting that secondary flows are required in order to keep the cellulose particles fully suspended. In addition, a simple scaling relationship was developed to model the entire torque signal of the slurry throughout settling and resuspension. Furthermore, qualitative and semi-quantitative agreement between the model and experimental results was observed.

09 BIOMASS FUELS↗
DOE OSTI2017

Molybdenum isotope fractionation and speciation in a euxinic lake—Testing ways to discern isotope fractionation processes in a sulfidic setting

The molybdenum (Mo) isotope composition in euxinic shales has been used as a proxy for the global distribution of anoxic conditions in ancient oceans, and since more recently also as a proxy for sulfide concentrations in depositional environments. However, there is currently no way to distinguish isotope fractionation at low bottom water sulfide concentrations in ‘local’ basins from ‘global’ secular isotope variations associated with changing seawater composition. This uncertainty is challenging the use of Mo isotopes for paleoceanographic reconstructions. To explore this further, we present new data from sediments deposited over the past ~ 9800 years in one of the best studied euxinic localities in the world: Lake Cadagno in Switzerland. The sample set allows us to test ways to discern isotope fractionation processes at play in a highly restricted euxinic basin. Most of our drill core samples (n = 18) show high δ98Mo values similar to previously studied shallow sediments, indicative of quantitative Mo removal from the water column (Dahl et al. 2010a). However, a few samples (n = 3) deposited between about 1200 and 3400 years ago carry low δ98Mo values and have been isotopically fractionated in the lake. Sedimentological and geochemical characterizations show that these δ98Mo-fractionated sediments formed during times of frequent injection of O2- and sediment-rich river water into the deep sulfidic water column. A positive correlation between δ98Mo and sedimentary Mo contents suggests that isotope fractionation occurred during times of non-quantitative Mo removal, although Mn-oxide cycling at the chemocline might also contribute a subordinate proportion of (98Mo-depleted) molybdenum into the sulfidic zone. Sedimentary Mo/U enrichments relative to oxic lake water further supports the hypothesis that a particulate Mo shuttle was most efficient during times of quantitative Mo removal. Therefore, periods with inefficient Mo capture are ascribed to incomplete conversion of molybdate to particle reactive Mo species when bottom water H2S levels were low or less stable than today. Using XAFS spectroscopy, we found that the two distinct Mo compounds predominating in the sediments (MoIV-S and MoVI-OS) are not diagnostic for isotope fractionation that has occurred in Lake Cadagno. Instead, we infer that δ98Mo-fractionated products (forming via a low-sulfide Mo pathway) can be subsequently altered with little or no isotopic imprint during remobilization and re-precipitation (e.g., at higher sulfide levels in the sediments) as well as during post-depositional oxidation. Future work could investigate local δ98Mo-fractionation processes expressed in other euxinic settings and explore other sedimentary metrics to constrain the steps involved in the euxinic burial pathway(s). One tantalizing prospect of this is to distinguish between local bottomwater sulfide levels and variations in the fraction of global seafloor anoxia from the Mo isotope composition in ancient euxinic mudrocks.

58 GEOSCIENCES↗
DOE OSTI2023

Domain-aware Control-oriented Neural Models for Autonomous Underwater Vehicles

Conventional physics-based modeling is a time-consuming bottleneck in control design for complex nonlinear systems like autonomous underwater vehicles (AUVs). In contrast, purely data-driven models, require a large number of observations and lack operational guarantees for safety-critical systems. Data-driven models leveraging available partially characterized dynamics have potential to provide reliable systems models in a typical data-limited scenario for high value complex systems, thereby avoiding months of expensive expert modeling time. In this work we explore this middle-ground between expert-modeled and pure data-driven modeling. We present control-oriented parametric models with varying levels of domain-awareness that exploit known system structure and prior physics knowledge to create constrained deep neural dynamical system models. We employ universal differential equations to construct data-driven blackbox and graybox representations of the AUV dynamics. In addition, we explore a hybrid formulation that explicitly models the residual error related to imperfect graybox models. We compare the prediction performance of the learned models for different distributions of initial conditions and control inputs to assess their suitability for control.

Shaw Cortez, Wenceslao E.↗