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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Effect of Tempering and Rolling on Fatigue Crack Growth Behavior of Modified 9Cr-1Mo Steel

In the present work, fatigue crack growth (FCG) behavior of modified 9Cr-1Mo steel has been investigated near fatigue threshold and Paris regime at room temperature (RT) and 650 °C. The modified 9Cr-1Mo steel was subjected to normalizing at 1050 °C followed by tempering (650-750 °C) and warm rolling at 550 °C followed by tempering in order to improve its precipitation state and refine microstructure. The microstructural parameters such as grain size, block size, lath width and precipitate size were examined and correlated with FCG and tensile properties of investigated samples. The results showed that decreasing tempering temperature significantly reduces the fatigue threshold at RT while increases the yield strength at elevated temperature (i.e., 650 °C). Increase in fatigue threshold is found to increase with the increase in block size, whereas high density of fine precipitates contributed the high temperature strength. The FCG rate in Paris regime was less affected by heat treatment and rolling process at RT. However, rolled sample tested at 650 °C shows a decrease in FCG rate as compared to heat-treated samples. In terms of higher FCG resistance and yield strength, warm rolling followed by tempering at low temperature (i.e., 700 °C) can be considered optimum. For the numerical simulation of FCG behavior, extended finite element method (XFEM) along with Paris law has been implemented in ABAQUS using python script. The simulated FCG behavior is found in good agreement with the experimental results.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Large Male Anthropomorphic Test Device (ATD) Finite Element Model (FEM) Correlation Improvement

As part of an activity for the Commercial Crew Program (CCP), a NASA Engineering and Safety Center (NESC) team conducted tests of small female, medium male, and large male Hybrid III anthropomorphic test device (ATD) head-neck complexes. The results were compared to simulation results from ATD finite element model (FEM) predictions to determine the validity and fidelity of ATD models for use in certifying CCP vehicles. The NESC was requested to help improve the large male ATD FEM to better represent physical responses. This report contains the outcome of the NESC assessment.

Kelly, Michael J.↗
DOE OSTI2013

A comparison of the lattice discrete particle method to the finite-element method and the K&C material model for simulating the static and dynamic response of concrete.

This report summarizes the work performed by the graduate student Jovanca Smith during a summer internship in the summer of 2012 with the aid of mentor Joe Bishop. The projects were a two-part endeavor that focused on the use of the numerical model called the Lattice Discrete Particle Model (LDPM). The LDPM is a discrete meso-scale model currently used at Northwestern University and the ERDC to model the heterogeneous quasi-brittle material, concrete. In the first part of the project, LDPM was compared to the Karagozian and Case Concrete Model (K&C) used in Presto, an explicit dynamics finite-element code, developed at Sandia National Laboratories. In order to make this comparison, a series of quasi-static numerical experiments were performed, namely unconfined uniaxial compression tests on four varied cube specimen sizes, three-point bending notched experiments on three proportional specimen sizes, and six triaxial compression tests on a cylindrical specimen. The second part of this project focused on the application of LDPM to simulate projectile perforation on an ultra high performance concrete called CORTUF. This application illustrates the strengths of LDPM over traditional continuum models.

Smith, Jovanca J.↗
DOE OSTI2014

An enriched finite element model with q-refinement for radiative boundary layers in glass cooling

Radiative cooling in glass manufacturing is simulated using the partition of unity finite element method. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary simplified P{sub 1} approximation for the radiation in non-grey semitransparent media. To integrate the coupled equations in time we consider a linearly implicit scheme in the finite element framework. A class of hyperbolic enrichment functions is proposed to resolve boundary layers near the enclosure walls. Using an industrial electromagnetic spectrum, the proposed method shows an immense reduction in the number of degrees of freedom required to achieve a certain accuracy compared to the conventional h-version finite element method. Furthermore the method shows a stable behaviour in treating the boundary layers which is shown by studying the solution close to the domain boundaries. The time integration choice is essential to implement a q-refinement procedure introduced in the current study. The enrichment is refined with respect to the steepness of the solution gradient near the domain boundary in the first few time steps and is shown to lead to a further significant reduction on top of what is already achieved with the enrichment. The performance of the proposed method is analysed for glass annealing in two enclosures where the simplified P{sub 1} approximation solution with the partition of unity method, the conventional finite element method and the finite difference method are compared to each other and to the full radiative heat transfer as well as the canonical Rosseland model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

