Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element Modeling”

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,603 records · Page 89

Simulating the effect of grain structure and porosity on creep for powder bed fusion 316H

This report details results from modeling studies focused on identifying the potential factors influencing creep performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF), with a focus on experimentally-observed and predicted differences between the response of the AM material and conventionally manufactured wrought 316H. The studies presented here systematically look at the role of grain structure, porosity and texture in the anisotropic creep behavior of the material under conditions expected in high temperature advanced nuclear reactors. A physics-based Crystal Plasticity Finite Element modeling approach has allowed us to look at these factors in isolation and study their role in the material's long-term deformation behavior. An update to the precipitation model used in the constitutive framework for 316H is also presented here, with the aim of improving the accuracy of our results from this modeling effort. This work is a step in the direction of gaining mechanistic understanding of the long-term deformation behavior of LPBF 316H, helping us model its long term performance and reducing the qualification time for the material.

36 MATERIALS SCIENCE↗

First in-core simultaneous measurements of nuclear heating and thermal neutron flux obtained with the innovative mobile calorimeter CALMOS inside the OSIRIS reactor

Nuclear heating inside a MTR reactor has to be known in order to design and to run irradiation experiments which have to fulfill target temperature constraints. This measurement is usually carried out by calorimetry. The innovative calorimetric system, CALMOS, has been studied and built in 2011 for the 70 MWth OSIRIS reactor operated by CEA. Thanks to a new type of calorimetric probe, associated to a specific displacement system, it provides measurements along the fissile height and above the core. This development required preliminary modelling and irradiation of mock-ups of the calorimetric probe in the ex-core area, where nuclear heating rate does not exceed 2 W.g{sup -1}. The calorimeter working modes, the different measurement procedures allowed with such a new probe, the main modeling and experimental results and expected advantages of this new technique have been already presented. However, these first in-core measurements were not performed beyond 6 W.g{sup -1}, due to an inside temperature limitation imposed by a safety authority requirement. In this paper, we present the first in-core simultaneous measurements of nuclear heating and conventional thermal neutron flux obtained by the CALMOS device at the 70 MW nominal reactor power. For the first time, this experimental system was operated in nominal in-core conditions, with nominal neutron flux up to 2.7 10{sup 14} n.cm{sup -2}.s{sup -1} and nuclear heating up to 12 W.g{sup -1}. A comprehensive measurement campaign carried out from 2013 to 2015 inside all accessible irradiation locations of the core, allowed to qualify definitively this new device, not only in terms of measurement ability but also in terms of reliability. After a brief reminder of the calorimetric cell configuration and displacement system specificities, first nuclear heating distributions at nominal power are presented and discussed. In order to reinforce the heating evaluation, a systematic comparison is made between results obtained by different methods, the probe calibration coefficient and the zero method. Thermal neutron flux evaluation from the SPND signal processing required a specific TRIPOLI-4 Monte Carlo calculation which has been performed with the precise CALMOS cell geometry. In addition, the Finite Element model for temperatures map prediction inside the calorimetric cell has been upgraded with the recent experimental data obtained up to 12 W.g{sup -1}. The Kc coefficient, taking into account nonlinearities with regard to the calibration, has been reevaluated so as to make relevant measurements up to the nominal reactor power. Finally, the experience feedback acquired until now with this first CALMOS version led us to improvement perspectives. A second device is currently under manufacturing and main technical options chosen for this second version are presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Finite-element analysis of transverse compressive and thermal loads on Nb 3 Sn wires with voids

High-field superconducting magnets play a very important role in many large-scale physics experiments, particularly particle colliders and fusion confinement devices such as Large Hadron Collider (LHC) and International Thermonuclear Experimental Reactor (ITER). The two most common superconductors used in these applications are NbTi and Nb 3 Sn. Nb 3 Sn wires are favored because of their significantly higher J c (critical current density) for higher field applications. The main disadvantage of Nb 3 Sn is that the superconducting performance of the wire is highly strain sensitive and it is very brittle. This strain sensitivity is strongly influenced by two factors: plasticity and cracked filaments. Cracks are induced by large stress concentrators that can be traced to the presence of voids in the wire. We develop detailed 2-D and 3-D finite-element models containing wire filaments and different possible distributions of voids in a bronze-route Nb 3 Sn wire. We apply compressive transverse loads for various cases of void distributions to simulate the stress and strain response of a Nb 3 Sn wire under the Lorentz force. Furthermore, this paper improves our understanding of the effect voids have on the Nb 3 Sn wire's mechanical properties, and in so, the connection between the distribution of voids and performance degradation such as the correlation between irreversible strain limit and the void-induced local stress concentrations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strain and lattice orientation distribution in SiN/Ge complementary metal–oxide–semiconductor compatible light emitting microstructures by quick x-ray nano-diffraction microscopy

