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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Algorithm for Stabilizing a POD-Based Dynamical System

This algorithm provides a new way to improve the accuracy and asymptotic behavior of a low-dimensional system based on the proper orthogonal decomposition (POD). Given a data set representing the evolution of a system of partial differential equations (PDEs), such as the Navier-Stokes equations for incompressible flow, one may obtain a low-dimensional model in the form of ordinary differential equations (ODEs) that should model the dynamics of the flow. Temporal sampling of the direct numerical simulation of the PDEs produces a spatial time series. The POD extracts the temporal and spatial eigenfunctions of this data set. Truncated to retain only the most energetic modes followed by Galerkin projection of these modes onto the PDEs obtains a dynamical system of ordinary differential equations for the time-dependent behavior of the flow. In practice, the steps leading to this system of ODEs entail numerically computing first-order derivatives of the mean data field and the eigenfunctions, and the computation of many inner products. This is far from a perfect process, and often results in the lack of long-term stability of the system and incorrect asymptotic behavior of the model. This algorithm describes a new stabilization method that utilizes the temporal eigenfunctions to derive correction terms for the coefficients of the dynamical system to significantly reduce these errors.

Kalb, Virginia L.↗

DOE OSTI2024

Weighted Composition Operators for Learning Nonlinear Dynamics

Operator theoretic methods in dynamical system have been dominated by the use of Koopman operators and their continuous time counterparts, such as Koopman Generators and Liouville Operators. The advantage gained from their use primarily stems from the ability to extract subspaces and eigenfunctions within a space of observables that are invariant with respect to the Koopman operator over that space. When this occurs, a dynamic mode decomposition of the systems state provides a linear model for the dynamical system. Not all Koopman operators have eigenfunctions that may be exploited in this manner. However, the framework can still be leveraged for approximations using other operators. In this setting, we present a different operator for the study of dynamical systems, the weighted composition operator. These operators are compact for a wide range of dynamics and spaces, and through their interactions with occupation kernels and vector valued kernels, they admit an estimation of the underlying dynamics. Here, this manuscript presents a new algorithm for the data driven study of dynamical systems from data, and also provides two numerical experiments where convergence is achieved as a proof of concept.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

A technique for designing active control systems for astronomical telescope mirrors

This paper considers the problem of designing a control system to achieve and maintain the required surface accuracy of the primary mirror of a large space telescope. Control over the mirror surface is obtained through the application of a corrective force distribution by actuators located on the rear surface of the mirror. The design procedure is an extension of a model control technique developed for distributed parameter plants with known eigenfunctions to include plants whose eigenfunctions must be approximated by numerical techniques. Instructions are given for constructing the mathematical model of the system, and a design procedure is developed for use with. typical numerical data in selecting the number and location of the actuators. Examples of actuator patterns and their effect on various errors are given.

Active optics↗

NASA NTRS1974

Electromagnetic and scalar diffraction by a right-angled wedge with a uniform surface impedance

The diffraction of an electromagnetic wave by a perfectly-conducting right-angled wedge with one surface covered by a dielectric slab or absorber is considered. The effect of the coated surface is approximated by a uniform surface impedance. The solution of the normally incident electromagnetic problem is facilitated by introducing two scalar fields which satisfy a mixed boundary condition on one surface of the wedge and a Neumann of Dirichlet boundary condition on the other. A functional transformation is employed to simplify the boundary conditions so that eigenfunction expansions can be obtained for the resulting Green's functions. The eigenfunction expansions are transformed into the integral representations which then are evaluated asymptotically by the modified Pauli-Clemmow method of steepest descent. A far zone approximation is made to obtain the scattered field from which the diffraction coefficient is found for scalar plane, cylindrical or sperical wave incident on the edge. With the introduction of a ray-fixed coordinate system, the dyadic diffraction coefficient for plane or cylindrical EM waves normally indicent on the edge is reduced to the sum of two dyads which can be written alternatively as a 2 X 2 diagonal matrix.

Hwang, Y. M.↗

NASA NTRS1973

Variational calculations for resonance oscillations of inhomogeneous plasmas

The electrostatic resonance properties of an inhomogeneous plasma column are reported by application of the Rayleigh-Ritz method. A description of the rf equation of motion and pressure term that expresses the system of equations in Euler-Lagrange form is presented. The Rayleigh-Ritz procedure is applied to the corresponding Lagrangian to obtain approximate resonance frequencies and eigenfunctions. An appropriate set of trial coordinate functions is defined, which leads to frequency and eigenfunction estimates.

Peng, Y. K. M.↗

NASA NTRS1974

Three-body Coulomb systems using generalized angular-momentum S states

