Results for “Dynamical mode decomposition” — Page 3

DOE OSTI2024

Weighted Composition Operators for Learning Nonlinear Dynamics

Operator theoretic methods in dynamical system have been dominated by the use of Koopman operators and their continuous time counterparts, such as Koopman Generators and Liouville Operators. The advantage gained from their use primarily stems from the ability to extract subspaces and eigenfunctions within a space of observables that are invariant with respect to the Koopman operator over that space. When this occurs, a dynamic mode decomposition of the systems state provides a linear model for the dynamical system. Not all Koopman operators have eigenfunctions that may be exploited in this manner. However, the framework can still be leveraged for approximations using other operators. In this setting, we present a different operator for the study of dynamical systems, the weighted composition operator. These operators are compact for a wide range of dynamics and spaces, and through their interactions with occupation kernels and vector valued kernels, they admit an estimation of the underlying dynamics. Here, this manuscript presents a new algorithm for the data driven study of dynamical systems from data, and also provides two numerical experiments where convergence is achieved as a proof of concept.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.

Aeroacoustics↗

DOE OSTI2024

Code for the manuscript "Mori-Zwanzig Modal Decomposition"

We would like to create an open source repository in LANL's github on code written in Julia, in which we implement and extend the data-driven Mori-Zwanzig method for extracting large-scale spatio-temporal structures from data, which we call MZMD. This method is an extension of Dynamic Mode Decomposition (DMD) in which Mori-Zwanzig memory kernels are included into the associated companion matrix. In the code we would like to release, we apply MZMD to a flow over a cylinder with Reynolds number 100 rather than the much larger data set used in the associated manuscript. DMD is used extensively in the fluid dynamics community mainly for extracting large scale spatio-temporal structures (patters) from flow data. This is useful for understanding the key mechanisms that generate certain complex dynamical process relevant in engineering design. In MZMD, we improve upon DMD by adding the Mori-Zwanzig memory kernels, and show this improvement is especially important in strongly nonlinear regions of the flow.

Woodward, Michael↗

NASA NTRS2019

Investigation of Transitional Shock-Wave/Boundary Layer Interactions Using Direct Numerical Simulations

Interaction of a transitional boundary layer with a normal-shock is investigated using unstructured tetrahedral meshes under the numerical framework of the space-time conservation element, solution element (CESE) method. The computations mimic recent experimental efforts at the University of Tennessee Space Institute, where a Mach 2.0 flow interacts with a tall cylinder attached to a flat plate. The location of the cylinder with respect to the flat plate leading edge determines if the incoming boundary layer is laminar, transitional or fully turbulent. Four representative flow conditions exemplifying laminar and transitional boundary layers are analyzed by direct numerical simulations. Similar to what was observed in the experiments for the case of transitional interaction, the computations reveal an intermittent upstream influence (UI) shock that repeatedly travels upstream from the lambda-foot toward the leading edge before vanishing. Through detailed unsteady flow analysis obtained using Fourier analysis and dynamic mode decomposition techniques, the presence of disturbances with similar frequencies as those measured in experiments were identified in the flow along with locations that appear to influence the dynamics of the flow.

Venkatachari, Balaji Shankar↗

DOE OSTI2022

Design of a supervisory control system for autonomous operation of advanced reactors

