Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive simulation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The predicted shapes of voids and Xe bubbles in UO 2

Morphology is a fundamental attribute when investigating voids and bubbles in UO 2 . Here we use molecular dynamics and Monte Carlo simulations to predict the lowest energy shapes for voids and bubbles in UO 2 . The energies of the {100}, {110} and {111} surfaces have been calculated and used to predict the equilibrium void shape from Wulff construction. This equilibrium shape is compared to low energy faceted voids exhibiting different relative proportions of each family of terminating surfaces. It is found that the equilibrium Wulff shape does not represent the lowest energy morphology for nm void sizes at temperatures between 300 K and 1200 K. Furthermore, the lowest energy faceted voids are slightly more energetically favourable than spherical voids, and as Xe is added, and bubble pressure increases, the faceted morphology becomes even more favourable than the spherical shape.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Do chaotic field lines cause fast reconnection in coronal loops?

Over the past decade, Boozer has argued that three-dimensional (3D) magnetic reconnection fundamentally differs from two-dimensional reconnection due to the fact that the separation between any pair of neighboring field lines almost always increases exponentially over distance in a 3D magnetic field. According to Boozer, this feature makes 3D field-line mapping chaotic and exponentially sensitive to small non-ideal effects; consequently, 3D reconnection can occur without intense current sheets. We test Boozer's theory via ideal and resistive reduced magnetohydrodynamic simulations of the Boozer–Elder coronal loop model driven by sub-Alfvénic footpoint motions [A. H. Boozer and T. Elder, Phys. Plasmas 28, 062303 (2021)]. In this work, our simulation results significantly differ from their predictions. The ideal simulation shows that Boozer and Elder under-predict the intensity of current density due to missing terms in their reduced model equations. Furthermore, resistive simulations of varying Lundquist numbers show that the maximal current density scales linearly rather than logarithmically with the Lundquist number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the Morphodynamics of a Wide Class of Large‐Scale Meandering Rivers: Insights Gained by Coupling LES With Sediment‐Dynamics

Abstract In meandering rivers, interactions between flow, sediment transport, and bed topography affect diverse processes, including bedform development and channel migration. Predicting how these interactions affect the spatial patterns and magnitudes of bed deformation in meandering rivers is essential for various river engineering and geoscience problems. Computational fluid dynamics simulations can predict river morphodynamics at fine temporal and spatial scales but have traditionally been challenged by the large scale of natural rivers. We conducted coupled large‐eddy simulation and bed morphodynamics simulations to create a unique database of hydro‐morphodynamic data sets for 42 meandering rivers with a variety of planform shapes and large‐scale geometrical features that mimic natural meanders. For each simulated river, the database includes (a) bed morphology, (b) three‐dimensional mean velocity field, and (c) bed shear stress distribution under bankfull flow conditions. The calculated morphodynamics results at dynamic equilibrium revealed the formation of scour and deposition patterns near the outer and inner banks, respectively, while the location of point bars and scour regions around the apexes of the meander bends is found to vary as a function of the radius of curvature of the bends to the width ratio. A new mechanism is proposed that explains this seemingly paradoxical finding. The high‐fidelity simulation results generated in this work provide researchers and scientists with a rich numerical database for morphodynamics and bed shear stress distributions in large‐scale meandering rivers to enable systematic investigation of the underlying phenomena and support a range of river engineering applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advancing AI-Driven Analysis in X-ray Absorption Spectroscopy: Spectral Domain Mapping and Universal Models

In recent years, rapid progress has been made in developing artificial intelligence (AI) and machine learning (ML) methods for X-ray absorption spectroscopy (XAS) analysis. Compared to traditional XAS analysis methods, AI/ML approaches offer dramatic improvements in efficiency and help eliminate human bias. To advance this field, we advocate an AI-driven XAS analysis pipeline that features several interconnected key building blocks: benchmarks, workflows, databases, and AI/ML models. Specifically, we present two case studies for XAS ML. In the first study, we demonstrate the importance of reconciling the discrepancies between simulation and experiment using spectral domain mapping (SDM). Our ML model, which is trained solely on simulated spectra, predicts an incorrect oxidation state trend for Ti atoms in a combinatorial zinc titanate film. After transforming the experimental spectra into a simulation-like representation using SDM, the same model successfully recovers the correct oxidation state trend. In the second study, we explore the development of universal XAS ML models that are trained on the entire periodic table, which enables them to leverage common trends across elements. Looking ahead, we envision that an AI-driven pipeline can unlock the potential of real-time XAS analysis to accelerate scientific discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2021

