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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Characterization of the arc in crossflow using a two-temperature nonequilibrium plasma flow model

Diverse industrial applications such as circuit breakers and wire arc spraying involve the interaction between an electric arc and a stream of gas impinging perpendicular to it, a configuration commonly referred to as the arc in crossflow. The arc in crossflow is simulated using a three-dimensional time-dependent two-temperature (heavy-species and electrons) plasma flow model to better capture plasma-gas interactions and deviations from Local Thermodynamic Equilibrium (LTE). The coupled fluid-electromagnetic flow model is solved in a monolithic manner using Variational Multiscale Finite Element Method. Simulation findings are validated with experimental findings and contrasted against results obtained with a LTE model. Results from the two-temperature model corroborate experimental observations while providing quantification of the deviation between heavy-species and electron temperatures. The model is used to characterize the arc in crossflow as a function of the Reynolds and Enthalpy dimensionless numbers, which encapsulate the inter-dependence among the main parameters total current, inflow velocity, and inter-electrode spacing. Furthermore, the characterization revealed the behavior of arc shape, voltage drop, arc power, the degree of nonequilibrium, as well as the characteristic plasma front thickness, with varying controlling parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Towards performance-portability of the Albany Finite Element analysis code using the Kokkos library of Trilinos

Performance portability on heterogeneous high-performance computing (HPC) systems is a major challenge faced today by code developers: parallel code needs to execute correctly as well as with high performance on machines with different architectures, operating systems, and software libraries. The Finite Element Method (FEM) is a popular and flexible method for discretizing partial differential equations arising in a wide variety of scientific, engineering, and industry applications that require HPC. This paper presents some preliminary results pertaining to our development of a performance portable implementation of the FEM-based Albany code. Performance portability is achieved using the Kokkos library of Trilinos. We present performance results for two different physics simulations modules in Albany: the Aeras global atmosphere dynamical code and the FELIX land-ice solver. As a result, numerical experiments show that our single code implementation gives reasonable performance across two multi-core/many-core architectures: NVIDIA GPUs and multi-core CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Modeling of the Heat-Affected and Thermomechanically Affected Zones in a Ti-6Al-4V Inertia Friction Weld

Inertia friction welding has been used across the aerospace, automotive, and power-generation industries for the fabrication of complex axisymmetric components for over forty years. The process involves one axisymmetric piece being held stationary and another piece being brought into contact set to rotate about its axis of symmetry by a flywheel with the system under an applied load across the joint. Plasticization at the joint interface through the frictional heating sees the two pieces bond together. The titanium alloy Ti-6Al-4V has been widely studied for inertia welding applications. A successful selection of processing parameters (flywheel energy and mass, applied load) allows an inertia welding process which produces a very high-integrity weld, with a minimal heat-affected zone (HAZ) and thermomechanically affected zone (TMAZ), formed as a narrow band at the interface and extending further into the material. The width of this narrow band of heated material is dependent upon the process parameters used. A series of experimental inertia friction welds were performed using Ti-6Al-4V, and a finite element (FE) modeling framework was developed using the FE code Deform in order to predict the widths of the HAZ and TMAZ at the weld interface. The experimentally observed HAZ boundaries were correlated with the thermal fields from the FE model, while TMAZ boundaries were correlated with the Von Mises plastic strain fields.

42 ENGINEERING↗
DOE OSTI2016

Seismic and LOCA safety testing and analysis program for Westinghouse PWR fuel under end of life conditions

The NRC issued Information Notice - IN 2012-09 in June 2012 which challenges approved assumptions based on the beginning of life (BOL) conditions in seismic and LOCA analysis. Westinghouse manufactures a variety of PWR fuel designs, including 14*14, 15*15, 16*16 and 17*17 array, also different fuel lengths from 12 feet to 14 feet. The fuels are supplied to the cores designed by Westinghouse and Combustion Engineering, and other cores around world. To address IN 2012-09 efficiently and economically, Westinghouse and PWR Owner Groups launched a testing and analysis program funded by PWROG. The program planed testing and analysis matrices based on types of fuel designs. The testing plan includes grid crush tests, fuel assembly mechanical tests and flowing water damping tests. Applying modelling technique, the program will reduce the number of test and cost. The ultimate goal of the program is to generate fuel assembly EOL models for all fuel designs manufactured by Westinghouse and to be ready for any plant specific seismic and LOCA analysis to address NRC IN for EOL effect. The phase I of PWROG for EOL conditions has been completed and several conclusions are obtained. First, the gapped grids reduce grid dynamic stiffness, local impact stiffness and crush strength. Also, due to gapped grids, both the fuel assembly natural frequency and stiffness are reduced. Secondly, applying higher fuel assembly damping due to the flowing water effect can significantly reduce fuel assembly impact forces. However, during a seismic event the reactor coolant pumps may trip resulting in a coast-down. As such, reduced flow rates that occur during the pump coast-down must be considered. Thirdly, the analytical model representation of a reactor core adequately characterizes the core response under seismic and LOCA excitations. Static and dynamic structural test results of a fuel assembly compared with analytical predictions indicate that the fuel assembly lateral and axial structural behavior can be fully determined with finite element models. Fourthly, analyses for EOL fuel assemblies under seismic and LOCA loads demonstrate EOL effect and flowing water damping effect. With flowing water damping benefit, significant margin can be recovered for fuel assemblies with EOL conditions and under seismic loads.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