This paper presents a study of the spatial distribution of strain and lattice orientation in CMOS-fabricated strained Ge microstripes using high resolution x-ray micro-diffraction. The recently developed model-free characterization tool, based on a quick scanning x-ray diffraction microscopy technique can image strain down to levels of 10{sup −5} (Δa/a) with a spatial resolution of ∼0.5 μm. Strain and lattice tilt are extracted using the strain and orientation calculation software package X-SOCS. The obtained results are compared with the biaxial strain distribution obtained by lattice parameter-sensitive μ-Raman and μ-photoluminescence measurements. The experimental data are interpreted with the help of finite element modeling of the strain relaxation dynamics in the investigated structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Advances in reduction techniques for tire contact problems

Some recent developments in reduction techniques, as applied to predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities, are reviewed. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of reduction techniques, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface. Also, the research topics which have high potential for enhancing the effectiveness of reduction techniques are outlined.

Noor, Ahmed K.↗

Reduction technique for tire contact problems

A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Noor, Ahmed K.↗

Reduced-Order Blade Mistuning Analysis Techniques Developed for the Robust Design of Engine Rotors

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo-Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using eigenfrequency curve veerings to identify "danger zones" in the operating conditions--ranges of rotational speeds and engine orders in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued. Several methods will be investigated, including the use of intentional mistuning patterns to mitigate the harmful effects of random mistuning, and the modification of disk stiffness to avoid reaching critical values of interblade coupling in the desired operating range. Recent research progress is summarized in the following paragraphs. First, significant progress was made in the development of the component mode mistuning (CMM) and static mode compensation (SMC) methods for reduced-order modeling of mistuned bladed disks (see the following figure). The CMM method has been formalized and extended to allow a general treatment of mistuning. In addition, CMM allows individual mode mistuning, which accounts for the realistic effects of local variations in blade properties that lead to different mistuning values for different mode types (e.g., mistuning of the first torsion mode versus the second flexural mode). The accuracy and efficiency of the CMM method and the corresponding Turbo-Reduce code were validated for an example finite element model of a bladed disk.

Min, James B.↗

MSC/NASTRAN Stress Analysis of Complete Models Subjected to Random and Quasi-Static Loads

Space payloads, such as those which fly on the Space Shuttle in Spacelab, are designed to withstand dynamic loads which consist of combined acoustic random loads and quasi-static acceleration loads. Methods for computing the payload stresses due to these loads are well known and appear in texts and NASA documents, but typically involve approximations such as the Miles' equation, as well as possible adjustments based on "modal participation factors." Alternatively, an existing capability in MSC/NASTRAN may be used to output exact root mean square [rms] stresses due to the random loads for any specified elements in the Finite Element Model. However, it is time consuming to use this methodology to obtain the rms stresses for the complete structural model and then combine them with the quasi-static loading induced stresses. Special processing was developed as described here to perform the stress analysis of all elements in the model using existing MSC/NASTRAN and MSC/PATRAN and UNIX utilities. Fail-safe and buckling analyses applications are also described.

Hampton, Roy W.↗

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗

Vertical Drop Testing and Analysis of the WASP Helicopter Skid Gear

Human occupant modeling and injury risk assessment have been identified as areas of research for improved prediction of rotorcraft crashworthiness within the NASA Aeronautics Program's Subsonic Rotary Wing Project. As part of this effort, an experimental program was conducted to assess the impact performance of a skid gear for use on the WASP kit-built helicopter, which is marketed by HeloWerks, Inc. of Hampton, Virginia. Test data from a drop test at an impact velocity of 8.4 feet-per-second were used to assess a finite element model of the skid gear test article. This assessment included human occupant analytic models developed for execution in LS-DYNA. The test article consisted of an aluminum skid gear mounted beneath a steel plate. A seating platform was attached to the upper surface of the steel plate, and two 95th percentile Hybrid III male Aerospace Anthropomorphic Test Devices (ATDs) were seated on the platform and secured using a four-point restraint system. The goal of the test-analysis correlation is to further the understanding of LS-DYNA ATD occupant models and responses in the vertical (or spinal) direction. By correlating human occupant experimental test data for a purely vertical impact with the LS-DYNA occupant responses, improved confidence in the use of these tools and better understanding of the limitations of the automotive-based occupant models for aerospace application can begin to be developed.

Fuchs, Yvonne T.↗

