Engineering Papers⌕ Search

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.

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

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Structural Analysis and Performance-Based Validation of a Composite Wing Spar

Electric-motor powered aircraft possess the ability to operate with efficient energy delivery, but lack the operational range of internal combustion engine powered aircraft. This range limitation requires the use of high aspect ratio, thin-chord wings to minimize aerodynamic drag losses, which results in highly loaded composite spar structures. High aspect ratio wings are required to increase mission durations for a NASA-developed experimental multi-rotor electric powered aircraft denoted as the Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) or “X-57”. This paper examines the structural performance of the composite main wing spars to validate spar strength using ply-based laminate finite element methods. Geometric scaling of a main spar test-section was initially proposed for proof-testing but sacrificed stability. Ply-based structures modeling with local structural features was implemented as a risk-reduction methodology. Ply-based modeling was selected to augment the conventional “building block” approach to reduce risk, and leverage a “performance-based” approval processes encouraged in Federal Aviation Administration (FAA) design guidance. Therefore, ply-based laminate modeling of the full-scale main spar and forward spar shear-web attachments were subsequently undertaken to determine load path complexity with predicted flight loads. Ply-based modeling included stress concentrations and interlaminate behavior at interface locations that can be obscured in traditional finite element sizing models. Analysis of the wing spar laminate ply-based models compared with bearing test coupon performance was used to reduce future wing assembly proof-testing burden and facilitate performance-based flight hardware safety for the X-57 experimental aircraft.

Moore, James B.↗
DOE OSTI2015

Interface profile optimization for planar stress wave attenuation in bi-layered plates

Stress waves scatter upon entering a new medium. This occurs due to the reflection and transmission of the waves, which depends on the impedance mismatch between the two materials and the angle of incidence. For a bi-layered structure with finite dimensions and constant impedance ratio, the scattering and intensity of the stress waves may be varied by changing the interface profile between the two layers. In this study, a methodology is proposed for optimizing the interface profile between the layers of a finite bi-layered plate for the objective of planar stress wave attenuation. The bi-layered plates are subjected at one end to highly impulsive loadings with various durations, and the geometry of the internal interface is optimized for the purpose of minimizing the amplitude of the maximum reaction force at the opposite fixed end. The optimization methodology is based on a genetic algorithm, which is coupled with a finite element method for analyzing the wave propagation behavior of the plates. Finally, it is observed that the interface profile and the amount of stress wave attenuation depend on the duration of the applied impulsive loading, with higher amounts of attenuation obtained when the wavelength associated with the impulsive load is small compared to the dimensions of the bi-layered plates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kernel methods for evolution of generalized parton distributions

Generalized parton distributions (GPDs) characterize the 3-dimensional structure of hadrons, combining information about their internal quark and gluon longitudinal momentum distributions and transverse position within the hadron. The dependence of GPDs on the factorization scale Q 2 allows one to connect hard exclusive processes involving GPDs at disparate energy and momentum scales, which is needed in global analyses of experimental data. Here, in this work, we explore how finite element methods can be used to construct fast and differentiable Q 2 evolution codes for GPDs in momentum space, which can be used in a machine learning framework. We show numerical benchmarks of the methods' accuracy, including a comparison to an existing evolution code from PARTONS/APFEL++, and provide a repository where the code can be accessed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to various forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.

Pan, Wenxiao↗
NASA NTRS1989

A finite difference treatment of Stokes-type flows: Preliminary report

The equations Laplacian operator omega = 0, (1.1a) and omega = Laplacian operator Chi, (1.1b) describe, in suitable units, 2-D Stokes flow of an incompressible fluid occupying a domain D in which omega is the vorticity and Chi is the stream function. The flow is uniquely determined by specifying the velocity on the boundary B of D, a condition which leads to specifying the stream function Chi and its normal derivative Chi sub n on B. A mathematically similar problem arises in describing the equilibrium of a flat plate in structural mechanics where a related 1-D problem by finite difference or finite element methods is to introduce effective methods for imposing the boundary conditions through which (1.1a) is coupled to (1.1b). These models thus provide a simple starting point for examining the general treatment of boundary conditions for more general time dependent Navier-Stokes incompressible flows. For the purpose of discussion it is assumed that D is a square domain. A standard finite difference method to solve (1.1) is to introduce a uniform grid and then use standard five point finite difference operators to express each equation in (1.1). At any point on the boundary B a value of Chi is specified by the boundary conditions but a value of omega at the same boundary mesh point will also be required to complete the computation. Methods are discussed which overcome the difficulty in solving these problems.

