Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the three-dimensional unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalizing procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the three-dimensional Euler equations.

Yokota, J. W.↗
NASA NTRS1989

Cloud field classification based on textural features

An essential component in global climate research is accurate cloud cover and type determination. Of the two approaches to texture-based classification (statistical and textural), only the former is effective in the classification of natural scenes such as land, ocean, and atmosphere. In the statistical approach that was adopted, parameters characterizing the stochastic properties of the spatial distribution of grey levels in an image are estimated and then used as features for cloud classification. Two types of textural measures were used. One is based on the distribution of the grey level difference vector (GLDV), and the other on a set of textural features derived from the MaxMin cooccurrence matrix (MMCM). The GLDV method looks at the difference D of grey levels at pixels separated by a horizontal distance d and computes several statistics based on this distribution. These are then used as features in subsequent classification. The MaxMin tectural features on the other hand are based on the MMCM, a matrix whose (I,J)th entry give the relative frequency of occurrences of the grey level pair (I,J) that are consecutive and thresholded local extremes separated by a given pixel distance d. Textural measures are then computed based on this matrix in much the same manner as is done in texture computation using the grey level cooccurrence matrix. The database consists of 37 cloud field scenes from LANDSAT imagery using a near IR visible channel. The classification algorithm used is the well known Stepwise Discriminant Analysis. The overall accuracy was estimated by the percentage or correct classifications in each case. It turns out that both types of classifiers, at their best combination of features, and at any given spatial resolution give approximately the same classification accuracy. A neural network based classifier with a feed forward architecture and a back propagation training algorithm is used to increase the classification accuracy, using these two classes of features. Preliminary results based on the GLDV textural features alone look promising.

Sengupta, Sailes Kumar↗
NASA NTRS1991

Partitioning of unstructured problems for parallel processing

Many large-scale computational problems are based on unstructured computational domains. Primary examples are unstructured grid calculations based on finite volume methods in computational fluid dynamics, or structural analysis problems based on finite element approximations. The question of how to distribute such unstructured computational domains over a large number of processors in a MIMD machine with distributed memory is addressed. A graph theoretical framework for these problems is established. Based on this framework three decomposition algorithms are introduced. In particular a new decomposition algorithm is discussed, which is based on the computation of an eigenvector of the Laplacian matrix associated with the graph. Numerical comparisons on large-scale two- and three-dimensional problems demonstrate the superiority of the new spectral bisection algorithm.

Simon, H. D.↗
NASA NTRS1992

Environmental effects on long term behavior of composite laminates

Model equations are presented for approximate methods simulating the long-term behavior of composite materials and structures in hot/humid service environments. These equations allow laminate property upgradings with time, and can account for the effects of service environments on creep response. These methodologies are illustrated for various individual and coupled temperature/moisture, longitudinal/transverse, and composite material type cases. Creep deformation is noted to rise dramatically for cases of matrix-borne, but not of fiber-borne, loading in hot, humid environments; the coupled influence of temperature and moisture is greater than a mere combination of their individual influences.

Singhal, S. N.↗
DOE OSTI2024

Development of creep resistant ferritic-martensitic steels for operation at 650°C

Creep resistant ferritic-martensitic steels capable of operating at 650°C have long been under development for electricity generation. This presentation will showcase the design and development of alloys CPJ-7 and JMP with approximately 420,000 h of in-house creep testing, the equivalent of ~50 years of cumulative creep tests. Approximately fifty different alloy compositions were manufactured as part of an ICME approach. CPJ-7 present superior creep properties compared to COST alloys for steam turbine and P91/92 for boiler applications. The prolonged creep life was attributed to slowing down the destabilization of fine precipitates at 650°C. The JMP steels were designed with higher Co for solid solution strengthening and increased W for matrix strength and stability. The creep life compared to CPJ-7 increased between 118 and 150% at 650°C depending on the stress. On a Larson-Miller plot, the performance of the JMP steels surpasses that of MARBN steel. Composition-microstructure-properties relationships will be discussed.

Detrois, Martin↗
NASA NTRS2012

Processing and Damage Tolerance of Continuous Carbon Fiber Composites Containing Puncture Self-Healing Thermoplastic Matrix

Research at NASA Langley Research Center (NASA LaRC) has identified several commercially available thermoplastic polymers that self-heal after ballistic impact and through-penetration. One of these resins, polybutadiene graft copolymer (PB(sub g)), was processed with unsized IM7 carbon fibers to fabricate reinforced composite material for further evaluation. Temperature dependent characteristics, such as the degradation point, glass transition (T(sub g)), and viscosity of the PBg polymer were characterized by thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and dynamic parallel plate rheology. The PBg resin was processed into approximately equal to 22.0 cm wide unidirectional prepreg tape in the NASA LaRC Advanced Composites Processing Research Laboratory. Data from polymer thermal characterization guided the determination of a processing cycle used to fabricate quasi-isotropic 32-ply laminate panels in various dimensions up to 30.5cm x 30.5cm in a vacuum press. The consolidation quality of these panels was analyzed by optical microscopy and acid digestion. The process cycle was further optimized based on these results and quasi-isotropic, [45/0/-45/90]4S, 15.24cm x 15.24cm laminate panels were fabricated for mechanical property characterization. The compression strength after impact (CAI) of the IM7/pBG composites was measured both before and after an elevated temperature and pressure healing cycle. The results of the processing development effort of this composite material as well as the results of the mechanical property characterization are presented in this paper.

