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 793 records · Page 44

Distributed and communication-efficient solutions to linear equations with special sparse structure

In this paper we report two distributed and communication-efficient algorithms based on the multi-agent system are proposed to solve a system of linear equations with the Laplacian sparse system matrix. One algorithm is based on the gradient descent method in optimization. In this algorithm, the agents only share partial information instead of all of their collective state vectors to save significant communication. The other algorithm is obtained by approximating Newton’s method for a faster convergence rate. Although it requires twice as much communication as the first one, it is still communication-efficient given the low dimension of the information shared among agents. The convergence at a linear rate is proved for both algorithms, and a comprehensive comparison of their convergence rate, communication burden, and computation costs is also performed. The proposed algorithms can be applied to various systems to solve those problems that can be modeled as a system of linear equations with a Laplacian sparse system matrix. Simulation results with the electric power system illustrate their effectiveness.

42 ENGINEERING↗

An Uncertainty Structure Matrix for Models and Simulations

Software that is used for aerospace flight control and to display information to pilots and crew is expected to be correct and credible at all times. This type of software is typically developed under strict management processes, which are intended to reduce defects in the software product. However, modeling and simulation (M&S) software may exhibit varying degrees of correctness and credibility, depending on a large and complex set of factors. These factors include its intended use, the known physics and numerical approximations within the M&S, and the referent data set against which the M&S correctness is compared. The correctness and credibility of an M&S effort is closely correlated to the uncertainty management (UM) practices that are applied to the M&S effort. This paper describes an uncertainty structure matrix for M&S, which provides a set of objective descriptions for the possible states of UM practices within a given M&S effort. The columns in the uncertainty structure matrix contain UM elements or practices that are common across most M&S efforts, and the rows describe the potential levels of achievement in each of the elements. A practitioner can quickly look at the matrix to determine where an M&S effort falls based on a common set of UM practices that are described in absolute terms that can be applied to virtually any M&S effort. The matrix can also be used to plan those steps and resources that would be needed to improve the UM practices for a given M&S effort.

Green, Lawrence L.↗

Carbon-rich Chondritic Clast PV1 from the Plainview H-chondrite Regolith Formation from H3 Chondrite Material by Possible Cometary Impact

Chondritic clast PV1 from the Plainview H-chondrite regolith breccia is a subrounded, 5-mm diameter unequilibrated chondritic fragment that contains 13 wt% C occurring mainly within irregularly shaped 30-400-micron-size opaque patches. The clast formed from H3 chondrite material as indicated by the mean apparent chondrule diameter (310 micron vs. approximately 300 micron in H3 chondrites), the mean Mg-normalized refractory lithophile abundance ratio (1.00 +/- 0.09 XH), the previously determined 0-isotopic composition (Delta O-17 = 0.66% vs. 0.68 +/- 0.04%0 in H3 chondrites and 0.73 +/- 0.09% in H4-6 chondrites), the heterogeneous olivine compositions in grain cores (with a minimum range of Fal-19), and the presence of glass in some chondrules. Although the clast lacks the fine-grained, ferroan silicate matrix material present in type 3 ordinary chondrites, PV1 contains objects that appear to be recrystallized clumps of matrix material. Similarly, the apparent dearth of radial pyroxene and cryptocrystalline chondrules in PV1 is accounted for by the presence of some recrystallized fragments of these chondrule textural types. All of the chondrules in PV1 are interfused indicating that temperatures must have briefly reached approximately 1100C (the approximate solidus temperature of H-chondrite silicate). The most likely source of this heating was by an impact. Some metal was lost during impact heating as indicated by the moderately low abundance of metallic Fe-Ni in PV1 (approximately 14 wt%) compared to that in mean H chondrites (approximately 18 wt%). The carbon enrichment of the clast may have resulted from a second impact event, one involving a cometary projectile, possibly a Jupiter-family comet. As the clast cooled, it experienced hydrothermal alteration at low water/rock ratios as evidenced by the thick rims of ferroan olivine around low-FeO olivine cores. The C-rich chondritic clast was later incorporated into the H-chondrite parent-body regolith and extensively fractured and faulted.

Rubin, Alan E.↗

Breccia Formation at a Complex Impact Crater: Slate Islands, Lake Superior, Ontario, Canada