Advanced reactors to be deployed in the coming decades will face deregulated energy markets, and may adopt flexible operation to boost profitability. To aid in the transition from baseload to flexible operation paradigm, autonomous operation is sought. This work focuses on the control aspect of autonomous operation. Specifically, a hierarchical control system is designed to support constraint enforcement during routine operational transients. Within the system, data-driven modeling, physics-based state observation, and classical control algorithms are integrated to provide an adaptable and robust solution. A 320 MW Fluoride-cooled High-temperature Pebble-bed Reactor is the design basis for demonstrating the proposed control system. The hierarchical control system consists of a supervisory layer and low-level layer. The supervisory layer receives requests to change the system's operating conditions (e.g., the current reactor power to meet a load -follow), and accepts or rejects them based on constraints that have been assigned. Constraints are issued to keep the plant within an optimal operating region. The low-level layer interfaces with the actuators of the system to fulfill requested changes, while maintaining tracking and regulation duties. Further, to accept requests at the supervisory layer, the Reference Governor algorithm was adopted. To model the dynamics of the reactor, a system identification algorithm, Dynamic Mode Decomposition, was utilized. To estimate the evolution of process variables that cannot be directly measured (e.g., the propagation of delayed neutron precursors), the Unscented Kalman Filter, incorporating a nonlinear model of nuclear dynamics, was adopted. The composition of these algorithms led to a numerical demonstration of constraint enforcement during a 40% power drop transient (at a rate of 5 %/min). Uncontrolled secondary-side temperatures were successfully constrained. Adaptability of the proposed system was demonstrated by modifying the constraint values, and enforcing them during the transient. Robustness was also demonstrated by enforcing constraints under noisy environments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Towards Improving luminosity using optics tuning and data-driven methods

The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.

Accelerator Physics↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

DOE OSTI2023

DRIPS: A framework for dimension reduction and interpolation in parameter space

Reduced-order models are often used to describe the behavior of complex systems, whose simulation with a full model is too expensive, or to extract salient features from the full model’s output. We introduce a new model-reduction framework DRIPS (dimension reduction and interpolation in parameter space) that combines the offline local model reduction with the online parameter interpolation of reduced-order bases (ROBs). The offline step of this framework relies on dynamic mode decomposition (DMD) to build a low-rank linear surrogate model, equipped with a local ROB, for quantities of interest derived from the training data generated by repeatedly solving the (nonlinear) high-fidelity model for multiple parameter points. The online step consists of the construction of a parametric reduced-order model for each target/test point in the parameter space, with the interpolation of ROBs done on a Grassman manifold and the interpolation of reduced-order operators done on a matrix manifold. The DMD component enables DRIPS to model (typically low-dimensional) quantities of interest directly, without having to access the (typically high-dimensional and possibly nonlinear) operators in a high-fidelity model that governs the dynamics of the underlying high-dimensional state variables, as required in projection-based reduced-order modeling. A series of numerical experiments suggests that DRIPS yields a model reduction, which is computationally more efficient than the commonly used projection-based proper orthogonal decomposition; it does so without requiring a prior knowledge of the governing equation for quantities of interest. Furthermore, for the nonlinear systems considered, DRIPS is more accurate than Gaussian-process interpolation (Kriging).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

Implementation of an Unsteady PSP System in the NASA TDT

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment. Note: this presentation is an MP4 video with sound, color with a run time of 10 minutes 37 seconds.

unsteady pressure-sensitive paint↗

DOE OSTI2025

Harnessing on-machine metrology data for prints with a surrogate model for laser powder directed energy deposition

In this study, we leverage the massive amount of multi-modal on-machine metrology data generated from Laser Powder Directed Energy Deposition (LP-DED) to construct a comprehensive surrogate model of the 3D printing process. By employing Dynamic Mode Decomposition with Control (DMDc), a data-driven technique, we capture the complex physics inherent in this extensive dataset. This physics-based surrogate model emphasizes thermodynamically significant quantities, enabling us to accurately predict key process outcomes. The model ingests 21 process parameters, including laser power, scan rate, and position, while providing outputs such as melt pool temperature, melt pool size, and other essential observables. Furthermore, it incorporates uncertainty quantification to provide bounds on these predictions, enhancing reliability and confidence in the results. We then deploy the surrogate model on a new, unseen part and monitor the printing process as validation of the method. Our experimental results demonstrate that the predictions align with actual measurements with high accuracy, confirming the effectiveness of our approach. Furthermore, this methodology not only facilitates real-time predictions but also operates at process-relevant speeds, establishing a basis for implementing feedback control in LP-DED.

Digital twins↗

DOE OSTI2024

Data-driven models of nonautonomous systems

