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

Mesoscale modeling of microstructure-dependent thermal conductivity in U-Zr fuels

In uranium-zirconium (U-Zr) based metallic fuels, different phases can form at different compositions and temperatures. Typically, lamellar δ-UZr 2 and α-U phases are the dominant microstructures in U-rich U-Zr alloys at temperatures below 880 K. In this work, a finite element method based mesoscale modeling technique is used to calculate the effective thermal conductivities of such heterogeneous microstructures, using the thermal conductivities of two individual phases and their interphase thermal resistance (Kapitza resistance) as input parameters. The Kapitza resistance between δ-UZr 2 and α-U is determined at different temperatures, which shows an approximately T 3 dependence in the temperature range between 300 and 800 K. In addition, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr 2 phase. Further, an analytical model is therefore developed to quantify the effects of both temperature and δ-UZr 2 aspect ratio on the Kapitza resistance. Using this newly developed Kapitza resistance model, the effective thermal conductivities of a number of δ-UZr 2 + α-U heterogeneous microstructures in U-Zr alloys, including non-lamellar microstructures, can be estimated accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

In-situ effect in cross-ply laminates under various loading conditions analyzed with hybrid macro/micro-scale computational models

In this article, multi-scale finite element analyses based on three-dimensional (3D) hybrid macro/micro-scale computational models subjected to various loading conditions are carried out to examine the in-situ effect imposed by the neighboring plies on the failure initiation and propagation of cross-ply laminates. A detailed comparative study on crack suppression mechanisms due to the effect of embedded laminar thickness and adjacent ply orientation is presented. Furthermore, we compare the results of in-situ transverse failure strain and strength between the computational models and analytical predictions. Good agreements are generally observed, indicating the constructed computational models are highly accurate to quantify the in-situ effect. Subsequently, empirical formulas for calculating the in-situ strengths as a function of embedded ply thickness and different ply angle between embedded and adjacent plies are developed, during which several material parameters are obtained using a reverse fitting method. Finally, a new set of failure criteria for σ 22 -τ 12 , σ 22 -τ 23 , and σ 11 -τ 12 accounting for the in-situ strengths are proposed to predict laminated composites failure under multi-axial stress states. This study demonstrates an effective and efficient computational technique towards the accurate prediction of the failure behaviors and strengths of cross-ply laminates by including the in-situ effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling Isotope Separation in Electrochemical Lithium Deposition

Naturally occurring Li consists of two stable isotopes, 6 Li with an abundance of about 7.5%, and 7 Li making up the remainder with 92.5%. The development of a 6 Li enrichment technique, in terms of technical reliability and environmental safety to reach 6 Li future requirements, represents a key step in the roadmap for nuclear fusion energy supply worldwide. This paper uses finite element analysis-based models to simulate electrochemical Li isotope separation, which is an attractive method in terms of simplicity, safety, and scalability. In the model, we quantitatively analyze how different electrochemical factors including thermodynamics, charge-transfer kinetics, and diffusivities affect the separation process (separation factor), together with cell parameters, such as cell length and current density. The maximum separation factor of 1.128 could be obtained with the cell under the optimal thermodynamic, kinetic, and diffusive conditions, which is among the highest separation factors ever reported. Furthermore, these results will assist in designing the actual isotope separation setup with large separation factor and appropriate timing for sample collection.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Hierarchical model reduction driven by a proper orthogonal decomposition for parametrized advection-diffusion-reaction problems

This work combines the Hierarchical Model (HiMod) reduction technique with a standard Proper Orthogonal Decomposition (POD) to solve parametrized partial differential equations for the modeling of advection-diffusion-reaction phenomena in elongated domains (e.g., pipes). This combination leads to what we define as HiPOD model reduction, which merges the reliability of HiMod reduction with the computational efficiency of POD. Two HiPOD techniques are presented and assessed by an extensive numerical verification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2025

Parameter extraction for a SPICE model of an hTron superconducting thermal switch

