POD-accelerated high-order multiscale method with corrected-basis reconstruction for transient multi-continuum problems
Abstract
Multi-continuum models provide an effective macroscopic description of transient flow in fractured and highly heterogeneous porous media by introducing several interacting state variables at each macroscale location. Although high-order multiscale (HOMS) methods reduce the cost of resolving microscopic heterogeneities while recovering local oscillatory features, repeated solution of the coupled homogenized system and repeated high-order reconstruction remain expensive in long-time and large-scale three-dimensional simulations. This work develops a proper orthogonal decomposition (POD) accelerated HOMS framework for transient multi-continuum problems. A vector-valued POD space is constructed from snapshots of the coupled homogenized solution, and the resulting block system is projected onto a low-dimensional space. Two recovery strategies are considered. The first reconstructs the reduced homogenized fields and performs full HOMS assembly at every output time. The second applies the spatial part of the multiscale recovery operator to retained POD modes offline and evaluates the online solution by direct linear combinations of corrected bases and time-derivative bases. An error estimate between the fully discrete homogenized solution and its POD-reduced approximation is established for the coupled system. A recovered POD-HOMS estimate is then obtained by combining the HOMS consistency error with POD truncation errors propagated through value, gradient, inter-continuum exchange, second-derivative, and time-derivative recovery terms. Numerical experiments for two- and three-dimensional heterogeneous configurations show that the reduced solutions retain the accuracy and local oscillatory behavior of the corresponding HOMS solution while substantially decreasing computational cost.
Related
- A macroscopic-shadow-corrected lattice Boltzmann method for fast, time-accurate simulation of low-Reynolds-number transient flows
- CLEAR: Complex Learned Explicit Analytical Regularization for Ultra-Accelerated 4D Flow CMR Reconstruction
- Accelerated Plug-and-Play Davis-Yin Splitting for Nonconvex Image Reconstruction
- Corrected L1 Scheme for the Fractional Subdiffusion Equation
- A multi-class kinetic traffic flow model: discrete-velocity formulation and diffusively-corrected macroscopic limits
BibTeX
@misc{peng2026podaccelerated,
title = {{POD-accelerated high-order multiscale method with corrected-basis reconstruction for transient multi-continuum problems}},
author = {Jiayuan Peng and Hao Dong and Jiale Linghu},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.17758},
eprint = {2610.17758},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.17758},
}
RIS
TY - UNPB
TI - POD-accelerated high-order multiscale method with corrected-basis reconstruction for transient multi-continuum problems
AU - Jiayuan Peng
AU - Hao Dong
AU - Jiale Linghu
PY - 2026
DA - 2026/10/13
PB - arXiv
DO - 10.48550/arxiv.2610.17758
UR - https://arxiv.org/abs/2610.17758
AB - Multi-continuum models provide an effective macroscopic description of transient flow in fractured and highly heterogeneous porous media by introducing several interacting state variables at each macroscale location. Although high-order multiscale (HOMS) methods reduce the cost of resolving microscopic heterogeneities while recovering local oscillatory features, repeated solution of the coupled homogenized system and repeated high-order reconstruction remain expensive in long-time and large-scale three-dimensional simulations. This work develops a proper orthogonal decomposition (POD) accelerated HOMS framework for transient multi-continuum problems. A vector-valued POD space is constructed from snapshots of the coupled homogenized solution, and the resulting block system is projected onto a low-dimensional space. Two recovery strategies are considered. The first reconstructs the reduced homogenized fields and performs full HOMS assembly at every output time. The second applies the spatial part of the multiscale recovery operator to retained POD modes offline and evaluates the online solution by direct linear combinations of corrected bases and time-derivative bases. An error estimate between the fully discrete homogenized solution and its POD-reduced approximation is established for the coupled system. A recovered POD-HOMS estimate is then obtained by combining the HOMS consistency error with POD truncation errors propagated through value, gradient, inter-continuum exchange, second-derivative, and time-derivative recovery terms. Numerical experiments for two- and three-dimensional heterogeneous configurations show that the reduced solutions retain the accuracy and local oscillatory behavior of the corresponding HOMS solution while substantially decreasing computational cost.
ER -