Diabatic and adiabatic representations for atomic collision problems.
Coupled equations for heavy particle motion expressed by introducing generalized matrix operator for effective momentum in atomic collision problems
Smith, F. T.↗
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.
Transient, compressible heat and mass transfer in porous media using the strongly implicit iteration procedure.
Coupled nonlinear partial differential equations describing heat and mass transfer in a porous matrix are solved in finite difference form with the aid of a new iterative technique (the strongly implicit procedure). Example numerical results demonstrate the characteristics of heat and mass transport in a porous matrix such as a charring ablator. It is emphasized that multidimensional flow must be considered when predicting the thermal response of a porous material subjected to nonuniform boundary conditions.
Curry, D. M.↗
Hypermatrix scheme for finite element systems on CDC STAR-100 computer
A study is made of the adaptation of the hypermatrix (block matrix) scheme for solving large systems of finite element equations to the CDC STAR-100 computer. Discussion is focused on the organization of the hypermatrix computation using Cholesky decomposition and the mode of storage of the different submatrices to take advantage of the STAR pipeline (streaming) capability. Consideration is also given to the associated data handling problems and the means of balancing the I/Q and cpu times in the solution process. Numerical examples are presented showing anticipated gain in cpu speed over the CDC 6600 to be obtained by using the proposed algorithms on the STAR computer.
Noor, A. K.↗
A study of the effects of state transition matrix approximations
The effects of using an approximate state transition matrix in orbit estimation are investigated. The approximate state transition matrix results when higher order geopotential terms in the equations of motion are ignored in the formation of the variational equations. Two methods of orbit estimation are considered: the differential correction procedure (DC) and the extended Kalman filter (EKF). The system used for the study is the Research & Development version of the Goddard Trajectory Determination System. The effects of the approximation are analyzed on a number of orbits. These orbits of various inclinations and semimajor axes. Other parameters studied include geopotential models and DC arc length.
May, J. A.↗
Buckling of rotating beams
The stability of a beam subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. The differential equations of motion are solved numerically using an integrating matrix method in combination with an eigenanalysis to determine the eigenvalues from which stability is assessed. The results clarify several differences which have been identified in the literature relating to the proper behavior of the critical rotational speed for buckling as the radius of rotation of the clamped end of the beam is reduced.
White, W. F., Jr.↗
Multiple scattering effects on wave propagation due to rain
The paper examines the multiple scattering of microwaves due to rain and the resulting incoherent intensity for a plane wave incident on a plane-parallel rain region. General formulations of the transfer equation using Stokes parameters are presented, and an extinction matrix which takes into account depolarization effects and nonspherical droplet shape is introduced. Scattering characteristics are calculated for spherical droplets using the Mie solution and the Laws-Parsons distributions. The ratios of incoherent to coherent received powers are calculated at 30 GHz for rain rates of 12.5, 50, 100, and 150 mm/h, rain thicknesses of 3 km and 10 km, and fields of view of 1.5, 5, and 15 deg.
Ishimaru, A.↗
Kinematic properties of rotary-wing and fixed-wing aircraft in steady coordinated high-g turns
An analytical approach to the study of flight dynamics of aircraft operating in a high-angle-of-attack flight regime and of helicopters operating in extreme thrust conditions is presented. Steady coordinated high-g turns are used to establish the initial equilibrium flight conditions near stall angles of attack. The kinematic properties of the aircraft in steady coordinated turns are examined: in high-g turns, pitch rate (independent of the angle of attack) is of a much larger magnitude than roll and yaw rate; a substantial roll rate is found to develop in steep turns for all angles of attack; the angle of attack also has a significant effect on the pitch attitude, with decreasing influence as the normal load factor increases. The exact small disturbance equations of motion of the aircraft in general steady turns are also developed for application to both rotary-wing and fixed-wing aircraft in extreme conditions. These equations are in a first-order, vector-matrix format, and are thus compatible with many efficient software packages developed in modern system theory.
Chen, R. T. N.↗
Determining 3-D motion and structure from image sequences
A method of determining three-dimensional motion and structure from two image frames is presented. The method requires eight point correspondences between the two frames, from which motion and structure parameters are determined by solving a set of eight linear equations and a singular value decomposition of a 3x3 matrix. It is shown that the solution thus obtained is unique.
Huang, T. S.↗
Kalman filter modeling
The formulation of appropriate state-space models for Kalman filtering applications is studied. The so-called model is completely specified by four matrix parameters and the initial conditions of the recursive equations. Once these are determined, the die is cast, and the way in which the measurements are weighted is determined foreverafter. Thus, finding a model that fits the physical situation at hand is all important. Also, it is often the most difficult aspect of designing a Kalman filter. Formulation of discrete state models from the spectral density and ARMA random process descriptions is discussed. Finally, it is pointed out that many common processes encountered in applied work (such as band-limited white noise) simply do not lend themselves very well to Kalman filter modeling.
Brown, R. G.↗
Calculation of compressible flow about three-dimensional inlets with auxiliary inlets, slats and vanes by means of a panel method
An efficient and user oriented method was constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: (1) the use of a panel method; (2) a technique of superposition for obtaining solutions at any inlet operating condition; and (3) employment of an advanced matrix iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that yield a more complete comprehension of the flowfield than was possible previously.
Hess, J. L.↗
Calculation of compressible flow about three-dimensional inlets with auxiliary inlets, slats and vanes by means of a panel method
An efficient and user-oriented method has been constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: the use of a panel method, a technique of superposition for obtaining solutions at any inlet operating condition, and employment of an advanced matrix-iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that, taken together, yield a more complete comprehension of the flowfield than had been possible previously. Examples of these features are presented for some complicated configurations, and where possible, comparisons are made between calculation and experiment.
Hess, J. L.↗
Development of improved gravity field estimates using existing altimeter data
The formulation of a Fourier series and geographic representation expressions for orbit error were reported. Addition activities include the study of undulation errors in the Lagrangian formulation, and the development and testing of all the necessary software for the above expression. Additionally, the matrix formulations that led to the formation of observed equations and covariance propagation were developed.
Source record↗
Compliance matrices for cracked bodies
An algorithm is developed to construct the compliance matrix for a cracked solid in the integral-equation formulation of two-dimensional linear-elastic fracture mechanics. The integral equation is reduced to a system of algebraic equations for unknown values of the dislocation-density function at discrete points on the interval from -1 to 1, using the numerical procedure described by Gerasoulis (1982). Sample numerical results are presented, and it is suggested that the algorithm is especially useful in cases where iterative solutions are required; e.g., models of fiber-reinforced concrete, rocks, or ceramics where microcracking, fiber bridging, and other nonlinear effects are treated as nonlinear springs along the crack surfaces (Ballarini et al., 1984).
Ballarini, R.↗
Determination of settings of a tilted head-cutter for generation of hypoid and spiral bevel gears
Kinematics of Gleason mechanisms of hypoid and spiral bevel cutting machines are considered. These mechanisms are designated to install the position and tilt of the head cutter. The tilt of the head cutter with standard blades provides the required pressure angle. The authors have developed the matrix presentation of kinematics of these meachanisms and basic equations for the required settings. An example is presented based on the developed computation procedure.
Litvin, F. L.↗
Krylov subspace methods on supercomputers
A short survey of recent research on Krylov subspace methods with emphasis on implementation on vector and parallel computers is presented. Conjugate gradient methods have proven very useful on traditional scalar computers, and their popularity is likely to increase as three-dimensional models gain importance. A conservative approach to derive effective iterative techniques for supercomputers has been to find efficient parallel/vector implementations of the standard algorithms. The main source of difficulty in the incomplete factorization preconditionings is in the solution of the triangular systems at each step. A few approaches consisting of implementing efficient forward and backward triangular solutions are described in detail. Polynomial preconditioning as an alternative to standard incomplete factorization techniques is also discussed. Another efficient approach is to reorder the equations so as to improve the structure of the matrix to achieve better parallelism or vectorization. An overview of these and other ideas and their effectiveness or potential for different types of architectures is given.
Saad, Youcef↗
Parallel processing techniques for finite element analysis of nonlinear large truss structures
Methods were developed for parallel processing of finite element solutions of large truss structures. The parallel processing techniques were implemented in two stages, i.e., the repeated forming of the nonlinear global stiffness matrix and the solving of the global system of equations. The Sequent Balance 21000 parallel computer was employed to demonstrate the procedures and the speed-up.
Chien, L. S.↗
A new technique for simulating composite material. Task 2: Analytical solutions with Generalized Impedance Boundary Conditions (GIBCs)
The diffraction problem associated with a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GSMF) in conjunction with the dual integral equation approach. The multilayer slab is replaced by a surface obeying a generalized impedance boundary condition (GIBC) to facilitate the computation of the pertinent Wiener Hopf split functions and their zeros. Both E(sub z) and H(sub z) polarizations are considered and a number of scattering patterns are presented, some of which are compared to exact results available for a homogeneous recessed slab.
Ricoy, M. A.↗
Electromagnetic scattering from two-dimensional thick material junctions
The problem of the plane wave diffraction is examined by an arbitrary symmetric two dimensional junction, where Generalized Impedance Boundary Conditions (GIBCs) and Generalized Sheet Transition Conditions (GSTCs) are employed to simulate the slabs. GIBCs and GSTCs are constructed for multilayer planar slabs of arbitrary thickness and the resulting GIBC/GSTC reflection coefficients are compared with exact counterparts to evaluate the GIBCs/GSTCs. The plane wave diffraction by a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GDMF) in conjunction with the dual integral equation approach. Various scattering patterns are computed and validated with exact results where possible. The diffraction by a material discontinuity in a thick dielectric/ferrite slab is considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of unknown constants is obtained, and these constants are evaluated for the recessed slab geometry by comparison with the solution obtained therein. Several other simplified cases are also presented and discussed. An eigenfunction expansion method is introduced to determine the unknown solution constants in the general case. This procedure is applied to the non-unique solution in terms of unknown constants; and scattering patterns are presented for various slab junctions and compared with alternative results where possible.
Ricoy, M. A.↗