Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 289 records · Page 16

R-Matrix Analysis of the Neutron-Induced Cross Sections on 143 Nd and 147,149 Sm Measured at LANSCE [Slides]

This presentation covers how an established used of R-Matrix codes is being implemented at LANSCE. This presentation displays how there are unique features accessible with detectors at LANSCE that can be used to perform more complete R-Matrix analysis. Additionally discussed is how spin separation is a powerful technique that can be applied in DANCE measurements. Unique method to measure transmission with DICER are discussed. Finally, transmission and capture data are being analyzed in this work. The study can be extended to other channels. The ongoing analysis of the 143 Nd data is presented in this work. The same method is being applied to analyze the 147,149 Sm data. The 143 Nd and 147,149 Sm work started in FY-23 funded by NCSP.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗

Stochastic microstructural analysis of ceramic matrix composites using a high-fidelity multiscale framework

This paper presents a high-fidelity multiscale framework to simulate the mechanical behavior of ceramic matrix composites (CMCs), accounting for the complex micromorphology captured using detailed material characterization. First, high-resolution micrographs are obtained for the specific carbon fiber silicon-carbide-nitride matrix (C/SiNC) CMC to characterize the variability of the architectural features and manufacturing-induced defects at the microscale. An image processing algorithm is then used to precisely estimate the size and distribution of all subscale features and defects from the micrographs. The information is then used to generate a three-dimensional stochastic representative volume element (SRVE) to reconstruct microscale constituents accounting for the variability. As a result, the generated SRVEs are simulated using the high-fidelity generalized method of cells (HFGMC) micromechanics theory to investigate the effects of defects on the elastic properties of C/SiNC CMCs.

Ceramic Matrix Composites↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

Scalable, Infiltration-Free Ceramic Matrix Composite Manufacturing

Manufacturing of ceramic matrix composites (CMCs) with carbon fiber and carbon matrix includes a time- and labor-intensive ceramic infiltration step that is responsible for more than half of the total manufacturing cost. To make CMCs cost-competitive in price-sensitive markets, like concentrated solar power, it is essential to develop a CMC manufacturing process that skips the ceramic infiltration process. In this work, we will show how a high-char-yield preceramic resin can eliminate the infiltration step while maintaining material density and evaluate the technoeconomic impact of our CMC manufacturing process. We monitor both morphology and porosity to determine the quality of CMCs made with different preceramic resins and evaluate the impact of polymer infiltration and pyrolyzing cycles.

matrix composite, polymer infiltration and pyrolys↗

Simplified, Infiltration-free Ceramic Matrix Composite Manufacturing

Manufacturing of ceramic matrix composites (CMCs) with carbon fiber and carbon matrix includes a time- and labor-intensive ceramic infiltration step that is responsible for more than half of the total manufacturing cost. To make CMCs cost-competitive in price-sensitive markets, like concentrated solar power, it is essential to develop a CMC manufacturing process that skips the ceramic infiltration process. In this work, we will show how a high-char-yield preceramic resin can eliminate the infiltration step while maintaining material density and evaluate the technoeconomic impact of our CMC manufacturing process. Here, we monitor both morphology and porosity to determine the quality of CMCs made with different preceramic resins and evaluate the impact of polymer infiltration and pyrolyzing cycles.

36 MATERIALS SCIENCE↗

Polymer-Infiltrated Metal–Organic Frameworks for Thin-Film Composite Mixed-Matrix Membranes with High Gas Separation Properties

Thin-film composite mixed-matrix membranes (TFC-MMMs) have potential applications in practical gas separation processes because of their high permeance (gas flux) and gas selectivity. In this study, we fabricated a high-performance TFC-MMM based on a rubbery comb copolymer, i.e., poly(2-[3-(2H-benzotriazol-2-yl)-4-hydroxyphenyl] ethyl methacrylate)-co-poly(oxyethylene methacrylate) (PBE), and metal–organic framework MOF-808 nanoparticles. The rubbery copolymer penetrates through the pores of MOF-808, thereby tuning the pore size. In addition, the rubbery copolymer forms a defect-free interfacial morphology with polymer-infiltrated MOF-808 nanoparticles. Consequently, TFC-MMMs (thickness = 350 nm) can be successfully prepared even with a high loading of MOF-808. As polymer-infiltrated MOF is incorporated into the polymer matrix, the PBE/MOF-808 membrane exhibits a significantly higher CO 2 permeance (1069 GPU) and CO 2 /N 2 selectivity (52.7) than that of the pristine PBE membrane (CO 2 permeance = 431 GPU and CO 2 /N 2 selectivity = 36.2). Therefore, the approach considered in this study is suitable for fabricating high-performance thin-film composite membranes via polymer infiltration into MOF pores.

