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 793 records · Page 44

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗

Modeling Space Radiation Induced Bone Changes in Rat Femurs through Finite Element Analysis

As the duration of manned missions outside of the Earth’s protective shielding increase, astronauts are at risk for exposure to space radiation. Various organ systems may be damaged due to exposure. This study investigates the bonestrength changes using finite element modeling of Long Evans rats (n=85) subjected to graded, head-only proton (0, 10, 25, and 100 cGy, 150 MeV/n) and 28silicon (0, 10, 25, and 50 cGy, 300 MeV/n) radiation. The strength of the femoral neck will be examined due its clinical relevance to hip fractures. It has been shown in previous studies that bone mineral density was not reduced at the site of fracture. These findings question whether measurements of bone mineral density may be used to assess risk of hip fracture. The mechanisms leading to the irregular relationship between bone density and strength are still uncertain within literature and investigated to greater extent in clinical applications. Finite element analysis within this study simulated physiological loading of the femoral neck. No significant changes in femoral neck strength were found across doses of proton or 28silicon head-only radiation. Future work includes performing mechanical testing of the bone samples. Moving from mouse to larger animal models may also provide the increased lifespan for assessing the long-term outcomes of radiation exposure.

Dale Johnson↗

Pressure/temperature fluid cell apparatus for the neutron powder diffractometer instrument: Probing atomic structure in situ

This contribution describes a new local structure compatible gas/liquid cell apparatus for probing disordered materials at high pressures and variable temperatures in the Neutron Powder Diffraction instrument at the Lujan Neutron Scattering Center, Los Alamos National Laboratory. The new sample environment offers choices for sample canister thickness and canister material type. Finite element modeling is utilized to establish maximum allowable working pressures of 414 MPa at 15 K and 121 MPa at 600 K. High quality atomic pair distribution function data extraction and modeling have been demonstrated for a calibration standard (Si powder) and for supercritical and subcritical CO 2 measurements. As a result, the new sample environment was designed to specifically target experimental studies of the local atomic structures involved in geologic CO 2 sequestration, but will be equally applicable to a wide variety of energy applications, including sorption of fluids on nano/meso-porous solids, clathrate hydrate formation, catalysis, carbon capture, and H 2 and natural gas uptake/storage.

36 MATERIALS SCIENCE↗

Quantifying Grain Level Stress-Strain Behavior for AM40 via Instrumented Microindentation

ABSTRACT Microindentation is performed on hot isostatic pressed (HIP) Mg-Al (AM40) alloy samples produced by high-pressure die cast (HPDC) process for the purpose of quantifying the mechanical properties of the α-Mg grains. The process of obtaining elastic modulus and hardness from indentation load-depth curves is well established in the literature. A new inverse method is developed to extract plastic properties in this study. The method utilizes empirical yield strength-hardness relationship reported in the literature together with finite element modeling of the individual indentation. Due to the shallow depth of the indentation, indentation size effect (ISE) is taken into account when determining plastic properties. The stress versus strain behavior is determined for a series of indents. The resulting average values and standard deviations are obtained for future use as input distributions for microstructure-based property prediction of AM40.

AM40↗

Finite element approximations in transient analysis

Alternative techniques for finite element models for transient analyses are explored, using as an example a pinned-pinned beam. Attention is given to the relationships between desired accuracy and the available computer resources. Stiffness, mass, damping, and loading components approximations are examined with note taken of optimization methods. Generation of regular grid results indicates a linear relationship between the number of elements or degrees of freedom and the number of eigenvalues for a given accuracy. Comparisons of the finite element models in terms of accuracy/element or accuracy/degree of freedom were found to be inaccurate for determining model efficiency. Higher accuracies were available from lumped mass and damping models than assumptions of consistent loading, which was preferred to statically equivalent loading.

Melosh, R. J.↗

Residual stresses in Cu matrix composite surface deposits after laser melt injection

Abstract Tungsten carbide particles reinforced metal matrix composite (MMC) coatings can significantly improve surface wear resistance owing to their increased surface hardness. However, the presence of macro‐ and micro‐residual stresses in MMC coatings can have detrimental effects, such as reducing service life. In this study, neutron diffraction was used to determine the residual stresses in spherical fused tungsten carbide (sFTC) reinforced Cu matrix composite surface deposits after laser melt injection. We also developed a thermo‐mechanical coupled finite element model to predict residual stresses. Our findings reveal that sFTC/Cu composite deposits produced with a preheating temperature of 400°C have low residual stresses, with a maximum tensile residual stress of 98 MPa in the Cu matrix on the top surface. In contrast, the sFTC/bronze (CuAl10Ni5Fe4) composite deposit exhibits very high residual stresses, with a maximum tensile residual stress in the Cu matrix on the top surface reaching 651 MPa. These results provide a better understanding of the magnitudes and distributions of residual stresses in sFTC‐reinforced Cu matrix composite surface deposits manufactured via laser melt injection.

