Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Visualizing Flutter Mechanism as Traveling Wave Through Animation of Simulation Results for the Semi-Span Super-Sonic Transport Wind-Tunnel Model

It has long been recognized that frequency and phasing of structural modes in the presence of airflow play a fundamental role in the occurrence of flutter. Animation of simulation results for the long, slender Semi-Span Super-Sonic Transport (S4T) wind-tunnel model demonstrates that, for the case of mass-ballasted nacelles, the flutter mode can be described as a traveling wave propagating downstream. Such a characterization provides certain insights, such as (1) describing the means by which energy is transferred from the airflow to the structure, (2) identifying airspeed as an upper limit for speed of wave propagation, (3) providing an interpretation for a companion mode that coalesces in frequency with the flutter mode but becomes very well damped, (4) providing an explanation for bursts of response to uniform turbulence, and (5) providing an explanation for loss of low frequency (lead) phase margin with increases in dynamic pressure (at constant Mach number) for feedback systems that use sensors located upstream from active control surfaces. Results from simulation animation, simplified modeling, and wind-tunnel testing are presented for comparison. The simulation animation was generated using double time-integration in Simulink of vertical accelerometer signals distributed over wing and fuselage, along with time histories for actuated control surfaces. Crossing points for a zero-elevation reference plane were tracked along a network of lines connecting the accelerometer locations. Accelerometer signals were used in preference to modal displacement state variables in anticipation that the technique could be used to animate motion of the actual wind-tunnel model using data acquired during testing. Double integration of wind-tunnel accelerometer signals introduced severe drift even with removal of both position and rate biases such that the technique does not currently work. Using wind-tunnel data to drive a Kalman filter based upon fitting coefficients to analytical mode shapes might provide a better means to animate the wind tunnel data.

Christhilf, David M.↗
DOE OSTI2022

Frictional properties of Opalinus Clay: influence of humidity, normal stress and grain size on frictional stability

SUMMARY The Opalinus Clay (OPA) is a clay-rich formation considered as a potential host rock for radioactive waste repositories and as a caprock for carbon storage in Switzerland. Its very low permeability (10−19 to 10−21 m2) makes it a potential sealing horizon, however the presence of faults that may be activated during the lifetime of a repository project can compromise the long-term hydrological confinement, and lead to mechanical instability. Here, we have performed laboratory experiments to test the effect of relative humidity (RH), grain size (g.s.) and normal stress on rate-and-state frictional properties and stability of fault laboratory analogues corresponding to powders of OPA shaly facies. The sifted host rock powders at different grain size fractions (<63 μm and 63 < g.s. < 125 μm), at room (∼25 per cent) and 100 per cent humidity, were slid in double-direct shear configuration, under different normal stresses (5–70 MPa). We observe that peak friction, μpeak and steady-state friction, μss, depend on water vapour content and applied normal stress. Increasing relative humidity from ∼25 per cent RH (room humidity) to 100 per cent RH causes a decrease of frictional coefficient from 0.41 to 0.35. The analysis of velocity-steps in the light of rate-and-state friction framework shows that the stability parameter (a–b) is always positive (velocity-strengthening), and it increases with increasing sliding velocity and humidity. The dependence of (a–b) on slip rate is lost as normal stress increases, for each humidity condition. By monitoring the variations of the layer thickness during the velocity steps, we observe that dilation (Δh) is directly proportional to the sliding velocity, decreases with normal stress and is unaffected by humidity. Microstructural analysis shows that most of the deformation is accommodated within B-shear zones, and the increase of normal stress (σn) promotes the transition from strain localization and grain size reduction to distributed deformation on a well-developed phyllosilicate network. These results suggest that: (1) the progressive loss of velocity dependence of frictional stability parameter (a–b) at σn > 35 MPa is dictated by a transition from localized to distributed deformation and (2) water vapour content does not affect the deformation mechanisms and dilation, whereas it decreases steady-state friction (μss), and enhances fault stability.

58 GEOSCIENCES↗
DOE OSTI2023

ANN-based ground motion model for Turkey using stochastic simulation of earthquakes