Efficiently simulating large circuits is crucial to the development of superconducting nanowire-based electronics. However, current simulation tools for this technology are not adapted to the scaling of circuit size and complexity. We focus on the multilayered heater-nanocryotron (hTron), a promising superconducting nanowire-based switch used in applications such as superconducting nanowire single-photon detector readout. Previously, the hTron was modeled using traditional finite-element methods, which fall short in simulating systems at a larger scale. An empirical-based method would be better adapted to this task, enhancing both simulation speed and agreement with experimental data. In this work, we perform switching current and activation delay measurements on 17 hTron devices. We then develop a method for extracting physical fitting parameters used to characterize the devices. We build a SPICE behavioral model that reproduces the static and transient device behavior using these parameters, and validate it by comparing its performance to a model developed in prior work, showing an improvement in simulation time by several orders of magnitude. Furthermore, our model provides circuit designers with a tool to help understand the hTron’s behavior during all design stages, thus enabling broader use of the hTron across various new areas of application.

Caloritronics↗
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2023

An extended three-field principle to scale-bridge the granular micromechanics of polymer-bonded particulate materials

Observed phenomena of damage, viscoelasticity, and viscoplasticity in polymer-bonded particulate materials (PBPM), like polymer-bonded explosive (PBX), are accommodated in a scale-bridging constitutive model based on granular micromechanics. An extended three-field variational principle with min - max optimization problem is introduced to homogenize the bonded aggregates’ grain-scale contact mechanics. Similarly to prior applications of a three-field principle, volume constraints are invoked to embed microscale-to-macroscale correspondences for kinematic and static volumetric measures, like mean strain and its energy-conjugate stress invariant. Unlike prior applications, our three-field principle realizes a strain-driven model for effective constitutive behavior due to bonded intergranular contact. In conclusion, the resultant micro-sphere integration-based homogenization is demonstrated in comparison to existing publicly-available measurements to capture the extreme tension/compression constitutive asymmetry exhibited by PBPMs (specifically PBX-9501) during material softening.

42 ENGINEERING↗
DOE OSTI2022

New Hybrid Model for Evaluating the Frequency-Dependent Leakage Inductance of a Variable Inductance Transformer (VIT)

Skin and proximity effects can cause a significant drop in the effective leakage inductance of a transformer when the operating frequency is increased. Although the magnetic image method-based double-2-D model can calculate the lowfrequency leakage inductance with sufficient accuracy, it is inherently a frequency-independent model. While Dowell’s 1-D model uses frequency-dependent relations to account for both skin and proximity effects, its accuracy is severely affected by the assumed winding geometry. In this paper, a hybrid model is proposed that uses superposition to combine a modified Dowell’s model with the double-2-D model. The proposed model is investigated on a variable inductance transformer (VIT)—a partially-filled transformer whose leakage inductance can be varied by moving one of the windings mechanically. The frequency-dependent leakage inductances of the VIT evaluated using the hybrid model are in excellent agreement with the corresponding finite element method (FEM) simulated and experimentally measured values, thereby validating the proposed hybrid model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Sensitivity of additively manufactured AA7075 to variation in feedstock composition and print parameters

Here we developed a combined finite element and CALPHAD based model of the Laser Powder Bed Fusion (LPBF) process for AA7075 alloy that considers the effect of feedstock composition and print parameters. A single-pass of a laser on a layer of AA7075 alloy powder has been considered. Sensitivity of temperature evolution and melt pool geometry to variation in the stoichiometry of the feedstock powder and laser source characteristics have been studied. Our results indicate that deviation (up to 10%) of the feedstock composition from the AA7075 raises the maximum temperature and increases melt pool size. Excess Cu content shows the largest melt pool width and depth among all the cases. The peak temperature is higher than the standard feedstock composition in all cases, except when the Cu concentration is reduced. Increasing the scan power also results in a higher peak temperature and a larger melt pool size. Furthermore, the temperature's rise time increases by lowering the scan speed.

42 ENGINEERING↗
DOE OSTI2022

Enhanced critical axial tensile strain limit of CORC ® wires: FEM and analytical modeling

