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 217 records · Page 12

LAW Simulant Recipes for Evaluation of Real-Time, In-Line Monitoring Instruments

Statistically designed recipes were generated for evaluation of in-line instrumentation for the Real-Time In-Line Monitoring (RTIM) program in support of the Hanford Direct-Feed Low Activity Waste (DFLAW) program. The composition matrices were prepared based on Low Activity Waste (LAW) feed vector, melter feed, and Effluent Management Facility (EMF) feed information provided by the flowsheet integration modelling for the DFLAW portion of the Waste Treatment Plant (WTP) mission. These composition matrices were converted to recipes for the LAW feed using the approaches for LAW simulant recipes developed by Savannah River National Laboratory (SRNL). Recipes for the glass forming chemicals (GFCs) were developed based on expected mineral forms to be used to for each element to be added to the recipe. The recipe matrices are shown in the Appendix.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Enhanced Hanford Low-Activity Waste Glass Property Data Development: Phase 5 and Phase 6

This report summarizes the data collected on two test matrices of low-activity waste (LAW) glass compositions intended to expand the composition-property database: Low-Activity Waste (LAW) Phase 5: Expansion of LAW Glass Composition Boundaries and LAW Phase 6: High PCT and VHT Response Glass Matrix. Both matrix glass compositions were statistically designed to expand the LAW glass composition region. The analyses performed on these glasses include chemical composition (for target compositional verification), density, viscosity, electrical conductivity, crystal fraction, container centerline cooling with crystal identification, the product consistency test (PCT) response, the vapor hydration test (VHT) response, and sulfur solubility. Because of the slightly different scope of the two matrices, not all methods were applied to both matrices. Specifically, the following measurements were taken only on the LAW Phase 5: Expansion of LAW Glass Composition Boundaries glasses: crystal fraction as a function of temperature, density (ρ), viscosity (η), and electrical conductivity (EC, e). Combined, these data contribute a significant amount, 51 glasses, to the database for high LAW loaded enhanced waste glasses. Most of these data are focused near the boundaries of acceptable PCT and VHT responses, where prediction uncertainties are most impactful.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data Driven Correlated Noise Simulation for the ICEBERG LArTPC

Accurate electronic-noise simulation is essential for low-energy physics in liquid-argon TPCs. More realistic noise modeling allows us to better tune reconstruction algorithms and more reliably assess and optimize signal-detection thresholds. We present a data-driven noise simulation framework developed for the ICEBERG test stand for DUNE that generates synthetic noise waveforms that reproduce both (i) the measured per-channel magnitude of the Fast Fourier Transform (FFT) and (ii) frequency-dependent channel-to-channel correlations observed in ICEBERG noise data. Using a dedicated noise-only dataset, we build a compact noise model containing per-channel FFT-magnitude targets together with a small set of band-wise cross-wire color matrices. White noise is generated in the frequency domain by drawing circular-symmetric complex Gaussian coefficients with random phases and scaling them to match the measured FFT-magnitude targets, and cross-wire correlations are subsequently imposed using the stored color matrices. The model and algorithm were integrated into the LArSoft + Wire-Cell Toolkit simulation chain and validated by comparing waveform structure, frequency-domain spectra, and band-limited correlation matrices from simulated noise and ICEBERG data. This approach can be extended to other LArTPC operating conditions.

Ghosh, Avik [Iowa State U.]↗

Genomic Metrics Applied to Rhizobiales ( Hyphomicrobiales ): Species Reclassification, Identification of Unauthentic Genomes and False Type Strains

