Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrices”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Pebble Tracking Transport Algorithm for Pebble Bed Reactor Analysis

Pebble bed reactor (PBR) designs pose unique challenges to modeling and simulation. The challenges originate from the multiple levels of heterogeneity partly due to random packing in and flow through the core of fuel pebbles. This paper proposes a novel pebble tracking transport (PTT) algorithm to address some of the challenges facing PBR analysis. Each individual pebble is tracked and modeled on an unstructured mesh of which the nodes are aligned with the pebble centroids. The neutron transport equation is solved on this pebble tracking mesh with pre-assembled elemental matrices. These elemental matrices are independent of the cross sections and the streaming directions, thus can be assembled before the calculation without introducing too much memory overhead especially when combined with spatial domain decomposition. Numerical results show that these matrices can be assembled efficiently. PTT is implemented in Rattlesnake, the INL MOOSE-based multigroup radiation transport application. Both the angular and spatial convergence are investigated and the results indicate that low order angular and spatial resolution is typically sufficient for analysis. A complete work flow for PBR analysis with PTT including, pebble packing, meshing, cross section generation, solution of the transport equation, postprocessing, etc. is demonstrated with a simplified reactor with 41,048 pebbles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Compressed Expanded Natural Graphite (CENG) Processing for PCM Composites

The use of phase change materials (PCMs) in thermal energy storage applications has received considerable attention in recent decades. Organic PCMs are popular due to their high latent heat of fusion, noncorrosive properties, and relative stability over many charge and discharge cycles. A primary limitation of these materials is their low thermal conductivity. This has led researchers to develop various methods to increase thermal conductivity by seeding PCM with or impregnating them into more conductive materials. One method is to impregnate PCMs into compressed expanded natural graphite (CENG) matrices, which can improve thermal conductivity by a factor of 100. CENG matrices have received particular interest due to their low cost, high porosity, small (nano/micro) pore size, high pore density, high thermal conductivity, and ability to be compressed into many geometries. PCM/CENG matrix composites have been extensively studied; however, the effect that CENG processing has on PCM saturation and the overall matrix thermal conductivity has not been well investigated. This processing includes four major steps including graphite intercalation, thermal shock, compression, and PCM saturation. Intercalation involves soaking graphite flakes in sulfuric and/or nitric acid to intercalate the acids between the graphene layers. The graphite flakes are then subjected to a high-temperature thermal shock, during which the intercalated acid is gasified rapidly, pushing the graphene layers apart, resulting in accordion-shaped graphite "worms". The "worms" are then compacted to a desired bulk density and then soaked with molten PCM until fully saturated. The properties of the produced CENG matrix, and its ability to allow PCM permeation, are sensitive to the processing parameters, namely, the thermal shock temperature and exposure time, as well as the matrix apparent density or porosity. Here, we study the effect of the thermal shock conditions necessary to expand intercalated graphite flakes on PCM saturation and the expanded graphite's thermal conductivity and morphology. We found that the thermal shock temperature exhibits the greatest influence. At greater shock temperatures, SEM images showed that expanded graphite worms exhibited greater density of pores, thus increasing total surface area within the matrices. Increasing thermal shock temperature yielded greater overall PCM saturation, as well as an increased rate of saturation. Improvements in PCM saturation rate and overall saturation are obtained as the shock temperature is increased. Longer exposure to thermal shock also improves initial saturation rates and is beneficial if a shortened impregnation time is needed. Thermal shock conditions did not impact thermal conductivity; however, conductivity was largely affected by matrix porosity. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to that observed in previous studies.

compressed expanded natural graphite↗
DOE OSTI2012

Comparison of homogenized and enhanced diffusion solutions of model PWR problems

