Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 685 records · Page 38

Mechanical Properties of Nanostructured Materials Determined Through Molecular Modeling Techniques

The potential for gains in material properties over conventional materials has motivated an effort to develop novel nanostructured materials for aerospace applications. These novel materials typically consist of a polymer matrix reinforced with particles on the nanometer length scale. In this study, molecular modeling is used to construct fully atomistic models of a carbon nanotube embedded in an epoxy polymer matrix. Functionalization of the nanotube which consists of the introduction of direct chemical bonding between the polymer matrix and the nanotube, hence providing a load transfer mechanism, is systematically varied. The relative effectiveness of functionalization in a nanostructured material may depend on a variety of factors related to the details of the chemical bonding and the polymer structure at the nanotube-polymer interface. The objective of this modeling is to determine what influence the details of functionalization of the carbon nanotube with the polymer matrix has on the resulting mechanical properties. By considering a range of degree of functionalization, the structure-property relationships of these materials is examined and mechanical properties of these models are calculated using standard techniques.

Clancy, Thomas C.↗

Fracture of Carbon Nanotube - Amorphous Carbon Composites: Molecular Modeling

Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

Study of X-ray diffraction on aging alloys with different degrees of volume discrepancy between the precipitation phase and the matrix

Study of the diffraction patterns of X rays diffracted by dispersion-hardened alloys with various degrees of difference between the specific atomic volumes of their precipitation phases and matrices. Experiments were performed on aging alloys KhN77TIuR (nimonic) and KhN67VMTIuO, which have the same types of gamma solid solutions and precipitating gamma prime phases, but differ considerably with respect to the volume factor. It is found that in the decay of supersaturated solid solutions for which the difference between the specific atomic volumes of the precipitation phase and the matrix is small (alloy KhN67VMTIuO) the nature of the intensity distribution of the diffraction lines remains practically unchanged regardless of the aging, while in cases where the volume discrepancy is large (alloy KhN77TIuR) significant changes in the appearance of the diffraction lines are recorded.-

Gitgarts, M. I.↗

Spectral form factor of a quantum spin glass

It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p-spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The S-matrix of 2D type 0B string theory. Part II. D-instanton effects

We study the effect of D-instantons on closed string scattering amplitudes in the two-dimensional type 0B string theory from the worldsheet perspective. We find that the contribution from a pair of D-instanton and anti-D-instanton to the closed string reflection amplitude, with a suitable contour prescription for the integration over the D-instanton moduli space, agrees with the corresponding leading non-perturbative corrections in the proposed dual matrix quantum mechanics. This analysis is further extended to thermal observables defined at finite temperature. The infrared divergence in charged processes is understood through the measure factor for charged D-instantons, and can be treated with spacetime dimensional regularization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparison of δ 13 C analyses of individual foraminifer ( Orbulina universa ) shells by secondary ion mass spectrometry and gas source mass spectrometry

Rationale: The use of secondary ion mass spectrometry (SIMS) to perform micrometer-scale in situ carbon isotope (δ 13 C) analyses of shells of marine microfossils called planktic foraminifers holds promise to explore calcification and ecological processes. The potential of this technique, however, cannot be realized without comparison to traditional whole-shell δ 13 C values measured by gas source mass spectrometry (GSMS). Methods: Paired SIMS and GSMS δ 13 C values measured from final chamber fragments of the same shell of the planktic foraminifer Orbulina universa are compared. The SIMS–GSMS δ 13 C differences (Δ 13 C SIMS-GSMS ) were determined via paired analysis of hydrogen peroxide-cleaned fragments of modern cultured specimens and of fossil specimens from deep-sea sediments that were either untreated, sonicated, and cleaned with hydrogen peroxide or vacuum roasted. After treatment, fragments were analyzed by a CAMECA IMS 1280 SIMS instrument and either a ThermoScientific MAT-253 or a Fisons Optima isotope ratio mass spectrometer (GSMS). Results: Paired analyses of cleaned fragments of cultured specimens (n = 7) yield no SIMS–GSMS δ 13 C difference. However, paired analyses of untreated (n = 18) and cleaned (n = 12) fragments of fossil shells yield average Δ 13 C SIMS-GSMS values of 0.8‰ and 0.6‰ (±0.2‰, 2 SE), respectively, while vacuum roasting of fossil shell fragments (n = 11) removes the SIMS–GSMS δ 13 C difference. Conclusions: The noted Δ 13 C SIMS-GSMS values are most likely due to matrix effects causing sample–standard mismatch for SIMS analyses but may also be a combination of other factors such as SIMS measurement of chemically bound water. The volume of material analyzed via SIMS is ~10 5 times smaller than that analyzed by GSMS; hence, the extent to which these Δ 13 C SIMS-GSMS values represent differences in analyte or instrument factors remains unclear.

