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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Investigation on Vibration Response of Aluminum Foam Beams Using Speckle Interferometry

This paper firstly presents a brief review on the representations of speckle fringe pattern in vibration measurement using time-averaged electronic speckle pattern interference (ESPI) method. Based on laser phase noise, a new representation of speckle fringe intensity generated by real-time subtraction ESPI method for vibrating measurement is then proposed. In the experimental performance, the real-time subtraction ESPI method is employed to inspect the vibrating response of a closed-cell aluminum foam cantilever beam. In contract with the results of the finite elemental method simulation, the vibration mode shapes obtained by the ESPI method are well agreed with numerical prediction. The optical vibration analysis is also carried out to determine the effective Young’s modulus of aluminum foam, and the results verify the validation of the ESPI method for investigation on mechanical properties of metal foam materials.

36 MATERIALS SCIENCE↗
NASA NTRS1995

The origin of spurious solutions in computational electromagnetics

The origin of spurious solutions in computational electromagnetics, which violate the divergence equations, is deeply rooted in a misconception about the first-order Maxwell's equations and in an incorrect derivation and use of the curl-curl equations. The divergence equations must be always included in the first-order Maxwell's equations to maintain the ellipticity of the system in the space domain and to guarantee the uniqueness of the solution and/or the accuracy of the numerical solutions. The div-curl method and the least-squares method provide rigorous derivation of the equivalent second-order Maxwell's equations and their boundary conditions. The node-based least-squares finite element method (LSFEM) is recommended for solving the first-order full Maxwell equations directly. Examples of the numerical solutions by LSFEM for time-harmonic problems are given to demonstrate that the LSFEM is free of spurious solutions.

Jiang, Bo-Nan↗
NASA NTRS2001

Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) for a 3-D Flexible Wing

The formulation and implementation of an optimization method called Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) are extended from single discipline analysis (aerodynamics only) to multidisciplinary analysis - in this case, static aero-structural analysis - and applied to a simple 3-D wing problem. The method aims to reduce the computational expense incurred in performing shape optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, Finite Element Method (FEM) structural analysis and sensitivity analysis tools. Results for this small problem show that the method reaches the same local optimum as conventional optimization. However, unlike its application to the win,, (single discipline analysis), the method. as I implemented here, may not show significant reduction in the computational cost. Similar reductions were seen in the two-design-variable (DV) problem results but not in the 8-DV results given here.

Gumbert, Clyde R.↗
NASA NTRS2001

Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) for a 3-D Flexible Wing

The formulation and implementation of an optimization method called Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) are extended from single discipline analysis (aerodynamics only) to multidisciplinary analysis - in this case, static aero-structural analysis and applied to a simple 3-D wing problem. The method aims to reduce the computational expense incurred in performing shape optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, Finite Element Method (FEM) structural analysis and sensitivity analysis tools. Results for this small problem show that the method reaches the same local optimum as conventional optimization. However, unlike its application to the rigid wing (single discipline analysis), the method, as implemented here, may not show significant reduction in the computational cost. Similar reductions were seen in the two-design-variable (DV) problem results but not in the 8-DV results given here.

Gumbert, Clyde R.↗
NASA NTRS2011

Effect of Boundary Conditions on the Axial Compression Buckling of Homogeneous Orthotropic Composite Cylinders in the Long Column Range

Buckling loads for long isotropic and laminated cylinders are calculated based on Euler, Fluegge and Donnell's equations. Results from these methods are presented using simple parameters useful for fundamental design work. Buckling loads for two types of simply supported boundary conditions are calculated using finite element methods for comparison to select cases of the closed form solution. Results indicate that relying on Donnell theory can result in an over-prediction of buckling loads by as much as 40% in isotropic materials.

Mikulas, Martin M., Jr.↗
DOE OSTI2016

Hybrid PN-SN Calculations with SAAF for the Multiscale Transport Capability in Rattlesnake

Two interface conditions, the Lagrange multiplier method and the upwinding method, for hybrid \pn-\sn calculations is proposed for the self-adjoint angular flux (SAAF) formulation of the transport equation using the continuous finite element method (FEM) for spatial discretization. These interface conditions are implemented in Rattlesnake, the radiation transport application built on MOOSE, for the on-going multiscale transport simulation effort at INL. For smoothing the solution at the interface for the Lagrange multiplier method, a method based on \sn Lagrange interpolation on the sphere is proposed. Numerical results indicate that the interface conditions give the expected convergence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
NASA NTRS1986

