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 487 records · Page 27

Apollo 14 breccias - General characteristics and classification.

The Apollo 14 breccias are complex and are characterized by a wide range of clast types and textures. The principal fragment types in the breccias are (1) hornfelsed noritic microbreccias; (2) micronorite hornfels; (3) devitrified glass; and (4) ophitic and other kinds of basalt. Based on fragment population, nature of the matrix, grain size and porosity, metamorphic history, and bulk chemical composition, the Apollo 14 breccias can be classified chiefly into regolith microbreccias, Fra Mauro breccias, and spherule-rich microbreccias. Further subdivisions are based principally on the interpretation of metamorphic history. Twenty-six Apollo 14 breccias have been so classified. The complex fragment types contained in Apollo 14 breccias represent multiple episodes of heating and fragmentation. They probably are pre-Imbrian and were not produced by the impact effects of a single event. The variability of their fragment populations is important in the interpretation of pre-Imbrian geologic history.

Chao, E. C. T.↗

Evaluating Crystallinity in Thermoplastic composites

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 manufacturing required in space. Most thermoplastic resins are semicrystalline polymers. The properties of semicrystalline thermoplastics are largely influenced by their crystallinity, which can vary due to many factors including thermal treatments, environmental conditions, and mechanical deformation. Monitoring the crystallinity of thermoplastic composites is key to ensuring these materials reliably meet the high demands required by space exploration. This talk discusses the use of multiple techniques such as Polarized Light Optical Microscopy and Fourier-Transform Infrared Spectroscopy to characterize the crystallinity in various thermoplastic composites, including carbon fiber reinforced PMCs and novel bio-based Martian and Lunar regolith composites designed for in-situ manufacturing. 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 space applications.

Thermoplastics↗

Flight-determined stability analysis of multiple-input-multiple-output control systems

Singular value analysis can give conservative stability margin results. Applying structure to the uncertainty can reduce this conservatism. This paper presents flight-determined stability margins for the X-29A lateral-directional, multiloop control system. These margins are compared with the predicted unscaled singular values and scaled structured singular values. The algorithm was further evaluated with flight data by changing the roll-rate-to-aileron-command-feedback gain by +/- 20 percent. Also presented are the minimum eigenvalues of the return difference matrix which bound the singular values. Extracting multiloop singular values from flight data and analyzing the feedback gain variations validates this technique as a measure of robustness. This analysis can be used for near-real-time flight monitoring and safety testing.

Burken, John J.↗

Layout optimization with algebraic multigrid methods

Finding the optimal position for the individual cells (also called functional modules) on the chip surface is an important and difficult step in the design of integrated circuits. This paper deals with the problem of relative placement, that is the minimization of a quadratic functional with a large, sparse, positive definite system matrix. The basic optimization problem must be augmented by constraints to inhibit solutions where cells overlap. Besides classical iterative methods, based on conjugate gradients (CG), we show that algebraic multigrid methods (AMG) provide an interesting alternative. For moderately sized examples with about 10000 cells, AMG is already competitive with CG and is expected to be superior for larger problems. Besides the classical 'multiplicative' AMG algorithm where the levels are visited sequentially, we propose an 'additive' variant of AMG where levels may be treated in parallel and that is suitable as a preconditioner in the CG algorithm.

Regler, Hans↗

Aging and Spaceflight: Catalase Targeted to Mitochondria Alters Skeletal Structure and Responses to Musculoskeletal Disuse

Microgravity and ionizing radiation in the spaceflight environment pose multiple challenges to homeostasis and may contribute to cellular stress. Effects may include increased generation of reactive oxygen species (ROS), DNA damage and repair error, cell cycle arrest, cell senescence or death. Our central hypothesis is that prolonged exposure to the spaceflight environment leads to excess production of ROS and oxidative damage, culminating in accelerated tissue degeneration which resembles aging. The main goal of this project is to determine the importance of cellular redox defense for physiological adaptations and tissue degeneration in the space environment. To accomplish this, we will use both wildtype (WT) mice and a well-established, genetically-engineered animal model (mCAT mice) which displays extended lifespan (Schriner et al. 2005). The animal model selected to test these ideas is engineered to quench ROS in mitochondria by targeted over-expression of the human catalase gene to the mitochondrial matrix. We showed previously that mCAT mice express the catalase transgene in skeletal tissues, bone forming osteoblasts, and bone resorbing osteoclasts. In addition, mCAT mice also display increased catalase activity in bone. Our findings revealed that exposure of adult, male, C57Bl/6J mice to simulated spaceflight (hindlimb unloading and gamma radiation) led to an increase in markers of oxidative damage (malondialdehyde, 4-hydroxynonenol) in skeletal tissue of WT mice but not mCAT mice. To extend our hypothesis to other, spaceflight-relevant tissues, we are performing a ground-based study simulating 30 days of spaceflight by hindlimb unloading to determine potential protective effects of mitochondrial catalase activity on aging of multiple tissues (cardiovascular, nervous and skeletal).

