Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Incompressible Navier-Stokes, GMRES, and Multi-Element Airfoils

The talk will discuss current work in developing and applying the INS2D and INS3D flow solvers to problems in high-lift aerodynamics, primarily multi-element airfoils. These flow solvers solve the incompressible Navier-Stokes equations using the method of artificial compressibility. High-lift system flowfields are perhaps the most challenging aeronautical configurations for CFD. Difficulties with convergence for multi-element configurations using fine-grid overset meshes has led to research in using Krylov-space iterative matrix solvers. In particular, a Generalized Minimum Residual (GMRES) solver has found to be a dramatic improvement over a previously used Gauss-Seidal line-relaxation solver.

Rogers, Stuart↗
NASA NTRS2000

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗
NASA NTRS1992

Environmental effects on long term behavior of composite laminates

Model equations are presented for approximate methods simulating the long-term behavior of composite materials and structures in hot/humid service environments. These equations allow laminate property upgradings with time, and can account for the effects of service environments on creep response. These methodologies are illustrated for various individual and coupled temperature/moisture, longitudinal/transverse, and composite material type cases. Creep deformation is noted to rise dramatically for cases of matrix-borne, but not of fiber-borne, loading in hot, humid environments; the coupled influence of temperature and moisture is greater than a mere combination of their individual influences.

Singhal, S. N.↗
NASA NTRS1996

Cyclic Load Effects on Long Term Behavior of Polymer Matrix Composites

A methodology to compute the fatigue life for different ratios, r, of applied stress to the laminate strength based on first ply failure criteria combined with thermal cyclic loads has been developed and demonstrated. Degradation effects resulting from long term environmental exposure and thermo-mechanical cyclic loads are considered in the simulation process. A unified time-stress dependent multi-factor interaction equation model developed at NASA Lewis Research Center has been used to account for the degradation of material properties caused by cyclic and aging loads. Effect of variation in the thermal cyclic load amplitude on a quasi-symmetric graphite/epoxy laminate has been studied with respect to the impending failure modes. The results show that, for the laminate under consideration, the fatigue life under combined mechanical and low thermal amplitude cyclic loads is higher than that due to mechanical loads only. However, as the thermal amplitude increases, the life also decreases. The failure mode changes from tensile under mechanical loads only to the compressive and shear at high mechanical and thermal loads. Also, implementation of the developed methodology in the design process has been discussed.

Shah, A. R.↗
NASA NTRS1975

JANNAF boundary layer integral matrix procedure

Accurate predictions of the thrust loss due to boundary layer effects and of the wall heat flux are very important to the design and performance evaluation of rocket nozzles. A method used in analytical procedures for liquid fuel rocket engine performance prediction and evaluation is presented. A computer program is examined that is a fast and accurate procedure for solving the set of boundary layer equation (momentum, energy, and species) for laminar or turbulent, chemically reacting flows with a wide variety of boundary conditions. Results of comparison of the various turbulent models are presented. A summary of the modifications and additions to the program is examined.

Evans, R. M.↗
NASA NTRS1973

Matrix Perturbation Techniques in Structural Dynamics

Matrix perturbation are developed techniques which can be used in the dynamical analysis of structures where the range of numerical values in the matrices extreme or where the nature of the damping matrix requires that complex valued eigenvalues and eigenvectors be used. The techniques can be advantageously used in a variety of fields such as earthquake engineering, ocean engineering, aerospace engineering and other fields concerned with the dynamical analysis of large complex structures or systems of second order differential equations. A number of simple examples are included to illustrate the techniques.

Caughey, T. K.↗
NASA NTRS1972

Development of programs for computing characteristics of ultraviolet radiation

Efficient programs were developed for computing all four characteristics of the radiation scattered by a plane-parallel, turbid, terrestrial atmospheric model. They were developed (FORTRAN 4) and tested on the IBM /360 computers with 2314 direct access storage facility. The storage requirement varies between 200K and 750K bytes depending upon the task. The scattering phase matrix (or function) is expanded in a Fourier series whose number of terms depend upon the zenith angles of the incident and scattered radiations, as well as on the nature of aerosols. A Gauss-Seidel procedure is used for obtaining the numerical solution of the transfer equation.

Dave, J. V.↗
NASA NTRS1984

Spatial variation analyses of Thematic Mapper data for the identification of linear features in agricultural landscapes