SUMMARY Turkey is characterized by a high level of seismic activity attributed to its complex tectonic structure. The country has a dense network to record earthquake ground motions; however, to study previous earthquakes and to account for potential future ones, ground motion simulations are required. Ground motion simulation techniques offer an alternative means of generating region-specific time-series data for locations with limited seismic networks or regions with seismic data gaps, facilitating the study of potential catastrophic earthquakes. In this research, a local ground motion model (GMM) for Turkey is developed using region-specific simulated records, thus constructing a homogeneous data set. The simulations employ the stochastic finite-fault approach and utilize validated input-model parameters in distinct regions, namely Afyon, Erzincan, Duzce, Istanbul and Van. To overcome the limitations of linear regression-based models, artificial neural network is used to establish the form of equations and coefficients. The predictive input parameters encompass fault mechanism (FM), focal depth (FD), moment magnitude (Mw), Joyner and Boore distance (RJB) and average shear wave velocity in the top 30 m (Vs30). The data set comprises 7359 records with Mw ranging between 5.0 and 7.5 and RJB ranging from 0 to 272 km. The results are presented in terms of spectral ordinates within the period range of 0.03–2.0 s, as well as peak ground acceleration and peak ground velocity. The quantification of the GMM uncertainty is achieved through the analysis of residuals, enabling insights into inter- and intra-event uncertainties. The simulation results and the effectiveness of the model are verified by comparing the predicted values of ground motion parameters with the observed values recorded during previous events in the region. The results demonstrate the efficacy of the proposed model in simulating physical phenomena.

Karimzadeh, Shaghayegh (ORCID:0000000337531676)↗
DOE OSTI2020

Additive Manufacturing of Energy Harvesting Material System for Active Wireless MEMS Sensors

Combustion systems make up a large sum of the world’s yearly energy production despite advances in renewable energy sources. Precise control of pressure and temperature in combustion chambers allows for lower carbon emissions, and overall higher efficiency. Wireless temperature and pressure sensors are critical for the operation of combustion chambers. Energy harvesters have emerged as a viable way to power wireless sensors when heat and mechanical energy are available. Lithium niobate (LNB), barium titanate (BTO), and lead zirconate titanate (PZT) pyroelectric ceramics were studied and fabricated through additive manufacturing (AM) for use as energy harvesting structures in combustion environments. The electric power generation potential of LNB and PZT ceramics was studied by varying the temperature from 50 to 60 °C and by introducing cyclic compression loads of 2000 N amplitude. It was found that thermal energy conversion had the highest output of 500 nW and that combined thermal and mechanical conversion did not increase the harvesting potential because of the competing contributions of both effects. Powder bed and slurry extrusion AM methods were used to fabricate BTO and PZT ceramic structures. Optimization of the powder-based binder jetting method produced BTO ceramics with a relative density of 36.77% a piezoelectric coefficient of 153 pC/N. The density of the ceramics increased up to 56% when increasing the saturation of binder in the powder to because of liquid phase sintering. Finally, PZT structures were manufactured through AM and fin features were added to the design to enhance the heat transfer along the material. The harvested power density of the flat samples was 3.643 μW/cm 3 and of the finned samples was 3.034 μW/cm3. This research paves the ways for the development of self-powered wireless sensors in critical areas of operation such as combustion.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Building a machine learning surrogate model for wildfire activities within a global Earth system model

Abstract. Wildfire is an important ecosystem process, influencing land biogeophysical and biogeochemical dynamics and atmospheric composition. Fire-driven loss of vegetation cover, for example, directly modifies the surface energy budget as a consequence of changing albedo, surface roughness, and partitioning of sensible and latent heat fluxes. Carbon dioxide and methane emitted by fires contribute to a positive atmospheric forcing, whereas emissions of carbonaceous aerosols may contribute to surface cooling. Process-based modeling of wildfires in Earth system land models is challenging due to limited understanding of human, climate, and ecosystem controls on fire counts, fire size, and burned area. Integration of mechanistic wildfire models within Earth system models requires careful parameter calibration, which is computationally expensive and subject to equifinality. To explore alternative approaches, we present a deep neural network (DNN) scheme that surrogates the process-based wildfire model with the Energy Exascale Earth System Model (E3SM) interface. The DNN wildfire model accurately simulates observed burned area with over 90 % higher accuracy with a large reduction in parameterization time compared with the current process-based wildfire model. The surrogate wildfire model successfully captured the observed monthly regional burned area during validation period 2011 to 2015 (coefficient of determination, R2=0.93). Since the DNN wildfire model has the same input and output requirements as the E3SM process-based wildfire model, our results demonstrate the applicability of machine learning for high accuracy and efficient large-scale land model development and predictions.

