Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix power”

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 433 records · Page 24

Optimizing Communication in 2D Grid-Based MPI Applications at Exascale

The new reality of exascale computing faces many challenges in achieving optimal performance on large numbers of nodes. A key challenge is the efficient utilization of the message-passing interface (MPI), a critical component for process communication. This paper explores communication optimization strategies to harness the GPU-accelerated architectures of these supercomputers. We focus on MPI applications where processors form a two-dimensional process grid, a common arrangement in applications involving dense matrix operations. This configuration offers a unique opportunity to implement innovative strategies to improve performance and maintain effective load distribution. We study two applications— Dist-FW (Apsp:all-pair-shortest-path) and HPL-MxP (LU factorization with Mixed precision)—on two accelerated systems: Summit (IBM Power and NVIDIA V100) and Frontier (AMD EPYC and MI250X). These supercomputers are operated by the Oak Ridge Leadership Computing Facility (OLCF) and are currently ranked #1 and #5 on the Top500 list. We show how to scale up both applications to exascale levels and tackle the MPI challenges related to implementation, synchronization, and performance. We also compare the performance of several communication strategies at an unprecedented scale. Accurately predicting application performance becomes crucial for cost reduction as the computational scale grows. To address this, we suggest a hyperbolic model as a better alternative to the traditional one-sided asymptotic model for predicting future application performance at such large scales.

Lu, Hao↗

A Dipolar Chiral Spin Liquid on the Breathed Kagome Lattice

Continuous control over lattice geometry, when combined with long-range interactions, offers a powerful yet underexplored tool to generate highly frustrated quantum spin systems. By considering long-range dipolar antiferromagnetic interactions on a breathed Kagome lattice, we demonstrate how these tools can be leveraged to stabilize a chiral spin liquid. We support this prediction with large-scale density-matrix renormalization group calculations and explore the surrounding phase diagram, identifying a route to adiabatic preparation via a locally varying magnetic field. At the same time, we identify the relevant low-energy degrees of freedom in each unit cell, providing a complementary language to study the chiral spin liquid. Finally, we carefully analyze its stability and signatures in finite-sized clusters, proposing direct, experimentally viable measurements of the chiral edge mode in both Rydberg atom and ultracold polar molecule arrays.

Machado↗

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗

Network reconfiguration and distributed energy resource scheduling for improved distribution system resilience

Electric utility companies work to restore as much load as possible after power outages caused by extreme weather events. In this paper, an outage management strategy is proposed to enhance distribution system resilience through network reconfiguration and distributed energy resources (DERs) scheduling. After a line fault, the proposed algorithm can identify radial network topology based on the rank of the incidence matrix. The reconfiguration is implemented by switching tie lines and sectionalizing lines. With the new network topology, an optimal DER scheduling problem is solved to minimize the accumulative cost for dispatchable DER operation and load reduction. Finally, the optimal topology that minimizes the accumulative cost is selected from all radial topologies. The computational workload is relatively low because only linear programming needs to be solved. Using the case studies of the IEEE 69-bus and IEEE 123-bus systems, we consider the worst-case scenarios in which faults occur in the upstream feeder. The simulation results demonstrate that the proposed strategy allows for a relatively high percentage of the load to remain in service after line faults. Furthermore, compared with microgrid-formation approaches, the proposed strategy has advantages when applied to the distribution systems with several normally-open tie lines and low DER penetration.

42 ENGINEERING↗

Holographic Quantum Simulation of Strongly Correlated Electron Systems

The project aimed to demonstrate a new holographic quantum simulation approach and co‐ designed quantum hardware to tackle three specific problems that fall within the broad umbrella of unraveling the physics of strongly correlated electron systems (SCES). These tasks were: (1) holographic preparation of ground‐ and thermal‐ states of correlated magnetic and electronic systems including quasi‐2d frustrated‐spin, Fermi‐Hubbard, and fractional quantum Hall (FQH) systems, (2) holographic‐simulation of long‐time out‐of‐equilibrium dynamics and (3) holographic analogs of embedding methods such as dynamical mean‐ field theory (DMFT) and density‐matrix embedding theory (DMET) to solve systems with complex structure or long‐range interactions. These tasks are prototypes for the kinds of material simulation problems of interest to BES, such as the simulation of multiferroic materials, perovskite photovoltaics and high‐temperature superconductors, that tax the capabilities of the most powerful classical supercomputers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Water Vapor Resistant SiC/SiC Composite for Hydrogen-Based Turbines (Final Scientific Technical Report)

