Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Optimization of multiple flyby trajectories

A procedure has been developed which minimizes total delta-V (instantaneous velocity change) for a multiple flyby trajectory with constraints on flyby altitude and orientation. The solution is found by varying the locations of maneuver points between each flyby to minimize the delta-Vs at the maneuver points. Each trajectory segment connecting consecutive maneuver points is found by solving an N-body analog to Lambert's problem. Multiconic techniques are used for trajectory propagation and for computation of the state transition matrix. The constrained parameter optimization problem is converted to an unconstrained problem by means of penalty functions and then solved with a quasi-Newton algorithm utilizing analytic first derivatives. This procedure has been successfully applied to Galileo satellite tour trajectories.

Damario, L. A.↗
NASA NTRS1989

Multiple boundary condition testing error analysis

Techniques for interpreting data from multiple-boundary-condition (MBC) ground tests of large space structures are developed analytically and demonstrated. The use of MBC testing to validate structures too large to stand alone on the ground is explained; the generalized least-squares mass and stiffness curve-fitting methods typically applied to MBC test data are reviewed; and a detailed error analysis is performed. Consideration is given to sensitivity coefficients, covariance-matrix theory, the correspondence between test and analysis modes, constraints and step sizes, convergence criteria, and factor-analysis theory. Numerical results for a simple beam problem are presented in tables and briefly characterized. The improved error-updating capabilities of MBC testing are confirmed, and it is concluded that reasonably accurate results can be obtained using a diagonal covariance matrix.

Glaser, R. J.↗
NASA NTRS2007

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" (DR) method is introduced for analyzing lightning field changes (DELTAEs) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an Overdetermined Fixed Matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all 8 parameters (x, y, z, Q, x', y', z', Q') are found, yet a numerical search of only 4 parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused DeltaEs derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the Lightning Detection And Ranging (LDAR) system at the KSC-ER, and from least squares error estimation theory, and the method is shown to be a useful "stand-alone" charge retrieval tool. Since more than one charge distribution describes a finite set of DELTAEs (i.e., solutions are non-unique), and since there can exist appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence, an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions, thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning in several Florida storms.

Koshak, W. J.↗
NASA NTRS2007

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" ("DR") method is introduced for analyzing lightning field changes ((Delta)Es) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an overdetermined fixed matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all eight parameters Ix. y, z, Q, x', y', z', Q') are found, yet a numerical search of only four parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused (Delta)Es derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and U.S. Air Force Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the lightning detection and ranging (LDAR) system at the KSC-ER. and from least squares error estimation theory, and the method is shown to be a useful "stand alone" charge retrieval tool. Since more than one charge distribution describes a finite set of (Delta)Es (i.e., solutions are nonunique), and since there can be appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence. an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions. thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning m several Florida storms.

Koshak, W. J.↗
NASA NTRS2009

Self-Assembling, Flexible, Pre-Ceramic Composite Preforms

In this innovation, light weight, high temperature, compact aerospace structures with increased design options are made possible by using self-assembling, flexible, pre-ceramic composite materials. These materials are comprised of either ceramic or carbon fiber performs, which are infiltrated with polymer precursors that convert to ceramics upon thermal exposure. The preform architecture can vary from chopped fibers formed into blankets or felt, to continuous fibers formed into a variety of 2D or 3D weaves or braids. The matrix material can also vary considerably. For demonstration purposes, a 2D carbon weave was infiltrated with a SiC polymer precursor. The green or unfired material is fabricated into its final shape while it is still pliable. It is then folded or rolled into a much more compact shape, which will occupy a smaller space. With this approach, the part remains as one continuous piece, rather than being fabricated as multiple sections, which would require numerous seals for eventual component use. The infiltrated preform can then be deployed in-situ. The component can be assembled into its final shape by taking advantage of the elasticity of the material, which permits the structure to unfold and spring into its final form under its own stored energy. The pre-ceramic composites are converted to ceramics and rigidized immediately after deployment. The final ceramic composite yields a high-temperature, high-strength material suitable for a variety of aerospace structures. The flexibility of the material, combined with its high-temperature structural capacity after rigidization, leads to a less complex component design with an increased temperature range. The collapsibility of these structures allows for larger components to be designed and used, and also offers the potential for increased vehicle performance. For the case of collapsible nozzle extensions, a larger nozzle, and thus a larger nozzle exit plane, is possible because interference with surrounding structures can be avoided in the collapsed state. The larger exit plane leads to an increase in expansion area ratio, which has the potential to increase thrust and overall rocket performance. In general, the use of advanced ceramic materials can lead to improved engine and vehicle performance. The ceramics can run hotter, so less cooling is required. Fuel to coolant ratios can be balanced more readily to reduce weight. Engine efficiency can also be increased with hotter combustion and exhaust temperatures. In addition, the ceramic composites themselves can reduce the component weight by as much as 50 percent, which can translate into greater payload for the vehicle

