Projection-Based Dimensional Reduction of Adaptively Refined Nonlinear Models

Abstract

Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computational models such as projection-based reduced-order models. This paper presents a complete framework for projection-based model order reduction (PMOR) of nonlinear problems in the presence of AMR that builds on elements from existing methods and augments them with critical new contributions. In particular, it proposes an analytical algorithm for computing a pseudo-meshless inner product between adapted solution snapshots for the purpose of clustering and PMOR. It exploits hyperreduction—specifically, the energy-conserving sampling and weighting hyperreduction method—to deliver for nonlinear and/or parametric problems the desired computational gains. Most importantly, the proposed framework for PMOR in the presence of AMR capitalizes on the concept of state-local reduced-order bases to make the most of the notion of a supermesh, while achieving computational tractability. Its features are illustrated with CFD applications grounded in AMR and its significance is demonstrated by the reported wall-clock speedup factors.

Keywords

Adaptive mesh refinement (AMR) / Computational fluid dynamics / Energy-conserving sampling and weighting (ECSW) / Model order reduction / Reduced-order model / Supermesh

References

[1.] Amsallem D, Zahr MJ, Farhat C. Nonlinear model order reduction based on local reduced-order bases. Int. J. Numer. Meth. Eng., 2012, 92(10): 891-916,

[2.] Bank RE, Sherman AH. An adaptive, multi-level method for elliptic boundary value problems. Computing, 1981, 26(2): 91-105,

[3.] Barnett J, Farhat C. Quadratic approximation manifold for mitigating the Kolmogorov barrier in nonlinear projection-based model order reduction. J. Comput. Phys., 2022, 464(111): 348,

[4.] Borker R, Huang D, Grimberg S, et al.. Mesh adaptation framework for embedded boundary methods for computational fluid dynamics and fluid-structure interaction. Int. J. Numer. Meth. Fluids, 2019, 90(8): 389-424,

[5.] Carlberg K, Bou-Mosleh C, Farhat C. Efficient non-linear model reduction via a least-squares Petrov-Galerkin projection and compressive tensor approximations. Int. J. Numer. Meth. Eng., 2011, 86(2): 155-181,

[6.] Carlberg K, Farhat C, Cortial J, et al.. The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows. J. Comput. Phys., 2013, 242: 623-647,

[7.] Chaturantabut S, Sorensen DC. Nonlinear model reduction via discrete empirical interpolation. SIAM J. Sci. Comput., 2010, 32(5): 2737-2764,

[8.] Eftang JL, Patera AT, Rønquist EM. An hp certified reduced basis method for parametrized elliptic partial differential equations. SIAM J. Sci. Comput., 2010, 32(6): 3170-3200,

[9.] Elkan, C.: Using the triangle inequality to accelerate K-means. In: Proceedings of the Twentieth International Conference on Machine Learning. AAAI Press, ICML 2003, pp. 147–153 (2003)

[10.] Farhat C, Avery P, Chapman T, et al.. Dimensional reduction of nonlinear finite element dynamic models with finite rotations and energy-based mesh sampling and weighting for computational efficiency. Int. J. Numer. Meth. Eng., 2014, 98(9): 625-662,

[11.] Farhat C, Chapman T, Avery P. Structure-preserving, stability, and accuracy properties of the energy-conserving sampling and weighting method for the hyper reduction of nonlinear finite element dynamic models. Int. J. Numer. Meth. Eng., 2015, 102(5): 1077-1110,

[12.] Farhat C, Geuzaine P, Brown G. Application of a three-field nonlinear fluid-structure formulation to the prediction of the aeroelastic parameters of an F-16 fighter. Comput. Fluids, 2003, 32(1): 3-29,

[13.] Farhat, C., Grimberg, S., Manzoni, A., et al.: Computational bottlenecks for PROMs: pre-computation and hyperreduction. In: Benner, P., Grivet-Talocia, S., Quarteroni, A., et al. (eds.) Model Order Reduction - Volume 2: Snapshot-Based Methods and Algorithms. De Gruyter, Berlin, chap. 5, pp. 181–244 (2020). https://doi.org/10.1515/9783110671490-005

[14.] Farhat Research Group: Aero-F (2022). frg.bitbucket.io

[15.] Farrell P, Maddison J. Conservative interpolation between volume meshes by local Galerkin projection. Comput. Methods Appl. Mech. Eng., 2011, 200(1–4): 89-100,

[16.] Geuzaine P, Brown G, Harris C, et al.. Aeroelastic dynamic analysis of a full F-16 configuration for various flight conditions. AIAA J., 2003, 41(3): 363-371,

[17.] Gräßle C, Hinze M. POD reduced order modeling for evolution equations utilizing arbitrary finite element discretizations. Adv. Comput. Math., 2018, 44(6): 1941-1978,

[18.] Grimberg S, Farhat C, Tezaur R, et al.. Mesh sampling and weighting for the hyperreduction of nonlinear Petrov-Galerkin reduced-order models with local reduced-order bases. Int. J. Numer. Meth. Eng., 2021, 122(7): 1846-1874,

[19.] Grimberg S, Farhat C, Youkilis N. On the stability of projection-based model order reduction for convection-dominated laminar and turbulent flows. J. Comput. Phys., 2020, 419(109): 681,

[20.] Haasdonk B, Ohlberger M. Reduced basis method for finite volume approximations of parametrized linear evolution equations. ESAIM: Math. Model. Num. Anal., 2008, 42(2): 277-302,

[21.] Huang W. Mathematical principles of anisotropic mesh adaptation. Commun. Comput. Phys., 2006, 1(2): 276-310

[22.] Jiao X, Heath MT. Common-refinement-based data transfer between non-matching meshes in multiphysics simulations. Int. J. Numer. Meth. Eng., 2004, 61(14): 2402-2427,

[23.] Leutenegger, S.T., Lopez, M.A., Edgington, J.: STR: a simple and efficient algorithm for R-tree packing. In: Proceedings of the 13th International Conference on Data Engineering, IEEE, pp. 497–506 (1997)

[24.] Mitchell, W.F.: Unified multilevel adaptive finite element methods for elliptic problems. University of Illinois at Urbana-Champaign (1988)

[25.] Remacle JF, Li X, Shephard MS, et al.. Anisotropic adaptive simulation of transient flows using discontinuous galerkin methods. Int. J. Numer. Meth. Eng., 2005, 62(7): 899-923,

[26.] Schmidt, A., Siebert, K.G.: ALBERT: an adaptive hierarchical finite element toolbox. Albert-Ludwigs-Univ., Math. Fak. (2000)

[27.] Ullmann S, Rotkvic M, Lang J. POD-Galerkin reduced-order modeling with adaptive finite element snapshots. J. Comput. Phys., 2016, 325: 244-258,

[28.] Yano M. A minimum-residual mixed reduced basis method: Exact residual certification and simultaneous finite-element reduced-basis refinement. ESAIM: Math. Model. Numer. Anal., 2016, 50(1): 163-185,

[29.] Yano M. A reduced basis method for coercive equations with an exact solution certificate and spatio-parameter adaptivity: Energy-norm and output error bounds. SIAM J. Sci. Comput., 2018, 40(1): A388-A420,

Funding

Air Force Office of Scientific Research(FA9550-22-1-0004); National Science Foundation(2019287888)