FY23 Simulation of Elastic-Plastic Failure Propagation

This memo summarizes the simulation of ductile failure propagation work conducted under the ASC project “V&V of Ductile Failure” conducted during FY 23. Physically, the failure propagation consists of crack propagation in the material. In the numerical setting—specifically in a finite element model—propagation can be accomplished through element death when critical conditions occur locally at an element that is then deleted from the simulation. The validation of the finite element models is evaluated by direct comparison between the experimental and simulation results regarding the rate of crack growth and its influence on the load-deflection response of the specimens tested. This work considers two geometries that display stable crack propagation under displacement-controlled conditions. The first geometry consists of hat specimens loaded in compression with nominally identical geometries but made with three different materials: Steel A286, Al 7075-T651 and 304L stainless steel. The three materials represent a range of ductility values that affect the response and crack propagation within the specimen. The crack induced propagates under an essentially mode-II type of deformation. The second geometry consists of a pre-cracked 304L stainless steel compact tension test specimen loaded so as to induce a mode-I deformation at the crack.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array

Purpose: Feasibility of targeted and volumetric hyperthermia (40–45 °C) delivery to the prostate with a commercial MR-guided endorectal ultrasound phased array system, designed specifically for thermal ablation and approved for ablation trials (ExAblate 2100, Insightec Ltd.), was assessed through computer simulations and tissue-equivalent phantom experiments with the intention of fast clinical translation for targeted hyperthermia in conjunction with radiotherapy and chemotherapy. Methods: The simulations included a 3D finite element method based biothermal model, and acoustic field calculations for the ExAblate ERUS phased array (2.3 MHz, 2.3 × 4.0 cm{sup 2}, ∼1000 channels) using the rectangular radiator method. Array beamforming strategies were investigated to deliver protracted, continuous-wave hyperthermia to focal prostate cancer targets identified from representative patient cases. Constraints on power densities, sonication durations and switching speeds imposed by ExAblate hardware and software were incorporated in the models. Preliminary experiments included beamformed sonications in tissue mimicking phantoms under MR temperature monitoring at 3 T (GE Discovery MR750W). Results: Acoustic intensities considered during simulation were limited to ensure mild hyperthermia (T{sub max} < 45 °C) and fail-safe operation of the ExAblate array (spatial and time averaged acoustic intensity I{sub SATA} < 3.4 W/cm{sup 2}). Tissue volumes with therapeutic temperature levels (T > 41 °C) were estimated. Numerical simulations indicated that T > 41 °C was calculated in 13–23 cm{sup 3} volumes for sonications with planar or diverging beam patterns at 0.9–1.2 W/cm{sup 2}, in 4.5–5.8 cm{sup 3} volumes for simultaneous multipoint focus beam patterns at ∼0.7 W/cm{sup 2}, and in ∼6.0 cm{sup 3} for curvilinear (cylindrical) beam patterns at 0.75 W/cm{sup 2}. Focused heating patterns may be practical for treating focal disease in a single posterior quadrant of the prostate and diffused heating patterns may be useful for heating quadrants, hemigland volumes or even bilateral targets. Treatable volumes may be limited by pubic bone heating. Therapeutic temperatures were estimated for a range of physiological parameters, sonication duty cycles and rectal cooling. Hyperthermia specific phasing patterns were implemented on the ExAblate prostate array and continuous-wave sonications (∼0.88 W/cm{sup 2}, 15 min) were performed in tissue-mimicking material with real-time MR-based temperature imaging (PRFS imaging at 3.0 T). Shapes of heating patterns observed during experiments were consistent with simulations. Conclusions: The ExAblate 2100, designed specifically for thermal ablation, can be controlled for delivering continuous hyperthermia in prostate while working within operational constraints.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1994

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗
NASA NTRS2010

