Search NASA⌕ Search

SEARCH · Search NASA

Results for “EIGENFUNCTION”

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 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 NTRS1986

Methods for the solution of radiative transfer equation

To obtain an exact solution of the radiative-transfer equation in media where both absorption and scattering are significant, the usual approach is to use a numerical method. Three methods are known in the literature: invariant imbedding, eigenvalue-eigenfunction, and matrix doubling. This paper examines the practical application of these methods to the problem of emission from an inhomogeneous (Rayleigh) layer, the effects of layer parameters on the stability. It is found that invariant imbedding is most suitable for computing emission from an inhomogeneous layer with a temperature profile but tends to be unstable as the optical thickness of the layer increases beyond 0.5. On the other hand, the matrix-doubling method is stable for arbitrary optical thickness but is not suitable for handling multilayers. The eigenvalue-eigenfunction method is more stable than the invariant imbedding as optical thickness increases up to 2.0. It also permits temperature profile in the layer, but the computation is much more complicated. It is less stable than the matrix-doubling method when optical thickness is larger than 2.0. In general, the choice of a method is dependent on the nature of the problem.

Chen, M. F.↗
NASA NTRS1987

A linear thermal stability analysis for the vertical structure of alpha model accretion disks

The linear equations that describe the thermal stability of the convective and radiative vertical structure of an alpha model accretion disk are derived. For a specific model chosen to be representative of a dwarf nova disk, it is found that the eigenfunction for the temperature shows large excursions in the zone of partial ionization, but that the eigenfunctions for the unstable fundamental mode perturbations of the vertical distance and flux are nearly constant over the vertical extent of the disk. The growth of any instability is thus nearly homologous and the stability criteria derived for vertically averaged structures represent an adequate approximation.

Saio, Hideyuki↗
NASA NTRS1987

Solar oscillations - A method for deriving nonlinear effects

The frequencies of solar oscillations are so closely spaced that nonlinear interactions among modes are probable. The rate of interaction is proportional to an integral involving the eigenfunctions of the interacting modes, which usually are known only numerically. An approximation in which the eigenfunctions are strictly sinusoidal functions of a suitably defined radial variable, and the numerical details of stellar structure are banished to a coefficient in the integrand, are here explored. The physical assumptions are the same as in the asymptotic approximation of p- or g-modes. This method should allow a general investigation as to likely nonlinear interactions and which modes may participate in such interactions. The coupling of two p-modes by possible large-scale internal magnetic fields is introduced as an anisotropic pressure response to a displacement. Pairs of modes differing in frequency by less than the fraction magnetic/thermal pressure are strongly coupled, and energy appears to oscillate slowly between the two associated spherical harmonics. Potentially, upper limits may be derived for internal magnetic fields.

Wentzel, Donat G.↗
NASA NTRS1988

Sampling inhomogeneous turbulent fields

The reconstruction of an inhomogeneous random process from a finite number of discrete samples can be performed in terms of the Karhunen-Loeve (KL) expansion for that process. The n(th) eigenfunction has n - 1 zero crossings which are the sampling points for the inhomogeneous process. The rapid variation of the KL eigenfunctions makes it unnecessary to have a high density of sampling (or grid points) near the wall. However, this result should not be construed to indicate that with spectral simulations significantly fewer grid points are required with the KL expansion as compared to other orthogonal expansions. Moin and Moser (1989) have shown that the advantage of the KL expansion over Chebychev expansion rapidly diminishes when high percentage (say 90 percent) energy recovery is demanded.

Adrian, R. J.↗
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 NTRS1989

The normal modes of a resonant cavity containing discrete inhomogeneities - The influence of fibril magnetic fields on the solar acoustic oscillations

