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,135 records · Page 63

Structure and creep rupture properties of directionally solidified eutectic gamma/gamma-prime-alpha alloy

A simple ternary gamma/gamma-prime-alpha alloy of nominal composition (wt-%) Ni-32Mo-6Al has been directionally solidified at 17 mm/h and tested in creep rupture at 1073, 1173, and 1273 K. A uniform microstructure consisting of square-shaped Mo fibers in a gamma + gamma-prime matrix was found despite some variation in the molybdenum and aluminum concentrations along the growth direction. Although the steady-state creep rate is well described by the normal stress temperature equation, the stress exponent (12) and the activation energy (580 kJ/mol) are high. The rupture behavior is best characterized by the Larson-Miller parameter where the constant equals 20.

Whittenberger, J. D.↗

Inviscid steady/unsteady flow calculations

The solution of the Euler equations using a flux splitting procedure is considered for low subsonic to high supersonic flows. Steady and unsteady, internal and external flow fields, are computed. For transient flows, a direct sparse matrix solver is applied to compute the flow field at each instant of time. Oscillation free normal and oblique shocks are captured. Unstart and restart of a simplified two-dimensional inlet is investigated.

Pordal, H. S.↗

Cratering and penetration experiments in teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility after the spacecraft spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments to reproduce such features and to understand the relationships between projectile size and the resulting crater or penetration hole diameter over a wide range of impact velocities. Such relationships are needed to derive the size and mass frequency distribution and flux of natural and man-made particles in low-earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres into pure Teflon targets at velocities ranging from 1 to 7 km/s. Target thickness varied over more than three orders of magnitude from finite halfspace targets to very thin films. Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration hole diameter can become larger than that of a standard crater. The crater diameter in infinite halfspace Teflon targets increases, at otherwise constant impact conditions, with encounter velocity by a factor of V (exp 0.44). In contrast, the penetration hole size in very thin foils is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations apply to all impacts in which the shock-pulse duration of the projectile is shorter than that assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Horz, Friedrich↗

Multi-Dimensional Asymptotically Stable 4th Order Accurate Schemes for the Diffusion Equation

An algorithm is presented which solves the multi-dimensional diffusion equation on co mplex shapes to 4th-order accuracy and is asymptotically stable in time. This bounded-error result is achieved by constructing, on a rectangular grid, a differentiation matrix whose symmetric part is negative definite. The differentiation matrix accounts for the Dirichlet boundary condition by imposing penalty like terms. Numerical examples in 2-D show that the method is effective even where standard schemes, stable by traditional definitions fail.

Abarbanel, Saul↗

Stochastic Optimal Control via Bellman's Principle

This paper presents a method for finding optimal controls of nonlinear systems subject to random excitations. The method is capable to generate global control solutions when state and control constraints are present. The solution is global in the sense that controls for all initial conditions in a region of the state space are obtained. The approach is based on Bellman's Principle of optimality, the Gaussian closure and the Short-time Gaussian approximation. Examples include a system with a state-dependent diffusion term, a system in which the infinite hierarchy of moment equations cannot be analytically closed, and an impact system with a elastic boundary. The uncontrolled and controlled dynamics are studied by creating a Markov chain with a control dependent transition probability matrix via the Generalized Cell Mapping method. In this fashion, both the transient and stationary controlled responses are evaluated. The results show excellent control performances.

Crespo, Luis G.↗

Engineering and Scientific Applications: Using MatLab(Registered Trademark) for Data Processing and Visualization

MatLab(TradeMark)(MATrix LABoratory) is a numerical computation and simulation tool that is used by thousands Scientists and Engineers in many countries. MatLab does purely numerical calculations, which can be used as a glorified calculator or interpreter programming language; its real strength is in matrix manipulations. Computer algebra functionalities are achieved within the MatLab environment using "symbolic" toolbox. This feature is similar to computer algebra programs, provided by Maple or Mathematica to calculate with mathematical equations using symbolic operations. MatLab in its interpreter programming language form (command interface) is similar with well known programming languages such as C/C++, support data structures and cell arrays to define classes in object oriented programming. As such, MatLab is equipped with most of the essential constructs of a higher programming language. MatLab is packaged with an editor and debugging functionality useful to perform analysis of large MatLab programs and find errors. We believe there are many ways to approach real-world problems; prescribed methods to ensure foregoing solutions are incorporated in design and analysis of data processing and visualization can benefit engineers and scientist in gaining wider insight in actual implementation of their perspective experiments. This presentation will focus on data processing and visualizations aspects of engineering and scientific applications. Specifically, it will discuss methods and techniques to perform intermediate-level data processing covering engineering and scientific problems. MatLab programming techniques including reading various data files formats to produce customized publication-quality graphics, importing engineering and/or scientific data, organizing data in tabular format, exporting data to be used by other software programs such as Microsoft Excel, data presentation and visualization will be discussed.

Sen, Syamal K.↗

Solar system applications of Mie theory and of radiative transfer of polarized light