WIRE: Resource-efficient Scaling with Online Prediction for DAG-based Workflows

This paper introduces WIRE that manages resources for the DAG-based workflows on IaaS clouds. WIRE predicts and plans resources over the MAPE (Monitor-Analyze-Plan-Execute) loops to: 1) Estimate task performance with online data, 2) Conduct simulations to predict the upcoming loads based on online estimates and workflow DAGs, 3) Apply a resource-steering policy to size cloud instance pools for the maximal parallelism that is consistent with low cost. We implement WIRE on Pegasus WMS/HTCondor and evaluate its performance on the ExoGENI network cloud. The results show that WIRE attains low resource cost with the performance that is typically within a factor of two of optimal.

Xie, Bing↗
DOE OSTI2021

End-to-End Simulations of a 3.4-Meter Detector Wall for Neutron-Diagnosed Subcritical Experiments

The Nevada National Security Site, together with Los Alamos National Laboratory and Lawrence Livermore National Laboratory, is developing a novel diagnostic to measure the reactivity of subcritical experiments. This capability is known as neutron-diagnosed subcritical experiments. The decay of the fission gamma rays from the neutron-interrogated subcritical experiment is measured as a function of time with a large (~3-meter diameter) detector wall consisting of 151 individual detector pixels. The data from this current mode measurement inform the neutron multiplication factor, keff, and thus the relative reactivity of the subcritical experiment. The Nevada National Security Site developed the Gamma Array Simulation Toolkit initially to help inform the design of the individual detector pixels and the 3.4-meter diameter detector wall to be fielded as part of neutron-diagnosed subcritical experiments. This toolkit is now being used to simulate and predict the performance of the final design of the individual detector pixels and the aggregate detector wall. Key detector characteristics evaluated from these simulations include impulse response, pulse height spectrum, number of photoelectrons per MeV, detector efficiency, and cross talk between detector pixels. Collectively, the results of these simulations inform how well the fission gamma ray die-off distribution from a neutron-diagnosed subcritical experiment measurement can be resolved. This is critical to determining the relative reactivity of the experiment. Furthermore, the Gamma Array Simulation Toolkit can be used to aid in the analysis of the experimental data from a neutron-diagnosed subcritical experiment measurement once the simulation has been benchmarked.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

End-to-End Simulations of a 3.4-Meter Detector Wall for Neutron-Diagnosed Subcritical Experiments

The Nevada National Security Site, together with Los Alamos National Laboratory and Lawrence Livermore National Laboratory, is developing a novel diagnostic to measure the reactivity of subcritical experiments. This capability is known as neutron-diagnosed subcritical experiments. The decay of the fission gamma rays from the neutron-interrogated subcritical experiment is measured as a function of time with a large (~3-meter diameter) detector wall consisting of 151 individual detector pixels. The data from this current mode measurement inform the neutron multiplication factor, keff, and thus the relative reactivity of the subcritical experiment. The Nevada National Security Site developed the Gamma Array Simulation Toolkit initially to help inform the design of the individual detector pixels and the 3.4-meter diameter detector wall to be fielded as part of neutron-diagnosed subcritical experiments. This toolkit is now being used to simulate and predict the performance of the final design of the individual detector pixels and the aggregate detector wall. Key detector characteristics evaluated from these simulations include impulse response, pulse height spectrum, number of photoelectrons per MeV, detector efficiency, and cross talk between detector pixels. Collectively, the results of these simulations inform how well the fission gamma ray die-off distribution from a neutron-diagnosed subcritical experiment measurement can be resolved. This is critical to determining the relative reactivity of the experiment. Furthermore, the Gamma Array Simulation Toolkit can be used to aid in the analysis of the experimental data from a neutrondiagnosed subcritical experiment measurement once the simulation has been benchmarked.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

End-to-End Simulations of a 3.4-Meter Detector Wall for Neutron-Diagnosed Subcritical Experiments