Model problem comparisons in slab geometry are made between two forms of homogenized diffusion theory and enhanced diffusion theory. The pin-cell discontinuity factors for homogenized diffusion calculations are derived from homogenized variational nodal P1 response matrices and from standard finite differencing. Enhanced diffusion theory consists of applying quasi-reflected interface conditions to reduce variational nodal Pn response matrices to one degree of freedom per interface, without homogenization within the cell. As expected both homogenized diffusion methods preserve reaction rates exactly if the discontinuity factors are derived from the P 11 reference solutions. If no reference lattice solution is available, discontinuity factors may be approximated from single cells with reflected boundary conditions; the computational effort is then comparable to calculating the enhanced diffusion response matrices. In this situation enhanced diffusion theory gives the most accurate results and finite difference discontinuity factors the least accurate. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

A Mobile Automated Tomographic Gamma Scanning System - 13231

Canberra Industries have recently designed and built a new automated Tomographic Gamma Scanning (TGS) system for mobile deployment. The TGS technique combines high-resolution gamma spectroscopy with low spatial resolution 3-dimensional image reconstruction to provide increased accuracy over traditional approaches for the assay of non-uniform source distributions in low-to medium-density, non-heterogeneous matrices. Originally pioneered by R. Estep at Los Alamos National Laboratory (LANL), the TGS method has been further developed and commercialized by Canberra Industries in recent years. The present system advances the state of the art on several fronts: it is designed to be housed in a standard cargo transport container for ease of transport, allowing waste characterization at multiple facilities under the purview of a single operator. Conveyor feed, drum rotator, and detector and collimator positioning mechanisms operated by programmable logic control (PLC) allow automated batch mode operation. The variable geometry settings can accommodate a wide range of waste packaging, including but not limited to standard 220 liter drums, 380 liter overpack drums, and smaller 20 liter cans. A 20 mCi Eu-152 transmission source provides attenuation corrections for drum matrices up to 1 g/cm{sup 3} in TGS mode; the system can be operated in Segmented Gamma Scanning (SGS) mode to measure higher density drums. To support TGS assays at higher densities, the source shield is sufficient to house an alternate Co-60 transmission source of higher activity, up to 250 mCi. An automated shutter and attenuator assembly is provided for operating the system with a dual intensity transmission source. The system's 1500 kg capacity rotator turntable can handle heavy containers such as concrete lined 380 liter overpack drums. Finally, data acquisition utilizes Canberra's Broad Energy Germanium (BEGE) detector and Lynx MCA, with 32 k channels, providing better than 0.1 keV/channel resolution to support both isotopic analysis with the MGA/MGAU software and a wide 3 MeV dynamic range. The calibration and verification of the system is discussed, and quantitative results are presented for a variety of drum types and matrices. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Facility Utilization and Risk Analysis for Remediation of Legacy Transuranic Waste at the Savannah River Site - 13572