58 GEOSCIENCES↗

Randomized Cholesky Preconditioning for Graph Partitioning Applications

Graph partitioning has emerged as an area of interest due to its use in various applications in computational research. One way to partition a graph is to solve for the eigenvectors of the corresponding graph Laplacian matrix. This project focuses on the eigensolver LOBPCG and the evaluation of a new preconditioner: Randomized Cholesky Factorization (rchol). This proconditioner was tested for its speed and accuracy against other well-known preconditioners for the method. After experiments were run on several known test matrices, rchol appears to be a better preconditioner for structured matrices. This research was sponsored by National Nuclear Security Administration Minority Serving Institutions Internship Program (NNSA-MSIIP) and completed at host facility Sandia National Laboratories. As such, after discussion of the research project itself, this report contains a brief reflection on experience gained as a result of participating in the NNSA-MSIIP.

97 MATHEMATICS AND COMPUTING↗

Heat pipe temperature control utilizing a soluble gas absorption reservior

A new gas-controlled heat pipe design is described which uses a liquid matrix reservior, or sponge, to replace the standard gas reservior. Reservior volume may be reduced by a factor of five to ten for certain gas-liquid combinations, while retaining the same level of temperature control. Experiments with ammonia, butane, and carbon dioxide control gases with methanol working fluid are discussed.

Saaski, E. W.↗

Accelerated fatigue durability of a high performance composite

The fatigue behavior of multidirectional graphite-epoxy laminates was analyzed theoretically and experimentally in an effort to establish an accelerated testing methodology. Analysis of the failure mechanism in fatigue of the laminates led to the determination of the failure mode governing fracture. The nonlinear, cyclic-dependent shear modulus was used to calculate the changing stress field in the laminate during the fatigue loading. Fatigue tests were performed at three different temperatures: 25 C, 74 C, and 114 C. The prediction of the S-N curves was made based on the artificial static strength artificial static strength at a reference temperature and the fatigue functions associated with them. The prediction of an S-N curve at other temperatures was performed using shifting factors determined for the specific failure mode. For multidirectional laminates, different S-N curves at different temperatures could be predicted using these shifting factors. Different S-N curves at different temperatures occur only when the fatigue failure mode is matrix dominated. It was found that whenever the fatigue failure mode is fiber dominated, temperature, over the range investigated, had no influence on the fatigue life. These results permit the prediction of long-time, low temperature fatigue behavior from data obtained in short time, high temperature testing, for laminates governed by a matrix failure mode.

Rotem, A.↗

SiC or B4C-B/low strategic element content composite

Advancements in materials technologies are needed to provide the aerospace industry with alternate material options in the event of future strategic metal shortages and to optimize performance of engines. Composite materials are promising candidates for such an application. The potential of SiC and B4C-B filament reinforced low strategic element iron-base alloy content composites or use at 760 to 870 C is considered. A limiting factor towards developing this material for high temperature use has been the reaction between the filament and matrix material during fabrication and service which degrades filament strength. A low temperature fabrication process to limit filament/matrix reaction is being developed which involves the use of hollow cathode sputtering to coat the filaments with iron-base alloys of various compositions. An investigation is being conducted to determine the interfacial reaction effects of SiC and B4C-B filaments with iron-base alloys to develop an understanding of filament/matrix alloy compatibility for 760 to 870 C service.

Petrasek, D. W.↗

Component mode damping assignment techniques

A relation between the system modal damping matrix and the component modal damping matrix is derived from First Principles. An optimization problem is then formulated to select all the component modes' damping ratios that best satisfy the above derived relation. A weighting matrix is used in the cost functional to stress the relative importance of the diagonal terms in the damping matrix. Inequality constraints are also added to the optimization problem to pick only nonnegative component modes' damping factors. The optimization problem may be solved algebraically or iteratively. The proposed techniques are successfully used on a high order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

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.↗