Nonautonomous dynamical systems are characterized by time-dependent inputs, which complicates the discovery of predictive models describing the spatiotemporal evolution of the state variables of quantities of interest from their temporal snapshots. When dynamic mode decomposition (DMD) is used to infer a linear model, this difficulty manifests itself in the need to approximate the time-dependent Koopman operators. Our approach is to approximate the original nonautonomous system with a modified system derived via a local parameterization of the time-dependent inputs. The modified system comprises a sequence of local parametric systems, which are subsequently approximated by a parametric surrogate model using the DRIPS (dimension reduction and interpolation in parameter space) framework. The offline step of DRIPS relies on DMD to build a linear surrogate model, endowed with reduced-order bases for the observables mapped from training data. The online step interpolates on suitable manifolds to construct a sequence of iterative parametric surrogate models; the target/test parameter points on these manifolds are specified by a local parameterization of the test time-dependent inputs. Here, we use numerical experimentation to demonstrate the robustness of our method and compare its performance with that of deep neural networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Online real-time learning of dynamical systems from noisy streaming data

Abstract Recent advancements in sensing and communication facilitate obtaining high-frequency real-time data from various physical systems like power networks, climate systems, biological networks, etc. However, since the data are recorded by physical sensors, it is natural that the obtained data is corrupted by measurement noise. In this paper, we present a novel algorithm for online real-time learning of dynamical systems from noisy time-series data, which employs the Robust Koopman operator framework to mitigate the effect of measurement noise. The proposed algorithm has three main advantages: (a) it allows for online real-time monitoring of a dynamical system; (b) it obtains a linear representation of the underlying dynamical system, thus enabling the user to use linear systems theory for analysis and control of the system; (c) it is computationally fast and less intensive than the popular extended dynamic mode decomposition (EDMD) algorithm. We illustrate the efficiency of the proposed algorithm by applying it to identify the Van der Pol oscillator, the chaotic attractor of the Henon map, the IEEE 68 bus system, and a ring network of Van der Pol oscillators.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2019

A Qualitative Investigation of Selected Infrared Flow Visualization Image Processing Techniques

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel testenvironment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is discussed as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity.This method makes use of dynamic mode decomposition and discrete-time filtering to separate the background reflections that dominate low emissivity surfaces from the aerodynamic driven surface temperature gradients.This process will be shown to resolve the onset of boundary layer transition on a bare metal wing as well as identify and resolve hidden features in the image data. While the implementation of this technique is very preliminary it demonstrates the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Image Processeing↗

NASA NTRSDate not supplied

Unsteady PSP in the NASA Transonic Dynamics Tunnel

For the first time, unsteady pressure sensitive paint (uPSP) has been applied in the NASA Langley Transonic Dynamics Tunnel. Obtaining global surface pressure measurements using the uPSP technique required the development of a new paint formulation for use in the low oxygen heavy gas atmosphere, as well as environmental enclosures to protect sensitive electro-optical components from the high temperature, low pressure environment present during tunnel operation. A high-speed datalink connecting the wind tunnel to Langley’s local high performance compute resource was also established for this test to enable near real time processing of the large datasets that were obtained throughout the campaign. A high-speed lifetime measurement technique was also utilized to yield steady state surface pressures at each condition using the same equipment that was used to provide unsteady measurements. Important metrics such as pressure time histories and power spectral density are compared against traditional unsteady pressure point measurements, and more advanced data products such as dynamic mode decomposition are also explored to provide insight into the underlying flow phenomena.

Daniel T. Reese↗

DOE OSTI2023

The Construction of Reduced Order Models for the HYBRID Repository

Reduced order models (ROMs) present an opportunity to reduce the computational burden of IES analysis using HYBRID models without significantly sacrificing data accuracy. Other applications of ROMs are digital twins, stochastic optimization studies, model exporting, and protection of proprietary information. In this paper the capability of creating ROMs of the models within the HYBRID repository will be demonstrated. Specifically, a linear regression ROM of a gas turbine and a dynamic mode decomposition with control (DMDc) ROM of a lithium-ion battery are constructed.

97 MATHEMATICS AND COMPUTING↗