transgenic mice↗

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗

Numerical Simulations of Single and Multiple Scattering by Fractal Ice Clusters

We consider the scattering model in the form of a vertically and horizontally homogeneous particulate slab of an arbitrary optical thickness composed of widely separated fractal aggregates built of small spherical ice monomers. The aggregates are generated by applying three different approaches, including simulated cluster-cluster aggregation (CCA) and diffusion-limited aggregation (DLA) procedures. Having in mind radar remote-sensing applications, we report and analyze the results of computations of the backscattering circular polarization ratio obtained using efficient superposition T-matrix and vector radiative-transfer codes. The computations have been performed at a wavelength of 12.6 cm for fractal aggregates with the following characteristics: monomer refractive index m=1.78+i0.003, monomer radius r=1 cm, monomer packing density p=0.2, overall aggregate radii R in the range 4<=R<=10 cm and fractal dimensions D(sub f) 2.5 and 3. We show that for aggregates generated with simulated CCA and DLA procedures, the respective values of the backscattering circular polarization ratio differ weakly for D(sub f) 2.5, but the differences can increase somewhat for D(sub f)3, especially in case of an optically semi-infinite medium. For aggregates with a spheroidal overall shape, the dependence of the circular polarization ratio on the cluster morphology can be quite significant and increases with increasing the aspect ratio of the circumscribing spheroid.

Dlugach, Janna M.↗

Boundary layer simulator improvement

Boundary Layer Integral Matrix Procedure (BLIMPJ) has been identified by the propulsion community as the rigorous boundary layer program in connection with the existing JANNAF reference programs. The improvements made to BLIMPJ and described herein have potential applications in the design of the future Orbit Transfer Vehicle engines. The turbulence model is validated to include the effects of wall roughness and a way is devised to treat multiple smooth-rough surfaces. A prediction of relaminarization regions is examined as is the combined effects of wall cooling and surface roughness on relaminarization. A turbulence model to represent the effects of constant condensed phase loading is given. A procedure is described for thrust decrement calculation in thick boundary layers by coupling the T-D Kinetics Program and BLIMPJ and a way is provided for thrust loss optimization. Potential experimental studies in rocket nozzles are identified along with the required instrumentation to provide accurate measurements in support of the presented new analytical models.

Praharaj, Sarat C.↗

Coupled 2-dimensional cascade theory for noise and unsteady aerodynamics of blade row interaction in turbofans. Volume 1: Theory development and parametric studies

Typical analytical models for interaction between rotor and stator in a turbofan analyze the effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a separate calculation. In those models, there is a one-to-one relationship between wake harmonics and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the BPF noise harmonic, etc. This report presents a more complete model in which flow tangency boundary conditions are satisfied on two cascades in relative motion for several harmonics simultaneously. By an extension of S.N. Smith's code for two dimensional flat plate cascades, the noise generation/frequency scattering/blade row reflection problem is solved in a single matrix inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy in the higher BPF harmonics due to relative motion between the blade rows. Furthermore, when swirl between the rotor and stator is modeled, a 'mode trapping' effect occurs which explains observations on fans operating at rotational speeds below BFP cuton: the BPF mode amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. However, energy scattered into higher harmonics does propagate and dominates the spectrum at two and three times BPF. This report presents the complete derivation of the theory, comparison with a previous (more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either upstream or downstream of the stator and to counter rotation stages. The theory has been coded in a FORTRAN program called CUP2D, documented in Volume 2 of this report. It is concluded that the new features of this analysis - unsteady coupling, frequency scattering, and flow turning between rotor and stator - have a profound effect on noise generation caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate for noise analysis and prediction.

Hanson, Donald B.↗

Processing and Characterization of Carbon Nanotube Composites