Zhang, Xingxing↗

Determination of crystal growth rates during rapid solidification of polycrystalline aluminum by nano-scale spatio-temporal resolution in situ transmission electron microscopy

In situ investigations of rapid solidification in polycrystalline Al thin films were conducted using nano-scale spatio-temporal resolution dynamic transmission electron microscopy. Differences in crystal growth rates and asymmetries in melt pool development were observed as the heat extraction geometry was varied by controlling the proximity of the laser-pulse irradiation and the associated induced melt pools to the edge of the transmission electron microscopy support grid, which acts as a large heat sink. Experimental parameters have been established to maximize the reproducibility of the material response to the laser-pulse-related heating and to ensure that observations of the dynamical behavior of the metal are free from artifacts, leading to accurate interpretations and quantifiable measurements with improved precision. Interface migration rate measurements revealed solidification velocities that increased consistently from ∼1.3 m s{sup −1} to ∼2.5 m s{sup −1} during the rapid solidification process of the Al thin films. Under the influence of an additional large heat sink, increased crystal growth rates as high as 3.3 m s{sup −1} have been measured. The in situ experiments also provided evidence for development of a partially melted, two-phase region prior to the onset of rapid solidification facilitated crystal growth. Using the experimental observations and associated measurements as benchmarks, finite-element modeling based calculations of the melt pool evolution after pulsed laser irradiation have been performed to obtain estimates of the temperature evolution in the thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Milestone 5 test report. Task 5, subtask 5.2: Tile to foam strength tests

This report summarizes work that has been performed to date on the strength of a cryotank insulation system using Rohacell foam and TUFI-coated AETB-12 ceramic tiles directly bonded to a simulated graphite-epoxy tank wall. Testing utilized a custom specimen design which consists of a long tensile specimen with eccentric loading to induce curvature similar to the curvature expected due to 'pillowing' of the tank when pressurized. A finite element model was constructed to predict the specific element strains in the test article, and to assist with design of the test specimen to meet the specific goals of curvature and laminate strain. The results indicate that the heat treated 3.25-pcf density Rohacell foam does not provide sufficient strength for the induced stresses due to curvature and stress concentration at the RTV bondline to the TUFI tile. The test was repeated using higher density non-heat treated Rohacell foam (6.9 pcf) without foam failure. The finite element model was shown to predict specimen behavior, and validation of the model was successful. It is pertinent to mention that the analyses described herein accurately predicted the failure of the heat treated foams and based on this analysis method it is expected that the untreated 3.25 pcf Rohacell foam will be successful.

Greenberg, H. S.↗

Behavior of Concrete Panels Reinforced with Synthetic Fibers, Mild Steel, and GFRP Composites Subjected to Blasts

The paper presents experimental data generated for calibrating finite element models to predict the performance of reinforced concrete panels with a wide range of construction details under blast loading. The specimens were 1.2 m square panels constructed using Normal Weight Concrete (NWC) or Fiber Reinforced Concrete (FRC). FRC consisted of macro-synthetic fibers dispersed in NWC. Five types of panels were tested: NWC panels with steel bars; FRC panels without additional reinforcement; FRC panels with steel bars; NWC panels with glass fiber reinforced polymer (GFRP) bars; and NWC panels reinforced with steel bars and external GFRP laminates on both faces. Each panel type was constructed with three thicknesses: 152 mm, 254 mm, and 356 mm. FRC panels with steel bars had the best performance for new construction. NWC panels reinforced with steel bars and external GFRP laminates on both faces had the best performance for strengthening or rehabilitation of existing structures. The performance of NWC panels with GFRP bars was strongly influenced by the bar spacing. The behavior of the panels is classified in terms of damage using immediate occupancy, life safety, and near collapse performance levels. Preliminary dynamic simulations are compared to the experimental results.

36 MATERIALS SCIENCE↗

Additive Manufacturing for Powering the Blue Economy Applications: A Tidal Turbine Blade Case Study: Preprint

As the marine renewable energy industry continues to expand, innovation in the manufacturing space must grow accordingly to reduce costs and ensure the economic feasibility of new technologies. Additive manufacturing, more commonly known as 3D printing, provides an alternative for rapid prototyping of marine hydrokinetic technologies, particularly supporting Powering the Blue Economy initiatives of the U.S. Department of Energy Water Power Technologies Office. This study explores the application of additive manufacturing in the development of marine hydrokinetic structures, focusing on material and printing method selection, design, and analysis of a 3D-printed spar for an axial-flow tidal turbine blade. Corrosion-resistant metals were deemed ideal due to the loads and harsh marine environment the blade would experience. Laser metal deposition methods were determined to be the most effective and scalable for the considered scale. The designed spar adapts its geometry to the blade - a feature uniquely suited to additive manufacturing - and is intended to serve as the blade's primary structural component. A finite element model was used to study stresses and deformations under loading conditions. The spar was manufactured using 316L stainless steel through direct energy deposition, and defects were assessed and recorded. Future efforts will include mechanical testing of the spar. This research establishes a benchmark process for using additive manufacturing in developing marine hydrokinetic structures, paving the way for future optimization and techno-economic analysis.

