[论文解读] On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems: Well-Posedness and Conservation

David A. Kopriva, Jan Nordström+1 more|arXiv (Cornell University)|May 10, 2021

Computational Fluid Dynamics and Aerodynamics参考文献 25被引用 5

一句话总结

本文通过能量法证明了双曲问题中重叠网格方法的适定性、守恒性和能量有界性，建立了其理论基础。该研究提出了一种新颖的罚函数方法，结合界面与内部重叠区域的罚项，即使在系数矩阵不可同时对角化的情况下，也能在多维情形下确保稳定性与原问题的等价性。

ABSTRACT

We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is ususally the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem.

研究动机与目标

• 建立双曲问题重叠网格方法的理论适定性，这是数值稳定性的先决条件。
• 解决现有文献中缺乏多维情形下稳定性证明的问题，特别是当系数矩阵不能同时对角化时。
• 构建一个框架，确保在多空间维情形下实现能量有界性、守恒性以及与原始单域问题解的等价性。
• 提出一种基于罚函数的耦合策略，使得无需对系数矩阵进行对角化即可实现适定性。
• 提供一种严格的基于能量的分析，考虑重叠区域的存在以及重叠区域中能量的重复计算。

提出的方法

• 使用能量法，推导一维和二维空间中重叠区域问题的能量估计。
• 在一维情形下应用特征界面条件，对标量方程和可对角化系统实现能量有界性。
• 在区域界面及重叠区域内引入多个内部罚函数项，以消除不被初始数据有界住的寄生项。
• 定义一个全局能量泛函，通过范数修正（T2）来校正重叠区域中重复计算的能量。
• 在重叠区域使用分部积分，关联内部与表面贡献（T1），从而实现对各个子区域的能量有界性估计。
• 施加特定的耦合矩阵条件（如式(128)），以确保罚项引入耗散，维持能量有界性与守恒性。

实验结果

研究问题

• RQ1 在使用特征界面条件时，双曲系统重叠区域问题在何种条件下是适定的？
• RQ2 当系数矩阵不能同时对角化时，为何能量法在多维情形下无法提供能量有界性？
• RQ3 是否可通过基于罚函数的方法，在无需矩阵对角化的情况下，恢复多维重叠问题中的能量有界性与适定性？
• RQ4 如何在具有重叠子区域的重叠区域问题中保证能量守恒？
• RQ5 耦合矩阵需满足何种条件，才能确保重叠问题的解与原始单域问题的解等价？

主要发现

• 在一维问题中采用特征耦合时，若系数矩阵可对角化，则重叠区域问题具有能量有界性且是适定的。
• 当系数矩阵不能同时对角化（如一般多维问题中），由于存在无界寄生项，能量法失效。
• 在界面及重叠区域内引入罚项后，若耦合矩阵满足条件(128)：(1−η)Σm_u = ηΣm_v，则可恢复能量有界性。
• 采用基于罚函数的耦合策略后，总能量被初始数据有界：E(T) ≤ ||ω₀||²_Ω，从而保证稳定性。
• 重叠问题的解与原始单域问题的解等价：在Ωu上有u = ω，在Ωv上有v = ω。
• 确保能量有界性的耦合条件同时保证了守恒性，即一个子区域损失的通量会被另一子区域所获得。