Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗
DOE OSTI2022

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Threshold Helium Leakage Detection Switch with Ultra Low Power Operation

Detecting helium leakage is important in many applications, such as in dry cask nuclear waste storage systems. This work develops a helium detection system based on the relative permittivity (dielectric constant) difference between air and helium. This difference changes the status of an electrostatic microelectromechanical system (MEMS) switch. The switch is a capacitive-based device and requires a very negligible amount of power. Exciting the switch’s electrical resonance enhances the MEMS switch sensitivity to detect low helium concentration. This work simulates two different MEMS switch configurations: a cantilever-based MEMS modeled as a single-degree-freedom model and a clamped-clamped beam MEMS molded using the COMSOL Multiphysics finite-element software. While both configurations demonstrate the switch’s simple operation concept, the clamped-clamped beam was selected for detailed parametric characterization due to its comprehensive modeling approach. The beam detects at least 5% helium concentration levels when excited at 3.8 MHz, near electrical resonance. The switch performance decreases at lower excitation frequencies or increases the circuit resistance. The MEMS sensor detection level was relatively immune to beam thickness and parasitic capacitance changes. However, higher parasitic capacitance increases the switch’s susceptibility to errors, fluctuations, and uncertainties.

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

Quantum speed limit for the out-of-time-ordered correlator from an open-system perspective

Scrambling, the delocalization of initially localized quantum information, is commonly characterized by the out-of-time-ordered correlator (OTOC). Employing the OTOC–Renyi-2 entropy theorem, we derive a quantum speed limit for the OTOC, which sets a lower bound for the rate with which information can be scrambled. This bound becomes particularly tractable by describing the scrambling of information in a closed quantum system as an effective decoherence process of an open system interacting with an environment. We prove that decay of the OTOC can be bounded by the strength of the system-environment coupling and two-point environmental correlation functions. We validate our analytic bound numerically using the nonintegrable transverse field Ising model. Furthermore, our results provide a universal and model-agnostic quantitative framework for understanding the dynamical limits of information spreading across quantum many-body physics, condensed matter systems, and engineered quantum platforms.

Fermions↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

DiscoFluxM version 1.x

DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.

Luscher, Darby↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

Here, in this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling strategy for progressive failure prediction in lithium-ion batteries under mechanical abuse

The prediction for the internal failure of lithium-ion batteries (LIBs) under external mechanical abuse loading remains a challenge for safe design. This paper systematically studies the modeling approach for progressive failure simulation and short-circuit prediction. Hemispherical indentation experiments are conducted on a battery specimen containing 5 representative sandwich (RS) layers to examine the deformation and progressive failure behavior of battery components. Three different models, viz, a high-fidelity detailed model, an intermediate homogenized model and a fully homogenized model, are developed and employed to simulate the progressive failure in a multi-layer battery specimen. The simulation results obtained using the three models all correlate well with experimental phenomena, capturing the overall stress-strain response and ultimate failure of the battery structure. The numerical results can help us understand the failure behavior and the mechanism for mechanical abuse-induced short circuit in LIBs. The feasibility of the three models for use in progressive failure prediction is compared by analyzing the capability of the models in predicting the critical failure events, the computational efficiency of the models, and the feasibility to implement coupled mechanical-electrochemical-thermal simulations. The results of this work provide useful insights on practical choices for the modeling strategy and safety design of LIBs under mechanical abuse conditions.

25 ENERGY STORAGE↗
DOE OSTI2023

Thermal-Hydraulic Analysis of the AFIP-7 Irradiation Test in the Advanced Test Reactor - Oxide Growth Prediction and Correlation

Knowing the thickness of the oxide layer on the surface of aluminum fuel cladding is vitally important for predicting fuel temperature due to the low thermal conductivity of the oxide layer. Several correlation models for predicting oxide growth can be found in the literature. In previous research, the correlations were combined with heat transfer simulations in Abaqus, a finite element analysis code, to forecast the oxide growth. However, this approach requires heat transfer coefficients for modeling heat exchanges with the external flow field, and such coefficients were obtained through empirical equations. Since different empirical equations yield varying heat transfer coefficients, the cladding temperature and predicted oxide thickness both carry a high degree of uncertainty. Here, this research develops a new approach that integrates the fluid flow, fluid and solid heat transfer, and oxide growth correlation(s) into a single computational fluid dynamics model. We demonstrate this approach’s ability to predict oxide development on the AFIP-7 plates during two Advanced Test Reactor (ATR) irradiation cycles. The projected oxide thickness falls within the experimental measurements taken during post irradiation examination.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