Jaskowiak, Martha H.↗
NASA NTRSDate not supplied

Evaluating Crystallinity in Thermoplastic Composites for Aerospace Applications

Polymer matrix composites (PMCs) offer many benefits for the aerospace industry due to their potential for weight reduction when compared to metal or ceramic based materials. Most PMCs currently in flight use thermoset matrices, however, thermoplastic resins are being explored as alternatives due to their ability to be remelted, which is of particular interest due to the potential for in-situ repair and faster production. Most thermoplastic resins are semicrystalline polymers. The properties of semicrystalline thermoplastics are directly influenced by their crystallinity, which can vary due to many factors including thermal treatments, environmental conditions, and mechanical deformation. Within thermoplastic PMC parts, crystallinity gradients can arise due to variations in part geometry, across part thicknesses, and along bonded joints. Monitoring the crystallinity of thermoplastic composites is key to ensuring these materials meet the high demands required for aerospace. Several different analytical techniques exist that can be used to characterize the bulk crystallinity of thermoplastic materials. However, many existing methods lack the specificity required to identify the subtle variations in crystallinity that may play a significant role in the performance and durability of PMC parts. Because of this, a significant amount of work is still required to fully characterize and understand the crystallinity profiles of thermoplastic PMCs and the resulting impact to material properties. This talk discusses the use of multiple techniques such as Differential Scanning Calorimetry, Polarized Light Optical Microscopy, and Fourier-Transform Infrared Spectroscopy to characterize the crystallinity in carbon fiber/thermoplastic composites. Samples of different crystallinity profiles were manufactured using various cooling procedures. This work aims to provide the fundamental data necessary to understand the effects of crystallinity on thermoplastic PMCs, which is key to advancing their use in aerospace applications.

Thermoplastics↗
NASA NTRS1982

Cure kinetics of epoxy matrix resin by differential scanning calorimetry

A study was made on the cure kinetics of an epoxy neat-resin (Narmco 5208) using Differential Scanning Calorimetry (DSC). Two interrelated analytical methods were applied to dynamic DSC data for evaluating the kinetic parameters, such as activation energy, E, the order of reaction, n, and the total heat of polymerization (or crosslinking), delta H sub t. The first method was proposed by Ellerstein (1968), and uses a thorough differential-integral analysis of a single DSC curve to evaluate the kinetic parameters. The second method was proposed by Kissinger (1957), and uses multiple DSC curves obtained at various heating rates to evaluate E regardless of n. Kinetic analysis of Narmco 5208 epoxy resin showed that the reaction order, n, is substantially affected by the rate of heating; i.e., n is approximately 2 at slow scan rates but is reduced to 1.5 at higher scan rates. The activation energy, E, is not affected by the scan rate, and the average value of E is 25.6 + or - 1.8 kcal/mole.

Cizmecioglu, M.↗
NASA NTRS1989

Internal object loading for multiple cooperating robot manipulators

For an object being rigidly grasped and manipulated by multiple robotic mechanisms, the internal loading characteristics at a common coordinate set within the object are considered. It is demonstrated that representation of internal forces and moments in these common coordinates gives insight into force and load distribution schemes developed previously. In particular, it is shown how internal loads may be created in some end-effector force distribution schemes even when no component in the null-space of the grasp matrix is included. It is further shown that a particular pseudoinverse of the grasp matrix, which can be shown to be consistent with the kinematic constraints, may be used to eliminate this situation. The case of a two-arm system is used to illustrate the concepts introduced.

