Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Matrix”

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 199 records · Page 11

Pareto-optimal target definition for multi-axis random vibration testing

In random vibration testing with multiple control channels, existing control laws require specification of a complete spectral density matrix at each control frequency. Spectral density matrices include autospectral densities on the diagonal and cross-spectral densities on the off-diagonal. In practice, the off-diagonal terms are often unknown, and recent vibration testing research has focused on fixing the diagonal and specifying the off-diagonal to minimize the required control energy, subject to a constraint that the target matrix is positive semidefinite. This paper shows that, even with a fixed diagonal, off-diagonal terms strongly affect control residuals. This overlooked effect occurs in both square and rectangular systems. By jointly considering input energy and control residuals, open-loop inputs are derived directly from the diagonal without specifying the off-diagonal terms. Vibration targets that can be used in closed-loop control are then derived using the optimal inputs, with positive semidefinite constraints applied during the derivation. The result is a set of Pareto-optimal control solutions. For each solution in the set, any other possible solution produces greater control error, greater input energy, or both. A balanced solution is selected automatically, though others can be chosen based on test needs. Simulations and experiments show that the proposed method outperforms state-of-the-art energy-minimizing approaches, achieving significant reductions in both control error and input energy.

Autospectral density↗

Grain Boundary Segregation Suppresses Local Short‐Range Ordering in Nanocrystalline High‐Entropy Alloys

Multi-principal-element alloys like high-entropy alloys (HEAs) have potential applications in many engineering fields due to their unique mechanical/functional properties. While HEAs are generally considered random solid solutions, recent studies revealed that they are prone to short-range-ordering (SRO) due to the complex multi-pair-wise interactions among the constituent elements. Meanwhile, SROs' evolution can sometimes be deleterious, and it is necessary to have control over their evolution. Examining the AlCoCrFe-Zr model alloy, long-range ordering occurs following the expectation of enthalpic predictions. Advanced characterization techniques—transmission electron microscopy, high-energy synchrotron X-ray diffraction/pair distribution function, and atom probe tomography, reveal that SRO is suppressed in as-milled and GB-decorated NC-(AlCoCrFe)100-xZrx (x = 0–1.5 atomic %). Warren-Cowley coefficient calculations are further used to validate the suppression of SRO. Besides the low segregation enthalpies of Cr, Fe, and Zr, and the high-mixing enthalpy of Cr and Fe, the short diffusion path to GBs due to high-GB density in the NC-HEAs and the higher energy state of the GBs than the matrix promotes GB-segregation that further alters the matrix chemistry and consequently disfavors SRO formation within the matrix. Despite the GB-segregation of Cr, Fe, and Zr, the matrices and GBs remain in a random solid solution.

36 MATERIALS SCIENCE↗

Morphology-Dependent Resonances of Spherical Droplets with Numerous Microscopic Inclusions

We use the recently extended superposition T-matrix method to study the behavior of a sharp Lorenz-Mie resonance upon filling a spherical micrometer-sized droplet with tens and hundreds of randomly positioned microscopic inclusions. We show that as the number of inclusions increases, the extinction cross-section peak and the sharp asymmetry-parameter minimum become suppressed, widen, and move toward smaller droplet size parameters, while ratios of diagonal elements of the scattering matrix exhibit sharp angular features indicative of a distinctly nonspherical particle. Our results highlight the limitedness of the concept of an effective refractive index of an inhomogeneous spherical particle.

resonance↗

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗

Parallel, iterative solution of sparse linear systems: Models and architectures

A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. A data transfer model predicting both the probability that data must be transferred between two tasks and the amount of data to be transferred is presented. This model is used to derive an execution time model for predicting parallel execution time and an optimal number of tasks given the dimension and sparsity of the coefficient matrix and the costs of computation, synchronization, and communication. The suitability of different parallel architectures for solving randomly sparse linear systems is discussed. Based on the complexity of task scheduling, one parallel architecture, based on a broadcast bus, is presented and analyzed.

Reed, D. A.↗

Parallel, iterative solution of sparse linear systems - Models and architectures

