Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Method (FEM)”

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 181 records · Page 10

Application of Finite Element Method to Analyze Inflatable Waveguide Structures

A Finite Element Method (FEM) is presented to determine propagation characteristics of deformed inflatable rectangular waveguide. Various deformations that might be present in an inflatable waveguide are analyzed using the FEM. The FEM procedure and the code developed here are so general that they can be used for any other deformations that are not considered in this report. The code is validated by applying the present code to rectangular waveguide without any deformations and comparing the numerical results with earlier published results.

Deshpande, M. D.↗

Application of Model Based Parameter Estimation for Fast Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using Hybrid FEM/MoM Technique

Model Based Parameter Estimation (MBPE) is presented in conjunction with the hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique for fast computation of the input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FENI/MoM technique is used to form an integro-partial- differential equation to compute the electric field distribution of a cavity-backed aperture antenna. In MBPE, the electric field is expanded in a rational function of two polynomials. The coefficients of the rational function are obtained using the frequency derivatives of the integro-partial-differential equation formed by the hybrid FEM/ MoM technique. Using the rational function approximation, the electric field is obtained over a frequency range. Using the electric field at different frequencies, the input characteristics of the antenna are obtained over a wide frequency range. Numerical results for an open coaxial line, probe-fed coaxial cavity and cavity-backed microstrip patch antennas are presented. Good agreement between MBPE and the solutions over individual frequencies is observed.

Reddy C. J.↗

User's Manual for FEMOM3DR

FEMoM3DR is a computer code written in FORTRAN 77 to compute radiation characteristics of antennas on 3D body using combined Finite Element Method (FEM)/Method of Moments (MoM) technique. The code is written to handle different feeding structures like coaxial line, rectangular waveguide, and circular waveguide. This code uses the tetrahedral elements, with vector edge basis functions for FEM and triangular elements with roof-top basis functions for MoM. By virtue of FEM, this code can handle any arbitrary shaped three dimensional bodies with inhomogeneous lossy materials; and due to MoM the computational domain can be terminated in any arbitrary shape. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computers on which the code is intended to run.

Reddy, C. J.↗

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗

A 3-D Magnetic Analysis of a Linear Alternator For a Stirling Power System

The NASA Glenn Research Center and the Department of Energy (DOE) are developing advanced radioisotope Stirling convertors, under contract with Stirling Technology Company (STC), for space applications. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. This paper presents a 3-D finite element method (FEM) approach for evaluating Stirling convertor linear alternators. Preliminary correlations with open-circuit voltage measurements provide an encouraging level of confidence in the model. Spatial plots of magnetic field strength (H) are presented in the region of the exciting permanent magnets. These plots identify regions of high H, where at elevated temperature and under electrical load, the potential to alter the magnetic moment of the magnets exists. This implies the need for further testing and analysis.

Geng, Steven M.↗

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.↗

Direct Calculation of the Scattering Amplitude Without Partial Wave Analysis

Two new developments in scattering theory are reported. We show, in a practical way, how one can calculate the full scattering amplitude without invoking a partial wave expansion. First, the integral expression for the scattering amplitude f(theta) is simplified by an analytic integration over the azimuthal angle. Second, the full scattering wavefunction which appears in the integral expression for f(theta) is obtained by solving the Schrodinger equation with the finite element method (FEM). As an example, we calculate electron scattering from the Hartree potential. With minimal computational effort, we obtain accurate and stable results for the scattering amplitude.

Shertzer, J.↗

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.↗

Simultaneous Aerodynamic and Structural Design Optimization (SASDO) for a 3-D Wing

The formulation and implementation of an optimization method called Simultaneous Aerodynamic and Structural Design Optimization (SASDO) is shown as an extension of the Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) method. It is extended by the inclusion of structure element sizing parameters as design variables and Finite Element Method (FEM) analysis responses as constraints. The method aims to reduce the computational expense. incurred in performing shape and sizing optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, FEM structural analysis and sensitivity analysis tools. SASDO is applied to a simple. isolated, 3-D wing in inviscid flow. Results show that the method finds the saine local optimum as a conventional optimization method with some reduction in the computational cost and without significant modifications; to the analysis tools.

Gumbert, Clyde R.↗

Structural Qualification Testing of the WindSat Payload Using Sine Bursts Near Structural Resonance

Sine burst tests are often used for structural qualification of space flight hardware. In most instances, the driving frequency of the shaker is specified far below the structure's first resonant mode, such that the entire test article sees uniform acceleration. For large structures, this limits qualification testing to lower parts of the structure, or else it over-tests the lower structure to achieve qualification of the upper structure. The WindSat payload, a 10.5 foot tall graphite/epoxy, titanium, and aluminum radiometer, experiences accelerations at the six foot diameter reflector nearly four times that at the spacecraft interface. Due to size of the payload, the number of bonded joints, and the lightweight reflector support structure design and construction, using static pull testing to qualify all of the bonded joints in the upper structure would result in large, expensive, and extensive test fixturing. Sine burst testing near the first two structural resonant modes was performed on the WindSat payload to achieve the correct load factor distribution up the stack for structural qualification. In this presentation, how finite element method (FEM) sine burst predictions were used in conjunction with low level random and sine burst tests to achieve correct qualification test load factor distribution on the WindSat payload is discussed. Also presented is the risk mitigation approach for using the uncorrelated FEM in this procedure.

Pontius, Jim↗

A 3-D Magnetic Analysis of a Stirling Convertor Linear Alternator Under Load

