A least squares iterative method for singular equations
Least squares iterative method for solving simultaneous linear equations having singular coefficient matrix
Tewarson, R. P.↗
SEARCH · Search NASA
Results for “Matrix equations”
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.
FCAP - A new tool for the performance and structural analysis for complex flexible aircraft with active control
A very general formulation has been developed for the analysis of complex flexible-aircraft configurations in presence of active control. The formulation is incorporated into a computer program, Flight Control Analysis Program (FCAP), which is designed in a modular fashion to incorporate aircraft dynamics and aerodynamics for complex configurations, as well as sensor, actuator and control logic element dynamics. Formulation of the total aircraft dynamic system is accomplished in matrix form by casting equations in state vector format. The system stability and performance are determined in either the frequency or time domain using classical analysis techniques. Fully unsteady aerodynamics is used for the evaluation of the flutter characteristics as well as the gust response of the aircraft.
Morino, L.↗
Modal cost analysis for linear matrix-second-order systems
Reduced models and reduced controllers for systems governed by matrix-second-order differential equations are obtained by retaining those modes which make the largest contributions to quadratic control objectives. Such contributions, expressed in terms of modal data, used as mode truncation criteria, allow the statement of the specific control objectives to influence the early model reduction from very high order models which are available, for example, from finite element methods. The relative importance of damping, frequency, and eigenvector in the mode truncation decisions are made explicit for each of these control objectives: attitude control, vibration suppression and figure control. The paper also shows that using modal cost analysis (MCA) on the closed loop modes of the optimally controlled system allows the construction of reduced control policies which feedback only those closed loop modal coordinates which are most critical to the quadratic control performance criterion. In this way, the modes which should be controlled (and hence the modes which must be observable by choice of measurements), are deduced from truncations of the optimal controller.
Skelton, R. E.↗
Linearized dynamical model for the NASA/IEEE SCOLE configuration
The linearized equation of motion for the NASA/IEEE SCOLE configuration are developed. The derivation is based on the method of Lagrange and the equations are assembled into matrix second order form.
Hotz, A. F.↗
A generalized analysis of dual-polarization radar measurements of rain
In the study of meteorological phenomena by radar, it has been shown that dual-polarization radar measurements can yield useful data about precipitation structures within the radar beam. The techniques available for dual-polarized radar measurement are differential reflectivity, or Z(DR), linear depolarization ratio (LDR), and circular depolarization ratio (CDR). This paper presents a theoretical analysis of backscattering from a volume filled with raindrops using Rayleigh scattering theory. A generalized radar equation is given in matrix form, and the dependence of backscatter measurements in terms of a rainstorm model including raindrop shape, size, and canting angle is examined. Calculated results which show the dependence of Z(DR) on beam elevation angle and polarization angle are presented. Some results for LDR and CDR are presented, showing how these quantities depend on raindrop parameters.
Stapor, D. P.↗
Response and characteristics of structures subjected to S-H waves
A study of the dynamic characteristics of a coupled translational-rotational system is given. The formulation of the problem considers the soil-structure interaction effects by utilizing the impedance functions at the foundation of a structure. Due to the fact that the coefficient matrix in the characteristic equation is frequency dependent in nature, iterations have to be performed to find the nature frequencies of the system. Examples and discussions are presented. Comparisons of the analytical results from various approaches are also given.
Wu, S. T.↗
A complete spherical harmonic approach to luni-solar tides
In this work a spherical harmonic theory of ocean tides is presented. The theory is based on Laplace tide equations modified to include turbulence with constant eddy viscosity, linearized bottom friction, and oceanic loading and self-gravitation. Variable bathymetry is also treated in harmonic terms, and no-flow boundary conditions are applied at continental coastlines. The tide and boundary constraint equations are reduced to matrix form and solved by a weighted least-squares procedure. Five zonal luni-solar tides, ranging in period from 14 days to 18.6 yr, are investigated using the theory; such tides have typically been difficult to compute using traditional numerical approaches. The polar motion and changes in the length of day induced by these long-period tides are calculated. Tidal solutions are compared extensively with results from other tidal theories and from recent satellite and sea-level observations. The greatest limitation to accurate prediction of zonal tides - for any theory - appears to be the marginal failure of all tide theories to conserve mass globally; the use of additional mass constraints may be warranted.
Dickman, S. R.↗
Inverse dynamics of adaptive structures used as space cranes
As a precursor to the real-time control of fast moving adaptive structures used as space cranes, a formulation is given for the flexibility induced motion relative to the nominal motion (i.e., the motion that assumes no flexibility) and for obtaining the open loop time varying driving forces. An algorithm is proposed for the computation of the relative motion and driving forces. The governing equations are given in matrix form with explicit functional dependencies. A simulator is developed to implement the algorithm on a digital computer. In the formulations, the distributed mass of the crane is lumped by two schemes, vz., 'trapezoidal' lumping and 'Simpson's rule' lumping. The effects of the mass lumping schemes are shown by simulator runs.
Das, S. K.↗
Central difference TVD and TVB schemes for time dependent and steady state problems
We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.
Jorgenson, P.↗
Development of a steady potential solver for use with linearized, unsteady aerodynamic analyses
A full potential steady flow solver (SFLOW) developed explicitly for use with an inviscid unsteady aerodynamic analysis (LINFLO) is described. The steady solver uses the nonconservative form of the nonlinear potential flow equations together with an implicit, least squares, finite difference approximation to solve for the steady flow field. The difference equations were developed on a composite mesh which consists of a C grid embedded in a rectilinear (H grid) cascade mesh. The composite mesh is capable of resolving blade to blade and far field phenomena on the H grid, while accurately resolving local phenomena on the C grid. The resulting system of algebraic equations is arranged in matrix form using a sparse matrix package and solved by Newton's method. Steady and unsteady results are presented for two cascade configurations: a high speed compressor and a turbine with high exit Mach number.
Hoyniak, Daniel↗
A procedure for utilization of a damage-dependent constitutive model for laminated composites
Described here is the procedure for utilizing a damage constitutive model to predict progressive damage growth in laminated composites. In this model, the effects of the internal damage are represented by strain-like second order tensorial damage variables and enter the analysis through damage dependent ply level and laminate level constitutive equations. The growth of matrix cracks due to fatigue loading is predicted by an experimentally based damage evolutionary relationship. This model is incorporated into a computer code called FLAMSTR. This code is capable of predicting the constitutive response and matrix crack damage accumulation in fatigue loaded laminated composites. The structure and usage of FLAMSTR are presented along with sample input and output files to assist the code user. As an example problem, an analysis of crossply laminates subjected to two stage fatigue loading was conducted and the resulting damage accumulation and stress redistribution were examined to determine the effect of variations in fatigue load amplitude applied during the first stage of the load history. It was found that the model predicts a significant loading history effect on damage evolution.
Lo, David C.↗
Central difference TVD and TVB schemes for time dependent and steady state problems
We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.
Jorgenson, P.↗
Electron impact polarization of atomic spectral lines. I - A general theoretical scheme
A suitable theoretical scheme able to describe, in a wide variety of astrophysical situations, the phenomenon of atomic line polarization by electron impact is developed. Starting from the general principles of quantum mechanics and assuming the Born approximation, the rate equations for the density matrix elements of a multilevel atomic system, interacting with a nonrelativistic electron beam having any kind of angular distribution, are derived in full generality. The resulting theory generalizes the previous ones by accounting for the collisional rates and the cross sections concerning both inelastic and superelastic collisions (in any geometrical situation), and, moreover, by taking into account the coherences among Zeeman sublevels split by a magnetic field. As an example of particular relevance, the general formulas derived in the first sections of the paper are subsequently particularized to the case of the electric dipole interaction.
Fineschi, Silvano↗
Robust root clustering for linear uncertain systems using generalized Lyapunov theory
Consideration is given to the problem of matrix root clustering in subregions of a complex plane for linear state space models with real parameter uncertainty. The nominal matrix root clustering theory of Gutman & Jury (1981) using the generalized Liapunov equation is extended to the perturbed matrix case, and bounds are derived on the perturbation to maintain root clustering inside a given region. The theory makes it possible to obtain an explicit relationship between the parameters of the root clustering region and the uncertainty range of the parameter space.
Yedavalli, R. K.↗
Computational strategies in the dynamic simulation of constrained flexible MBS
This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.
Amirouche, F. M. L.↗
Mobility of Electron in DNA Crystals by Laser Radiation
The mobility of electrons in laser radiated DNA is closed to the energy transfer and energy migration of a biological molecule. Arrhenius has studied the conductivity of the electrons in a biological molecule. But his result is far from the experimental result and meanwhile the relation between some parameters in his theory and the micro-quantities in DNA is not very clear. In this paper, we propose a new phonon model of electron mobility in DNA and use Lippman-Schwinger equation and S-matrix theory to study the mobility of electrons in DNA crystal. The result is relatively close to the experiment result and some parameters in Arrhenius theory are explained in our work.
Zhang, Kaixi↗
Generalization of the Euler Angles
It is shown that the Euler angles can be generalized to axes other than members of an orthonormal triad. As first shown by Davenport, the three generalized Euler axes, hereafter: Davenport axes, must still satisfy the constraint that the first two and the last two axes be mutually perpendicular if these axes are to define a universal set of attitude parameters. Expressions are given which relate the generalized Euler angles, hereafter: Davenport angles, to the 3-1-3 Euler angles of an associated direction-cosine matrix. The computation of the Davenport angles from the attitude matrix and their kinematic equation are presented. The present work offers a more direct development of the Davenport angles than Davenport's original publication and offers additional results.
Bauer, Frank H.↗
Many-Body Effects in Quantum-Well Intersubband Transitions
Intersubband polarization couples to collective excitations of the interacting electron gas confined in a semiconductor quantum well (Qw) structure. Such excitations include correlated pair excitations (repellons) and intersubband plasmons (ISPs). The oscillator strength of intersubband transitions (ISBTs) strongly varies with QW parameters and electron density because of this coupling. We have developed a set of kinetic equations, termed the intersubband semiconductor Bloch equations (ISBEs), from density matrix theory with the Hartree-Fock approximation, that enables a consistent description of these many-body effects. Using the ISBEs for a two-conduction-subband model, various many-body effects in intersubband transitions are studied in this work. We find interesting spectral changes of intersubband absorption coefficient due to interplay of the Fermi-edge singularity, subband renormalization, intersubband plasmon oscillation, and nonparabolicity of bandstructure. Our results uncover a new perspective for ISBTs and indicate the necessity of proper many-body theoretical treatment in order for modeling and prediction of ISBT line shape.
Li, Jian-Zhong↗