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 523 records · Page 29

Using the NASTRAN Thermal Analyzer to simulate a flight scientific instrument package

The NASTRAN Thermal Analyzer has proven to be a unique and useful tool for thermal analyses involving large and complex structures where small, thermally induced deformations are critical. Among its major advantages are direct grid point-to-grid point compatibility with large structural models; plots of the model that may be generated for both conduction and boundary elements; versatility of applying transient thermal loads especially to repeat orbital cycles; on-line printer plotting of temperatures and rate of temperature changes as a function of time; and direct matrix input to solve linear differential equations on-line. These features provide a flexibility far beyond that available in most finite-difference thermal analysis computer programs.

Lee, H.-P.↗

Digital explicit model following with unstable model dynamics

In the optimal regulator formulation of discrete explicit model following the Riccati equation may fail to reach a steady-state value for model dynamics which are not asymptotically stable. Such conditions often arise in aircraft applications when flying quality criteria based on step inputs are used to define the model equations. Mathematical conditions are presented which insure a steady-state value of the model-following gain matrix regardless of the behavior of the underlying Riccati equation. These results are applied to the design of a model-following controller for the lateral motion of a typical fighter aircraft using unstable model equations.

Armstrong, E. S.↗

Rocket engine thrust chamber heat transfer calculation and analysis

A parametric study of the heat transfer rate along the wall of a rocket nozzle is presented. The influences of different parameters; laminar and turbulent Lewis number, mixture ratio, initial wall temperature distribution, and eddy viscosity, were considered. The numerical evaluation of these influences on heat transfer rate was done by using three different compressible, reacting laminar and turbulent boundary layer computer programs; MABL (Mass Addition Boundary Layer Program), MABL-KE (MABL program is modified to include turbulent kinetic energy equation), and BLIMP (Boundary Layer Integral Matrix Procedure). This study also provided an excellent opportunity to evaluate the efficiencies of these three computer programs and to suggest one of them for future computational purposes.

Saha, H.↗

Discretized partial differential equations - Examples of control systems defined on modules

The purpose of this paper is to show how the important problems of linear system theory can be solved concisely for a particular class of linear systems, namely block circulant systems, by exploiting the algebraic structure. This type of system arises in lumped approximations to linear partial differential equations. The computation of the transition matrix, the variation of constants formula, observability, controllability, pole allocation, realization theory, stability and quadratic optimal control are discussed. In principle, all questions which are solved here could also be solved by standard methods; the present paper clearly exposes the structure of the solution, and thus permits various savings in computational effort.

Brockett, R. W.↗

Rotor dynamic state and parameter identification from simulated forward flight transients, part 1

State and parameter identifications from simulated forward flight blade flapping measurements are presented. The transients were excited by progressing cyclic pitch stirring or by hub stirring with constant stirring acceleration. Rotor dynamic inflow models of varying degree of sophistication were used from a one parameter inflow model (equivalent Lock number) to an eight parameter inflow model. The maximum likelihood method with assumed fixed measurement error covariance matrix was applied. The rotor system equations for both fixed hub and tilting hub are given. The identified models were verified by comparing true responses with predicted responses. An optimum utilization of the simulated measurement data can be defined. From the numerical results it can be anticipated that brief periods of either accelerated cyclic pitch stirring or of hub stirring are sufficient to extract with adequate accuracy up to 8 rotor dynamic inflow parameters plus the blade Lock number from the transients.

Hohenemser, K. H.↗

Method of fan sound mode structure determination computer program user's manual: Microphone location program

A computer user's manual describing the operation and the essential features of the microphone location program is presented. The Microphone Location Program determines microphone locations that ensure accurate and stable results from the equation system used to calculate modal structures. As part of the computational procedure for the Microphone Location Program, a first-order measure of the stability of the equation system was indicated by a matrix 'conditioning' number.

Pickett, G. F.↗

Dynamic modeling of structures from measured complex modes

A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.

Ibrahim, s. R.↗

Accelerated testing of composites

It is shown that the Zhurkov method for testing the strength of solids can be applied to dynamic tension and to cyclic loading and provides a viable approach to accelerated testing of composites. Data from the literature are used to demonstrate a straightforward application of the method to dynamic tension of glass fiber and cyclic loading for glass/polymer, metal matrix, and graphite/epoxy composites. Zhurkov's equation can be used at relatively high loads to obtain failure times at any temperature of interest. By taking a few data points at one or two other temperatures the spectrum of failure times can be expanded to temperatures not easily accessible.

Papazian, H. A.↗

Computer programs for smoothing and scaling airfoil coordinates

