SEARCH · Engineering Papers

Results for “Finite element methods”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

NASA NTRS1995

P-Finite-Element Program For Analysis Of Plates

BUCKY is p-finite-element computer program for highly accurate analysis of structures. Used to analyze buckling, bending, and in-plane stress-and-strain behaviors of plates. Provides elastic-plastic solutions for isotropic plates in states of plane stress, and axisymmetric solution sequence used to treat three-dimensional problems. Computes response of plate to variety of loading and boundary conditions by use of higher-order displacement function in p-finite-element method. Enables user to obtain results more accurate than obtained by use of traditional h-finite elements. Written in FORTRAN 77.

Smith, James P.↗

NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗

NASA NTRS1977

A finite element simulation of sound attenuation in a finite duct with a peripherally variable liner

Using multimodal analysis, a variational finite element method is presented for analyzing sound attenuation in a three-dimensional finite duct with a peripherally variable liner in the absence of flow. A rectangular element, with cubic shaped functions, is employed. Once a small portion of a peripheral liner is removed, the attenuation rate near the frequency where maximum attenuation occurs drops significantly. The positioning of the liner segments affects the attenuation characteristics of the liner. Effects of the duct termination are important in the low frequency ranges. The main effect of peripheral variation of the liner is a broadening of the attenuation characteristics in the midfrequency range. Because of matrix size limitations of the presently available computer program, the eigenvalue equations should be solved out of core in order to handle realistic sources.

Watson, W. R.↗

DOE OSTI2015

Temperature Profile in Fuel and Tie-Tubes for Nuclear Thermal Propulsion Systems

A finite element method to calculate temperature profiles in heterogeneous geometries of tie-tube moderated LEU nuclear thermal propulsion systems and HEU designs with tie-tubes is developed and implemented in MATLAB. This new method is compared to previous methods to demonstrate shortcomings in those methods. Typical methods to analyze peak fuel centerline temperature in hexagonal geometries rely on spatial homogenization to derive an analytical expression. These methods are not applicable to cores with tie-tube elements because conduction to tie-tubes cannot be accurately modeled with the homogenized models. The fuel centerline temperature directly impacts safety and performance so it must be predicted carefully. The temperature profile in tie-tubes is also important when high temperatures are expected in the fuel because conduction to the tie-tubes may cause melting in tie-tubes, which may set maximum allowable performance. Estimations of maximum tie-tube temperature can be found from equivalent tube methods, however this method tends to be approximate and overly conservative. A finite element model of heat conduction on a unit cell can model spatial dependence and non-linear conductivity for fuel and tie-tube systems allowing for higher design fidelity of Nuclear Thermal Propulsion.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2016

Energy dissipation via acoustic emission in ductile crack initiation

In this article, we present a modeling approach to estimate the energy release due to ductile crack initiation in conjunction to the energy dissipation associated with the formation and propagation of transient stress waves typically referred to as acoustic emission. To achieve this goal, a ductile fracture problem is investigated computationally using the finite element method based on a compact tension geometry under Mode I loading conditions. To quantify the energy dissipation associated with acoustic emission, a crack increment is produced given a pre-determined notch size in a 3D cohesive-based extended finite element model. The computational modeling methodology consists of defining a damage initiation state from static simulations and linking such state to a dynamic formulation used to evaluate wave propagation and related energy redistribution effects. The model relies on a custom traction separation law constructed using full field deformation measurements obtained experimentally using the digital image correlation method. The amount of energy release due to the investigated first crack increment is evaluated through three different approaches both for verification purposes and to produce an estimate of the portion of the energy that radiates away from the crack source in the form of transient waves. The results presented herein propose an upper bound for the energy dissipation associated to acoustic emission, which could assist the interpretation and implementation of relevant nondestructive evaluation methods and the further enrichment of the understanding of effects associated with fracture.

42 ENGINEERING↗

DOE OSTI2015

A Novel Transverse Flux Machine for Vehicle Traction Applications

A novel transverse flux machine topology for electric vehicle traction application using ferrite magnets is presented in this paper. The proposed transverse flux topology utilizes novel magnet arrangements in the rotor that are similar to Halbach-array to boost flux linkage; on the stator side, cores are alternately arranged around a pair of ring windings in each phase to make use of the entire rotor flux that eliminates end windings. Analytical design considerations and finite element methods are used for an optimized design of a scooter in-wheel motor. Simulation results from Finite Element Analysis (FEA) show the motor achieved comparable torque density to conventional rare-earth permanent magnet machines. This machine is a viable candidate for direct drive applications with low cost and high torque density.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗

DOE OSTI2015

Novel Transverse Flux Machine for Vehicle Traction Applications: Preprint

A novel transverse flux machine topology for electric vehicle traction applications using ferrite magnets is presented in this paper. The proposed transverse flux topology utilizes novel magnet arrangements in the rotor that are similar to the Halbach array to boost flux linkage; on the stator side, cores are alternately arranged around a pair of ring windings in each phase to make use of the entire rotor flux that eliminates end windings. Analytical design considerations and finite-element methods are used for an optimized design of a scooter in-wheel motor. Simulation results from finite element analysis (FEA) show that the motor achieved comparable torque density to conventional rare-earth permanent magnet (PM) machines. This machine is a viable candidate for direct-drive applications with low cost and high torque density.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2014

