Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

Multidisciplinary Approach to Aerospike Nozzle Design

A model of a linear aerospike rocket nozzle that consists of coupled aerodynamic and structural analyses has been developed. A nonlinear computational fluid dynamics code is used to calculate the aerodynamic thrust, and a three-dimensional finite-element model is used to determine the structural response and weight. The model will be used to demonstrate multidisciplinary design optimization (MDO) capabilities for relevant engine concepts, assess performance of various MDO approaches, and provide a guide for future application development. In this study, the MDO problem is formulated using the multidisciplinary feasible (MDF) strategy. The results for the MDF formulation are presented with comparisons against separate aerodynamic and structural optimized designs. Significant improvements are demonstrated by using a multidisciplinary approach in comparison with the single-discipline design strategy.

Korte, J. J.↗

Status of the KTH-NASA Wind-Tunnel Test for Acquisition of Transonic Nonlinear Aeroelastic Data

This paper presents a status report on the collaboration between the Royal Institute of Technology (KTH) in Sweden and the NASA Langley Research Center regarding the design, fabrication, modeling, and testing of a full-span lighter configuration in the Transonic Dynamics Tunnel (TDT). The goal of the test is to acquire transonic limit-cycle- oscillation (LCO) data, including accelerations, strains, and unsteady pressures. Finite element models (FEMs) and aerodynamic models are presented and discussed along with results obtained to date.

Silva, Walter A.↗

Electrical four-point probing of spherical metallic thin films coated onto micron sized polymer particles

Micron-sized metal-coated polymer spheres are frequently used as filler particles in conductive composites for electronic interconnects. However, the intrinsic electrical resistivity of the spherical thin films has not been attainable due to deficiency in methods that eliminate the effect of contact resistance. In this work, a four-point probing method using vacuum compatible piezo-actuated micro robots was developed to directly investigate the electric properties of individual silver-coated spheres under real-time observation in a scanning electron microscope. Poly(methyl methacrylate) spheres with a diameter of 30 μm and four different film thicknesses (270 nm, 150 nm, 100 nm, and 60 nm) were investigated. By multiplying the experimental results with geometrical correction factors obtained using finite element models, the resistivities of the thin films were estimated for the four thicknesses. These were higher than the resistivity of bulk silver.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chromium Poisoning ERMINE Model Data

This is the data from the paper "Systematic and Predictive Trends to Chromium Poisoning in Solid Oxide Fuel Cell Cathodes" by Hokon Kim et al., appearing in the Journal of Power Sources. The files here are the microstructures, mesh files, and outputs from the ERMINE finite element model of chromium poisoning in SOFC cathodes. Further description can be found in Readme.rtf and in the paper..

Chromium,SOFC↗

Precise Determination of the Zero-Gravity Surface Figure of a Mirror without Gravity-Sag Modeling

The zero-gravity surface figure of optics used in spaceborne astronomical instruments must be known to high accuracy, but earthbound metrology is typically corrupted by gravity sag. Generally, inference of the zero-gravity surface figure from a measurement made under normal gravity requires finite-element analysis (FEA), and for accurate results the mount forces must be well characterized. We describe how to infer the zero-gravity surface figure very precisely using the alternative classical technique of averaging pairs of measurements made with the direction of gravity reversed. We show that mount forces as well as gravity must be reversed between the two measurements and discuss how the St. Venant principle determines when a reversed mount force may be considered to be applied at the same place in the two orientations. Our approach requires no finite-element modeling and no detailed knowledge of mount forces other than the fact that they reverse and are applied at the same point in each orientation. If mount schemes are suitably chosen, zero-gravity optical surfaces may be inferred much more simply and more accurately than with FEA.

measurements↗

Finite element meshing of ANSYS (trademark) solid models

A large scale, general purpose finite element computer program, ANSYS, developed and marketed by Swanson Analysis Systems, Inc. is discussed. ANSYS was perhaps the first commercially available program to offer truly interactive finite element model generation. ANSYS's purpose is for solid modeling. This application is briefly discussed and illustrated.

Kelley, F. S.↗