58 GEOSCIENCES↗
DOE OSTI2019

Adaptive harnessing damping in hydrokinetic energy conversion by two rough tandem-cylinders using flow-induced vibrations

Alternating lift technology (ALT) harnesses hydrokinetic energy from currents and tides and is different from conventional steady lift technologies (turbines). ALT is still under development with an important research issue being the mitigation of wake effect on downstream cylinders in Flow Induced Oscillations (FIO) of multi-cylinder Current Energy Converters (CEC). Rather than adjust the configuration/parameters of the converter, a more direct and effective way is to actively adjust the harnessed power. In this paper, a hydrokinetic energy converter using FIO of two cylinders in tandem in one-degree-of-freedom oscillators, with nonlinear adaptive damping and linear spring stiffness, is tested experimentally. FIO include Vortex Induced Vibrations (VIV), galloping, and their coexistence. Introducing adaptive damping into the tandem cylinders overcomes the shielding effect, which previous research on two tandem cylinders has shown to reduce power output by the downstream cylinder. Shielding along with fixed damping result either in ceasing motion due to excessive damping, or in low harnessed energy due to insufficient damping. In this experimental study, damping-to-velocity rate, linear spring-stiffness, cylinder spacing, and flow-velocity are the parameters with Reynolds number 30,000 ≤ Re ≤ 120,000. Comparison to linear-oscillators in FIO shows that this nonlinear converter, with velocity-proportional damping coefficient, is more effective over the entire FIO range but especially in galloping, where both flow and cylinder speeds are higher. Experimental results for energy harvesting, efficiency and instantaneous energy of the converter are presented and discussed supported by amplitude and frequency response data. Here the results show that the nonlinear, adaptive, velocity-proportional damping coefficient is an effective way to increase the overall harnessed power and the power of downstream cylinder. The harnessed power in the VIV to galloping transition increases by up to 94%. The most significant improvement is in the galloping region, where the increase in harnessed power and efficiency is around 33%.

13 HYDRO ENERGY↗
NASA NTRS2020

Model Validation for Bigelow Expandable Activity Module (BEAM) with Stowage

The Bigelow Expandable Activity Module (BEAM) was a collaboration between NASA and Bigelow Aerospace to develop an expandable habitat technology that could be used for future space exploration missions such as the Artemis program. On April 26, 2016, the BEAM was launched and berthed to the International Space Station (ISS) on the aft port of Node 3. After the originally planned two-year mission life, BEAM was approved for use as a stowage module for a 109-cargotransfer bag equivalent (CTBE) to alleviate stowage limitations onboard the ISS. The present work builds on the previous 2017 model validation to flight data that demonstrated that BEAM could generally accommodate cargo stowage without adversely affecting BEAM thermal control, primarily for condensation management. Once the stowage was configured in early 2019, model validation efforts resumed. The thermal model revisions included the CTBE stowage and updated CFD-derived heat transfer coefficients for various intra-module ventilation (IMV) flow rates. The coefficients were provided by the Crew and Thermal Systems Division (EC) at NASA JSC and Jacobs Technology. The model is now revalidated and able to predict minimum temperatures with good agreement to flight temperatures (within 1 °C) for both historic and recurring minimum values. Additionally, the model is partially validated for reduced flow rates. With this model, the team was able to provide new flight rule recommendations for condensation and IMV operations management.

Sydney Taylor↗
DOE OSTI2025

EVALUATION OF OPERATION TEMPERATURES UNDER NATURAL CONVECTION HELIUM FLOW IN A CONFINED CAVITY

