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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Separation of crack extension modes in orthotropic delamination models

In the analysis of an interface crack between dissimilar elastic materials, the mode of crack extension is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. This behavior currently limits the applicability of interfacial fracture mechanics as a means to predict composite delamination. The Virtual Crack Closure Technique (VCCT) is a method used to extract mode 1 and mode 2 energy release rates from numerical fracture solutions. The mode of crack extension extracted from an oscillatory solution using the VCCT is not unique due to the dependence of mode on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent crack extension modes for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior. Instead, it is argued that physically reasonable, Delta-independent modes of crack extension can be extracted from oscillatory solutions. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent modes of crack extension result. The modified technique has potential as a consistent method for extracting crack extension modes from numerical solutions. The Delta-independent modes extracted using this technique can also serve as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.

Beuth, Jack L.↗
NASA NTRS1984

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗
NASA NTRS1985

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗
NASA NTRS1986

Optimal application of Morrison's iterative noise removal for deconvolution

Morrison's iterative method of noise removal can be applied for both noise removal alone and noise removal prior to deconvolution. This method is applied to noise of various noise levels added to determine the optimum use of the method. The phase shift method of migration and modeling is evaluated and the results are compared to Stolt's approach. A method is introduced by which the optimum iterative number for deconvolution can be found. Statistical computer simulation is used to describe the optimum use of two convergent iterative techniques for seismic data. The Always-Convergent deconvolution technique was applied to data recorded during the quantitative analysis of materials through NonDestructive Evaluation (NDE) in which ultrasonic signals were used to detect flaws in substances such as composites.

Ioup, George E.↗
NASA NTRS1988

The N2-P3 complex of the evoked potential and human performance

The N2-P3 complex and other endogenous components of human evoked potential provide a set of tools for the investigation of human perceptual and cognitive processes. These multidimensional measures of central nervous system bioelectrical activity respond to a variety of environmental and internal factors which have been experimentally characterized. Their application to the analysis of human performance in naturalistic task environments is just beginning. Converging evidence suggests that the N2-P3 complex reflects processes of stimulus evaluation, perceptual resource allocation, and decision making that proceed in parallel, rather than in series, with response generation. Utilization of these EP components may provide insights into the central nervous system mechanisms modulating task performance unavailable from behavioral measures alone. The sensitivity of the N2-P3 complex to neuropathology, psychopathology, and pharmacological manipulation suggests that these components might provide sensitive markers for the effects of environmental stressors on the human central nervous system.

Odonnell, Brian F.↗
NASA NTRS1991

A geometric nonlinear degenerated shell element using a mixed formulation with independently assumed strain fields

A mixed formulation is chosen to overcome deficiencies of the standard displacement-based shell model. Element development is traced from the incremental variational principle on through to the final set of equilibrium equations. Particular attention is paid to developing specific guidelines for selecting the optimal set of strain parameters. A discussion of constraint index concepts and their predictive capability related to locking is included. Performance characteristics of the elements are assessed in a wide variety of linear and nonlinear plate/shell problems. Despite limiting the study to geometric nonlinear analysis, a substantial amount of additional insight concerning the finite element modeling of thin plate/shell structures is provided. For example, in nonlinear analysis, given the same mesh and load step size, mixed elements converge in fewer iterations than equivalent displacement-based models. It is also demonstrated that, in mixed formulations, lower order elements are preferred. Additionally, meshes used to obtain accurate linear solutions do not necessarily converge to the correct nonlinear solution. Finally, a new form of locking was identified associated with employing elements designed for biaxial bending in uniaxial bending applications.

Graf, Wiley E.↗
NASA NTRS1989

A least-squares finite element method for incompressible Navier-Stokes problem