Rose, M. E.↗
DOE OSTI2024

Adaptive Data-Driven Deep-Learning Surrogate Model for Frontal Polymerization in Dicyclopentadiene

Frontal polymerization (FP) is a self-sustaining curing process that enables rapid and energy-efficient manufacturing of thermoset polymers and composites. Computational methods conventionally used to simulate the FP process are time-consuming, and repeating simulations are required for sensitivity analysis, uncertainty quantification, or optimization of the manufacturing process. Here, in this work, we develop an adaptive surrogate deep-learning model for FP of dicyclopentadiene (DCPD), which predicts the evolution of temperature and degree of cure orders of magnitude faster than the finite-element method (FEM). The adaptive algorithm provides a strategy to select training samples efficiently and save computational costs by reducing the redundancy of FEM-based training samples. The adaptive algorithm calculates the residual error of the FP governing equations using automatic differentiation of the deep neural network. A probability density function expressed in terms of the residual error is used to select training samples from the Sobol sequence space. The temperature and degree of cure evolution of each training sample are obtained by a 2D FEM simulation. The adaptive method is more efficient and has a better prediction accuracy than the random sampling method. With the well-trained surrogate neural network, the FP characteristics (front speed, shape, and temperature) can be extracted quickly from the predicted temperature and degree-of-cure fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

SIAM Conference on Parallel Processing for Scientific Computing, 4th, Chicago, IL, Dec. 11-13, 1989, Proceedings

Attention is given to such topics as an evaluation of block algorithm variants in LAPACK and presents a large-grain parallel sparse system solver, a multiprocessor method for the solution of the generalized Eigenvalue problem on an interval, and a parallel QR algorithm for iterative subspace methods on the CM2. A discussion of numerical methods includes the topics of asynchronous numerical solutions of PDEs on parallel computers, parallel homotopy curve tracking on a hypercube, and solving Navier-Stokes equations on the Cedar Multi-Cluster system. A section on differential equations includes a discussion of a six-color procedure for the parallel solution of elliptic systems using the finite quadtree structure, data parallel algorithms for the finite element method, and domain decomposition methods in aerodynamics. Topics dealing with massively parallel computing include hypercube vs. 2-dimensional meshes and massively parallel computation of conservation laws. Performance and tools are also discussed.

Dongarra, Jack↗
NASA NTRS1989

Procedures for shape optimization of gas turbine disks

Two procedures, the feasible direction method and sequential linear programming, for shape optimization of gas turbine disks are presented. The objective of these procedures is to obtain optimal designs of turbine disks with geometric and stress constraints. The coordinates of the selected points on the disk contours are used as the design variables. Structural weight, stress and their derivatives with respect to the design variables are calculated by an efficient finite element method for design senitivity analysis. Numerical examples of the optimal designs of a disk subjected to thermo-mechanical loadings are presented to illustrate and compare the effectiveness of these two procedures.

Cheu, Tsu-Chien↗
DOE OSTI2017

Overset meshing coupled with hybridizable discontinuous Galerkin finite elements

We introduce the use of hybridizable discontinuous Galerkin (HDG) finite element methods on overlapping (overset) meshes. Overset mesh methods are advantageous for solving problems on complex geometrical domains. We also combine geometric flexibility of overset methods with the advantages of HDG methods: arbitrarily high-order accuracy, reduced size of the global discrete problem, and the ability to solve elliptic, parabolic, and/or hyperbolic problems with a unified form of discretization. This approach to developing the ‘overset HDG’ method is to couple the global solution from one mesh to the local solution on the overset mesh. We present numerical examples for steady convection–diffusion and static elasticity problems. The examples demonstrate optimal order convergence in all primal fields for an arbitrary amount of overlap of the underlying meshes.

42 ENGINEERING↗
NASA NTRS2004

Integrated Force Method Solution to Indeterminate Structural Mechanics Problems