The Spallation Neutron Source (SNS) is a high-power accelerator-based pulsed neutron source led by Oak Ridge National Laboratory (ORNL) to achieve high fluxes of neutrons for scientific experiments. Active and passive cooling of the systems and parts forming the SNS have been considered to warrant the safe operation of the facility. The diverse cooling systems make use of conjugated heat transfer mechanisms to provide a stable operation temperature for all components in the machine. Thermal power deposited into stainless-steel piping lines due to particle radiation may reach values of up to 1.2 W/cc in the regions located closer to the center of the lower IRP. These energy deposition levels, in not actively cooled components, such as the transfer-line-outer-vacuum-layer may increase the temperature of the component beyond design requirement limitations. The evaluation of the operation temperatures for the former components relies in the assumption that a low-pressure helium atmosphere provides enough heat removal capacity based on natural convection phenomena. In this work, the evaluation of steady state temperatures in components such as the CMS transfer lines has been evaluated using computational fluid dynamics (CFD), analytical correlations and experimental measurements. The companion experiments were conducted in a closed helium system at pressures varying from 1.1 to 1.5 bar. A copper rod was affixed horizontally between viewing windows and heated at constant power, and measurements were made of both the rod temperature and ambient temperature via a system fiberoptic distributed temperature sensors and RTDs. It was found that the measured heat transfer coefficients agree well with the predictions of Churchill and Chu correlations across the range of cases considered. Additionally, the ambient helium volume above the rod was imaged via background oriented schlieren (BOS), and these data was used to determine the line-averaged density gradients in this region. These gradients were compared to simulation data to validate the predictions of natural convection simulations.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗
DOE OSTI2016

Magma intrusion near Volcan Tancitaro: Evidence from seismic analysis

Between May and June 2006, an earthquake swarm occurred near Volcan Tancítaro in Mexico, which was recorded by a temporary seismic deployment known as the MARS network. We located ~1000 events from this seismic swarm. Previous earthquake swarms in the area were reported in the years 1997, 1999 and 2000. We relocate and analyze the evolution and properties of the 2006 earthquake swarm, employing a waveform cross-correlation-based phase repicking technique. Hypocenters from 911 events were located and divided into eighteen families having a correlation coefficient at or above 0.75. 90% of the earthquakes provide at least sixteen phase picks. We used the single-event location code Hypo71 and the P-wave velocity model used by the Jalisco Seismic and Accelerometer Network to improve hypocenters based on the correlation-adjusted phase arrival times. We relocated 121 earthquakes, which show clearly two clusters, between 9–10 km and 3–4 km depth respectively. The average location error estimates are <1 km epicentrally, and <2 km in depth, for the largest event in each cluster. Depths of seismicity migrate upward from 16 to 3.5 km and exhibit a NE-SW trend. The swarm first migrated toward Paricutin Volcano but by mid-June began propagating back toward Volcán Tancítaro. In addition to its persistence, noteworthy aspects of this swarm include a quasi-exponential increase in the rate of activity within the first 15 days; a b-value of 1.47; a jug-shaped hypocenter distribution; a shoaling rate of ~5 km/month within the deeper cluster, and a composite focal mechanism solution indicating largely reverse faulting. As a result, these features of the swarm suggest a magmatic source elevating the crustal strain beneath Volcan Tancítaro.

58 GEOSCIENCES↗
DOE OSTI2017

Monovalent versus Divalent Cation Diffusion in Thiospinel Ti2S4

Diffusion coefficients (D) for both Li+ and Mg2+ in Ti2S4 were measured using the galvanic intermittent titration technique (GITT) as a function of both ion concentration (x) and temperature. During discharge at 60 degrees C, D-Li descends gradually from 2 X 10(-8) cm(2)/s at x(Li) approximate to 0 to 2 X 10(-9) cm(2)/s at x(Li) approximate to 1.9. In contrast, D-Mg decreases sharply from 2 X 10(-8) to 1 x 10(-12) cm(2)/s by x(Mg) approximate to 0.8. This kinetic factor limits the maximum practical discharge capacity of MgxTi2S4. The difference in behavior vis a vis Li+ implies that either increasing Mg2+ occupation of the tetrahedral site at x(Mg) > 0.6 and/or interactions between diffusing cations play a larger role in mediating the diffusion of divalent compared to monovalent cations. Diffusion activation energies (E-a) extracted from the temperature-dependent data revealed that E-a,E-Mg (540 +/- 80 meV) is about twice that of E-a,E-Li (260 +/- 50 meV), explaining the poorer electrochemical performance of MgxTi2S4 at room temperature.

Bonnick, Patrick J.↗
NASA NTRS2020