Surface Remeshing with Robust High-Order Reconstruction

Remeshing is an important problem in variety of applications, such as finite element methods and geometry processing. Surface remeshing poses some unique challenges, as it must deliver not only good mesh quality but also good geometric accuracy. For applications such as finite elements with high-order elements (quadratic or cubic elements), the geometry must be preserved to high-order (third-order or higher) accuracy, since low-order accuracy may undermine the convergence of numerical computations. The problem is particularly challenging if the CAD model is not available for the underlying geometry, and is even more so if the surface meshes contain some inverted elements. We describe remeshing strategies that can simultaneously produce high-quality triangular meshes, untangling mildly folded triangles and preserve the geometry to high-order of accuracy. Our approach extends our earlier works on high-order surface reconstruction and mesh optimization by enhancing its robustness with a geometric limiter for under-resolved geometries. We also integrate high-order surface reconstruction with surface mesh adaptation techniques, which alter the number of triangles and nodes. We demonstrate the utilization of our method to meshes for high-order finite elements, biomedical image-based surface meshes, and complex interface meshes in fluid simulations.

Ray, Navamita↗

NASA NTRS1986

Unstructured multigrid methods for the compressible Euler equations

A two-step explicit finite element based solution procedure for the compressible Euler equations is described. Convergence to steady state can be accelerated by using a multigrid technique. The geometric flexibility of the finite element method is retained by using a sequence of unnested grids. The viability of this approach is demonstrated for elliptic problems and initial experiences with the technique applied to the Euler equations are reported.

Loehner, R.↗

NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

DOE OSTI2020

Toward large-scale simulation of residual stress and distortion in wire and arc additive manufacturing

In this study, we aim to advance the structural analysis of wire and arc additive manufacturing (WAAM) by considering the thermomechanical features inherent in direct energy deposition. Simulation approaches including the iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit finite element method (FEM) were developed for accelerating additive manufacturing stress analysis that is very time consuming by conventional numerical methods. ISM and DMRM take advantage of the strong nonlinearity phenomenon near the moving heat source by reducing the global iterations and refining the local mesh, respectively. In addition to the use of GPUs, the explicit FEM is accelerated by a time scaling technique based on the inherent strain concept. The residual stress and distortion of two large builds were analyzed, showing very consistent numerical results and good agreement with experiments. Compared with the commercial software Abaqus, the novel approaches reduced the computational cost substantially without compromising accuracy. The GPU code showed the highest computational efficiency (30~70×), while DMRM and ISM had acceleration factors of 9× and 3×, respectively. Finally, such high-fidelity modeling approaches will be very useful for building up a digital twin of WAAM to reduce development time and cost.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗

DOE OSTI2019

High-Frequency Acoustic Mode Identification of Unstable Combustors

Abstract High frequency thermoacoustic instabilities are becoming increasingly problematic in modern combustion systems. Understanding which acoustic mode is being excited is important for understanding potential mechanisms and control approaches — for example, influence of a helical shear layer mode on the flame has profoundly different effects on the first tangential acoustic mode, than a higher order axial mode of similar frequency. Nonetheless, the modal density increases with frequency and it becomes increasingly difficult to determine which acoustic mode is self-excited, based upon frequency calculations alone. Moreover, access issues and cost usually limit the number of pressure probes that can be distributed axially and azimuthally in the combustor. This paper presents a methodology for identifying the acoustic mode by using high temperature pressure transducers flush mounted in a combustion chamber. Modal identification is demonstrated with a siren under non-reacting conditions. The siren is mounted on the chamber to excite longitudinal and azimuthal waves. Five acoustic sensors at different axial and azimuthal locations measure the pressure fluctuations simultaneously. Given the forcing frequency and the speed of sound, the pressure distribution in the combustor is reconstructed in the time domain from the measured data by using a least squares method to determine its mode shapes. In addition, the finite element method (FEM) solver is used to provide the eigenfrequencies and corresponding mode shapes. The test results demonstrate that the mode shapes from the reconstructed data and corresponding frequencies are consistent with those predicted from the FEM, which validates the methodology in this study. In addition, the methodology is extended to practical reacting cases without the siren to determine the acoustic mode shapes of naturally occurring instabilities. In these cases, the modal features have strong stochastic features, such as what appear to be stochastic variations in overall amplitude and relative amplitudes of clockwise and counterclockwise waves.

Kim, J.↗

NASA NTRS1983

Design, analysis and test of composite curved frames for helicopter fuselage structure

Curved beam effects in composite frame structures representative of a light helicopter airframe are examined, and currently available analytical methods for studying these effects are briefly reviewed. A finite element study of curved composite frames is then reported, and it is shown that the curved frame effects must be accurately accounted for to avoid premature fracture. The finite element method is shown to be accurate to within 10 percent in accounting for the curved beam effects in composite structures.

Rich, M. J.↗

NASA NTRS1986

Iterative methods for elliptic finite element equations on general meshes

Iterative methods for arbitrary mesh discretizations of elliptic partial differential equations are surveyed. The methods discussed are preconditioned conjugate gradients, algebraic multigrid, deflated conjugate gradients, an element-by-element techniques, and domain decomposition. Computational results are included.

Nicolaides, R. A.↗

NASA NTRS1988

Spectral element methods: Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered: the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