additive manufacturing↗

Acoustic wave propagation in heterogeneous structures including experimental validation

A finite element model was developed to solve for the acoustic pressure and energy fields in a heterogeneous suppressor. The derivations from the governing equations assumed that the material properties could vary with position resulting in a heterogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between different materials. For a two media region consisting of part air and part bulk absorber, a model was used to describe the bulk absorber properties in two directions. Complex metallic structures inside the air duct are simulated by simply changing element properties from air to the structural material in a pattern to describe the desired shapes. To verify the numerical theory, experiments were conducted without flow in a rectangular duct with a single folded cavity mounted above the duct and absorbing material mounted inside a cavity. Changes in a nearly plane wave sound field were measured on the wall opposite the absorbing cavity. Fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber, as a function of distance along the duct. The finite element model provides a convenient method for evaluating the acoustic properties of bulk absorbers.

Baumeister, Kenneth J.↗

Acoustic wave propagation in heterogeneous structures including experimental validation

A finite element model was developed to solve for the acoustic pressure and energy fields in a heterogeneous suppressor. The derivations from the governing equations assumed that the material properties could vary with position resulting in a heterogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between different materials. For a two-media region consisting of part air and part bulk absorber, a model was used to describe the bulk absorber properties in two directions. Complex metallic structures inside the air duct are simulated by simply changing element properties from air to the structural material in a pattern to describe the desired shapes. To verify the numerical theory, experiments were conducted without flow in a rectangular duct with a single folded cavity mounted above the duct and absorbing material mounted inside a cavity. Changes in a nearly plane wave sound field were measured on the wall opposite the absorbing cavity. Fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber, as a function of distance along the duct. The finite element model provides a convenient method for evaluating the acoustic properties of bulk absorbers.

Baumeister, Kenneth J.↗

A finite element overlay technique for modeling pinned composite joints

A finite element technique using overlaid plane stress elements is presented for modeling a pinned composite joint. The technique allows for discrete modeling of the pin and web regions of the joint in a two-dimensional (2-D) finite element analysis. Thus, predictive capability is substantially increased without the added complexity of a three-dimensional (3-D) analysis. The overlay technique requires a thorough evaluation of the model constraints between the joint components. Bearing panels are used to account for a reduced compressive modulus under bearing loads. Application of the plane stress overlay technique to the Filament Wound Case (FWC) Program is discussed. The experimental approach to determine bearing panel moduli using a double lap shear test is detailed, and finite element results are correlated with full-scale FWC joint test results.

Colvin, G. E.↗

Improved performance of HgCdTe infrared detector focal plane arrays by modulating light field based on photonic crystal structure

An HgCdTe long-wavelength infrared focal plane array photodetector is proposed by modulating light distributions based on the photonic crystal. It is shown that a promising prospect of improving performance is better light harvest and dark current limitation. To optimize the photon field distributions of the HgCdTe-based photonic crystal structure, a numerical method is built by combining the finite-element modeling and the finite-difference time-domain simulation. The optical and electrical characteristics of designed HgCdTe mid-wavelength and long-wavelength photon-trapping infrared detector focal plane arrays are obtained numerically. The results indicate that the photon crystal structure, which is entirely compatible with the large infrared focal plane arrays, can significantly reduce the dark current without degrading the quantum efficiency compared to the regular mesa or planar structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance assessment of 700-bar compressed hydrogen storage for light duty fuel cell vehicles

In this study, type 4 700-bar compressed hydrogen storage tanks were modeled using ABAQUS. The finite element model was first calibrated against data for 35-L subscale test tanks to obtain the composite translation efficiency, and then applied to full sized tanks. Two variations of the baseline T700/epoxy composite were considered in which the epoxy was replaced with a low cost vinyl ester resin and low cost resin with an alternate sizing. The results showed that the reduction in composite weight was attributed primarily to the lower density of the resin and higher fiber volume fraction in the composite due to increased squeeze-out with the lower viscosity vinyl ester resin. The system gravimetric and volumetric capacities for the onboard storage system that holds 5.6 kg H 2 are 4.2 wt% (1.40 kWh/kg) and 24.4 g-H 2 /L (0.81 kWh/L), respectively. The system capacities increase and carbon fiber requirement decreases if the in-tank amount of unrecoverable hydrogen is reduced by lowering the tank "empty" pressure. Models of an alternate tank design showed potential 4-7% saving in composite usage for tanks with a length-to-diameter (L/D) ratio of 2.8-3.0 but no saving for L/D of 1.7. Lastly, a boss with smaller opening and longer flange does not appear to reduce the amount of helical windings.