59 BASIC BIOLOGICAL SCIENCES↗

Graph Partitioning and Sparse Matrix Ordering using Reinforcement Learning and Graph Neural Networks

We present a novel method for graph partitioning, based on reinforcement learning and graph convolutional neural networks. Our approach is to recursively partition coarser representations of a given graph. The neural network is implemented using SAGE graph convolution layers, and trained using an advantage actor critic (A2C) agent. We present two variants, one for finding an edge separator that minimizes the normalized cut or quotient cut, and one that finds a small vertex separator. The vertex separators are then used to construct a nested dissection ordering to permute a sparse matrix so that its triangular factorization will incur less fill-in. The partitioning quality is compared with partitions obtained using METIS and SCOTCH, and the nested dissection ordering is evaluated in the sparse solver SuperLU. Our results show that the proposed method achieves similar partitioning quality as METIS and SCOTCH. Furthermore, the method generalizes across different classes of graphs, and works well on a variety of graphs from the SuiteSparse sparse matrix collection.

97 MATHEMATICS AND COMPUTING↗

$\mathcal{R}$ -matrix analysis of η + nat Cl reactions relevant to molten-salt reactor designs [Abstract]

Performed with support from the U.S. Nuclear Criticality Safety Program in an effort to provide improved chlorine cross section and corresponding covariance data, the R-matrix analysis of neutron induced reactions for two stable chlorine isotopes, 35,37 Cl, was performed in the energy range of thermal up to 1.2 MeV. Starting from the repository of the ENDF/B-VIII.0 library and following recent measurement campaigns, this work represents a significant improvement in the evaluation of the ( η , p ) reaction channel. The evaluation methodology uses the $\mathcal{R}$ -matrix code SAMMY to generate a set of Reich-Moore resonance parameters. As predicted by recently measured 35 Cl( η , p ) data, the presented evaluation features a dramatic increase in the magnitude of the ( η , p ) reaction channel over the ENDF/B-VIII.0 and previous nuclear data ENDF/B released libraries. Together with details of the evaluation procedure, the impact of the updates in the outgoing proton emissions on reactivity coefficients for different molten-salt reactor designs is presented and discussed.

charged-particle↗

Perovskite Grown in Gallium Nitride Nanowire Matrix for Stable and High‐Efficiency X‐Ray Detection

High-quality quasi-2D perovskites in a GaN nano-wire matrix are grown to build a 3D hetero-structure for high-performance X-ray sensing. In the 3D hetero-structure, GaN nano-wire matrix serves as an n-type charge collector that can rapidly extract carriers through the bulk film of the perovskite layer. Together with a p-type top electrode, a p–i–n diode with the 3D hetero-structure is built, that exhibits a rectified current–voltage characteristic. After analyzing the interface energy alignment, it is found that the fermi levels of the perovskite and GaN are aligned in the dark, and a quasi-fermi level splits upon illumination, introducing a built-in electrical field at the interface. As a result, strong photo-induced current is observed from the diode without an external field. Finally, the 3D diode for X-ray detection demonstration is used, revealing a sensitivity of 308.9 µC Gy air -1 cm -2 at an exceptionally low applied field of 0.125 V µm -1 . The X-ray-induced signal from the 3D diode is stable after 155 cycles of X-ray irradiation under a constant electric field. This demonstration informs a new 3D architecture for high-performance X-ray sensing, and it shows that GaN is a robust n-type interface for perovskite optoelectronic devices.

36 MATERIALS SCIENCE↗

BAD2matrix: Phylogenomic matrix concatenation, indel coding, and more

