SEARCH · Engineering Papers

Results for “Finite element”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Finite element thermal analysis of convectively-cooled aircraft structures

The design complexity and size of convectively-cooled engine and airframe structures for hypersonic transports necessitate the use of large general purpose computer programs for both thermal and structural analyses. Generally thermal analyses are based on the lumped-parameter finite difference technique, and structural analyses are based on the finite element technique. Differences in these techniques make it difficult to achieve an efficient interface. It appears, therefore, desirable to conduct an integrated analysis based on a common technique. A summary is provided of efforts by NASA concerned with the development of an integrated thermal structural analysis capability using the finite element method. Particular attention is given to the development of conduction/forced-convection finite element methodology and applications which illustrate the capabilities of the developed concepts.

Wieting, A. R.↗

NASA NTRS1991

A three-dimensional finite-element thermal/mechanical analytical technique for high-performance traveling wave tubes

Current research in high-efficiency, high-performance traveling wave tubes (TWT's) has led to the development of novel thermal/ mechanical computer models for use with helical slow-wave structures. A three-dimensional, finite element computer model and analytical technique used to study the structural integrity and thermal operation of a high-efficiency, diamond-rod, K-band TWT designed for use in advanced space communications systems. This analysis focused on the slow-wave circuit in the radiofrequency section of the TWT, where an inherent localized heating problem existed and where failures were observed during earlier cold compression, or 'coining' fabrication technique that shows great potential for future TWT development efforts. For this analysis, a three-dimensional, finite element model was used along with MARC, a commercially available finite element code, to simulate the fabrication of a diamond-rod TWT. This analysis was conducted by using component and material specifications consistent with actual TWT fabrication and was verified against empirical data. The analysis is nonlinear owing to material plasticity introduced by the forming process and also to geometric nonlinearities presented by the component assembly configuration. The computer model was developed by using the high efficiency, K-band TWT design but is general enough to permit similar analyses to be performed on a wide variety of TWT designs and styles. The results of the TWT operating condition and structural failure mode analysis, as well as a comparison of analytical results to test data are presented.

Bartos, Karen F.↗

DOE OSTI2016

Partition of unity finite element method for quantum mechanical materials calculations

The current state of the art for large-scale quantum-mechanical simulations is the planewave (PW) pseudopotential method, as implemented in codes such as VASP, ABINIT, and many others. However, since the PW method uses a global Fourier basis, with strictly uniform resolution at all points in space, it suffers from substantial inefficiencies in calculations involving atoms with localized states, such as first-row and transition-metal atoms, and requires significant nonlocal communications, which limit parallel efficiency. Real-space methods such as finite-differences (FD) and finite-elements (FE) have partially addressed both resolution and parallel-communications issues but have been plagued by one key disadvantage relative to PW: excessive number of degrees of freedom (basis functions) needed to achieve the required accuracies. In this paper, we present a real-space partition of unity finite element (PUFE) method to solve the Kohn–Sham equations of density functional theory. In the PUFE method, we build the known atomic physics into the solution process using partition-of-unity enrichment techniques in finite element analysis. The method developed herein is completely general, applicable to metals and insulators alike, and particularly efficient for deep, localized potentials, as occur in calculations at extreme conditions of pressure and temperature. Full self-consistent Kohn–Sham calculations are presented for LiH, involving light atoms, and CeAl, involving heavy atoms with large numbers of atomic-orbital enrichments. We find that the new PUFE approach attains the required accuracies with substantially fewer degrees of freedom, typically by an order of magnitude or more, than the PW method. As a result, we compute the equation of state of LiH and show that the computed lattice constant and bulk modulus are in excellent agreement with reference PW results, while requiring an order of magnitude fewer degrees of freedom to obtain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Nonconforming Mesh Refinement for High-Order Finite Elements

In this study, we propose a general algorithm for nonconforming adaptive mesh refinement (AMR) of unstructured meshes in high-order finite element codes. Our focus is on $h$-refinement with a fixed polynomial order. The algorithm handles triangular, quadrilateral, hexahedral, and prismatic meshes of arbitrarily high-order curvature, for any order finite element space in the de Rham sequence. We present a flexible data structure for meshes with hanging nodes and a general procedure to construct the conforming interpolation operator, both in serial and in parallel. The algorithm and data structure allow anisotropic refinement of tensor product elements in two and three dimensions and support unlimited refinement ratios of adjacent elements. We report numerical experiments verifying the correctness of the algorithms and perform a parallel scaling study to show that we can adapt meshes containing billions of elements and run efficiently on 393,000 parallel tasks. Finally, we illustrate the integration of dynamic AMR into a high-order Lagrangian hydrodynamics solver.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2018

Finite Element Simulations of Two Vertical Drop Tests of F-28 Fuselage Sections

In March 2017, a vertical drop test of a forward fuselage section of a Fokker F-28 MK4000 aircraft was conducted as part of a joint NASA/FAA project to investigate the performance of transport aircraft under realistic crash conditions. In June 2017, a vertical drop test was conducted of a wing-box fuselage section of the same aircraft. Both sections were configured with two rows of aircraft seats, in a triple-double configuration. A total of ten Anthropomorphic Test Devices (ATDs) were secured in seats using standard lap belt restraints. The forward fuselage section was also configured with luggage in the cargo hold. Both sections were outfitted with two hat racks, each with added ballast mass. The drop tests were performed at the Landing and Impact Research facility located at NASA Langley Research Center in Hampton, Virginia. The measured impact velocity for the forward fuselage section was 346.8-in/s onto soil. The wing-box section was dropped with a downward facing pitch angle onto a sloping soil surface in order to create an induced forward acceleration in the airframe. The vertical impact velocity of the wing-box section was 349.2-in/s. A second objective of this project was to assess the capabilities of finite element simulations to predict the test responses. Finite element models of both fuselage sections were developed for execution in LS-DYNA(Registered Trademark), a commercial explicit nonlinear transient dynamic code. The models contained accurate representations of the airframe structure, the hat racks and hat rack masses, the floor and seat tracks, the luggage in the cargo hold for the forward section, and the detailed under-floor structure in the wing-box section. Initially, concentrated masses were used to represent the inertial properties of the seats, restraints, and ATD occupants. However, later simulations were performed that included finite element representations of the seats, restraints, and ATD occupants. These models were developed to more accurately replicate the seat loading of the floor and to enable prediction of occupant impact responses. Models were executed to generate analytical predictions of airframe responses, which were compared with test data to validate the model. Comparisons of predicted and experimental structural deformation and failures were made. Finally, predicted and experimental soil deformation and crater depths were also compared for both drop test configurations.

Jackson, Karen E.↗

DOE OSTI2022

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1985

Application of the finite element method in the calculation of transmission loss of flat and curved panels

This investigation represents an extension of a study of Roussos (1985) who considered the noise transmission loss of a rectangular plate in an infinite baffle. Roussos, who employed an analytical formulation, considered an unstiffened plate. While it is difficult to consider stiffeners by means of analytical methods, the difficulties can be avoided by employing a finite element procedure. For this reason, the present study is concerned with the implementation of a finite element method. The representation of the panel transmission loss is discussed, and the determination of the panel motion by means of the finite element technique is described, taking into account an isotropic flat panel, the exciting force, an eigenvalue problem, the radiation pressure, a plate element, and a cylindrical shell element. Numerical results are considered for a flat panel, a curved panel, and a stiffened flat panel.

Koval, L. R.↗

DOE OSTI2018

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1977

Finite element stress analysis of the human left ventricle whose irregular shape is developed from single plane cineangiocardiogram

The three-dimensional left ventricular chamber geometrical model is developed from single plane cineangiocardiogram. This left ventricular model is loaded by an internal pressure monitored by cardiac catheterization. The resulting stresses in the left ventricular model chamber's wall are determined by computerized finite element procedure. For the discretization of this left ventricular model structure, a 20-node, isoparametric finite element is employed. The analysis and formulation of the computerised procedure is presented in the paper, along with the detailed algorithms and computer programs. The procedure is applied to determine the stresses in a left ventricle at an instant, during systole. Next, a portion (represented by a finite element) of this left ventricular chamber is simulated as being infarcted by making its active-state modulus value equal to its passive-state value; the neighbouring elements are shown to relieve the 'infarcted' element of stress by themselves taking on more stress.

Ghista, D. N.↗

DOE OSTI2016

Generalization of mixed multiscale finite element methods with applications

Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixed multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. The algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

DOE OSTI2015

Parallel domain decomposition method for finite element approximation of 3D steady state non‐Newtonian fluids

Summary We introduce a stabilized finite element method for the 3D non‐Newtonian Navier–Stokes equations and a parallel domain decomposition method for solving the sparse system of nonlinear equations arising from the discretization. Non‐Newtonian flow problems are, generally speaking, more challenging than Newtonian flows because the nonlinearities are not only in the convection term but also in the viscosity term, which depends on the shear rate. Many good iterative methods and preconditioning techniques that work well for the Newtonian flows do not work well for the non‐Newtonian flows. We employ a Galerkin/least squares finite element method, with stabilization parameters adjusted to count the non‐Newtonian effect, to discretize the equations, and the resulting highly nonlinear system of equations is solved by a Newton–Krylov–Schwarz algorithm. In this study, we apply the proposed method to some inelastic power‐law fluid flows through the eccentric annuli with inner cylinder rotation and investigate the robustness of the method with respect to some physical parameters, including the power‐law index and the Reynolds number ratios. We then report the superlinear speedup achieved by the domain decomposition algorithm on a computer with up to 512 processors. Copyright © 2015 John Wiley & Sons, Ltd.

Shiu, Wen‐Shin↗

DOE OSTI2018

Determination of Thermoelastic Properties of Carbon Nanotube/Epoxy Composites Using Finite Element Method

In the present study, three-dimensional finite element (FE) analysis has been performed for determination of thermoelastic properties of carbon nanotube (CNT)/epoxy composites considering both linear and nonlinear stress–strain relations of CNT. A square representative volume element has been considered for the FE analysis (FEA), and the effect of considering nonlinear stress–strain relation of CNT on the prediction of thermoelastic properties has been determined by comparing the results with the theoretical results. Results show that in order to obtain realistic values of thermoelastic properties, the nonlinear stress–strain relation of CNT needs to be used in the FEA. Not considering nonlinear stress–strain relation of CNTs in the FEA may lead to estimating erroneous values of thermoelastic properties of CNT composites which may have adverse effect in the design of structural components made of such composites subjected to thermomechanical loading.

36 MATERIALS SCIENCE↗

NASA NTRS1981

Adaptive Finite Element Technology in Integrated Design and Analysis

An assessment of the potential impact of adaptive finite element technology on the analysis part of the aircraft structural synthesis process is presented. The main conclusion is that adaptive application of the p-version of the finite element method based on indirect error estimation procedures results in substantial cost reduction and increased reliability of the computed data. Adaptivity based on direct a posteriori error estimation has the potential for additional savings.

Szabo, B. A.↗

NASA NTRS1987

Design of Software for Design of Finite Element for Structural Analysis

The concepts of software engineering which allow a user of the finite element method to describe a model, to collect and to check the model data in a data base as well as to form the matrices required for a finite element calculation are examined. Next the components of the model description are conceived including the mesh tree, the topology, the configuration, the kinematic boundary conditions, the data for each element, and the loads. The possibilities for description and review of the data are considered. The concept of the segments for the modularization of the programs follows the components of the model description. The significance of the mesh tree as a globular guiding structure will be understood in view of the principle of the unity of the model, the mesh tree, and the data base. The user-friendly aspects of the software system will be summarized: the principle of language communication, the data generators, error processing, and data security.

Helfrich, Reinhard↗

NASA NTRS1989

A finite element approach for large motion dynamic analysis of multibody structures in space

A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Chang, Che-Wei↗

DOE OSTI2016

Investigation of grain-scale microstructural variability in tantalum using crystal plasticity-finite element simulations

In this research, a crystal plasticity-finite element (CP-FE) model is used to investigate the effects of microstructural variability at a notch tip in tantalum single crystals and polycrystals. It is shown that at the macroscopic scale, the mechanical response of single crystals is sensitive to the crystallographic orientation while the response of polycrystals shows relatively small susceptibility to it. However, at the microscopic scale, the local stress and strain fields in the vicinity of the crack tip are completely determined by the local crystallographic orientation at the crack tip for both single and polycrystalline specimens with similar mechanical field distributions. Variability in the local metrics used (maximum von Mises stress and equivalent plastic strain at 3% deformation) for 100 different realizations of polycrystals fluctuates by up to a factor of 2–7 depending on the local crystallographic texture. Comparison with experimental data shows that the CP model captures variability in stress–strain response of polycrystals that can be attributed to the grain-scale microstructural variability. In conclusion, this work provides a convenient approach to investigate fluctuations in the mechanical behavior of polycrystalline materials induced by grain morphology and crystallographic orientations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