Wind Tunnel Testing of High Efficiency Low Power (HELP) Actuation for Active Flow Control

A High Efficiency Low Power (HELP) methodology has been successfully developed and tested on an active flow control (AFC) version of the 10%-scale high-lift Common Research Model (CRM-HL) at the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel (14x22). The AFC variant of the CRM-HL, designated as CRM-SHL-AFC, was integrated with modular HELP actuator cartridges on the shoulder of its highly deflected (≥50°) simplehinged flaps. A reference conventional CRM-HL configuration equipped with Fowler flaps and a nacelle chine was chosen to provide the targeted lift enhancement goals and for comparison with the results of the CRM-SHL-AFC. The current data are presented with the Transonic Wall Interference Correction System (TWICS) method applied. The HELP actuators, which use a combination of unsteady sweeping jets and steady discrete jets in tandem, were designed to overcome strong adverse pressure gradients, while minimizing the pneumatic power usage. The power coefficient (C(sub π)), which takes account of both supply air pressure and mass flow usage for the AFC actuators, is a useful parameter for judging the actuators’ performance efficiency and perhaps for scaling-up of the AFC system. Full HELP actuation coverage with a constant spanwise nozzle pressure ratio (NPR) was able to achieve the lift enhancement goals for the entire lift curve and was found to be the most effective AFC case. For the most effective case, mass flow rates greater than ~0.91 lbm/s and NPR values greater than ~1.8 (corresponding to C(sub π) = 0.18) are needed to achieve the lift enhancement goal at α = 9°, whereas mass flow rates greater than ~1.23 lbm/s and NPR values greater than ~2.3 (corresponding to C(sub π) = 0.3) are needed to achieve the lift enhancement goal at α = 17°. Surface pressure distributions indicate that the HELP actuation on the flap shoulder increased the suction pressures (and flow circulation) globally in both the streamwise and the spanwise directions, and thereby enhanced the lift over the entire high-lift system. The CRM-SHLAFC configuration equipped with HELP actuation was able to match or exceed the lift performance of the reference conventional CRM-HL, thus meeting the objective of the research.

John C Lin↗
DOE OSTI2022

Organic Evaporation and Oxidation Testing in Support of Hanford Sample-and-Send

The Hanford site has approximately 56 million gallons of radioactive mixed waste stored in 177 underground storage tanks. The Hanford Waste Treatment and Immobilization Plant (WTP) is being built to treat and immobilize the tank waste. The baseline method for immobilization of Low Activity Waste (LAW) through the WTP is vitrification, but additional immobilization capacity is needed to supplement the initial LAW melters. An alternative cementitious waste form is being investigated for that future immobilization method to supplement vitrification. However, one impediment to a cementitious waste form is the presence of Land Disposal Restricted (LDR) organic chemicals in tank waste. This work evaluates potential avenues to eliminate that impediment to permit possible use of a cementitious waste form and work towards a decision whether additional LDR organic pretreatment would be required. Savannah River National Laboratory (SRNL) performed testing using simulants to examine evaporation as a method to remove some prevalent organics from LAW. Spiking the caustic LAW simulant with selected regulated organic chemicals found one that clearly decomposes because of caustic instability. Oxidation testing of other organic chemicals found some LDR organics degrade as desired and others are stable in the presence of peroxide and permanganate. In addition to studies with simulants, a literature review was performed to evaluate radiological stability of LDR organics. Descriptions of the experimental details, equipment, and results are included in this report. Evaporation testing consisted of preparing the LAW simulant, spiking that simulant with organic chemicals, and evaporating the mixture via differential distillation. The apparatus was a laboratory-scale vacuum evaporator operated at 60 ±5 torr absolute (vacuum evaporation) and also at atmospheric pressure. The LAW simulant represented the liquid expected to be retrieved from the Hanford tank farms at approximately 4.0 M [Na + ] total sodium ion concentration. The concentration of the organic chemicals added was significantly higher than typically found in the tank waste samples since the higher levels were necessary to assist in analytical measurement and tracking of the spiked species. Organic chemicals were chosen for the work with a consideration of how their volatility compares with that of methanol. This was done by comparing the ratio of the pure water Henry’s law coefficient (K h ) of methanol to that of the compound in question (hereafter termed the K h ratio), where ratios above unity indicated less volatility than methanol. Methanol was chosen because it is a common regulated chemical with relatively low volatility but which has been removed by evaporation in previous laboratory work. While organic separation results depend on evaporator design, laboratory experiments verified that organic partitioning to the overhead condensate stream by evaporation is a practical process. The work reported here found difficulties in quantitative analysis of the organic chemicals in aqueous samples. Most of the time there was insufficient analysis to close a mass balance for evaporator runs, but qualitative evidence of carryover was obtained. The methods were also able to show whether organic chemicals were susceptible or resistant to solution oxidation in permanganate or hydrogen peroxide tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Overview of the Salt at Wall Thermal Exchanges (SWATH) Experiment

