Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Koopmanism”

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 73 records · Page 4

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗

Nonunique fraction of Fock exchange for defects in two-dimensional materials

By investigating the vacancy and substitutional defects in monolayer WS 2 with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band-gap and band-edge positions. Fixing the mixing parameter of Fock exchange based on the band gap can lead to qualitatively incorrect defect physics in two-dimensional materials. Instead, excellent agreement is achieved with both experiment and many-body perturbation theory within the GW approximation once the mixing parameters are tuned individually for the defect species and the band edges. Here, we show the departure from a unique optimized mixing parameter is inherent to two-dimensional systems as the band edges experience a reduced screening while the localized defects are subject to bulklike screening.

2-dimensional systems↗

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗

pnnl/grid_prediction

Two datadriven predictive approaches, namely, {\em Koopman Operator Theoretic (KOT)-based} model and {\em Graph Neural Network (GNN)-based model}, to enable effective power system state predictions. The KOT-based approaches (Robust DMD, deepDMD) capture the power system evolution as a linear dynamical system on an abstract space. The GNNs model the spatio-temporal correlations using graph convolutional network and are called Spatio-Temporal Graph Convolutional Network (STGCN). These predictive models are trained, tested and compared rigorously based on their predictions of frequencies in the IEEE 68 bus system when subjected to a disturbance. GridSTAGE framework developed at Pacific Northwest National Laboratory is leveraged to generate multiple datasets (in the form of PMU measurements) for training and testing by strategically creating load changes across the spatial locations of the network

Nandanoor, Sai Pushpak↗

Simple self-consistent method for excited states in density functional theory to characterize defect-derived behavior in wide-band-gap-based microelectronic materials

This final report summarizes the results of the Laboratory Direct Research and Development (LDRD) Project Number 229740. Wide band gap semiconductors such as gallium nitride (GaN) have features highly desirable for multiple mission electronic applications. Realization of their potential requires atomic-scale understanding of electronic behavior. The principal experimental tools for electronically probing defects in GaN are chemically undifferentiating and lack a practical theoretical counterpart needed to identify and characterize specific defects. This project investigated whether a simple idea for modeling defect excited states and their associated photoluminescence (PL) energies is viable, as a path to accelerate the understanding of defect behavior and gain valuable insights into engineering new electronic materials and devices. The research implemented a non-self-consistent total-energy evaluation of a Koopmans-type estimation of an excited electronic state energy in density functional theory (DFT) calculations, and proceeded to design, implement, and assess a self-consistent method for computing excited states based upon an OCcupation-Constrained-DFT (occ-DFT). The occ-DFT was verified in test calculations of defect excited states and validated against well-characterized PL data for 3d transition metal defects in GaN. The method proved stable and robust in computing excited states and gave accurate predictions compared to experimental PL data. The combined ground state/excited-state capability proved capable of chemically differentiating defect species in GaN. In application to 3d dopants in GaN, we reinterpreted extensive experimental literature, proposed new defects as prospective candidates for use in quantum information applications, and outlined design strategies to create and exploit these potentially useful functional defects in GaN.

36 MATERIALS SCIENCE↗

Low Purge and Mercury Recovery Testing with Sludge Batch 10 Tank 40 Simulant

Researchers at the Savannah River National Laboratory (SRNL) have completed testing as requested by Savannah River Mission Completion (SRMC) to perform experiments to determine the impact of using a lower or inert purge in Sludge Batch (SB) 10 processing under the Nitric-Glycolic Acid (NGA) flowsheet. A key objective of this testing was also to determine the mercury speciation and recovery during each experiment. The testing was performed as part of the SB10 Technical Task Request (TTR) and Task Technical and Quality Assurance Plan (TTQAP). Two sets of tests were performed, and a Run Plan was approved prior to each set of experiments to document the planned testing. Three initial experiments were completed to determine whether a low air purge would be beneficial to CPC processing at higher acid stoichiometry (110%) based on the Koopman minimum acid (KMA) equation (116% Hsu). One of the tests, an inert nitrogen purge experiment, was also completed to demonstrate that excluding oxygen did not introduce any new hazards. The experiments were designed to be identical except for the change in purge gas and purge flowrate from run to run. After reviewing the results from the initial three experiments, six additional tests were proposed by SRNL to support the lower purge study and to look for processing alternatives for improving mercury recovery. These additional six tests were all completed at a very low acid stoichiometry to mimic the pH experienced during processing in DWPF (~7). DWPF is processing SB10 sludge at an acid stoichiometry of 90% based on the Hsu equation. The additional SRNL experiments were performed at an acid stoichiometry of 62.5% KMA stoichiometry (66.3% Hsu) to produce a Sludge Receipt and Adjustment Tank (SRAT) product with a pH of about 7. All experiments used simulants of both SWPF streams, although no entrained solvent was added during any of the experiments

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data Analysis with Graphical Models: Software Tools

Probabilistic graphical models (directed and undirected Markov fields, and combined in chain graphs) are used widely in expert systems, image processing and other areas as a framework for representing and reasoning with probabilities. They come with corresponding algorithms for performing probabilistic inference. This paper discusses an extension to these models by Spiegelhalter and Gilks, plates, used to graphically model the notion of a sample. This offers a graphical specification language for representing data analysis problems. When combined with general methods for statistical inference, this also offers a unifying framework for prototyping and/or generating data analysis algorithms from graphical specifications. This paper outlines the framework and then presents some basic tools for the task: a graphical version of the Pitman-Koopman Theorem for the exponential family, problem decomposition, and the calculation of exact Bayes factors. Other tools already developed, such as automatic differentiation, Gibbs sampling, and use of the EM algorithm, make this a broad basis for the generation of data analysis software.

