Search NASA⌕ Search

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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗
NASA NTRS1976

Derivatives of a statically reduced stiffness matrix with respect to sizing variables

An expression is obtained for the first derivatives with respect to the sizing variables of a statically reduced stiffness matrix that is a nonlinear function of the sizing variables, where the unreduced stiffness matrix is a linear function of the sizing variables. An accepted procedure to reduce the number of degrees of freedom is to eliminate a number of nodal displacements from the degrees of freedom such that the accuracy of the flutter analysis is not significantly affected. In a typical optimization procedure with flutter constraints, the derivative of the stiffness matrix may be used in a form that contains the characteristic vector of the flutter matrix equation and the transpose of the characteristic vector of the adjoint flutter matrix equation corresponding to a particular solution of the flutter equation.

Oconnell, R. F.↗
NASA NTRS1983

An analytical investigation of shape control of large space structures by applied temperatures

An analytical procedure for the static shape control of flexible space structures subjected to thermal distortions is developed which is based on prescribing temperatures in control elements having much higher coefficients of thermal expansion than the main structure. The temperatures at the control elements are defined so as to minimize the overall thermal distortion of the structure from its ideal shape, and a matrix equation is obtained which can be solved for the set of optimum control temperatures. A formulation of the procedure for continuous structures governed by differential equations and a formulation for discrete (finite element modeled) structures governed by matrix equations are presented. The equations from the continuous formulation are employed for the shape control of a simple beam distorted by nonuniform heating, and the discrete formulation is applied in a general purpose finite-element structural analysis computer program for the shape control of a 750 m radiometer antenna reflector dish subjected to orbital heating. A reduction in thermal distortion by a factor of nearly 50 was obtained with the use of only seven control elements. Results for four different sets of control locations for the antenna are presented in which reductions in distortion of up to a factor of four were obtained.

Haftka, R. T.↗
NASA NTRS1992

Scattering from dielectric structures above impedance surfaces and resistive sheets

Interest in understanding of electromagnetic interaction with rough surfaces has prompted the study of scattering from typical dielectric humps over impedance surfaces. It is shown that the Green's function of the problem for a resistive sheet resembles that of the impedance surface. Hence both problems are considered here. A numerical solution for the scattered field of a two-dimensional dielectric object, possibly inhomogeneous, with arbitrary cross section above the impedance surface or resistive sheet is sought. First the Green's function of the problem is derived based on the exact image theory. This form of the Green's function is amenable to numerical computation. Then the induced polarization currents are calculated by casting the integral equations into a matrix equation via the method of moments. Numerical problems in calculation of the Green's function when both source and observation points are close to the surface are discussed. Comparison of numerical results with a perturbation solution shows excellent agreement between the two methods.

Sarabandi, Kamal↗
NASA NTRS1994

Scattering From A Thick Dichroic Microwave Reflector

Thick Frequency Selective Surface with Rectangular Apertures, TFSSRA, computer program developed to calculate scattering parameters of thick frequency-selective surface with rectangular apertures on skew grid at oblique angle of incidence. Integral equation transformed into matrix equation, which is then solved. Written in FORTRAN 77 with single precision.

Chen, Jacqueline C.↗
NASA NTRS1985

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗
NASA NTRS1977

Application of symbolic/numeric matrix solution techniques to the NASTRAN program

The matrix solving algorithm of any finite element algorithm is extremely important since solution of the matrix equations requires a large amount of elapse time due to null calculations and excessive input/output operations. An alternate method of solving the matrix equations is presented. A symbolic processing step followed by numeric solution yields the solution very rapidly and is especially useful for nonlinear problems.

Buturla, E. M.↗
NASA NTRS1996

Asymptotic Waveform Evaluation (AWE) Technique for Frequency Domain Electromagnetic Analysis