Savannah River Nuclear Solutions (SRNS) completed the Accelerated TRU Project for remediating legacy waste at the Savannah River Site with significant cost and schedule efficiencies due to early identification of resources and utilization of risk matrices. Initial project planning included identification of existing facilities that could be modified to meet the technical requirements needed for repackaging and remediating the waste. The project schedule was then optimized by utilization of risk matrices that identified alternate strategies and parallel processing paths which drove the overall success of the project. Early completion of the Accelerated TRU Project allowed SRNS to pursue stretch goals associated with remediating very difficult TRU waste such as concrete casks from the hot cells in the Savannah River National Laboratory. Project planning for stretch goals also utilized existing facilities and the risk matrices. The Accelerated TRU project and stretch goals were funded under the American Recovery and Reinvestment Act (ARRA). (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Using Green's Functions to Compute 3 Dimensional Neutron Flux Maps - 16016

The knowledge of the radionuclide content of radioactive waste is of utmost importance for safety and waste management reasons. Numerical simulations are used by EDF-DP2D to anticipate dismantling and radioactive waste management. The main numerical activation schemes include the computing of a 3 dimensional neutron flux map, obtained on the basis of a neutron propagation calculation. The codes used solve the transport equation called the Boltzmann equation without simplification. The Monte Carlo method is used to solve the Boltzmann equation. Random series of numbers are used to simulate the lives of millions of neutrons. The codes follow each neutron individually, from its birth to its disappearance by leakage, absorption, or fission, and then the neutron flux map is calculated at the nominal power rating conditions. The main drawback of the Monte Carlo codes is the large computing time needed to converge within acceptable statistical criteria. This extended time period is incompatible with the multiple fuel managements used by EDF to operate its NPPs because it is necessary to make a Monte Carlo calculation for each fuel loading pattern. The solution developed by EDF-DP2D is to use the importance factors of Green's Functions to simplify the computing of the 3 dimensional neutron flux maps. Physicists talk about Green's Functions where mathematicians talk about elementary solutions of linear equations with constant coefficients. With an activation scheme based on the Green's Functions approach the 3D neutron flux is obtained on the basis of the convolution of the normalized source contributions (i.e. normalized response matrices) with the neutrons emitted by the fuel assemblies. The response matrices to unitary sources are computed with the Monte Carlo code MCNP for different tallies of interest and the neutron sources emitted by the fuel are computed with a core code for different fuel loading patterns. The principal interest of this approach is the calculation of the response matrices to unitary sources that can be used for all power distributions, thus limiting the Monte Carlo computing time. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

QuadSync: Quadrifocal tensor synchronization via Tucker decomposition

In structure from motion, quadrifocal tensors capture more information than their pairwise counterparts (essential matrices), yet they have often been thought of as impractical and only of theoretical interest. In this work, we challenge such beliefs by providing a new framework to recover n cameras from the corresponding collection of quadrifocal tensors. We form the block quadrifocal tensor and show that it admits a Tucker decomposition whose factor matrices are the stacked camera matrices, and which thus has a multilinear rank of (4,4,4,4) independent of n. We develop the first synchronization algorithm for quadrifocal tensors, using Tucker decomposition, alternating direction method of multipliers, and iteratively reweighted least squares. We further establish relationships between the block quadrifocal, trifocal, and bifocal tensors, and introduce an algorithm that jointly synchronizes these three entities. Numerical experiments demonstrate the effectiveness of our methods on modern datasets, indicating the potential and importance of using higher-order information in synchronization.

Miao, Daniel [University of Minnesota]↗
DOE OSTI2026

Two-Tower Quantum Matrix Chain Multiplication: Trading Qubits for Depth

Matrix chain multiplication -- computing $\mathcal{W} = M^{(0)}\cdots M^{(K-1)}$ where $M^{(k)} \in \mathbb{R}^{P_k \times P_{k+1}}$-- arises in scientific computing, machine learning, and graph analysis. Despite the importance of this problem, for chains of distinct matrices, the classical number of operations grows linearly with the chain length $K$ and polynomially in the matrix dimensions. We present \emph{Two-Tower Matrix Multiplication}, a quantum subroutine that encodes the product $\mathcal{W}$ of the $K$ matrices into a quantum state in circuit depth $\mathcal{O}(\max_{k} \mathrm{polylog} (P_k P_{k+1}))$, which is independent of~$K$ within the QRAM-based state-preparation model, whereas the qubit count is $\mathcal{O}\bigl(\sum_{k} \log P_k \bigr)$; the total gate count remains linear in $K$, so the gain is in the circuit depth. The construction interleaves state-preparation operators across two layers; within each layer, all operators act on disjoint registers and execute in parallel. This subroutine can be specialized for the chain-vector case, which computes the product of $K-1$ matrices applied to a vector. We prove the correctness of the subroutine for all $K$ and provide two implementations using the Qiskit and QCLAB frameworks. The subroutine is applicable to any downstream quantum algorithm that operates on a matrix encoded in the statevector, including norm estimation, graph-matrix powers, linear system solving, and quantum machine learning kernels.

Antonioli, Giacomo [Pisa U.] (ORCID:00090000668703↗
DOE OSTI2013

Combinatorial theory of the semiclassical evaluation of transport moments. I. Equivalence with the random matrix approach

To study electronic transport through chaotic quantum dots, there are two main theoretical approaches. One involves substituting the quantum system with a random scattering matrix and performing appropriate ensemble averaging. The other treats the transport in the semiclassical approximation and studies correlations among sets of classical trajectories. There are established evaluation procedures within the semiclassical evaluation that, for several linear and nonlinear transport moments to which they were applied, have always resulted in the agreement with random matrix predictions. We prove that this agreement is universal: any semiclassical evaluation within the accepted procedures is equivalent to the evaluation within random matrix theory. The equivalence is shown by developing a combinatorial interpretation of the trajectory sets as ribbon graphs (maps) with certain properties and exhibiting systematic cancellations among their contributions. Remaining trajectory sets can be identified with primitive (palindromic) factorisations whose number gives the coefficients in the corresponding expansion of the moments of random matrices. The equivalence is proved for systems with and without time reversal symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Optimal system size for complex dynamics in random neural networks near criticality

In this article, we consider a model of dynamical agents coupled through a random connectivity matrix, as introduced by Sompolinsky et al. [Phys. Rev. Lett. 61(3), 259–262 (1988)] in the context of random neural networks. When system size is infinite, it is known that increasing the disorder parameter induces a phase transition leading to chaotic dynamics. We observe and investigate here a novel phenomenon in the sub-critical regime for finite size systems: the probability of observing complex dynamics is maximal for an intermediate system size when the disorder is close enough to criticality. We give a more general explanation of this type of system size resonance in the framework of extreme values theory for eigenvalues of random matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Nanostructured Metal Oxide Catalysts via Building Block Syntheses

A broadly applicable methodology has been developed to prepare new single site catalysts on silica supports. This methodology requires of three critical components: a rigid building block that will be the main structural and compositional component of the support matrix; a family of linking reagents that will be used to insert active metals into the matrix as well as cross link building blocks into a three dimensional matrix; and a clean coupling reaction that will connect building blocks and linking agents together in a controlled fashion. The final piece of conceptual strategy at the center of this methodology involves dosing the building block with known amounts of linking agents so that the targeted connectivity of a linking center to surrounding building blocks is obtained. Achieving targeted connectivities around catalytically active metals in these building block matrices is a critical element of the strategy by which single site catalysts are obtained. This methodology has been demonstrated with a model system involving only silicon and then with two metal-containing systems (titanium and vanadium). The effect that connectivity has on the reactivity of atomically dispersed titanium sites in silica building block matrices has been investigated in the selective oxidation of phenols to benezoquinones. 2-connected titanium sites are found to be five times as active (i.e. initial turnover frequencies) than 4-connected titanium sites (i.e. framework titanium sites).

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

Convexity Conditions and the Legendre-Fenchel Transform for the Product of Finitely Many Positive Definite Quadratic Forms

While the product of finitely many convex functions has been investigated in the field of global optimization, some fundamental issues such as the convexity condition and the Legendre-Fenchel transform for the product function remain unresolved. Focusing on quadratic forms, this paper is aimed at addressing the question: When is the product of finitely many positive definite quadratic forms convex, and what is the Legendre-Fenchel transform for it? First, we show that the convexity of the product is determined intrinsically by the condition number of so-called 'scaled matrices' associated with quadratic forms involved. The main result claims that if the condition number of these scaled matrices are bounded above by an explicit constant (which depends only on the number of quadratic forms involved), then the product function is convex. Second, we prove that the Legendre-Fenchel transform for the product of positive definite quadratic forms can be expressed, and the computation of the transform amounts to finding the solution to a system of equations (or equally, finding a Brouwer's fixed point of a mapping) with a special structure. Thus, a broader question than the open 'Question 11' in Hiriart-Urruty (SIAM Rev. 49, 225-273, 2007) is addressed in this paper.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the Kräuter–Seifter Theorem on Permanent Divisibility

The paper investigates the divisibility of the permanent function of (1,−1)-matrices by different powers of 2. It is shown that the Kräuter–Seifter bound is the best possible one for generic(1,−1)-matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the Determinantal Range of Matrix Products

Let matrices A,C ∈ M{sub n} have eigenvalues α{sub 1}, . . ., α{sub n} and γ{sub 1}, . . . , γ{sub n}, respectively. The set of complex numbers D{sub C}(A) = {det(A−UCU*) : U ∈ M_n, U*U = I_n} is called the C-determinantal range of A. The paper studies various conditions under which the relation D{sub C}(R S) = D{sub C}(S R) holds for some matrices R and S.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On New Structures in the Theory of Fully Nonlinear Equations

We describe the current state of the theory of equations with m-Hessian stationary and evolution operators. It is quite important that new algebraic and geometric notions appear in this theory. In the present work, a list of those notions is provided. Among them, the notion of m-positivity of matrices is quite important; we provide a proof of an analog of Sylvester’s criterion for such matrices. From this criterion, we easily obtain necessary and sufficient conditions for existence of classical solutions of the first initial boundary-value problem for m-Hessian evolution equations. The asymptotic behavior of m-Hessian evolutions in a semibounded cylinder is considered as well.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Birational Darboux Coordinates on Nilpotent Coadjoint Orbits of Classical Complex Lie Groups, the Case of 2 × 2 Jordan Blocks

We consider the problem of constructing birational Darboux coordinates on nilpotent coadjoint orbits of the complex Lie groups SO(N, ℂ) and Sp(N, ℂ). The nilpotent case is the most difficult one. Difficulties arise if the Jordan form of matrices from the orbit under consideration contains Jordan blocks of sizes of different parity. The desired coordinates have been found on orbits consisting of matrices with 1 × 1 and 2 × 2 Jordan blocks. Explicit formulas for them are given in the paper.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The Effects of Sc Addition on the Microstructure and Mechanical Properties of Be-Al Alloy Fabricated by Induction Melting

This study investigated the alloying effects of 0.4 wt.% Sc on the microstructure, mechanical properties, and fracture behavior of a cast Be-Al alloy prepared by vacuum induction melting. The micromechanics of Sc-containing secondary-phase particles (SPs) were characterized by nanoindentation, and the fracture behavior of the alloy was investigated via tensile fractography and optical metallography. The addition of 0.4 wt.% Sc significantly modified the alloy microstructure by decreasing the secondary dendritic arm spacing (SDAS), weakening the dendritic structure of Be, and producing an equiaxed/cellular alloy microstructure with refined Be grains. Sc alloying increased the Vickers hardness by 36%, elastic modulus by 28%, yield strength by 35%, ultimate strength by 52%, but decreased the elongation from 3.6 to 2.9%. Sc-containing SPs had a higher elastic modulus, yield strength, microhardness, and a lower plasticity index than the Be and Al matrices, resulting in interface separation during the fracture process. SPs located at the Be/Al interfaces did not act as microcrack sources for the matrices, whereas highly faceted SPs within the Be grains induced microcrack formation and decreased the elongation. The failure mode of the cast Be-Al-0.4Sc alloy still involved the brittle cleavage failure of Be regions and ductile dimple failure of Al regions, which was not the case with the cast Be-Al alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact Damage Behavior of Vinyl Ester-, Epoxy-, and Nylon 6-Based Basalt Fiber Composites

Composite laminates reinforced by woven basalt fibers and based on different polymer matrices were prepared and mechanically characterized, in static and dynamic conditions, to compare the results and evaluate the influence of the matrix type. In this regard, a thermoplastic resin, polyamide 6 (PA6), and two thermosetting matrices, vinyl ester and epoxy grades, commonly used for high-performance applications, were considered. Indentation depth measurements carried out by a confocal microscope and ultrasound analyses performed on specimens subjected to falling dart tests at levels of impact energies lower than the perforation one permitted to evaluate the extent of plastic deformations and of subjected delaminations, respectively. These data have provided useful information about the involved damage mechanisms and highlighted interesting perspective of use for the eco-friendly thermoplastic sample.

36 MATERIALS SCIENCE↗