The NASA Glenn Research Center (GRC), the Department of Energy (DOE), and the Stirling Technology Company (STC) are developing Stirling convertors for Stirling Radioisotope Power Systems (SRPS) to provide electrical power for future NASA deep space missions. STC is developing the 55-We Technology Demonstration Convertor (TDC) under contract to DOE. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. This paper presents a 3-dimensional finite element method (FEM) approach for evaluating Stirling convertor linear alternators. The model extends a magnetostatic analysis previously reported at the 35th Intersociety Energy Conversion Engineering Conference (IECEC) to include the effects of the load current. STC's 55-We linear alternator design was selected to validate the model. Spatial plots of magnetic field strength (H) are presented in the region of the exciting permanent magnets. The margin for permanent magnet demagnetization is calculated at the expected magnet operating temperature for the near earth environment and for various average magnet temperatures. These thermal conditions were selected to represent a worst-case condition for the planned deep space missions. This paper presents plots that identify regions of high H where the potential to alter the magnetic moment of the magnets exists.

Geng, Steven M.↗

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↗

Space Radiation Transport Methods Development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be utilized in the final design as verification of the deterministic method optimized design.

Wilson, J. W.↗

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.↗

Thermal-Structural Optimization of Integrated Cryogenic Propellant Tank Concepts for a Reusable Launch Vehicle

A next generation reusable launch vehicle (RLV) will require thermally efficient and light-weight cryogenic propellant tank structures. Since these tanks will be weight-critical, analytical tools must be developed to aid in sizing the thickness of insulation layers and structural geometry for optimal performance. Finite element method (FEM) models of the tank and insulation layers were created to analyze the thermal performance of the cryogenic insulation layer and thermal protection system (TPS) of the tanks. The thermal conditions of ground-hold and re-entry/soak-through for a typical RLV mission were used in the thermal sizing study. A general-purpose nonlinear FEM analysis code, capable of using temperature and pressure dependent material properties, was used as the thermal analysis code. Mechanical loads from ground handling and proof-pressure testing were used to size the structural geometry of an aluminum cryogenic tank wall. Nonlinear deterministic optimization and reliability optimization techniques were the analytical tools used to size the geometry of the isogrid stiffeners and thickness of the skin. The results from the sizing study indicate that a commercial FEM code can be used for thermal analyses to size the insulation thicknesses where the temperature and pressure were varied. The results from the structural sizing study show that using combined deterministic and reliability optimization techniques can obtain alternate and lighter designs than the designs obtained from deterministic optimization methods alone.

Johnson, Theodore F.↗

Development of a Fatigue Crack Growth Coupon for Highly Plastic Stress Conditions

This paper presents an analytical approach used to develop a novel fatigue crack growth coupon for a highly plastic 3-D stress field condition. The flight hardware investigated in this paper is a large separation bolt that fractures using pyrotechnics at the appointed time during the flight sequence. The separation bolt has a deep notch that produces a severe stress concentration and a large plastic zone when highly loaded. For this geometry, linear-elastic fracture mechanics (LEFM) techniques are not valid due to the large nonlinear stress field. Unfortunately, industry codes that are generally available for fracture mechanics analysis and fatigue crack growth (e.g. NASGRO (11) are limited to LEFM and are available for only a limited number of geometries. The results of LEFM based codes are questionable when used on geometries with significant plasticity. Therefore elastic-plastic fracture mechanics (EPFM) techniques using the finite element method (FEM) were used to analyze the bolt and test coupons. scale flight hardware is very costly in t e r n of assets, laboratory resources, and schedule. Therefore to alleviate some of these problems, a series of novel test coupons were developed to simulate the elastic-plastic stress field present in the bolt.

Allen, Phillip A.↗

Fast Computation of Frequency Response of Cavity-Backed Apertures Using MBPE in Conjunction with Hybrid FEM/MoM Technique

The hybrid Finite Element Method(FEM)/Method of Moments(MoM) technique has become popular over the last few years due to its flexibility to handle arbitrarily shaped objects with complex materials. One of the disadvantages of this technique, however, is the computational cost involved in obtaining solutions over a frequency range as computations are repeated for each frequency. In this paper, the application of Model Based Parameter Estimation (MBPE) method[1] with the hybrid FEM/MoM technique is presented for fast computation of frequency response of cavity-backed apertures[2,3]. In MBPE, the electric field is expanded in a rational function of two polynomials. The coefficients of the rational function are obtained using the frequency-derivatives of the integro-differential equation formed by the hybrid FEM/MoM technique. Using the rational function approximation, the electric field is calculated at different frequencies from which the frequency response is obtained.

Reddy, C. J.↗

Static Aeroelastic Analysis of Transonic Wind Tunnel Models Using Finite Element Methods

A computational method for accurately predicting the static aeroelastic deformations of typical transonic transport wind tunnel models is described. The method utilizes a finite element method (FEM) for predicting the deformations. Extensive calibration/validation of this method was carried out using a novel wind-off wind tunnel model static loading experiment and wind-on optical wing twist measurements obtained during a recent wind tunnel test in the National Transonic Facility (NTF) at NASA LaRC. Further validations were carried out using a Navier-Stokes computational fluid dynamics (CFD) flow solver to calculate wing pressure distributions about several aeroelastically deformed wings and comparing these predictions with NTF experimental data. Results from this aeroelastic deformation method are in good overall agreement with experimentally measured values. Including the predicted deformations significantly improves the correlation between CFD predicted and experimentally measured wing & pressures.

Hooker, John R.↗