Most finite element schemes for solving the Navier-Stokes equations can be categorized into the Galerkin mixed method and the penalty method. The mixed method leads to a saddle-point problem. In order to guarantee the existence of a solution, the combination of velocity and pressure interpolations requires satisfaction of the Ladyzhenskaya Babuska Brezzi (LBB) consistency condition which precludes the use of equal order interpolations and many seemingly natural pairs of velocity and pressure elements. In a previous paper a least-squares finite element method based on the first order velocity-pressure-vorticity formulation for the Stokes problem was proposed. This method leads to a minimization problem. The choice of combination of elements is thus not subject to the LBB condition. The numerical experiments exhibit the optimal rate of convergence for all variables with equal order interpolations. A theoretical error analysis supports the numerical results. In this paper the least-squares finite element method is extended to solving the incompressible Navier-Stokes problem.

Jiang, Bonan↗
NASA NTRS2004

Jupiter's Atmospheric Temperatures: From Voyager IRIS to Cassini CIRS

Retrievals run on Cassini Composite Infrared Spectrometer data obtained during the distant Jupiter flyby have been used to generate global temperature maps of the planet in the troposphere and stratosphere. Similar retrievals were performed on Voyager 1 IRIS data and have provided the first detailed IRIS map of the stratosphere. In both data sets, high latitude troposphere temperatures are presented for the first time, and the meridional gradients indicate the presence of circumpolar jets. Thermal winds were calculated for each data set and show strong vertical shears in the zonal winds at low latitudes. The temperatures retrieved from the two spacecraft were also compared with yearly ground-based data obtained over the intervening two decades. Tropospheric temperatures reveal gradual changes at low latitudes, with little obvious seasonal or short-term variation (Orton et al. 1994). Stratospheric temperatures show much more complicated behavior over short timescales, consistent with quasi-quadrennial oscillations at low latitudes, as suggested in prior analyses of shorter intervals of ground- based data (Orton et al. 1991, Friedson 1999). A scaling analysis indicates that meridional motions, mechanically forced by wave or eddy convergence, play an important role in modulating the temperatures and winds in the upper troposphere and stratosphere on seasonal and shorter time scales. At latitudes away from the equator, the mechanical forcing can be derived simply from a temporal record of temperature and its vertical derivative. Ground-based observations with improved vertical resolution and/or long-term monitoring from spacecraft are required for this purpose.

Simon-Miller, Amy A.↗
NASA NTRS1998

Computational Analysis of Shock Layer Emission Measurements in an Arc-Jet Facility

This paper reports computational analysis of radiation emission experiments in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. Recently, as part of ongoing arc-jet characterization work, spectroscopic radiation emission experiments have been conducted at the 20 MW NASA Ames arc-jet facility. The emission measurements were obtained from the arc-jet freestream and from a shock layer formed in front of flatfaced models. Analysis of these data is expected to provide valuable information about the thermodynamic state of the gas in the arc-jet freestream and in the shock layer as well as thermochemical equilibration processes behind the shock in arc-jet flows. Knowledge of the thermodynamic state of the gas in arc-jet test flows and especially within the shock layer is essential to interpret the heat transfer measurements such as in surface catalysis experiments. The present work is a continuation of previous work and focuses on analysis of the emission data obtained at relatively low-pressure conditions for which the arc-jet shock layer is expected to be in thermal and chemical nonequilibrium. Building blocks of the present computational analysis are: (1) simulation of nonequilibrium expanding flow in the converging-diverging conical nozzle and supersonic jet; (2) simulation of nonequilibrium shock layer formed in front of the flat-faced cylinder model; and (3) prediction of line-of-sight radiation from the computed flowfield. For computations of the nonequilibrium flow in the conical nozzle and shock layer, multi-temperature nonequilibrium codes with the axisymmetric formulation are used. For computations of line-of-sight radiation. a nonequilibrium radiation code (NEQAIR) is used to predict emission spectra from the computed flowfield. Computed line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer will be compared with those deduced from the experimental spectra. Detailed comparisons of computational and experimental spectra will also be presented.

Gokcen, Tahir↗
NASA NTRS1992

