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 271 records · Page 15

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

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

Owocki, Stanley P.↗

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

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

Puls, J.↗

Theorems on symmetrics and flux conservation in radiative transfer using the matrix operator theory

The matrix operator approach to radiative transfer is shown to be a very powerful technique in establishing symmetry relations for multiple scattering in inhomogeneous atmospheres. Symmetries are derived for the reflection and transmission operators using only the symmetry of the phase function. These results will mean large savings in computer time and storage for performing calculations for realistic planetary atmospheres using this method. The results have also been extended to establish a condition on the reflection matrix of a boundary in order to preserve reciprocity. Finally energy conservation is rigorously proven for conservative scattering in inhomogeneous atmospheres.

Kattawar, G. W.↗

Theorems on symmetries and flux conservation in radiative transfer using the matrix operator theory.

The matrix operator approach to radiative transfer is shown to be a very powerful technique in establishing symmetry relations for multiple scattering in inhomogeneous atmospheres. Symmetries are derived for the reflection and transmission operators using only the symmetry of the phase function. These results will mean large savings in computer time and storage for performing calculations for realistic planetary atmospheres using this method. The results have also been extended to establish a condition on the reflection matrix of a boundary in order to preserve reciprocity. Finally energy conservation is rigorously proven for conservative scattering in inhomogeneous atmospheres.

Kattawar, G. W.↗

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Interaction

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, N. V.↗

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Migration

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, Nelson V.↗

Investigating and implementing enhancements to the simulation of short-term collision hazards

A software tool, the Relative Collision Matrix (RCM), has been developed to provide a quick-look representation of the shortterm collision hazard to space systems from a fragmentation event in Earth orbit. The software performs multiple fragmentation simulations of space objects to quantify the probability of collision for a satellite or a constellation of satellites nearby. Previously, the results were displayed in a color matrix format which showed the relative hazard of each constellation. The RCM can be used for scientific research and operational assessments even though it was designed for test and evaluation applications. Because of its successful use as an analytical tool, the capabilities of RCM are being extended by enhancing the orbital hazard analysis routines, developing ballistic trajectory hazard analysis routines, and expanding the breakup modeling. Improvements are also being made to the RCM's usability and presentation quality by developing a graphical user interface and by providing graphical animated and nonanimated output.

Huth, Jeffrey S.↗

Unsteady airloading panel method for prop-fans

An unsteady, linear, compressible, three dimensional lifting surface panel method is presented for use with propellers and prop-fans. Both gust and blade vibration problems can be treated. Derivation of the theory is based on the acceleration potential method and can be considered a generalization of the analogous wing theory to include rotation and multiple blades. As with the wing methods, the integral equation is discretized and solved by matrix inversion. Sample calculations are use to explore the difference between 2D and 3D results and the significance of unsteadiness in angular inflow calculations for propellers. It is shown at the once-per-revolution frequency of propellers, that quasi-steady methods should be adequate. At the high frequencies of interest in noise analysis, 2D strip methods give adequate results away from blade tips. For flutter, frequencies lie in an intermediate range where a full 3D unsteady technique seems to be required.

Hanson, D. B.↗

Method for Fabricating Composite Structures Including Continuous Press Forming and Pultrusion Processing

A method for fabricating composite structures at a low-cost, moderate-to-high production rate is disclosed. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply, and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform, and debulking the partially debulked preform to form a netshape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length, and electron-beam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates, resulting in lower cost and high structural performance.

Farley, Gary L.↗

Method for Fabricating Composite Structures Using Continuous Press Forming

A method for fabricating composite structures at a low-cost. moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a net-shape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electron-beam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates. resulting in lower cost and high structural performance.

Farley, Gary L.↗

Method for Fabricating Composite Structures Using Pultrusion Processing

A method for fabricating composite structures at a low-cost, moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a net-shape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electron-beam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates, resulting in lower cost and high structural performance.

Farley, Gary L.↗

Method for Fabricating Composite Structures Using Pultrusion Processing

