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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Study of Aircraft Structural Response and Occupant Loading During a Ditching Event Utilizing LS-DYNA Simulation

In this study, the capability of finite element model (FEM) analysis to assess aircraft structural response and occupant injury risk during aircraft water ditching conditions was evaluated. A FEM of the Fokker F28 aircraft, previously validated against full-scale crash test data, was simulated in a realistic water impact environment utilizing recorded conditions from the ditching of US Airways Flight 1549 into the Hudson River. The FEM of the F28 aircraft, was developed from a NASTRAN loads model provided by the vehicle manufacturer, calibrated using component fuselage drop-tests, and validated through test-analysis comparisons from a full-scale crash test of the F28 onto soil. Simulation results are presented in several categories: vehicle kinematics, structural accelerations, damage, occupant responses, and injury risk predictions. Results were consistent with the general outcomes of Flight 1549, indicating that the developed vehicle FEM, validated against land impact test data, may provide an effective tool for evaluating water ditching events.

Jacob Putnam↗
NASA NTRS1989

A NASTRAN/TREETOPS solution to a flexible, multi-body dynamics and controls problem on a UNIX workstation

Demands for nonlinear time history simulations of large, flexible multibody dynamic systems has created a need for efficient interfaces between finite-element modeling programs and time-history simulations. One such interface, TREEFLX, an interface between NASTRAN and TREETOPS, a nonlinear dynamics and controls time history simulation for multibody structures, is presented and demonstrated via example using the proposed Space Station Mobile Remote Manipulator System (MRMS). The ability to run all three programs (NASTRAN, TREEFLX and TREETOPS), in addition to other programs used for controller design and model reduction (such as DMATLAB and TREESEL, both described), under a UNIX Workstation environment demonstrates the flexibility engineers now have in designing, developing and testing control systems for dynamically complex systems.

Benavente, Javier E.↗
DOE OSTI2015

Improving the Accuracy of PCMI Simulations with More Realistic Geometry and Material Models

Assessment of fuel performance codes is an endeavour that requires re-evaluation as new simulation capabilities become available. BISON, the finite element-based fuel performance code in development at Idaho National Laboratory, has been used to simulate a variety of experiments. Calculations from these simulations are compared to experimental measurements for the purpose of evaluating the progress of code development. One such experiment, from the Risø 3 Fission Gas Project, provides measurements of fission gas release and clad diameter. A subset of the Risø3 Fission Gas Project, which is referred to here, consisted of a base-irradiation period of a full-length fuel rod in the Quad-Cities-1 BWR. The base irradiation test was followed by a “bump test”, in which a sub-section of the original rod was subjected to a strongly non-uniform axial power profile. The bump test took place in Risø test reactor DR3 in a water-cooled HP1 rig under BWR conditions. This experiment was chosen because it allows for an evaluation of several aspects of BISON, including fully coupled thermo-mechanics, contact, and several nonlinear material models. A summary of simulation results and comparison to the experiment data has been published previously, however, this paper shows new results with recently added features to BISON. These new features include representing the pellets as discrete entities rather than a smeared column, including primary creep in the clad creep model, and the coupling of fuel gaseous swelling to fission gas release. Comparing to previous results, the addition of theses features has narrowed the gap between BISON calculations and experimental measurements for final clad diameter and fission gas release.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Coupled nonlinear elasticity, plastic slip, twinning, and phase transformation in single crystal titanium for plate impact loading

