Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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

Improvements in High Temperature Gas Cooled Reactor Modeling Capabilities in the Pronghorn Code

This report details the improvement of pebble bed reactor modeling capabilities in the Pronghorn code in fiscal year 2022. The following accomplishments are reported: Deployment of weakly compressible finite volume formulation to the HTR- PM reference plan model; Enable modeling of stagnant gas gaps in the finite volume formulation; Enable using all Pronghorn correlations available in the finite element version in the finite volume version; Modeling of decay heat in pebble bed reactors; Simplifying the input for multiphysics equilibrium core calculations and significant reduction of execution time; Implementation of advanced correlations developed by the Center of Excellence for Thermal-Fluids Applications in Nuclear Energy . In addition, this report includes a development plan for Pronghorn and associated NEAMS tools for prismatic gas-cooled reactors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Verification and Regression Testing of a Physically Stabilized Layered Solid Element Formulation in DYNA3D/Paradyn

Recent modeling of filament-wound composite structures drove the need for a new element type within the finite element code DYNA3D/Paradyn. This new layered solid element was implemented and designed to capture the kinematics and constitutive be havior of various lamina layers defined with arbitrary orientations and volume fractions within the element to accurately model a laminated composite material. The new el ement uses a single integration point in each of the lamina layers defined to capture the constitutive behavior, and a novel physical stabilization routine is used to prevent hourglassing while mitigating shear and volumetric locking. This report documents the verification testing and the regression tests added to the DYNA3D/Paradyn software quality assurance test suite to assess the proper implementation of the new layered solid element. There were a total of 27 tests added the DYNA3D test suite, which consist of problems using the layered solid element with a single layer or multiple layers for isotropic and orthotropic material models. These tests are a combination of simple kinematically driven patch tests, beam bending tests, plate bending tests, and more complicated problems. Ultimately, the testing done verifies that the element behaves as expected and is implemented correctly.

42 ENGINEERING↗
DOE OSTI2022

Pronghorn Porous Media Model Validation with Pressure Drop Measurements

The verification and validation (V&V) of Pronghorn is imperative to assert its accuracy when predicting the fluid velocity, temperature, and pressure in high temperature gas-cooled reactors. Pronghorn is a coarse-mesh, intermediate-fidelity, and multidimensional thermal-hydraulics (TH) code developed by the Idaho National Laboratory (INL). New pebble bed experiments are used to observe the details of the fluid motion and pressure drop in the porous bed under the reactor normal operation. This paper focuses on the validation of the Pronghorn compressible and incompressible Navier-Stokes equations using the pressure drop measurements performed at the engineering-scale pebble bed facility at the Texas A&M university (TAMU). Various pressure drop correlations and porosity functions are implemented in both Pronghorn and STAR-CCM+ to compare the pressure drop due to the combined viscous and inertial resistances in the porous bed. The correlations accounting for the near-wall effect are also utilized to observe if the pressure drop estimates can be improved. Pronghorn porous media models predict the pressure drop well relative to the STAR-CCM+ simulation results and 1D correlations, and both the finite element method (FEM) and finite volume method (FVM) perform accurately. Pronghorn models are also validated with the experimental measurements given the different Reynolds number ranges and specific aspect ratios. The likelihood of the statistical significance between the pressure drop measurements and specific correlations or simulations is low provided that the overlap of their confidence intervals is more than the half of a single arm. Several validation metrics are reasonable in regard to the similar studies from other literature. The precise average pebble bed porosity estimation has much impact on the pressure drop, and the Foumeny and Montillet (dense packing) models carry out the accurate pressure drop prediction by considering the near-wall effect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2021

Synchrotron experiment and simulation studies of magnesium-steel interface manufactured by impact welding

The effective weight reduction in the automotive industry by the wide adoption of lightweight magnesium (Mg) alloys demands high-quality joint between magnesium alloys and massively-used steels in order to wring the excess weight with strength and safety assurance. However, Mg-steel joint is difficult to achieve because there is no mutual solubility between magnesium and steel and huge disparity in physical properties. An impact-based welding method recently showed successful Mg-steel joining. In this work, the characteristics of Mg-steel interface joined by the impact welding method were investigated. Synchrotron high-energy X-ray computed tomography and diffraction were applied to characterize the microstructure across Mg-steel interface. Results revealed a deposit layer formed at the joint interface where Fe-rich particles spread deep into the Mg matrix. High-resolution 3D morphology of Mg-steel interface demonstrated the trapped pores and cracks inside the deposit layer. Finally, the formation of the deposit layer and the void/cracking evolution were analyzed by using finite element models. These findings provide insights into the immiscible Mg-steel joining process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal conductivity of crack-containing media: A numerical study

