Reduction of nonlinear embedded boundary models for problems with evolving interfaces

Embedded boundary methods alleviate many computational challenges, including those associated with meshing complex geometries and solving problems with evolving domains and interfaces. Developing model reduction methods for computational frameworks based on such methods seems however to be challenging. Indeed, most popular model reduction techniques are projection-based, and rely on basis functions obtained from the compression of simulation snapshots. In a traditional interface-fitted computational framework, the computation of such basis functions is straightforward, primarily because the computational domain does not contain in this case a fictitious region. This is not the case however for an embedded computational framework because the computational domain typically contains in this case both real and ghost regions whose definitions complicate the collection and compression of simulation snapshots. The problem is exacerbated when the interface separating both regions evolves in time. This paper addresses this issue by formulating the snapshot compression problem as a weighted low-rank approximation problem where the binary weighting identifies the evolving component of the individual simulation snapshots. The proposed approach is application independent and therefore comprehensive. It is successfully demonstrated for the model reduction of several two-dimensional, vortex-dominated, fluid-structure interaction problems.

[1]  Pei Chen,et al.  Optimization Algorithms on Subspaces: Revisiting Missing Data Problem in Low-Rank Matrix , 2008, International Journal of Computer Vision.

[2]  C. Farhat,et al.  Interpolation Method for Adapting Reduced-Order Models and Application to Aeroelasticity , 2008 .

[3]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[4]  David A. Ham,et al.  POD reduced-order unstructured mesh modeling applied to 2D and 3D fluid flow , 2013, Comput. Math. Appl..

[5]  E. Candès The restricted isometry property and its implications for compressed sensing , 2008 .

[6]  Earl H. Dowell,et al.  Stabilization of projection-based reduced order models of the Navier–Stokes , 2012 .

[7]  T. Belytschko,et al.  Robust and provably second‐order explicit–explicit and implicit–explicit staggered time‐integrators for highly non‐linear compressible fluid–structure interaction problems , 2010 .

[8]  Siep Weiland,et al.  Missing Point Estimation in Models Described by Proper Orthogonal Decomposition , 2004, IEEE Transactions on Automatic Control.

[9]  Raphaël Loubère,et al.  ReALE: A reconnection-based arbitrary-Lagrangian-Eulerian method , 2010, J. Comput. Phys..

[10]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[11]  C. Farhat,et al.  Dimensional reduction of nonlinear finite element dynamic models with finite rotations and energy‐based mesh sampling and weighting for computational efficiency , 2014 .

[12]  Xianyi Zeng,et al.  A systematic approach for constructing higher-order immersed boundary and ghost fluid methods for fluid-structure interaction problems , 2012, J. Comput. Phys..

[13]  Juan Du,et al.  Non-linear model reduction for the Navier-Stokes equations using residual DEIM method , 2014, J. Comput. Phys..

[14]  Takayuki Okatani,et al.  On the Wiberg Algorithm for Matrix Factorization in the Presence of Missing Components , 2007, International Journal of Computer Vision.

[15]  Charbel Farhat,et al.  Reduced-order fluid/structure modeling of a complete aircraft configuration , 2006 .

[16]  Karen Willcox,et al.  Parametric reduced-order models for probabilistic analysis of unsteady aerodynamic applications , 2007 .

[17]  Charbel Farhat,et al.  The GNAT method for nonlinear model reduction: Effective implementation and application to computational fluid dynamics and turbulent flows , 2012, J. Comput. Phys..

[18]  Danny C. Sorensen,et al.  Nonlinear Model Reduction via Discrete Empirical Interpolation , 2010, SIAM J. Sci. Comput..

[19]  Ivan Markovsky,et al.  Software for weighted structured low-rank approximation , 2014, J. Comput. Appl. Math..

[20]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[21]  Charbel Farhat,et al.  FIVER: A finite volume method based on exact two-phase Riemann problems and sparse grids for multi-material flows with large density jumps , 2012, J. Comput. Phys..

[22]  Jacob K. White,et al.  Model order reduction for nonlinear dynamical systems based on trajectory piecewise-linear approximations , 2006 .

[23]  Earl H. Dowell,et al.  Low-dimensional modelling of high-Reynolds-number shear flows incorporating constraints from the Navier–Stokes equation , 2013, Journal of Fluid Mechanics.

[24]  Charbel Farhat,et al.  A three-dimensional torsional spring analogy method for unstructured dynamic meshes , 2002 .

[25]  Kevin G. Wang,et al.  Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems , 2011 .

[26]  Andrew W. Fitzgibbon,et al.  Damped Newton algorithms for matrix factorization with missing data , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[27]  Ionel M. Navon,et al.  A Dual-Weighted Approach to Order Reduction in 4DVAR Data Assimilation , 2008, Monthly Weather Review.

[28]  A. Patera,et al.  Certified real‐time solution of the parametrized steady incompressible Navier–Stokes equations: rigorous reduced‐basis a posteriori error bounds , 2005 .

[29]  C. Farhat,et al.  Efficient non‐linear model reduction via a least‐squares Petrov–Galerkin projection and compressive tensor approximations , 2011 .

[30]  Prateek Jain,et al.  Low-rank matrix completion using alternating minimization , 2012, STOC '13.

[31]  Ivan Markovsky,et al.  Low Rank Approximation - Algorithms, Implementation, Applications , 2018, Communications and Control Engineering.

[32]  Marco Fossati,et al.  An ALE mesh movement scheme for long‐term in‐flight ice accretion , 2012 .

[33]  C. Peskin The immersed boundary method , 2002, Acta Numerica.

[34]  C. Farhat,et al.  Mixed explicit/implicit time integration of coupled aeroelastic problems: Three‐field formulation, geometric conservation and distributed solution , 1995 .

[35]  Emmanuel J. Candès,et al.  Exact Matrix Completion via Convex Optimization , 2008, Found. Comput. Math..

[36]  Bart Vandereycken,et al.  Low-Rank Matrix Completion by Riemannian Optimization , 2013, SIAM J. Optim..

[37]  Earl H. Dowell,et al.  Reduced-Order Modeling of Flutter and Limit-Cycle Oscillations Using the Sparse Volterra Series , 2012 .

[38]  C. Farhat,et al.  Torsional springs for two-dimensional dynamic unstructured fluid meshes , 1998 .

[39]  J. Halleux,et al.  An arbitrary lagrangian-eulerian finite element method for transient dynamic fluid-structure interactions , 1982 .

[40]  Tommi S. Jaakkola,et al.  Weighted Low-Rank Approximations , 2003, ICML.

[41]  Rama Chellappa,et al.  Large-Scale Matrix Factorization with Missing Data under Additional Constraints , 2010, NIPS.