Grimsley, Brian W.↗
DOE OSTI2020

Thermal shock resistance of multilayer silicon carbide receiver tubes for 800°C molten salt concentrating solar power application

CSP power tower receiver systems during rapid transient weather periods can be vulnerable to thermal shock conditions from rain that which can facilitate the onset of leaks and failures that can have catastrophic consequences. Silicon carbide (SiC) materials have attractive receiver application characteristics for being light weight, having high-strength and excellent thermal shock resistance performance which make them a particularly good fit for receiver absorber materials in CSP. In this investigation, the performance characteristics of Ceramic Tubular Products (CTP) SiC ceramic matrix composite (CMC), multilayered tubes were explored with respect to thermal shock performance for solar receiver applications in next generation CSP plants. Here, thermal shock testing was performed at the Sandia National Laboratories (SNL) Solar Furnace facility using a dynamic stage and thermal shock tube test setup. The tubes tested under incident solar heat flux of 100 W/cm2 were heated with inner tube temperatures reaching approximately 800 °C, with outer temperatures exceeding or just reaching 1000 El for the multilayer and monolithic SiC tubes respectively. The tubes were then quenched with simulated rain. The tubes were then cooled and subjected to hoop stress analysis using an Instron device to assess their subsequent mechanical strength. Furthermore, the on-sun study experimental results indicate an average of 24.2% and 97% higher hoop strength for the CMC tubes than those composed of monolithic SiC and aluminum oxide (A1203) respectively.

14 SOLAR ENERGY↗
NASA NTRS1983

EGRET - The high energy gamma ray telescope for NASA's Gamma Ray Observatory

The EGRET high energy gamma-ray telescope under development for NASA's Gamma Ray Observatory will have an energy range of approximately 12 to 30,000 MeV, energy resolution of about 15 percent FWHM over most of that range, an effective area of about 2000 sq cm at high energies, and single photon angular accuracy of approximately 2 deg at 100 MeV, less than 0.1 deg above 5 GeV. This instrument can locate strong sources to an accuracy of about 5 arc min. The instrument utilizes a set of digital spark chambers interleaved with tantalum foils for detection and identification of gamma-ray events, and a large NaI(Tl) scintillator for energy determination. The system is triggered by a coincidence matrix using two arrays of plastic scintillation counters and a large plastic scintillator anticoincidence dome that rejects incident charged particles.

Fichtel, C. E.↗
NASA NTRS1989

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗
NASA NTRS1991

Deformation and failure mechanisms in metal matrix composites

An investigation was undertaken to determine the key deformation mechanisms and their interaction leading to failure of both 0 degree and 90 degree Ti 15-3/SCS-6 laminae under monotonic loading. The experimental results suggest that inelastic deformation in the 0-degree lamina is dominated by plastic deformation and that in the 90-degree lamina is dominated by both fiber-matrix debonding and plasticity. The loading-unloading response, monitoring of Poisson's ratio and microscopy were utilized to identify the key deformation mechanisms. The sequence of deformation mechanisms leading to failure are identified for both the 0 and the 90-degree specimens. The threshold strains for plasticity or damage which are referred to as 'microdeformation' in the 0 deg and 90 deg laminae are approximately 0.004 and 0.002, respectively, at room temperature. These strain levels may be considered critical in initiation based structural design with these composites.

Newaz, G.↗
NASA NTRS2004

Inducer Hydrodynamic Forces in a Cavitating Environment

