Model order reduction of HDG discretizations for harmonic wave propagation
Abstract
We propose and analyze a projection-based model reduction method for the rapid and reliable solution of time-harmonic acoustic wave propagation problems over a frequency band, with particular emphasis on seismic applications. The high-fidelity model is based on a hybridizable discontinuous Galerkin (HDG) discretization, formulated as a generalized saddle-point system coupling the volume and trace unknowns. The basis construction combines residual-based greedy sampling of the frequency parameter, enrichment with frequency derivatives of the HDG solution, and proper orthogonal decomposition (POD) compression. We establish the well posedness of both the high-fidelity and reduced problems and characterize their stability with inf-sup constants. We investigate the interaction between static condensation, inherent to the HDG discretization, and projection onto the reduced spaces, characterizing its implications for the reduced formulation and approximation errors. Since these operations do not commute in general, we derive a priori error estimates for both the trace and volume variables to quantify the gap between the two operations. We further identify a lifting condition on the volume reduced space that eliminates these gaps and yields a trace error estimate independent of the volume approximation error. Numerical experiments on the two-dimensional Marmousi model and a three-dimensional heterogeneous medium, with more than 10 7 unknowns, illustrate the theoretical estimates and show that the reduced model achieves speedups of several orders of magnitude compared to the high-fidelity solver.
Related
- Mixed Hybrid High-Order Methods for Sixth-Order Problems: Fully Non-conforming and C0-Conforming Discretizations
- Two-level additive Schwarz preconditioners with GDSW coarse spaces for high-order discretizations in 2D
- High-Order Fully Discrete Entropy-Stable DG Methods with Compact Explicit Single-Step Time Discretizations
- Two-level additive Schwarz preconditioners with GDSW coarse spaces for high-order discretizations in 3D
- Targeted Entropy-Stable Nonlinear Artificial Dissipation for High-Order Spectral-Element Summation-By-Parts Discretizations
BibTeX
@misc{barucq2026model,
title = {{Model order reduction of HDG discretizations for harmonic wave propagation}},
author = {Hélène Barucq and Florian Faucher and Ha Pham and Ishak Tifouti},
howpublished = {HAL},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.20145},
eprint = {2610.20145},
archivePrefix = {arXiv},
url = {https://inria.hal.science/hal-05805324v1},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Model order reduction of HDG discretizations for harmonic wave propagation
AU - Hélène Barucq
AU - Florian Faucher
AU - Ha Pham
AU - Ishak Tifouti
PY - 2026
DA - 2026/10/14
PB - HAL
DO - 10.48550/arxiv.2610.20145
UR - https://inria.hal.science/hal-05805324v1
AB - We propose and analyze a projection-based model reduction method for the rapid and reliable solution of time-harmonic acoustic wave propagation problems over a frequency band, with particular emphasis on seismic applications. The high-fidelity model is based on a hybridizable discontinuous Galerkin (HDG) discretization, formulated as a generalized saddle-point system coupling the volume and trace unknowns. The basis construction combines residual-based greedy sampling of the frequency parameter, enrichment with frequency derivatives of the HDG solution, and proper orthogonal decomposition (POD) compression. We establish the well posedness of both the high-fidelity and reduced problems and characterize their stability with inf-sup constants. We investigate the interaction between static condensation, inherent to the HDG discretization, and projection onto the reduced spaces, characterizing its implications for the reduced formulation and approximation errors. Since these operations do not commute in general, we derive a priori error estimates for both the trace and volume variables to quantify the gap between the two operations. We further identify a lifting condition on the volume reduced space that eliminates these gaps and yields a trace error estimate independent of the volume approximation error. Numerical experiments on the two-dimensional Marmousi model and a three-dimensional heterogeneous medium, with more than 10 7 unknowns, illustrate the theoretical estimates and show that the reduced model achieves speedups of several orders of magnitude compared to the high-fidelity solver.
ER -