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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
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

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗
DOE OSTI2024

Simplified inelastic constitutive models for ASME Section III, Division 5 design by inelastic analysis

This report describes the development of simplified, universal constitutive model that captures the high temperature monotonic and cyclic behavior of a range of commonly-used high temperature materials. The goal of the work is to provide a simple, universal constitutive model to replace the current bespoke models for Grade 91, 316H, and Alloy 617 included in Nonmandatory Appendix HBB-Z of the ASME Boiler & Pressure Vessel Code, and to extend this model to cover Alloy 800H. We initiated this work in response to feedback from reactor vendors and other Code users requesting simplified models, compared to the current models, that are easier to implement and use in commercial finite element analysis software. This report describes the completion of this effort by developing a model to correct the defects in standard model forms presently used for high temperature material modeling, described in past work, developing and implementing new numerical methods to train this model against test data, and then actually training the model for the four materials. The report provides a complete mathematical description of the model along with the tabulated material coefficients for the four materials. The final step will be to formulate an ASME Code change to introduce the new models into the Code.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A general framework for substructuring-based domain decomposition methods for models having nonlocal interactions [including apppendix B]

A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

BISON TRISO Modeling Advancements and Validation to AGR-1 Data

BISON is a finite element-based nuclear fuel performance code. Among its unique characteristics are its ability to model 1D, 2D, and 3D geometries and its applicability to a wide variety of nuclear fuels. For the last eight years, BISON has included a beginning capability to model tri-structural isotropic (TRISO) fuel. Recently, interest in TRISO fuel has grown, and a significant effort has been made to improve BISON’s capabilities in this area. Capability development has occurred for each material present in TRISO fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. New elastic, creep, swelling, thermal expansion, thermal conductivity, and fission gas release (FGR) models are available. New models for the graphite matrix are also now available. Another important addition is the ability to perform statistical failure analysis of large samples of fuel particles. This new capability, which continues to grow, enables evaluation of failure due to pressure or crack formation by analyzing many thousands of particles. This enables realistic calculations of fission product release from the many particles in a TRISO-fueled reactor. These capabilities were checked via regression and verification tests. A large number of code benchmarking problems were also run, showing that BISON’s results closely match those of other software tools. Finally, a significant validation effort was completed in which fission product release, measured as part of the AGR-1 capsule experiments, was compared to BISON outputs. BISON outputs compared very well to the experimental data and to PARFUME results. Interest in BISON’s TRISO capabilities is growing, with the U.S. Nuclear Regulatory Commission (NRC) and Westinghouse Electric Company receiving training during the past year. Multiple other entities have expressed interest in or are actively using BISON. Kairos Power, LLC, has a strong partnership with Idaho National Laboratory (INL) regarding the use of BISON for TRISO analysis. While its capabilities still continue to grow, BISON has already become a powerful tool for TRISO analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗
DOE OSTI2022

Deep Learning-Assisted Real-Time Forward Modeling of Electromagnetic Logging in Complex Formations

Higher dimensional (i.e., 2-D and 3-D) modeling is indispensable to correctly evaluate the responses of electromagnetic logging tools in complex formation environments. However, limited by the high computational cost of rigorous modeling such as finite difference method and the finite element method, the real-time applications in the well logging industry primarily rely on the 1-D forward solver, which would result in erroneous formation evaluation for complex scenarios. As a result, aiming at realizing fast modeling for electromagnetic logging tools in complex formations, this paper proposes a general framework assisted by deep neural networks (DNNs). The framework consists of three modules: earth model classification, parameter extraction, and surrogate construction. Separate DNNs are trained and tested for different modules. The accuracy and efficiency of the DNN assisted fast modeling are validated by several experiments. Here, this study finds that the fast modeling assisted by DNNs is able to calculate the tool responses and reconstruct the subsurface formations in real-time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Linear viscoelastic characterization of electrically conductive adhesives used as interconnect in photovoltaic modules

Electrically conductive adhesives (ECAs) are incorporated into recent designs of photovoltaic (PV) modules and replace the traditional metallic solders as interconnects. This transition depicts a significant material change, and a proper understanding of the interconnects' mechanical response has not yet been established. However, such an understanding is necessary to (a) identify the driving forces for module degradation and failure, (b) allow for module design optimization, and (c) enable accurate lifetime predictions. This study summarizes the framework for the mechanical materials characterization and modeling of ECAs for PV applications. Only high-fidelity material models are able to capture the rate and temperature dependency of the ECA interconnect and allow for accurate modeling of the materials response. Furthermore, a linear viscoelastic representation is found to describe the mechanical response of the ECAs sufficiently well. The effects of curing conditions and environmental exposure are investigated, and material models for a variety of ECAs are reported and prepared for the use in numerical simulations. Overall, a finite element simulation of a generic submodel of a shingled cell module is used to highlight the need for high-fidelity material models and demonstrates the error made in the predicted stress states by using less sophisticated models.

