Compositional maps for registration in complex geometries

Abstract

We develop and analyze a parametric registration procedure for manifolds associated with the solutions to parametric partial differential equations in two-dimensional domains. Given the domain $\Omega \subset \mathbb{R}^2$ and the manifold $M=\{ u_{\mu} : \mu\in P\}$ associated with the parameter domain $P \subset \mathbb{R}^P$ and the parametric field $\mu\mapsto u_{\mu} \in L^2(\Omega)$, our approach takes as input a set of snapshots from $M$ and returns a parameter-dependent mapping $\Phi: \Omega \times P \to \Omega$, which tracks coherent features (e.g., shocks, shear layers) of the solution field and ultimately simplifies the task of model reduction. We consider mappings of the form $\Phi=\texttt{N}(\mathbf{a})$ where $\texttt{N}:\mathbb{R}^M \to {\rm Lip}(\Omega; \mathbb{R}^2)$ is a suitable linear or nonlinear operator; then, we state the registration problem as an unconstrained optimization statement for the coefficients $\mathbf{a}$. We identify minimal requirements for the operator $\texttt{N}$ to ensure the satisfaction of the bijectivity constraint; we propose a class of compositional maps that satisfy the desired requirements and enable non-trivial deformations over curved (non-straight) boundaries of $\Omega$; we develop a thorough analysis of the proposed ansatz for polytopal domains and we discuss the approximation properties for general curved domains. We perform numerical experiments for a parametric inviscid transonic compressible flow past a cascade of turbine blades to illustrate the many features of the method.