An Adaptive Geometry-Free Thermo-Mechanical Model for Directed Energy Deposition Process Modeling

This presentation describes a novel, geometry-free thermo-mechanical model with adaptive subdomain con- struction to accurately predict the thermal conditions, distortions, and residual stresses throughout the directed energy deposition (DED) process. A novel finite element workflow is designed to con- duct the numerical analysis, based on the multi-app and data transfer capabilities in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE). Unlike with traditional methods, the part geometry in this model is not predefined. Instead, it is a combined effect of the processing parameters and material properties. At each time step, the model utilizes a subdomain construction paradigm to model the material deposition. A specialized mesh adaptivity scheme is incorporated to provide an accurate prediction while reducing the overall computational cost. The results generated by the proposed model show general agreement with the experimental measurements for the single track scan with varying processing parameters and demonstrate reasonable predictions for higher material buildups.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A framework for data-driven solution and parameter estimation of PDEs using conditional generative adversarial networks

We employ and adapt the image-to-image translation concept based on conditional generative adversarial networks (cGAN) for learning a forward and an inverse solution operator of partial differential equations (PDEs). We focus on steady-state solutions of coupled hydromechanical processes in heterogeneous porous media and present the parameterization of the spatially heterogeneous coefficients, which is exceedingly difficult using standard reduced-order modeling techniques. We show that our framework provides a speed-up of at least 2,000 times compared to a finite-element solver and achieves a relative root-mean-square error (r.m.s.e.) of less than 2% for forward modeling. For inverse modeling, the framework estimates the heterogeneous coefficients, given an input of pressure and/or displacement fields, with a relative r.m.s.e. of less than 7%, even for cases where the input data are incomplete and contaminated by noise. The framework also provides a speed-up of 120,000 times compared to a Gaussian prior-based inverse modeling approach while also delivering more accurate results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anomaly Detection in Materials Digital Twins with Multiscale ICME for Additive Manufacturing

Detecting anomaly in fatigue and fracture experimental materials science is an interesting yet challenging topic. The reasons are threefold. First, the anomalous microstructure feature that gives rise to structural failure is small, sometimes in the order of 10 -7 of the interrogated volume. This, in turn, results in a highly imbalanced classification problem in machine learning (ML). Second, the consequence is high, in the sense that the test specimen is destructed in such case. Third, the convolution between microstructure stochasticity and the small probability of void nucleation, growth, and coalescence makes failure and fracture a hard-to-predict and challenging problem in materials science due to its irreproducibility, even experimentally. In this paper, we developed a materials digital twin and applied anomaly detection methods to detect voids and anomaly in additive manufacturing (AM). The materials digital twin is driven by two integrated computational materials engineering (ICME) models, which are kinetic Monte Carlo (kMC) and crystal plasticity finite element method (CPFEM). In conclusion, we demonstrated that by using anomaly detection, it is possible to detect voids and other defects in materials digital twin, which paves way for future research in integrating materials digital twin with its physical counterpart.

ICME↗
DOE OSTI2020

Transient magnetic field diffusion considerations relevant to magnetically assisted indirect drive inertial confinement fusion

We report application of a magnetic field to an indirect drive inertial confinement fusion target requires diffusion of the field through the high-Z and electrically conducting Hohlraum. The onset of the external field generates eddy currents in the Hohlraum wall that result in (1) a reduction of the peak field at the capsule, (2) heating of the Hohlraum wall through Ohmic dissipation, and (3) wall movement due to the inward force from the eddy current interacting with the field. Heating of the wall causes an increase in blackbody radiation which can preheat the capsule and frozen deuterium–tritium fuel, while wall motion leads to potential misalignment of the lasers at the Hohlraum wall. Limiting these detrimental effects sets requirements on the tolerable magnitude of each effect. We present a nonlinear model for B-field diffusion through an infinitely long thin-walled cylinder with a temperature dependent resistivity, to show that a 15 μm thick wall of pure gold fails to meet these requirements. A new Hohlraum material made from an alloy of Au and Ta has a measured resistivity of ≥60 times that of Au and is shown with the nonlinear model to meet the requirements for magnetization. We compare the nonlinear model to simulations of the actual Hohlraum target using a finite element code which includes temperature-dependent Hohlraum resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