Strength of materials problems have been classified into determinate and indeterminate problems. Determinate analysis primarily based on the equilibrium concept is well understood. Solutions of indeterminate problems required additional compatibility conditions, and its comprehension was not exclusive. A solution to indeterminate problem is generated by manipulating the equilibrium concept, either by rewriting in the displacement variables or through the cutting and closing gap technique of the redundant force method. Compatibility improvisation has made analysis cumbersome. The authors have researched and understood the compatibility theory. Solutions can be generated with equal emphasis on the equilibrium and compatibility concepts. This technique is called the Integrated Force Method (IFM). Forces are the primary unknowns of IFM. Displacements are back-calculated from forces. IFM equations are manipulated to obtain the Dual Integrated Force Method (IFMD). Displacement is the primary variable of IFMD and force is back-calculated. The subject is introduced through response variables: force, deformation, displacement; and underlying concepts: equilibrium equation, force deformation relation, deformation displacement relation, and compatibility condition. Mechanical load, temperature variation, and support settling are equally emphasized. The basic theory is discussed. A set of examples illustrate the new concepts. IFM and IFMD based finite element methods are introduced for simple problems.

Patnaik, Surya N.↗
NASA NTRS2002

Calculation of Scattering Amplitude Without Partial Analysis: Inclusion of Exchange - II

There was a method for calculating the whole scattering amplitude, f(Omega(sub k)), directly. The idea was to calculate the complete wave function Psi numerically, and use it in an integral expression for f, which can be reduced to a 2 dimensional quadrature. The original application was for e-H scattering without exchange. There the Schrodinger reduces a 2-d partial differential equation (pde), which was solved using the finite element method (FEM). Here we extend the method to the exchange approximation. The S.E. can be reduced to a pair of coupled pde's, which are again solved by the FEM. The formal expression for f(Omega(sub k)) consists two integrals, f+/- = f(sub d) +/- f(sub e); f(sub d) is formally the same integral as the no-exchange f. We have also succeeded in reducing f(sub e) to a 2-d integral. Results will be presented at the meeting.

Temkin, Aaron↗
DOE OSTI2022

Computational Methods to Characterize Panel Loading Conditions for Accelerated Testing

Panel cracking and degradation due to wind loading are known to have a detrimental effect on power output. Recreating these damaging conditions in a controlled experimental setting requires an understanding of the loads being generated at the panel surface as a function of both wind speed, wind direction, and panel orientation. To better understand these relationships, a computational fluid dynamics (CFD) simulation package was constructed using an open-source Python library for solving partial differential equations with the finite element method. This CFD package allows the simulated wind speed and panel orientation of a multi-panel array to be easily changed and provides the capability to measure the traction forces at discrete points along the sun-facing and ground-facing surface of an interior panel residing in the wake of one or more upstream panels. A parameter exploration was performed in which the panel angle was varied in even increments from -70 deg to 70 deg and the wind speed was varied from 2 to 30 m/s. During post processing, the measured traction along the panel surface was averaged spatially and interpreted as a time-varying signal, where further processing of this signal yielded the root mean square amplitude and a characteristic frequency associated with the loading. This study found that higher wind speeds are generally associated with increased amplitude of loading and that the panel orientation angle can significantly exacerbate or mitigate this loading. These outcomes are presented along with current work on higher-fidelity verification simulations and recommendations for performing accelerated experimental testing.

loading↗
DOE OSTI2016

A three-dimensional meso-macroscopic model for Li-Ion intercalation batteries

Through this study, we present a three-dimensional computational formulation for electrode-electrolyte-electrode system of Li-Ion batteries. The physical consistency between electrical, thermal and chemical equations is enforced at each time increment by driving the residual of the resulting coupled system of nonlinear equations to zero. The formulation utilizes a rigorous volume averaging approach typical of multiphase formulations used in other fields and recently extended to modeling of supercapacitors [1]. Unlike existing battery modeling methods which use segregated solution of conservation equations and idealized geometries, our unified approach can model arbitrary battery and electrode configurations. The consistency of multi-physics solution also allows for consideration of a wide array of initial conditions and load cases. The formulation accounts for spatio-temporal variations of material and state properties such as electrode/void volume fractions and anisotropic conductivities. The governing differential equations are discretized using the finite element method and solved using a nonlinearly consistent approach that provides robust stability and convergence. The new formulation was validated for standard Li-ion cells and compared against experiments. Finally, its scope and ability to capture spatio-temporal variations of potential and lithium distribution is demonstrated on a prototypical three-dimensional electrode problem.