The present study aims to evaluate the effective thermal conduction behavior of crack-containing media through micromechanical modeling. Numerical analyses were performed using the finite element method, first using periodic arrays of elliptical pores and then reducing the elliptical minor axis to reach the crack limit. This parametric approach was seen to be able to capture the effect of discontinuity across the crack face in a straightforward manner. The overall thermal conductivity of the model structures, encompassing a range of crack densities and spatial distributions, was investigated. It was found that, for cracks with mixed orientations, the effective thermal conductivity is insensitive to the actual crack alignment and can be uniquely represented by the crack density. The mixed-orientation periodic crack configurations are also representative of random sizes, orientations, and spatial distributions of cracks, thus exhibiting a high degree of generality. The modeling analyses can also be used to determine simple empirical expressions for two- and three-dimensional media containing random cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Digital force prediction for milling

This paper presents a fully digital, integrated approach for milling force prediction with arbitrary end mill-work material combinations. The approach includes: 1) structured light scanning to identify the end mill’s cutting edge macro-geometry along the tool axis; 2) structured light scanning to measure the cutting edge cross-sectional rake and relief profiles; 3) finite element analysis of orthogonal cutting to determine the force model coefficients that relate the force to chip area using the work material’s constitutive model and tool’s rake and relief profiles; and 4) time domain simulation with inputs that include the measured cutting edge macro-geometry, finite element-based force model, and measured structural dynamics. Milling force predictions are compared to in-process measurements to validate the method.

42 ENGINEERING↗
DOE OSTI2022

Development of a Subchannel Capability for Liquid-Metal Fast Reactors in Pronghorn

This report details the development and demonstration of an entirely new capability in Pronghorn, namely the ability to model liquid-metal fast reactor (LMFR) flow conditions on the engineering scale. We developed two modeling approaches for LMFR that can be used separately or be combined into a hybrid simulation: (1) a modern subchannel capability called Pronghorn-Subchannel for square and hexagonal lattices, and (2) a porous flow capability for LMFR geometries. The report emphasizes the novel aspects of the developed subchannel capability and the interoperability of the subchannel capability, porous flow capability, and multiphysics tools within the multiphysics object oriented simulation environment (MOOSE). Here we demonstrate the ability to: (1) Accurately model subchannel flow in hexagonal lattices; (2) Couple the subchannel flow model to multidimensional finite-element method (FEM) or finite-volume method (FVM) heat conduction models; (3) Model LMFRs using Pronghorn’s porous media FVM approach; (4) Couple porous flow FVM and subchannel models in a single simulation; (5) Explicitly model inter-wrapper flows along with conjugate heat transfer from the intra-element flow; and (6) Demonstrate the numerical robustness of the subchannel algorithm by simulating intra-element flow recirculation in a high-buoyancy, low-flow fuel element.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Analysis of Element Size and Element Type on SS304L Material Model Performance

Finite element solutions are only valid if the information inputted into the models is a valid assumption. This is why it is important to determine if the model behaves as expected especially in the material response. Material models that consider deformation in the elastic and plastic regions are desirable. As a material deforms plastically, hardening will occur until the ultimate strength is reached due to the reduction of the cross sectional area in a uniaxial tensile test environment. A phenomenological approximation of this behavior which is known as power law hardening (plh) can be expressed as σ = σ y + K $ϵ^{n}_{p}$ where K and n are empirically determined constants, and σ y is the proportional limit of the material. The first objective of this study is to verify the behavior for a calibrated plh material model for selective deviatoric hexahedral elements when applied to a different material stock. The Sierra SM input for this calibrated model can be found in the appendix. The second objective of this study is to determine the differences in failure characteristics for considerably large shell elements versus that of fine hexahedral elements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cladding Profilometry Analysis of Experimental Breeder Reactor-II Metallic Fuel Pins with HT9, D9, and SS316 Cladding