Abstract Conductor on Round Core (CORC ® ) cables and wires are composed of spiraled high-temperature superconducting (HTS) rare-earth barium copper oxide (REBCO) tapes, wound in multiple layers, and can carry very high currents in background magnetic fields of more than 20 T. They combine isotropic flexibility and high resilience to electromagnetic and thermal loads. The brittle nature of HTS tapes limits the maximum allowable axial tensile strain in superconducting cables. An intrinsic tensile strain above about 0.45% will introduce cracks in the REBCO layer of straight HTS tapes resulting in irreversible damage. The helical fashion at which the REBCO tapes are wound around the central core allows tapes to experience only a fraction of the total axial tensile strain applied to the CORC ® wire. As a result, the critical strain limit of CORC ® wires can be increased by a factor of more than 10 that of REBCO tapes. Finite element (FE) and analytical models are developed to predict the performance of CORC ® wires under axial tensile strain. A parametric analysis is carried out by varying the winding angle, the Poisson’s ratio of the CORC ® wire core, the core diameter, and the tape width. The results show that a small variation in winding angle can have a significant impact on the cable’s axial tensile strain tolerance. While the radial contraction of the helically wound tapes in a CORC ® wire under axial tensile strain depends on its winding angle, it is mostly driven by the Poisson’s ratio of the central core, affecting the tape strain state and thus its performance. Contact pressure from multiple layers within the CORC ® wire also affects the CORC ® wire performance. The FE model can be used to optimize the cable design for specific application conditions, resulting in an irreversible strain limit of CORC ® cables and wires as high as 7%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine-Learned Linear Structural Dynamics

The tension between accuracy and computational cost is a common thread throughout computational simulation. One such example arises in the modeling of mechanical joints. Joints are typically confined to a physically small domain and yet are computationally expensive to model with a high-resolution finite element representation. A common approach is to substitute reduced-order models that can capture important aspects of the joint response and enable the use of more computationally efficient techniques overall. Unfortunately, such reduced-order models are often difficult to use, error prone, and have a narrow range of application. In contrast, we propose a new type of reduced-order model, leveraging machine learning, that would be both user-friendly and extensible to a wide range of applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Generating MCNP Input Files for Unstructured Mesh Geometries

The Los Alamos National Laboratory’s (LANL) Monte Carlo N-Particle (MCNP)1 transport code version 6.3 (also known as MCNP6.3) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is time-consuming and error-prone as the complexities of geometries increase. A UM geometry model is a collection of finite elements representing a solid geometry. The first step of the MCNP UM calculation is using other software packages to create a finite element mesh representation of a solid 3D geometry because the MCNP code cannot be used to generate a UM model. Computer-aided design (CAD) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. Some mesh generation software packages may also be used to create solid geometries and thus CAD files are not needed. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software suite. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with a 6.3 version, the MCNP code can process HDF5 mesh input files. We only focus on the UM models formatted as Abaqus input files in this report since currently no external software can be used to generate HDF5 mesh input files for MCNP UM calculations. The MCNP code version 6.3 can be used to convert the Abaqus mesh input files into the HDF5 mesh input files, but this option is typically used by the MCNP code development team to test the HDF5 mesh input file feature. Several software packages (such as Abaqus, Attila4MC, or Cubit) can be used to create the Abaqus input files for MCNP UM calculations. An MCNP UM calculation using an Abaqus model requires two input file types: MCNP and Abaqus input files. The Abaqus input files needed for MCNP UM calcu lations must have the correct Abaqus syntax and meet the additional requirements by the MCNP code. The MCNP code can process only Abaqus input files that make use of part and assembly definitions, where elements in each part must be grouped into one or more element sets (i.e., elset) using *Elset keyword lines with specified naming formats. The MCNP and Abaqus input files required for MCNP UM simulations must be related; pseudo-cells in an MCNP input file must be constructed from mesh model data from an Abaqus input file. For large complex UM models, it is tedious to manually create MCNP UM input files. The um pre op (unstructured mesh pre operations) program with the -m option can be used to create a skeleton MCNP input file from an Abaqus input file [6]. Since the um pre op program was written in Fortran and was not written for optimized performance, this program is a deprecated feature in the MCNP code version 6.3 and may be removed in the next release of the code. To improve calculation flow of multiphysics calculations, a Python3 code called write mcnp um input has been developed to generate an MCNP input file instead of using the um_pre_op -m option. This Python code was initially released to the public in 2020. We have updated this Python code for MCNP6.3 and it was used to generate the MCNP input files used to verify the MCNP6.3 code. The write_mcnp_um_input code is included with the MCNP6.3 code package which will be released to the public through the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory. This report is a revision of LA-UR-20-27139 report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