Vertical Drop Testing and Analysis of the WASP Helicopter Skid Gear

Human occupant modeling and injury risk assessment have been identified as areas of research for improved prediction of rotorcraft crashworthiness within the NASA Aeronautics Program s Subsonic Rotary Wing Project. As part of this effort, an experimental program was conducted to assess the impact performance of a skid gear for use on the WASP kit-built helicopter, which is marketed by HeloWerks, Inc. of Hampton, Virginia. Test data from a drop test at an impact velocity of 8.4 feet-per-second were used to assess a finite element model of the skid gear test article. This assessment included human occupant analytic models developed for execution in LS-DYNA. The test article consisted of an aluminum skid gear mounted beneath a steel plate. A seating platform was attached to the upper surface of the steel plate, and two 95th percentile Hybrid III male Aerospace Anthropomorphic Test Devices (ATDs) were seated on the platform and secured using a four-point restraint system. The goal of the test-analysis correlation is to further the understanding of LS-DYNA ATD occupant models and responses in the vertical (or spinal) direction. By correlating human occupant experimental test data for a purely vertical impact with the LS-DYNA occupant responses, improved confidence in the use of these tools and better understanding of the limitations of the automotive-based occupant models for aerospace application can begin to be developed.

Fuchs, Yvonne T.↗

Characterization and Modeling of Cell Wall Imperfections in Aluminum Honeycomb Cores using X-ray CT Imaging

Honeycomb cell imperfections may significantly affect the failure of sandwich composites in core crushing and shear modes. Identification and quantification of the different types of imperfections and evaluation of their effects on the mechanical response of honeycomb cores in flatwise compression is presented. X-Ray computed tomography (CT) scans of honeycomb cores in co-cured composite sandwich panels were performed. Image analysis methods are applied to characterize and quantify core cell geometry and cell wall waviness. An approach is developed for extracting cell wall mid surface planes and waviness for constructing finite element models using shell elements. Analysis results from models with the as-manufactured imperfections are compared to results from models with idealized geometries to evaluate the effect of the as-manufactured imperfections on the initiation of cell wall buckling and subsequent core crushing.

aluminum honeycomb↗

Colocalized Raman spectroscopy – scanning electrochemical microscopy investigation of redox flow battery dialkoxybenzene redoxmer degradation pathways

Non-aqueous redox flow batteries offer high voltages for grid-level energy storage technologies. However, decomposition of the redoxmers - redox-active molecules that make up the anolyte and catholyte in the negative and positive cell compartments - is an important challenge to overcome for long-term storage. Here, we present a spectroelectrochemical study of the catholyte candidate 2,3-dimethyl-1,4-dialkoxybenzene (C7) and its decomposition mechanisms in the presence of a model nucleophilic base, pyridine. We utilize colocalized Raman microscopy and scanning electrochemical microscopy (Raman-SECM) to quantify the chemical rates of decom-position and qualitatively identify the reaction intermediates and products. A detailed study on how the Raman-SECM parameters (electrode distance, substrate electrode, and laser focal height) influence the Raman signal of the charged catholyte C7 ∙+ is presented. Using optimized conditions, we monitored the deprotonation of C7 ∙+ via Raman spectroscopy and the subsequent hydrogen abstraction from solvent molecules to regenerate C7 via electrochemistry. Finite element modeling was used to fit electrochemical and spectroscopic data, quantifying the deprotonation rates as k dep = 2000 and 700 L mol -1 s -1 and abstraction rates as k abs = 0.5 and 0.2 s - 1 in propylene carbonate and acetonitrile solvents, respectively. Our results show the value of spatiotemporal reso-lution in evaluating the chemical and electrochemical behavior of materials for redox flow batteries.

25 ENERGY STORAGE↗

The New NASA Orbital Debris Engineering Model ORDEM2000

The NASA Orbital Debris Program Office at Johnson Space Center has developed a new computer-based orbital debris engineering model, ORDEM2000, which describes the orbital debris environment in the low Earth orbit region between 200 and 2000 km altitude. The model is appropriate for those engineering solutions requiring knowledge and estimates of the orbital debris environment (debris spatial density, flux, etc.). ORDEM2000 can also be used as a benchmark for ground-based debris measurements and observations. We incorporated a large set of observational data, covering the object size range from 10 mm to 10 m, into the ORDEM2000 debris database, utilizing a maximum likelihood estimator to convert observations into debris population probability distribution functions. These functions then form the basis of debris populations. We developed a finite element model to process the debris populations to form the debris environment. A more capable input and output structure and a user-friendly graphical user interface are also implemented in the model. ORDEM2000 has been subjected to a significant verification and validation effort. This document describes ORDEM2000, which supersedes the previous model, ORDEM96. The availability of new sensor and in situ data, as well as new analytical techniques, has enabled the construction of this new model. Section 1 describes the general requirements and scope of an engineering model. Data analyses and the theoretical formulation of the model are described in Sections 2 and 3. Section 4 describes the verification and validation effort and the sensitivity and uncertainty analyses. Finally, Section 5 describes the graphical user interface, software installation, and test cases for the user.

