AI summary · arXiv: 2209.02244v1 · math.NA

The mpEDMD Algorithm for Data-Driven Computations of Measure-Preserving Dynamical Systems

Matthew J. Colbrook
Original paper Submitted to arXiv on 06 Sep. 2022 License: NONEXCLUSIVE-DISTRIB 1.0

How to read this summary? Illustrated guide

The essentials

  • Study focuses on analysis and decomposition of nonlinear dynamical systems using spectral information of Koopman operators
  • Koopman operators have ability to globally linearize systems, but computing their spectral information is challenging due to infinite-dimensional nature
  • Authors propose new algorithm called measure-preserving extended dynamic mode decomposition (mpEDMD)
  • mpEDMD is data-driven algorithm that enforces measure-preserving truncations of Koopman operators using orthogonal Procrustes problem
  • Operates with general dictionary of observables, making it flexible and compatible with various types of data and pre-existing DMD-type methods
  • Key contribution of mpEDMD is its ability to converge to spectral quantities of Koopman operators for general measure-preserving dynamical systems
  • Theoretical proofs provided for convergence of mpEDMD in terms of projection-valued and scalar-valued spectral measures, spectra, and Koopman mode decompositions
  • Investigation into convergence rates of approximating spectral measures when using delay embedding (Krylov subspaces) with increasing dictionary size
  • Demonstrated effectiveness through several challenging examples, showing increased robustness to noise compared to other DMD-type methods
  • Successfully captures energy conservation and cascade phenomena in experimental measurements of turbulent boundary layer flow with high Reynolds number (>6×10^4) and large state-space dimension (>10^5)
  • Overall, mpEDMD offers promising approach for analyzing nonlinear dynamical systems by providing accurate approximations of Koopman operator spectral information.

Key points written by AI: they may contain errors, check them against the paper.

AI-generated keywords: Koopman operators mpEDMD spectral information measure-preserving dynamical systems data-driven algorithm

Read further

  • Comprehensive Summary
  • Layman's Summary
  • Blog article
  • Original abstract
The study focuses on the analysis and decomposition of nonlinear dynamical systems using the spectral information of Koopman operators. These operators have the ability to globally linearize such systems, but their infinite-dimensional nature makes computing their spectral information challenging. To address this issue, the authors propose a new algorithm called measure-preserving extended dynamic mode decomposition (mpEDMD). mpEDMD is a data-driven algorithm that utilizes an orthogonal Procrustes problem to enforce measure-preserving truncations of Koopman operators. It operates with a general dictionary of observables, making it flexible and compatible with various types of data and pre-existing DMD-type methods. The key contribution of mpEDMD is its ability to converge to the spectral quantities of Koopman operators for general measure-preserving dynamical systems. The authors provide theoretical proofs for the convergence of mpEDMD in terms of projection-valued and scalar-valued spectral measures, spectra, and Koopman mode decompositions. They also investigate the convergence rates of approximating spectral measures when using delay embedding (Krylov subspaces) with increasing dictionary size. To demonstrate the effectiveness of mpEDMD, several challenging examples are presented. The algorithm shows increased robustness to noise compared to other DMD-type methods. Additionally, it successfully captures energy conservation and cascade phenomena in experimental measurements of a turbulent boundary layer flow with high Reynolds number (>6×10^4) and large state-space dimension (>10^5). Overall, mpEDMD offers a promising approach for analyzing nonlinear dynamical systems by providing accurate approximations of Koopman operator spectral information. Its flexibility, ease-of-use, and demonstrated performance make it a valuable tool for researchers in various fields dealing with complex dynamical systems.

Was this synthesis useful?

Score: 0

Why do we need votes?

Votes are used to determine whether we need to re-run our summarizing tools. If the count reaches -10, our tools can be restarted.

Figures and plots

Ask about a figure by number or use its question button to request visual analysis. Up to two figures and four image panels per question. Values read from plots are approximate; check the original data.

AI assistant instructions?

Tools for this paper

Explain it differently 3 free rewrites left today
Save to a collection Group papers in collections to compare them, chat with them and export them. Create a free account Log in

Advanced tools Premium

For researchers who follow a field closely: citation alerts, journal-club slides, version comparison and Zotero export. Discover Premium

Zotero
Citation alerts Citation alerts guide

Papers around this one

See them on a map

Ranked by how close their title and abstract are to this paper’s. Select some to read them together.