A theory of the multiple scattering of polarized light is discussed using the doubling method of van de Hulst. The concept of the Stokes parameters is derived and used to develop the form of the scattering phase matrix of a single particle. The diffuse reflection and transmission matrices of a single scattering plane parallel atmosphere are expressed as a function of the phase matrix, and the symmetry properties of these matrices are examined. Four matrices are required to describe scattering and transmission. The scattering matrix that results from the addition of two identical layers is derived. Using the doubling method, the scattering and transmission matrices of layers of arbitrary optical thickness can be derived. The doubling equations are then rewritten in terms of their Fourier components. Computation time is reduced since each Fourier component doubles independently. Computation time is also reduced through the use of symmetry properties.

Whitehill, L. P.↗

A math model for high velocity sensoring with a focal plane shuttered camera.

A new mathematical model is presented which describes the image produced by a focal plane shutter-equipped camera. The model is based upon the well-known collinearity condition equations and incorporates both the translational and rotational motion of the camera during the exposure interval. The first differentials of the model with respect to exposure interval, delta t, yield the general matrix expressions for image velocities which may be simplified to known cases. The exposure interval, delta t, may be replaced under certain circumstances with a function incorporating blind velocity and image position if desired. The model is tested using simulated Lunar Orbiter data and found to be computationally stable as well as providing excellent results, provided that some external information is available on the velocity parameters.

Morgan, P.↗

On-orbit Systems Identification of Flexible Spacecraft Technol., 1984, Pt. 2 p 465-481

Future spacecraft include configurations which are too flexible to be adequately tested prior to flight and which will require on-orbit systems identification to ensure safe operation of the flight control system. The structural dynamics model will evolve and its accuracy will improve in stages as ground tests of full-scale components and replica-scale models are performed. State Space Modeling and Conditional Maximum Likelihood Parameter Estimation methodology can provide the formal probability-based framework for the process of upgrading a model as additional test results are obtained. Although the number of unknown parameters can be reduced by the use of canonical forms for the stability matrix, the number of unknown model parameters quickly becomes unmanageable unless advantage is taken of the relationship of a much fewer number of global model parameters. Distributed parameter systems or partial differential equation models are one way to take advantage of such global parameters to reduce the number of unknown model parameters.

Taylor, L.↗

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗

Methods for the solution of radiative transfer equation

To obtain an exact solution of the radiative-transfer equation in media where both absorption and scattering are significant, the usual approach is to use a numerical method. Three methods are known in the literature: invariant imbedding, eigenvalue-eigenfunction, and matrix doubling. This paper examines the practical application of these methods to the problem of emission from an inhomogeneous (Rayleigh) layer, the effects of layer parameters on the stability. It is found that invariant imbedding is most suitable for computing emission from an inhomogeneous layer with a temperature profile but tends to be unstable as the optical thickness of the layer increases beyond 0.5. On the other hand, the matrix-doubling method is stable for arbitrary optical thickness but is not suitable for handling multilayers. The eigenvalue-eigenfunction method is more stable than the invariant imbedding as optical thickness increases up to 2.0. It also permits temperature profile in the layer, but the computation is much more complicated. It is less stable than the matrix-doubling method when optical thickness is larger than 2.0. In general, the choice of a method is dependent on the nature of the problem.

Chen, M. F.↗

Cratering and penetration experiments in Teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility (LDEF) after the spacecraft had spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments in an effort to reproduce such features and to -- hopefully -- understand the relationships between projectile size and the resulting crater or penetration-hole diameter over a wide range of impact velocity. Such relationships are needed to derive the size- and mass-frequency distribution and flux of natural and man-made particles in low-Earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres of 3.175 mm (1/8 inch) nominal diameter (Dp) into pure Teflon FEP targets at velocities ranging from 1 to 7 km/s. Target thickness (T) was varied over more than three orders of magnitude from infinite halfspace targets (Dp/T less than 0.1) to very thin films (Dp/T greater than 100). Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration-hole diameter (Dh) can become larger than that of a standard crater (Dc) at relative target thicknesses of Dp/T = 0.6-0.9. The crater diameter is infinite halfspace Teflon targets increases -- at otherwise constant impact conditions -- with encounter velocity by a factor of V0.44. In contrast, the penetration-hole size is very thin foils (Dp/T greater than 50) is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for the velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations dominate all impacts in which the shock-pulse duration of the projectile (tp) is shorter than that of the target (tt). We also demonstrate that each penetration hole from space-retrieved surfaces may be assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Hoerz, Friedrich↗

Development of Implicit Methods in CFD NASA Ames Research Center 1970's - 1980's

The focus here is on the early development (mid 1970's-1980's) at NASA Ames Research Center of implicit methods in Computational Fluid Dynamics (CFD). A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The emphasis will be on the Euler equations. A review of material pertinent to the solution of the Euler equations within the framework of implicit methods will be presented. The eigensystem of the equations will be used extensively in developing a framework for various methods applied to the Euler equations. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented.

Pulliam, Thomas H.↗

