Abstract:This article is concerned with quadrature/cubature rules able to deal with multiple subspaces of functions, in such a way that the integration points are common for all the subspaces, yet the (nonnegative) weights are tailored to each specific subspace. These subspace‐adaptive weights cubature rules can be used to accelerate computational mechanics applications requiring efficiently evaluating spatial integrals whose integrand function dynamically switches between multiple pre‐computed subspaces. One of such applications is local hyperreduced‐order modeling (HROM), in which the solution manifold is approximately represented as a collection of basis matrices, each basis matrix corresponding to a different region in parameter space. The proposed optimization framework is discrete in terms of the location of the integration points, in the sense that such points are selected among the Gauss points of a given finite element mesh, and the target subspaces of functions are represented by orthogonal basis matrices constructed from the values of the functions at such Gauss points, using the singular value decomposition (SVD). This discrete framework allows us to treat also problems in which the integrals are approximated as a weighted sum of the contribution of each finite element, as in the energy‐conserving sampling and weighting method of C. Farhat and co‐workers. Two distinct solution strategies are examined. The first one is a greedy strategy based on an enhanced version of the empirical cubature method (ECM) developed by the authors elsewhere (we call it the subspace‐adaptive weights ECM, SAW‐ECM for short), while the second method is based on a convexification of the cubature problem so that it can be addressed by linear programming algorithms. We show in a toy problem involving integration of polynomial functions that the SAW‐ECM clearly outperforms the other method both in terms of computational cost and optimality. On the other hand, we illustrate the performance of the SAW‐ECM in the construction of a local HROMs in a highly nonlinear equilibrium problem (large strains regime). We demonstrate that, provided that the subspace‐transition errors are negligible, the error associated to hyperreduction using adaptive weights can be controlled by the truncation tolerances of the SVDs used for determining the basis matrices. We also show that the number of integration points decreases notably as the number of subspaces increases, and that, in the limiting case of using as many subspaces as snapshots, the SAW‐ECM delivers rules with a number of integration points only dependent on the intrinsic dimensionality of the solution manifold and the degree of overlapping required to avoid subspace‐transition errors. The Python source codes of the proposed SAW‐ECM are openly accessible in the public repository https://github.com/Rbravo555/localECM.

Kernel-based active subspaces with application to computational fluid dynamics parametric problems using the discontinuous Galerkin method

Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems

Model Order Reduction by means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics

Friedrichs' systems discretized with the Discontinuous Galerkin method: domain decomposable model order reduction and Graph Neural Networks approximating vanishing viscosity solutions

Comparative study of inner-outer Krylov solvers for linear systems in structured and high-order unstructured CFD problems

Active-Set Reduced-Space Methods with Nonlinear Elimination for Two-Phase Flow Problems in Porous Media.

Multi-fidelity data fusion through parameter space reduction with applications to automotive engineering

Multi-Level Solution Strategies for Implicit Time Discontinuous Galerkin Discretizations

An efficient shape parametrisation by free-form deformation enhanced by active subspace for hull hydrodynamic ship design problems in open source environment

Accelerating high order discontinuous Galerkin solvers through a clustering-based viscous/turbulent-inviscid domain decomposition

Discovering Artificial Viscosity Models for Discontinuous Galerkin Approximation of Conservation Laws using Physics-Informed Machine Learning

Two-level dynamic load-balanced p-adaptive discontinuous Galerkin methods for compressible CFD simulations

Enhancing Low-Order Discontinuous Galerkin Methods with Neural Ordinary Differential Equations for Compressible Navier--Stokes Equations

Two Numerical Schemes Via an Auxiliary Function for Nonlocal Problems With Integrable Kernels

KFVM-WENO: A high-order accurate kernel-based finite volume method for compressible hydrodynamics

Kernel Multi-Grid on Manifolds

Discontinuous Galerkin schemes for hyperbolic systems in non-conservative variables: quasi-conservative formulation with subcell finite volume corrections

An optimisation-based domain-decomposition reduced order model for parameter-dependent non-stationary fluid dynamics problems

A subspace-adaptive weights cubature method with application to the local hyperreduction of parameterized finite element models

Augmented Flexible Krylov Subspace methods with applications to Bayesian inverse problems

A subspace‐adaptive weights cubature method with application to the local hyperreduction of parameterized finite element models