Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference 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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗
NASA NTRS1999

Open-Mode Debonding Analysis of Curved Sandwich Panels Subjected to Heating and Cryogenic Cooling on Opposite Faces

Increasing use of curved sandwich panels as aerospace structure components makes it vital to fully understand their thermostructural behavior and identify key factors affecting the open-mode debonding failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich panels with different radii of curvature. The curved sandwich panels are either simply supported or clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the concave side. The finite-element method was used to study the effects of panel curvature and boundary condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich panels. The critical stress point, where potential debonding failure could initiate, was found to be at the midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and the associated sandwich core depth stretching.

Ko, William L.↗
NASA NTRS2018

Analysis of the Influence of Layer Shifting on the Elastic Response and Damage Nucleation and Growth in Woven Composite Laminates

The influence of relative layer shifting on the elastic and damage response of plain weave composite laminates is analyzed using a continuum damage mechanics approach in combination with the finite element method. First, the homogenized properties of the woven composite as a function of the number of layers and of layer shifting are presented. Next, the damage development in various shifting configurations is studied using different damage constitutive models for the matrix and the fiber bundles. It is shown that the impact of layer shifting on both the elastic response and the nonlinear damage response is significant. Most notably, the model captures changes in the damage mechanisms within the woven composite that occur due to layer shifting, resulting in stiffer, more brittle behavior, which has been shown experimentally in the literature. Model results in the linear and nonlinear regimes are shown to be consistent with both an independent analytical model and reported experiments.

layer shifting↗
NASA NTRS1986

Numerical prediction of cold turbulent flow in combustor configurations with different centerbody flame holders

A numerical study of cold turbulent flow in a combustor containing a centerbody flame holder is presented. The axisymmetric Navier-Stokes equations incorporating a k-epsilon turbulence model were solved in a nonorthogonal curvilinear coordinate system. The finite volume method, applied to a staggered grid system, was used to discretize the differential equations. These finite differenced equations were then solved iteratively utilizing the SIMPLE algorithm described by Patankar (1980). Solutions were obtained for two combustor duct geometries with various centerbody flame holders which included a disk, a cone, and a sphere. The extent of mixing due to these bodies was evaluated. A comparison with previously obtained experimental data yielded moderately good agreement.

Yung, C. N.↗
NASA NTRS2006

User-Defined Material Model for Progressive Failure Analysis

An overview of different types of composite material system architectures and a brief review of progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model (or UMAT) for use with the ABAQUS/Standard1 nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details and use of the UMAT subroutine are described in the present paper. Parametric studies for composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented.

Knight, Norman F. Jr.↗
NASA NTRS1979

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗
NASA NTRS1993

Recent advances and progress towards an integrated interdisciplinary thermal-structural finite element technology

An integrated finite element approach is presented for interdisciplinary thermal-structural problems. Of the various numerical approaches, finite element methods with direct time integration procedures are most widely used for these nonlinear problems. Traditionally, combined thermal-structural analysis is performed sequentially by transferring data between thermal and structural analysis. This approach is generally effective and routinely used. However, to solve the combined thermal-structural problems, this approach results in cumbersome data transfer, incompatible algorithmic representations, and different discretized element formulations. The integrated approach discussed in this paper effectively combines thermal and structural fields, thus overcoming the above major shortcomings. The approach follows Lax-Wendroff type finite element formulations with flux and stress based representations. As a consequence, this integrated approach uses common algorithmic representations and element formulations. Illustrative test examples show that the approach is effective for integrated thermal-structural problems.

Namburu, Raju R.↗
NASA NTRS2002

Computational AeroAcoustics for Fan Noise Prediction