Solving large, sparse, linear systems of equations is a fundamental problem in large scale scientific and engineering computation. A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. A data transfer model predicting both the probability that data must be transferred between two tasks and the amount of data to be transferred is presented. This model is used to derive an execution time model for predicting parallel execution time and an optimal number of tasks given the dimension and sparsity of the coefficient matrix and the costs of computation, synchronization, and communication. The suitability of different parallel architectures for solving randomly sparse linear systems is discussed. Based on the complexity of task scheduling, one parallel architecture, based on a broadcast bus, is presented and analyzed.

Reed, D. A.↗

Metal Decorated Multi-Walled Carbon Nanotube/Polyimide Composites with High Dielectric Constants and Low Loss Factors

The measurement of observable electromagnetic phenomena in materials and their derived intrinsic electrical material properties are of prime importance in the discovery and development of material systems for electronic and aerospace applications. Nanocomposite materials comprised of metal decorated multi-walled carbon nanotubes (MWCNTs) were prepared by a facile method and characterized. Metal particles such as silver(Ag), platinum(Pt) and palladium(Pd) with diameters ranging from less than 5 to over 50 nanometers were distributed randomly on the MWCNTs. The present study is focused on silver decorated MWCNTs dispersed in a polyimide matrix. The Ag-containing MWCNTs were melt mixed into Ultem(TradeMark) and the mixture extruded as ribbons. The extruded ribbons exhibited a moderate to high degree of MWCNT alignment as determined by HRSEM. These ribbons were then fabricated into test specimens while maintaining MWCNT alignment and subsequently characterized for electrical and electromagnetic properties at 8-12 GHz. The results of the electromagnetic characterization showed that certain sample configurations exhibited a decoupling of the permittivity (epsilon ) and loss factor (epsilon") indicating that these properties could be tailored within certain limits. The decoupling and independent control of these fundamental electrical material parameters offers a new class of materials with potential applications in electronics, microwave engineering and optics.

Elliott, Holly A.↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A probabilistic formulation for fracture energy of continuous fibre-matrix composites

Since the fracture energy is a key design variable in brittle-matrix composites, the calculation of the energy and its probabilistic distribution becomes a prominent issue. A demonstration is presented of the way in which the statistical moments of the fracture-energy random variable may be directly calculated. The case considered is that in which continuous fibers possess periodic points of weakness along their length; the shear strength at the interface caused by the reentry of the fibers into the matrix will produce fiber-debonding.

Wetherhold, Robert C.↗

Comprehensive Thematic T-Matrix Reference Database: A 2014-2015 Update

The T-matrix method is one of the most versatile and efficient direct computer solvers of the macroscopic Maxwell equations and is widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper is the seventh update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2013. It also lists a number of earlier publications overlooked previously.

electromagnetic scattering↗

FT-IR microscopic mappings of early mineralization in chick limb bud mesenchymal cell cultures

Chick limb bud mesenchymal cells differentiate into chondrocytes and form a cartilaginous matrix in culture. In this study, the mineral formed in different areas within cultures supplemented with 4 mM inorganic phosphate, or 2.5, 5.0, and 10 mM beta-glycerophosphate (beta GP), was characterized by Fourier-transform infrared (FT-IR) microscopy. The relative mineral-to-matrix ratios, and distribution of crystal sizes at specific locations throughout the matrix were measured from day 14 to day 30. The only mineral phase detected was a poorly crystalline apatite. Cultures receiving 4 mM inorganic phosphate had smaller crystals which were less randomly distributed around the cartilage nodules than those in the beta GP-treated cultures. beta GP-induced mineral consisted of larger, more perfect apatite crystals. In cultures receiving 5 or 10 mM beta GP, the relative mineral-to-matrix ratios (calculated from the integrated intensities of the phosphate and amide I bands, respectively) were higher than in the cultures with 4 mM inorganic phosphate or in the in vivo calcified chick cartilage.

NASA Discipline Musculoskeletal↗

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗

Window convolution of the galaxy clustering bispectrum

In galaxy survey analysis, the observed clustering statistics do not directly match theoretical predictions but rather have been processed by a window function that arises from the survey geometry including the sky footprint, redshift-dependent background number density and systematic weights. While window convolution of the power spectrum is well studied, for the bispectrum with a larger number of degrees of freedom, it poses a significant numerical and computational challenge. In this work, we consider the effect of the survey window in the tripolar spherical harmonic decomposition of the bispectrum and lay down a formal procedure for their convolution via a series expansion of configuration-space three-point correlation functions, which was first proposed by Sugiyama et al. (2019). We then provide a linear algebra formulation of the full window convolution, where an unwindowed bispectrum model vector can be directly premultiplied by a window matrix specific to each survey geometry. To validate the pipeline, we focus on the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) luminous red galaxy (LRG) sample in the South Galactic Cap (SGC) in the redshift bin 0.4 ≤ z ≤ 0.6. We first perform convergence checks on the measurement of the window function from discrete random catalogues, and then investigate the convergence of the window convolution series expansion truncated at a finite of number of terms as well as the performance of the window matrix. This work highlights the differences in window convolution between the power spectrum and bispectrum, and provides a streamlined pipeline for the latter for current surveys such as DESI and the Euclid mission.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learning–Augmented Laser-Induced Breakdown Spectroscopy for Spectral Discrimination of Iron Oxalates