MODELING A SPENT NUCLEAR FUEL CASK SEISMIC TEST

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a canister dry storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental US to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. The earthquake conditions will represent seismic hazards in the 2,000-to-20,000-year return period range. The test will use instrumented pressurized water reactor fuel assemblies with surrogate mass inside the fuel rods instead of radioactive fuel pellets. Dummy assemblies with similar mass and dimensions to actual fuel assemblies will occupy the rest of the fuel assembly locations in a SNF canister that holds 32 fuel assemblies. The canister will be located inside a mockup vertical concrete cask. Instrumentation will record the motion of the major components of the complex dynamic system, and strain gauges will be used to record the cladding strain at select locations. Preparations for the test require modeling predictions to identify the range of response of the system and to help select specific earthquake cases to be simulated on the shake table from a large set of potential cases. This paper describes the pretest nonlinear finite element modeling efforts that have been completed to date, including cask system level modeling and fuel assembly modeling in LS-DYNA. The cask system level models are critical for anticipating sliding or tipping of an unanchored cask during the test. The fuel assembly model is needed to estimate the range of cladding strain response and fuel assembly structural response to be expected during the test and cladding strain measurements will be key metrics for model validation and the development of modeling best practices after the test is concluded.

Klymyshyn, Nicholas A.↗

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗

Thermal-vacuum testing and in-situ optical alignment measurements of the HALOE telescope/sun sensor assembly

Design features, test data and projected performance levels of the telescope/sun sensor assembly for the Halogen Occultation Experiment (HALOE) on board the Upper Atmosphere Research Satellite are described. HALOE will gather data on stratospheric chemical species with an IR telescope pointed through the atmosphere at the sun during occultation periods. The pointing accuracy will need to be 1 arcmin elevation and 1.6 arcmin azimuthal. Adjustments of the gimbaled instrument are guided by a sun sensor and its associated electronics. The components were subjected to thermal-vacuum, optical boresight stability tests and the data generated were compared with a finite element model of the telescope and sun sensor. The tests consisted of exposure to various thermal gradients while the telescope and sun sensor alignments were tracked by optical deflectometry. The thermal behavior was compared with predictions made with a finite difference model. Alignments were within tolerable ranges and the thermal behavior model was concluded valid for predicting the thermal behavior of orbiting instruments.

Foss, R. A.↗

Modeling used fuel response to normal conditions of transportation package drops to access geometric sensitivities

The US Department of Energy Used Fuel Disposition program is performing research to determine the loading conditions that used nuclear fuel (UNF) experiences during normal conditions of transportation (NCT). This loading environment includes mechanical shock and vibration as the UNF is transported by rail or truck, and it also includes potential package free drop conditions, which are the topic of this study. NCT drop test requirements for UNF package designs are defined in 10 CFR 71.71. The free drop test height is 30 cm for the UNF packages considered in this study. Explicit dynamic finite element models of generic rail and truck UNF packages containing a highly resolved 17x17 PWR fuel assembly are used to estimate the potential range of loading experienced by UNF during NCT horizontal package drops onto an unyielding surface. There are many variables and uncertainties that can affect the loading of UNF during this type of free drop scenario. This study focuses on geometric variables, such as gap conditions and small adjustments to the angle of impact, in a number of analysis cases. The LS-DYNA finite element analyses calculate detailed results for all the components of the fuel assembly. The results show that the peak localized axial strain in the fuel cladding is estimated to be 0.006421 (nominal) with an upper bound of 0.014511, based on the peak instantaneous bending moment of 63.65 N-m and other uncertainty factors. For the assumed low burnup condition of the UNF, the cladding yield strain is 0.011027 and the failure strain is 0.015574, which means that the calculated upper bound strain approaches the cladding failure limit without exceeding it. The peak shear force on the cladding was calculated to be 7879 N. This could be significant to cladding integrity if applied to the UNF in a situation where the fuel pellets do not support the shear transfer. Further evaluation is needed to determine if this is a credible threat to cladding integrity for these 30 cm free drop scenarios.

42 ENGINEERING↗

Modeling the Effect of Damage on the Electrical Resistivity of Melt-Infiltrated SiC/SiC Composites