The Nevada National Security Site, together with Los Alamos National Laboratory and Lawrence Livermore National Laboratory, is developing a novel diagnostic to measure the reactivity of subcritical experiments. This capability is known as neutron-diagnosed subcritical experiments. The decay of the fission gamma rays from the neutron-interrogated subcritical experiment is measured as a function of time with a large (~3-meter diameter) detector wall consisting of 151 individual detector pixels. The data from this current mode measurement inform the neutron multiplication factor, keff, and thus the relative reactivity of the subcritical experiment. The Nevada National Security Site developed the Gamma Array Simulation Toolkit initially to help inform the design of the individual detector pixels and the 3.4-meter diameter detector wall to be fielded as part of neutron-diagnosed subcritical experiments. This toolkit is now being used to simulate and predict the performance of the final design of the individual detector pixels and the aggregate detector wall. Additionally, key detector characteristics evaluated from these simulations include impulse response, pulse height spectrum, number of photoelectrons per MeV, detector efficiency, and cross talk between detector pixels. Collectively, the results of these simulations inform how well the fission gamma ray die-off distribution from a neutron-diagnosed subcritical experiment measurement can be resolved. This is critical to determining the relative reactivity of the experiment. Furthermore, the Gamma Array Simulation Toolkit can be used to aid in the analysis of the experimental data from a neutron-diagnosed subcritical experiment measurement once the simulation has been benchmarked.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks

Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.

Sun, Yixuan (ORCID:0000000311093380)↗
DOE OSTI2021

Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules

Advances in machine learning provide the ability to leverage data from expensive simulations of high-energy-density experiments to significantly cut down on computational time and costs associated with the search for optimal target designs. This study presents an application of cutting-edge Bayesian optimization methods to the one-dimensional (1D) design optimization of double shell graded layer targets for inertial confinement fusion experiments. This investigation attempts to reduce hydrodynamic instabilities while retaining high yields for future NIF experiments. Machine learning methods can use predictive physics simulations to identify graded layer designs from within the vast design space that demonstrate high predicted performance, including novel designs with high uncertainty in performance that may hold unexpected promise. By applying machine learning tools to the simulation design, we map the trade-off between 1D yield and instability, specifically isolating parameter ranges, which maintain high performance while showing significantly improved Rayleigh–Taylor stability over the point design. Furthermore, the groundwork laid in this study will be a useful design tool for future NIF experiments with graded layer targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Early Solar Wind and Dynamo Magnetic Field Topology Predictions for (16) Psyche and Other Asteroids

Abstract Asteroid (16) Psyche is a metal‐rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field topology for several distinct pathways: (i) an early solar wind‐induced magnetization imparted after a larger body was impacted, forming the present‐day asteroid, (ii) a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile, and (iii) magnetization in the turbulent solar nebula disk. For pathway (i) we find the field to be predominantly dipolar and spin axis‐aligned. For pathway (ii) we find the field to be either dipolar and spin axis‐misaligned, or highly multipolar. We also find that (iii) a field produced earlier before the solar nebula cleared, would be highly multipolar. In cases (i) and (ii) we also place constraints on the field strength. Simple detection of a magnetic field without constraining its topology and temporal variability would be insufficient to confirm a remanent source, due to the influence of the present‐day solar wind, electromagnetic induction, and (16) Psyche's high obliquity. For sufficiently strong fields however, the field topology and orientation may reveal key observable consequences of the nature and history of (16) Psyche. Our framework is also broadly applicable to the study of magnetic fields from other asteroids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

We report Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

13 HYDRO ENERGY↗
DOE OSTI2020

Recreating Giants Impacts in the Laboratory: Shock Compression of MgSiO 3 Bridgmanite to 14 Mbar

Understanding giant impacts requires accurate description of how extreme pressures and temperatures affect the physical properties of the constituent materials. Here, we report shock experiments on two polymorphs of MgSiO 3 : enstatite and bridgmanite (perovskite) crystals.We obtain pressure-density shock equation of state to 14 Mbar and more than 9 g/cm 3 , a 40% increase in density from previous data on MgSiO 3 . Density-functional-theory molecular dynamics (DFT-MD) simulations provide predictions for the shock Hugoniot curves for bridgmanite and enstatite and suggest that the Grüneisen parameter decreases with increasing density. The good agreement between the simulations and the experimental data, including for the shock temperature along the enstatite Hugoniot reveals that DFT-MD simulations reproduce well the behavior of dense fluid MgSiO 3 .We also reveal a high optical reflectance indicative of a metal-like electrical conductivity which supports the hypothesis that magma oceans may contribute to planetary magnetic field generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2021

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