Liou, Jer-Chyi↗

Mechanical Performance of the First Two Prototype 4.5 m Long Nb$_3$Sn Low-$β$ Quadrupole Magnets for the Hi-Lumi LHC Upgrade

The U.S. High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) team is collaborating with CERN in the design and fabrication of the first 4.5 m long MQXFA magnets, a 150 mm aperture high-field Nb$_3$Sn quadrupole magnet that uses the aluminum shell-based bladder-and-key technology. The first two prototype magnets, MQXFAP1 and MQXFAP2, were assem-bled and tested while the first pre-series structure (MQXFA03) was in fabrication. This paper summarizes the mechanical perfor-mance of these prototype structures based on the comparison of measured strain gauge data with finite element model analyses from all load steps to powering. The MQXFAP1 magnet almost reached ultimate current before a short to ground was detected and the test was stopped. The MQXFAP2 magnet experienced a low training performance due to a fractured aluminum shell. MQXFAP1b was rebuilt with a new replacement coil, but an old coil limited the magnet from achieving the ultimate current. The mitigations and analyses of these prototype magnets are discussed in the context of the transition to pre-series production.

43 PARTICLE ACCELERATORS↗

Using Numerical Modeling to Simulate Space Capsule Ground Landings

Experimental work is being conducted at the National Aeronautics and Space Administration s (NASA) Langley Research Center (LaRC) to investigate ground landing capabilities of the Orion crew exploration vehicle (CEV). The Orion capsule is NASA s replacement for the Space Shuttle. The Orion capsule will service the International Space Station and be used for future space missions to the Moon and to Mars. To evaluate the feasibility of Orion ground landings, a series of capsule impact tests are being performed at the NASA Langley Landing and Impact Research Facility (LandIR). The experimental results derived at LandIR provide means to validate and calibrate nonlinear dynamic finite element models, which are also being developed during this study. Because of the high cost and time involvement intrinsic to full-scale testing, numerical simulations are favored over experimental work. Subsequent to a numerical model validated by actual test responses, impact simulations will be conducted to study multiple impact scenarios not practical to test. Twenty-one swing tests using the LandIR gantry were conducted during the June 07 through October 07 time period to evaluate the Orion s impact response. Results for two capsule initial pitch angles, 0deg and -15deg , along with their computer simulations using LS-DYNA are presented in this article. A soil-vehicle friction coefficient of 0.45 was determined by comparing the test stopping distance with computer simulations. In addition, soil modeling accuracy is presented by comparing vertical penetrometer impact tests with computer simulations for the soil model used during the swing tests.

Heymsfield, Ernie↗

Surface orientation ambiguity for single molecules at dielectric interfaces

Fluorescent molecules emit light in a dipole radiation pattern that can be used to infer their orientation through defocused fluorescence microscopy. Proper measurement of the orientation requires mathematical modeling of the radiation pattern expected for a dipole in the geometry of interest and subsequent comparison against experimental data. We point out an ambiguity in common calculations of these patterns that appears to compromise orientation measurements for molecules that are especially near dielectric surfaces. This results in a rotation of the measured emission dipole toward the surface for near-interface molecules, which can be mistaken for a preferentially horizontal orientation among the emitters. The proper treatment for on-surface emitters requires consideration of finite-sized current elements between two dielectric media, and we show that the theoretical ambiguity can be lifted via finite-element modeling. A prescription is provided for correcting measured orientations at arbitrary interfaces.

Dey, E. [University of Texas, Arlington, TX (Unite↗

Semiempirical Model Would Control Czochralski Process

Semiempirical mathematical model proposed for control of growth of single crystals of silicon by Czochralski process. Expresses dependence of pulling rate and shape of liquid/solid interface upon important process variables; radius of growing crystal, temperature of crucible, level of melt, and height of exposed portion of crucible wall. Necessary to control shape of interface in manner consistent with other variables, to maintain radially uniform concentration of dopant, and reduce thermally induced stresses in vicinity of interface. Used to simulate complete growth cycles without requiring excessive computer time consumed by rigorous finite-element modeling.

Dudukovic, M. P.↗