Radioprotection of Human Cell Nuclear DNA by Polyamines: Radiosensitivity of Chromatin is Influenced by Tightly Bound Spermine

The polyamines putrescine (PUT) and spermine (SPM) were examined for their ability to protect human cell Deoxyribonucleic Acid (DNA) against the formation of radiation-induced double-strand breaks (DSBs). As observed previously, under conditions where polyamines were shown to be almost completely absent, association with nuclear matrix protein into a nucleoid, and organization into chromatin structure, protected DNA from induction of DSBs by factors of 4.5 and 95, respectively. At concentrations below 1 mM, PUT or SPM provided equivalent levels of protection to deproteinized nuclear DNA, consistent with their capacity to scavenge radiation-induced radicals. At constant ionic strength, 5 mM SPM protected deproteinized DNA and nucleoid DNA and DNA in nuclear chromatin by factors of 100 and 26, respectively. At 5 mM, SPM provided 15 times greater protection of deproteinized DNA than did PUT. Under physiologically relevant conditions, 5 mM SPM protected DNA in the intact nucleus from the induction of DSBs by a factor of 2 relative to DNA in the absence of SPM. Studies of SPM binding during cellular fractionation revealed that a significant fraction of the cellular SPM is tightly bound in the nucleus but can be removed by extended washing. Thus the association of SPM with nuclear chromatin appears to be a significant contributor to the resistance of the cell's DNA to the induction of DSBs.

Warters, Raymond L.↗

Examination of Surface Roughness on Light Scattering by Long Ice Columns by Use of a Two-Dimensional Finite-Difference Time-Domain Algorithm

Natural particles such as ice crystals in cirrus clouds generally are not pristine but have additional micro-roughness on their surfaces. A two-dimensional finite-difference time-domain (FDTD) program with a perfectly matched layer absorbing boundary condition is developed to calculate the effect of surface roughness on light scattering by long ice columns. When we use a spatial cell size of 1/120 incident wavelength for ice circular cylinders with size parameters of 6 and 24 at wavelengths of 0.55 and 10.8 mum, respectively, the errors in the FDTD results in the extinction, scattering, and absorption efficiencies are smaller than similar to 0.5%. The errors in the FDTD results in the asymmetry factor are smaller than similar to 0.05%. The errors in the FDTD results in the phase-matrix elements are smaller than similar to 5%. By adding a pseudorandom change as great as 10% of the radius of a cylinder, we calculate the scattering properties of randomly oriented rough-surfaced ice columns. We conclude that, although the effect of small surface roughness on light scattering is negligible, the scattering phase-matrix elements change significantly for particles with large surface roughness. The roughness on the particle surface can make the conventional phase function smooth. The most significant effect of the surface roughness is the decay of polarization of the scattered light.

Sun, W.↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Penny-shaped crack in a fiber-reinforced matrix

Using a slender inclusion model developed earlier, the elastostatic interaction problem between a penny-shaped crack and elastic fibers in an elastic matrix is formulated. For a single set and for multiple sets of fibers oriented perpendicularly to the plane of the crack and distributed symmetrically on concentric circles, the problem was reduced to a system of singular integral equations. Techniques for the regularization and for the numerical solution of the system are outlined. For various fiber geometries numerical examples are given, and distribution of the stress intensity factor along the crack border was obtained. Sample results showing the distribution of the fiber stress and a measure of the fiber-matrix interface shear are also included.

Narayanan, T. V.↗

Penny-shaped crack in a fiber-reinforced matrix

Using the slender inclusion model developed earlier the elastostatic interaction problem between a penny-shaped crack and elastic fibers in an elastic matrix is formulated. For a single set and for multiple sets of fibers oriented perpendicularly to the plane of the crack and distributed symmetrically on concentric circles the problem is reduced to a system of singular integral equations. Techniques for the regularization and for the numerical solution of the system are outlined. For various fiber geometries numerical examples are given and distribution of the stress intensity factor along the crack border is obtained. Sample results showing the distribution of the fiber stress and a measure of the fiber-matrix interface shear are also included.

Narayanan, T. V.↗