Motivated by considerations of the interaction between fibril magnetic fields and solar p-modes, the acoustic spectrum of a cylindrical cavity filled with ideal gas in which a number of magnetic flux tubes are embedded is studied. A formalism, based on the T-matrix approach to acoustic scattering, is developed which can be used to determine the eigenfrequencies and eigenfunctions for any arbitrary distribution of flux tubes. For weak scatterers, the frequency shifts and velocity eigenfunctions are calculated using perturbation theory for the cases of a single flux tube and a random distribution of up to 100 flux tubes. The results of this 'exact' approach are used to give a critical appraisal of the predictions of theories based on some form of averaging, such as the one discussed recently by Bogdan and Zweibel (1987).

Bogdan, Thomas J.↗
NASA NTRS1991

The inviscid compressible Goertler problem

The growth rate is studied of Goertler vortices in a compressible flow in the inviscid limit of large Goertler number. Numerical solutions are obtained for 0(1) wavenumbers. The further limits of large Mach number and large wavenumber with 0(1) Mach number are considered. It is shown that two different types of disturbance modes can appear in this problem. The first is a wall layer mode, so named as it has its eigenfunctions trapped in a thin layer away from the wall and termed a trapped layer mode for large wavenumbers and an adjustment layer mode for large Mach numbers, since then this mode has its eigenfunctions concentrated in the temperature adjustment layer. The near crossing of the modes which occurs in each of the limits mentioned is investigated.

Dando, Andrew↗
NASA NTRS1992

Diskoseismology: Probing accretion disks. II - G-modes, gravitational radiation reaction, and viscosity

A scalar potential is used to derive a single partial differential equation governing the oscillation of a disk. The eigenfunctions and eigenfrequencies of a variety of disk models are found to fall into two main classes which are analogous to the p-modes and g-modes in the sun. Specifically, the eigenfunctions and eigenfrequencies of isothermal disks are computed, and the way in which these results can be generalized to other disk models is indicated. The (assumed) relatively small rates of growth or damping of the modes due to various mechanisms, in particular gravitational radiation reaction and parameterized models of viscosity are also computed. It is found that for certain parameters the p-modes are unstable to gravitational radiation reaction (CFS instability), while both the p-modes and g-modes are unstable to viscosity unless highly anisotropic viscosity models are considered.

Nowak, Michael A.↗
NASA NTRS1992

Diskoseismology - Signatures of black hole accretion disks

General relativity requires the existence of a spectrum of oscillations which are trapped near the inner edge of accretion disks around black holes. We have developed a general formalism for analyzing the normal modes of such acoustic perturbations of arbitrary thin disk models, approximating the dominant relativistic effects via a modified Newtonian potential (these modes do not exist in Newtonian gravity). The eigenfunctions and eigenfrequencies of a variety of disk models are found to fall in to two main classes, which are analogous to the p-modes and g-modes in the sun. In this work, we compute the eigenfunctions and eigenfrequencies of isothermal disks. The (relatively small) rates of growth or damping of these oscillations due to gravitational radiation and parameterized models of viscosity are also computed.

Nowak, Michael↗
NASA NTRS1993

Instability of rectangular jets

The instability of rectangular jets is investigated using a vortex-sheet model. It is shown that such jets support four linearly independent families of instability waves. Within each family there are infinitely many modes. A way to classify these modes according to the characteristics of their mode shapes or eigenfunctions is proposed. It is demonstrated that the boundary element method can be used to calculate the dispersion relations and eigenfunctions of these instability wave modes. The method is robust and efficient. A parametric study of the instability wave characteristics has been carried out. A sample of the numerical results is reported here. It is found that the first and third modes of each instability wave family are corner modes. The pressure fluctuations associated with these instability waves are localized near the corners of the jet. The second mode, however, is a center mode with maximum fluctuations concentrated in the central portion of the jet flow. The center mode has the largest spatial growth rate. It is anticipated that as the instability waves propagate downstream the center mode would emerge as the dominant instability of the jet.

Tam, Christopher K. W.↗
NASA NTRS1992

The form of spin orbitals for open-shell restricted Hartree-Fock reference functions