Recent advances in the synthesis of large-scale quantities of carbon nanotubes (CNT) have provided the opportunity to study the mechanical properties of polymer matrix composites using these novel materials as reinforcement. Nanocomp Technologies, Inc. currently supplies large sheets with dimensions up to 122 cm x 244 cm containing both single-wall and few-wall CNTs. The tubes are approximately 1 mm in length with diameters ranging from 8 to 12 nm. In the present study being conducted at NASA Langley Research Center (LaRC), single and multiple layers of CNT sheets were infused or coated with various polymer solutions that included commercial toughened-epoxies and bismaleimides, as well as a LaRC developed polyimide. The resulting CNT composites were tested in tension using a modified version of ASTM D882-12 to determine their strength and modulus values. The effects of solvent treatment and mechanical elongation/alignment of the CNT sheets on the tensile performance of the composite were determined. Thin composites (around 50 wt% CNT) fabricated from acetone condensed and elongated CNT sheets with either a BMI or polyimide resin solution exhibited specific tensile moduli approaching that of toughened epoxy/ IM7 carbon fiber unidirectional composites.

Can, Roberto J.↗

Development of advanced fuel cell system, phase 3

A multiple task research and development program was performed to improve the weight, life, and performance characteristics of hydrogen-oxygen alkaline fuel cells for advanced power systems. Gradual wetting of the anode structure and subsequent long-term performance loss was determined to be caused by deposition of a silicon-containing material on the anode. This deposit was attributed to degradation of the asbestos matrix, and attention was therefore placed on development of a substitute matrix of potassium titanate. An 80 percent gold 20 percent platinum catalyst cathode was developed which has the same performance and stability as the standard 90 percent gold - 10 percent platinum cathode but at half the loading. A hybrid polysulfone/epoxy-glass fiber frame was developed which combines the resistance to the cell environment of pure polysulfone with the fabricating ease of epoxy-glass fiber laminate. These cell components were evaluated in various configurations of full-size cells. The ways in which the baseline engineering model system would be modified to accommodate the requirements of the space tug application are identified.

Handley, L. M.↗

Systems Maintenance Automated Repair Tasks (SMART)

SMART is a uniform automated discrepancy analysis and repair-authoring platform that improves technical accuracy and timely delivery of repair procedures for a given discrepancy (see figure a). SMART will minimize data errors, create uniform repair processes, and enhance the existing knowledge base of engineering repair processes. This innovation is the first tool developed that links the hardware specification requirements with the actual repair methods, sequences, and required equipment. SMART is flexibly designed to be useable by multiple engineering groups requiring decision analysis, and by any work authorization and disposition platform (see figure b). The organizational logic creates the link between specification requirements of the hardware, and specific procedures required to repair discrepancies. The first segment in the SMART process uses a decision analysis tree to define all the permutations between component/ subcomponent/discrepancy/repair on the hardware. The second segment uses a repair matrix to define what the steps and sequences are for any repair defined in the decision tree. This segment also allows for the selection of specific steps from multivariable steps. SMART will also be able to interface with outside databases and to store information from them to be inserted into the repair-procedure document. Some of the steps will be identified as optional, and would only be used based on the location and the current configuration of the hardware. The output from this analysis would be sent to a work authoring system in the form of a predefined sequence of steps containing required actions, tools, parts, materials, certifications, and specific requirements controlling quality, functional requirements, and limitations.

Schuh, Joseph↗

Piloted Simulation Assessment of the Impact of Flexible Structures on Handling Qualities of Generic Supersonic Aircraft

The NASA Langley Research Center Cockpit Motion Facility (CMF) was used to conduct a piloted simulation assessment of the impact of flexible structures on flying qualities. The CMF was used because of its relatively high bandwidth, six degree-of-freedom motion capability. Previous studies assessed and attempted to mitigate the effects of multiple dynamic aeroservoelastic modes (DASE). Those results indicated problems existed, but the specific cause and effect was difficult to ascertain. The goal of this study was to identify specific DASE frequencies, damping ratios, and gains that cause degradation in handling qualities. A generic aircraft simulation was developed and designed to have Cooper-Harper Level 1 handling qualities when flown without DASE models. A test matrix of thirty-six DASE modes was implemented. The modes had frequencies ranging from 1 to 3.5 Hz and were applied to each axis independently. Each mode consisted of a single axis, frequency, damping, and gain, and was evaluated individually by six subject pilots with test pilot backgrounds. Analysis completed to date suggests that a number of the DASE models evaluated degrade the handling qualities of this class of aircraft to an uncontrollable condition.

Stringer, Mary T.↗

Basic failure mechanisms in advanced composites

Fundamental failure mechanisms in carbon-epoxy composites were studied for more reliable prediction of the performance of these materials. Single and multiple fiber specimens were tested under tensile loads, and the sequence of failure events was observed. Parameters such as resin crack sensitivity, fiber surface treatment and variations in fibers from batch to batch are being evaluated. The analysis of bulk composite fracture processes using acoustic emission techniques is being studied in order to correlate microscopic observations with bulk composite behavior. Control of the fracture process through matrix and interface modification is being attempted, and study of failure processes in composite/metal specimens is being conducted. Most of the studies involved DEN 438 epoxy novolac as the matrix, but some experiments are now underway using the higher temperature resin ERLA 4617.