Taxonomic decisions within the order Rhizobiales have relied heavily on the interpretations of highly conserved 16S rRNA sequences and DNA–DNA hybridizations (DDH). Currently, bacterial species are defined as including strains that present 95–96% of average nucleotide identity (ANI) and 70% of digital DDH (dDDH). Thus, ANI values from 520 genome sequences of type strains from species of Rhizobiales order were computed. From the resulting 270,400 comparisons, a ≥95% cut-off was used to extract high identity genome clusters through enumerating maximal cliques. Coupling this graph-based approach with dDDH from clusters of interest, it was found that: (i) there are synonymy between Aminobacter lissarensis and Aminobacter carboxidus, Aurantimonas manganoxydans and Aurantimonas coralicida, “Bartonella mastomydis,” and Bartonella elizabethae, Chelativorans oligotrophicus, and Chelativorans multitrophicus, Rhizobium azibense, and Rhizobium gallicum, Rhizobium fabae, and Rhizobium pisi, and Rhodoplanes piscinae and Rhodoplanes serenus; (ii) Chelatobacter heintzii is not a synonym of Aminobacter aminovorans; (iii) “Bartonella vinsonii” subsp. arupensis and “B. vinsonii” subsp. berkhoffii represent members of different species; (iv) the genome accessions GCF_003024615.1 (“Mesorhizobium loti LMG 6125T”), GCF_003024595.1 (“Mesorhizobium plurifarium LMG 11892T”), GCF_003096615.1 (“Methylobacterium organophilum DSM 760T”), and GCF_000373025.1 (“R. gallicum R-602 spT”) are not from the genuine type strains used for the respective species descriptions; and v) “Xanthobacter autotrophicus” Py2 and “Aminobacter aminovorans” KCTC 2477T represent cases of misuse of the term “type strain”. Aminobacter heintzii comb. nov. and the reclassification of Aminobacter ciceronei as A. heintzii is also proposed. To facilitate the downstream analysis of large ANI matrices, we introduce here ProKlust (“Prokaryotic Clusters”), an R package that uses a graph-based approach to obtain, filter, and visualize clusters on identity/similarity matrices, with settable cut-off points and the possibility of multiple matrices entries.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

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

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↗

Efficient Unitary Designs from Random Sums and Permutations

A unitary k-design is an ensemble of unitaries that matches the first k moments of the Haar measure. In this work, we provide two efficient constructions of k-designs on n-qubits using new random matrix theory techniques. Our first construction is based on exponentiating sums of random i.i.d. Hermitian matrices and uses O(k2n2)-many gates. In the spirit of central limit theorems, we show that this random sum approximates the Gaussian Unitary Ensemble (GUE). We then show that the product of just two exponentiated GUE matrices is already approximately Haar random. Our second construction is based on products of exponentiated sums of random permutations and uses Õ(k poly (n)) many gates. The k dependence is optimal (up to polylogarithmic factors) and is inherited from the efficiency of existing k-wise independent permutations. Furthermore, replacing random permutations with quantum-secure pseudorandom permutations (PRPs), we also obtain a pseudorandom unitary (PRU) ensemble that is secure under nonadaptive queries. A central feature of both proofs is a new connection between the polynomial method in quantum query complexity and the large-dimension (N) expansion in random matrix theory. In particular, the first construction uses the polynomial method to control high moments of certain random matrix ensembles without requiring delicate Weingarten calculations. In doing so, we define and solve a moment problem on the unit circle, asking whether a finite number of equally weighted points can reproduce a given set of moments. In our second construction, the key step is to exhibit an orthonormal basis for irreducible representations of the partition algebra that has a low-degree large-N expansion. This allows us to show that the distinguishing probability is a low-degree rational polynomial of the dimension N.

algebra↗

Communication-Avoiding and Memory-Constrained Sparse Matrix-Matrix Multiplication at Extreme Scale

Sparse matrix-matrix multiplication (SpGEMM) is a widely used kernel in various graph, scientific computing and machine learning algorithms. In this paper, we consider SpGEMMs performed on hundreds of thousands of processors generating trillions of nonzeros in the output matrix. Distributed SpGEMM at this extreme scale faces two key challenges: (1) high communication cost and (2) inadequate memory to generate the output. Furthermore, we address these challenges with an integrated communication-avoiding and memory-constrained SpGEMM algorithm that scales to 262,144 cores (more than 1 million hardware threads) and can multiply sparse matrices of any size as long as inputs and a fraction of output fit in the aggregated memory. As we go from 16,384 cores to 262,144 cores on a Cray XC40 supercomputer, the new SpGEMM algorithm runs 10x faster when multiplying large-scale protein-similarity matrices.

