Results for “Matrix models” — Page 68

DOE OSTI2021

Subleading corrections to the S 3 free energy of necklace quiver theories dual to massive IIA

We investigate the S 3 free energy of N = 3 Chern-Simons-matter quiver gauge theories with gauge group U(N) r (r ≥ 2) where the sum of Chern-Simons levels does not vanish, beyond the leading order in the large-N limit. We take two different approaches to explore the sub-leading structures of the free energy. First we evaluate the matrix integral for the partition function in the ’t Hooft limit using a saddle point approximation. Second we use an ideal Fermi-gas model to compute the same partition function, but in the limit of fixed Chern-Simons levels. The resulting expressions for the free energy F = – log Z are then compared in the overlapping parameter regime. The Fermi-gas approach also hints at a universal 1/6 log N correction to the free energy. Since the quiver gauge theories we consider are dual to massive Type IIA theory, we expect the sub-leading correction of the planar free energy in the large ’t Hooft parameter limit to match higher-derivative corrections to the tree-level holographic dual free energy, which have not yet been fully investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗

DOE OSTI2022

Natural orbitals for the ab initio no-core configuration interaction approach

Ab initio no-core configuration interaction (NCCI) calculations for the nuclear many-body problem have traditionally relied upon an antisymmetrized product (Slater determinant) basis built from harmonic oscillator orbitals. The accuracy of such calculations is limited by the finite dimensions which are computationally feasible for the truncated many-body space. We therefore seek to improve the accuracy obtained for a given basis size by optimizing the choice of single-particle orbitals. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which maximizes the occupation of low-lying orbitals, thus accelerating convergence in a configuration-interaction basis, while also possibly providing physical insight into the single-particle structure of the many-body wave function. We describe the implementation of natural orbitals in the NCCI framework and examine the nature of the natural orbitals thus obtained, the properties of the resulting many-body wave functions, and the convergence of observables. After taking 3 He as an illustrative testbed, we explore aspects of NCCI calculations with natural orbitals for the ground state of the p-shell neutron halo nucleus 6 He .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

PAO 1.0: A Python Library for Adversarial Optimization

PAO is a Python-based package for Adversarial Optimization. The goal of this package is to provide a general modeling and analysis capability for bilevel, trilevel and other multilevel optimization forms that express adversarial dynamics. PAO integrates two different modeling abstractions: 1. Algebraic models extend the modeling concepts in the Pyomo algebraic modeling language to express problems with an intuitive algebraic syntax. Thus, we expect that this modeling abstraction will commonly be used by PAO end-users. 2. Compact models express objective and constraints in a manner that is typically used to express the mathematical form of these problems (e.g. using vector and matrix data types). PAO denes custom Multilevel Problem Representations (MPRs) that simplify the implementation of solvers for bilevel, trilevel and other multilevel optimization problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Enhanced Hanford Low-Activity Waste Glass Property Data Development: Phase 2

This report summarizes and analyzes the data collected on a second test matrix of 42 low-activity waste (LAW) glass compositions intended to expand the LAW glass composition region over which glass property-composition models are valid. The 42 LAW glass compositions were selected using a statistical layered experimental design approach to explore a new glass composition region that overlaps and expands beyond the previously explored LAW glass composition region. The layered design of 42 LAW glasses consisted of 21 outer-layer glasses, 15 inner-layer glasses, a center glass (three replicates), and two glasses previously tested at another laboratory. The outer-layer glass compositions include extrema based on current projections of the LAW feeds and melter operating temperatures. One of the outer-layer compositions was eliminated due to the inability to form a glass. Therefore, only 41 glasses were analyzed. The analyses performed on these glasses include chemical composition (for target compositional verification), density, viscosity, electrical conductivity, crystal fraction, canister centerline cooling with crystal identification, the Product Consistency Test, sulfur solubility, and the Vapor Hydration Test. This report discusses the results obtained from this testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Using optimal estimation to retrieve winds from velocity-azimuth display (VAD) scans by a Doppler lidar