Development of a Aerothermoelastic-Acoustics Simulation Capability of Flight Vehicles

A novel numerical, finite element based analysis methodology is presented in this paper suitable for accurate and efficient simulation of practical, complex flight vehicles. An associated computer code, developed in this connection, is also described in some detail. Thermal effects of high speed flow obtained from a heat conduction analysis are incorporated in the modal analysis which in turn affects the unsteady flow arising out of interaction of elastic structures with the air. Numerical examples pertaining to representative problems are given in much detail testifying to the efficacy of the advocated techniques. This is a unique implementation of temperature effects in a finite element CFD based multidisciplinary simulation analysis capability involving large scale computations.

Gupta, K. K.↗
NASA NTRS1991

Microfracture in high temperature metal matrix crossply laminates

Microfracture (fiber/matrix fracture, interphase debonding and inter-ply delamination) in high temperature metal matrix composites (HTMMC), subjected to both mechanical and thermal loading, is computationally simulated. A crossply 0.3 fiber volume ratio SiC/Ti15 composite with 0/90/0 layup is evaluated for microfracture using a multicell finite element model. A computational simulation procedure based on strain energy release rates is used to predict the fracture process and establish the hierarchy of fracture modes. Microfracture results for various loading cases are presented and discussed.

Mital, Subodh K.↗
DOE OSTI2018

Thermal Wear Sensing System: Proof of the Concept

A technique for wear sensing using a sacrificial layer with specific material properties has been developed successfully. This is achieved by recording the temperature evolution underneath the sacrificial organic layer during sliding tests, and by correlating the recorded temperature with a wear depth. Finite element modeling was used to calculate the contact temperature which spreads through the sacrificial layer under different sliding test conditions. Numerical simulation was used to determine the optimal parameters for the best correlation between temperature and wear depth. These parameters were used to select the appropriate sacrificial layer for wear sensing. Thus, an epoxy layer was implemented and experimental measurements were carried out to check the correlation between the measured temperature underneath the epoxy layer and wear depth. It was found that a low thermal conductivity of the sacrificial layer of less than 1 W/m K and a relative high thermal conductivity of the substrate, namely more than 4 W/m K are the optimal parameters for our thermal wear sensing system. The thickness of the sacrificial layer is also a crucial parameter and has to be taken into account when designing such a sensing system.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Micro-taper as focusing or scattering optical element

We consider micro-taper (narrow refractive axicon) as optical element which is focusing or scattering in dependence on axicon’s cone angle. The diffraction of laser beam by micro-taper is simulated by two methods: multiply internal ray reflections using geometrical approach and Helmholtz equation solving using finite elements method. Based on ray optics we derive analytic formulas for conical angles values which provide focusing or scattering features of micro-taper. Numerical simulation by finite elements method verifies theoretical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Room Temperature Stamping of High-Strength Aluminum for Lightweight Structural Automotive Components

The use of aluminum (Al) 7xxx series in automobile structural applications is limited by factors such as its low room-temperature formability in T6 temper and the need for off-line artificial aging to regain high strength if it is formed in soft temper. In this study, two tempers were established that enable Al 7075 blanks to be room temperature stamped to create side impact beams that achieve near-T6 strength without post-forming artificial aging. It is shown that natural aging at room temperature, after solution treatment and before stamping, provides an alternative to high temperature pre-aging and hot stamping techniques in 7xxx alloys. Both tempers enabled room temperature stamping of the beams without cracking, indicating good formability. The Al beams were ~38% lighter than the benchmark boron steel beam. Further, the hardness in the low-deformation regions of the stamped Al beams were within 88% of a conventional 7075-T6 sheet. Paint-bake treatment of the beams appeared to cause over aging, lowering the ultimate tensile strength and ductility. Finite element three point bending simulations of the beams showed higher bending force and energy absorption before failure for the Al 7075 showing good promise for 7xxx series Al as a light-weight alternative to steel beams.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Correction to the Stress–Strain Curve During Multistage Hot Deformation of 7150 Aluminum Alloy Using Instantaneous Friction Factors

