An optimization-based registration approach to geometry reduction

Abstract: We develop and assess an optimization-based approach to parametric geometry reduction. Given a family of parametric domains, we aim to determine a parametric diffeomorphism $\Phi$ that maps a fixed reference domain $\Omega$ into each element of the family, for different values of the parameter; the ultimate goal of our study is to determine an effective tool for parametric projection-based model order reduction of partial differential equations in parametric geometries. For practical problems in engineering, explicit parameterizations of the geometry are likely unavailable: for this reason, our approach takes as inputs a reference mesh of $\Omega$ and a point cloud ${y_i^{\rm raw}}{i=1}^Q$ that belongs to the boundary of the target domain $V$ and returns a bijection $\Phi$ that approximately maps $\Omega$ in $V$. We propose a two-step procedure: given the point clouds ${x_j}{j=1}^N\subset \partial \Omega$ and ${y_i^{\rm raw}}{i=1}^Q \subset \partial V$, we first resort to a point-set registration algorithm to determine the displacements ${ v_j }{j=1}^N$ such that the deformed point cloud ${y_j:= x_j+v_j }_{j=1}^N$ approximates $\partial V$; then, we solve a nonlinear non-convex optimization problem to build a mapping $\Phi$ that is bijective from $\Omega$ in $\mathbb{R}^d$ and (approximately) satisfies $\Phi(x_j) = y_j$ for $j=1,\ldots,N$.We present a rigorous mathematical analysis to justify our approach; we further present thorough numerical experiments to show the effectiveness of the proposed method.