14 SOLAR ENERGY↗
DOE OSTI2022

Cubit v.16.X

SAND2021-3051 O CUBIT is a full-featured software toolkit for robust generation of two- and three-dimensional finite element meshes (grids) and geometry preparation. Its main goal is to reduce the time required to generate meshes - particularly large hex meshes of complicated, interlocking assemblies. It is a solid-modeler-based preprocessor that meshes volumes and surfaces for finite element analysis. CUBIT also includes state-of-the-art smoothing algorithms and provides an extensive suite of tools for geometry decomposition and mesh generation. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Stimpson, Clint↗
DOE OSTI2022

Generating a debris model for a structure

A computer generates a computer model based on the output data of a finite element analysis (FEA) performed on a structure that has experienced a break up event. The output data is processed to generate a debris model. The debris model comprises data defining the structure, as well as any fragments or pieces of debris caused by forces exerted on the structure during the break up event. The debris model can then be utilized as input data to generate other computer models.

Laverty, Richard Rion↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Room Temperature Stamping of High-Strength Aluminum for Lightweight Structural Automotive Components

The use of aluminum (Al) 7xxx series in automobile structural applications is limited by factors such as its low room-temperature formability in T6 temper and the need for off-line artificial aging to regain high strength if it is formed in soft temper. In this study, two tempers were established that enable Al 7075 blanks to be room temperature stamped to create side impact beams that achieve near-T6 strength without post-forming artificial aging. It is shown that natural aging at room temperature, after solution treatment and before stamping, provides an alternative to high temperature pre-aging and hot stamping techniques in 7xxx alloys. Both tempers enabled room temperature stamping of the beams without cracking, indicating good formability. The Al beams were ~38% lighter than the benchmark boron steel beam. Further, the hardness in the low-deformation regions of the stamped Al beams were within 88% of a conventional 7075-T6 sheet. Paint-bake treatment of the beams appeared to cause over aging, lowering the ultimate tensile strength and ductility. Finite element three point bending simulations of the beams showed higher bending force and energy absorption before failure for the Al 7075 showing good promise for 7xxx series Al as a light-weight alternative to steel beams.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Novel processing route for design and manufacturing of metal toughened nanoceramics: Al–Al 2 O 3 nanocermets

Monolithic alumina is a highly brittle ceramic with a very low fracture toughness, making it an unreliable material for mechanical and structural applications. Ceramic-metal composites (cermets) are the subject of intensive interest to obtain materials which can address the intrinsic brittleness of ceramics without compromising their high strength and light weight. In the present investigation, aluminium-alumina (Al–Al 2 O 3 ) nanocermets were developed. The nanocermets possessed the toughness of Al phase and the strength of reinforced nano Al 2 O 3 particles. Here, the nanocermets green bodies were prepared such that upon sintering each adjacent Al 2 O 3 grain was conformally coated with a continuous Al film. Mechanical milling via a combination of cryo-milling and ball-milling was used to prepare Al–Al 2 O 3 nano-powder mixture. Spark plasma sintering was used to consolidate the green nanocermets. Results show that specimens prepared from 35 h milled powders exhibited superior mechanical properties as compared to other nanocermets, with an improvement in microhardness, elastic modulus, compressive strength, and flexural strength by 88%, 78%, 69% and 40% with respect to sintered pure Al. The indentation toughness of 35 h milled nanocermet showed an increment of 12.5% with respect to sintered pure Al 2 O 3 . Characterization techniques such as XRD, SEM, HR-TEM, X-ray tomography were employed to investigate the phase evolution and structural morphology, while FEM simulations were performed for understanding the crack propagation mechanisms in nanocermets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plasticity and fracture behavior of Inconel 625 manufactured by laser powder bed fusion: Comparison between as-built and stress relieved conditions

In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of control blade insertion on the deformation behavior of SiC-SiC channel boxes in BWRs

This article describes the analysis of distortion of a silicon carbide fiber-reinforced, silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of the temperature and neutron flux dependent irradiation-induced swelling of SiC, the SiC-SiC composite channel box can undergo distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermal-hydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that large gradients in fast neutron flux (up to 35–40% across a single axial level) will develop across the channel box. Due to these gradients, the channel box will undergo time-dependent bending for all the control blade positions in the assembly. The time-dependent bowing behavior is dominated by the transient swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The bending will cause temporary interference between the channel box and control blade, and the interference is expected to be most severe for the fully inserted control blade position. The developed stresses due to differential swelling in the channel box exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the tensile strength of the material, and therefore, development of a full, through-thickness crack in the channel box is not expected. Finally, further work is recommended to explore and evaluate the mitigation strategies.

42 ENGINEERING↗