An overview of the current state-of-the-art in computational aeroacoustics as applied to fan noise prediction at NASA Glenn is presented. Results from recent modeling efforts using three dimensional inviscid formulations in both frequency and time domains are summarized. In particular, the application of a frequency domain method, called LINFLUX, to the computation of rotor-stator interaction tone noise is reviewed and the influence of the background inviscid flow on the acoustic results is analyzed. It has been shown that the noise levels are very sensitive to the gradients of the mean flow near the surface and that the correct computation of these gradients for highly loaded airfoils is especially problematic using an inviscid formulation. The ongoing development of a finite difference time marching code that is based on a sixth order compact scheme is also reviewed. Preliminary results from the nonlinear computation of a gust-airfoil interaction model problem demonstrate the fidelity and accuracy of this approach. Spatial and temporal features of the code as well as its multi-block nature are discussed. Finally, latest results from an ongoing effort in the area of arbitrarily high order methods are reviewed and technical challenges associated with implementing correct high order boundary conditions are discussed and possible strategies for addressing these challenges ore outlined.

Envia, Ed↗
NASA NTRS1974

Lunar paleointensity from three Apollo 15 crystalline rocks using an A.R.M. method

An anhysteretic remanent magnetization method described by Banerjee and Mellema (1974) is used in the lunar paleointensity studies reported. The range of the difference in the paleointensity values obtained by the new method is not quite as great as the range of differences ordinarily found in values determined with the conventional method reported by Thellier and Thellier (1959). The results of the investigation show that the three Apollo 15 rocks studied acquired their natural remanent magnetization in a finite-sized magnetic field which was two or three orders of magnitude greater than the present ambient interplanetary field at the moon.

Banerjee, S. K.↗
NASA NTRS1983

A chimera grid scheme

A mesh system composed of multiple overset body-conforming grids is described for adapting finite-difference procedures to complex aircraft configurations. In this so-called 'chimera mesh,' a major grid is generated about a main component of the configuration and overset minor grids are used to resolve all other features. Methods for connecting overset multiple grids and modifications of flow-simulation algorithms are discussed. Computational tests in two dimensions indicate that the use of multiple overset grids can simplify the task of grid generation without an adverse effect on flow-field algorithms and computer code complexity.

Steger, J. L.↗
NASA NTRS1985

Application of a viscous-inviscid interaction method to predict transonic separated flows

A viscous-inviscid interaction numerical procedure is used to predict turbulent transonic separated flows over axisymmetric bodies. The viscous flow solutions are obtained by solving the boundary-layer equations inversely in a coupled manner by a finite-difference scheme; the inviscid flow is determined by a finite-difference solution to the full-potential equation in conservative form. Parameters affecting the convergence rate of the interaction procedure are discussed. The predictions of a new turbulence model proposed by Johnson and King are compared with experimental data and numerical solutions to the Navier-Stokes equations for flows with large regions of separation.

Lee, D. S.↗
NASA NTRS1989

Navier-Stokes Simulation Of Turbine Rotor/Stator Interaction

Patched-grid method used to effect three-dimensional Navier-Stokes numerical simulation of flow interacting with both rotor and stator in turbine. Involves use of two or more grids or sets of grids to solve finite-difference versions of Navier-Stokes or other equations of fluid flow. Ability to perform such calculations contributes to understanding of complicated turbomachinery flows and to consequent improvements in designs of turbomachinery.

Rai, Man Mohan↗
NASA NTRS1989

A method for solution of the Euler-Bernoulli beam equation in flexible-link robotic systems

An efficient numerical method for solving the partial differential equation (PDE) governing the flexible manipulator control dynamics is presented. A finite-dimensional model of the equation is obtained through discretization in both time and space coordinates by using finite-difference approximations to the PDE. An expert program written in the Macsyma symbolic language is utilized in order to embed the boundary conditions into the program, accounting for a mass carried at the tip of the manipulator. The advantages of the proposed algorithm are many, including the ability to (1) include any distributed actuation term in the partial differential equation, (2) provide distributed sensing of the beam displacement, (3) easily modify the boundary conditions through an expert program, and (4) modify the structure for running under a multiprocessor environment.

Tzes, Anthony P.↗
NASA NTRS1993

Estimation of Kalman filter gain from output residuals

