SEARCH · Search NASA
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.
An extensive examination of ellipsometric equations for anisotropic surfaces and films is reported. It is shown that the reflection matrix can be calculated by writing Maxwell's equations in 6 x 6 matrix form and by applying appropriate boundary conditions at proper points in the development of the formalism, reducing the equation for the propagation of light to an eigenvalue problem using a 4 x 4 matrix.
Two methods for natural mode vibration analysis are discussed. The first consists of a direct approach based on a finite element representation of the complete structure as an entity. The mass and stiffness matrices for the complete structure are assembled by properly combining the mass and stiffness matrices of the individual elements into which the structure has been divided. The second approach is that of component mode synthesis. This method is based on the concept of synthesizing the natural modes of the complete structure from modes of conveniently difined substructures, or components, into which the structure has been partitioned. In this way the expedient of reducing the system degrees of freedom, and thus the size of the eigenvalue problem, can be introduced by partial modal synthesis.
The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.
Quantum mechanical calculations require the repeated solution of a Schrödinger equation for the wavefunctions of the system, from which materials properties follow. Recent work has shown the effectiveness of enriched finite element type Galerkin methods at significantly reducing the degrees of freedom required to obtain accurate solutions. However, time to solution has been adversely affected by the need to solve a generalized rather than standard eigenvalue problem and the ill-conditioning of associated system matrices. In this work, we address both issues by proposing a stable and efficient orbital-enriched partition of unity method to solve the Schrödinger boundary-value problem in a parallelepiped unit cell subject to Bloch-periodic boundary conditions. In the proposed partition of unity method, the three-dimensional domain is covered by overlapping patches, with a compactly-supported weight function associated with each patch. A key ingredient in our approach is the use of non-negative weight functions that possess the flat-top property, i.e., each weight function is identically equal to unity over some finite subset of its support. This flat-top property provides a pathway to devise a stable approximation over the whole domain. On each patch, we use $p$th degree orthogonal (Legendre) polynomials that ensure $p$th order completeness, and in addition include eigenfunctions of the radial Schrödinger equation. Furthermore, we adopt a variational lumping approach to construct a (block-)diagonal overlap matrix that yields a standard eigenvalue problem for which there exist efficient eigensolvers. The accuracy, stability, and efficiency of the proposed method is demonstrated for the Schrödinger equation with a harmonic potential as well as a localized Gaussian potential. In conclusion, we show that the proposed approach delivers optimal rates of convergence in the energy, and the use of orbital enrichment significantly reduces the number of degrees of freedom for a given desired accuracy in the energy eigenvalues while the stability of the enriched approach is fully maintained.
A general approach to analyzing eigenmodes in anisotropic and gyrotropic 3D photonic crystals based on dielectric and magnetic media is proposed. This approach is based on the representation of stationary macroscopic Maxwell equations in the operator form corresponding to the quantum-mechanical equation for a photon with spin s = 1. In these equations, the strengths of electric and magnetic fields are put into correspondence with state vectors in the complex Hilbert space. The permittivity and permeability serve as operators acting on these vectors. It is shown that the problem of determining the eigenmodes of a photonic crystal is reduced to searching for eigenvectors and eigenvalues of the Hermitian operator characterizing the spin–orbit interaction of a photon in the periodic anisotropic medium under study. It is proposed to use photon states with a certain wave vector (certain momentum) and a certain linear or circular spin polarization as a basis for the representation operator equations. One-dimensional photonic crystals are considered as an example. The influence of anisotropy and gyrotropy on the dispersion of eigenmodes in these crystals is investigated. The group velocity of the eigenmodes, momentum transferred by them, and spin angular momentum are analyzed for the case of gyrotropic media.
The structure-control system optimization problem is formulated with constraints on the closed-loop eigenvalues and the efficiency of the reduced order system. The feasibility of the approach is illustrated by designing the ACOSS-FOUR structure with a reduced order system and improving the efficiency characteristics of the structures-control system.
A very simple modification is presented for the Chebyshev tau method which can eliminate spurious eigenvalues, proceeding from a consideration of the vorticity-streamfunction reformulation of the Chebyshev tau method and the Chebyshev-Galerkin method, which have no spurious modes. Consideration of a model problem indicates that these two approaches are equivalent, and that they reduce to the present modification of the tau method. This modified tau method also eliminates spurious eigenvalues from the Orr-Sommerfeld equation.
The stability of a self-gravitating, nonrotating, isothermal gas layer threaded by a one-dimensional equipartition magnetic field, immersed in a rigid isothermal layer of stars, is considered with respect to waves with motions in and perpendicular to the X plane, where X is the equilibrium magnetic and gravitational field vectors. When motions are perpendicular to the X plane, the magnetic field hinders gravitational instability, increasing the minimum length necessary to produce instability by the factor (1 + alpha) to the 1/2 power, where alpha is the ratio of magnetic pressure to gas pressure. When motions are in the plane, no such simple analytical solution is found. However, the resulting system of equations reduces to a single fourth order differential equation in the perturbed gas potential that defines an eigenvalue problem for the marginally unstable state when the four appropriate boundary conditions are considered.
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.
The Subsonic Ultra Green Aircraft Research (SUGAR) Phase III was led by Dr. Rakesh K. Kapania and Dr. Joseph A. Schetz at the Multidisciplinary Analysis and Design Center for Advanced Vehicles, Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg VA. The research was performed from December 2014 to December 2015. Three major areas were investigated: Multidisciplinary Design Optimization (MDO) studies of truss braced wing (TBW) and strut braced wing (SBW) vehicles at cruise Mach numbers of 0.7 and 0.8 for a flight mission similar to current market single aisle configurations. The performance and the characteristics of the optimized vehicles were compared to the SUGAR Phase II TBW vehicle. These results were obtained without applying any of the extended transonic aerodynamic and aeroelastic tools that will be discussed later. It was found that the cruise Mach number has a large effect on the best “truss” configuration. At Mach 0.7, an SBW has a better fuel consumption and better take-off gross weight (TOGW). However, at Mach 0.8, the TBW is superior because the jury strut aids in satisfying the flutter constraint; Two-dimensional, steady, transonic aerodynamic analysis of the Boeing Airfoil J (BACJ) airfoil was performed for a range of thickness ratios, Mach numbers and lift coefficients. Reynolds-averaged Navier-Stokes (RANS) equations were solved to obtain the lift-curve slope, wave drag coefficient, the location of the center of pressure and to predict the separation at the trailing edge, which may lead to buffeting. One of the goals was to develop a database of lift-curve slope and the location of center of pressure, which could be used in a transonic aeroelastic analysis. Another goal was to compare the wave drag coefficients to those predicted by Lock’s fourth-power law and also to compare the transonic effects obtained from RANS simulations to those predicted by the Korn equations. A third goal was to develop a buffet boundary, which can be integrated into the MDO framework to prevent the optimized designs from probable buffeting; A state-space transonic aeroelastic analysis tool was developed, which can incorporate the nonlinear transonic effects in the unsteady aerodynamics but is yet computationally cheap when used within the VT MDO framework. The aeroelastic analysis uses Leishman- Beddoes (LB) indicial functions, which generated a state-space representation of the aeroelastic system. The indicial functions allow the incorporation of data for steady lift-curve slope and location of the center of pressure. Thus, the steady transonic effects are included, and the unsteady aerodynamic responses are a linearization about the steady results. The aeroelastic approach discretizes the wing into numerous strips, which results in a large eigenvalue problem as each strip has eight augmented aerodynamic states as per the LB theory. Thus, to reduce the computation expense, a reduced order model (ROM) was developed. The approach was validated using a few examples.
An extension of the Tridiagonal Reduction (FEER) method to complex eigenvalue analysis in NASTRAN is described. As in the case of real eigenvalue analysis, the eigensolutions closest to a selected point in the eigenspectrum are extracted from a reduced, symmetric, tridiagonal eigenmatrix whose order is much lower than that of the full size problem. The reduction process is effected automatically, and thus avoids the arbitrary lumping of masses and other physical quantities at selected grid points. The statement of the algebraic eigenvalue problem admits mass, damping and stiffness matrices which are unrestricted in character, i.e., they may be real, complex, symmetric or unsymmetric, singular or non-singular.
Partial eigensystem assignment with output feedback can lead to an unstable closed-loop system. However, output feedback with passive linear time-invariant systems, such as flexible space structures, is guaranteed to be stable if the controller is dissipative. This paper presents a novel approach for synthesis of dissipative output feedback gain matrices for assigning a selected number of closed-loop poles. Dissipativity of a gain matrix is known to be equivalent to positive semidefiniteness of the symmetric part of the matrix. A sequential procedure is presented to assign one self-conjugate pair of closed-loop eigenvalues at each step using dissipative output feedback gain matrices, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The problem of assigning one closed-loop pair is reduced to a constrained solution of a system of quadratic equations, and necessary and sufficient conditions for the existence of a solution are presented. A minimax approach is presented for determining parameters which satisfy these conditions. This method can assign as many closed-loop system poles as the number of control inputs. A numerical example of damping enhancement for a flexible structure is presented to demonstrate the approach.
This paper presented a novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once, the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices, which reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.
A NASTRAN analysis of the external tank (ET) substructure of the 1/8-scale space shuttle structural dynamics model is described. The NASTRAN hydroelastic procedures were used to form a model of the liquid oxygen portion of the ET. Large computer storage requirements and running times were required unless these procedures were modified. Several possibilities were demonstrated including the substitution of the real for the complex eigenvalue routine and the use of the OMIT capability to reduce the number of fluid coordinates. A NASTRAN model of the complete ET was formed and reduced to 252 degrees of freedom. A review of the eigenvector extracted, using the unsymmetrical inverse power method, indicated that the structural OMIT's resulted in unsatisfactory modal deflections. Modifications to the NASTRAN hydroelastic capability appeared necessary to generate a successful ET model. An approach is described which, by assuming incompressibility, reduces the fluid representation to a symmetric mass matrix which can be added to the structural mass. The problem is solved using faster and more efficient eigenvalue routines. Using this approach the ET NASTRAN model was analyzed for three separate weight conditions.
A study was made of certain special cases of a generalized eigenvalue problem. Let A and B be nxn matrics. One may construct a certain polynomial, P(A,B, lambda) which specializes to the characteristic polynomial of B when A equals I. In particular, when B is hermitian, that characteristic polynomial, P(I,B, lambda) has real roots, and one can ask: are the roots of P(A,B, lambda) real when B is hermitian. We consider the case where A is positive definite and show that when N equals 3, the roots are indeed real. The basic tools needed in the proof are Shur's theorem on majorization for eigenvalues of hermitian matrices and the interlacing theorem for the eigenvalues of a positive definite hermitian matrix and one of its principal (n-1)x(n-1) minors. The method of proof first reduces the general problem to one where the diagonal of B has a certain structure: either diag (B) = diag (1,1,1) or diag (1,1,-1), or else the 2 x 2 principal minors of B are all 1. According as B has one of these three structures, we use an appropriate method to replace A by a positive diagonal matrix. Since it can be easily verified that P(D,B, lambda) has real roots, the result follows. For other configurations of B, a scaling and a continuity argument are used to prove the result in general.
A new design procedure for a class of multivariable integral plus proportional servo control systems is presented. The procedure utilizes a null space formulation of the eigensystem assignment problem assuming complete state feedback. Reduced order observers are used to estimate unavailable states. The resulting system realizes arbitrary closed-loop eigenvalues, approximates specified modal behavior and achieves low eigensystem sensitivity with respect to plant parameter variations. A tutorial example to illustrate the design procedure is included.
The method of simultaneous iteration with shift is extended to extraction of m-eigenpairs of a general eigenvalue problem of large order n in a parallel processing environment. The algorithm combines the power method and the Jacobi technique, and reduces to performing four basic operations. Parallel implementation of the algorithm is discussed in detail. The analysis accounts for computation and communication costs, and utilizes a parallel processing architecture of the ensemble type. Expressions for the computational efficiency and speedup are defined as a function of the problem and hardware parameters. Selected representative problems exhibit efficiencies ranging from 60 to 98 percent.
IA novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices with appropriate sparsity that reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one conjugate pair of eigenvalues at each step using symmetric output feedback gain matrices, and the eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley Research Center (Controls-Structures-Interaction Evolutionary model) is presented to demonstrate the feasibility of the proposed approach.