New material innovations are needed for the extreme environments of next generation, hydrogen-fueled turbine engines. Ceramic matrix composites (CMCs) offer high-temperature capability but are susceptible to degradation in high-temperature water vapor. Pratt & Whitney’s (P&W) baseline silicon carbide (SiC)/SiC CMC was initially applied to design a component in a hydrogen-fueled engine, verifying the existence of a design space and extracting boundary conditions to inform testing parameters. Several material innovations, including two fibers, interface coatings (IFC), and self-healing matrices (SHM), were investigated to improve high-temperature performance in water vapor. A boron-doped pyrocarbon (B-PyC) IFC and a SHM matrix with layers of zirconium nitride (ZrN) or zirconium diboride (ZrB 2 ) were developed to fabricate minicomposites with either standard Hi-Nicalon™ Type S (HNS) fibers or new Tyranno® SA4 fibers (SA4). The B-PyC IFC is functional but lacks in providing improved durability. A matrix consisting of thick layers of SiC and thin layers of ZrB 2 shows promise as a SHM that can effectively seal matrix cracks in this extreme environment. Minicomposites with HNS fibers generally outperform those with SA4 fibers.

08 HYDROGEN↗

Discovery of Signatures, Anomalies, and Precursors in Synchrophasor Data with Matrix Profile and Deep Recurrent Neural Networks (Final Project Report)

The widespread deployment of phasor measurement unit (PMU) across the U.S. together with the burgeoning machine learning technology made it possible to develop data-driven PMU data analytics to improve grid security and reliability in a more insightful and effective manner. Although PMU applications have been explored for over a decade, the representative PMU usage is limited to the bulk power system monitoring mainly due to the data integrity issues associated with PMUs (typically missing, fragmented, and wrongly amplified data). To forge a breakthrough on this stalemate and embrace PMUs for power system control and protection as well, we applied various advanced machine learning and big data analysis technology to the power system event detection and classification as the first step toward the power system control and protection pertaining to grid security enhancement.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. Here, we present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_ℓ$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of $C_ℓ$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{IA}$, defined as $S_8D_{IA}=σ_8(Ω_m/0.3)^{0.5}D_{IA}$, where $D_{IA}=1-0.11(A_{IA}-1)$. We find $S_8D_{IA}=0.8954_{-0.039}^{+0.054}$, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Network Based Estimation of Wind Farm Power and Velocity Data Under Changing Wind Direction

This paper describes an estimation algorithm for velocity and power output signals in a wind farm under changing wind direction. A graph-theoretic definition describes the wind farm as a collection of nodes (turbines) and time-varying weighted edges (inter-turbine wake propagation) that change as a function of incoming wind direction. The velocity at each turbine is determined through a discrete input-output model. Changes in wind direction serve as the input and the output is defined in terms of a time-varying weighted adjacency matrix that depends on the time-delay of information propagation between turbines. These delays, which are defined in terms of the advection speed of the wind and the distance between the turbines, capture the delayed effect of wind direction changes on the inter-connectivity of the graph as the wind conditions at the farm inlet propagate through the turbine array. An event-based update framework is employed to capture time-dependent topology changes due to shifts in wind direction. Simulation results for dynamically changing wind inlet directions to a circular wind farm are compared to predictions from both the static and dynamic versions of the FLOw Redirection and Induction in Steady State (FLORIS) model. The approach is shown to enable real-time tracking of dynamic changes to wind farm power output within a framework that can be easily integrated into real-time, horizon-based, control strategies that typically do not account for wind direction changes.

distributed↗

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring Material Solutions for Supercritical CO 2 Applications above 800 °C

There has been recent interest in exploring revolutionary supercritical CO 2 (sCO 2 ) power cycles, and this exploratory investigation was seeking materials with CO 2 compatibility at up to 1200 °C. Here, initial exposures were conducted at 0.1 and 2 MPa CO 2 for up to 1000 h at 900–1200 °C. As expected, specimens of Mo and W that might be used as matrix materials in cermets were rapidly attacked under these conditions. Even an alumina-forming FeCrAlMo alloy showed high mass gains in less than 100 h at 1200 °C due to the formation of Fe-rich oxide. However, at 900–1100 °C, more protective behavior was observed for FeCrAlMo specimens, with or without pre-oxidation, in 0.1 MPa CO 2 , but increased attack was observed in 2 MPa CO 2 . In contrast, most Ni-based alloys exposed at 900–1100 °C showed higher mass gains and thicker reaction products than formed in air. Thus, Ni-based alloys appear less compatible with CO 2 environments above 800 °C compared to lower temperatures. Low mass gains were observed for CVD SiC at 900–1200 °C, but MoSi 2 and Mo(Si,Al) 2 specimens did not form protective scales under these conditions at 1000 and 1100 °C.

36 MATERIALS SCIENCE↗

Pair density wave in the doped three-band Hubbard model on two-leg square cylinders

