Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Supplement of: The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El NioSouthern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium.An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records.Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the amount effect. During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher↗
NASA NTRS1990

Free wake analysis of hover performance using a new influence coefficient method

A new approach to the prediction of helicopter rotor performance using a free wake analysis was developed. This new method uses a relaxation process that does not suffer from the convergence problems associated with previous time marching simulations. This wake relaxation procedure was coupled to a vortex-lattice, lifting surface loads analysis to produce a novel, self contained performance prediction code: EHPIC (Evaluation of Helicopter Performance using Influence Coefficients). The major technical features of the EHPIC code are described and a substantial amount of background information on the capabilities and proper operation of the code is supplied. Sample problems were undertaken to demonstrate the robustness and flexibility of the basic approach. Also, a performance correlation study was carried out to establish the breadth of applicability of the code, with very favorable results.

Quackenbush, Todd R.↗
NASA NTRS1988

Modal cost analysis for simple continua

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations, it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for simple continua such as beam-like structures. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗
NASA NTRS1974

Thermal stress analysis of reusable surface insulation for shuttle

An iterative procedure for accurately determining tile stresses associated with static mechanical and thermally induced internal loads is presented. The necessary conditions for convergence of the method are derived. An user-oriented computer program based upon the present method of analysis was developed. The program is capable of analyzing multi-tiled panels and determining the associated stresses. Typical numerical results from this computer program are presented.

Ojalvo, I. U.↗
NASA NTRS1998

Vertical and Horizontal Wind Structure Prior to and During the 1997/1998 ENSO Event

Analyses of ECMWF and satellite wind products are used to investigate low frequency fluctuations in the wind fields immediately preceding and during the 1997/1998 ENSO event. An extended EOF analysis of the horizontal and vertical wind structure that allows for the study of propagating features indicates that the dominant mode of variability established in the lower levels of the atmosphere is a convergence in the central tropical Pacific. At mid-levels, the westerlies increase with time over most of the equatorial Pacific, and there are closed cyclonic features on either side of the equator. At upper levels, the flow tends to become divergent across the equatorial Pacific, and there is a strong flow from the northwest Pacific across Indonesia that then turns to the west along the equator in the Indian Ocean. The addition of satellite data for improved surface wind estimates reveal that the Indian Ocean winds were anomalously strong during the winter preceding the initiation of the ENSO event. An analysis of the momentum convergence flux at the surface indicates a very strong anomaly 9 months prior to the ENSO event that begins in the western Indian Ocean and propagates to the east, arriving over the Pacific warm pool area coincident with the beginning of the El Nino. Changes in the near surface thermal and water vapor structure over the Indian Ocean during this period will also be discussed.

Adamec, David↗
NASA NTRS2009

New Flutter Analysis Technique for CFD-based Unsteady Aeroelasticity

This paper presents a flutter analysis technique for the transonic flight regime. The technique uses an iterative approach to determine the critical dynamic pressure for a given mach number. Unlike other CFD-based flutter analysis methods, each iteration solves for the critical dynamic pressure and uses this value in subsequent iterations until the value converges. This process reduces the iterations required to determine the critical dynamic pressure. To improve the accuracy of the analysis, the technique employs a known structural model, leaving only the aerodynamic model as the unknown. The aerodynamic model is estimated using unsteady aeroelastic CFD analysis combined with a parameter estimation routine. The technique executes as follows. The known structural model is represented as a finite element model. Modal analysis determines the frequencies and mode shapes for the structural model. At a given mach number and dynamic pressure, the unsteady CFD analysis is performed. The output time history of the surface pressure is converted to a nodal aerodynamic force vector. The forces are then normalized by the given dynamic pressure. A multi-input multi-output parameter estimation software, ERA, estimates the aerodynamic model through the use of time histories of nodal aerodynamic forces and structural deformations. The critical dynamic pressure is then calculated using the known structural model and the estimated aerodynamic model. This output is used as the dynamic pressure in subsequent iterations until the critical dynamic pressure is determined. This technique is demonstrated on the Aerostructures Test Wing-2 model at NASA's Dryden Flight Research Center.