Walker, Ian D.↗
NASA NTRS1989

Comparison of joint space versus task force load distribution optimization for a multiarm manipulator system

It is often proposed that the redundancy in choosing a force distribution for multiple arms grasping a single object should be handled by minimizing a quadratic performance index. The performance index may be formulated in terms of joint torques or in terms of the Cartesian space force/torque applied to the body by the grippers. The former seeks to minimize power consumption while the latter minimizes body stresses. Because the cost functions are related to each other by a joint angle dependent transformation on the weight matrix, it might be argued that either method tends to reduce power consumption, but clearly the joint space minimization is optimal. A comparison of these two options is presented with consideration given to computational cost and power consumption. Simulation results using a two arm robot system are presented to show the savings realized by employing the joint space optimization. These savings are offset by additional complexity, computation time and in some cases processor power consumption.

Soloway, Donald I.↗
NASA NTRS1991

Polarimetric signatures of a layer of random nonspherical discrete scatterers overlying a homogeneous half-space based on first- and second-order vector radiative transfer theory

Complete polarimetric signatures of a layer of random, nonspherical discrete scatterers overlying a homogeneous half space are studied with the first- and second-order solutions of the vector radiative transfer theory. Some of the salient features of the numerical results are as follows: (1) the inclusion of the nondiagonal extinction matrix in the vector radiative transfer theory accounts for an appreciable phase difference between vv and hh polarizations, particularly for aligned scatterers; (2) the ensemble-averaged scattered Stokes vector is generally partially polarized, with the degree of polarization less than unity; (3) there generally exists a pedestal in the copolarization return when plotted as a function of ellipticity and orientation angles, which may be due to heterogeneity of scattering objects and/or multiple scattering effects; and (4) multiple scattering effects generally enhance the pedestal in copolarization return, decrease the degree of polarization, affect phase difference, and also enhance the depolarization return.

Tsang, Leung↗
NASA NTRS2013

The Effect of Scale Dependent Discretization on the Progressive Failure of Composite Materials Using Multiscale Analyses

A multiscale modeling methodology, which incorporates a statistical distribution of fiber strengths into coupled micromechanics/ finite element analyses, is applied to unidirectional polymer matrix composites (PMCs) to analyze the effect of mesh discretization both at the micro- and macroscales on the predicted ultimate tensile (UTS) strength and failure behavior. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a PMC tensile specimen that initiates at the repeating unit cell (RUC) level. Three different finite element mesh densities were employed and each coupled with an appropriate RUC. Multiple simulations were performed in order to assess the effect of a statistical distribution of fiber strengths on the bulk composite failure and predicted strength. The coupled effects of both the micro- and macroscale discretizations were found to have a noticeable effect on the predicted UTS and computational efficiency of the simulations.

Ricks, Trenton M.↗
NASA NTRS1991

A 640-MHz 32-megachannel real-time polyphase-FFT spectrum analyzer

A polyphase fast Fourier transform (FFT) spectrum analyzer being designed for NASA's Search for Extraterrestrial Intelligence (SETI) Sky Survey at the Jet Propulsion Laboratory is described. By replacing the time domain multiplicative window preprocessing with polyphase filter processing, much of the processing loss of windowed FFTs can be eliminated. Polyphase coefficient memory costs are minimized by effective use of run length compression. Finite word length effects are analyzed, producing a balanced system with 8 bit inputs, 16 bit fixed point polyphase arithmetic, and 24 bit fixed point FFT arithmetic. Fixed point renormalization midway through the computation is seen to be naturally accommodated by the matrix FFT algorithm proposed. Simulation results validate the finite word length arithmetic analysis and the renormalization technique.

Zimmerman, G. A.↗
NASA NTRS1993

Theoretical calculations on the electron absorption spectra of selected Polycyclic Aromatic Hydrocarbons (PAH) and derivatives