The Slate Islands impact structure is the eroded remnant of a approximately 30-32 km-diameter complex impact structure located in northern Lake Superior, Ontario, Canada. Target rocks are Archean supracrustal and igneous rocks and Proterozoic metavolcanics, metasediments, and diabase. A wide variety of breccias occurs on the islands, many of which contain fragments exhibiting shock metamorphic features. Aphanitic, narrow and inclusion-poor pseudotachylite veins, commonly with more or less parallel boundaries and apophyses branching off them, represent the earliest breccias formed during the compression stage of the impact process. Coarse-grained, polymictic elastic matrix breccias form small to very large, inclusion-rich dikes and irregularly shaped bodies that may contain altered glass fragments. These breccias have sharp contacts with their host rocks and include a wide range of fragment types some of which were transported over minimum distances of approximately 2 km away from the center of the structure. They cut across pseudotachylite veins and contain inclusions of them. Field and petrographic evidence indicate that these polymictic breccias formed predominantly during the excavation and central uplift stages of the impact process. Monomictic breccias, characterized by angular fragments and transitional contacts with their host rocks, occur in parautochthonous target rocks, mainly on the outlying islands of the Slate Islands archipelago. A few contain fragmented and disrupted, coarse-grained, polymictic clastic matrix breccia dikes. This is an indication that at least some of these monomictic breccias formed late in the impact process and that they are probably related to a late crater modification stage. A small number of relatively large occurrences of glass-poor, suevitic breccias occur at the flanks of the central uplift and along the inner flank of the outer ring of the Slate Islands complex crater. A coarse, glass-free, allogenic breccia, containing shatter-coned fragments derived from Proterozoic target rocks (upper target strata), observed at two locations may be analogous to the 'Bunt Breccia' of the Ries crater in Germany. At one of these locations this breccia lies close to a crater suevite deposit. At the other, it overlies parautochthonous, monomictic breccia. The State Islands impact breccias are superbly exposed, much better than breccias in most other terrestrial impact structures. Observations, including those indicative of multiple and and sequential processes, provide insight on how impact breccias form and how they relate to the various phases of the impact process. Eventually they will lead to an improved understanding of planetary impact processes.

Dressler, B. O.↗

A preliminary characterization of the tensile and fatigue behavior of tungsten-fiber/Waspaloy-matrix composite

A microstructural study and a preliminary characterization of the room temperature tensile and fatigue behavior of a continuous, tungsten fiber, Waspaloy-matrix composite was conducted. A heat treatment was chosen that would allow visibility of planar slip if it occurred during deformation, but would not allow growth of the reaction zone. Tensile and fatigue tests showed that the failed specimens contained transverse cracks in the fibers. The cracks that occurred in the tensile specimen were observed at the fracture surface and up to approximately 4.0 mm below the fracture surface. The crack spacing remained constant along the entire length of the cracked fibers. Conversely, the cracks that occurred in the fatigue specimen were only observed in the vicinity of the fracture surface. In instances where two fiber cracks occurred in the same plane, the matrix often necked between the two cracked fibers. Large groups of slip bands were generated in the matrix near the fiber cracks. Slip bands in the matrix of the tensile specimen were also observed in areas where there were no fiber cracks, at distances greater than 4 mm from the fracture surface. This suggests that the matrix plastically flows before fiber cracking occurs.

Corner, Ralph E.↗

The performance of additively manufactured Haynes 282 in supercritical CO 2

The use of supercritical carbon dioxide (sCO 2 ) as a working fluid is garnering interest in next generation power production systems due to the possibility of increased operational efficiencies and lower associated costs. Its implementation requires alloys with excellent high temperature strength and corrosion resistance, which makes Haynes 282® (H282) a suitable candidate. With increasing adoption of additive manufacturing (AM) in the energy industry, there is a need to investigate long-term sCO2 exposure on AM produced materials. Additively manufactured H282 (AM-H282) samples built using laser powder bed fusion in two different orientations were analyzed through imaging, gravimetric analysis, and room-temperature tensile testing before and after 1000-h exposure to CO 2 at 750 °C and 20 MPa. Performed imaging included observation of the sample surfaces and of the material bulk (cross sections) through means of Scanning Electron Microscopy (SEM) and optical microscopy. In comparison to wrought H282, it was found that the exposed material exhibited a similar Cr 2 O 3 oxide protective layer about 2 μm in thickness with additional top-surface TiO 2 oxide and carbon-rich precipitates. It was also similarly observed that internal oxidation was present but limited to a depth of 20 μm, and a 1–3 μm γ’ denuded region appeared to surround all internal and surface oxidation. Thermal aging effects were noted with the precipitation of needle-like structures in the γ/γ’ matrix and a coarsening of the γ’ precipitates from 28 nm to 73 nm. Mass measurements analyzing oxide precipitation revealed a larger increase compared to wrought, representing an approximate 20% difference. Tensile testing results showed similar behavior to wrought with a slight increase in yield strength and a large reduction in elongation in the exposed samples.