A new set of spin orbitals is proposed for single-reference correlation procedures based on an open-shell restricted Hartree-Fock (RHF) reference function. The spatial parts of the new spin orbitals are from an open-shell RHF calculation. These spin orbitals form a complete, orthonormal set. The RHF single determinant based on these spin orbitals is an S sub 2 eigenfunction, but a linear combination of S sub z eigenfunctions. The advantage of these spin orbitals is that equations for perturbation theory or coupled-cluster theory are symmetric in the spin indices and hence the number of independent parameters to be determined is dramatically reduced.

Jayatilaka, Dylan↗
NASA NTRS1994

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗
NASA NTRS1995

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗
NASA NTRS1994

Gravitational instabilities in a proto-planetary disk including the effects of magnetic fields

We investigate the gravitational instability of a thin, Keplerian protoplanetary disk including the effects of a largely azimuthal magnetic field. The model follows that of our previous work (Noh, Vishniac, & Cochran 1991) except for the inclusion of a magnetic field. The disk is assumed to consist of neutral and ionized gas and neutral dust which are coupled by gravity and friction. The growth rates and eigenfunctions are calculated numerically using nonaxisymmetric linear perturbation methods. The results show that the growth rate has a maximum at some intermediate azimuthal number m, but for each value of m it is reduced relative to the unmagnetized case. The effects of the magnetic field appear more strongly on small scales. As the strength of the equilibrium magnetic field increases the growth rates decrease, and the maximum instability occurs at a lower value on m due to the increasing magnetic pressure. The response of each component to the magnetic field is discussed using the behavior of the eigenfunctions in the radial direction. With the inclusion of the magnetic field, the effects of the ionization fraction and friction on the growth rates also appear to be important for high m modes. Increasing the ionization fraction or the friction suppresses instability, but only slightly changes the maximally unstable azimuthal scales. The enhanced growth rates due to a dust component for which thermal pressure is negligible are somewhat reduced by the inclusion of a magnetic field. The effects of different boundary conditions (reflecting and transmitting) on the growth rates are also shown.

Noh, Hyerim↗
NASA NTRS1994

Low-dimensional description of the dynamics in separated flow past thick airfoils

Results are presented for the numerical simulation of unsteady viscous incompressible flow past thick airfoils. Specifically, flow past a NACA 4424 at an angle of attack of 2.5 deg and Reynolds numbers in the range of 1700-4000 has been simulated using the spectral element method. At these conditions the flow is separatedd and an unsteady wake is formed. Application of the method of empirical eigenfunction reveals the structure of the most energetic components of the flow. These are found to occur in pairs that, through phase exchange, are responsible for the vortex shedding. A set of ordinary differential equations is obtained for the amplitudes of these eigenfunctions by a Galerkin projection of the Navier-Stokes equations. The solutions of the model system are compared with the full simulation. The work is of relevance to the transition process and observed routes to chaos in airfoil wakes.

Deane, Anil E.↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
NASA NTRS1998

A High Frequency Model of Cascade Noise

Closed form asymptotic expressions for computing high frequency noise generated by an annular cascade in an infinite duct containing a uniform flow are presented. There are two new elements in this work. First, the annular duct mode representation does not rely on the often-used Bessel function expansion resulting in simpler expressions for both the radial eigenvalues and eigenfunctions of the duct. In particular, the new representation provides an explicit approximate formula for the radial eigenvalues obviating the need for solutions of the transcendental annular duct eigenvalue equation. Also, the radial eigenfunctions are represented in terms of exponentials eliminating the numerical problems associated with generating the Bessel functions on a computer. The second new element is the construction of an unsteady response model for an annular cascade. The new construction satisfies the boundary conditions on both the cascade and duct walls simultaneously adding a new level of realism to the noise calculations. Preliminary results which demonstrate the effectiveness of the new elements are presented. A discussion of the utility of the asymptotic formulas for calculating cascade discrete tone as well as broadband noise is also included.

Envia, Edmane↗