Closest papers

  1. Related
    DelayDiffEq: Generating Delay Differential Equation Solvers via Recursive Embedding of Ordinary Differential Equation Solvers
    David Widmann, Chris Rackauckas · 2022 · Numerical Analysis · AI summary available
    Abstract

    Traditional solvers for delay differential equations (DDEs) are designed around only a single method and do not effectively use the infrastructure of their more-developed ordinary differential equation (ODE) counterparts. In this work we present DelayDiffEq, a Julia package for numerically solving delay differential equations (DDEs) which leverages the multitude of numerical algorithms in OrdinaryDiffEq for solving both stiff and non-stiff ODEs, and manages to solve challenging stiff DDEs. We describe how compiling the ODE integrator within itself, and accounting for discontinuity propagation, leads to a design that is effective for DDEs while using all of the ODE internals. We highlight som…

  2. Related
    Low-complexity computation of plate eigenmodes with Vekua approximations and the Method of Particular Solutions
    Gilles Chardon, Laurent Daudet · 2013 · Numerical Analysis · AI summary available
    Abstract

    This paper extends the Method of Particular Solutions (MPS) to the computation of eigenfrequencies and eigenmodes of plates. Specific approximation schemes are developed, with plane waves (MPS-PW) or Fourier-Bessel functions (MPS-FB). This framework also requires a suitable formulation of the boundary conditions. Numerical tests, on two plates with various boundary conditions, demonstrate that the proposed approach provides competitive results with standard numerical schemes such as the Finite Element Method, at reduced complexity, and with large flexibility in the implementation choices.

  3. Related
    A robust solution strategy for the Cahn-Larché equations
    Erlend Storvik, Jakub Wiktor Both, Jan Martin Nordbotten et al. · 2022 · Numerical Analysis · AI summary available
    Abstract

    In this paper we propose a solution strategy for the Cahn-Larch\'e equations, which is a model for linearized elasticity in a medium with two elastic phases that evolve subject to a Ginzburg-Landau type energy functional. The system can be seen as a combination of the Cahn-Hilliard regularized interface equation and linearized elasticity, and is non-linearly coupled, has a fourth order term that comes from the Cahn-Hilliard subsystem, and is non-convex and nonlinear in both the phase-field and displacement variables. We propose a novel semi-implicit discretization in time that uses a standard convex-concave splitting method of the nonlinear double-well potential, as well as special treatment…

  4. Related
    A Computationally Efficient Vectorized Implementation of Localizing Gradient Damage Method in MATLAB
    Subrato Sarkar · 2023 · Numerical Analysis · AI summary available
    Abstract

    In this work, a recently developed fracture modeling method called localizing gradient damage method (LGDM) is implemented in MATLAB. MATLAB is well-known in the computational research community for its simple and easy-to-learn coding interface. As a result, MATLAB is generally preferred for the initial development (prototyping) of computational models by researchers. However, MATLAB-developed codes are seldom used for large-scale simulations (after initial development is complete) due to their computational inefficiency. Hence, a computationally efficient implementation of LGDM using MATLAB vectorization is presented in this work. The choice of LGDM (as the fracture modeling method) is base…

  5. Related
    An overview on deep learning-based approximation methods for partial differential equations
    Christian Beck, Martin Hutzenthaler, Arnulf Jentzen et al. · 2020 · Discrete Contin. Dyn. Syst. Ser. B 28 (2023), no. 6, 3697-3746 · AI summary available
    Abstract

    It is one of the most challenging problems in applied mathematics to approximatively solve high-dimensional partial differential equations (PDEs). Recently, several deep learning-based approximation algorithms for attacking this problem have been proposed and tested numerically on a number of examples of high-dimensional PDEs. This has given rise to a lively field of research in which deep learning-based methods and related Monte Carlo methods are applied to the approximation of high-dimensional PDEs. In this article we offer an introduction to this field of research by revisiting selected mathematical results related to deep learning approximation methods for PDEs and reviewing the main ide…

  6. Related
    Compatible finite element methods for geophysical fluid dynamics
    Colin J. Cotter · 2023 · Numerical Analysis · AI summary available
    Abstract

    This article surveys research on the application of compatible finite element methods to large scale atmosphere and ocean simulation. Compatible finite element methods extend Arakawa's C-grid finite difference scheme to the finite element world. They are constructed from a discrete de Rham complex, which is a sequence of finite element spaces which are linked by the operators of differential calculus. The use of discrete de Rham complexes to solve partial differential equations is well established, but in this article we focus on the specifics of dynamical cores for simulating weather, oceans and climate. The most important consequence of the discrete de Rham complex is the Hodge-Helmholtz d…

  7. Related
    Provably stable numerical method for the anisotropic diffusion equation in toroidally confined magnetic fields
    Dean Muir, Kenneth Duru, Matthew Hole et al. · 2023 · Numerical Analysis · AI summary available
    Abstract

    We present a novel numerical method for solving the anisotropic diffusion equation in toroidally confined magnetic fields which is efficient, accurate and provably stable. The continuous problem is written in terms of a derivative operator for the perpendicular transport and a linear operator, obtained through field line tracing, for the parallel transport. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing …

  8. Related
    Convergence Analysis for Nonlinear GMRES
    Yunhui He · 2025 · Numerical Analysis · AI summary available
    Abstract

    In this work, we revisit nonlinear generalized minimal residual method (NGMRES) applied to nonlinear problems. NGMRES is used to accelerate the convergence of fixed-point iterations, which can substantially improve the performance of the underlying fixed-point iterations. We consider NGMRES with a finite window size $m$, denoted as NGMRES($m$). However, there is no convergence analysis for NGMRES($m$) applied to nonlinear systems. We prove that for general $m>0$, the residuals of NGMRES($m$) converge r-linearly under some conditions. For $m=0$, we prove that the residuals of NGMRES(0) converge q-linearly.

How are these papers found?

We compare the title and abstract of this paper with those of every paper summarized on this site, and with the papers Semantic Scholar recommends, those it cites and those citing it (for a much-cited paper, its most recent citations and, from OpenAlex, the most cited ones). Among the prior work and the follow-ups, the papers that use each other in depth, a method it uses and the much-cited follow-ups come a little higher. The levels (very close, close, related) come from the similarity of the texts, not from a judgment of quality; * marks a paper without an abstract, compared by its title. The sentence under a paper is written by an AI from the two abstracts: check it in the papers. A paper marked not relevant by three readers leaves the list.

For now only among the papers summarized on this site: search again to add the prior work, the follow-ups and Semantic Scholar’s recommendations.

The previous summary was created more than a year ago and can be re-run (if necessary) by clicking on the Run button below.

Disclaimer: The AI-based summarization tool and virtual assistant provided on this website may not always provide accurate and complete summaries or responses. We encourage you to carefully review and evaluate the generated content to ensure its quality and relevance to your needs.