Abstract. Low-powered commercially available coherent Doppler lidar (CDL) wind profilers provide continuous measurement of vertical profiles of wind in the lower troposphere, usually close to or up to the top of the planetary boundary layer. The vertical extent of these wind profiles is limited by the availability of scatterers and thus varies substantially throughout the day and from one day to the next. This makes it challenging to develop continuous products that rely on CDL-observed wind profiles. In order to overcome this problem, we have developed a new method for wind profile retrievals from CDL that combines the traditional velocity-azimuth display (VAD) technique with optimal estimation (OE) to provide continuous wind profiles up to 3 km. The new method exploits the level-to-level covariance present in the wind profile to fill in the gaps where the signal-to-noise ratio of the CDL return is too low to provide reliable results using the traditional VAD method. Another advantage of the new method is that it provides the full error covariance matrix of the solution and profiles of information content, which more easily facilitates the assimilation of the observed wind profiles into numerical weather prediction models. This method was tested using yearlong CDL measurements at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Central Facility in 2019. Comparison with the ARM operational CDL wind profile product and collocated radiosonde wind measurements shows excellent agreement (R2>0.99) with no degradation in results where the traditional VAD provided a valid solution. In the region where traditional VAD does not provide results, the OE wind speed and wind vector have uncertainties of 3.44 and 4.33 m s−1, respectively. As a result, the new method provides additional information over the standard technique and increases the effective range of existing CDL systems without the need for additional hardware.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Enhanced Diffusive Transport in Fluctuating Porous Media

Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

