Results for “Coherent Mode Decomposition”

DOE OSTI2022

Memory and CPU efficient coherent mode decomposition of partially coherent synchrotron radiation with subtraction of common quadratic phase terms

Application examples of a memory and CPU efficient coherent mode decomposition (CMD) method for wave-optics based simulation of the partially coherent undulator radiation propagation through a hard X-ray beamline in a 3rd generation synchrotron radiation source are presented. The high efficiency of the method is achieved thanks to the analytical treatment of the common quadratic phase terms that are developed in the phase of cross-spectral density (CSD) of partially coherent radiation at a distance from source. This treatment allows for a considerable, several orders of magnitude, reduction of the 4D CSD mesh density (and the memory occupied by the CSD) required for ensuring sufficient accuracies of wavefront propagation simulations with the modes produced by the CMD at a beamline entrance. This method, implemented in the “Synchrotron Radiation Workshop” open-source software, dramatically increases the feasibility of the CMD of 4D CSD for producing 2D coherent modes for a large variety of applications at storage rings and other types of radiation sources.

36 MATERIALS SCIENCE↗

DOE OSTI2020

CPU and memory efficient coherent mode decomposition for the partially coherent x-ray simulations

In this work, the method of the Coherent Mode Decomposition (CMD) is applied to numerical wave propagation calculations for partially-coherent X-rays, using the Fourier optics and compatible methods. Its CPU and memory efficiency is discussed in various cases of the wavefront at the source and the beam waist. With the absence of the quadratic phase terms, the required sampling density of the electric fields is effectively reduced. The problem size is thus moderate and the method is feasible to be implemented on a single-node CPU server. In other cases, the same argument holds with proper treatments of the quadratic phase terms. Tests on CMD and the modes propagation are done for the case of the Coherent Hard X-ray beamline of the National Synchrotron Light Source II, using the Synchrotron Radiation Workshop software. We observe a few hundred or less dominant decomposed modes that resemble the electric fields converge to the wavefront intensity at a high accuracy of over 99%.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Coherent Mode Decomposition for Kolmogorov Optical Turbulence in a Finite Aperture

An analysis of the coherent mode decomposition of an optical field after propagation through atmospheric turbulence is presented. The coherent modes represent an ideal basis by which to decompose the field for design of mode-limited optical systems. Using the rotational symmetry of the Fredholm integral operator for Kolmogorov optical turbulence, it is shown that the coherent modes exhibit separable solutions classified by radial and azimuthal quantum numbers. The study of the coherent modes is then reduced that of the radial functions determined by the aperture-coherence diameter ratio and obscuration ratio. Analysis of the spectrum of eigenvalues yields sharp bounds on the efficiency of receivers using incoherent or coherent combining with mode-limited photonic devices in the presence of Kolmogorov turbulence. The effective number of modes needed to represent Kolmogorov optical turbulence is studied via the von Neumann entropy, purity, and largest eigenvalue, and the differences in the different definitions is discussed. The similarity of the coherent modes to linearly polarized (LP) fiber modes is quantified yielding a precise characterization of the maximum gain that can be achieved in mode-limited systems via mode shaping techniques. As a final application, a mode sorting technique is presented for optimally splitting power from atmospherically degraded light into a finite number of modes simultaneously maximizing total coupling efficiency and minimizing the average and instantaneous power ratio between channels.

free-space optical communication↗

NASA NTRSDate not supplied

Coherent Mode Decomposition for Kolmogorov Optical Turbulence in a Finite Aperture

An analysis of the coherent mode decomposition of an optical field after propagation through atmospheric turbulence is presented. The coherent modes represent an ideal basis by which to decompose the field for design of mode-limited optical systems. Due to rotational symmetry of the Fredholm integral operator for Kolmogorov optical turbulence, the coherent modes exhibit separable solutions classified by a radial and azimuthal mode index. The study of the coherent modes is then reduced to that of the radial functions determined by the coherence ratio and obscuration ratio. Analysis of the spectrum of eigenvalues yields sharp bounds on the efficiency of receivers using incoherent or coherent combining with mode-limited photonic devices in the presence of Kolmogorov turbulence. The effective number of modes needed to represent Kolmogorov optical turbulence is studied via the von Neumann entropy, purity, and largest eigenvalue, and the differences in the different definitions is discussed. The similarity of the coherent modes to linearly polarized (LP) fiber modes is quantified yielding a precise characterization of the maximum gain that can be achieved in mode-limited systems via mode shaping techniques. As a final application, a mode sorting technique is presented for optimally splitting power from atmospherically degraded light into a finite number of modes simultaneously maximizing total coupling efficiency and minimizing the difference in average power between channels.

free-space optical communication↗

DOE OSTI2024

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Developments in SRW Code and Sirepo Framework Supporting Simulation of Time-Dependent Coherent X-ray Scattering Experiments