Electrical resistivity (ER) measurements are a possible health monitoring technique for CMC components in future aerospace applications. In order to use ER measurements to detect and identify damage, it is necessary to understand how each specific damage state will affect the ER response. In this paper, finite element models are developed and applied to quantify the effect of specific damage states on the electrical resistivity response in a melt-infiltrated silicon carbide (SiC) fiber-reinforced silicon carbide composite. The electrical resistivity of several damage states are calculated by simulating the electric current flow through the damaged micro-structure. This is achieved by performing the numerical solution of the steady-state conservation of charge density equation. Numerical results reveal that cracking of the inter-tow matrix has the most profound effect on the composite electrical resistivity. Also, fiber/matrix debonding at matrix cracks in the 0° tows (tows aligned with the loading direction) may cause a significant increase in the electrical resistivity, but only if the fiber coating resistivity is 1000 Ω-cm or less. Cracks in the 90° tows and the crack opening displacement have very little effect on the composite electrical resistivity.

finite element analysis↗

Evaluation of Transverse Thermal Stresses in Composite Plates Based on First-Order Shear Deformation Theory

A postprocessing procedure is presented for the evaluation of the transverse thermal stresses in laminated plates. The analytical formulation is based on the first-order shear deformation theory and the plate is discretized by using a single-field displacement finite element model. The procedure is based on neglecting the derivatives of the in-plane forces and the twisting moments, as well as the mixed derivatives of the bending moments, with respect to the in-plane coordinates. The calculated transverse shear stiffnesses reflect the actual stacking sequence of the composite plate. The distributions of the transverse stresses through-the-thickness are evaluated by using only the transverse shear forces and the thermal effects resulting from the finite element analysis. The procedure is implemented into a postprocessing routine which can be easily incorporated into existing commercial finite element codes. Numerical results are presented for four- and ten-layer cross-ply laminates subjected to mechanical and thermal loads.

Rolfes, R.↗

Frequency dependence of the acoustic field generated from a spherical cavity transducer with open ends

Resolution of high intensity focused ultrasound (HIFU) focusing is limited by the wave diffraction. We have developed a spherical cavity transducer with two open ends to improve the focusing precision without sacrificing the acoustic intensity (App Phys Lett 2013; 102: 204102). This work aims to theoretically and experimentally investigate the frequency dependence of the acoustic field generated from the spherical cavity transducer with two open ends. The device emits high intensity ultrasound at the frequency ranging from 420 to 470 kHz, and the acoustic field is measured by a fiber optic probe hydrophone. The measured results shows that the spherical cavity transducer provides high acoustic intensity for HIFU treatment only in its resonant modes, and a series of resonant frequencies can be choosen. Furthermore, a finite element model is developed to discuss the frequency dependence of the acoustic field. The numerical simulations coincide well with the measured results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nondestructively Visualizing and Understanding the Mechano-Electro-chemical Origins of “Soft Short” and “Creeping” in All-Solid-State Batteries

All-solid-state Li-metal batteries (ASLMBs) represent a significant breakthrough in the quest to overcome limitations associated with traditional Li-ion batteries, particularly in energy density and safety aspects. However, widespread implementation is stymied due to a lack of profound understanding of the complex mechano-electro-chemical behavior of Li metal in the ASLMBs. Herein, operando neutron imaging and X-ray computed tomography (XCT) are leveraged to nondestructively visualize Li behaviors within ASLMBs. This approach offers real-time observations of Li evolutions, both pre- and post- occurrence of a “soft short”. The coordination of 2D neutron radiography and 3D neutron tomography enables charting of the terrain of Li metal deformation operando. Concurrently, XCT offers a 3D insight into the internal structure of the battery following a “soft short”. Despite the manifestation of a “soft short”, the persistence of Faradaic processes is observed. To study the elusive “soft short”, phase field modeling is coupled with electrochemistry and solid mechanics theory. The research unravels how external pressure curbs dendrite growth, potentially leading to dendrite fractures and thus uncovering the origins of both “soft” and “hard” shorts in ASLMBs. Furthermore, by harnessing finite element modeling, it dive deeper into the mechanical deformation and the fluidity of Li metal.

25 ENERGY STORAGE↗

Sensor Development for Liquid Water Detection in Dry Storage Casks: FY 2023 Status