A need exists for digitized information pertaining to linear features such as roads, streams, water bodies and agricultural field boundaries as component parts of a data base. For many areas where this data may not yet exist or is in need of updating, these features may be extracted from remotely sensed digital data. This paper examines two approaches for identifying linear features, one utilizing raw data and the other classified data. Each approach uses a series of data enhancement procedures including derivation of standard deviation values, principal component analysis and filtering procedures using a high-pass window matrix. Just as certain bands better classify different land covers, so too do these bands exhibit high spectral contrast by which boundaries between land covers can be delineated. A few applications for this kind of data are briefly discussed, including its potential in a Universal Soil Loss Equation Model.

Pelletier, R. E.↗
NASA NTRS1989

State-space solutions to standard H2 and H(infinity) control problems

Simple state-space formulas are derived for all controllers solving the standard H(infinity) problem of finding, for a given number gamma greater than 0, all controllers such that the H(infinity) norm of the closed-loop transfer function is (strictly) less than gamma. It is known that a controller exists if and only if the unique stabilizing solutions to two algebraic Riccati equations are positive definite and the spectral radius of their product is less than gamma squared. Under these conditions, a parameterization of all controllers solving the problem is given as a linear fractional transformation (LFT) on a contractive stable free parameter. The state dimension of the coefficient matrix for the LFT, constructed using the two Riccati solutions, equals that of the plant and has a separation structure reminiscent of classical LQG (i.e., H2) theory. A standard H2 solution is developed in parallel.

Doyle, John C.↗
NASA NTRS1992

Reference equations of motion for automatic rendezvous and capture

The analysis presented in this paper defines the reference coordinate frames, equations of motion, and control parameters necessary to model the relative motion and attitude of spacecraft in close proximity with another space system during the Automatic Rendezvous and Capture phase of an on-orbit operation. The relative docking port target position vector and the attitude control matrix are defined based upon an arbitrary spacecraft design. These translation and rotation control parameters could be used to drive the error signal input to the vehicle flight control system. Measurements for these control parameters would become the bases for an autopilot or feedback control system (FCS) design for a specific spacecraft.

Henderson, David M.↗
NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗
NASA NTRS1993

A multigrid solver for the semiconductor equations

We present a multigrid solver for the exponential fitting method. The solver is applied to the current continuity equations of semiconductor device simulation in two dimensions. The exponential fitting method is based on a mixed finite element discretization using the lowest-order Raviart-Thomas triangular element. This discretization method yields a good approximation of front layers and guarantees current conservation. The corresponding stiffness matrix is an M-matrix. 'Standard' multigrid solvers, however, cannot be applied to the resulting system, as this is dominated by an unsymmetric part, which is due to the presence of strong convection in part of the domain. To overcome this difficulty, we explore the connection between Raviart-Thomas mixed methods and the nonconforming Crouzeix-Raviart finite element discretization. In this way we can construct nonstandard prolongation and restriction operators using easily computable weighted L(exp 2)-projections based on suitable quadrature rules and the upwind effects of the discretization. The resulting multigrid algorithm shows very good results, even for real-world problems and for locally refined grids.

Bachmann, Bernhard↗
NASA NTRS1995

A new technique for ground simulation of hypervelocity debris

A series of hypervelocity damage experiments were preformed on spacecraft materials. These experiments employed a technique which accelerates micro flyer plates simulating space debris traveling at 3 to 8 km/sec. The apparatus used to propel the micro flyer plates was compact and fit well into a space environmental chamber equipped with instrumentation capable of analyzing the vapor ejected from the sample. Mechanical damage to the sample was also characterized using optical and scanning electron microscpopy. Data for this work was obtained from hypervelocity impacts on a polysulfone resin and a graphite polysulfone composite. Polysulfone was selected because it was flown on the Long Duration Exposure Facility (LDEF) which spent several years in low earth orbit (LEO). Chemistry of the vapor produced by the impact was analyzed with a time of flight mass spectrometer, (TOFMS). This represents the first time that ejected vapors from hypervelocity collisions were trapped and analyzed with a mass spectrometer. With this approach we are able to study changes in the vapor chemistry as a function of time after impact, obtain a velocity measurement of the vapor, and estimate a temperature of the surface at time of impact using dynamic gas equations. Samples of the vapor plume may be captured and examined by transmission electron microscopy. Studies were also conducted to determine mechanical damage to a graphite polysulfone composite and a polysulfone resin. Impact craters were examined under optical and scanning electron microscopes. The collision craters in the matrix were typical of those shown in conventional shock experiments. However, the hypervelocity collisions with the graphite polysulfone composite were remarkably different than those with the resin.