High purity single crystal titanium (Ti) under shockwave loading in plate impact experiment is modelled and simulated, with help of the finite element. A thermodynamically consistent system of equations is formulated in the frame of large deformation to consider material anisotropy, rate dependence, and multi-physics including nonlinear elasticity, dislocation based plastic slip, deformation twinning, and phase transformation. A novel kinematics is proposed to consider phase transformation in the α twin variants and also to consider the dislocation based slip in all of components of parent material, primary twins and the ω phase. The thermodynamically consistent driving forces for plastic slip, twinning, and phase transformation are developed based upon the second law of thermodynamics. For nonlinear elasticity, the stiffness matrix is dependent on the volume fractions of all components, and the volumetric part of Cauchy stress is obtained from the nonlinear equation of state. Dislocation based plastic slip in multi-variant and multiphase heterogeneous materials is developed and interactions of dislocations among all slip modes are taken into account. A mechanism for dislocation density evolution during twinning and phase transformation is proposed. The shock loading along the [0001] and [10$\bar{1}1$ directions of single crystal high purity Ti is investigated computationally. Very good correspondence between simulation and experiment is obtained, which includes pole figures, volume fractions of components, free surface velocity, peak pressure, and phase transformation pressure. Multiple experimental phenomena are interpreted based upon the progression of dislocation slip, deformation twinning, and phase transformation. In compression with the [10$\bar{1}1$ crystal, a higher volume fraction in the primary twins but a lower secondary twin volume fraction in the [0001] crystal in experiment was observed. Finally, the higher propensity for phase transformation occurs in the [0001] crystal is reproduced. In addition to material texture, distributions of temperature, stresses, and plastic strains dependent on the impact loading directions are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Radiation and viscous dissipation effect on square porous annulus

The present study is carried out to investigate the effect of radiation and viscous dissipation in a square porous annulus subjected to outside hot T{sub h} and inside cold T{sub c} temperature. The square annulus has a hollow section of dimension D×D at the interior of annulus. The flow is assumed to obey Darcy law. The governing equations are non-dimensionalised and solved with the help of finite element method. Results are discussed with respect to viscous dissipation parameter, radiation parameter and size of the hollow section of annulus.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

An Integrated Numerical Method to Predict the Strain Sensing Behavior of Flexible Conductive Fibers under Electric–Thermal Conditions

Flexible conductive textiles have attracted an increasing amount of attention due to being flexible and stretchable with multi-functions, in addition to the strain sensing function. In this paper, an integrated numerical method is established to predict the strain sensing behavior of PPy-coated flexible conductive fibers. Based on our previous investigation on the strain sensing behavior of the PPy-coated elastic conductive fibers, it has been known that the main strain sensing mechanisms of the conductive fibers are micro-cracks opening under stretching and closing under unloading. Therefore, in the numerical approach, extended finite element method is used to simulate the multi-crack initiation and propagation. Then, by establishing a representative element model including a single crack, the resistance variation with the increase in crack width/depth can be obtained using the electric–thermal package in ABAQUS software. Combining the relationship between crack width and strain, the relationship between fractional increment in resistance (ΔR/R{sub 0}) and the strain can be obtained. The prediction of the fractional increment in resistance (ΔR/R{sub 0}) versus strain is consistent with the experimental results of PPy-coated Lycra fibers. The new integrated method is further verified by experimental results of PPy-coated XLA fibers, showing that the integrated numerical approach can capture the complicated strain sensing mechanisms of PPy-coated elastic fibers and predict the strain sensing behavior of flexible conductive fibers. This approach will make it possible to shorten the period, and reduce the cost in R&D of the PPy-coated electrically conductive fibers, comparing with the corresponding experimental work.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Computational simulation of intermingled-fiber hybrid composite behavior

Three-dimensional finite-element analysis and a micromechanics based computer code ICAN (Integrated Composite Analyzer) are used to predict the composite properties and microstresses of a unidirectional graphite/epoxy primary composite with varying percentages of S-glass fibers used as hydridizing fibers at a total fiber volume of 0.54. The three-dimensional finite-element model used in the analyses consists of a group of nine fibers, all unidirectional, in a three-by-three unit cell array. There is generally good agreement between the composite properties and microstresses obtained from both methods. The results indicate that the finite-element methods and the micromechanics equations embedded in the ICAN computer code can be used to obtain the properties of intermingled fiber hybrid composites needed for the analysis/design of hybrid composite structures. However, the finite-element model should be big enough to be able to simulate the conditions assumed in the micromechanics equations.

Mital, Subodh K.↗
DOE OSTI2024

NEAMS-Multiphysics MOOSE End Year Framework Activities FY20

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Stochastic-Strength-Based Damage Simulation of Ceramic Matrix Composite Laminates

The Finite Element Analysis-Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to characterize and predict the progressive damage response of silicon-carbide-fiber-reinforced reaction-bonded silicon nitride matrix (SiC/RBSN) composite laminate tensile specimens. Studied were unidirectional laminates [0] (sub 8), [10] (sub 8), [45] (sub 8), and [90] (sub 8); cross-ply laminates [0 (sub 2) divided by 90 (sub 2),]s; angled-ply laminates [plus 45 (sub 2) divided by -45 (sub 2), ]s; doubled-edge-notched [0] (sub 8), laminates; and central-hole laminates. Results correlated well with the experimental data. This work was performed as a validation and benchmarking exercise of the FEAMAC/CARES program. FEAMAC/CARES simulates stochastic-based discrete-event progressive damage of ceramic matrix composite and polymer matrix composite material structures. It couples three software programs: (1) the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), (2) the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program (CARES/Life), and (3) the Abaqus finite element analysis program. MAC/GMC contributes multiscale modeling capabilities and micromechanics relations to determine stresses and deformations at the microscale of the composite material repeating-unit-cell (RUC). CARES/Life contributes statistical multiaxial failure criteria that can be applied to the individual brittle-material constituents of the RUC, and Abaqus is used to model the overall composite structure. For each FEAMAC/CARES simulation trial, the stochastic nature of brittle material strength results in random, discrete damage events that incrementally progress until ultimate structural failure.

composite structures↗
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
DOE OSTI2018

Control rod treatment for FEM based radiation transport methods

This paper presents a novel control rod treatment for finite element method (FEM) based reactor physics simulations. The new method falls into the class of homogenization methods which treat material interfaces within mesh elements by defining a representative, homogenized nuclear cross section set. In contrast to existing homogenization based cusping treatments that attempt to approximate flux-volume weighting of the cross sections, we propose to exactly integrate terms appearing in the variational form of the transport or diffusion equation over mixed material elements. This goal is achieved by defining smoothly varying cross sections in elements that straddle material boundaries. It is shown that the smoothly varying cross sections can be obtained by projecting the piecewise set of cross sections onto the set of Legendre polynomials up to a finite order. Implementation of smoothly varying cross sections is supported within the Rattlesnake radiation multiphysics code that is currently under development at the Idaho National Laboratory. Homogenization methods traditionally suffer from cusping effects that originate from overestimation of absorption cross sections in strongly neutron-absorbing materials as found in control rods. We demonstrate that the presented method is (1) highly effective in removing cusping effects and (2) provides accurate results based on the C5G7-TD benchmark, a time-dependent version of the C5G7 benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

The semidiscrete Galerkin finite element modeling of compressible viscous flow past an airfoil

The primary project was the numerical simulation, by a finite element/finite difference method, of the viscous flow about an airfoil. The secondary project involved the numerical simulation of the three-dimensional separated and vortex-dominated flow about a hemispherically capped cylinder in the transonic regime. Preliminary calculations were started for the hemisphere-cylinder at 0 and 5 degree angle of attack. The solution of the flow field about airfoils and wings is required to determine the important parameters of lift, moment, and drag. Viscous effects must be accounted for if the drag is to be accurately calculated. At present there are basically two approaches to the numerical simulation of the flow field, the use of fully viscous models and the inviscid/viscous models. The fully viscous models require the solution of an approximation of the Navier-Stokes equations and therefore should simulate most of the physical mechanisms. A fast, accurate, and computationally efficient inviscid flow solver was recently developed by Hartwich. It is thought that Hartwich's program coupled to a fast, accurate, and computationally efficient boundary layer code, will make an excellent tool for airfoil design. The purpose of the primary project was to develop a compressible boundary layer code using the semidiscrete Galerkin finite element method. The numerical scheme employed used the combination of a Dorodnitsyn formulation of the boundary layer equations, with a finite difference/finite element procedure (semidiscrete Galerkin method), in the solution of the compressible two-dimensional boundary layer equations. A laminar compressible boundary layer code was developed and tested for a NACA 0012 airfoil at a Mach number of 0.5, a Reynolds number of 5000, and zero angle of attack. At present the boundary layer program solves up to, but not beyond, separation.

Meade, Andrew J.↗
DOE OSTI2018

A thermal-mechanical finite element workflow for directed energy deposition additive manufacturing process modeling

This work proposes a finite element (FE) analysis workflow to simulate directed energy deposition (DED) additive manufacturing at a macroscopic length scale (i.e. part length scale) and to predict thermal conditions during manufacturing, as well as distortions, strength and residual stresses at the completion of manufacturing. The proposed analysis method incorporates a multi-step FE workflow to elucidate the thermal and mechanical responses in laser engineered net shaping (LENS) manufacturing. For each time step, a thermal element activation scheme captures the material deposition process. Then, activated elements and their associated geometry are analyzed first thermally for heat flow due to radiation, convection, and conduction, and then mechanically for the resulting stresses, displacements, and material property evolution. Finally, simulations agree with experimentally measured in situ thermal measurements for simple cylindrical build geometries, as well as general trends of local hardness distribution and plastic strain accumulation (represented by relative distribution of geometrically necessary dislocations).

42 ENGINEERING↗
DOE OSTI2020

Effect of heat treatment on microstructural evolution and hardness homogeneity in laser powder bed fusion of alloy 718

The thermal history developed in laser powder bed fusion (LPBF) processes has been shown to be complex resulting in equally complex microstructures and mechanical properties. Microscopic observations and Vickers micro-hardness mapping measurements were carried out on diff ;erent section planes of LPBF alloy 718 cuboids. Three-dimensional finite element analysis was used to simulate thermal history and to predict the residual stress distribution in the as-built material. Computational thermodynamics was used to predict the micro-segregation and nucleation driving force of various phases in the bulk and in segregated regions. Varied heat-treatments such as simulated hot isostatic pressing, and double aging were applied. Their influence on the microstructure, microsegregation, precipitate formation, and micro-hardness variations of LPBF alloy 718 were investigated. Hardness map results showed heterogeneous micro-hardness on the xy- and xz-planes of the as-built parts where the bottom plane and center regions had larger hardness of ~315 HV 0.5 while the top plane and contours showed hardness of ~300 HV 0.5 . It was found that the aging treatment increased the overall hardness of the as-built condition from ~310 HV 0.5 to 470 HV 0.5 but also increased the hardness gradient throughout the coupon. After simulated hot isostatic pressing process (i.e., without applied pressure) at 1020 °C for 4 h followed by water quench (HIPWQ), the hardness gradient and hardness was minimized (~210 HV 0.5 ) as the microstructure transitioned from heterogeneous columnar grains in the as-built condition to more uniform recrystallized grains. A double aging treatment was applied to enhance hardness from ~210 HV 0.5 to ~440 HV 0.5 . Finally, HIPWQ followed by double aging produced a homogeneous microstructure and more uniform hardness map with enhanced mechanical properties in LPBF alloy 718 coupons.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions

A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.

42 ENGINEERING↗
DOE OSTI2026

On the roles of welding residual stresses in determination of fracture toughness in austenitic stainless steel SUS 304 pipeline girth welds

Welding residual stresses especially the high tensile stresses are proved to have negative impacts on the fatigue and fracture behaviors of welded structures. In this study, a virtual fabrication of test specimens from welding process to specimen preparation was carried out by numerical simulation. An austenitic stainless steel multi-pass pipe welding was simulated by transient thermal–mechanical finite element analysis, the residual stresses were then mapped into the test specimen to evaluate fracture toughness. The findings in this study confirmed that, residual stress can be high in a sub-sized compact tensile specimen, which may accelerate or hinder the crack propagation during actual fatigue and fracture tests as reported in recent years. The influence of the cutting location and orientation of the specimen on fracture performance was investigated systematically to provide a fundamental understanding of welding residual stress and necessary insights into the specimen preparation procedure. Considering the limitation of measuring techniques and the complexity of the stress distribution, the developed numerical model can be a very useful tool to elucidate the stress evolution and quantify the effect of remaining welding stress on fracture toughness.

Fracture behavior↗
DOE OSTI2016

Mitigating cutting-induced plasticity in the contour method, Part 2: Numerical analysis

Cutting-induced plasticity can have a significant effect on the measurement accuracy of the contour method. The present study examines the benefit of a double-embedded cutting configuration that relies on self-restraint of the specimen, relative to conventional edge-crack cutting configurations. A series of finite element analyses are used to simulate the planar sectioning performed during double-embedded and conventional edge-crack contour cutting configurations. The results of numerical analyses are first compared to measured results to validate the cutting simulations. The simulations are then used to compare the efficacy of different cutting configurations by predicting the deviation of the residual stress profile from an original (pre-cutting) reference stress field, and the extent of cutting-induced plasticity. Comparisons reveal that while the double-embedded cutting configuration produces the most accurate residual stress measurements, the highest levels of plastic flow are generated in this process. As a result, this cutting-induced plastic deformation is, however, largely confined to small ligaments formed as a consequence of the sample sectioning process, and as such it does not significantly affect the back-calculated residual stress field.

36 MATERIALS SCIENCE↗