Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigenvalues”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Separation analysis, a tool for analyzing multigrid algorithms

The separation of vectors by multigrid (MG) algorithms is applied to the study of convergence and to the prediction of the performance of MG algorithms. The separation operator for a two level cycle algorithm is derived. It is used to analyze the efficiency of the cycle when mixing of eigenvectors occurs. In particular cases the separation analysis reduces to Fourier type analysis. The separation operator of a two level cycle for a Schridubger eigenvalue problem, is derived and analyzed in a Fourier basis. Separation analysis gives information on how to choose performance relaxations and inter-level transfers. Separation analysis is a tool for analyzing and designing algorithms, and for optimizing their performance.

Costiner, Sorin↗
NASA NTRS1994

Control theory and splines, applied to signature storage

In this report the problem we are going to study is the interpolation of a set of points in the plane with the use of control theory. We will discover how different systems generate different kinds of splines, cubic and exponential, and investigate the effect that the different systems have on the tracking problems. Actually we will see that the important parameters will be the two eigenvalues of the control matrix.

Enqvist, Per↗
NASA NTRS1995

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗
NASA NTRS1995

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗
NASA NTRS1995

Aeroelastic stability of wind turbine blade/aileron systems

Aeroelastic stability analyses have been performed for the MOD-5A blade/aileron system. Various configurations having different aileron torsional stiffness, mass unbalance, and control system damping have been investigated. The analysis was conducted using a code recently developed by the General Electric Company - AILSTAB. The code extracts eigenvalues for a three degree of freedom system, consisting of: (1) a blade flapwise mode; (2) a blade torsional mode; and (3) an aileron torsional mode. Mode shapes are supplied as input and the aileron can be specified over an arbitrary length of the blade span. Quasi-steady aerodynamic strip theory is used to compute aerodynamic derivatives of the wing-aileron combination as a function of spanwise position. Equations of motion are summarized herein. The program provides rotating blade stability boundaries for torsional divergence, classical flutter (bending/torsion) and wing/aileron flutter. It has been checked out against fixed-wing results published by Theodorsen and Garrick. The MOD-5A system is stable with respect to divergence and classical flutter for all practical rotor speeds. Aileron torsional stiffness must exceed a minimum critical value to prevent aileron flutter. The nominal control system stiffness greatly exceeds this minimum during normal operation. The basic system, however, is unstable for the case of a free (or floating) aileron. The instability can be removed either by the addition of torsional damping or mass-balancing the ailerons. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Strain, J. C.↗
NASA NTRS1995

Combustion-acoustic stability analysis for premixed gas turbine combustors

Lean, prevaporized, premixed combustors are susceptible to combustion-acoustic instabilities. A model was developed to predict eigenvalues of axial modes for combustion-acoustic interactions in a premixed combustor. This work extends previous work by including variable area and detailed chemical kinetics mechanisms, using the code LSENS. Thus the acoustic equations could be integrated through the flame zone. Linear perturbations were made of the continuity, momentum, energy, chemical species, and state equations. The qualitative accuracy of our approach was checked by examining its predictions for various unsteady heat release rate models. Perturbations in fuel flow rate are currently being added to the model.

Darling, Douglas↗
NASA NTRS1990

Partitioning sparse matrices with eigenvectors of graphs

The problem of computing a small vertex separator in a graph arises in the context of computing a good ordering for the parallel factorization of sparse, symmetric matrices. An algebraic approach for computing vertex separators is considered in this paper. It is shown that lower bounds on separator sizes can be obtained in terms of the eigenvalues of the Laplacian matrix associated with a graph. The Laplacian eigenvectors of grid graphs can be computed from Kronecker products involving the eigenvectors of path graphs, and these eigenvectors can be used to compute good separators in grid graphs. A heuristic algorithm is designed to compute a vertex separator in a general graph by first computing an edge separator in the graph from an eigenvector of the Laplacian matrix, and then using a maximum matching in a subgraph to compute the vertex separator. Results on the quality of the separators computed by the spectral algorithm are presented, and these are compared with separators obtained from other algorithms for computing separators. Finally, the time required to compute the Laplacian eigenvector is reported, and the accuracy with which the eigenvector must be computed to obtain good separators is considered. The spectral algorithm has the advantage that it can be implemented on a medium-size multiprocessor in a straightforward manner.