The Asymptotic Waveform Evaluation (AWE) technique is applied to a generalized frequency domain electromagnetic problem. Most of the frequency domain techniques in computational electromagnetics result in a matrix equation, which is solved at a single frequency. In the AWE technique, the Taylor series expansion around that frequency is applied to the matrix equation. The coefficients of the Taylor's series are obtained in terms of the frequency derivatives of the matrices evaluated at the expansion frequency. The coefficients hence obtained will be used to predict the frequency response of the system over a frequency range. The detailed derivation of the coefficients (called 'moments') is given along with an illustration for electric field integral equation (or Method of Moments) technique. The radar cross section (RCS) frequency response of a square plate is presented using the AWE technique and is compared with the exact solution at various frequencies.

Cockrell, C. R.↗
NASA NTRS1980

Collisional redistribution of radiation. II - The effects of degeneracy on the equations of motion for the density matrix. III - The equation of motion for the correlation function and the scattered spectrum

The effect of correlations between an absorber atom and perturbers in the binary-collision approximation are applied to degenerate atomic systems. A generalized absorption profile which specifies the final state of the atom after an absorption event is related to the total intensities of Rayleigh scattering and fluorescence from the atom. It is suggested that additional dynamical information to that obtainable from ordinary absorption experiments is required in order to describe redistributed atomic radiation. The scattering of monochromatic radiation by a degenerate atom is computed in a binary-collision approximation; an equation of motion is derived for the correlation function which is valid outside the quantum-regression regime. Solutions are given for the weak-field conditions in terms of generalized absorption and emission profiles that depend on the indices of the atomic multipoles.

Burnett, K.↗
NASA NTRS2003

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗
NASA NTRS2004

A Model for the Oxidation of Carbon Silicon Carbide Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of carbon silicon carbide (C/SiC) composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The numerical method is demonstrated by utilizing the method to model the carbon oxidation and weight loss behavior of C/SiC specimens during thermogravimetric experiments. The numerical method is used to study the physics of carbon oxidation in carbon silicon carbide composites.

Sullivan, Roy M.↗
NASA NTRS2008

Oxidation Behavior of Carbon Fiber-Reinforced Composites

OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.

Sullivan, Roy M.↗
NASA NTRS1987

Partitioned Matrices for Combined Linear Systems

Multiple-input, multiple-output mathematical models combined. Report presents mathematical background for combination of linear multiple-input, multiple-output equations of subsystems into overall matrix equations for whole system. Development of such system mathematical models essential for design of complicated control systems.

Duke, Eugene L.↗
NASA NTRS1992

Active structural control for the mini-MAST and ACES structures

Two experiments based on the maximum entropy/optimal projection design approach are discussed which were conducted by Harris Corporation as part of the NASA Controls Structures Interaction Guest Investigator program. The design equations consist of four coupled matrix equations which specialize to the standard linear-quadratic Guassian Riccati equations, when the plant is known perfectly and a full order controller is desired. These experiments demonstrate successful control system design and implementation for flexible structures.

Collins, E. G., Jr.↗
NASA NTRS1979

Energy conserving norms for the solution of hyperbolic systems of partial differential equations

The hyperbolic system of partial differential equations with a real constant square coefficient matrix A is considered. The problem of finding an energy conserving norm for the solution of the system is reduced to the problem of characterizing those matrices appearing in the boundary conditions which satisfy two specific matrix equations. Necessary and sufficient conditions on the coefficient matrix A and the matrices appearing in boundary conditions are derived for an energy conserving norm. The conditions serve as criteria on a given system which determine whether or not the solution will have its energy conserved in some norm. Examples of specific systems and boundary conditions are also provided.

Gunzburger, M. D.↗
NASA NTRS1969

Finite element formulation for linear thermoviscoelastic materials

Report presents the finite difference equations in time and finite element matrix equations in space for general linear thermovisoelastic problems. The equations are derived for a general three-dimensional body but are applicable to one- and two-dimensional configurations with minor changes.

Chen, J. C.↗