Buntine, Wray L.↗

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Grand Challenges and Opportunities in Stimulated Dynamic and Resonant Catalysis

Traditional heterogeneous catalysis is constrained by kinetic and thermodynamic limits, such as the Sabatier principle and reaction equilibrium. Dynamic and resonant catalysts hold promise to overcome these limitations by actively oscillating a catalyst’s physical or electronic structure at the time scale of the catalytic cycle, allowing programmable control over reaction pathways, and leading to improved rate and selectivity. External stimuli such as temperature swing, mechanical strain, electric charge, and light can perturb catalyst surfaces in different ways, altering adsorbate coverage, binding energies, and transition states beyond what steady-state catalysis allows. This work surveys the current state of dynamic catalysis, introduces the concept of “ stimulando ” characterization for observing transient dynamics, and outlines key modeling, mechanistic, and benchmarking strategies to advance the field toward improved chemical transformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A compendium of human gene functions derived from evolutionary modelling

A comprehensive, computable representation of the functional repertoire of all macromolecules encoded within the human genome is a foundational resource for biology and biomedical research. The Gene Ontology Consortium has been working towards this goal by generating a structured body of information about gene functions, which now includes experimental findings reported in more than 175,000 publications for human genes and genes in experimentally tractable model organisms 1,2 . Here, we describe the results of a large, international effort to integrate all of these findings to create a representation of human gene functions that is as complete and accurate as possible. Specifically, we apply an expert-curated, explicit evolutionary modelling approach to all human protein-coding genes. This approach integrates available experimental information across families of related genes into models that reconstruct the gain and loss of functional characteristics over evolutionary time. The models and the resulting set of 68,667 integrated gene functions cover approximately 82% of human protein-coding genes. The functional repertoire reveals a marked preponderance of molecular regulatory functions, and the models provide insights into the evolutionary origins of human gene functions. We show that our set of descriptions of functions can improve the widely used genomic technique of Gene Ontology enrichment analysis. The experimental evidence for each functional characteristic is recorded, thereby enabling the scientific community to help review and improve the resource, which we have made publicly available.

59 BASIC BIOLOGICAL SCIENCES↗

Atacama Cosmology Telescope measurements of a large sample of candidates from the Massive and Distant Clusters of WISE Survey

Galaxy clusters are an important tool for cosmology, and their detection and characterization are key goals for current and future surveys. Using data from the Wide-field Infrared Survey Explorer (WISE), the Massive and Distant Clusters of WISE Survey (MaDCoWS) located 2839 significant galaxy overdensities at redshifts 0.7 ≲ z ≲ 1.5, which included extensive follow-up imaging from the Spitzer Space Telescope to determine cluster richnesses. Concurrently, the Atacama Cosmology Telescope (ACT) has produced large area millimeter-wave maps in three frequency bands along with a large catalog of Sunyaev-Zeldovich (SZ)-selected clusters as part of its Data Release 5 (DR5). Aims. We aim to verify and characterize MaDCoWS clusters using measurements of, or limits on, their thermal SZ effect signatures. Here, we also use these detections to establish the scaling relation between SZ mass and the MaDCoWS-defined richness.

79 ASTRONOMY AND ASTROPHYSICS↗

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗

State-resolved ultrafast charge and spin dynamics in [Co/Pd] multilayers

We use transient absorption spectroscopy with circularly polarized x rays to detect laser-excited hole states below the Fermi level and compare their dynamics with that of unoccupied states above the Fermi level in ferromagnetic [Co/Pd] multilayers. While below the Fermi level, an instantaneous and significantly stronger demagnetization is observed, above the Fermi level, the demagnetization is delayed by 35 ± 10 fs. This provides a direct visualization of how ultrafast demagnetization proceeds via initial spin-flip scattering of laser-excited holes to the subsequent formation of spin waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Simons Observatory: Cryogenic half wave plate rotation mechanism for the small aperture telescopes

We present the requirements, design, and evaluation of the cryogenic continuously rotating half-wave plate (CHWP) for the Simons Observatory (SO). SO is a cosmic microwave background polarization experiment at Parque Astronómico de Atacama in northern Chile that covers a wide range of angular scales using both small (Ø0.42 m) and large (Ø6 m) aperture telescopes. In particular, the small aperture telescopes (SATs) focus on large angular scales for primordial B-mode polarization. To this end, the SATs employ a CHWP to modulate the polarization of the incident light at 8 Hz, suppressing atmospheric 1/f noise and mitigating systematic uncertainties that would otherwise arise due to the differential response of detectors sensitive to orthogonal polarizations. The CHWP consists of a 505 mm diameter achromatic sapphire HWP and a cryogenic rotation mechanism, both of which are cooled down to ~50 K to reduce detector thermal loading. Under normal operation, the HWP is suspended by a superconducting magnetic bearing and rotates with a constant 2 Hz frequency, controlled by an electromagnetic synchronous motor. We find that the number of superconductors and the number of magnets that make up the superconducting magnetic bearing are important design parameters, especially for the rotation mechanism’s vibration performance. The rotation angle is detected through an angular encoder with a noise level of 0.07 μrad √ $\overline{s}$. During a cooldown process, the rotor is held in place by a grip-and-release mechanism that serves as both an alignment device and a thermal path. In this paper, we provide an overview of the SO SAT CHWP: its requirements, hardware design, and laboratory performance.

47 OTHER INSTRUMENTATION↗