A pair density wave (PDW) is a superconducting (SC) state with a spatially modulated order parameter. Although much is known about the properties of the PDW state, its realization in microscopic models with divergent susceptibility has been challenging. Here, in this work, we report a density-matrix renormalization group study of a three-band Hubbard model (also known as the Emery model) for cuprates on long two-leg square cylinders. Upon light doping, we find that the ground state of the system is consistent with that of a PDW state with mutually commensurate and power-law SC, charge (CDW), and spin (SDW) density wave correlations. The SC correlations are dominant between neighboring Cu sites with d-wave pairing symmetry. Interestingly, we find that the near-neighbor interactions, especially the near-neighbor attractive V pd interaction between neighboring Cu and oxygen sites, can notably enhance the SC correlations while simultaneously suppressing the CDW correlations. For a modestly strong attractive V pd , the SC correlations become quasi-long-ranged with a divergent PDW susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bridging hadronic and vacuum structure by heavy quarkonia

We discuss the central and, mostly, spin-dependent potentials in heavy quarkonia $\overline{b}$$b$,$\overline{c}$$c$, with two goals in mind. The first is phenomenological: using the splitting between the 1S and 2S pairs, as well as the 1P and 2P quartet masses, we obtain pertinent matrix elements of all these potentials. The second is theoretical: using standard wave functions, we compute these matrix elements, quantifying the perturbative, nonperturbative and “string” contributions. The model for nonperturbative effects is a “dense instanton liquid model,” in which the QCD vacuum is made mostly of (strongly color correlated) instanton-anti-instanton pairs or “molecules.” Their contributions are calculated via standard Wilson lines, with or without extra powers of gauge fields. We find a reasonable description of all central, spin-spin, and spin-orbit forces at distances 𝑟 = 0–0.7 fm, relevant for $\overline{b}$$b$ and $\overline{c}$$c$ quarkonia.

Bound states↗

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cyclic bend over sheave performance of a steel tension member belt for wave energy conversion

Rotary power take-offs (PTOs) promise many benefits for commercial wave energy converters (WECs), including longer power strokes/reduced end-stop concerns, deployment advantages, and controllability. However, these winch-type PTOs can impose millions of bending cycles on the belt or rope element every year, which may lead to a failure mode known as Cyclic Bend Over Sheave, or “CBOS”. Synthetic rope, though attractive due to its high strength, low elongation, and durability, suffers from comparatively poor CBOS performance and struggles to reach a million load cycles without failing. Using belts instead of ropes for winching elements should improve CBOS performance due to having a smaller thickness that undergoes bending. Unfortunately, belts have their own host of problems. Maintaining even tension across the width of a belt subject to twist or fleet angles is a challenge, especially for WECs which can move in all 6 degrees of freedom in the ocean. This study sought to quantify the CBOS performance of a belt constructed of steel tension members in a polyurethane matrix containing longitudinal profiles. Belt samples were tested to failure or up to two million bend-unbend cycles and subject to a range of fleet angles. This paper summarizes the utilized testing methodology and the resulting belt performance.

16 TIDAL AND WAVE POWER↗

Uranium Doped Gallium Nitride Epitaxial Thin Films

Gallium nitride (GaN) is near ubiquitous in modern day technologies, forming the backbone of solid-state lighting and high-power electronics. Engineering the physical properties of GaN has been investigated to some degree by the incorporation or doping of most of the elements of the periodic table, but the actinides remain unexplored. Molecular beam epitaxy is used to demonstrate uranium doping of GaN single crystals. High structural quality of the host matrix is maintained despite partial elemental segregation of the uranium dopant into 1D structures at the levels presented here. Electronic transport measurements reveal relatively high conductivity, which persists down to cryogenic temperature and is characterized by the formation of narrow gaps in the electronic band structures very close to the Fermi level. Photoluminescence measurements reveal that the U-doped GaN exhibits optical behavior similar to that of the GaN substrate. The addition of actinide materials to a non-centrosymmetric, optically active, radiation-hard, and electronically tunable host matrix opens a world of possibilities for investigating and leveraging elements with high electron correlations in the pursuit of novel devices.

36 - MATERIALS SCIENCE↗

Cooperative Research and Development Agreement between National Energy Technology Laboratory and Electric Power Research Institute (EPRI) [Abstract]

Carbon dioxide (CO 2 ) capture from flue gas generated by fossil fuel-fired power plants has been proposed as an efficient approach to limit CO 2 emissions to the atmosphere. Due to the high cyclic CO 2 sorption capacity, well-tuned adsorption chemistry and non-volatility, solid sorbents are widely studied for CO 2 capture processes such as pressure swing adsorption (PSA) and temperature swing adsorption (TSA). Realistic applications of traditional solid sorbent systems face many challenges since moisture and heat management are problematic and solid-solid heat exchange is inefficient. Recently, a novel solid sorbent system, the sorbent polymer composite (SPC), has been developed to overcome those challenges in an energy-saving CO 2 capture process using TSA membrane contactors. A SPC material is comprised of a powdered solid sorbent embedded into a hydrophobic and porous polymer matrix. It allows gases to permeate through and achieve full contact with the sorbents while rejecting water moisture in the CO 2 capture process. In this collaborative research, NETL will conduct performance testing of a novel SPC material and EPRI will conduct a cost analysis based on that performance data.

20 FOSSIL-FUELED POWER PLANTS↗