Molten Salt Reactors have attracted increased attention in recent years because of the design and safety possibilities offered by the use of a liquid fuel. The European project SAMOFAR1 (2015-2019) is currently coordinating the research efforts on the Molten Fast Salt Reactor (MSFR) concept in Europe. The MSFR is a fast-spectrum breeder reactor with a large negative power coefficient that can be operated in a Thorium fuel cycle. The fast spectrum allows the reduction of the reprocessing requirements and a better reactor breeding ratio. Other advantages of the MSFR design are the possibility for actinide burning and extending fuel resources, on-line fuel loading and reprocessing and the use of novel passive safety systems such as the fuel salt draining system. The reference MSFR concept is a 3000 MW(th) reactor with three different circuits: the fuel circuit, the intermediate circuit and the power conversion system. The main components of the fuel circuit are the fuel salt which serves as fuel and coolant, the core cavity, the inlet and outlet pipes, the gas injection system (not shown in the figure), the salt-bubble separators, the fuel heat exchangers and the pumps. A lithium fluoride salt has been selected as the fuel matrix. The initial composition (non-irradiated) of the MSFR fuel salt is a mixture of a lithium fluoride, thorium fluoride salts and actinides fluoride (LiF-ThF{sub 4}-{sup 233}UF{sub 4} or LiF-ThF{sub 4}- {sup enr}UF{sub 4}-(Pu-MA)F{sub 3}), with the proportion of LiF fixed at about Development of accurate numerical models for the MSFR can only be accomplished by the use of a multiphysics approach which allows taking into account the coupling between relevant neutronics, thermal hydraulics and thermomechanics phenomena. As it is discussed in the next section, such a reactor model requires further progress on the knowledge of molten salt thermal hydraulics, in particular regarding the accuracy of the fuel salt Computational Fluid Dynamics (CFD) models. The Salt at Wall: Thermal Exchanges (SWATH) experiment is being designed to contribute to this particular point. SWATH is one of the experimental activities of the European H2020 SAMOFAR project and includes both experimental and numerical developments. The main objectives of SWATH are to: (i) Improve molten salt numerical models used for design and safety studies and (ii) Demonstrate the working principles of the cold plug device. The cold plug device is a key safety component of the salt draining system of the MSFR. The cold plug is intended to passively melt during accidental conditions thus allowing the flow of the nuclear fuel from the core cavity to the fuel salt draining tanks by gravity forces. The fuel salt draining tanks are designed to ensure that adequate margins are obtained regarding the fuel cooling capabilities and the reactivity control. The next sections provide an overview of the current challenges in molten salt thermal hydraulics modeling and the SWATH experiment design. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources

3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modernizing the active neutron collar poison rod correction using simulations to enhance PWR fresh fuel verification