Verifying the interactive convergence clock synchronization algorithm using the Boyer-Moore theorem prover

The application of formal methods to the analysis of computing systems promises to provide higher and higher levels of assurance as the sophistication of our tools and techniques increases. Improvements in tools and techniques come about as we pit the current state of the art against new and challenging problems. A promising area for the application of formal methods is in real-time and distributed computing. Some of the algorithms in this area are both subtle and important. In response to this challenge and as part of an ongoing attempt to verify an implementation of the Interactive Convergence Clock Synchronization Algorithm (ICCSA), we decided to undertake a proof of the correctness of the algorithm using the Boyer-Moore theorem prover. This paper describes our approach to proving the ICCSA using the Boyer-Moore prover.

Young, William D.↗
NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. III - Analysis of multiheats

An outline is presented of approaches for treating multiheats on the basis of the focal point convergence method. The method has been employed in the case of two highly characterized multiheats, including a 304 stainless steel studied in Japan and a low alloy carbon steel studied in England. The method makes use of the same functional form for all materials. Only the constants are varied for each multiheat. Completely computerized procedures are employed for the determination of the constants. Once the basic analysis has been performed, the representation of various members in the same system is achieved by adding linear expressions of log stress, changing only two constants in the equations to represent a selected heat.

Manson, S. S.↗
NASA NTRS1986

Exact and approximate gas dynamics using the tangent gas

For the determination of aerodynamic characteristics such as lift, drag, and moment coefficients, it is crucial to compute the properties of steady flow past an airfoil. This investigation provides a set of flow dependent grid systems and initial flowfield guesses which substantially improve convergence rates when applied to the Euler equations for flows past an airfoil. The basic equations are examined, taking into account nonlinear equations which are difficult to solve. A good approximation to these equations under certain conditions can be obtained by introducing the so-called 'tangent gas approximation' considered by Woods (1961), in which the isentropic relation between rho and p is replaced by a tangent to a curve. Attention is given to the solution procedure, the analysis (direct) problem, and a comparison of the tangent gas solution with the converged Euler solution.

Daripa, P. K.↗
NASA NTRS1987

Extrapolation methods for vector sequences

This paper derives, describes, and compares five extrapolation methods for accelerating convergence of vector sequences or transforming divergent vector sequences to convergent ones. These methods are the scalar epsilon algorithm (SEA), vector epsilon algorithm (VEA), topological epsilon algorithm (TEA), minimal polynomial extrapolation (MPE), and reduced rank extrapolation (RRE). MPE and RRE are first derived and proven to give the exact solution for the right 'essential degree' k. Then, Brezinski's (1975) generalization of the Shanks-Schmidt transform is presented; the generalized form leads from systems of equations to TEA. The necessary connections are then made with SEA and VEA. The algorithms are extended to the nonlinear case by cycling, the error analysis for MPE and VEA is sketched, and the theoretical support for quadratic convergence is discussed. Strategies for practical implementation of the methods are considered.

Smith, David A.↗
NASA NTRS1993

Uniform convergence of multigrid V-cycle iterations for indefinite and nonsymmetric problems

In this paper, we present an analysis of a multigrid method for nonsymmetric and/or indefinite elliptic problems. In this multigrid method various types of smoothers may be used. One type of smoother which we consider is defined in terms of an associated symmetric problem and includes point and line, Jacobi, and Gauss-Seidel iterations. We also study smoothers based entirely on the original operator. One is based on the normal form, that is, the product of the operator and its transpose. Other smoothers studied include point and line, Jacobi, and Gauss-Seidel. We show that the uniform estimates for symmetric positive definite problems carry over to these algorithms. More precisely, the multigrid iteration for the nonsymmetric and/or indefinite problem is shown to converge at a uniform rate provided that the coarsest grid in the multilevel iteration is sufficiently fine (but not depending on the number of multigrid levels).