Preconditioned conjugate residual methods for the solution of spectral equations

Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.

Wong, Y. S.↗
NASA NTRS1987

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semianalytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved from a rigorous mathematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggest that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both of the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.

Chaudhuri, Reaz A.↗
NASA NTRS1999

Experimentally Determined Crack Location and Mixed Mode Stress Intensity Factors

This presentation first examines five strain field equations: the two modes of the Westergaard equations, and a mixed mode; the generalized coordinate system, and the significance testing of coefficients. It also discusses the weighted multiple linear least squares method, which shows why it is difficult to find the crack tip. The presentation then discusses the experiments to examine stress intensity. Lastly, comparison of finite element method results and experimental results is reviewed.

Zachary, Loren W.↗
NASA NTRS1989

Global functions in global-local finite-element analysis of localized stresses in prismatic structures

An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.

Dong, Stanley B.↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2014

An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains

Many practical applications require the analysis of elastic wave propagation in a homogeneous isotropic media in an unbounded domain. One widely used approach for truncating the infinite domain is the so called method of perfectly matched layers (PMLs). Most existing PML formulations are developed for finite difference methods based on the first-order velocity-stress form of the elasticity equations, and they are not straight-forward to implement using standard finite element methods (FEMs) on unstructured meshes. Some of the problems with these formulations include the application of boundary conditions in half-space problems and in the treatment of edges and/or corners for time-domain problems. Several PML formulations, which do work with FEMs have been proposed, although most of them still have some of these problems and/or they require a large number of auxiliary nodal history/memory variables. Here in this work, we develop a new PML formulation for time-domain elastodynamics on a spherical domain, which reduces to a two-dimensional formulation under the assumption of axisymmetry. Our formulation is well-suited for implementation using FEMs, where it requires lower memory than existing formulations, and it allows for natural application of boundary conditions. We solve example problems on two-dimensional and three-dimensional domains using a high-order discontinuous Galerkin (DG) discretization on unstructured meshes and explicit time-stepping. We also study an approach for stabilization of the discrete equations, and we show several practical applications for quality factor predictions of micromechanical resonators along with verifying the accuracy and versatility of our formulation.

42 ENGINEERING↗
DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2017

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

In this work, a phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equation and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
NASA NTRS1985

Shape design sensitivity analysis using domain information

A numerical method for obtaining accurate shape design sensitivity information for built-up structures is developed and demonstrated through analysis of examples. The basic character of the finite element method, which gives more accurate domain information than boundary information, is utilized for shape design sensitivity improvement. A domain approach for shape design sensitivity analysis of built-up structures is derived using the material derivative idea of structural mechanics and the adjoint variable method of design sensitivity analysis. Velocity elements and B-spline curves are introduced to alleviate difficulties in generating domain velocity fields. The regularity requirements of the design velocity field are studied.

Seong, Hwal-Gyeong↗
NASA NTRS1996

Radiation Characteristics of Cavity Backed Aperture Antennas in Finite Ground Plane Using the Hybrid FEM/MoM Technique and Geometrical Theory of Diffraction

A technique using hybrid Finite Element Method (FEM)/Method of Moments (MoM), and Geometrical Theory of Diffraction (GTD) is presented to analyze the radiation characteristics of cavity fed aperture antennas in a finite ground plane. The cavity which excites the aperture is assumed to be fed by a cylindrical transmission line. The electromagnetic (EM) fields inside the cavity are obtained using FEM. The EM fields and their normal derivatives required for FEM solution are obtained using (1) the modal expansion in the feed region and (2) the MoM for the radiating aperture region(assuming an infinite ground plane). The finiteness of the ground plane is taken into account using GTD. The input admittance of open ended circular, rectangular, and coaxial line radiating into free space through an infinite ground plane are computed and compared with earlier published results. Radiation characteristics of a coaxial cavity fed circular aperture in a finite rectangular ground plane are verified with experimental results.

Reddy, C. J.↗