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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A controls engineering approach for analyzing airplane input-output characteristics

An engineering approach for analyzing airplane control and output characteristics is presented. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled coordinates. The scaling is accomplished by applying various transformations to the system to employ prior engineering knowledge of the airplane physics. Two different analysis techniques are then explained. Modal analysis techniques calculate the influence of each system input on each fundamental mode of motion and the distribution of each mode among the system outputs. The optimal steady state response technique computes the blending of steady state control inputs that optimize the steady state response of selected system outputs. Analysis of an example airplane model is presented to demonstrate the described engineering approach.

Arbuckle, P. Douglas↗
NASA NTRS1991

Finite element solution of transient fluid-structure interaction problems

A finite element approach using NASTRAN is developed for solving time-dependent fluid-structure interaction problems, with emphasis on the transient scattering of acoustic waves from submerged elastic structures. Finite elements are used for modeling both structure and fluid domains to facilitate the graphical display of the wave motion through both media. For the liquid, the use of velocity potential as the fundamental unknown results in a symmetric matrix equation. The approach is illustrated for the problem of transient scattering from a submerged elastic spherical shell subjected to an incident tone burst. The use of an analogy between the equations of elasticity and the wave equation of acoustics, a necessary ingredient to the procedure, is summarized.

Everstine, Gordon C.↗
NASA NTRS1991

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

CRAY 2 and CRAY Y-MP evaluation runs are undertaken for a three-dimensional, elastoplastic FEM code implementation of the 'autotasking' parallel-processing technique, which is used in all code components except the matrix-equation solver. For a typical example problem in the case of a four-processor CRAY 2, a speedup factor of 2.1 was found achievable in a dedicated environment, and 1.7 in a multiuser environment; for a three-processor CRAY Y-MP, speedup in a multiuser environment was about 2.4. Only minimal effort was required to implement autotasking in this implementation code.

Shivakumar, K. N.↗
NASA NTRS1990

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗
NASA NTRS1991

Homotopy approach to optimal, linear quadratic, fixed architecture compensation

Optimal linear quadratic Gaussian compensators with constrained architecture are a sensible way to generate good multivariable feedback systems meeting strict implementation requirements. The optimality conditions obtained from the constrained linear quadratic Gaussian are a set of highly coupled matrix equations that cannot be solved algebraically except when the compensator is centralized and full order. An alternative to the use of general parameter optimization methods for solving the problem is to use homotopy. The benefit of the method is that it uses the solution to a simplified problem as a starting point and the final solution is then obtained by solving a simple differential equation. This paper investigates the convergence properties and the limitation of such an approach and sheds some light on the nature and the number of solutions of the constrained linear quadratic Gaussian problem. It also demonstrates the usefulness of homotopy on an example of an optimal decentralized compensator.

Mercadal, Mathieu↗
NASA NTRS1992

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

A CRAY parallel processing technique called autotasking was implemented in a three-dimensional elasto-plastic finite-element code. The technique was evaluated on two CRAY supercomputers, a CRAY 2 and a CRAY Y-MP. Autotasking was implemented in all major portions of the code, except the matrix equations solver. Compiler directives alone were not able to properly multitask the code; user-inserted directives were required to achieve better performance. It was noted that the connect time, rather than wall-clock time, was more appropriate to determine speedup in multiuser environments. For a typical example problem, a speedup of 2.1 (1.8 when the solution time was included) was achieved in a dedicated environment and 1.7 (1.6 with solution time) in a multiuser environment on a four-processor CRAY 2 supercomputer. The speedup on a three-processor CRAY Y-MP was about 2.4 (2.0 with solution time) in a multiuser environment.

Shivakumar, K. N.↗
NASA NTRS1993

Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Owocki, Stanley P.↗
NASA NTRS1993

Modifying real convolutional codes for protecting digital filtering systems

A novel method is proposed for protecting digital filters from temporary and permanent failures that are not easily detected by conventional fault-tolerant computer design principles, on the basis of the error-detecting properties of real convolutional codes. Erroneous behavior is detected by externally comparing the calculated and regenerated parity samples. Great simplifications are obtainable by modifying the code structure to yield simplified parity channels with finite impulse response structures. A matrix equation involving the original parity values of the code and the polynomial of the digital filter's transfer function is formed, and row manipulations separate this equation into a set of homogeneous equations constraining the modifying scaling coefficients and another set which defines the code parity values' implementation.

Redinbo, G. R.↗
NASA NTRS1994

Evaluation of several non-reflecting computational boundary conditions for duct acoustics

Several non-reflecting computational boundary conditions that meet certain criteria and have potential applications to duct acoustics are evaluated for their effectiveness. The same interior solution scheme, grid, and order of approximation are used to evaluate each condition. Sparse matrix solution techniques are applied to solve the matrix equation resulting from the discretization. Modal series solutions for the sound attenuation in an infinite duct are used to evaluate the accuracy of each non-reflecting boundary conditions. The evaluations are performed for sound propagation in a softwall duct, for several sources, sound frequencies, and duct lengths. It is shown that a recently developed nonlocal boundary condition leads to sound attenuation predictions considerably more accurate for short ducts. This leads to a substantial reduction in the number of grid points when compared to other non-reflecting conditions.

Watson, Willie R.↗
NASA NTRS1993

On the synthesis of resonance lines in dynamical models of structured hot-star winds

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Puls, J.↗
NASA NTRS1997

Implicit Kalman filtering