Enhanced characterization and phase identification of post-PUREX Pu Oxalates (PuOXA) are pivotal for nonproliferation and pre-detonation nuclear forensics. Despite significant advances in the characterization of PuO 2 samples, little is known about the impact of both the chemical structure and oxidation states of PuOXA (i.e., Pu(III) and Pu(IV)) have on optical emission signatures. Here, we demonstrate the analytical capabilities of laser-induced breakdown spectroscopy (LIBS) applied to Fe(II) and Fe(III) oxalate samples as surrogates for PuOXA, highlighting the discriminating features in the LIBS emission spectra arising from differences in the oxidation states within mixed FeOXA samples. We report the enhancement of spectral feature selection using Principal Component Analysis (PCA), which enables the analytical superiority of machine learning algorithms such as Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest Regression (RFR) over conventional univariate techniques for phase discrimination and chemometric analysis. Cluster analysis revealed how both matrix effects and laser ablation influence cluster separability by introducing spectral artifacts that misdirect the maximization of variance. PCA-selected emission lines were used in the regression models, demonstrating that both univariate and multivariate linear regression models (i.e., PLSR and SVR) can achieve acceptable performance, with machine learning models outperforming conventional calibration regressions. Furthermore, the application of non-linearly activated PCA-selected emission lines illustrates how simplifying the data while retaining captured variance enables the use of less complex and more computationally efficient models. Furthermore, this is particularly evident in the underperformance of RFR, which suffers from increased computational costs and overfitting owing to its high complexity.

Oxalates↗

Thermal expansion behavior of graphite/glass and graphite/magnesium

The thermal expansion behavior of n (+/- 8)s graphite fiber reinforced magnesium laminate and four graphite reinforced glass-matrix laminates (a unidirectional laminate, a quasi-isotropic laminate, a symmetric low angle-ply laminate, and a random chopped-fiber mat laminate) was determined, and was found, in all cases, to not be significantly affected by thermal cycling. Specimens were cycled up to 100 times between -200 F and 100 F, and the thermal expansion coefficients determined for each material as a function of temperature were found to be low. Some dimensional changes as a function of thermal cycling, and some thermal-strain hysteresis, were observed.

Tompkins, Stephen S.↗

Optics of Water Microdroplets with Soot Inclusions: Exact Versus Approximate Results

We use the recently generalized version of the multi-sphere superposition T-matrix method (STMM) to compute the scattering and absorption properties of microscopic water droplets contaminated by black carbon. The soot material is assumed to be randomly distributed throughout the droplet interior in the form of numerous small spherical inclusions. Our numerically-exact STMM results are compared with approximate ones obtained using the Maxwell-Garnett effective-medium approximation (MGA) and the Monte Carlo ray-tracing approximation (MCRTA). We show that the popular MGA can be used to calculate the droplet optical cross sections, single-scattering albedo, and asymmetry parameter provided that the soot inclusions are quasi-uniformly distributed throughout the droplet interior, but can fail in computations of the elements of the scattering matrix depending on the volume fraction of soot inclusions. The integral radiative characteristics computed with the MCRTA can deviate more significantly from their exact STMM counterparts, while accurate MCRTA computations of the phase function require droplet size parameters substantially exceeding 60.

cloud droplets↗

Comprehensive Thematic T-Matrix Reference Database: a 2017-2019 Update

Since its inception in the mid-1960s, the T-matrix method has been one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It has been widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This compilation is the ninth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2017. It also includes a few earlier publications that have so far been overlooked.

Michael I Mishchenko↗