Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element simulation”

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

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

At least 901 records · Page 50

Calibration of Airframe and Occupant Models for Two Full-Scale Rotorcraft Crash Tests

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. Accelerations and kinematic data collected from the crash tests were compared to a system integrated finite element model of the test article. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the second full-scale crash test. This combination of heuristic and quantitative methods was used to identify modeling deficiencies, evaluate parameter importance, and propose required model changes. It is shown that the multi-dimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and co-pilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. This approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and test planning guidance. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, thereby reducing overall development costs.

Annett, Martin S.↗

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↗

Design, fabrication, simulation, and testing of additively manufactured lattice-based copper heat sinks

Here, this study investigates the potential advantages of lattice structure-based bound metal material extrusion (MEX) 3D printing for fabricating high-performance copper heat sinks. Copper powder-filled polymer filaments, with a copper content of > 90 wt.%, were developed specifically for the bound metal MEX 3D printing process. Three types of structures—planar, strut, and surface lattices—were 3D printed to facilitate efficient heat transfer pathways within the heat sinks. Subsequent post-processing steps, including polymer removal and sintering, were performed to achieve dense copper parts. Hot isostatic pressing was further employed to enhance the sintered density from 93 to 98%. Finite element analysis (FEA) simulations were conducted to assess the heat transfer efficiency of the designs, and heat transfer experiments were performed using a custom setup to validate the simulation results. Additionally, this research explores the use of extended hold times during pre-sintering and a reduced atmosphere to enhance the %IACS values (electrical conductivity) and thermal performance of the bound metal MEX 3D printed copper heat sinks. The investigation combines experimental analysis, including simulations and heat transfer experiments, to gain insights into the structure-material property relationships and optimize the thermal performance of the printed heat sinks.

