J. Nathan Kutz Department of Applied Mathematics University of - - PowerPoint PPT Presentation

j nathan kutz
SMART_READER_LITE
LIVE PREVIEW

J. Nathan Kutz Department of Applied Mathematics University of - - PowerPoint PPT Presentation

DATA-DRIVEN MODEL DISCOVERY AND COORDINATE EMBEDDINGS FOR PHYSICAL SYSTEMS ICERM 2019 J. Nathan Kutz Department of Applied Mathematics University of Washington Email: kutz@uw.edu Web: faculty.washington.edu/kutz Model Discovery Steven


slide-1
SLIDE 1

DATA-DRIVEN MODEL DISCOVERY AND COORDINATE EMBEDDINGS FOR PHYSICAL SYSTEMS ICERM 2019

  • J. Nathan Kutz

Department of Applied Mathematics University of Washington Email: kutz@uw.edu Web: faculty.washington.edu/kutz

slide-2
SLIDE 2

Model Discovery

Steven Brunton

Mechanical Engineering University of Washington

Joshua Proctor

Institute for Disease Modeling

slide-3
SLIDE 3

Mathematical Framework

x

Dynamics Measurement

State-space Parameters Stochastic effects Measurement model Measurement noise Dynamics

slide-4
SLIDE 4

What Could the Right Side Be?

Limited by your imagination 2nd degree polynomials

slide-5
SLIDE 5

Sparse Identification of Nonlinear Dynamics (SINDy)

slide-6
SLIDE 6
slide-7
SLIDE 7

Identifying ROMs

slide-8
SLIDE 8

Discovery Innovations

Schaeffer et al -- corrupt data, PDEs, integral formulation, convergence Dongbin Xiu & co-workers (2018) – Sampling strategies Guang Lin & co-workers (2018) -- Uncertainty Metrics Hod Lipson and co-workers (2006) — Symbolic/genetic regression Karniadakis, Raissai, Perdikaris …. — Neural Nets

Zheng, Askham, Brunton, Kutz & Aravkin (2018) – SR3 sparse relaxed regularized regression (for SINDy, LASSO, CS, TV, Matrix Completion …)

slide-9
SLIDE 9

Manifolds and Embeddings Observables & Coordinates

slide-10
SLIDE 10

Bernard Koopman 1931

Mezic (2004) Coifman, Kevrekidis, co-workers

slide-11
SLIDE 11

Koopman Invariant Subspaces

Brunton, Proctor & Kutz, PLOS ONE (2018)

slide-12
SLIDE 12

Burgers’ Equation

Cole-Hopf

Kutz, Proctor & Brunton, Complexity (2018)

slide-13
SLIDE 13

Dynamic Mode Decomposition

Travis Askham

Askham & Kutz, SIADS (2018)

slide-14
SLIDE 14

Approximate Dynamical Systems

Linear dynamics (equation-free) Eigenfunction expansion Least-square fit

slide-15
SLIDE 15

DMD with Control

Input Input Snapshots DMD generalization

Proctor, Brunton & Kutz, SIADS (2016)

slide-16
SLIDE 16

Koopman vs DMD: All about Observables!

slide-17
SLIDE 17

Neural Nets

slide-18
SLIDE 18

The Zoo NN Zoo

slide-19
SLIDE 19

NNs for Koopman Embedding

Bethany Lusch

Lusch et al. Nat. Comm (2018)

slide-20
SLIDE 20

NNs for Koopman Embedding

Bethany Lusch

Lusch et al. Nat. Comm (2018)

Handling the Continuous Spectra

Lusch, Kutz & Brunton, arxiv (2017)

slide-21
SLIDE 21

The Pendulum

slide-22
SLIDE 22

Flow Around a Cylinder

slide-23
SLIDE 23

Autoencoder + SINDy

Kathleen Champion

slide-24
SLIDE 24

Fast Learning

Charles Delahunt

slide-25
SLIDE 25

Moth Olfactory System

slide-26
SLIDE 26

Learning New Odors

slide-27
SLIDE 27

Sparsity for Learning

slide-28
SLIDE 28

Rapid Learning in NNs

slide-29
SLIDE 29

Comparisons

slide-30
SLIDE 30

Decoder Networks

slide-31
SLIDE 31

Structure of Mapping

Linear Maps: SVD (left singular vector) defines layer Erichson, Mathelin, Brunton, Mahoney & Kutz (2019)

slide-32
SLIDE 32

Measurements State space Mapping Approximate the full state space from limited measurements Optimization

Mathematical Framework

slide-33
SLIDE 33

Linear Maps

Singular value decomposition Data Linear measurements H Optimize (least-squares)

slide-34
SLIDE 34

Optimal Placement

Point measurements Optimal Sensors via QR pivots Manohar, Kutz & Brunton (2018) IEEE Control Systems Magazine Clark, Askham, Brunton & Kutz (2019) IEEE Sensors

slide-35
SLIDE 35

Nonlinear Mapping

General Form: Compositional Layers Universal Approximators: Hornik 1990

slide-36
SLIDE 36

Shallow Layer Mapping

Two Layers Composition Erichson, Mathelin, Brunton, Mahoney & Kutz (2019) SIAM

slide-37
SLIDE 37

Activation Functions

slide-38
SLIDE 38

Linear vs Nonlinear Maps

Improved Interpretability of Modes

slide-39
SLIDE 39

Improved Performance

slide-40
SLIDE 40

Robustness to Noise

slide-41
SLIDE 41

Improved Performance

slide-42
SLIDE 42

Super Resolution Analysis

slide-43
SLIDE 43

Conclusion

Model Discovery: Sparse regression provides parsimonious dynamical models Coordinates: Learning Koopman embeddings can provide optimal basis for dynamics Neural Networks: Structure and function matter

  • Connect discovery and coordinates
  • Structure can lead to fast (one-shot) learning with limited data
  • Discovery of improved coordinate embeddings through decoding

COMING SOON: A multi scale physics challenge set