A Well-Posed Two Phase Flow Model and its Numerical Solutions for Reactor Thermal-Fluids Analysis

A 7-equation two-phase flow model and its numerical implementation is presented for reactor thermal-fluids applications. The equation system is well-posed and treats both phases as compressible flows. The numerical discretization of the equation system is based on the finite element formalism. The numerical algorithm is implemented in the next generation RELAP-7 code (Idaho National Laboratory (INL)’s thermal-fluids code) built on top of an other INL’s product, the massively parallel multi-implicit multi-physics object oriented code environment (MOOSE). Some preliminary thermal-fluids computations are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Mechanics of cantilever beam: Implementation and comparison of FEM and MLPG approach

Two weak form solution approaches for partial differential equations, the well known meshbased finite element method and the newer meshless local Petrov Galerkin method are described and compared on a standard test case - mechanics of cantilever beam. The implementation, solution accuracy and calculation complexity are addressed for both approaches. We found out that FEM is superior in most standard criteria, but MLPG has some advantages because of its flexibility that results from its general formulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

Numerical Modeling of Ocular Dysfunction in Space

Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.

intraocular pressure↗
NASA NTRS1987

Initial dynamic load estimates during configuration design

This analysis includes the structural response to shock and vibration and evaluates the maximum deflections and material stresses and the potential for the occurrence of elastic instability, fatigue and fracture. The required computations are often performed by means of finite element analysis (FEA) computer programs in which the structure is simulated by a finite element model which may contain thousands of elements. The formulation of a finite element model can be time consuming, and substantial additional modeling effort may be necessary if the structure requires significant changes after initial analysis. Rapid methods for obtaining rough estimates of the structural response to shock and vibration are presented for the purpose of providing guidance during the initial mechanical design configuration stage.

Schiff, Daniel↗
DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗
DOE OSTI2025

Quantum State Transfer between Heterogeneous 3D QPUs

We present detailed analyses and microwave simulations of a quantum interconnect architecture that couples two distinct three-dimensional quantum processing units (QPUs): (i) a two-cell superconducting radio-frequency (SRF) cavity coupled to multiple transmon qubits and (ii) a waveguide-based cavity hosting a planar circuit with transmons and resonators. The modules are linked via a low-loss superconducting coaxial cable, which serves as a quantum bus. Finite-element electromagnetic simulations are performed to optimize coupling strengths and mode frequencies for high-fidelity quantum state transfer between cavities. These simulation results provide a comprehensive blueprint for the experimental realization of a flexible quantum interconnect architecture, coupling heterogeneous QPUs.

Yaker, Joey [Fermilab]↗
DOE OSTI2016

Laser-excited optical emission response of CdTe quantum dot/polymer nanocomposite under shock compression

Laser-driven shock compression experiments and corresponding finite element method simulations are carried out to investigate the blueshift in the optical emission spectra under continuous laser excitation of a dilute composite consisting of 0.15% CdTe quantum dots by weight embedded in polyvinyl alcohol polymer. This material is a potential candidate for use as internal stress sensors. The analyses focus on the time histories of the wavelength blue-shift for shock loading with pressures up to 7.3 GPa. The combined measurements and calculations allow a relation between the wavelength blueshift and pressure for the loading conditions to be extracted. It is found that the blueshift first increases with pressure to a maximum and subsequently decreases with pressure. This trend is different from the monotonic increase of blueshift with pressure observed under conditions of quasistatic hydrostatic compression. Additionally, the blueshift in the shock experiments is much smaller than that in hydrostatic experiments at the same pressure levels. The differences in responses are attributed to the different stress states achieved in the shock and hydrostatic experiments and the time dependence of the mechanical response of the polymer in the composite. The findings offer a potential guide for the design and development of materials for internal stress sensors for shock conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Thermal behavior characterization of a kilowatt-power-level cryogenically cooled Yb:YAG active mirror laser amplifier

We present an analysis of the thermal behavior of a high-energy kilowatt-average-power diode-pumped cryogenically cooled Yb:YAG active mirror laser amplifier based on measurements and simulations. Maps of the temperature distribution of the laser material at pump powers up to 1 kW were obtained for the first time by spatially and spectrally resolving the fluorescence induced by a scanning probe beam. The wavefront distortion resulting from the front surface deformation and the overall deformation of the gain medium assembly were measured using a Mach–Zehnder interferometer. The measured deformations agree well with the results of thermomechanical modeling using finite element method simulations, and with the results of focal length shift measurements. The relative contributions to the optical path difference (OPD) of the mechanical deformations, refractive index changes, and electronic contribution are discussed. We propose that the Cr 4+ :YAG cladding plays a significant role in both the temperature distribution and the overall OPD changes. Moreover, the pump-induced mechanical deformations of the assembly dominate the OPD changes in this kilowatt-average-pump-power cryogenically cooled Yb:YAG active mirror laser.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Simulation of Aircraft Landing Gears with a Nonlinear Dynamic Finite Element Code