Bramble, James H.↗
NASA NTRS2012

Influence of Primary Gage Sensitivities on the Convergence of Balance Load Iterations

The connection between the convergence of wind tunnel balance load iterations and the existence of the primary gage sensitivities of a balance is discussed. First, basic elements of two load iteration equations that the iterative method uses in combination with results of a calibration data analysis for the prediction of balance loads are reviewed. Then, the connection between the primary gage sensitivities, the load format, the gage output format, and the convergence characteristics of the load iteration equation choices is investigated. A new criterion is also introduced that may be used to objectively determine if the primary gage sensitivity of a balance gage exists. Then, it is shown that both load iteration equations will converge as long as a suitable regression model is used for the analysis of the balance calibration data, the combined influence of non linear terms of the regression model is very small, and the primary gage sensitivities of all balance gages exist. The last requirement is fulfilled, e.g., if force balance calibration data is analyzed in force balance format. Finally, it is demonstrated that only one of the two load iteration equation choices, i.e., the iteration equation used by the primary load iteration method, converges if one or more primary gage sensitivities are missing. This situation may occur, e.g., if force balance calibration data is analyzed in direct read format using the original gage outputs. Data from the calibration of a six component force balance is used to illustrate the connection between the convergence of the load iteration equation choices and the existence of the primary gage sensitivities.

Ulbrich, Norbert Manfred↗
NASA NTRS1999

Chemometric analysis for extraction of individual fluorescence spectrum and lifetimes from a target mixture

The present invention is a system for chemometric analysis for the extraction of the individual component fluorescence spectra and fluorescence lifetimes from a target mixture. The present invention combines a processor with an apparatus for generating an excitation signal to transmit at a target mixture and an apparatus for detecting the emitted signal from the target mixture. The present invention extracts the individual fluorescence spectrum and fluorescence lifetime measurements from the frequency and wavelength data acquired from the emitted signal. The present invention uses an iterative solution that first requires the initialization of several decision variables and the initial approximation determinations of intermediate matrices. The iterative solution compares the decision variables for convergence to see if further approximation determinations are necessary. If the solution converges, the present invention then determines the reduced best fit error for the analysis of the individual fluorescence lifetime and the fluorescence spectrum before extracting the individual fluorescence lifetime and fluorescence spectrum from the emitted signal of the target mixture.

Hallidy, William H.↗
NASA NTRS1999

Urban Heat Islands and Summertime Convective Thunderstorms in Atlanta: Three Case Studies

Data from both 27 sites in the Atlanta mesonet surface meteorological network and eight National Weather Service sites were analyzed for the period from 26 July to 3 August 1996. Analysis of the six precipitation events over the city during the period (each on a different day) showed that its urban heat island (UHI) induced a convergence zone that initiated three of the storms at different times of the day, i.e., 0630,0845, and 1445 EDT. Previous analysis has shown that New York City (NYC) effects summer daytime thunderstorm formation and/or movement. That study found that during nearly calm regional flow conditions the NYC UHI initiates convective activity. Moving thunderstorms, however, tended to bifurcate and to move around the city, due to its building barrier effect. The current Atlanta results thus agree with the NYC results with respect to thunderstorm initiation.

Bornstein, Robert↗
NASA NTRS1993

A coarse-grid correction/nonlinear relaxation algorithm for the three-dimensional, compressible Navier-Stokes equations

A combined coarse-grid correction/upwind relaxation strategy to provide rapid convergence for 3D high-speed viscous flowfields is discussed an evaluated. The construction and analysis of a simple two-grid acceleration procedure based on 'hyperbolic' multigrid concepts is presented. Numerical simulations of a 2D compression-corner flowfield, a 3D crossing shock/turbulent boundary layer interaction, and a 3D scramjet inlet flowfield are presented to illustrate the benefits of the approach. Results indicate that the procedure generally converges two or more times faster than the baseline algorithm.

Edwards, Jack R.↗