skip to main content


Title: On Occupation Kernels, Liouville Operators, and Dynamic Mode Decomposition
Using the newly introduced ``occupation kernels,'' the present manuscript develops an approach to dynamic mode decomposition (DMD) that treats continuous time dynamics, without discretization, through the Liouville operator. The technical and theoretical differences between Koopman based DMD for discrete time systems and Liouville based DMD for continuous time systems are highlighted, which includes an examination of these operators over several reproducing kernel Hilbert spaces.  more » « less
Award ID(s):
2028001
NSF-PAR ID:
10281591
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
2021 American Control Conference (ACC)
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Koopman decomposition is a nonlinear generalization of eigen-decomposition, and is being increasingly utilized in the analysis of spatio-temporal dynamics. Well-known techniques such as the dynamic mode decomposition (DMD) and its linear variants provide approximations to the Koopman operator, and have been applied extensively in many fluid dynamic problems. Despite being endowed with a richer dictionary of nonlinear observables, nonlinear variants of the DMD, such as extended/kernel dynamic mode decomposition (EDMD/KDMD) are seldom applied to large-scale problems primarily due to the difficulty of discerning the Koopman-invariant subspace from thousands of resulting Koopman eigenmodes. To address this issue, we propose a framework based on a multi-task feature learning to extract the most informative Koopman-invariant subspace by removing redundant and spurious Koopman triplets. In particular, we develop a pruning procedure that penalizes departure from linear evolution. These algorithms can be viewed as sparsity-promoting extensions of EDMD/KDMD. Furthermore, we extend KDMD to a continuous-time setting and show a relationship between the present algorithm, sparsity-promoting DMD and an empirical criterion from the viewpoint of non-convex optimization. The effectiveness of our algorithm is demonstrated on examples ranging from simple dynamical systems to two-dimensional cylinder wake flows at different Reynolds numbers and a three-dimensional turbulent ship-airwake flow. The latter two problems are designed such that very strong nonlinear transients are present, thus requiring an accurate approximation of the Koopman operator. Underlying physical mechanisms are analysed, with an emphasis on characterizing transient dynamics. The results are compared with existing theoretical expositions and numerical approximations. 
    more » « less
  2. Abstract

    Uvula‐induced snoring and associated obstructive sleep apnea is a complex phenomenon characterized by vibrating structures and highly transient vortex dynamics. This study aimed to extract signature features of uvula wake flows of different pathological origins and develop a linear reduced‐order surrogate model for flow control. Six airway models were developed with two uvula kinematics and three pharynx constriction levels. A direct numerical simulation (DNS) flow solver based on the immersed boundary method was utilized to resolve the wake flows induced by the flapping uvula. Key spatial and temporal responses of the flow to uvula kinematics and pharynx constriction were investigated using continuous wavelet transform (CWT), proper orthogonal decomposition (POD), and dynamic mode decomposition (DMD). Results showed highly complex patterns in flow topologies. CWT analysis revealed multiscale correlations in both time and space between the flapping uvular and wake flows. POD analysis successfully separated the flows among the six models by projecting the datasets in the vector space spanned by the first three eigenmodes. Perceivable differences were also captured in the time evolution of the DMD modes among the six models. A linear reduced‐order surrogate model was constructed from the predominant eigenmodes obtained from the DMD analysis and predicted vortex patterns from this surrogate model agreed well with the corresponding DNS simulations. The computational and analytical platform presented in this study could bring a variety of applications in breathing‐related disorders and beyond. The computational efficiency of surrogate modeling makes it well suited for flow control, forecasting, and uncertainty analyses.

     
    more » « less
  3. null (Ed.)
    Abstract For $$\gamma \in (0,2)$$ γ ∈ ( 0 , 2 ) , we define a weak $$\gamma $$ γ - Liouville quantum gravity ( LQG ) metric to be a function $$h\mapsto D_h$$ h ↦ D h which takes in an instance of the planar Gaussian free field and outputs a metric on the plane satisfying a certain list of natural axioms. We show that these axioms are satisfied for any subsequential limits of Liouville first passage percolation. Such subsequential limits were proven to exist by Ding et al. (Tightness of Liouville first passage percolation for $$\gamma \in (0,2)$$ γ ∈ ( 0 , 2 ) , 2019. ArXiv e-prints, arXiv:1904.08021 ). It is also known that these axioms are satisfied for the $$\sqrt{8/3}$$ 8 / 3 -LQG metric constructed by Miller and Sheffield (2013–2016). For any weak $$\gamma $$ γ -LQG metric, we obtain moment bounds for diameters of sets as well as point-to-point, set-to-set, and point-to-set distances. We also show that any such metric is locally bi-Hölder continuous with respect to the Euclidean metric and compute the optimal Hölder exponents in both directions. Finally, we show that LQG geodesics cannot spend a long time near a straight line or the boundary of a metric ball. These results are used in subsequent work by Gwynne and Miller which proves that the weak $$\gamma $$ γ -LQG metric is unique for each $$\gamma \in (0,2)$$ γ ∈ ( 0 , 2 ) , which in turn gives the uniqueness of the subsequential limit of Liouville first passage percolation. However, most of our results are new even in the special case when $$\gamma =\sqrt{8/3}$$ γ = 8 / 3 . 
    more » « less
  4. Abstract This work presents a rigorous theory for topological photonic materials in one dimension. The main focus is on the existence of interface modes that are induced by topological properties of the bulk structure. For a general 1D photonic structure with time-reversal symmetry, we investigate the existence of an interface mode that is induced by a Dirac point upon perturbation. Specifically, we establish conditions on the perturbation which guarantee the opening of a band gap around the Dirac point and the existence of an interface mode. For a periodic photonic structure with both time-reversal and inversion symmetry, the Zak phase is quantized, taking only two values 0 , π . We show that the Zak phase is determined by the parity (even or odd) of the Bloch modes at the band edges. For a photonic structure consisting of two semi-infinite systems on the two sides of an interface with distinct topological indices, we show the existence of an interface mode inside the common gap. The stability of the mode under perturbations is also investigated. Finally, we study resonances for finite topological structures. Our results are based on the transfer matrix method and the oscillation theory for Sturm–Liouville operators. The methods and results can be extended to general topological Sturm–Liouville systems in one dimension. 
    more » « less
  5. Microsphere photolithography (MPL) is a promising technique for cost-effective fabrication of large-scale metasurfaces. This approach generates an array of photonic jets by the collimated illumination of self-assembled microspheres. The photonic jets can be precisely steered within the unit cell defined by each microsphere by changing the angle of incidence. This allows for the creation of complex metasurface element geometries. Computer controlled articulation of the substrate relative to a static UV source allows the direct-write of different metasurface elements. However, this is time-consuming and requires registration between each exposure for complex features. This paper investigates a single exposure method with the dynamic continuous angle of incidence control provided by a Digital Micromirror Device (DMD) in the front Fourier plane of the projection system. The grayscale values of the DMD pixels can be adjusted to provide optical proximity correction. Larger patterns can be achieved by scanning the substrate relative to the exposure beam. This approach is demonstrated with the creation of hierarchical patterns. This work greatly simplifies the MPL exposure process for complex resonators and provides potential for full light field control.

     
    more » « less