Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigenvalues and eigenfunctions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Improved Finite Element Modeling of the Turbofan Engine Inlet Radiation Problem

Improvements have been made in the finite element model of the acoustic radiated field from a turbofan engine inlet in the presence of a mean flow. The problem of acoustic radiation from a turbofan engine inlet is difficult to model numerically because of the large domain and high frequencies involved. A numerical model with conventional finite elements in the near field and wave envelope elements in the far field has been constructed. By employing an irrotational mean flow assumption, both the mean flow and the acoustic perturbation problem have been posed in an axisymmetric formulation in terms of the velocity potential; thereby minimizing computer storage and time requirements. The finite element mesh has been altered in search of an improved solution. The mean flow problem has been reformulated with new boundary conditions to make it theoretically rigorous. The sound source at the fan face has been modeled as a combination of positive and negative propagating duct eigenfunctions. Therefore, a finite element duct eigenvalue problem has been solved on the fan face and the resulting modal matrix has been used to implement a source boundary condition on the fan face in the acoustic radiation problem. In the post processing of the solution, the acoustic pressure has been evaluated at Gauss points inside the elements and the nodal pressure values have been interpolated from them. This has significantly improved the results. The effect of the geometric position of the transition circle between conventional finite elements and wave envelope elements has been studied and it has been found that the transition can be made nearer to the inlet than previously assumed.

Roy, Indranil Danda↗
NASA NTRS1988

Multigrid method for stability problems

The problem of calculating the stability of steady state solutions of differential equations is addressed. Leading eigenvalues of large matrices that arise from discretization are calculated, and an efficient multigrid method for solving these problems is presented. The resulting grid functions are used as initial approximations for appropriate eigenvalue problems. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a nonstandard way in which the right-hand side of the coarse grid equations involves unknown parameters to be solved on the coarse grid. This leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem are presented which demonstrate the effectiveness of the method.

Ta'asan, Shlomo↗
NASA NTRS1988

Multigrid method for stability problems

The problem of calculating the stability of steady state solutions of differential equations is treated. Leading eigenvalues (i.e., having maximal real part) of large matrices that arise from discretization are to be calculated. An efficient multigrid method for solving these problems is presented. The method begins by obtaining an initial approximation for the dominant subspace on a coarse level using a damped Jacobi relaxation. This proceeds until enough accuracy for the dominant subspace has been obtained. The resulting grid functions are then used as an initial approximation for appropriate eigenvalue problems. These problems are being solved first on coarse levels, followed by refinement until a desired accuracy for the eigenvalues has been achieved. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a non-standard way in which the right hand side of the coarse grid equations involves unknown parameters to be solved for on the coarse grid. This in particular leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem demonstrate the effectiveness of the method proposed. Using an FMG algorithm a solution to the level of discretization errors is obtained in just a few work units (less than 10), where a work unit is the work involved in one Jacobi relization on the finest level.

Taasan, Shlomo↗
NASA NTRS1990

On the instabilities of supersonic mixing layers - A high-Mach-number asymptotic theory

The stability of a family of tanh mixing layers is studied at large Mach numbers using perturbation methods. It is found that the eigenfunction develops a multilayered structure, and the eigenvalue is obtained by solving a simplified version of the Rayleigh equation (with homogeneous boundary conditions) in one of these layers which lies in either of the external streams. This analysis leads to a simple hypersonic similarity law which explains how spatial and temporal phase speeds and growth rates scale with Mach number and temperature ratio. Comparisons are made with numerical results, and it is found that this similarity law provides a good qualitative guide for the behavior of the instability at high Mach numbers. In addition to this asymptotic theory, some fully numerical results are also presented (with no limitation on the Mach number) in order to explain the origin of the hypersonic modes (through mode splitting) and to discuss the role of oblique modes over a very wide range of Mach number and temperature ratio.

Balsa, Thomas F.↗
NASA NTRS1999

The Application of a Boundary Integral Equation Method to the Prediction of Ducted Fan Engine Noise

The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) is considered. Governing equations for the BIEM are based on linearized acoustics and describe the scattering of incident sound by a thin, finite-length cylindrical duct in the presence of a uniform axial inflow. A classical boundary value problem (BVP) is derived that includes an axisymmetric, locally reacting liner on the duct interior. Using potential theory, the BVP is recast as a system of hypersingular boundary integral equations with subsidiary conditions. We describe the integral equation derivation and solution procedure in detail. The development of the computationally efficient ducted fan noise prediction program TBIEM3D, which implements the BIEM, and its utility in conducting parametric noise reduction studies are discussed. Unlike prediction methods based on spinning mode eigenfunction expansions, the BIEM does not require the decomposition of the interior acoustic field into its radial and axial components which, for the liner case, avoids the solution of a difficult complex eigenvalue problem. Numerical spectral studies are presented to illustrate the nexus between the eigenfunction expansion representation and BIEM results. We demonstrate BIEM liner capability by examining radiation patterns for several cases of practical interest.

Dunn, M. H.↗
NASA NTRS1985

Sum rules and other properties involving resonance projection operators