As a theoretical component of the joint effort with the laboratory of Dr. Lou Allamandola to search for potential candidates for interstellar organic carbon compound that are responsible for the visible diffuse interstellar absorption bands (DIB's), quantum mechanical calculations were performed on the electron absorption spectra of selected polycyclic aromatic hydrocarbons (PAH) and derivatives. In the completed project, 15 different species of naphthalene, its hydrogen abstraction and addition derivatives, and corresponding cations and anions were studied. Using semiempirical quantum mechanical method INDO/S, the ground electronic state of each species was evaluated with restricted Hartree-Fock scheme and limited configuration interaction. The lowest energy spin state for each species was used for electron absorption calculations. Results indicate that these calculations are accurate enough to reproduce the spectra of naphthalene cation and anion observed in neon matrix. The spectral pattern of the hydrogen abstraction and addition derivatives predicted based on these results indicate that the electron configuration of the pi orbitals of these species is the dominant determinant. A combined list of 19 absorptions calculated from 4500 A to 10,400 A were compiled and suggested as potential candidates that are relevant for the DIB's absorptions. Continued studies on pyrene and derivatives revealed the ground state symmetries and multiplicities of its neutral, anionic, and cationic species. Spectral calculations show that the cation (B(sub 3g)-2) and the anion (A(sub u)-2) are more likely to have low energy absorptions in the regions between 10 kK and 20 kK, similar to naphthalene. These absorptions, together with those to be determined from the hydrogen abstraction and addition derivatives of pyrene, can be used to provide additional candidates and suggest experimental work in the search for interstellar compounds that are responsible for DIB's.

Du, Ping↗
NASA NTRS1993

Internal strain analysis of ceramics using scanning laser acoustic microscopy

Quantitative studies of material behavior characteristics are essential for predicting the functionality of a material under its operating conditions. A nonintrusive methodology for measuring the in situ strain of small dimeter (to 11 microns) ceramic fibers under uniaxial tensile loading and the local internal strains of ceramics and ceramic composites under flexural loading is introduced. The strain measurements and experimentally observed mechanical behavior are analyzed in terms of the microstructural development and fracture behavior of each test specimen evaluated. Measurement and analysis of Nicalon silicon carbide (SiC) fiber (15 microns diameter) indicate that the mean elastic modulus of the individual fiber is 185.3 GPa. Deviations observed in the experimentally determined elastic modulus values between specimens were attributed to microstructural variations which occur during processing. Corresponding variations in the fracture surface morphology were also observed. The observed local mechanical behavior of a lithium alumino-silicate (LAS) glass ceramic, a LAS/SiC monofilament composite, and a calcium alumino-silicate (CAS)/SiC fully reinforced composite exhibits nonlinearities and apparent hysteresis due to the subcritical mechanical loading. Local hysteresis in the LAS matrices coincided with the occurrence of multiple fracture initiation sites, localized microcracking, and secondary cracking. The observed microcracking phenomenon was attributed to stress relaxation of residual stresses developed during processing, and local interaction of the crack front with the microstructure. The relaxation strain and stress predicted on apparent mechanical hysteresis effects were defined and correlated with the magnitude of the measured fracture stress for each specimen studied. This quantitative correlation indicated a repeatable measure of the stress at which matrix microcracking occurred for stress relief of each material system. Stress relaxation occurred prior to the onset of steady state cracking conditions. The relaxation stress occurred at 18.5 percent of the fracture stress in LAS and 11.0 percent of the yield stress in CAS/SiC. The relaxation stress ratio was dependent upon the dominant fracture mode of the LAS/SiC specimens. Relaxation stress ratios greater than 0.30 were observed for specimens which fractured due to shear at the fiber matrix interface; specimens which fracture due to tensile cracking had relaxation stress ratios less than 0.30. The stress relaxation ratio appeared to be a specific characteristic of the glass ceramic material. The measured stress relaxation for LAS indicated a measure of the inherent residual stresses in the material due to processing and suggested localized toughening mechanisms for brittle material structures.

Kent, Renee M.↗
NASA NTRS2007

Quaternion Averaging

Many applications require an algorithm that averages quaternions in an optimal manner. For example, when combining the quaternion outputs of multiple star trackers having this output capability, it is desirable to properly average the quaternions without recomputing the attitude from the the raw star tracker data. Other applications requiring some sort of optimal quaternion averaging include particle filtering and multiple-model adaptive estimation, where weighted quaternions are used to determine the quaternion estimate. For spacecraft attitude estimation applications, derives an optimal averaging scheme to compute the average of a set of weighted attitude matrices using the singular value decomposition method. Focusing on a 4-dimensional quaternion Gaussian distribution on the unit hypersphere, provides an approach to computing the average quaternion by minimizing a quaternion cost function that is equivalent to the attitude matrix cost function Motivated by and extending its results, this Note derives an algorithm that deterniines an optimal average quaternion from a set of scalar- or matrix-weighted quaternions. Rirthermore, a sufficient condition for the uniqueness of the average quaternion, and the equivalence of the mininiization problem, stated herein, to maximum likelihood estimation, are shown.

Markley, F. Landis↗
NASA NTRSDate not supplied

Coordinated analyses of hydrated Kuiper Belt dust

We report on coordinated mineralogical, chemical, and isotopic analyses of hydrated carbon-rich interplanetary dust particles (IDPs). We focus on hydrated IDPs that contain high solar flare track densities consistent with an origin in outer Solar System bodies such as Edgeworth-Kuiper belt objects (EKBOs) [1]. We are searching for evidence of the heavy water reservoir predicted by self-shielding models [2]. The IDPs were embedded in sulfur and microtome thin sections were obtained for Fourier-transform infrared (FTIR), transmission electron microscopy (TEM), and nanoSIMS analyses. The IDPs were only partly microtomed and the remainder was extracted from the sulfur and pressed into In or Au for high precision O isotopic analyses. The hydrated IDPs in this study are compact low-porosity objects with a mineralogy dominated by fine-grained saponite and lesser serpentine, Fe-Ni sulfides, Mg-Fe carbonates, rare crystaline silicates, and abundant organic matter. The C in the IDPs is abundant (typically >4X CI) and occurs as isolated nanoglobules and finely disseminated within the phyllosilicate matrix of the IDPs. TEM analyses of the nanoglobules show O-rich, S-depleted rims with higher carbonyl functionality than their cores. FTIR spectra show strong aliphatic C-H stretches consistent with the presence of short-chain aliphatic hydrocarbons in these IDPs. The oxygen isotopic compositions measured in 6 IDPs all plot within error of the CCAM line [3]. For several of the IDPs we were able to obtain multiple analyses. The ∆^(17)O values range from -23‰ up to +3.4‰ with most analyses clustering above the terrestrial fractionation line. The O isotopic compositions and the high C contents of these IDPs are distinct from known hydrated carbonaceous chondrites. The track-rich hydrated IDPs record parent body aqueous alteration effects through interactions with 16O-poor H2O on a EKBO parent body. The aqueous alteration activity likely resulted from the heat generated by collisional processes among EKBOs given that these objects are believed to have accreted late after the decay of 26Al [4].

Lindsay P. Keller↗
NASA NTRS1992

A comparison between theoretical prediction and experimental measurement of the dynamic behavior of spur gears

A comparison was made between computer model predictions of gear dynamics behavior and experimental results. The experimental data were derived from the NASA gear noise rig, which was used to record dynamic tooth loads and vibration. The experimental results were compared with predictions from the DSTO Aeronautical Research Laboratory's gear dynamics code for a matrix of 28 load speed points. At high torque the peak dynamic load predictions agree with the experimental results with an average error of 5 percent in the speed range 800 to 6000 rpm. Tooth separation (or bounce), which was observed in the experimental data for light torque, high speed conditions, was simulated by the computer model. The model was also successful in simulating the degree of load sharing between gear teeth in the multiple tooth contact region.

Rebbechi, Brian↗
NASA NTRS1981

The extension of the thermal-vacuum test optimization program to multiple flights

The thermal vacuum test optimization model developed to provide an approach to the optimization of a test program based on prediction of flight performance with a single flight option in mind is extended to consider reflight as in space shuttle missions. The concept of 'utility', developed under the name of 'availability', is used to follow performance through the various options encountered when the capabilities of reflight and retrievability of space shuttle are available. Also, a 'lost value' model is modified to produce a measure of the probability of a mission's success, achieving a desired utility using a minimal cost test strategy. The resulting matrix of probabilities and their associated costs provides a means for project management to evaluate various test and reflight strategies.

Williams, R. E.↗