Recent advances in computational speed have made aircraft and spacecraft crash simulations using an explicit, nonlinear, transient-dynamic, finite element analysis code more feasible. This paper describes the development of a simple landing gear model, which accurately simulates the energy absorbed by the gear without adding substantial complexity to the model. For a crash model, the landing gear response is approximated with a spring where the force applied to the fuselage is computed in a user-written subroutine. Helicopter crash simulations using this approach are compared with previously acquired experimental data from a full-scale crash test of a composite helicopter.

Lyle, Karen H.↗
DOE OSTI2019

Strain localization and dynamic recrystallization in polycrystalline metals: Thermodynamic theory and simulation framework

We describe a theoretical and computational framework for adiabatic shear banding (ASB) and dynamic recrystallization (DRX) in polycrystalline materials. The Langer-Bouchbinder-Lookman (LBL) thermodynamic theory of polycrystalline plasticity, which we recently reformulated to describe DRX via the inclusion of the grain boundary density or the grain size as an internal state variable, provides a convenient and self-consistent way to represent the viscoplastic and thermal behavior of the material, with minimal ad-hoc assumptions regarding the initiation of yielding or onset of shear banding. Here, we implement the LBL-DRX theory in conjunction with a finite-element computational framework. Favorable comparison to experimental measurements on a top-hat AISI 316L stainless steel sample compressed with a split-Hopkinson pressure bar suggests the accuracy and usefulness of the LBL-DRX framework, and demonstrates the crucial role of DRX in strain localization.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Analysis of the static magnetic field-dependent optical transmission of Ni nanorod colloidal suspensions

The magnetic field-dependent optical transmission of dilute Ni nanorod aqueous suspensions was investigated. A series of four samples of nanorods were synthesized using the AAO template method and processed to stable colloids. The distributions of their length and diameter were characterized by analysis of TEM images and revealed average diameters of ∼25 nm and different lengths in the range of 60 nm–1100 nm. The collinear magnetic and optical anisotropy was studied by static field-dependent transmission measurements of linearly polarized light parallel and perpendicular to the magnetic field direction. The experimental results were modelled assuming the field-dependent orientation distribution function of a superparamagnetic ensemble for the uniaxial ferromagnetic nanorods in liquid dispersion and extinction cross sections for longitudinal and transversal optical polarization derived from different approaches, including the electrostatic approximation and the separation of variables method, both applied to spheroidal particles, as well as finite element method simulations of spheroids and capped cylindrical particles. The extinction cross sections were compared to reveal the differences associated with the approximations of homogeneous polarization and/or particle shape. The consequences of these approximations for the quantitative analysis of magnetic field-dependent optical transmission measurements were investigated and a reliable protocol derived. Furthermore, the changes in optical cross sections induced by electromagnetic interaction between two nanorods in parallel end-to-end and side-by-side configuration as a function of their separation were studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

Efficient Modal Basis Selection Criteria for Reduced-Order Nonlinear Simulation

A modal basis selection technique for a reduced-order nonlinear numerical simulation with application to two-dimensional structures is presented as a two-step procedure. A system identification analysis is first performed using proper orthogonal decomposition. Using these results, a set of load-invariant bases consisting of the normal modes is next selected. Two criteria for making the basis selection are offered; one using the modal assurance criterion and the other using the modal expansion theorem. The quality of the subsequent reduced-order analyses are examined through comparison with computationally intensive finite element nonlinear simulations in physical degrees-of-freedom. A clamped flat isotropic plate under a random acoustic loading is considered to demonstrate the procedure. It is found that the subject procedure enables formation of an accurate and computationally efficient reduced-order system applicable to a broad range of loading conditions.

Przekop, Adam↗
DOE OSTI2019

Influence of Single-Point Incremental Force Process Parameters on Forming Characteristics and Microstructure Evolution of AA-6061 Alloy Sheet

Single-point incremental forming (SPIF) process comprises a set of process variables, such as forming tool diameter, vertical pitch, spindle rotation speed, and tool velocity, which may affect the forming behavior of the sheet to be deformed. The objective of this work is to study the effect of SPIF process parameters on forming characteristics and microstructure development for the AA-6061 (T6) aluminum alloy sheet. The SPIF experiments and finite element (FE) simulations were performed at different process parameters to achieve conical shapes from the AA-6061 blanks. The effect of process parameters on forming forces, thickness uniformity in wall region of the cone and surface roughness of the blank was analyzed. A detailed microstructure study was performed to analyze the effect of process variables on microstructure and texture evolution during the SPIF process. This study reveals that the process parameters are likely to influence the texture development especially at high tool diameter and vertical pitch values. Therefore, suitability and consequences of using different combinations of tool diameter and vertical pitch values are discussed for AA-6061 alloy.

36 MATERIALS SCIENCE↗