Engineering and Scientific Applications: Using MatLab(Registered Trademark) for Data Processing and Visualization

MatLab(R) (MATrix LABoratory) is a numerical computation and simulation tool that is used by thousands Scientists and Engineers in many cou ntries. MatLab does purely numerical calculations, which can be used as a glorified calculator or interpreter programming language; its re al strength is in matrix manipulations. Computer algebra functionalities are achieved within the MatLab environment using "symbolic" toolbo x. This feature is similar to computer algebra programs, provided by Maple or Mathematica to calculate with mathematical equations using s ymbolic operations. MatLab in its interpreter programming language fo rm (command interface) is similar with well known programming languag es such as C/C++, support data structures and cell arrays to define c lasses in object oriented programming. As such, MatLab is equipped with most ofthe essential constructs of a higher programming language. M atLab is packaged with an editor and debugging functionality useful t o perform analysis of large MatLab programs and find errors. We belie ve there are many ways to approach real-world problems; prescribed methods to ensure foregoing solutions are incorporated in design and ana lysis of data processing and visualization can benefit engineers and scientist in gaining wider insight in actual implementation of their perspective experiments. This presentation will focus on data processing and visualizations aspects of engineering and scientific applicati ons. Specifically, it will discuss methods and techniques to perform intermediate-level data processing covering engineering and scientifi c problems. MatLab programming techniques including reading various data files formats to produce customized publication-quality graphics, importing engineering and/or scientific data, organizing data in tabu lar format, exporting data to be used by other software programs such as Microsoft Excel, data presentation and visualization will be discussed. The presentation will emphasize creating practIcal scripts (pro grams) that extend the basic features of MatLab TOPICS mclude (1) Ma trix and vector analysis and manipulations (2) Mathematical functions (3) Symbolic calculations & functions (4) Import/export data files (5) Program lOgic and flow control (6) Writing function and passing parameters (7) Test application programs

Sen, Syamal K.↗

Probabilistic simulation of long term behavior in polymer matrix composites

A methodology to compute cumulative probability distribution functions (CDF) of fatigue life for different ratios, r of applied stress to the laminate strength based on first ply failure criteria has been developed and demonstrated. Degradation effects due to long term environmental exposure and 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/aging of material properties due to cyclic loads. Fast probability integration method is used to perform probabilistic simulation of uncertainties. Sensitivity of fatigue life reliability to uncertainties in the primitive random variables are computed and their significance in the reliability based design for maximum life is discussed. The results show that the graphite/epoxy (0/+45/90) deg laminate with ply thickness 0.125 in. has 500,000 cycles life for applied stress to laminate strength ratio of 0.6 and a reliability of 0.999. Also, the fatigue life reliability has been found to be most sensitive to the ply thickness and matrix tensile strength. Tighter quality controls must therefore be enforced on ply thickness and matrix strength in order to achieve high reliability of the structure.

Shah, A. R.↗

On the factorization and fitting of molecular scattering information

The reported analysis is based on the factored IOS T-matrix. It is shown that line shape measurements may be used over a range of temperatures to evaluate inelastic scattering cross sections. Basic factorization or parameterization relations are derived by considering the wavefunction equations. The parameterization of cross sections is considered, taking into account the differential scattering amplitude and cross section, integral cross sections, phenomenological cross sections for general relaxation processes, and viscosity and diffusion cross sections. Thermal averages and rates are discussed, giving attention to integral cross sections and rates, and general phenomenological cross sections. The results of computational studies are also presented.

Goldflam, R.↗

Vibration analysis of rotor blades with pendulum absorbers

A comprehensive vibration analysis of rotor blades with spherical pendulum absorbers is presented. Linearized equations of motion for small oscillations about the steady-state deflection of a spherical pendulum on elastic rotor blades undergoing coupled flapwise bending, chordwise bending, and torsional vibrations are obtained. A transmission matrix formulation is given to determine the natural vibrational characteristics of rotor blades with spherical or simple flapping pendulum absorbers. The natural frequencies and mode shapes of a hingeless rotor blade with a spherical pendulum are computed.

Murthy, V. R.↗

Errors incurred in a plane-wave-type expansion of a Gaussian beam

The multipole expansion obtained by Morita et al. (1968) of the Gaussian laser beam used to levitate an aerosol particle in order that its complete phase matrix may be measured is compared with that of Tsai and Pogorzelski (1975) in order to demonstrate the effect of the incorrect expansion used by Morita. Errors incurred by the use of an equation in which one side satisfies the scalar wave equation while the other side does not and can be reduced to a plane wave amplitude are calculated as functions of the inverse of the wave number times the beam waist, the wave number times the radial spherical coordinate and the angular spherical coordinate. Errors on the order of a few percent, considered undetectable are obtained in the squared-field amplitudes due to the expansion, however, they are found to become significant (several tens of percent) when the angle is zero. It is concluded that the expansion of Morita should only be used in the regions where the spherical angle is less than 0.01 and its product with the wave number and the radial spherical coordinate is less than unity.

Kattawar, G. W.↗