Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 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↗
NASA NTRS1999

Development of a Three-Dimensional PSE Code for Compressible Flows: Stability of Three-Dimensional Compressible Boundary Layers

A program is developed to investigate the linear stability of three-dimensional compressible boundary layer flows over bodies of revolutions. The problem is formulated as a two dimensional (2D) eigenvalue problem incorporating the meanflow variations in the normal and azimuthal directions. Normal mode solutions are sought in the whole plane rather than in a line normal to the wall as is done in the classical one dimensional (1D) stability theory. The stability characteristics of a supersonic boundary layer over a sharp cone with 50 half-angle at 2 degrees angle of attack is investigated. The 1D eigenvalue computations showed that the most amplified disturbances occur around x(sub 2) = 90 degrees and the azimuthal mode number for the most amplified disturbances range between m = -30 to -40. The frequencies of the most amplified waves are smaller in the middle region where the crossflow dominates the instability than the most amplified frequencies near the windward and leeward planes. The 2D eigenvalue computations showed that due to the variations in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the most amplified disturbances are shifted to 120 degrees compared to 90 degrees for the parallel theory. It is also observed that the nonparallel amplification rates are smaller than that is obtained from the parallel theory.

Balakumar, P.↗
NASA NTRS1997

Analysis, Repair, and Management of the Total Ozone Mapping Spectrometer Database

In the ensuing period we were able to demonstrate that the origin of these filamentous patterns resulted from the action of synoptic-scale vortical velocity field on the global-scale background gradient of ozone concentration in the meridional direction. Hyperbolic flow patterns between long-lived atmospheric vortices bring together air parcels from different latitudes, thus creating large gradients along the separatrices leaving the hyperbolic (stagnation) point. This result is further confirmed by the KL analysis of the ozone field in the equatorial region, where the background concentration gradient vanishes. The spectral slope in this region has been found to lie close to -1, in agreement with Batchelor's prediction. Another outcome of this result is that it at least provides indirect evidence about the kinetic energy spectrum of the atmospheric turbulence in the range of scales approximately 200 to 2000 km. Namely, Batchelor's analysis is based on the assumption that the velocity field is large-scale, that is the kinetic energy spectrum decays as O(k(sup -3)) or steeper. Since the scalar spectrum is confirmed, this also supports this form of the kinetic energy spectrum. The study of equatorial regions of TOMS data revealed the efficiency of the KL method is in detecting and separating a wave-like measurement artifact inherently present in the dataset due to the non-perfect correction for cross-track bias. Just two to three eigenfunctions represent the error, which makes it possible to enhance the data by reconstituting it from the data by eliminating the subspace of artifactual eigenfunctions. This represents a highly efficient means for achieving an improved rendering of the data. This has been implemented on the database. A wide range of techniques and algorithms have been developed for the repair and extension of the TOMS database.

Sirovich, Lawrence↗
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 NTRSDate not supplied

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