Dynamic mode decomposition for data-driven analysis and reduced-order modeling of E × B plasmas: I. Extraction of spatiotemporally coherent patterns

The advent of data-driven/machine-learning based methods and the increase in data available from high-fidelity simulations and experiments has opened new pathways toward realizing reduced-order models for plasma systems that can aid in explaining the complex, multi-dimensional phenomena and enable forecasting and prediction of the systems’ behavior. In this two-part article, we evaluate the utility and the generalizability of the dynamic mode decomposition (DMD) algorithm for data-driven analysis and reduced-order modeling of plasma dynamics in cross-field E × B configurations. The DMD algorithm is an interpretable data-driven method that finds a best-fit linear model describing the time evolution of spatiotemporally coherent structures (patterns) in data. We have applied the DMD to extensive high-fidelity datasets generated using a particle-in-cell (PIC) code based on the cost-efficient reduced-order PIC scheme. In this part, we first provide an overview of the concept of DMD and its underpinning proper orthogonal and singular value decomposition methods. Two of the main DMD variants are next introduced. We then present and discuss the results of the DMD application in terms of the identification and extraction of the dominant spatiotemporal modes from high-fidelity data over a range of simulation conditions. We demonstrate that the DMD variant based on variable projection optimization (OPT-DMD) outperforms the basic DMD method in identification of the modes underlying the data, leading to notably more reliable reconstruction of the ground-truth. Furthermore, we show in multiple test cases that the discrete frequency spectrum of OPT-DMD-extracted modes is consistent with the temporal spectrum from the fast Fourier transform of the data. This observation implies that the OPT-DMD augments the conventional spectral analyses by being able to uniquely reveal the spatial structure of the dominant modes in the frequency spectra, thus, yielding more accessible, comprehensive information on the spatiotemporal characteristics of the plasma phenomena.