97 MATHEMATICS AND COMPUTING↗

Atmospheric and Soil Parameters in Five Urban Sites in Knoxville, Tennessee. 2024

This dataset, which contains ten csv files, reports several meteorological conditions such as temperature, wind speed, solar radiation, and net radiation measured in several urban parks in Knoxville, Tennessee, USA: Cumberland Estates Park (CEP), Socially Equal Energy Efficient Development (SEEED), West View Park (WVP), Victor Ashe Park (VAP), and West Hills Park (WHP). Also, hourly average of soil parameters such soil temperature, moisture, and matric potential measured in the five sites are included in this dataset. Air temperature, wind speed and direction, and solar radiation data were obtained from a METER ATMOS 41 All-in-One Weather Station (Pullman, Washington, USA). Soil moisture and temperature parameters were measured the METER Teros 12 probes embedded to a depth of 5 cm into the soil, while soil matric potential was measured with a METER Teros 21 probe. Incident and emitted radiation (shortwave and longwave) measurements were made using an Apogee (Logan, Utah, USA) net radiometer (Model SN-500-SS). This work is a part of a larger study which investigates the impact of soil moisture and plant evapotranspiration on ambient temperature and relative humidity in several city parks in Knoxville, Tennessee.

54 ENVIRONMENTAL SCIENCES↗

Polymer Zwitterions for Stabilization of CsPbBr 3 Perovskite Nanoparticles and Nanocomposite Films

Abstract Functional polymers with sulfobetaine or phosphorylcholine zwitterions as pendent groups are demonstrated as both ligands and host matrices for CsPbBr 3 perovskite nanoparticles (PNPs). These polymers produce nanocomposite films with excellent NP dispersion, optical transparency, and impressive resistance to NP degradation upon exposure to water. Multidentate interactions of the zwitterion‐containing copolymers with the PNPs induce dispersed or weakly aggregated nanocomposite morphologies, depending on the extent of zwitterionic functionality in the polymer. Incorporating additional functionality into the polymers, such as benzophenone pendent groups, yields lithographically patternable films, while time‐resolved photoluminescence measurements provide insight into the electronic impact of PNPs in zwitterionic polymer matrices.

Kim, Hyunki↗

Catalytically Active Site Mapping Realized through Energy Transfer Modeling

Abstract The demands of a sustainable chemical industry are a driving force for the development of heterogeneous catalytic platforms exhibiting facile catalyst recovery, recycling, and resilience to diverse reaction conditions. Homogeneous‐to‐heterogeneous catalyst transitions can be realized through the integration of efficient homogeneous catalysts within porous matrices. Herein, we offer a versatile approach to understanding how guest distribution and evolution impact the catalytic performance of heterogeneous host–guest catalytic platforms by implementing the resonance energy transfer (RET) concept using fluorescent model systems mimicking the steric constraints of targeted catalysts. Using the RET‐based methodology, we mapped condition‐dependent guest (re)distribution within a porous support on the example of modular matrices such as metal–organic frameworks (MOFs). Furthermore, we correlate RET results performed on the model systems with the catalytic performance of two MOF‐encapsulated catalysts used to promote CO 2 hydrogenation and ring‐closing metathesis. Guests are incorporated using aperture‐opening encapsulation, and catalyst redistribution is not observed under practical reaction conditions, showcasing a pathway to advance catalyst recyclability in the case of host–guest platforms. These studies represent the first generalizable approach for mapping the guest distribution in heterogeneous host–guest catalytic systems, providing a foundation for predicting and tailoring the performance of catalysts integrated into various porous supports.

Thompson, William J.↗

Polymer Zwitterions for Stabilization of CsPbBr 3 Perovskite Nanoparticles and Nanocomposite Films