Pothen, Alex↗
NASA NTRS1994

Nonlinear equilibration of localized instabilities on a baroclinic jet

Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Delsole, Timothy M.↗
NASA NTRS1994

Nonlinear stability of Halley comethosheath with transverse plasma motion

Weakly nonlinear Magneto Hydrodynamic (MHD) stability of the Halley cometosheath determined by the balance between the outward ion-neutral drag force and the inward Lorentz force is investigated including the transverse plasma motion as observed in the flanks with the help of the method of multiple scales. The eigenvalues and the eigenfunctions are obtained for the linear problem and the time evolution of the amplitude is obtained using the solvability condition for the solution of the second order problem. The diamagnetic cavity boundary and the adjacent layer of about 100 km thickness is found unstable for the travelling waves of certain wave numbers. Halley ionopause has been observed to have strong ripples with a wavelength of several hundred kilometers. It is found that nonlinear effects have stabilizing effect.

Srivastava, Krishna M.↗
NASA NTRS1994

Large-angle cosmic microwave background anisotropies in an open universe

If the universe is open, scales larger than the curvature scale may be probed by observation of large-angle fluctuations in the cosmic microwave background (CMB). We consider primordial adiabatic perturbations and discuss power spectra that are power laws in volume, wavelength, and eigenvalue of the Laplace operator. Such spectra may have arisen if, for example, the universe underwent a period of `frustated' inflation. The resulting large-angle anisotropies of the CMB are computed. The amplitude generally increases as Omega is decreased but decreases as h is increased. Interestingly enough, for all three Ansaetze, anisotropies on angular scales larger than the curvature scale are suppressed relative to the anisotropies on scales smaller than the curvature scale, but cosmic variance makes discrimination between various models difficult. Models with 0.2 approximately less than Omega h approximately less than 0.3 appear compatible with CMB fluctuations detected by Cosmic Background Explorer Satellite (COBE) and the Tenerife experiment and with the amplitude and spectrum of fluctuations of galaxy counts in the APM, CfA, and 1.2 Jy IRAS surveys. COBE normalization for these models yields sigma(sub 8) approximately = 0.5 - 0.7. Models with smaller values of Omega h when normalized to COBE require bias factors in excess of 2 to be compatible with the observed galaxy counts on the 8/h Mpc scale. Requiring that the age of the universe exceed 10 Gyr implies that Omega approximately greater than 0.25, while requiring that from the last-scattering term in the Sachs-Wolfe formula, large-angle anisotropies come primarily from the decay of potential fluctuations at z approximately less than 1/Omega. Thus, if the universe is open, COBE has been detecting temperature fluctuations produced at moderate redshift rather than at z approximately 1300.

Kamionkowski, Marc↗
NASA NTRS1995

The use of optimal polarizations for studying the microphysics of precipitation: Nonattenuating wavelengths

The objective of this work is to explore relationships between the microphysical properties of precipitation and optimal polarizations. The dependence of three optimal polarization parameters (asymmetry ratio A, optimal tilt tau(sub op), and optimal ellipticity epsilon(sub op) on the reflectivity-weighted mean drop shape, mean canting angle, and standard deviation of a Gaussian canting angle distribution is studied. This is accomplished by using computer simulations that provide the rms scattering matrix for an ensemble of canted drops with a prescribed two-parameter canting angle distribution. Also examined are the effects of propagation on the polarization parameters for nonattenuating wavelengths. The asymmetry ratio A is simply the ratio of the maximal to minimal total backscattered energy (ratio of the largest and smallest eigenvalue of the Graves power matrix G is identically equal to S dagger S). Similar to Z(sub DR), this ratio decreases with increasing mean axial ratio, but unlike Z(sub DR), it is not affected by canting (for a single drop). The dependence of A on the reflectivity-weighted mean drop shape is examined, and a power-law relationship similar to that which exists for Z(sub DR) is established. The asymmetry ratio A can be regarded as a generalization of Z(sub DR) because it requires only a measurement of linear depolarization ratio (in addition to Z(sub DR)), is independent of the propagation phase, and is less sensitive to canting. In a similar manner, the dependence of optimal ellipticity and tilt on the microphysical parameters is studied. In particular, it appears that the rms tilt of the optimal polarization ellipse is proportional to the variance of the canting angle distribution. Several other promising relationships between optimal polarizations and the microphysical variables of an ensemble of hydrometeors are also discussed.

Kwiatkowski, John M.↗
NASA NTRS1993

Linear static structural and vibration analysis on high-performance computers

Parallel computers offer the oppurtunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on massively-parallel computers hereafter referred to as Scalable High-Performance Computers (SHPC), for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. models for High-Speed Civil Transport). The goal of this research is to develop a new, efficient technique which extends structural analysis to SHPC and makes large-scale structural analyses tractable.

Baddourah, M. A.↗
NASA NTRS1994

Signatures of bulge triaxiality from kinematics in Baade's window

We study a sample of 62 Baade's Window, (l,b) = (1, -4)deg, K giants that have published proper motions, radial velocity, and metallicity. Using R(sub 0) = 8 kpc, we construct the velocity ellipsoids, namely the 3x3 velocity dispersion tensors, for the metal rich stars ((Fe/H) greater than or equal to 0) and metal poor stars ((Fe/H) less than or equal to -0.2). After diagonalizing the tensor, we find a vertex deviation characteristic of a nonaxisymmetric system. Eigenvalues for the two velocity ellipsoids (sigma(sub 1), sigma(sub 2), sigma(sub 3)) are (126, 89, 65) +/- 13 km/s for the metal rich sample and (154, 77, 83) +/- 25 km/s for the metal poor sample with their long axes pointing to two nearly perpendicular directions (l(sub v), b(sub v)) = (-65 +/- 9 deg, +14 +/- 9 deg) and (l(sub v), b(sub v)) = (25 +/- 14 deg, -11 +/- 14 deg), respectively. The vertex deviations of the velocity ellipsoids cannot be consistently explained by any oblate model. We are able to reject the hypothesis that the metal poor and metal rich populations are drawn from the same distribution at better than the 97% confidence level. We populate orbits in a realistic bar potential with a Gaussian velocity distribution, allowing us to simulate and interpret observations. We conclude that the data are consistent with a triaxial bulge pointing towards (l,b) with l less than 0 deg and b = 0 deg as suggested by earlier work on gas dynamics and the observed light distribution. We also predict that low latitude (absolute value of b less than or equal to 4 deg) bulge fields should show the vertex deviation more strongly and would therefore be the best locations for future proper motion studies. In the classification scheme of Athanassoula et al. (1983) the metal rich stars appear to occupy the B-family orbits which rotate in the prograde sense in the rest frame and have boxy shapes that are aligned with and supporting the bar. The metal poor stars in the sample lag behind the metal rich bulge and appear to occupy R-family orbits which rotate in the retrograde sense in the rest frame. They have nearly round loop shapes and are aligned perpendicularly to the bar, hence limit the triaxiality of the bar potential. The correlations between the metallicity and the orbit families can develop if the bulge forms dissipatively on a sufficiently long time scale. However, it is difficult to explain such correlations if most stars in the inner Galaxy form during the violent relaxation phase.

Zhao, Hongsheng↗
NASA NTRS1994

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗
NASA NTRS1995

Temporal variation of meandering intensity and domain-wide lateral oscillations of the Gulf Stream

The path of the Gulf Stream exhibits two modes of variability: wavelike spatial meanders associated with instability processes and large-sale lateral shifts of the path presumably due to atmospheric forcing. The objectives of this study are to examine the temporal variation of the intensity of spatial meandering in the stream, to characterize large-scale lateral oscillations in the stream's path, and to study the correlation betwen these two dynamically distinct modes of variability. The data used for this analysis are path displacemets ofthe Gulf Stream between 75 deg and 60 deg W obtained from AVHRR-derived (Advanced Very High Resolution Radiometer) infrared images for the period April 1982 through December 1989. Meandering intensity, measured by the spatial root-mean-sqaure displacement of the stream path, displays a 9-month dominant periodicity which is persistent through the study period. The 9-month fluctuation in meandering intensity may be related to the interaction of Rosseby waves with the stream. Interannual variation of meandering intensity is also found to be significant, with meandering being mich more intense during 1985 than it was in 1987. Annual variation, however,is weak and not well-defined.The spatially averaged position of the stream, which reflects nonmeandering large-scale lateral oscillations of the stream path, is dominated by an annual cycle. On average, the mean position is farthest north in November and farthest south in April. The first empirical orthogonal function mode of the space-time path displacements represents lateral oscillatins that are in-phase over the space-time domain. Interannual oscillations are also observed and are found to be weaker than the annual oscillation. The eigenvalue of the first mode indicates that about 21.5% of the total space-time variability of the stream path can be attibuted to domain-wide lateral oscillation. The correlation between meandering intensity and domain-wide lateral oscillations is very weak.

Lee, Tong↗
NASA NTRS1995

Asymptotic Spectra Of Banded Quasi-Toeplitz Matrices

Paper presents theoretical and numerical study of asymptotic spectra of eigenvalues of banded Toeplitz and quasi-Toeplitz matrices. Emphasis in study on non-normal banded Toeplitz and quasi-Toeplitz matrices of arbitrarily large order and relatively small bandwidth.

Beam, Richard↗
NASA NTRS1995

Optimization of composite box-beam structures including effects of subcomponent interactions

Minimum mass designs are obtained for a simple box beam structure subject to bending, torque and combined bending/torque load cases. These designs are obtained subject to point strain and linear buckling constraints. The present work differs from previous efforts in that special attention is payed to including the effects of subcomponent panel interaction in the optimal design process. Two different approaches are used to impose the buckling constraints. When the global approach is used, buckling constraints are imposed on the global structure via a linear eigenvalue analysis. This approach allows the subcomponent panels to interact in a realistic manner. The results obtained using this approach are compared to results obtained using a traditional, less expensive approach, called the local approach. When the local approach is used, in-plane loads are extracted from the global model and used to impose buckling constraints on each subcomponent panel individually. In the global cases, it is found that there can be significant interaction between skin, spar, and rib design variables. This coupling is weak or nonexistent in the local designs. It is determined that weight savings of up to 7% may be obtained by using the global approach instead of the local approach to design these structures. Several of the designs obtained using the linear buckling analysis are subjected to a geometrically nonlinear analysis. For the designs which were subjected to bending loads, the innermost rib panel begins to collapse at less than half the intended design load and in a mode different from that predicted by linear analysis. The discrepancy between the predicted linear and nonlinear responses is attributed to the effects of the nonlinear rib crushing load, and the parameter which controls this rib collapse failure mode is shown to be the rib thickness. The rib collapse failure mode may be avoided by increasing the rib thickness above the value obtained from the (linear analysis based) optimizer. It is concluded that it would be necessary to include geometric nonlinearities in the design optimization process if the true optimum in this case were to be found.

Ragon, Scott A.↗
NASA NTRS1995

Predictions of thermal buckling strengths of hypersonic aircraft sandwich panels using minimum potential energy and finite element methods

Thermal buckling characteristics of hypersonic aircraft sandwich panels of various aspect ratios were investigated. The panel is fastened at its four edges to the substructures under four different edge conditions and is subjected to uniform temperature loading. Minimum potential energy theory and finite element methods were used to calculate the panel buckling temperatures. The two methods gave fairly close buckling temperatures. However, the finite element method gave slightly lower buckling temperatures than those given by the minimum potential energy theory. The reasons for this slight discrepancy in eigensolutions are discussed in detail. In addition, the effect of eigenshifting on the eigenvalue convergence rate is discussed.

Ko, William L.↗