36 MATERIALS SCIENCE↗

Neutron elastic scattering on calcium isotopes from chiral nuclear optical potentials

We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Portable MPI Implementation of the SPAI Preconditioner in ISIS++

A parallel MPI implementation of the Sparse Approximate Inverse (SPAI) preconditioner is described. SPAI has proven to be a highly effective preconditioner, and is inherently parallel because it computes columns (or rows) of the preconditioning matrix independently. However, there are several problems that must be addressed for an efficient MPI implementation: load balance, latency hiding, and the need for one-sided communication. The effectiveness, efficiency, and scaling behavior of our implementation will be shown for different platforms.

Barnard, Stephen T.↗

Nucleation, growth, and strain relaxation of lattice-mismatched 3-5 semiconductor epitaxial layers

We have investigated the early stages of evolution of highly strained 2-D InAs layers and 3-D InAs islands grown by metal-organic chemical vapor deposition (MOCVD) on (100) and (111)B GaAs substrates. The InAs epilayer/GaAs substrate combination has been chosen because the lattice-mismatch is severe (approximately 7.2 percent), yet these materials are otherwise very similar. By examining InAs-on-GaAs composites instead of the more common In(x)Ga(1-x)As alloy we remove an additional degree of freedom (x) and thereby simplify data interpretation. A matrix of experiments is described in which the MOCVD growth parameters - susceptor temperature, Thin flux, and AsH3 flux - have been varied over a wide range. Scanning electron microscopy, atomic force microscopy, transmission electron microscopy, and electron microprobe analysis have been employed to observe the thin film surface morphology. In the case of 3-D growth, we have extracted activation energies and power-dependent exponents that characterize the nucleation process. As a consequence, optimized growth conditions have been identified for depositing approximately 250 A thick (100) and (111)B oriented InAs layers with relatively smooth surfaces. Together with preliminary data on the strain relaxation of these layers, the above results on the evolution of thin InAs films indicate that the (111)B orientation is particularly promising for yielding lattice-mismatched films that are fully relaxed with only misfit dislocations at the epilayer/substrate interface.

Welser, R. E.↗

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗

Aspects of the history of 66095 based on trace elements in clasts and whole rock

Halogens, P, U and Na are reported in anorthositic and basaltic clasts and matrix from rusty rock 66095. Large fractions of Cl and Br associated with the separated phases from 66095 are soluble in H2O. Up to two orders of magnitude variation in concentrations of these elements in the breccia components and varying H2O-soluble Cl/Br ratios indicate different sources of volatiles. An approximately constant ratio of the H2O- to 0.1 M HNO3-soluble Br in the various components suggests no appreciable alteration in the original distributions of this element in the breccia forming processes. Up to 50% or more of the phosphorus and of the non-H2O-soluble Cl was dissolved from most of the breccia components by 0.1 M HNO3. Clast and matrix residues from the leaching steps contain, in most cases, the Cl/P2O5 ratio found in 66095 whole rock and in a number of other Apollo 16 samples. Evidence that phosphates are the major P-phases in the breccia is based on the 0.1 M acid solubility of Cl and P in the matrix sample and on elemental concentrations which are consistent with those of KREEP.

Jovanovic, S.↗

Cell shape, cytoskeletal mechanics, and cell cycle control in angiogenesis

Capillary endothelial cells can be switched between growth and differentiation by altering cell-extracellular matrix interactions and thereby, modulating cell shape. Studies were carried out to determine when cell shape exerts its growth-regulatory influence during cell cycle progression and to explore the role of cytoskeletal structure and mechanics in this control mechanism. When G0-synchronized cells were cultured in basic fibroblast growth factor (FGF)-containing defined medium on dishes coated with increasing densities of fibronectin or a synthetic integrin ligand (RGD-containing peptide), cell spreading, nuclear extension, and DNA synthesis all increased in parallel. To determine the minimum time cells must be adherent and spread on extracellular matrix (ECM) to gain entry into S phase, cells were removed with trypsin or induced to retract using cytochalasin D at different times after plating. Both approaches revealed that cells must remain extended for approximately 12-15 h and hence, most of G1, in order to enter S phase. After this restriction point was passed, normally 'anchorage-dependent' endothelial cells turned on DNA synthesis even when round and in suspension. The importance of actin-containing microfilaments in shape-dependent growth control was confirmed by culturing cells in the presence of cytochalasin D (25-1000 ng ml-1): dose-dependent inhibition of cell spreading, nuclear extension, and DNA synthesis resulted. In contrast, induction of microtubule disassembly using nocodazole had little effect on cell or nuclear spreading and only partially inhibited DNA synthesis. Interestingly, combination of nocodazole with a suboptimal dose of cytochalasin D (100 ng ml-1) resulted in potent inhibition of both spreading and growth, suggesting that microtubules are redundant structural elements which can provide critical load-bearing functions when microfilaments are partially compromised. Similar synergism between nocodazole and cytochalasin D was observed when cytoskeletal stiffness was measured directly in living cells using magnetic twisting cytometry. These results emphasize the importance of matrix-dependent changes in cell and nuclear shape as well as higher order structural interactions between different cytoskeletal filament systems for control of capillary cell growth during angiogenesis.