An expansion of the three-body Coulomb potential in generalized angular-momentum eigenfunctions developed earlier by one of the authors is used to compute energy eigenvalues and eigenfunctions of bound S states of three-body Coulomb systems. The results for He, H(-), e(-)e(+)e(-), and pmu(-)p are compared with the results of other computational approaches.

Whitten, R. C.↗

NASA NTRS1975

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗

NASA NTRS1975

Variational calculations for resonance oscillations of inhomogeneous plasmas

The electrostatic resonance properties of an inhomogeneous plasma column are treated by the Rayleigh-Ritz method. In contrast to Parker, Nickel & Gould (1964), who carried out an exact computation, the present treatment uses a description of the RF equation of motion and pressure term that allows one to express the system of equations in Euler-Lagrange form. The Rayleigh-Ritz procedure is then applied to the corresponding Lagrangian, to obtain approximate resonance frequences and eigenfunctions. An appropriate set of trial coordinate functions is defined, which leads to frequency and eigenfunction estimates in excellent agreement with the work of Parker et al. (1964).

Peng, Y.-K. M.↗

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 NTRS1977

Theoretical intensity parameters for the vibration-rotation bands of ClO

Intensity parameters for the vibration-rotation bands of ClO were calculated by evaluating matrix elements of the vibrational eigenfunctions using the theoretical dipole moment curve. These vibrational eigenfunctions are based on realistic Klein-Dunham potentials deduced from experimental values (Coxon et al, 1976) for certain vibrational energy levels. The data can be applied to the interpretation of remote or in situ experiments to deduce stratospheric ClO concentrations from vibration-rotation absorption spectra.

Langhoff, S. R.↗

NASA NTRS1983

Elasticity solutions for a class of composite laminate problems with stress singularities

A study on the fundamental mechanics of fiber-reinforced composite laminates with stress singularities is presented. Based on the theory of anisotropic elasticity and Lekhnitskii's complex-variable stress potentials, a system of coupled governing partial differential equations are established. An eigenfunction expansion method is introduced to determine the orders of stress singularities in composite laminates with various geometric configurations and material systems. Complete elasticity solutions are obtained for this class of singular composite laminate mechanics problems. Homogeneous solutions in eigenfunction series and particular solutions in polynomials are presented for several cases of interest. Three examples are given to illustrate the method of approach and the basic nature of the singular laminate elasticity solutions. The first problem is the well-known laminate free-edge stress problem, which has a rather weak stress singularity. The second problem is the important composite delamination problem, which has a strong crack-tip stress singularity. The third problem is the commonly encountered bonded composite joints, which has a complex solution structure with moderate orders of stress singularities.

Wang, S. S.↗

NASA NTRS1985

Numerical Descriptions of Cosmic-Ray Transport

The behavior of energetic particles in the solar system is described by a well known Fokker-Planck equation. Although analytic methods yield insight into the nature of its solutions, especially in the diffusion regime, calculations that go beyond diffusion are very complicated. The reliability of these calculations is of concern, because numerical methods are notorious for their errors and artifacts. The well known Milne problem of classical transport theory was analyzed with the aid of three different numerical methods. These are: (1) The method of eigenfunctions in which the distribution function is approximated by a sum of eigenfunctions of the scattering operator, (2) Numerical solutions of a finite difference aquation; and (3) Direct simulation of the scattering and streaming of individual particles with the aid of Monte Carlo methods.

Earl, J. A.↗

NASA NTRS1985

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, 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 NTRS1985

Application of satellite data to tropic/subtropic moisture coupling

The objective is to develop analysis tools for use of satellite data to interpret synoptic-scale systems in data-void regions. Interim goals are to: (1) quantify the synoptic information content of satellite data; and (2) utilize these data in the diagnosis of moisture bursts in the eastern tropical Pacific Ocean. Researchers developed and implemented a statistical procedure for using TIROS N microwave data to infer infrared channel data for overcast conditions; they used the same procedure for deducing full TIROS N channel radiance profiles from NOAA 5 VTPR channel data over regions where the TIROS N data are missing. An empirical orthogonal function analysis of twice-daily channel radiance fields over the tropical eastern Pacific was completed. The vertically oriented eigenfunctions were interpreted in terms of typical meteorological events. The horizontal distribution of the eigenfunction amplitudes relates these meteorological signals to moisture bursts. A pair of moisture burst climatologies is complete: one of four years using infrared imagery (including the highly anomalous 1982 to 83 cold season); the other implementing 850 to 200 mb wind analyses in conjunction with GOES imagery. A number of different evaluations of the synoptic evolution of moisture fields (enhanced infrared imagery, moisture channel data, FGGE humidity analysis, and in situ station and sounding observations) are compared. All have limitations; all can be utilized together; all together are still less than adequate in the tropical Pacific.

Mcguirk, J. P.↗