Multistage stress–strain curve correction based on an instantaneous friction factor was studied for axisymmetric uniaxial hot compression of 7150 aluminum alloy. Experimental friction factors were calculated based on continuous isothermal axisymmetric uniaxial compression tests at various deformation parameters. Then, an instantaneous friction factor equation was fitted by mathematic analysis. After verification by comparing single-pass flow stress correction with traditional average friction factor correction, the instantaneous friction factor equation was applied to correct multistage stress–strain curves. The corrected results were reasonable and validated by multistage relative softening calculations. This research provides a broad potential for implementing axisymmetric uniaxial compression in multistage physical simulations and friction optimization in finite element analysis.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Microfracture in high temperature metal matrix crossply laminates

Microfracture (fiber/matrix fracture, interphase debonding, and inter-ply delamination) in high temperature metal matrix composites (HTMMC), subjected to both mechanical and thermal loading, is computationally simulated. A crossply 0.3 fiber volume ratio SiC/Ti15 composite with 0/90/0 lay-up is evaluated for microfracture using a multicell finite element model. A computational simulation procedure based on strain energy release rates is used to predict the fracture process and establish the hierarchy of fracture modes. Microfracture results for various loading cases are presented and discussed.

Mital, Subodh K.↗
NASA NTRS2014

Simulating the Impact Response of Composite Airframe Components

In 2010, NASA Langley Research Center obtained residual hardware from the US Army's Survivable Affordable Repairable Airframe Program (SARAP). The hardware consisted of a composite fuselage section that was representative of the center section of a Black Hawk helicopter. The section was fabricated by Sikorsky Aircraft Corporation and designated the Test Validation Article (TVA). The TVA was subjected to a vertical drop test in 2008 to evaluate a tilting roof concept to limit the intrusion of overhead mass items, such as the rotor transmission, into the fuselage cabin. As a result of the 2008 test, damage to the hardware was limited primarily to the roof. Consequently, when the post-test article was obtained in 2010, the roof area was removed and the remaining structure was cut into six different types of test specimens including: (1) tension and compression coupons for material property characterization, (2) I-beam sections, (3) T-sections, (4) cruciform sections, (5) a large subfloor section, and (6) a forward framed fuselage section. In 2011, NASA and Sikorsky entered into a cooperative research agreement to study the impact responses of composite airframe structures and to evaluate the capabilities of the explicit transient dynamic finite element code, LS-DYNA®, to simulate these responses including damage initiation and progressive failure. Finite element models of the composite specimens were developed and impact simulations were performed. The properties of the composite material were represented using both a progressive in-plane damage model (Mat 54) and a continuum damage mechanics model (Mat 58) in LS-DYNA. This paper provides test-analysis comparisons of time history responses and the location and type of damage for representative I-beam, T-section, and cruciform section components.

Jackson, Karen E.↗
DOE OSTI2014

Revealing origin of quasi-one dimensional current transport in defect rich two dimensional materials

The presence of defects in graphene have for a long time been recognized as a bottleneck for its utilization in electronic and mechanical devices. We recently showed that micro four-point probes may be used to evaluate if a graphene film is truly 2D or if defects in proximity of the probe will lead to a non-uniform current flow characteristic of lower dimensionality. In this work, simulations based on a finite element method together with a Monte Carlo approach are used to establish the transition from 2D to quasi-1D current transport, when applying a micro four-point probe to measure on 2D conductors with an increasing amount of line-shaped defects. Clear 2D and 1D signatures are observed at low and high defect densities, respectively, and current density plots reveal the presence of current channels or branches in defect configurations yielding 1D current transport. A strong correlation is found between the density filling factor and the simulation yield, the fraction of cases with 1D transport and the mean sheet conductance. The upper transition limit is shown to agree with the percolation threshold for sticks. Finally, the conductance of a square sample evaluated with macroscopic edge contacts is compared to the micro four-point probe conductance measurements and we find that the micro four-point probe tends to measure a slightly higher conductance in samples containing defects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