A method for fabricating composite structures at a low-cost, moderate-to-high production rate. A first embodiment of the method includes employing a continuous press forming fabrication process. A second embodiment of the method includes employing a pultrusion process for obtaining composite structures. The methods include coating yarns with matrix material, weaving the yarn into fabric to produce a continuous fabric supply and feeding multiple layers of net-shaped fabrics having optimally oriented fibers into a debulking tool to form an undebulked preform. The continuous press forming fabrication process includes partially debulking the preform, cutting the partially debulked preform and debulking the partially debulked preform to form a netshape. An electron-beam or similar technique then cures the structure. The pultrusion fabric process includes feeding the undebulked preform into a heated die and gradually debulking the undebulked preform. The undebulked preform in the heated die changes dimension until a desired cross-sectional dimension is achieved. This process further includes obtaining a net-shaped infiltrated uncured preform, cutting the uncured preform to a desired length and electronbeam curing (or similar technique) the uncured preform. These fabrication methods produce superior structures formed at higher production rates, resulting in lower cost and high structural performance.

Farley, Gary L.↗

Fabrication of Polyimide-Matrix/Carbon and Boron-Fiber Tape

The term HYCARB denotes a hybrid composite of polyimide matrices reinforced with carbon and boron fibers. HYCARB and an improved process for fabricating dry HYCARB tapes have been invented in a continuing effort to develop lightweight, strong composite materials for aerospace vehicles. Like other composite tapes in this line of development, HYCARB tapes are intended to be used to build up laminated structures having possibly complex shapes by means of automated tow placement (ATP) - a process in which a computer-controlled multiaxis machine lays down prepreg tape or tows. The special significance of the present process for making dry HYCARB for ATP is that it contributes to the reduction of the overall cost of manufacturing boron-reinforced composite-material structures while making it possible to realize increased compression strengths. The present process for making HYCARB tapes incorporates a "wet to dry" process developed previously at Langley Research Center. In the "wet to dry" process, a flattened bundle of carbon fiber tows, pulled along a continuous production line between pairs of rollers, is impregnated with a solution of a poly(amide acid) in N-methyl-2-pyrrolidinone (NMP), then most of the NMP is removed by evaporation in hot air. In the present case, the polyamide acid is, more specifically, that of LaRC. IAX (or equivalent) thermoplastic polyimide, and the fibers are, more specifically, Manganite IM7 (or equivalent) polyacrylonitrile- based carbon filaments that have a diameter of 5.2 m and are supplied in 12,000-filament tows. The present process stands in contrast to a prior process in which HYCARB tape was made by pressing boron fibers into the face of a wet carbon-fiber/ poly(amide acid) prepreg tape . that is, a prepreg tape from which the NMP solvent had not been removed. In the present process, one or more layer(s) of side-by-side boron fibers are pressed between dry prepreg tapes that have been prepared by the aforementioned gwet to dry h process. The multilayer tape is then heated to imidize the matrix material and remove most of the remaining solvent, and is pressed to consolidate the multiple layers into a dense tape. For tests, specimens of HYCARB tapes and laminated composite panels made from HYCARB tape were prepared as follows: HYCARB tapes were fabricated as described above. Each panel was made by laying down ten layers of tape, containing, variously, one, two, or three boron-fiber plies and the remainder carbon- fiber-only plies (see figure). Each panel was made by laying down ten layers of tape. Each panel was then cured by heating to a temperature of 225 C for 15 minutes, then pressing at 200 psi (A1.4 MPa) while heating to 371 C, holding at 371 C for 1 hour, then continuing to hold pressure during cooling. Control specimens that were otherwise identical except that they did not contain boron fibers also were prepared. In room-temperature flexural tests, the HYCARB specimens performed comparably to the control specimens; in room-temperature, open-hole compression tests, the HYCARB specimens performed slightly better, by amounts that increased with boron content.

Belvin, Harry L.↗

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

Efficient Preconditioning of a High-Order Solver for Multiple Physics

This work addresses preconditioning approaches for an implicit high-order solver frame-work applied to multiple physics. The solver is based on a space-time spectral element method and matrix-free Newton-Krylov solver developed at NASA over the recent years. Within this context, most preconditioning methods are impractical, as the computational time and memory requirements scale poorly with increasing polynomial orders. To improve computational efficiency, we first describe a novel entity-based Block Jacobi preconditioner for the continuous-Galerkin solution of the linear-elasticity and linear-shell equations. Second, we introduce a multigrid algorithm to further reduce time-to-solution on stiff cases arising from continuous-and discontinuous-Galerkin discretizations. Results obtained on relevant single-physics reference solutions, demonstrate the feasibility of the methods, paving the way for high-order solutions of fully coupled multi-physics problems.

STMD↗