A sum rule is derived for the auxiliary eigenvalues of an equation whose eigenspectrum pertains to projection operators which describe electron scattering from multielectron atoms and ions. The sum rule's right-hand side depends on an integral involving the target system eigenfunctions. The sum rule is checked for several approximations of the two-electron target. It is shown that target functions which have a unit eigenvalue in their auxiliary eigenspectrum do not give rise to well-defined projection operators except through a limiting process. For Hylleraas target approximations, the auxiliary equations are shown to contain an infinite spectrum. However, using a Rayleigh-Ritz variational principle, it is shown that a comparatively simple aproximation can exhaust the sum rule to better than five significant figures. The auxiliary Hylleraas equation is greatly simplified by conversion to a square root equation containing the same eigenfunction spectrum and from which the required eigenvalues are trivially recovered by squaring.

Berk, A.↗
NASA NTRS1977

Extraction of the large eddy structure of a turbulent boundary layer

Much attention is currently being focused on turbulence modeling by separation of large-scale and small-scale turbulence. One major difficulty with this method is the problem of separating and defining the large-scale structure. In this paper, a method is presented which does this separation by analyzing time-averaged, two-point double velocity correlation measurements. The basis of the method is a maximization of the inner product of a candidate velocity field with the measured turbulent velocity field. An eigenvalue problem results, whose solution provides a set of orthogonal eigenfunctions associated with specific modes of turbulent motion. The large-scale eddy, containing most of the turbulent energy, is associated with the first (largest) mode. The turbulent velocity field is then represented as a sum of the eigenfunctions. A short description of the equations is given, the methods of solution are discussed, and some results for a two-dimensional fully developed turbulent boundary layer are presented.

Lemmerman, L. A.↗
NASA NTRS1981

Existence of axisymmetric modes of thin cylindrical shells with axial displacement restraint

A necessary and sufficient condition is derived for the existence of the axisymmetric mode of cylindrical shells with a radial displacement having one half wave along the axis and axial displacement restrained at both ends. The condition is purely geometric and consists of an upper bound on the mean radius of thickness ratio for a given fixed value of the length to mean radius ratio, above which the mode with one half wave in radial displacement along the axis ceases to exist. The proof is based on enforcing two basic lemmas concerned with the simplicity of the eigenvalues of the shell and the uniform ordering of these eigenvalues with the number of nodes of their corresponding radial displacement eigenfunctions.

El-Raheb, M.↗
NASA NTRS1984

MMIC devices for active phased array antennas

The study of printed circuit discontinuities is necessary in order to design, for example, transitions between rectangular waveguides and printed circuits. New developments with respect to the analytical approaches to this problem are discussed. A summary of the progress in the experimental approach is presented. The accurate solution for the modes in various millimeter-wave waveguides is essential in the analysis of many integrated circuit components, such as filters and impedance transformers. Problems associated with the numerical computation of these modes in two frequently used waveguide forms, namely, the finline and microstrip, are presented. The spectral domain method of formulation, with a moment method solution, is considered. This approach can be readily extended to analyze an arbitrary configuration of dielectric and metallized regions in a shielded enclosure. Galerkin's method is used, where the testing and basic functions are the same. It is shown that the mode functions, or eigenfunctions, are more sensitive to errors than the phase constants, or eigenvalues. The approximate mode functions do not satisfy the orthogonality relationship well, resulting in difficulties when these modal solutions are used to form an approximate Green's function or are used in a mode matching analysis.

Mittra, R.↗
NASA NTRS1996

Dynamics and Control of Articulated Anisotropic Timoshenko Beams

The paper illustrates the use of continuum models in control design for stabilizing flexible structures. A 6-DOF anisotropic Timoshenko beam with discrete nodes where lumped masses or actuators are located provides a sufficiently rich model to be of interest for mathematical theory as well as practical application. We develop concepts and tools to help answer engineering questions without having to resort to ad hoc heuristic ("physical") arguments or faith. In this sense the paper is more mathematically oriented than engineering papers and vice versa at the same time. For instance we make precise time-domain solutions using the theory of semigroups of operators rather than formal "inverse Laplace transforms." We show that the modes arise as eigenvalues of the generator of the semigroup, which are then related to the eigenvalues of the stiffness operator. With the feedback control, the modes are no longer orthogonal and the question naturally arises as to whether there is still a modal expansion. Here we prove that the eigenfunctions yield a biorthogonal Riesz basis and indicate the corresponding expansion. We prove mathematically that the number of eigenvalues is nonfinite, based on the theory of zeros of entire functions. We make precise the notion of asymptotic modes and indicate how to calculate them. Although limited by space, we do consider the root locus problem and show for instance that the damping at first increases as the control gain increases but starts to decrease at a critical value, and goes to zero as the gain increases without bound. The undamped oscillatory modes remain oscillatory and the rigid-body modes go over into deadbeat modes. The Timoshenko model dynamics are translated into a canonical wave equation in a Hilbert space. The solution is shown to require the use of an "energy" norm which is no more than the total energy: potential plus kinetic. We show that, under an appropriate extension of the notion of controllability, rate feedback with a collocated sensor can stabilize the structure in the sense that all modes are damped and the energy decays to zero. An example, non-numeric, is worked out in some detail illustrating the concepts and theory developed.