The declared inventory of 235 U in fresh, low enriched uranium, nuclear fuel assemblies is routinely verified by nuclear safeguards inspectorates using the UNCL (Uranium Neutron Collar – Light Water Reactor Fuel) nondestructive assay instrument. The trend in modern fuels is towards higher initial enrichment, which in turn requires a larger number of burnable poison pins with higher Gadolinia concentration to hold down the initial reactivity. UNCL assay error for modern fuel assemblies typically exceeds 10%. The traditional algorithm used to correct the response for burnable poison content needs revision to achieve accuracy comparable to the $\sigma_\text{R}$ of 4.5% expected for poison-free assemblies. Here we review the problem and use a large set of 287 Monte Carlo simulations based on an experimentally benchmarked model corresponding to the 16×16 PWR array used in Brazilian Angra type II and III fuel. Simulated relative responses are used to evaluate the functional form of the poison correction and determine parameters optimal for 16×16 PWR assemblies. This update is conducted using between 4–24 Gadolinia pins ranging 2–11 wt% Gd 2 O 3 with assemblies having mean enrichments ranging 2.5%–5%. Experimental validation of the updated coefficients using seven measurements shows bias is reduced by an order of magnitude and $\sigma_\text{R}$R reduced from 10% to less than 2%. Two updated sets of coefficients are provided, the most accurate of which is directly useable in the existing analysis code (INCC) in use by inspectorate at fuel fabrication facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Kinetic and temperature sensitivity properties of soil exoenzymes through the soil profile down to one-meter depth at a temperate coniferous forest (Blodgett, CA)

This dataset contains data on kinetic and temperature sensitivity parameters of the exoenzymes β-glucosidase (BG), leucine/leucyl aminopeptidase (LAP) and acid phosphatase (AP) across a soil profile down to 90 cm depth at Blodgett forest, CA, USA. These data were generated to determine if kinetic and thermal properties of microbial exoenzymes involved in organic matter decomposition varied with soil depth, following variation in soil properties and microbial communities, as part of the study: Alves et al. (2021). Kinetic Properties of Microbial Exoenzymes Vary with Soil Depth but Have Similar Temperature Sensitivities Through the Soil Profile. Frontiers in Microbiology 12:3618. https://doi.org/10.3389/fmicb.2021.735282. This research was performed within the framework of the TES Belowground Biogeochemistry SFA project, in particular association with a long-term field warming experiment of the whole soil profile at Blodgett forest. Samples for this work were collected from locations representative of the field experimental plots. Potential enzyme activity rates were measured using laboratory fluorometric assays with soils collected at 0-10, 10-20, 30-40, 50-60, 60-70 and 80-90 cm deep in biological triplicates (i.e., three soil cores collected at different representative locations). Assays with each soil were conducted over a gradient of eight substrate concentrations per enzyme, and incubated at 4, 10, 16, 25, 35 or 50°C. Enzyme Michaelis-Menten kinetics were modeled over the eight substrate concentrations at each temperature, and the temperature sensitivity of the kinetic parameters was modeled over the six temperatures using linear Arrhenius/Q10 and non-linear Macromolecular Rate Theory (MMRT) models. The dataset includes the fully processed enzyme activity rate data used to model Michaelis-Menten kinetics, calculated kinetic and temperature sensitivity parameters, and basic soil and microbial biomass chemistry for each sample. All data is provided for each individual biological replicate, and kinetic and temperature sensitivity parameters are also provided as means of the biological replicates. The dataset also includes all raw measurement data and code used to parse, combine and perform the analyses described by Alves et al. (2021). For file descriptions, see the file-level metadata files: “enzymes_dataProcessed_flmd.xlsx” (processed data); “enzymes_dataRaw_flmd.xlsx” (compressed raw data and metadata); and “enzymes_code_flmd.xlsx” (compressed code). The experimental design, list of parameters measured, soil and microbial biomass chemistry data, and means of biological replicates for kinetic and temperature sensitivity parameters are also provided as human-readable tables in file “enzymeTraits_design_results_tables.xlsx”.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

Test of a model for trace element partition during closed-system solidification of a silicate liquid

Trace-element concentrations in separated minerals and rock fragments from a coarsely crushed and sieved sample of medium-grained lunar basalt are determined by high-precision instrumental neutron activation analysis for three main purposes. These are: (1) to test a previously proposed model for describing trace-element behavior during solidification of a silicate liquid under conditions of a closed system; (2) to compare trace-element concentrations among different size fractions of a comminuted basalt; and (3) to investigate small-scale heterogeneity within a single medium-grained basalt. It is found that an excellent mass balance for the whole rock is attained in terms of the trace-element concentrations of the component minerals and mesostasis, that the mixing model describing this mass balance yields a superior modal analysis for the whole rock, and that different size fractions of the coarsely crushed basalt vary in mineral composition. The closed-system model is shown to account properly for the average behavior of the trace elements during solidification of the basalt, and values of distribution coefficients are obtained for incompatible elements in that rock.

Haskin, L. A.↗