This paper presents a procedure to estimate the Kalman filter gain from input-output measurement data with a given system model. The system model can be a finite element model or an experimental model from any identification method. The procedure consists of three basic steps. First, the stochastic portion related to the residuals of the response is computed. Second, the coefficients of a linear difference model for the stochastic portion are estimated by a least-squares solution that minimizes the filter residual. Third, the Kalman filter gain is computed from these model coefficients. Experimental results are presented to illustrate the usefulness of the developed procedure.

Juang, Jer-Nan↗
NASA NTRS2001

Design, Generation and Tooth Contact Analysis (TCA) of Asymmetric Face Gear Drive With Modified Geometry

A new type of face gear drive for application in transmissions, particularly in helicopters, has been developed. The new geometry differs from the existing geometry by application of asymmetric profiles and double-crowned pinion of the face gear mesh. The paper describes the computerized design, simulation of meshing and contact, and stress analysis by finite element method. Special purpose computer codes have been developed to conduct the analysis. The analysis of this new type of face gear is illustrated with a numerical example.

Litvin, Faydor L.↗
NASA NTRS1995

Supersonic Transport Analysis on the IBM Parallel System SP2

Several studies of supersonic transport (SST) configurations have been undertaken by members of the High Speed Aerodynamics branch at NASA Ames. These computational investigations involved the analysis of shapes to study the sonic boom signatures, aerodynamic performance characteristics, as well as studies of nacelle/airframe integration. A variety of different computer codes were employed including both structured and unstructured codes. The AIRPLANE code has been used extensively in these investigations. This computer code solves the Euler equations for inviscid flow by exploiting an explicit finite volume method on a mesh of tetrahedral cells. AIRPLANE is capable of handling complete aircraft configurations including nacelles and diverters. An example of a generic SST configuration is shown and a comparison of computed and experimental force coefficients is presented. Most of the computations in support of the SST investigations have been run on the YMP and C-90 computers currently installed at NASA Ames. Additional information is contained in the original extended abstract.

Jameson, Antony↗
NASA NTRS2010

Single-Receiver GPS Phase Bias Resolution

Existing software has been modified to yield the benefits of integer fixed double-differenced GPS-phased ambiguities when processing data from a single GPS receiver with no access to any other GPS receiver data. When the double-differenced combination of phase biases can be fixed reliably, a significant improvement in solution accuracy is obtained. This innovation uses a large global set of GPS receivers (40 to 80 receivers) to solve for the GPS satellite orbits and clocks (along with any other parameters). In this process, integer ambiguities are fixed and information on the ambiguity constraints is saved. For each GPS transmitter/receiver pair, the process saves the arc start and stop times, the wide-lane average value for the arc, the standard deviation of the wide lane, and the dual-frequency phase bias after bias fixing for the arc. The second step of the process uses the orbit and clock information, the bias information from the global solution, and only data from the single receiver to resolve double-differenced phase combinations. It is called "resolved" instead of "fixed" because constraints are introduced into the problem with a finite data weight to better account for possible errors. A receiver in orbit has much shorter continuous passes of data than a receiver fixed to the Earth. The method has parameters to account for this. In particular, differences in drifting wide-lane values must be handled differently. The first step of the process is automated, using two JPL software sets, Longarc and Gipsy-Oasis. The resulting orbit/clock and bias information files are posted on anonymous ftp for use by any licensed Gipsy-Oasis user. The second step is implemented in the Gipsy-Oasis executable, gd2p.pl, which automates the entire process, including fetching the information from anonymous ftp

Bertiger, William I.↗
NASA NTRS2019

Comparison of Design Tools for Stress Analysis of Adhesively Bonded Joints

Analytical or semi-analytical models for stress analysis have long been a part of initial adhesive bond design and sizing. Even with the rise of general finite element software and methods, these design models have still remained a preferred method for fast and simple joint analysis. While these methods can yield fairly representative results, the models are constructed on the foundation of geometrical and material simplifications or assumptions that allow closed form or semi-closed form solutions. This study outlines the major differences in basic assumptions for three common design software packages under use at NASA, and shows the ramifications these assumptions in a few exemplar bonded joints.

Composites↗