Modeling efforts were undertaken in fiscal year (FY) 2023 to evaluate the feasibility and capability for sensing of liquid water inside canister-based dry cask storage systems (DCSSs). The focus was on the development of full-scale finite element models (FEMs) of ultrasound propagation through DCSS canister components. For these initial investigations, simulation of the baseplate component was targeted, envisioning water collection at the bottom of a vertically oriented canister. The effects of internal canister components (i.e., fuel basket, fuel assemblies) on ultrasound propagation in the baseplate component were evaluated in these efforts. The environment inside a DCSS confinement is intended to be inert and free of water to prevent potential corrosion of used fuel cladding or other internal hardware. Spent fuel assemblies are dried, after storage in water pools, to make sure water has been removed from assembly cavities. However, there is some uncertainty about the amount of residual water potentially left behind in a DCSS after drying processes, because water can become trapped in cavities or other small crevices in the surfaces formed by the fuel cladding, fuel assemblies, and other internal hardware components. Considering the complex space- and time-dependent temperature profiles in DCSSs, water may be in a liquid or gas phase depending on its location in the cask and how long the cask has been in storage. Evacuating most water and oxidizing agents contained within a canister is recommended by NUREG-1536 (NRC 2010), which covers DCSSs. As summarized by Salazar et al. (2020), existing guidance typically relies on achievable vacuum pumping pressures sustained over a hold period as a signal of dryness and water removal. However, time series data about water removal from full-scale commercial drying procedures are lacking (Hanson and Alsaed 2019). A review of drying specifications from several vendors led to the conclusion that if the specifications are followed correctly, the residual moisture in DCSSs should present an insignificant risk of cladding degradation (Knoll and Gilbert 1987). A more recent analysis concluded that much larger quantities of residual water could remain in DCSSs, but the amount would still not be expected to lead to significant corrosion of fuel cladding or other internal components (Jung et al. 2013). Industry drying procedures are mostly prescriptive in nature, and operational issues arising during the process could result in incomplete drying. A summary of operational issues and potential negative effects of residual water is provided by Salazar et al. (2019). Experimental efforts are ongoing to validate the extent of water removal in a DCSS based on drying procedures used at nuclear power plants through well-designed investigations of drying process efficacy and water retention (Durbin et al. 2021; Pulido et al. 2022a, Pulido et al. 2022b). This has been approached by simulating limited portions of a DCSS internal volume and fuel assemblies. So far, these efforts have focused on the effects of potential water trapping in the dashpot region of control rod guide tubes but are expected to be expanded to include the effects of other internal hardware features and failed fuel rod cladding. A method for detection and measurement of liquid water (Meyer et al. 2022; Meyer et al. 2021), in tandem with the drying information gained through experimental investigations, provides comprehensive bases for understanding the internal conditions of DCSSs and provides technical information to support licensing decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High Fidelity Modeling of Evolutionary Structures in IMOS

A methodology is proposed for testing, modal identification, and model correlation of structures which are assembled in an evolutionary manner. The methodology prescribes initial testing of a base structure without any attached components. The experimental modal properties are used to update the physical properties of the base structure finite element model using a Bayesian Estimation Technique (BET). The BET is implemented in METLAB using IMOS functions. As components are added to the base structure, additional tests are performed. Three testing options are proposed, all of which enable estimation of the physical parameters of the components alone. Furthermore, improvements to the original BET implementation have been made, which contribute to higher fidelity models and more efficient computational times.

structures testing assembling modeling BET Bayesia↗

Mixed Models and Reduction Techniques for Large-Rotation, Nonlinear Analysis of Shells of Revolution with Application to Tires

An effective computational strategy is presented for the large-rotation, nonlinear axisymmetric analysis of shells of revolution. The three key elements of the computational strategy are: (1) use of mixed finite-element models with discontinuous stress resultants at the element interfaces; (2) substantial reduction in the total number of degrees of freedom through the use of a multiple-parameter reduction technique; and (3) reduction in the size of the analysis model through the decomposition of asymmetric loads into symmetric and antisymmetric components coupled with the use of the multiple-parameter reduction technique. The potential of the proposed computational strategy is discussed. Numerical results are presented to demonstrate the high accuracy of the mixed models developed and to show the potential of using the proposed computational strategy for the analysis of tires.

Noor, A. K.↗

A variable stiffness strut using strain induced bending for vibration suppression in the evolutionary model

A new piezoceramic strut has been developed. This strut utilizes strain induced bending as a means of actuation. The significance of this design is that the displacement of the strut's motion arm is independent of the expansion of the piezoceramic material. Design criteria was developed using continuous and finite element modeling techniques. Experimental and analytical results were compared. Finally, a numerical analysis was conducted that demonstrated the benefits of the piezoceramic strut configuration.

Pokines, Brett↗