Mullin, J. V.↗

Accelerated Adaptive MGS Phase Retrieval

The Modified Gerchberg-Saxton (MGS) algorithm is an image-based wavefront-sensing method that can turn any science instrument focal plane into a wavefront sensor. MGS characterizes optical systems by estimating the wavefront errors in the exit pupil using only intensity images of a star or other point source of light. This innovative implementation of MGS significantly accelerates the MGS phase retrieval algorithm by using stream-processing hardware on conventional graphics cards. Stream processing is a relatively new, yet powerful, paradigm to allow parallel processing of certain applications that apply single instructions to multiple data (SIMD). These stream processors are designed specifically to support large-scale parallel computing on a single graphics chip. Computationally intensive algorithms, such as the Fast Fourier Transform (FFT), are particularly well suited for this computing environment. This high-speed version of MGS exploits commercially available hardware to accomplish the same objective in a fraction of the original time. The exploit involves performing matrix calculations in nVidia graphic cards. The graphical processor unit (GPU) is hardware that is specialized for computationally intensive, highly parallel computation. From the software perspective, a parallel programming model is used, called CUDA, to transparently scale multicore parallelism in hardware. This technology gives computationally intensive applications access to the processing power of the nVidia GPUs through a C/C++ programming interface. The AAMGS (Accelerated Adaptive MGS) software takes advantage of these advanced technologies, to accelerate the optical phase error characterization. With a single PC that contains four nVidia GTX-280 graphic cards, the new implementation can process four images simultaneously to produce a JWST (James Webb Space Telescope) wavefront measurement 60 times faster than the previous code.

Lam, Raymond K.↗

Particle Scattering in the Resonance Regime: Full-Wave Solution for Axisymmetric Particles with Large Aspect Ratios

Reliable descriptions of the optical properties of clouds and aerosols are essential for studies of radiative transfer in planetary atmospheres. The scattering algorithms provide accurate estimates of these properties for spherical particles with a wide range of sizes and refractive indices, but these methods are not valid for non-spherical particles (e.g., ice crystals, mineral dust, and smoke). Even though a host of methods exist for deriving the optical properties of nonspherical particles that are very small or very large compared with the wavelength, only a few methods are valid in the resonance regime, where the particle dimensions are comparable with the wavelength. Most such methods are not ideal for particles with sharp edges or large axial ratios. We explore the utility of an integral equation approach for deriving the single-scattering optical properties of axisymmetric particles with large axial ratios. The accuracy of this technique is shown for spheres of increasing size parameters and an ensemble of randomly oriented prolate spheroids of size parameter equal to 10.079368. In this last case our results are compared with published results obtained with the T-matrix approach. Next we derive cross sections, single-scattering albedos, and phase functions for cylinders, disks, and spheroids of ice with dimensions extending from the Rayleigh to the geometric optics regime. Compared with those for a standard surface integral equation method, the storage requirement and the computer time needed by this method are reduced, thus making it attractive for generating databases to be used in multiple-scattering calculations. Our results show that water ice disks and cylinders are more strongly absorbing than equivalent volume spheres at most infrared wavelengths. The geometry of these particles also affects the angular dependence of the scattering. Disks and columns with maximum linear dimensions larger than the wavelength scatter much more radiation in the forward and backward directions and much less radiation at intermediate phase angles than equivalent volume spheres.

Zuffada, Cinzia↗

Nonlinear Attitude Filtering Methods

The extended Kalman filter (EKF) is the workhorse of real-time spacecraft attitude estimation. Since the group SO3 of rotation matrices has dimension three, most attitude determination EKFs use lower- dimensional attitude parameterizations than the nine-parameter attitude matrix itself. The fact that all three- parameter representations of SO3 are singular or discontinuous for certain attitudes has led to extended discussions of constraints and attitude representations in EKFs. The most successful EKF uses a nonsingular parameterization for the global attitude, which necessarily has more than three parameters, while employing a three-component representation for the attitude errors. This filter has become known as the Multiplicative Extended Kalman Filter. These issues are now well understood, however, and the EKF has performed admirably in the vast majority of attitude determination applications. Nevertheless, poor performance or even divergence arising from the linearization implicit in the EKF has led to the development of nonlinear filters, most recently sigma point or unscented filters and particle filters.

Markley, F. Landis↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar↗