Roybal, R.↗
NASA NTRS1995

Implicit method for the computation of unsteady flows on unstructured grids

An implicit method for the computation of unsteady flows on unstructured grids is presented. Following a finite difference approximation for the time derivative, the resulting nonlinear system of equations is solved at each time step by using an agglomeration multigrid procedure. The method allows for arbitrarily large time steps and is efficient in terms of computational effort and storage. Inviscid and viscous unsteady flows are computed to validate the procedure. The issue of the mass matrix which arises with vertex-centered finite volume schemes is addressed. The present formulation allows the mass matrix to be inverted indirectly. A mesh point movement and reconnection procedure is described that allows the grids to evolve with the motion of bodies. As an example of flow over bodies in relative motion, flow over a multi-element airfoil system undergoing deployment is computed.

Venkatakrishnan, V.↗
NASA NTRS2005

A Novel Porosity Model for Use in Hydrocode Simulations

Introduction: Numerical modeling of impact cratering has reached a high degree of sophistication; however, the treatment of porous materials still poses a large problem in hydrocode calculations. Porosity plays only a minor role in the formation of large craters on most planetary objects, but impacts on comets are believed to be highly affected by the presence of porosity, which may be as much as 80%. The upcoming Deep Impact Mission (launched January 2005) will provide more detailed data about the composition of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of its nucleus. The numerical simulations of such impact events requires an appropriate model for how pore space in the comet is crushed out during the violent initial stage of the impact event. Most hydro-codes compute the pressure explicitly using an "equation of state" (EOS) for each material, which relates changes in density and internal energy to changes in pressure. The added complication introduced by porosity is that changes in a material s density are due to both the closing of pore space (compaction) and compression of the matrix. The amount of resistance to volume change and the amount of irreversible work done during these two processes is very different; it is far easier to compact a porous material sample than to compress a non-porous sample of the same material. As an alternative to existing porosity models, like the Pdot(alpha) model [1], we present a novel approach for dealing with the compaction of porosity in hydrocode calculations.

Wuennemann, K.↗
NASA NTRS2007

Composite Nanomechanics: A Mechanistic Properties Prediction

A unique mechanistic theory is described to predict the properties of nanocomposites. The theory is based on composite micromechanics with progressive substructuring down to a nanoscale slice of a nanofiber where all the governing equations are formulated. These equations hav e been programmed in a computer code. That computer code is used to predict 25 properties of a mononanofiber laminate. The results are pr esented graphically and discussed with respect to their practical sig nificance. Most of the results show smooth distributions. Results for matrix-dependent properties show bimodal through-the-thickness distr ibution with discontinuous changes from mode to mode.

Chamis, Christos C.↗
NASA NTRS1988

A cumulative damage model for continuous fiber composite laminates with matrix cracking and interply delaminations

The present cumulative damage model for the prediction of stiffness loss in graphite/epoxy laminates applies a thermomechanical constitutive theory for elastic composites with distributed damage. The model proceeds from a continuum mechanics and thermodynamics approach in which the distributed damage is characterized by a set of second-order tensor-valued internal state variables. A set of damage-dependent laminated plate equations is obtained; this is developed by modifying classical Kirchhoff plate theory.

Allen, David H.↗
NASA NTRS1985

Parametric Stiffness Control of Flexible Structures

An unconventional method for control of flexible space structures using feedback control of certain elements of the stiffness matrix is discussed. The advantage of using this method of configuration control is that it can be accomplished in practical structures by changing the initial stress state in the structure. The initial stress state can be controlled hydraulically or by cables. The method leads, however, to nonlinear control equations. In particular, a long slender truss structure under cable induced initial compression is examined. both analytical and numerical analyses are presented. Nonlinear analysis using center manifold theory and normal form theory is used to determine criteria on the nonlinear control gains for stable or unstable operation. The analysis is made possible by the use of the exact computer algebra system MACSYMA.

Moon, F. C.↗