BISON finite element method fuel performance simulations were conducted using an existing automated process that couples the Fuels Irradiation & Physics Database (FIPD) and the Integral Fast Reactor Materials Information System database by writing input files and comparing the BISON output to post-irradiation fuel pin profilometry measurements contained within the databases. The importance of this work is to demonstrate the ability to benchmark fuel performance metallic fuel models within BISON using Experimental Breeder Reactor-II fuel pin data for a number of similar pins, while building off previous modeling efforts. Changes to the generic BISON input file include implementing pin specific axial power and flux profiles, pin specific fluences, frictional contact, and irradiation-induced volumetric swelling models for cladding. A statistical analysis of irradiation-induced volumetric swelling models for HT9, D9, and SS316 was performed for experiments X421/X421A, X441/X441A, and X486. Between these three experiments, there were 174 post-irradiation examination (PIE) profilometries used for validating the swelling models presented using a standard error of the estimate (SEE) method. Implementation of the volumetric swelling models for D9 and SS316 claddings was found to have a significant impact on the BISON profilometry simulated, where HT9 clad pins had an insignificant change due to low fluence values. BISON profilometry simulated for HT9, D9, and SS316 fuel pins agreed with PIE profilometry measurements, with assembly SEE values being 4.4 × 10−3 for X421A, 2.0 × 10−3 for X441A, and 2.8 × 10−3 for X486. D9 clad pins in X421/X421A had the highest SEE values, which is due to the BISON simulated profilometry being shifted axially. While this work accomplished its purpose to demonstrate the modeling of multiple fuel pins from the databases to help validate models, the results suggest that the continued development of metallic fuel models is necessary for qualifying new metallic fuel systems to better capture some physical performance phenomena, such as the hot pressing of U-Pu-Zr and the fuel cladding chemical interaction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2025

Prelude to Simulation Diagnostics in MOOSE: Diagnosing Meshes

As finite element simulations increase in complexity, so does the need for high quality meshes. Without proper diagnostics, issues pertaining to poor mesh quality or improper mesh construction may go undetected, causing simulations to fail unexpectedly. This paper details the new mesh diagnostic system available in MOOSE, which provides users with the tools to inspect mesh quality for both externally and internally generated meshes. This system can detect issues arising from element volume size, non-conformal elements, intersecting edges, missing sidesets, as well as other common mesh issues. As meshes increase in size and complexity, such issues become more difficult to diagnose through visual inspection alone.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

High-Order Mesh hr-adaptivity for Surface Fitting to Implicit Geometries

We present an ℎ𝑟-adaptivity framework for morphing a given mesh to fit a target surface prescribed as the zero isocontour of a discrete function. In this framework, high-order meshing is posed as a variational minimization problem that depends on the mesh quality prescribed via the target matrix optimization paradigm (TMOP) and position of a subset of mesh nodes with respect to the target surface. The proposed formulation ensures that the variational problem is converged and mesh quality degradation near the surface is limited, even when the mesh topology is incompatible with the target surface. Additionally, a mesh subset-based approach and ℎ-refinement is introduced to efficiently increase fitting accuracy while reducing the computational cost of the mesh morphing problem. The ℎ𝑟-adaptivity technique extends to different element types in two- and three-dimensions, and can be used in existing finite element and spectral element frameworks to obtain high-order body-fitted meshes. Various numerical experiments demonstrate the robustness and accuracy of the fitting approach for problems of practical interest such as Lagrangian hydrodynamics and topology optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nondestructive Evaluation of Stress Corrosion Cracking in a Welded Steel Plate Using Guided Ultrasonic Waves

Stress corrosion cracking (SCC) had occurred in early-generation high-level nuclear waste tanks constructed by welding carbon steel. This paper describes an ultrasonic inspection system and its fundamental ability to detect and quantify the length of SCC on thick welded steel plates. The finite element method (FEM) was applied to simulate the welding process to estimate the welding residual stress field. Growth of stress corrosion cracks is driven by crack stress intensities exceeding the subcritical cracking threshold intensity. The subject plate was experimentally inspected with ultrasonic nondestructive evaluation (NDE) techniques to characterize the extent of SCC. The NDE system uses a piezoelectric transducer to generate guided waves in the thick steel plate, and a scanning laser Doppler vibrometer (SLDV) to measure multidimensional time–space wavefield data over a user-defined scanning area in the plate surface. The measured wavefield data can show wave interactions in a localized area in the plate due to the presence of the discontinuities of the SCC. To generate an inspection image that can precisely show the crack’s location and/or the dimension, the wavefield data are further processed to generate inspection image that maps the entire sample plate so the crack can be clearly identified in the plate while its length can be readily estimated. In conclusion, the ultrasonic test results for crack length agree well with the visually estimated length and are close to that predicted by the FEM for cracks in the weld residual stress field.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