Balakrishnan, A. V.↗
NASA NTRS1978

Modal content of noise generated by a coaxial jet in a pipe

Noise generated by air flow through a coaxial obstruction in a long, straight pipe was investigated with concentration on the modal characteristics of the noise field inside the pipe and downstream of the restriction. Two measurement techniques were developed for separation of the noise into the acoustic duct modes. The instantaneous mode separation technique uses four microphones, equally spaced in the circumferential direction, at the same axial location. The time-averaged mode separation technique uses three microphones mounted at the same axial location. A matrix operation on time-averaged data produces the modal pressure levels. This technique requires the restrictive assumption that the acoustic modes are uncorrelated with each other. The measured modal pressure spectra were converted to modal power spectra and integrated over the frequency range 200-6000 Hz. The acoustic efficiency levels (acoustic power normalized by jet kinetic energy flow), when plotted vs. jet Mach number, showed a strong dependence on the ratio of restriction diameter to pipe diameter. The acoustic energy flow analyses based on the thermodynamic energy equation and on the results of Mohring both resulted in orthogonality properties for the eigenfunctions of the radial mode shape equation. These orthogonality relationships involve the eigenvalues and derivatives of the radial mode shape functions.

Kerschen, E. J.↗
NASA NTRS1986

Wave-packet evolution and quantization

A general relationship is established between wave-packet evolution and quantization in classically integrable systems. Because of wave-packet spreading, one cannot simply take the Fourier transform of the time evolution of a wave-packet. Instead, one must propagate the wave-packet using the actions as Hamiltonians. The energy eigenvalues which result are the Einstein-Brillouin-Keller values, and new forms for the eigenfunctions appear. These are free of caustic singularities, and represent averages of wave packets over the invariant torus.

Littlejohn, R. G.↗
NASA NTRS1987

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗
NASA NTRS1989

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗
NASA NTRS1987

The characteristic Rossby frequency

The characteristic Rossby frequency is defined for a fixed zonal wavenumber perturbation as the variational integral of the Rayleigh-Ritz method. It is a measure of the time scale of the disturbance. For a disturbance which locally has the shape of an eigenfunction but is not global in extent, the characteristic Rossby frequency is very close to the true eigenvalue, and additionally remains unchanged under linear inviscid dynamics. Results are presented for the shallow water equations, both with and without a mean zonal wind. The characteristic Rossby frequency of a wavenumber 1 perturbation having the shape of the second symmetric Rossby mode but confined to the Northern Hemisphere is close to the corresponding Rossby frequency. This finding is helpful in understanding the behavior of the observed wavenumber 1 pattern of January 1979, which propagated westward with nearly the pure Rossby frequency but was discernible only in the Northern Hemisphere (as discussed by Daley and Williamson).

Straus, David M.↗
NASA NTRS1986

Eigenfunction methods in magnetospheric radial-diffusion theory

Complete sets of orthonormal basis functions constructed according to a generalization of the quantum-mechanical WKB approximation can be used to generate a nearly-diagonal matrix representation of the radial-transport operator for ring-current ions in the presence of radial diffusion and charge exchange. The resulting eigenfunctions (constructed by weighting the basis functions in proportion to the respective components of the eigenvectors of the matrix representation) and eigenvalues provide a spatial and temporal description of the evolving phase-space density during and following a magnetospheric disturbance (e.g., a magnetic storm). A linear superposition of the basis functions can also be used to eliminate any discrepancy between the steady-state solution of the transport equation and the appropriate WKB approximation of this steady-state solution.

Schulz, Michael↗
NASA NTRS1985

Numerical Techniques for Scattering from Submerged Objects

To represent the final results in terms of matrices, one expands all appropriate physical quantities in terms of partial wave basis states. This includes expansions for the incident and scattered fields and the surface quantities. The method then utilizes the Huygen-Poincare integral representation for both the exterior and interior solutions, leading to the required matrix equations. One thus deals with matrix equations, the complexity of which depends on the nature of the problem. It is shown that in general a transition matrix T can be obtained relating the incident field A with the scattered field f having the form T = PQ(-1), where f = TA. The structure of Q can be quite complicated and can itself be composed of other matrix inversions such as arise from layered objects. Recent improvements in this method appropriate for a variety of physical problems are focused on, and on their implementation. Results are outlined from scattering simulations for very elongated submerged objects and resonance scattering from elastic solids and shells. The final improvement concerns eigenfunction expansions of surface terms, arising from solution of the interior problem, obtained via a preconditioning technique. This effectively reduces the problem to that of obtaining eigenvalues of a Hermitian operator. This formalism is reviewed for scattering from targets that are rigid, sound-soft, acoustic, elastic solids, elastic shells, and elastic layered objects. Two sets of the more interesting results are presented. The first concerns scattering from elongated objects, and the second to thin elastic spheroids.

Werby, M. F.↗