Functional polymers with sulfobetaine or phosphorylcholine zwitterions as pendent groups are demonstrated as both ligands and host matrices for CsPbBr 3 perovskite nanoparticles (PNPs). These polymers produce nanocomposite films with excellent NP dispersion, optical transparency, and impressive resistance to NP degradation upon exposure to water. Multidentate interactions of the zwitterion-containing copolymers with the PNPs induce dispersed or weakly aggregated nanocomposite morphologies, depending on the extent of zwitterionic functionality in the polymer. Incorporating additional functionality into the polymers, such as benzophenone pendent groups, yields lithographically patternable films, while time-resolved photoluminescence measurements provide insight into the electronic impact of PNPs in zwitterionic polymer matrices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytically Active Site Mapping Realized through Energy Transfer Modeling

Abstract The demands of a sustainable chemical industry are a driving force for the development of heterogeneous catalytic platforms exhibiting facile catalyst recovery, recycling, and resilience to diverse reaction conditions. Homogeneous‐to‐heterogeneous catalyst transitions can be realized through the integration of efficient homogeneous catalysts within porous matrices. Herein, we offer a versatile approach to understanding how guest distribution and evolution impact the catalytic performance of heterogeneous host–guest catalytic platforms by implementing the resonance energy transfer (RET) concept using fluorescent model systems mimicking the steric constraints of targeted catalysts. Using the RET‐based methodology, we mapped condition‐dependent guest (re)distribution within a porous support on the example of modular matrices such as metal–organic frameworks (MOFs). Furthermore, we correlate RET results performed on the model systems with the catalytic performance of two MOF‐encapsulated catalysts used to promote CO 2 hydrogenation and ring‐closing metathesis. Guests are incorporated using aperture‐opening encapsulation, and catalyst redistribution is not observed under practical reaction conditions, showcasing a pathway to advance catalyst recyclability in the case of host–guest platforms. These studies represent the first generalizable approach for mapping the guest distribution in heterogeneous host–guest catalytic systems, providing a foundation for predicting and tailoring the performance of catalysts integrated into various porous supports.

Thompson, William J.↗

On‐site surrogates for large‐scale calibration

Abstract Motivated by a computer model calibration problem from the oil and gas industry, involving the design of a honeycomb seal, we develop a new Bayesian methodology to cope with limitations in the canonical apparatus stemming from several factors. We propose a new strategy of on‐site design and surrogate modeling for a computer simulator acting on a high‐dimensional input space that, although relatively speedy, is prone to numerical instabilities, missing data, and nonstationary dynamics. Our aim is to strike a balance between data‐faithful modeling and computational tractability in a calibration framework—tailoring the computer model to a limited field experiment. Situating our on‐site surrogates within the canonical calibration apparatus requires updates to that framework. We describe a novel yet intuitive Bayesian setup that carefully decomposes otherwise prohibitively large matrices by exploiting the sparse blockwise structure. Empirical illustrations demonstrate that this approach performs well on toy data and our motivating honeycomb example.

Huang, Jiangeng↗

A Nonstationary and Non-Gaussian Moving Average Model for Solar Irradiance

Historically, power has flowed from large power plants to customers. Increasing penetration of distributed energy resources such as solar power from rooftop photovoltaic has made the distribution network a two-way-street with power being generated at the customer level. The incorporation of renewables introduces additional uncertainty and variability into the power grid. Distribution network operation studies are being adapted to include renewables; however, such studies require high quality solar irradiance data that adequately reflect realistic meteorological variability. Data from satellite-based products are spatially complete, but temporally coarse, whereas solar irradiances exhibit high frequency variation at very fine timescales. We propose a new stochastic method for temporally downscaling global horizontal irradiance (GHI) to 1 min resolution, but we do not consider the spatial aspect due to limited availability of the in situ irradiance measurements. Solar irradiance's first and second-order structures vary diurnally and seasonally, and our model adapts to such nonstationarity. Empirical irradiance data exhibits highly non-Gaussian behavior; we develop a nonstationary and non-Gaussian moving average model that is shown to capture realistic solar variability at multiple timescales. We also propose a new estimation scheme based on Cholesky factors of empirical autocovariance matrices, bypassing difficult and inaccessible likelihood-based approaches. The model is demonstrated for a case study of three locations that are located in diverse climates through the United States. The model is compared against competitors from the literature and is shown to provide better uncertainty and variability quantification on testing data.

Cholesky factor↗