Detailed descriptions are given of the theoretical methods and associated computer codes of a program to smooth and a program to scale arbitrary airfoil coordinates. The smoothing program utilizes both least-squares polynomial and least-squares cubic spline techniques to smooth interatively the second derivatives of the y-axis airfoil coordinates with respect to a transformed x-axis system which unwraps the airfoil and stretches the nose and trailing-edge regions. The corresponding smooth airfoil coordinates are then determined by solving a tridiagonal matrix of simultaneous cubic-spline equations relating the y-axis coordinates and their corresponding second derivatives. A technique for computing the camber and thickness distribution of the smoothed airfoil is also discussed. The scaling program can then be used to scale the thickness distribution generated by the smoothing program to a specific maximum thickness which is then combined with the camber distribution to obtain the final scaled airfoil contour. Computer listings of the smoothing and scaling programs are included.

Morgan, H. L., Jr.↗

Dynamic Modeling of Structures from Measured Complex Modes

A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.

Ibrahim, S. R.↗

Interactive application of quadratic expansion of chi-square statistic to nonlinear curve fitting

This report contains a detailed theoretical description of an all-purpose, interactive curve-fitting routine that is based on P. R. Bevington's description of the quadratic expansion of the Chi-Square statistic. The method is implemented in the associated interactive, graphics-based computer program. Taylor's expansion of Chi-Square is first introduced, and justifications for retaining only the first term are presented. From the expansion, a set of n simultaneous linear equations is derived, then solved by matrix algebra. A brief description of the code is presented along with a limited number of changes that are required to customize the program of a particular task. To evaluate the performance of the method and the goodness of nonlinear curve fitting, two typical engineering problems are examined and the graphical and tabular output of each is discussed. A complete listing of the entire package is included as an appendix.

Badavi, F. F.↗

Vector solution for the mean electromagnetic fields in a layer of random particles

The mean electromagnetic fields are found in a layer of randomly oriented particles lying over a half space. A matrix-dyadic formulation of Maxwell's equations is employed in conjunction with the Foldy-Lax approximation to obtain equations for the mean fields. A two variable perturbation procedure, valid in the limit of small fractional volume, is then used to derive uncoupled equations for the slowly varying amplitudes of the mean wave. These equations are solved to obtain explicit expressions for the mean electromagnetic fields in the slab region in the general case of arbitrarily oriented particles and arbitrary polarization of the incident radiation. Numerical examples are given for the application to remote sensing of vegetation.

Lang, R. H.↗

Control of a slow moving space crane as an adaptive structure

Assuming that the space crane is an adaptive structure with length-adjustable bars and taking as controls the length-adjustments of these bars, the computation of the incremental controls corresponding to the motion of a payload along its minimum-energy trajectory is given in terms of the inverse-transpose of matrix B of the joint equilibrium equations Bs = p, where s lists the bar forces and p lists the nodal loads. The compensation of the controls for elastic deformations and support movements are shown. It is also shown that the computations may be done automatically and in real time by an attached processor once the characteristics of the crane's maneuver are keyed in.

Utku, S.↗

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗

A flux-split solution procedure for unsteady inlet flows

The unstart and restart of an axisymmetric inlet is investigated using a flux-split procedure applied to the Euler and Reduced Navier Stokes (RNS) equations. A time consistent direct sparse matrix solver is applied to compute the transient flow field both internal and external to the inlet. Time varying oblique and normal shocks are captured. The code is quite general and is applicable for subsonic, transonic and supersonic free streams. The current analysis is concerned with supersonic flight conditions.

Pordal, H. S.↗

A diagonal implicit scheme for computing flows with finite-rate chemistry

A new algorithm for solving steady, finite-rate chemistry, flow problems is presented. The new scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The source Jacobian matrix is replaced by a diagonal matrix which is tailored to account for the fastest reactions in the chemical system. A point-implicit procedure is discussed and then the algorithm is included into the LU-SGS scheme. Solutions are presented for hypervelocity reentry and Hydrogen-Oxygen combustion. For the LU-SGS scheme a CFL number in excess of 10,000 has been achieved.

Eberhardt, Scott↗

Transient behavior of supersonic flow through inlets

A solution technique to compute inlet flow behavior is presented. The phenomena of inlet unstart and restart are investigated using a flux-split procedure applied to the Euler and Reduced Navier Stokes (RNS) equations. A time consistent direct sparse matrix solver in conjunction with a domain decomposition strategy is applied to compute the transient flow behavior both internal and external to the inlet. Time varying shocks and time varying recirculation regions are efficiently analyzed. The code is quite general and is suitable for the computation of flow for a wide variety of geometries and over a wide range of Mach and Reynolds numbers.

Pordal, H. S.↗

A manual for PARTI runtime primitives

Primitives are presented that are designed to help users efficiently program irregular problems (e.g., unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial differential equations solvers) on distributed memory machines. These primitives are also designed for use in compilers for distributed memory multiprocessors. Communications patterns are captured at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