08 HYDROGEN↗

Error Localization Examples: Looking for a Needle in a Haystack

Finite element models (FEM) are routinely developed and used during fabrication of high dollar-value hardware. NASA, as part of the pre-flight certification of launch vehicles, routinely conducts vibration and static tests to calibrate models used for flight-risk assessments. During model calibration, certain areas of the model are modified, using engineering judgment and sensitivity analysis, to match the test results. Unfortunately, tools to identify problem areas in the FEM using test data directly are scarce and infrequently applied. Over the years, error localization algorithms have been proposed with very limited success. Recently, the Analytical Dynamics Model Improvement (ADMI) algorithm, which computes closed-form mass and stiffness corrections to match the test data exactly, have been shown to be an effective Error Localization Algorithm (ELA). The paper discusses three examples where ELA is used with simulated test data to locate problem areas. To gain confidence in the approach, the exact answer is shown along with ELA results. Results show that ELA is able to identify general problem areas consistent with known problem areas. In all examples, the ELA identified area is larger than the exact problem area. Nonetheless, with proper optimization tools, calibration results using the ELA identified areas provide excellent results.

model calibration↗

The Influence of PV Module Materials and Design on Solder Joint Thermal Fatigue Durability

Finite element model (FEM) simulations have been performed to elucidate the effect of flat plate photovoltaic (PV) module materials and design on PbSn eutectic solder joint thermal fatigue durability. The statistical method of Latin Hypercube sampling was employed to investigate the sensitivity of simulated damage to each input variable. Variables of laminate material properties and their thicknesses were investigated. Using analysis of variance, we determined that the rate of solder fatigue was most sensitive to solder layer thickness, with copper ribbon and silicon thickness being the next two most sensitive variables. By simulating both accelerated thermal cycles (ATCs) and PV cell temperature histories through two characteristic days of service, we determined that the acceleration factor between the ATC and outdoor service was independent of the variables sampled in this study. This result implies that an ATC test will represent a similar time of outdoor exposure for a wide range of module designs. This is an encouraging result for the standard ATC that must be universally applied across all modules.

14 SOLAR ENERGY↗

A full-ring variable-aperture cadmium zinc telluride system for whole-body single photon emission computed tomography: realistic simulations with phantoms

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. Purpose We have designed a cadmium zinc telluride (CZT) variable-aperture full-ring SPECT scanner instrumented with a broad-energy tungsten collimator intended to fill this technological gap. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT scanner breaks the paradigm of the standard dual- and triple-head rotational SPECT systems, utilizing a larger detector area in each scan increasing the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized 99mTc phantoms. Methods We developed a realistic model of our detector by using a combination of a Geant4 Monte Carlo simulation and a CZT detector response model based on a finite element model. Our approach models the characteristic low-energy tail effect in CZT that noticeably affects the sensitivity and the quality of the scatter correction in CZT detectors. We implement a modified dual energy window scatter correction adapted to include the CZT low-energy tail effect. A dedicated correction is also developed to eliminate the undesirable truncation observed in images given the presence of detector edges and gaps between detectors, due to the non-rotational nature of our device. Corrections for the attenuation, detector response and the presence of collimators are also included. The images are reconstructed using the maximum-likelihood expectation-maximization algorithm implemented in the reconstruction open software STIR. Detector and reconstruction performance are characterized with a Derenzo phantom and a body-sized National Electrical Manufacturers Association (NEMA) Image Quality (IQ) phantom containing 99mTc. Results Our SPECT design can resolve 6.4mm rods in a Derenzo phantom and obtain a good image contrast with the IQ phantom. Explicit testing of the gap and edge correction is provided, showing an excellent performance in eliminating the image truncation artifacts. Our modified scatter correction shows no overestimation of the contrast-recovery ratio for our realistic CZT detector model, as opposed to the cases without correction and with a standard dual-energy window scatter correction. Conclusions In this paper, we further demonstrate the performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99mTc imaging for brain perfusion and 123I imaging for dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design performs similarly to the VERITON SPECT from Spectrum Dynamics, a commercial full-ring CZT SPECT camera, with the potential advantage of the broader energy range of application given by our custom-design tungsten collimators. Here, our device combines high sensitivity and image resolution with a broad-energy imaging application for the purpose of clinical imaging and theranostics of emerging radionuclides.

62 RADIOLOGY AND NUCLEAR MEDICINE↗