Pak, Chan-gi↗
NASA NTRS1998

Asymptotic Spreading Rate of Initially Compressible Jets: Experiment and Analysis

Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub J). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub J), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Zaman, K. B. M. Q.↗
NASA NTRS1998

Asymptotic Spreading Rate of Initially Compressible Jets-Experiment and Analysis

Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub j). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub j), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Zaman, K. B. M. Q.↗
NASA NTRS1987

Error analysis for spectral approximation of the Korteweg-De Vries equation

The conservation and convergence properties of spectral Fourier methods for the numerical approximation of the Korteweg-de Vries equation are analyzed. It is proved that the (aliased) collocation pseudospectral method enjoys the same convergence properties as the spectral Galerkin method, which is less effective from the computational point of view. This result provides a precise mathematical answer to a question raised by several authors in recent years.

Maday, Y.↗
NASA NTRS1994

Structural/aerodynamic Blade Analyzer (SAB) User's Guide, Version 1.0

The structural/aerodynamic blade (SAB) analyzer provides an automated tool for the static-deflection analysis of turbomachinery blades with aerodynamic and rotational loads. A structural code calculates a deflected blade shape using aerodynamic loads input. An aerodynamic solver computes aerodynamic loads using deflected blade shape input. The two programs are iterated automatically until deflections converge. Currently, SAB version 1.0 is interfaced with MSC/NASTRAN to perform the structural analysis and PROP3D to perform the aerodynamic analysis. This document serves as a guide for the operation of the SAB system with specific emphasis on its use at NASA Lewis Research Center (LeRC). This guide consists of six chapters: an introduction which gives a summary of SAB; SAB's methodology, component files, links, and interfaces; input/output file structure; setup and execution of the SAB files on the Cray computers; hints and tips to advise the user; and an example problem demonstrating the SAB process. In addition, four appendices are presented to define the different computer programs used within the SAB analyzer and describe the required input decks.

Morel, M. R.↗
NASA NTRS1993

Acceleration of convergence and spectrum transformation of implicit finite difference operators associated with Navier-Stokes equations

Implicit difference operator spectra are presently computed by applying eigensystem analysis techniques to finite-difference formulations of 2D Euler and Navier-Stokes equations, and attention is given to these iterative methods' convergence and stability characteristics by taking into account the effects of grid geometry, time-step, numerical viscosity, and boundary conditions. On the basis of the eigenvalue distributions for various flow configurations, the feasibility of applying such convergence-acceleration techniques as eigenvalue annihilation and relaxation is discussed. Spectrum-shifting is applied to NASA-Ames' ARC2D flow code, achieving a 20-33 percent efficiency.

Saleem, M.↗
NASA NTRS1979

An investigation of several factors involved in a finite difference procedure for analyzing the transonic flow about harmonically oscillating airfoils and wings

Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. Since sinusoidal motion is assumed, the unsteady equation is independent of time. Three finite difference investigations are discussed including a new operator for mesh points with supersonic flow, the effects on relaxation solution convergence of adding a viscosity term to the original differential equation, and an alternate and relatively simple downstream boundary condition. A method is developed which uses a finite difference procedure over a limited inner region and an approximate analytical procedure for the remaining outer region. Two investigations concerned with three-dimensional flow are presented. The first is the development of an oblique coordinate system for swept and tapered wings. The second derives the additional terms required to make row relaxation solutions converge when mixed flow is present. A finite span flutter analysis procedure is described using the two-dimensional unsteady transonic program with a full three-dimensional steady velocity potential.

Ehlers, F. E.↗
NASA NTRS1986

Divergence study of a high-aspect ratio, forward-swept wing

