Toby van Gastelen, Wouter Edeling, Benjamin Sanderse
{"title":"Modeling Advection-Dominated Flows with Space-Local Reduced-Order Models","authors":"Toby van Gastelen, Wouter Edeling, Benjamin Sanderse","doi":"arxiv-2409.08793","DOIUrl":null,"url":null,"abstract":"Reduced-order models (ROMs) are often used to accelerate the simulation of\nlarge physical systems. However, traditional ROM techniques, such as those\nbased on proper orthogonal decomposition (POD), often struggle with\nadvection-dominated flows due to the slow decay of singular values. This\nresults in high computational costs and potential instabilities. This paper\nproposes a novel approach using space-local POD to address the challenges\narising from the slow singular value decay. Instead of global basis functions,\nour method employs local basis functions that are applied across the domain,\nanalogous to the finite element method. By dividing the domain into subdomains\nand applying a space-local POD within each subdomain, we achieve a\nrepresentation that is sparse and that generalizes better outside the training\nregime. This allows the use of a larger number of basis functions, without\nprohibitive computational costs. To ensure smoothness across subdomain\nboundaries, we introduce overlapping subdomains inspired by the partition of\nunity method. Our approach is validated through simulations of the 1D advection\nequation discretized using a central difference scheme. We demonstrate that\nusing our space-local approach we obtain a ROM that generalizes better to flow\nconditions which are not part of the training data. In addition, we show that\nthe constructed ROM inherits the energy conservation and non-linear stability\nproperties from the full-order model. Finally, we find that using a space-local\nROM allows for larger time steps.","PeriodicalId":501162,"journal":{"name":"arXiv - MATH - Numerical Analysis","volume":"117 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - MATH - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.08793","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Reduced-order models (ROMs) are often used to accelerate the simulation of
large physical systems. However, traditional ROM techniques, such as those
based on proper orthogonal decomposition (POD), often struggle with
advection-dominated flows due to the slow decay of singular values. This
results in high computational costs and potential instabilities. This paper
proposes a novel approach using space-local POD to address the challenges
arising from the slow singular value decay. Instead of global basis functions,
our method employs local basis functions that are applied across the domain,
analogous to the finite element method. By dividing the domain into subdomains
and applying a space-local POD within each subdomain, we achieve a
representation that is sparse and that generalizes better outside the training
regime. This allows the use of a larger number of basis functions, without
prohibitive computational costs. To ensure smoothness across subdomain
boundaries, we introduce overlapping subdomains inspired by the partition of
unity method. Our approach is validated through simulations of the 1D advection
equation discretized using a central difference scheme. We demonstrate that
using our space-local approach we obtain a ROM that generalizes better to flow
conditions which are not part of the training data. In addition, we show that
the constructed ROM inherits the energy conservation and non-linear stability
properties from the full-order model. Finally, we find that using a space-local
ROM allows for larger time steps.
降阶模型(ROM)通常用于加速大型物理系统的模拟。然而,传统的 ROM 技术(如基于适当正交分解(POD)的技术)由于奇异值衰减缓慢,在处理以对流为主的流动时往往力不从心。这导致了高计算成本和潜在的不稳定性。本文提出了一种使用空间局部 POD 的新方法,以解决奇异值衰减缓慢带来的挑战。我们的方法采用了局部基函数,而不是全局基函数,这种方法适用于整个域,类似于有限元方法。通过将域划分为子域,并在每个子域中应用空间局部 POD,我们实现了稀疏的表达,并能在训练区外更好地泛化。这样就可以使用更多的基函数,而不会产生过高的计算成本。为了确保子域边界的平滑性,我们引入了受unity 分区方法启发的重叠子域。我们通过模拟使用中心差分方案离散化的一维平流方程,验证了我们的方法。我们证明,利用我们的空间局部方法,我们得到的 ROM 能够更好地泛化不属于训练数据的流动条件。此外,我们还证明所构建的 ROM 继承了全阶模型的能量守恒和非线性稳定性特性。最后,我们发现使用空间局部 ROM 可以实现更大的时间步长。