25 ENERGY STORAGE↗
NASA NTRS1980

Application of finite-element-techniques to the interaction of conduction and radiation in an absorbing, scattering and emitting medium

In this paper, the authors demonstrate that a Galerkin finite element method of analysis, utilizing isoparametric elements, offers a viable means of solving continuum thermal radiation problems with conduction in a participating medium. The participating medium was considered to be a gray radiation medium exhibiting isotropic absorption, emission, and scattering characteristics and optical properties that are independent of temperature. The medium was considered to be bounded by infinite parallel opaque, gray surfaces with diffuse emission and reflection characteristics. In solving this problem, a finite element formulation was developed to describe a system in radiative equilibrium. Then the results of this first analysis were linked with a second finite element model which incorporated conduction into the analysis. The results of this study were found to be in good agreement with existing published data. The model offers the following advantageous features: geometric generality, a computational algorithm which is 'convenient' and 'computable', and a functional basis for extension of the radiation model to higher order approximation.

Wu, S. T.↗
NASA NTRS1985

Contributions to the finite element solution of the fan noise radiation problem

The radiation of fan generated noise to the far field from a nacelle of realistic geometry is investigated using the finite element method. Several innovations have been introduced to minimize the computational requirements and create a highly efficient numerical scheme. The innovations include: (1) formulation of the problem in terms of velocity potential and density in such a way that no inlet mean flow velocity derivatives are required in the field equations, (2) the use of 'wave envelope' elements in an outer region permitting a grid much coarser than would be used for conventional finite elements, (3) the use of a mesh which deforms with an increase of forward flight speed so that mesh lines are always lines of constant phase and rays for a point source, permitting the use of wave envelope elements and simple boundary conditions for any case of forward velocity, (4) an efficient scheme for introducing the noise source via modal amplitude coefficients, and (5) the use of a frontal solution technique which for physically realistic problems drastically reduces the active storage requirements. The finite element scheme is outlined, as are the specific details of the innovations. Results are given for cases where comparable experimental data are available.

Eversman, W.↗
NASA NTRS1983

Laser microspeckle technique in displacement measurement near a crack tip

The laser speckle method has been found quite useful in obtaining in-plane displacement measurements. It is especially useful if the displacement field in a small region can be effectively determined. By obtaining directly the speckle patterns at higher magnifications, a better distinction of the displacements in the vicinity of a crack tip is possible. In this short report some results are obtained and compared with those calculated from linear fracture mechanics and the finite element method for an aluminum crack.

Tang, Z. Q.↗
NASA NTRS1995

Analysis of three-dimensional-cavity-backed aperture antennas using a Combined Finite Element Method/Method of Moments/Geometrical Theory of Diffraction technique

A combined finite element method (FEM) and method of moments (MoM) technique is presented to analyze the radiation characteristics of a cavity-fed aperture in three dimensions. Generalized feed modeling has been done using the modal expansion of fields in the feed structure. Numerical results for some feeding structures such as a rectangular waveguide, circular waveguide, and coaxial line are presented. The method also uses the geometrical theory of diffraction (GTD) to predict the effect of a finite ground plane on radiation characteristics. Input admittance calculations for open radiating structures such as a rectangular waveguide, a circular waveguide, and a coaxial line are shown. Numerical data for a coaxial-fed cavity with finite ground plane are verified with experimental data.

Reddy, C. J.↗
DOE OSTI2017

Pressure Self-focusing Effect and Novel Methods for Increasing the Maximum Pressure in Traditional and Rotational Diamond Anvil Cells

The main principles of producing a region near the center of a sample, compressed in a diamond anvil cell (DAC), with a very high pressure gradient and, consequently, with high pressure are predicted theoretically. The revealed phenomenon of generating extremely high pressure gradient is called the pressure self-focusing effect. Initial analytical predictions utilized generalization of a simplified equilibrium equation. Then, the results are refined using our recent advanced model for elastoplastic material under high pressures in finite element method (FEM) simulations. The main points in producing the pressure self-focusing effect are to use beveled anvils and reach a very thin sample thickness at the center. We find that the superposition of torsion in a rotational DAC (RDAC) offers drastic enhancement of the pressure self-focusing effect and allows one to reach the same pressure under a much lower force and deformation of anvils.

42 ENGINEERING↗