Marshall Space Flight Center has developed and demonstrated a measurement device for sensing and resolving the hydrodynamic loads on fluid machinery. The device - a derivative of the six-component wind tunnel balance - senses the forces and moments on the rotating device through a weakened shaft section instrumented with a series of strain gauges. This rotating balance was designed to directly measure the steady and unsteady hydrodynamic loads on an inducer, thereby defining the amplitude and frequency content associated with operating in various cavitation modes. The rotating balance was calibrated statically using a dead-weight load system in order to generate the 6 x 12 calibration matrix later used to convert measured voltages to engineering units. Structural modeling suggested that the rotating assembly first bending mode would be significantly reduced with the balance s inclusion. This reduction in structural stiffness was later confirmed experimentally with a hammer-impact test. This effect, coupled with the relatively large damping associated with the rotating balance waterproofing material, limited the device s bandwidth to approximately 50 Hertz Other pre-test validations included sensing the test article rotating assembly built-in imbalance for two configurations and directly measuring the assembly mass and buoyancy while submerged under water. Both tests matched predictions and confirmed the device s sensitivity while stationary and rotating. The rotating balance was then demonstrated in a water test of a full-scale Space Shuttle Main Engine high-pressure liquid oxygen pump inducer. Experimental data was collected a scaled operating conditions at three flow coefficients across a range of cavitation numbers for the single inducer geometry and radial clearance. Two distinct cavitation modes were observed symmetric tip vortex cavitation and alternate-blade cavitation. Although previous experimental tests on the same inducer demonstrated two additional cavitation modes at lower inlet pressures, these conditions proved unreachable with the rotating balance installed due to the intense dynamic environment. The sensed radial load was less influenced by flow coefficient than by cavitation number or cavitation mode although the flow coefficient range was relatively narrow. Transition from symmetric tip vortex to alternate-blade cavitation corresponded to changes in both radial load magnitude and radial load orientation relative to the inducer. Sensed moments indicated that the effective load center moved downstream during this change in cavitation mode. An occurrence of "higher+rdex cavitation" was also detected in both the stationary pressures and the rotating balance data although the frequency of the phenomena was well above the reliable bandwidth of the rotating balance. In summary the experimental tests proved both the concept and device s capability despite the limitations and confirmed that hydrodynamically-induced forces and moments develop in response to the unbalanced pressure field, which is, in turn, a product of the cavitation environment.

Skelley, Stephen E.↗
DOE OSTI2020

Stochastically Realized Observables for Excitonic Molecular Aggregates

In this work, we show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}$($N$ 3 ) for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}$($N$ log($N$)), enabling quick study of systems with a million of coupled molecules on the micrometer size scale. We calculate several important experimental observables, including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Gravitational wave modes in matter

Here, a general linear gauge-invariant equation for dispersive gravitational waves (GWs) propagating in matter is derived. This equation describes, on the same footing, both the usual tensor modes and the gravitational modes strongly coupled with matter. It is shown that the effect of matter on the former is comparable to diffraction and therefore negligible within the geometrical-optics approximation. However, this approximation is applicable to modes strongly coupled with matter due to their large refractive index. GWs in ideal gas are studied using the kinetic average-Lagrangian approach and the gravitational polarizability of matter that we have introduced earlier. In particular, we show that this formulation subsumes the kinetic Jeans instability as a collective GW mode with a peculiar polarization, which is derived from the dispersion matrix rather than assumed a priori. This forms a foundation for systematically extending GW theory to GW interactions with plasmas, where symmetry considerations alone are insufficient to predict the wave polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
NASA NTRS1974

Analysis of twisting stiffness for a multifiber composite layer

The twisting stiffness of a rectangular cross section consisting of a single row of solid circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Dirichlet torsion problem of a multielement region and solved by the boundary-point least-squares method. Numerical results for a single-fiber square cross section compare favorably with previous relaxation-method results. New numerical results for three and five-fiber composites suggest that the torsional rigidity of a multifiber composite can be approximated from the torsional rigidities of single and three-fiber models.

Bert, C. W.↗
NASA NTRS1979

Completion of the Viking Labeled Release experiment on Mars

The final Labeled Release (LR) cycle on each Viking lander tested a surface sample that had been stored for several months at approximately 10 C prior to the onset of the active sequence. At each lander site, activity was strongly diminished. This thermal sensitivity of the active agent on the surface of Mars is consistent with a biological explanation of the LR experiment. At the end of one of these cycles, the incubation mixture was heated to 50 C to release any radioactive gas trapped in the sample matrix. The results suggest that more than one carbon substrate is involved in the LR reaction on Mars. The thermal data from the stored samples, coupled with data from previous cycles, have formed the basis for evaluation of the thermal decomposition of the Mars active agent. The slope of the resulting Arrhenius plot has been used to test the fit of other flight data and to calculate the activation energy for thermal decomposition of the Mars agent. The results and their interpretation still leave unresolved the question of whether the Mars LR data were generated by biological or chemical activity.

Levin, G. V.↗
NASA NTRS1990

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗
NASA NTRS1991

Mechanical properties of monolithic and particulate composites of L12 forms of Al3Ti

Preliminary studies of simple compressive properties of particulate composites with L1(2)-modified Al3Ti's as the matrices were performed as functions of matrix composition, amount of TiB2, processing technique, temperature, and strain rate. At and below about 900 K the introduction of TiB2 particles overwhelms minor chemical effects, and extremely high-yield strengths (about 2000 MPa at the 20 vol pct TiB2 level) were obtained. As L1(2)-modified Al3Ti has about half the density of the Ni-base superalloys, such values translate to strength/density ratios approximately four times those for superalloys. This advantage is not maintained by the Al3Ti-based composites above about 1000 K, particularly under slow strain rate plastic flow conditions, and it is proposed that this behavior could be due to the dissolution of intermetallic second phases as well as grain boundary creep mechanisms.

Whittenberger, J. D.↗