ATOMIC Simulations and Experimental Data for CaCO3 Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1] and are used to develop a statistical characterization of matrix effects. Our overall suite of simulations includes contains several sets of simulations: training and validation sets of simulations for three and four element mixtures of calcium, carbon, oxygen, and nitrogen (included to account for atmosphere) along with simulations of the individual elements. The 4-element simulations include the mixture of all four elements mentioned and for each of the four individual elements. The 3-element simulations include output for the mixture of calcium, carbon, oxygen and for these three individual elements. The training data were produced using a 600-run design, shown in Figure 1, that varies input parameters temperature (T), electron density (Ne), and proportion of the elements calcium, carbon, oxygen, and nitrogen (Ca; C; O; N) for the 4 element output. The 3-element output includes all parameters except for the proportion of nitrogen. The element proportions (all the variables but T and Ne) sum to one and are unused in the single-element simulations. The validation data was produced with a 80-run design shown in Figure 2. The training and validation simulation outputs for the 4-element simulations for the mixture and for the single element calcium are shown as sample simulations in Figures 3 and 4 respectively. The simulations produce spectra over a range of 190nm - 950nm that roughly mimics the range collected by the SciAps Z-300 LIBS instrument that was used for the experimental data. The measured spectra for a CaCO3 (which may include contribution from Earth's atmosphere) in the experiment is shown in in Figure 5. All files are kept in directories whose names indicate the elemental composition (CaCO3, Ca, C, O, or N), number of elements (3 or 4), and purpose (training, which is not labeled in the file name, or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. [1] J Colgan, EJ Judge, DP Kilcrease, and JE Barefield II. Ab-initio modeling of an iron laser-induced plasma: Comparison between theoretical and experimental atomic emission spectra. Spectrochimica Acta Part B: Atomic Spectroscopy, 97:65{73}, 2014.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Case Study of Using Kokkos and SYCLs Performance-Portable Frameworks for Milc-Dslash Benchmark on NVIDIA, AMD and Intel GPUs

Six of the top ten supercomputers in the TOP500 list from June 2021 rely on NVIDIA GPUs to achieve their peak compute bandwidth. With the announcement of Aurora, Frontier, and El Capitan, Intel and AMD have also entered the domain of providing GPUs for scientific computing. A consequence of the increased diversity in the GPU landscape is the emergence of portable programming models such as Kokkos, SYCL, OpenCL, and OpenMP, which allow application developers to maintain a single-source code across a diverse range of hardware architectures. While the portable frameworks try to optimize the compute resource usage on a given architecture, it is the programmers responsibility to expose parallelism in an application that can take advantage of thousands of processing elements available on GPUs. In this paper, we introduce a GPU-friendly parallel implementation of Milc-Dslash that exposes multiple hierarchies of parallelism in the algorithm. Milc-Dslash was designed to serve as a benchmark with highly optimized matrix-vector multiplications to measure the resource utilization on the GPU systems. The parallel hierarchies in the Milc-Dslash algorithm are mapped onto a target hardware using Kokkos and SYCL programming models. We present the performance achieved by Kokkos and SYCL implementations of Milc-Dslash on NVIDIA A100 GPU, AMD MI100 GPU, and Intel Gen9 GPU. Additionally, we compare the Kokkos and SYCL performances with those obtained from the versions written in CUDA and HIP programming models on NVIDIA A100 GPU and AMD MI100 GPU, respectively.

Dufek, Amanda S↗

DOE OSTI2021

Excitonic wave-packet evolution in a two-orbital Hubbard model chain: A real-time real-space study

Motivated by experimental developments introducing the concept of spin-orbit separation, we study the real-space real-time evolution of an excitonic wave packet using a two-orbital Hubbard model in a chain. The exciton is created by exciting an electron from a lower-energy half-filled orbital to a higher-energy empty orbital. Here we carry out the real-time dynamics of the resulting excitonic wave packet using the time-dependent density matrix renormalization group. We find clear evidence of charge-spin and spin-orbit separation in real space, by tracking the time evolution of local observables. We show that the velocity of the orbiton can be tuned varying the interorbital interactions. We also present a comparative study of the dynamics of a hole in one-orbital and two-orbital Hubbard models. Moreover, we analyze the dynamics of an exciton with spin-flip excitation, where we observe fractionalized spinons induced by Hund's interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2020

Optimal D-FACTS Placement in Moving Target Defense Against False Data Injection Attacks

Moving target defense (MTD) is a defense strategy to detect stealthy false data injection (FDI) attacks against the power system state estimation using distributed flexible AC transmission system (D-FACTS) devices. However, existing studies neglect to address a fundamental yet critical issue, i.e., the D-FACTS placement, by assuming that all lines are equipped with D-FACTS devices. Here, to tackle this problem, we first derive analytical necessary conditions and requirements on the D-FACTS placement for a complete MTD. Further, we propose sufficient conditions using a graph theory-based topology analysis to ensure that the MTD under the proposed D-FACTS placement has the maximum rank of its composite matrix, which is indicative of the MTD effectiveness. Based on the analytical conditions, we design D-FACTS placement algorithms by using the minimum number of D-FACTS devices to achieve the maximum MTD effectiveness. A novel MTD-based ACOPF model, in which the reactance of D-FACTS lines is introduced as decision variables, is proposed to find a trade-off between the system loss and the MTD effectiveness. Numerical results on 6-bus, IEEE 14-bus, and IEEE 118-bus systems show the efficacy of MTDs using the proposed D-FACTS placement algorithms in maximizing the composite matrix rank and detecting FDI attacks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

A General Materials Data Science Framework for Quantitative 2D Analysis of Particle Growth from Image Sequences

Abstract Phase transformations are a challenging problem in materials science, which lead to changes in properties and may impact performance of material systems in various applications. We introduce a general framework for the analysis of particle growth kinetics by utilizing concepts from machine learning and graph theory. As a model system, we use image sequences of atomic force microscopy showing the crystallization of an amorphous fluoroelastomer film. To identify crystalline particles in an amorphous matrix and track the temporal evolution of the particle dispersion, we have developed quantitative methods of 2D analysis. 700 image sequences were analyzed using a neural network architecture, achieving 0.97 pixel-wise classification accuracy as a measure of the correctly classified pixels. The growth kinetics of isolated and impinged particles were tracked throughout time using these image sequences. The relationship between image sequences and spatiotemporal graph representations was explored to identify the proximity of crystallites from each other. The framework enables the analysis of all image sequences without the requirement of sampling for specific particles or timesteps for various materials systems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

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 Ljungströ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↗