Non-NASA Center↗

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗

First-principles study of electron transport in ScN

We investigate the conduction-band structure and electron mobility in rocksalt ScN based on density functional theory. The first-principles band structure allows us to obtain band velocities and effective masses as a function of energy. Electron-phonon scattering is assessed by explicitly computing the q-dependent electron-phonon matrix elements, with the inclusion of the long-range electrostatic interaction. The influence of free-carrier screening on the electron transport is assessed using the random phase approximation. We find a notable enhancement of electron mobility when the carrier concentration exceeds 10 20 cm -3 . We calculate the room-temperature electron mobility in ScN to be 587 cm2 /Vs at low carrier concentrations. When the carrier concentration is increased, the electron mobility starts to decrease significantly around n = 10 19 cm -3 , and drops to 240 cm 2 /Vs at n = 10 21 cm -3 . We also explore the influence of strain in (111)- and (100)-oriented ScN films. For (111) films, we find that a 1.0% compressive epitaxial strain increases the in-plane mobility by 72 cm 2 /Vs and the out-of-plane mobility by 50 cm 2 /Vs. For (100) films, a 1.0% compressive epitaxial strain increases the out-of-plane mobility by as much as 172 cm2 /Vs, but has a weak impact on the in-plane mobility. Our study sheds light on electron transport in ScN at different electron concentrations and shows how strain engineering could increase the electron mobility.

36 MATERIALS SCIENCE↗

Reduced-Order Models of Static Power Grids based on Spectral Clustering

For large-scale interconnected power systems that cover large geographical areas, certain electrical studies are required so that appropriate decisions ensure system reliability and low cost. For such studies, it is often neither practical nor necessary to model in detail the entire power system, which is increasingly complex due to a more diverse range of grid assets to choose from in both short and long-term planning. The goal of this paper is to present a methodology to reduce the order of large-scale power networks based on spectral graph theory given that current methods for static network reduction are not scalable. A brief analysis of some spectral clustering properties to determine which graph Laplacian matrix should be used and why is included. The analysis shows that the utilization of the normalized graph Laplacian is more advantageous for clustering purposes. Techniques are proposed to approximate cost functions for the aggregated generators. This is done via linear regression. The reduced-order model obtained with the proposed methodology has an accuracy above 94% and solves the scalability issue commonly present in other reduction methods. If the utilization of the reduced-order model is either constrained to load levels above mid-peak demand, or cost functions of aggregated units are approximated via a piecewise quadratic approach, then the error distribution is in the order of 10^-3. .

Baquedano-Aguilar, Mario D.↗

Fabrication of graphite/epoxy cases for orbit insertion motors

The fabrication procedures are described for filament-wound rocket motor cases, approximately 26.25 inches long by 25.50 inches diameter, utilizing graphite fibers. The process utilized prepreg tape which consists of Fortafil 4-R fibers in the E-759 epoxy resin matrix. This fabrication effect demonstrated an ability to fabricate high quality graphite/epoxy rocket motor cases in the 26.25 inch by 25.50 inch size range.

Schmidt, W. W.↗

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗

The effect of particle size on the electrical conductivity of CuCl (Al2O3) composites

The conductivity of CuCl containing Al2O3 of 0.06, 0.3, 1, 3, 8, and 15 micron sized particles was measured between 25 and 390 C. Conductivity was enhanced for the 0.06 and 0.3 sized Al2O3 particles for temperatures below approximately 250 C. The maximum enhancement occurred at 10 m/o of 0.06 micron Al2O3 at 25 C. Uncertain degrees of agglomeration as well as the grain size of the matrix were found to be significant.

Chang, M. R.-W.↗