A study has been conducted in the NASA-Langley Transonic Dynamics wind tunnel to determine the divergence characteristics of a high aspect ratio forward swept wing of rectangular planform, for sweep angles of zero, -15, -30, -45, and -60 deg. In addition to a rectangular wingtip, a tip whose geometry lay parallel with the freestream flow was tested at a sweep angle of -45 deg. The primary objective was to obtain data that could be used to verify the divergence prediction capabilities of an aeroelastic analysis code employing kernel function aerodynamics. The analytical predictions of convergence were found to be conservative for all forward sweep angles, and the effect of the two tip shapes on divergence dynamic pressure were accurately predicted.

Cole, S. R.↗
NASA NTRS1989

Comparison of measured induced velocities with results from a closed-form finite state wake model in forward flight

An unsteady, three-dimensional induced-flow model is applied to compute the induced-flow distribution of a rotor in forward flight, and numerical results are compared against LDV measurements for both time-averaged and unsteady induced flow at the disk. The former flow shows good agreement with measured data, except just behind the pylon at the lowest advanced ratio and near the blade tips for rectangular blades at high advance ratios. The method performs as well or better than other codes that have been applied to these data, takes less computing time, and is better suited to aeroelastic analysis. Results with only four harmonics and 15 state variables converge to all fundamental characteristics of the time-averaged flow.

Peters, David A.↗
NASA NTRS1990

Development of load-dependent Ritz vector method for structural dynamic analysis of large space structures

The development and preliminary assessment of a method for dynamic structural analysis based on load-dependent Ritz vectors are presented. The vector basis is orthogonalized with respect to the mass and structural stiffness in order that the equations of motion can be uncoupled and efficient analysis of large space structure performed. A series of computer programs was developed based on the algorithm for generating the orthogonal load-dependent Ritz vectors. Transient dynamic analysis performed on the Space Station Freedom using the software was found to provide solutions that require a smaller number of vectors than the modal analysis method. Error norm based on the participation of the mass distribution of the structure and spatial distribution of structural loading, respectively, were developed in order to provide an indication of vector truncation. These norms are computed before the transient analysis is performed. An assessment of these norms through a convergence study of the structural response was performed. The results from this assessment indicate that the error norms can provide a means of judging the quality of the vector basis and accuracy of the transient dynamic solution.

Ricles, James M.↗
NASA NTRS1993

Increased depth-diameter ratios in the Medusae Fossae Formation deposits of Mars

Depth to diameter ratios for fresh impact craters on Mars are commonly cited as approximately 0.2 for simple craters and 0.1 for complex craters. Recent computation of depth-diameter ratios in the Amazonis-Memnonia region of Mars indicates that craters within the Medusae Fossae Formation deposits found in this region display greater depth-diameter ratios than expected for both simple and complex craters. Photoclinometric and shadow length techniques have been used to obtain depths of craters within the Amazonis-Memnonia region. The 37 craters in the 2 to 29 km diameter range and displaying fresh impact morphologies were identified in the area of study. This region includes the Amazonian aged upper and middle members of the Medusae Fossae Formation and Noachian aged cratered and hilly units. The Medusae Fossae Formation is characterized by extensive, flat to gently undulating deposits of controversial origin. These deposits appear to vary from friable to indurated. Early analysis of crater degradation in the Medusae Fossae region suggested that simple craters excavated to greater depths than expected based on the general depth-diameter relationships derived for Mars. However, too few craters were available in the initial analysis to estimate the actual depth-diameter ratios within this region. Although the analysis is continuing, we are now beginning to see a convergence towards specific values for the depth-diameter ratio depending on geologic unit.

Barlow, N. G.↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS1996

Schwarz Methods: To Symmetrize or not to Symmetrize

A preconditioning theory for Schwarz methods is presented. The theory establishes sufficient conditions for multiplicative and additive Schwarz algorithms to yield self-adjoint positive definite preconditioners. It allows for the analysis and use of non-variational and non-convergent linear methods as preconditioners for conjugate gradient methods, and it is applied to domain decomposition and multigrid. This paper illustrates why symmetrizing may be a bad idea for linear methods. Numerical examples are presented for a test problem.

Holst, Michael↗