Common steps in phylogenomic matrix production include biological sequence concatenation, morphological data concatenation, insertion/deletion (indel) coding, gene content (presence/absence) coding, removing uninformative characters for parsimony analysis, recording with reduced amino acid alphabets, and occupancy filtering. Existing software does not accomplish these tasks on a phylogenomic scale using a single program. BAD2matrix is a Python script that performs the above-mentioned steps in phylogenomic matrix construction for DNA or amino acid sequences as well as morphological data. The script works in UNIX-like environments (e.g., LINUX, MacOS, Windows Subsystem for LINUX).

59 BASIC BIOLOGICAL SCIENCES↗

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗

Sparse Matrix-Based HPC Tomography

Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Furthermore, flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. Furthermore, the solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide.

97 MATHEMATICS AND COMPUTING↗

New oscillation and scattering constraints on the tau row matrix elements without assuming unitarity

The tau neutrino is the least well measured particle in the Standard Model. Most notably, the tau neutrino row of the lepton mixing matrix is quite poorly constrained when unitarity is not assumed. In this paper, we identify data sets involving tau neutrinos that improve our understanding of the tau neutrino part of the mixing matrix, in particular ν ${\tau}$ appearance in atmospheric neutrinos. We present new results on the elements of the tau row leveraging existing constraints on the electron and muon rows for the cases of unitarity violation, with and without kinematically accessible steriles. We also show the expected sensitivity due to upcoming experiments and demonstrate that the tau neutrino row precision may be comparable to the muon neutrino row in a careful combined fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The soft $ \mathcal{S} $-matrix in gravity

The gravitational $ \mathcal{S}$-matrix defined with an infrared (IR) cutoff factorizes into hard and soft factors. The soft factor is universal and contains all the IR and collinear divergences. Here we show, in a momentum space basis, that the intricate expression for the soft factor is fully reproduced by two boundary currents, which live on the celestial sphere. The first of these is the supertranslation current, which generates spacetime supertranslations. The second is its symplectic partner, the Goldstone current for spontaneously broken supertranslations. The current algebra has an off-diagonal level structure involving the gravitational cusp anomalous dimension and the logarithm of the IR cutoff. It is further shown that the gravitational memory effect is contained as an IR safe observable within the soft $ \mathcal{S}$-matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Emergent unitarity in de Sitter from matrix integrals

We study Jackiw-Teitelboim gravity with positive cosmological constant as a model for de Sitter quantum gravity. We focus on the quantum mechanics of the model at past and future infinity. There is a Hilbert space of asymptotic states and an infinite-time evolution operator between the far past and far future. This evolution is not unitary, although we find that it acts unitarily on a subspace up to non-perturbative corrections. These corrections come from processes which involve changes in the spatial topology, including the nucleation of baby universes. There is significant evidence that this 1+1 dimensional model is dual to a 0+0 dimensional matrix integral in the double-scaled limit. So the bulk quantum mechanics, including the Hilbert space and approximately unitary evolution, emerge from a classical integral. We find that this emergence is a robust consequence of the level repulsion of eigenvalues along with the double scaling limit, and so is rather universal in random matrix theory.

2D gravity↗

Distributed non-negative matrix factorization with determination of the number of latent features

The holistic analysis and understanding of the latent (that is, not directly observable) variables and patterns buried in large datasets is crucial for data-driven science, decision making and emergency response. Such exploratory analyses require devising unsupervised learning methods for data mining and extraction of the latent features, and non-negative matrix factorization (NMF) is one of the prominent such methods. NMF is based on compute-intense non-convex constrained minimization, which, for large datasets requires fast and distributed algorithms. However, current parallel implementations of NMF fail to estimate the number of latent features. In practice, identifying these features is both difficult and significant for pattern recognition and latent feature analysis, especially for large dense matrices. Here, we introduce a distributed NMF algorithm coupled with distributed custom clustering followed by a stability analysis on dense data, which we call DnMFk, to determine the number of latent variables. The results on synthetic data and the classical Swimmer data set demonstrate the accuracy of model determination while scaling nearly linearly across multiple processors for large data. Further, we employ DnMFk to determine the number of hidden features from a terabyte matrix.

97 MATHEMATICS AND COMPUTING↗