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 955 records · Page 53

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

OpenTurbine [SWR-23-07]

OpenTurbine is an open-source code designed for the simulation of wind turbine structural dynamics. It is a flexible multi-body dynamics solver wherein the wind turbine is modeled as a collection of beam finite elements, rigid bodies, and constraints. OpenTurbine is equipped with an application programming interface primarily designed for coupling to computational fluid dynamics codes, but it is also suited for coupling to aerodynamics solvers based on blade element momentum theory.

Sprague, MichaelA.↗

A Monte Carlo Rectangle Packing Algorithm for Identifying Likely Spatial Distributions of Final Closure Cap Subsidence in the E-Area Low-Level Waste Facility

A Monte Carlo rectangle packing algorithm, summarized in Figure ES-1, has been developed to provide insights into the likely spatial distribution of subsided areas or holes in the final closure cap for the E-Area Low-Level Waste Facility (ELLWF) resulting from the placement of non-crushable waste containers within slit and engineered trench disposal units. Using historical records provided by Solid Waste Management as an indicator of future waste streams, a probability distribution of unique waste containers was constructed, and Monte Carlo sampling was used to sequentially select waste containers for placement in a rectangular region bounded by the disposal unit footprint. Within this Monte Carlo approach, the probability of any container type being selected is weighted by the fractional abundance relative to all other container types. For each Monte Carlo realization, the frequency of non-crushable waste containers appearing at each location on a spatially discretized mesh is updated. After many realizations, the location(s) with the highest frequencies of non-crushable waste containers are regarded as the most likely locations to experience subsidence (i.e., for a hole to form). Utilizing this approach has enabled the down-selection of a computationally intractable set of subsidence hole combinations to be explored in three-dimensional finite element groundwater flow and contaminant transport models. Such an approach can be applied to a number of different packaging problems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Impact of Control Blade Insertion on the Deformation Behavior of SiC-SiC Channel Boxes in Boiling Water Reactors

This report 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 temperature and neutron flux dependent irradiation induced swelling of SiC, the SiC-SiC composite channel box can undergo significant 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 thermalhydraulics 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 the channel box will undergo time dependent lateral bowing for all the control blade positions in the assembly. The bowing behavior is dominated by the swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The lateral bending will cause a temporary interference between the channel box and the control blade. The developed stresses exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the ultimate tensile strength of the material, and therefore, development of a full, through thickness crack in the channel box is not expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Initial Assessment of the Design Margins of Different ASME Section III, Division 5 Design Rules

This report develops a method for assessing the design margin of the ASME Section III, Division 5 Class A design rules, focusing on a definition of margin as the ratio between the actual, expected component life and the ASME design life. Full inelastic finite element simulations with a complete damage model capturing the available creep, creep-fatigue, fatigue, and tension test failure data produce the expected component lives, which can then be compared to corresponding ASME design calculations. The focus of this particular work is on Alloy 617 and the ASME creep-fatigue design rules, but the approach is general and could be applied to other design limits and to other materials. A margin assessment of a representative Alloy 617 component shows the ASME creep-fatigue design rules applied to this component have a margin of 10 on design life. The report describes several challenges in completing a more comprehensive margin assessment of the ASME rules, in particular challenges in developing probabilistic assessment methods for high temperature structural components. The report details these challenge and discusses the immediate possibilities of the deterministic margin assessment method developed here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental Study of Mechanical Response of Polyurethane Rubber at Varying Strain Rates for W80-4 System

This report describes an experimental study to determine the mechanical behavior of the polyurethane rubber material that was used in the W80-4 MCC Shock/Breach Phase 1 test series. Compression experiments were conducted on cylindrical specimens over a wide range of loading rates to characterize the material over the range of strain rates that were experienced in the shock/breach testing. Additionally, specimen diameter was varied to determine the effect of confinement on the material response and was found to be significant. This data is used to populate a material model to enable accurate analyses and finite element simulations of the shock/breach test series.

36 MATERIALS SCIENCE↗

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗

Power Electronics Materials and Bonded Interfaces - Reliability and Lifetime

The emergence of wide-bandgap (WBG) devices in power electronics packages has brought significant attention to the reliability and thermal performance of the bonded interfaces. These interfaces, such as die-attach and substrate-attach, need to withstand the severe stress induced on them by the high-temperature (> 200°C) operation of the WBG devices. Sintered silver and transient liquid phase alloys are key bonded interfaces that can enable the high-temperature operation of WBG devices, however it is essential to evaluate their reliability and study the possible failure mechanisms before they can be used in power electronics packages. Sintered silver is already being adopted in commercial packages as a die-attach, but large-area substrate-attach offers additional reliability challenges. In this project, we investigate the reliability of these interface materials under accelerated thermal cycling and correlate these results with finite element simulations to develop lifetime prediction models. The failure mechanisms of these bonded interfaces are also studied.

47 OTHER INSTRUMENTATION↗

Adaptive time scaling for accelerating explicit finite element analysis

A method for accelerating an explicit finite element analysis (FEA) simulation of a modeled system or process includes performing an initial iteration of the FEA simulation according to a baseline time interval via an FEA computing network, and calculating a criteria ratio of a predetermined set of scaling criteria for the modeled system or process. The method includes determining a time-scaling factor using the criteria ratio via the FEA computing network as a function of the criteria ratio, and then applying the time-scaling factor to the baseline time interval to generate a scaled time interval. The scaled time interval accelerates simulation time of the FEA simulation. The method includes performing a subsequent iteration of the explicit FEA simulation at the scaled time interval using the FEA computing network. The process continues for subsequent iterations, with the time-scaling factor adapting with each iteration.

Chen, Jian↗