Physical optics simulations for beamlines and experiments are essential for the effective use of synchrotron light source facilities such as NSLS-II at BNL. The SRW software package supports such source-to-detector simulations for coherent X-ray scattering and imaging experiments through its Python interface and Sirepo browser-based graphical user interface. This allows one to define custom sample models, assess the feasibility of an experiment, and estimate most appropriate beamline settings before using valuable beamtime. We discuss the recent use of general-purpose GPU resources and coherent mode decomposition algorithms in SRW to accelerate physical optics simulations with partially coherent X-rays. To illustrate these new capabilities, we describe simulations of typical time series of partially coherent scattering images used in X-ray Photon Correlation Spectroscopy (XPCS) experiments; aiming to characterize the nanoscale dynamics of a disordered sample, representing a solution of nanoparticles undergoing Brownian diffusion.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Physical optics simulations for synchrotron radiation sources

Here we describe approaches to high-accuracy physical optics calculations used for the development of x-ray beamlines at synchrotron radiation sources, as well as simulation of experiments and processing of experimental data at some of these beamlines. We pay special attention to the treatment of the partial coherence of x rays, a topic of high practical importance for modern low-emittance high-brightness synchrotron radiation facilities. The approaches are based, to a large extent, on the works of Emil Wolf and co-authors, including the basic scalar diffraction theory and the coherent mode decomposition method. The presented simulation examples are related to the case of the novel Coherent Diffractive Imaging beamline that is currently under development at the National Synchrotron Light Source II at the Brookhaven National Laboratory.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Three-dimensional phase evolution and stress-induced non-uniform Li intercalation behavior in lithium iron phosphate

A lithium iron phosphate (LiFePO 4 ) olivine cathode displays rich phase transition behavior during charging/discharging. Despite extensive studies, the evolution of the LiFePO 4 /FePO 4 two-phase coexistence structure and the underlying mechanism are still a point of debate due to the intricate and often disparate phase morphologies observed. Here we apply phase-field simulations to provide detailed insights into the three-dimensional phase evolution process in LiFePO 4 upon delithiation under a large driving force. The simulation reveals that the coherency stress arising from the LiFePO 4 /FePO 4 lattice mismatch destabilizes the initially flat delithiation front and induces filamentary growth of the FePO 4 phase along [010], which is in excellent agreement with a recent operando X-ray imaging study. The combined simulation and experimental results provide the first illustration of the phenomenon of stress-induced instability of the lithium (de)intercalation front in battery electrode compounds. We show that this phenomenon is facilitated by the surface mode of coherent spinodal decomposition in LiFePO 4 and influenced by the anisotropies of misfit strain, elasticity and lithium diffusion. It is detrimental to battery performance and life by causing stress concentration and reducing the lithium (de)intercalation kinetics and thus should be mitigated.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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↗

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 OSTI2021

Transverse mode analysis for free-space laser beams using Bayesian analysis

Reliable and rapid assessment of the transverse mode quality of a free-space laser beam has a wide range of applications in laser development, research, and utilization. It has become even more important with recent advances in developing orbital angular momentum photon beams across a broad spectral region. In this work, a general modal analysis method for a free-space multimode laser beam has been developed based on Bayesian analysis. After transforming mode decomposition into a linear system problem, a Gaussian probabilistic model is used to find a closed-form solution. The method is found to be robust with the presence of Gaussian noise. Prior knowledge about the mode content can be incorporated into the method to improve the solution for situations when coherent disturbances or contamination are present in the laser beam. This method can be used to analyze the mode content for laser beams in different bases, such as Hermite–Gaussian (HG) modes and Laguerre–Gaussian (LG) modes. Three applications of this method are presented: a detailed modal analysis of the beam image from the incoherent intensity addition of HG modes and two examples of mode decomposition using the complex wavefront from the coherent superposition of HG and LG modes. The feasibility of this method is demonstrated using various simulation results. Based on digital images of a laser beam recorded without complex wavelength-limiting optics, in principle, this method can be used in a wide wavelength range from infrared to ultraviolet, and possibly x ray.

Liu, Peifan (ORCID:0000000273816153)↗

DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1994

Global-scale modes of surface temperature variability on interannual to century timescales

Using 100 years of global temperature anomaly data, we have performed a singluar value decomposition of temperature variations in narrow frequency bands to isolate coherent spatio-temporal modes of global climate variability. Statistical significance is determined from confidence limits obtained by Monte Carlo simulations. Secular variance is dominated by a globally coherent trend; with nearly all grid points warming in phase at varying amplitude. A smaller, but significant, share of the secular variance corresponds to a pattern dominated by warming and subsequent cooling in the high latitude North Atlantic with a roughly centennial timescale. Spatial patterns associated with significant peaks in variance within a broad period range from 2.8 to 5.7 years exhibit characteristic El Nino-Southern Oscillation (ENSO) patterns. A recent transition to a regime of higher ENSO frequency is suggested by our analysis. An interdecadal mode in the 15-to-18 years period and a mode centered at 7-to-8 years period both exhibit predominantly a North Atlantic Oscillation (NAO) temperature pattern. A potentially significant decadal mode centered on 11-to-12 years period also exhibits an NAO temperature pattern and may be modulated by the century-scale North Atlantic variability.

Mann, Michael E.↗

NASA NTRSDate not supplied

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗

DOE OSTI2020

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