Ajjarapu, Kameswara Pavan Kumar [Oak Ridge Nationa↗

The Shifted Interface Method: A flexible approach to embedded interface computations

We propose a new embedded finite element method to simulate partial differential equations over domains with internal interfaces. Our approach belongs to the family of surrogate/approximate interface methods and relies on the idea of shifting the location and value of jump interface conditions. This choice has the goal of preserving optimal convergence rates while avoiding small cut cells and related problematic issues, typical of traditional embedded methods. The proposed approach is accurate, robust, efficient, and simple to implement. We apply this concept to internal interface computations in the context of the mixed Poisson problem, also known as the Darcy flow problem. An extensive set of numerical tests is provided to demonstrate the performance of the proposed approach.

42 ENGINEERING↗

On residual stress analysis and microstructural evolution for stainless steel type 304 spent nuclear fuel canisters weld joint: Numerical and experimental studies

Chloride-induced stress corrosion cracking (CISCC) in the weldments of austenitic stainless steel canisters is one of the primary safety concerns during the dry storage of used nuclear fuel at Independent Spent Fuel Storage Installations in coastal areas. In order to evaluate the CISCC behavior in the canister, a 3D sequentially coupled thermo-mechanical finite element model was built to simulate the residual stresses induced by two intersecting longitudinal and circumferential multi-pass welds in austenitic stainless steel mockup canister. Weld-induced residual stresses from simulation were in good agreement with the experimental measurement results by deep-hole drilling and contour methods. Through-wall high tensile axial and hoop stresses were observed in both the fusion zone and heat affected zone (HAZ) in the longitudinal and circumferential weld, respectively. On the other hand, the microstructure within the weldments of 304 stainless steel mockup was characterized using optical microscopy, scanning electron microscopy, electron back-scattered diffraction and energy dispersive X-ray spectroscopy methods. A few amount of MnS inclusions in both the base metal and HAZ, and a lot of (Mn,Si) x O y inclusions in the fusion zone were found. Meanwhile, there is some α'-martensite in the base metal, while the martensite was reversed in the HAZ. In the fusion zone, there is a lot of δ-ferrite. The martensite reversion phenomenon was also verified by the combination of Feritscope test and thermal simulation, which could predict a martensite reversion temperature as 400C°. The welding residual stresses and microstructure analysis results could offer the instructive information for the prediction of CISCC behavior of the spent nuclear fuel canisters.

36 MATERIALS SCIENCE↗

Thermocouple-tip-exposing temperature assessment technique for evaluating photothermal conversion efficiency of plasmonic nanoparticles at low laser power density

A new thermocouple (TC) tip-exposing temperature assessment technique that combines experimental temperature measurements with a numerical model of the photothermal conversion efficiency η is presented. The proposed technique is designed to evaluate η for a gold-coated superparamagnetic iron oxide nanoparticle (SPIO-Au NP) solution (26 nm, 12–70 ppm) at low continuous wave laser power (103 mW, 532 nm) irradiation in a convenient manner under ambient conditions. The TC tip temperature is measured during the first 30 s of the laser exposure, and the results are combined with a finite element model to simulate the temperature rise of the NP solution for a given concentration. The value of η is adjusted in the model until the model agrees with the measured transient TC temperature rise. Values of η = 1.00 were observed for all concentrations. Theoretical predictions of η derived by Mie theory confirmed the near unity conversion efficiency of the as-synthesized SPIO-Au NPs. Advantages of the current technique include co-locating the TC tip in the geometric center of the laser-heated region, rather than outside of this region. In addition, the technique can be done under ambient room conditions using unmodified commercially available hardware.

Numerical methods↗

Electromagnetic Modeling of Wind tunnel Magnetic Suspension and Balance Systems

An analytical framework for the open-loop behavior of a permanent magnet element levitated within an applied magnetic field are shown, in the context of application to wind tunnel Magnetic Suspension and Balance Systems (MSBS). Various modes of motion are identified, which correlate with observed behavior of the NASA/ODU 6-inch MSBS. This system is being developed as a technology demonstrator for dynamic stability testing of atmospheric entry capsules, with possible application to a supersonic wind tunnel. The analytic equations are populated with field information from experimental measurements and a finite element model. Finally, a system simulation can provide alternative estimates of system dynamic behavior via model linearization, with preliminary comparisons presented. Taken together, the analytical framework, experimental measurement, finite element analysis, and dynamic simulation provide a complete understanding of the characteristics and behavior of the MSBS.

magnetic suspension↗

A continuum dislocation-based model of wedge microindentation of single crystals

Recent Electron Backscatter Diffraction (EBSD) experiments have revealed the emergence of heterogeneous dislocation microstructures forming under a wedge indenter in fcc crystals, where micro-meter dislocation patterns challenge the predictions of traditional models of plasticity. In order to explain the formation of these features and develop a relationship between the force-displacement curve and the dislocation substructure, here we present a model of wedge indentation based on the continuum theory of dislocations. The model accounts for large deformation kinematics through the multiplicative split of the deformation gradient tensor, where the incompatible plastic component of deformation results from the flux of dislocations on different and interacting slips systems. Constitutive equations for dislocation fluxes are determined from a dissipative variational principle. As a result, each dislocation density satisfies an initial-boundary value problem with convective-diffusive character, which is coupled to the macroscopic stress and displacement fields governing the deformation process. Solution to the self-consistent continuum formulation is found using the finite element method. Computer simulations mimic the experimental conditions of wedge micro-indentation experiments into Ni single-crystals used by Kysar et al. (2010a). A comparison of overall dislocation density distribution and macroscopic mechanical response shows good overall agreement with the experimental results in terms of the detailed features of dislocation patterns and lattice rotations as well as the macroscopic force-displacement response.

36 MATERIALS SCIENCE↗

Dynamic behavior of chemically tunable mechano-responsive hydrogels

Here, using theory and simulation, we model the mechanical behavior of gels that encompass loops and dangling chain ends. If the loops remain folded and dangling ends are chemically inert, then these topological features just serve as defects. If, however, the loops unfold to expose the hidden (“cryptic”) binding sites and the ends of the dangling chains are reactive, these moieties can form bonds that improve the gel's mechanical properties. For gels with a lower critical solubility temperature (LCST), we systematically switch on the possible unfolding and binding events. To quantify the resulting effects, we derive equations for the gel's equilibrium and dynamic elastic moduli. We also use a finite element approach to simulate the gel's response to deformation and validate the analytic calculations. Herein, we show that the equilibrium moduli are highly sensitive to the presence of unfolding and binding transitions. The dynamical moduli are sensitive not only to these structural changes, but also to the frequency of deformation. For example, when reactive ends bind to exposed cryptic sites at T = 29 °C and relatively high frequency, the storage shear modulus is 119% greater than the corresponding equilibrium value, while the storage Young's modulus is 109% greater than at equilibrium. These findings provide guidelines for tuning the chemical reactivity of loops and dangling ends and the frequency of deformation to tailor the mechano-responsive behavior of polymer networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated simulation of vertical fracture propagation induced by water injection and its borehole electromagnetic responses in shale gas systems

We analyze fracture propagation induced by hydraulic fracturing with water injection and examine their detectability with crosswell electromagnetic (EM) geophysical methods. For rigorous 3D coupled flow-geomechanical modeling, we employ a numerical method that can model failure by tensile and shear stresses, dynamic nonlinear permeability, dual continuum approach, and thermo-poro-mechanical effects. From numerical simulation, we find that the fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics whereas water saturation is determined by multiphase flow. At early times, the water front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with the injected water. However, at late times, movement of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, filled with reservoir gas. During fracture propagation, the coupled flow-geomechanical models are transformed via a rock-physics model into electrical conductivity models. We employ a full 3D finite-element EM geophysical simulator to evaluate the sensitivity of the crosswell EM method to fracture propagation. It is shown that anomalous distribution of electrical conductivity is closely related to the injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to electrical conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for monitoring the injected fluids during hydraulic fracturing operations.

58 GEOSCIENCES↗

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗

Identification of the dynamics of a two-dimensional grid structure using least square lattice filters

The basic theory of least square lattice filters and their use in identification of structural dynamics systems is summarized. Thereafter, this theory is applied to a two-dimensional grid structure made of overlapping bars. Previously, this theory has been applied to an integral beam. System identification results are presented for both simulated and experimental tests and they are compared with those predicted using finite element modelling. The lattice filtering approach works well for simulated data based on finite element modelling. However, considerable discrepancy exists between estimates obtained from experimental data and the finite element analysis. It is believed that this discrepancy is the result of inadequacies in the finite element modelling to represent the damped motion of the laboratory apparatus.

Montgomery, R. C.↗

3D Finite Element Analysis of Particle-Reinforced Aluminum

Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.

Shen, H.↗

Stability Analysis of Intertank Formed Skin/Stringer Compression Panel with Simulated Damage

The External Tank (ET) is a component of the Space Shuttle launch vehicle that contains fuel and oxidizer. During launch, the ET supplies the space shuttle main engines with liquid hydrogen and liquid oxygen. In addition to supplying fuel and oxidizer, it is the backbone structural component of the Space Shuttle. It is comprised of a liquid hydrogen (LH2) tank and a liquid oxygen (LOX) tank, which are separated by an Intertank. The Intertank is a stringer-stiffened cylindrical structure with hat-section stringers that are roll formed from aluminum-lithium alloy Al-2090. Cracks in the Intertank stringers of the STS-133 ET were noticed after a November 5, 2010 launch attempt. The cracks were approximately nine inches long and occurred on the forward end of the Intertank (near the LOX tank), along the fastener line, and were believed to have occurred while loading the ET with the cryogenic propellants. These cracks generated questions about the structural integrity of the Intertank. In order to determine the structural capability of the Intertank with varying degrees of damage, a finite element model (FEM) simulating a 1995 compression panel test was analyzed and correlated to test data. Varying degrees of damage were simulated in the FEM, and non-linear stability analyses were performed. The high degree of similarity between the compression panel and the Intertank provided confidence that the ET Intertank would have similar capabilities.

Harper, David W.↗

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlocal Analysis of Natural Vibrations of Carbon Nanotubes

When the atomic level is considered, the dynamic properties have a major influence on nanostructures behavior because of their ultrahigh or very high natural frequencies. Vibrations are essential in the analysis of resonators, oscillators, and sensors for carbon nanotube-based devices. The present paper is dedicated to the eigenfrequencies analysis of single-walled carbon nanotubes. In the analysis, the small-scale coefficient is introduced, and the nonlocal elasticity theory is applied in the modal analysis of the carbon nanotubes. The correlation between nonlocal small-scale parameter and the vibrational behavior of the carbon nanotubes is studied. The determination of nonlocal parameter is based on the structural finite element model of the nanotube. A detailed parametric study is realized to investigate the effects of the length-to-diameter ratio L/D, nonlocal coefficient, the influence of boundary conditions on eigenfrequencies of carbon nanotubes. The main contribution of this work is the evaluation of the nonlocal parameter for single-walled carbon nanotubes based on the finite element approximation in modeling the dynamic behavior instead of commonly used MD simulations. It was clearly displayed that the application of the nonlocal beam models involves analysis of the small-scale parameter, whose value determines the final eigenfrequencies. The highest values of eigenfrequencies were observed for the local beam model. It was presented that the eigenfrequencies for carbon nanotubes decrease when the nanotube length increases for all analyzed modes. On the other hand, when the diameter of carbon nanotubes increases, the eigenfrequency values rise. The numerical validation of the nonlocal parameter revealed the high dependence on the boundary conditions and the mode number. However, the influence of the mode number and the boundary conditions is vanishing for long single-walled carbon nanotubes.

36 MATERIALS SCIENCE↗

Effect of annealing on the laser induced damage of polished and CO{sub 2} laser-processed fused silica surfaces

We investigate the effect of different heat treatments on the laser-induced damage probabilities of fused silica samples. Isothermal annealing in a furnace is applied, with different temperatures in the range 700–1100 °C and 12 h annealing time, to super-polished fused silica samples. The surface flatness and laser damage probabilities at 3 ns, 351 nm are measured before and after the different annealing procedures. We have found a significant improvement of the initial laser damage probabilities of the silica surface after annealing at 1050 °C for 12 h. A similar study has been conducted on CO{sub 2} laser-processed sites on the surface of the samples. Before and after annealing, we have studied the morphology of the sites, the evolution of residual stress, and the laser-induced damage threshold measured at 351 nm, 3 ns. In this case, we observe that the laser damage resistance of the laser created craters can reach the damage level of the bare fused silica surface after the annealing process, with a complete stress relieve. The obtained results are then compared to the case of local annealing process by CO{sub 2} laser irradiation during 1 s, and we found similar improvements in both cases. The different results obtained in the study are compared to numerical simulations made with a thermo-mechanical model based on finite-element method that allows the simulation of the isothermal or the local annealing process, the evolution of stress and fictive temperature. The simulation results were found to be very consistent with experimental observations for the stresses evolution after annealing and estimation of the heat affected area during laser-processing based on the density dependence with fictive temperature. Following this work, the temperature for local annealing should reach 1330–1470 °C for an optimized reduction of damage probability and be below the threshold for material removal, whereas furnace annealing should be kept below the annealing point to avoid sample deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