SEARCH · Engineering Papers
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.
Focal Area(s): (2) Predictive modeling through the use of AI techniques and AI-derived model components: This white paper focuses on using artificial intelligence (AI) to discover connections across processes and scales and to represent those processes in predictive models of anthropogenic impacts on the water cycle.
The anomalous magnetic moment of muon has been used as an indication of physics beyond the Standard Model. Fermilab Muon g-2 experiment has published the most precise measurement of the anomalous magnetic moment of muon with an uncertainty of 460 part-per-billion which showed 4.2 sigma discrepancy between the experiment average and the 2020 g-2 Theory Initiative Standard Model prediction. This talk will show the recent results from the Fermilab Muon g-2 experiment while also briefly cover the status of the Standard Model prediction on calculating the anomalous magnetic moment of muon.
The deformation behavior of commercially pure rolled titanium subjected to strain path changes is studied using experiments and a crystal plasticity model. Four different loading combinations are performed at room temperature to study the activation of slip, twinning, de-twinning and double-twinning in hexagonal closed packed titanium. The strain paths considered are: rolling direction compression (RDC) followed by normal direction compression (NDC), RDC followed by transverse direction compression (TDC), NDC followed by RDC, and NDC followed by TDC. An EBSD-based analysis of the distribution of boundary misorientation angles before and after reload was developed to analyze the evolution of {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twins. This analysis supports the model results concerning the treatment of twin reorientation. A de-twinning and double-twinning model accounting for back stress effects, an important feature of strain path changes, is implemented within the framework of the visco-plastic self-consistent (VPSC) model along with a dislocation density (DD) based hardening scheme. In the model, plasticity is accommodated by prismatic $\langle$a$\rangle$, basal $\langle$a$\rangle$ and pyramidal $\langle$c+a$\rangle$ slip modes, and {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twinning modes. The VPSC model predicts the evolution of twinning, de-twinning and double-twinning processes for both tensile and compression twinning modes under strain path change. The model predicts macroscopic stress-strain response, texture evolution, and twin volume fraction that are in agreement with experimental observations. The evolution of texture is investigated in detail by separately analyzing the twinned domains, rather than the evolution of the global texture.
The Advance Sensors and Instrumentation fiscal year report on the direct application of the HTIR-TC drift model on other commercial thermocouples. The Type N and K thermocouples are useful in irradiation environments up to 1100°C (or 1290°C for short bursts and/or single use), but over temperature and high fluence can cause the TCs to drift. The HTIR-TC drift model predicts the end use environment of temperature and irradiation and the total time scale to show the amount of drift a thermocouple will undergo. The model predicts that HTIR-TCs at full temperature range and fluence of 10^21 N/cm^2 will exhibit a drift of under -1%. The type N and K TCs will see a drift of -0.5% at the same fluence, but temperatures below 1100°C. Above this temperature the Type N and K TCs quickly de-calibrate and will give erroneous temperatures measurements.
Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.
Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. In this paper, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model’s accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal–metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.
ABSTRACT Using the GAEA semi-analytic model, we analyse the connection between Damped Ly α systems (DLAs) and H i in galaxies. Our state-of-the-art semi-analytic model is tuned to reproduce the local galaxy H i mass function, and that also reproduces other important galaxy properties, including the galaxy mass–gas metallicity relation. To produce catalogues of simulated DLAs we throw 105 random lines of sight in a composite simulated volume: dark matter haloes with log$(\frac{M_{200}}{ {\rm M}_{\odot }}) \ge 11.5$ are extracted from the Millennium Simulation, while for $9.2 \le \log (\frac{M_{200}}{ \mathrm{M}_{\odot }})\lt 11.5$ we use the Millennium II, and for $8 \le \log (\frac{M_{200}}{\mathrm{M}_{\odot }}) \lt 9.2$ a halo occupation distribution model. At 2 < z < 3, where observational data are more accurate, our fiducial model predicts the correct shape of the column density distribution function, but its normalization falls short of the observations, with the discrepancy increasing at higher redshift. The agreement with observations is significantly improved increasing both the H i masses and the disc radii of model galaxies by a factor of 2, as implemented ‘a posteriori’ in our 2M−2R model. In the redshift range of interest, haloes with $M_{200} \ge {10}^{11} \, \mathrm{M}_{\odot }$ give the major contribution to ΩDLA, and the typical DLA host halo mass is $\sim \!{10}^{11} \, \mathrm{M}_{\odot }$. The simulated DLA metallicity distribution is in relatively good agreement with observations, but our model predicts an excess of DLAs at low metallicities. Our results suggest possible improvements for the adopted modelling of the filtering mass and metal ejection in low-mass haloes.
Recently, heat transport was investigated using a directly driven beryllium sphere [Farmer et al., Phys. Plasmas 27, 082701 (2020)]. Models that overly restrict heat transport were rejected. This paper extends work to directly driven gold spheres where radiation loss is more important. Here, gold coated spheres are directly driven at the OMEGA laser facility at intensities of 5×10 14 W/cm 2 . Plasma conditions, laser coupling, and x-ray flux are all measured. Additionally, comparisons to 2D radiation-hydrodynamic simulations are performed. Simulations use three common heat transport models: local transport with flux limiters of f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. It is shown that both the SNB model and f = 0.15 match the measured plasma conditions with the SNB model better capturing the temporal evolution of electron temperature. The f = 0.03 model predicts too low of an electron density and too hot of a temperature. The measured scattered light is roughly 6% of the incident energy, the f = 0.15 and SNB models predict 0.5% uncoupled light, and f = 0.03, 38% uncoupled light. The x-ray fluxes in the f = 0.15 and SNB simulations rise too quickly and are just outside the measurement's error, while the x-ray flux in the f = 0.03 simulation is low by a factor of two-three. For these reasons, the f = 0.03 model is rejected.
The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.
The rapid changes occurring in the polar regions require an improved understanding of the processes that are driving these changes. At the same time, increased human activities such as marine navigation, resource exploitation, aviation, commercial fishing, and tourism require reliable and relevant weather information. One of the primary goals of the World Meteorological Organization's Year of Polar Prediction (YOPP) project is to improve the accuracy of numerical weather prediction (NWP) at high latitudes. During YOPP, two Canadian “supersites” were commissioned and equipped with new ground-based instruments for enhanced meteorological and system process observations. Additional pre-existing supersites in Canada, the United States, Norway, Finland, and Russia also provided data from ongoing long-term observing programs. These supersites collected a wealth of observations that are well suited to address YOPP objectives. In order to increase data useability and station interoperability, novel Merged Observatory Data Files (MODFs) were created for the seven supersites over two Special Observing Periods (February to March 2018 and July to September 2018). All observations collected at the supersites were compiled into this standardized NetCDF MODF format, simplifying the process of conducting pan-Arctic NWP verification and process evaluation studies. This paper describes the seven Arctic YOPP supersites, their instrumentation, data collection and processing methods, the novel MODF format, and examples of the observations contained therein. MODFs comprise the observational contribution to the model intercomparison effort, termed YOPP site Model Intercomparison Project (YOPPsiteMIP). All YOPPsiteMIP MODFs are publicly accessible via the YOPP Data Portal (Whitehorse: https://doi.org/10.21343/a33e-j150, Huang et al., 2023a; Iqaluit: https://doi.org/10.21343/yrnf-ck57, Huang et al., 2023b; Sodankylä: https://doi.org/10.21343/m16p-pq17, O'Connor, 2023; Utqiagvik: https://doi.org/10.21343/a2dx-nq55, Akish and Morris, 2023c; Tiksi: https://doi.org/10.21343/5bwn-w881, Akish and Morris, 2023b; Ny-Ålesund: https://doi.org/10.21343/y89m-6393, Holt, 2023; and Eureka: https://doi.org/10.21343/r85j-tc61, Akish and Morris, 2023a), which is hosted by MET Norway, with corresponding output from NWP models.
In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.
A credible simulation of disposal room porosity at the Waste Isolation Pilot Plant (WIPP) requires a tenable compaction model for the 55-gallon waste containers within the room. A review of the legacy waste material model, however, revealed several out-of-date and untested assumptions that could affect the model’s compaction behavior. For example, the legacy model predicted non-physical tensile out-of-plane stresses under plane strain compression. (Plane strain compression is similar to waste compaction in the middle of a long drift.) Consequently, a suite of new compaction experiments were performed on containers filled with surrogate, non-degraded, waste. The new experiments involved uniaxial, triaxial, and hydrostatic compaction tests on quarter-scale and full-scale containers. Special effort was made to measure the volume strain during uniaxial and triaxial tests, so that the lateral strain could be inferred from the axial and volume strain. These experimental measurements were then used to calibrate a pressure dependent, viscoplastic, constitutive model for the homogenized compaction behavior of the waste containers. This new waste material model’s predictions agreed far better with the experimental measurements than the legacy model’s predictions, especially under triaxial and hydrostatic conditions. Under plane strain compression, the new model predicted reasonable compressive out-of-plane stresses, instead of tensile stresses. Moreover, the new model’s plane strain behavior was substantially weaker for the same strain, yet substantially stronger for the same porosity, than the legacy model’s behavior. Although room for improvement exists, the new model appears ready for prudent engineering use.
Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.
Electromagnetic Ion Cyclotron (EMIC) wave scattering has been proved to be responsible for the fast loss of both radiation belt (RB) electrons and ring current (RC) protons. However, its role in the concurrent dropout of these two co–located populations remains to be quantified. In this work, we study the effect of EMIC wave scattering on both populations during the 27 February 2014 storm by employing the global physics–based RAM–SCB model. Throughout this storm event, MeV RB electrons and 100s keV RC protons experienced simultaneous dropout following the occurrence of intense EMIC waves. By implementing data–driven initial and boundary conditions, we perform simulations for both populations through the interplay with EMIC waves and compare them against Van Allen Probes observations. Notably, the results indicate that by including EMIC wave scattering loss, especially by the He–band EMIC waves, the model aligns closely with data for both populations. Additionally, we investigate the simulated pitch angle distributions (PADs) for both populations. Including EMIC wave scattering in our model predicts a 90° peaked PAD for electrons with stronger losses at lower pitch angles, while protons exhibit an isotropic PAD with enhanced losses at pitch angles above 40°. Furthermore, our model predicts considerable precipitation of both particle populations, predominantly confined to the afternoon to midnight sector (12 hr < MLT < 24 hr) during the storm's main phase, corresponding closely with the presence of EMIC waves.
Applied potentials have been demonstrated as a powerful tool to promote heterogeneous Brønsted acid catalysis by orders of magnitude, leveraging interfacial electric fields to stabilize protonated intermediates. However, the use of flat two-dimensional electrodes with inherently low active site densities limits the application of conventional thermochemical characterization techniques that can probe the nature of catalytic active sites. Here, we use kinetic analyses with an electrostatics-based model to elucidate the intricacies of potential-induced rate promotion, employing liquid-phase dehydration of 1-methylcyclopentanol catalyzed by carboxylic acid groups on carbon nanotubes as a probe system. By using a basket electrode to directly polarize catalyst powder, we demonstrate that thermocatalytic reaction rates can be promoted by 100,000-fold, exhibiting a log–linear dependence on applied potential with rate-potential scalings as high as 125 ± 4 mV per 10-fold rate increase. In agreement with model predictions, we show that lower ionic strengths attenuate potential sensitivity, resulting from a weakening of the interfacial electric field that interacts with the acidic proton. Furthermore, we experimentally confirm the model-predicted “isokinetic potential” (at ∼0.6 V vs Ag/AgCl)─the potential at which all rate scaling lines at various ionic strengths intersect, making the rate independent of ionic strength. Base titrations reveal that only ∼8% of the carboxylic acid sites are catalytically active, yet these same active sites are operational at the highest and lowest potentials. Collectively, our results provide a key methodology for modeling catalytic effects of electric fields, quantifying active sites under applied potential, and demonstrating fundamental principles of electric field-induced rate promotion.
Predictive maintenance (PdM) has shown great potential for achieving substantial cost savings and enhancing the economic competitiveness of nuclear power plants (NPPs) in today's energy market. Among the different modeling approaches that exist, machine learning (ML) tools in particular have a demonstrated ability to handle high dimensional and multivariate data and to extract hidden relationships within data in industrial environments. While ML methods show great potential, their lack of explainability---especially for black-box models---is a major hurdle to their adoption. Moreover, considering the supposed trade-off between explainability and performance challenges, careful consideration must be made as to which of these quality aspects takes precedence in light of multiple modeling options, resource availability, and domain characteristics. The present work evaluates the performance of six ML models, each with a different degree of explainability, in classifying the conditions of circulating water pumps (CWPs) by utilizing sensor data from nuclear power plants. To determine the drivers behind the trade-offs presented by this array of models, this work also tests different combinations of CWP units as the training and testing data, degrees of data imbalance, and objective functions for hyperparameter tuning. It was found that black-box models tend to afford superior performance in cases where there are far more instances of one type of labeled data than of any other type. It is recommended that a guided procedure be followed for designing and delivering an ML system that is sufficiently explainable to all involved stakeholders.
This paper presents results from system -level modeling of a water -based reactor cavity cooling system using RELAP5-3D. The computational model is benchmarked with experimental data from a half -scale RCCS test facility at Argonne National Laboratory. The model prediction is first compared with a two-phase oscillatory baseline experimental case where mixed accuracy is obtained. The model shows reasonable prediction of mass flow rate, pressure, and temperature but significant overprediction of void fraction. The model prediction is then compared with a fault case where the inlet of the risers is gradually reduced using a throttling valve. As the valve is closed, the model is able to predict some major flow phenomena observed in the experiment such as the dampening of oscillations, the reintroduction of oscillations, as well as boiling, flashing, and geysering in the risers. However, the timeline of these events are not well captured by the model. The model is also used to investigate the evolution of flow regime in the chimney. This work highlights that the semiempirical constitutive relations used in RELAP-3D could have a strong influence on the accuracy of the model in two-phase oscillatory flows.
Abstract Developing prediction models for emerging infectious diseases from relatively small numbers of cases is a critical need for improving pandemic preparedness. Using COVID-19 as an exemplar, we propose a transfer learning methodology for developing predictive models from multi-modal electronic healthcare records by leveraging information from more prevalent diseases with shared clinical characteristics. Our novel hierarchical, multi-modal model ( $${\textsc {TransMED}}$$ T R A N S MED ) integrates baseline risk factors from the natural language processing of clinical notes at admission, time-series measurements of biomarkers obtained from laboratory tests, and discrete diagnostic, procedure and drug codes. We demonstrate the alignment of $${\textsc {TransMED}}$$ T R A N S MED ’s predictions with well-established clinical knowledge about COVID-19 through univariate and multivariate risk factor driven sub-cohort analysis. $${\textsc {TransMED}}$$ T R A N S MED ’s superior performance over state-of-the-art methods shows that leveraging patient data across modalities and transferring prior knowledge from similar disorders is critical for accurate prediction of patient outcomes, and this approach may serve as an important tool in the early response to future pandemics.