For an implicitly defined discrete system, a new algorithm for Kalman filtering is developed and an efficient numerical implementation scheme is proposed. Unlike the traditional explicit approach, the implicit filter can be readily applied to ill-conditioned systems and allows for generalization to descriptor systems. The implementation of the implicit filter depends on the solution of the congruence matrix equation (A1)(Px)(AT1) = Py. We develop a general iterative method for the solution of this equation, and prove necessary and sufficient conditions for convergence. It is shown that when the system matrices of an implicit system are sparse, the implicit Kalman filter requires significantly less computer time and storage to implement as compared to the traditional explicit Kalman filter. Simulation results are presented to illustrate and substantiate the theoretical developments.

Non-NASA Center↗
NASA NTRS2006

General Solution for Theoretical Packet Data Loss Rate

Communications systems which transfer blocks ("frames") of data must use a marker ("frame synchronization pattern") for identifying where a block begins. A technique ("frame synchronization strategy") is used to locate the start of each frame and maintain synchronization as additional blocks are processed. A device which strips out the frame synchronization pattern [FSP] and provides an "end of frame" pulse is called a frame synchronizer. As clock and data errors are introduced into the system, the start-of-block marker becomes displaced and/or corrupted. The capability of the frame synchronizer to stay locked to the pattern under these conditions is a figure of merit for the frame synchronization strategy. It is important to select a strategy which will stay locked nearly all the time at bit error rates where the data is usable. ("Bit error rate" [BER] is the fraction of binary bits which are inverted by passage through a communication system.) The fraction of frames that are discarded because the frame synchronizer is not locked is called "Percent Data Loss" or "Packet Data Loss rate" [PDL]. A general approach for accurately predicting PDL given BER was developed in Theoretical Percent Data Loss Calculation and Measurement Accuracy, T. P. Kelly, LESC-30554, December 1992. Kelly gave a solution in terms of matrix equations, and only addressed "level" channel encoding. This paper goes on to give a closed-form polynomial solution for the most common class of frame synchronizer strategies, and will also address "mark" and "space" (differential) channel encoding, and burst error environments. The paper is divided into four sections and follows a logically ordered presentation, with results developed before they are evaluated. However, most readers will derive the greatest benefit from this paper by treating the results as reference material. The result developed for differential encoding can be extended to other applications (like block codes) where the probability is needed that a block contains only a certain number of errors.

Lansdowne, Chatwin↗
NASA NTRSDate not supplied

Cathode Modeling of Solid-State Batteries

The search for safe, reliable, and compact high-capacity energy storage devices has led to increased interest in all-solid-state battery research. The use of solid electrolytes provides enhanced safety and durability due to their reduced flammability and increased mechanical strength compared to organic liquid electrolytes. Still, the use of solid electrolytes remains challenging. Computational modeling plays a substantial role in addressing these challenges. A particle dynamics electromechanical model for simulating electrochemical processes in a solid-state battery cathode will be presented. The model presents cathode microstructure at the particle level as a mixture of ionically conductive solid electrolyte particles, electrically conductive carbon additives, and cathodic reactant particles. After densification, the particle connectivity is analyzed to reconstruct the complex electric network connecting reactant particles with an anodic and cathodic current collectors through the electrolyte and carbon particles. The Kirchhoff’s matrix equation describing this electric network, is solved to obtain values of various critical parameters, such as the overall conductivity of the cathode for lithium ions and electrons, cathodic reactant material utilization, and the distribution of the electric current and voltages within the cathode. In addition, by representing the reactant particles as electrolyte or galvanic microcells governed by the Butler-Volmer electrochemical equation, the overall performance of battery cells during charge or discharge processes, respectively, can be predicted for a given cathodic powder composition. The presented model, executed on a high-performance computing architecture, essentially provides a valuable guidance in designing and developing future solid-state batteries.

solid-state battery↗
NASA NTRS1987

Solution of the two-dimensional Navier-Stokes equations using sparse matrix solvers

The use of direct sparse matrix solvers in the solution of the Navier-Stokes equations is investigated. The Yale Sparse Matrix Package and its implementation in the solution algorithm is described. The streamfunction-vorticity form of the Navier-Stokes equations are discretized and linearized and the resulting system of equations are solved using this package. Several viscous flow problems are investigated, including flow in a cavity and flow around a NACA0012 airfoil. Massively separated flow around a sine wave airfoil is investigated and high Reynolds number solutions are obtained. A solution of the unsteady flow around a Joukowski airfoil at high angle of attack is presented.

Bender, Erich E.↗
NASA NTRS2017

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Rose, Geoffrey K.↗
NASA NTRS1949

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗
NASA NTRS1992

Towards the development of micromechanics equations for ceramic matrix composites via fiber substructuring

A generic unit cell model which includes a unique fiber substructuring concept is proposed for the development of micromechanics equations for continuous fiber reinforcement ceramic composites. The unit cell consists of three constituents: fiber, matrix, and an interphase. In the present approach, the unit cell is further subdivided into several slices and the equations of micromechanics are derived for each slice. These are subsequently integrated to obtain ply level properties. A stand alone computer code containing the micromechanics model as a module is currently being developed specifically for the analysis of ceramic matrix composites. Towards this development, equivalent ply property results for a SiC/Ti-15-3 composite with 0.5 fiber volume ratio are presented and compared with those obtained from customary micromechanics models to illustrate the concept. Also, comparisons with limited experimental data for the ceramic matrix composite, SiC/RBSN (Reaction Bonded Silicon Nitride) with a 0.3 fiber volume ratio are given to validate the concepts.

Murthy, P. L. N.↗