Gradient-preserving hyper-reduction of nonlinear dynamical systems via discrete empirical interpolation

by   Cecilia Pagliantini, et al.

This work proposes a hyper-reduction method for nonlinear parametric dynamical systems characterized by gradient fields such as (port-)Hamiltonian systems and gradient flows. The gradient structure is associated with conservation of invariants or with dissipation and hence plays a crucial role in the description of the physical properties of the system. Traditional hyper-reduction of nonlinear gradient fields yields efficient approximations that, however, lack the gradient structure. We focus on Hamiltonian gradients and we propose to first decompose the nonlinear part of the Hamiltonian, mapped into a suitable reduced space, into the sum of d terms, each characterized by a sparse dependence on the system state. Then, the hyper-reduced approximation is obtained via discrete empirical interpolation (DEIM) of the Jacobian of the derived d-valued nonlinear function. The resulting hyper-reduced model retains the gradient structure and its computationally complexity is independent of the size of the full model. Moreover, a priori error estimates show that the hyper-reduced model converges to the reduced model and the Hamiltonian is asymptotically preserved. Whenever the nonlinear Hamiltonian gradient is not globally reducible, i.e. its evolution requires high-dimensional DEIM approximation spaces, an adaptive strategy is performed. This consists in updating the hyper-reduced Hamiltonian via a low-rank correction of the DEIM basis. Numerical tests demonstrate the applicability of the proposed approach to general nonlinear operators and runtime speedups compared to the full and the reduced models.


page 1

page 2

page 3

page 4


Rank-adaptive structure-preserving reduced basis methods for Hamiltonian systems

This work proposes an adaptive structure-preserving model order reductio...

Adaptive symplectic model order reduction of parametric particle-based Vlasov-Poisson equation

High-resolution simulations of particle-based kinetic plasma models typi...

Reduced models with nonlinear approximations of latent dynamics for model premixed flame problems

Efficiently reducing models of chemically reacting flows is often challe...

Dynamical reduced basis methods for Hamiltonian systems

We consider model order reduction of parameterized Hamiltonian systems d...

Structure-preserving model order reduction of Hamiltonian systems

We discuss the recent developments of projection-based model order reduc...

Model order reduction strategies for weakly dispersive waves

We focus on the numerical modelling of water waves by means of depth ave...

S-OPT: A Points Selection Algorithm for Hyper-Reduction in Reduced Order Models

While projection-based reduced order models can reduce the dimension of ...