[论文解读] Wave propagation across acoustic / Biot's media: a finite-difference method

Guillaume Chiavassa, Bruno Lombard|arXiv (Cornell University)|Sep 15, 2011

Lattice Boltzmann Simulation Studies参考文献 40被引用 30

一句话总结

本论文提出了一种高阶有限差分方法，用于基于Biot方程模拟流体与多孔弹性介质之间的波传播，包含开放、密封及非理想孔隙界面条件。通过结合四阶ADER时间积分、时空网格局部加密以及具有亚格子分辨率的浸入界面法，该方法能够精确捕捉扩散型慢纵波，并实现高精度，其结果经由所有类型界面的解析解验证。

ABSTRACT

Numerical methods are developed to simulate the wave propagation in heterogeneous 2D fluid / poroelastic media. Wave propagation is described by the usual acoustics equations (in the fluid medium) and by the low-frequency Biot's equations (in the porous medium). Interface conditions are introduced to model various hydraulic contacts between the two media: open pores, sealed pores, and imperfect pores. Well-possedness of the initial-boundary value problem is proven. Cartesian grid numerical methods previously developed in porous heterogeneous media are adapted to the present context: a fourth-order ADER scheme with Strang splitting for time-marching; a space-time mesh-refinement to capture the slow compressional wave predicted by Biot's theory; and an immersed interface method to discretize the interface conditions and to introduce a subcell resolution. Numerical experiments and comparisons with exact solutions are proposed for the three types of interface conditions, demonstrating the accuracy of the approach.

研究动机与目标

• 开发一种鲁棒的数值方法，用于模拟由声学方程与Biot方程控制的异质流体/多孔介质中的波传播。
• 准确模拟三种类型的水力界面条件：开放孔隙、密封孔隙以及具有有限渗透率的非理想孔隙。
• 解决模拟慢纵波的挑战，该波具有扩散特性，且在界面附近需要精细的空间分辨率。
• 确保流体-多孔弹性耦合系统的初边值问题具有适定性。
• 在不使用非结构化网格的结构化笛卡尔网格上，实现复杂几何形状下的高阶精度与稳定性。

提出的方法

• 采用带有Strang分裂的四阶ADER格式进行时间积分，确保最优CFL稳定性条件。
• 该方法对流体（声学）与多孔弹性（Biot方程）区域分别使用独立求解器，并通过浸入界面法进行耦合。
• 利用亚格子分辨率技术对水力接触界面条件（开放、密封、非理想孔隙）进行离散化，以处理非对齐界面。
• 在界面附近采用时空网格局部加密策略，以解析局部化且具有扩散特性的慢纵波。
• 采用基于泰勒展开的重构方法结合奇异值分解，计算界面点处的修正有限差分模板，确保一致性和精度。
• 在计算中使用归一化物理参数以改善条件数，数值实验中采用归一化因子 𝒩 = 1000。

实验结果

研究问题

• RQ1 如何在具有复杂界面条件的流体与多孔弹性介质中实现波传播模拟的高阶精度？
• RQ2 不同类型的水力接触（开放孔隙、密封孔隙、非理想孔隙）对波的透射与反射有何影响？
• RQ3 基于笛卡尔网格的有限差分方法能否在不使用非结构化网格的情况下准确解析Biot理论预测的慢纵波？
• RQ4 ADER时间积分、时空网格局部加密与浸入界面法的结合如何提升稳定性与精度？
• RQ5 所提出的方法在任意形状界面与不同材料属性下，其精度能保持到何种程度？

主要发现

• 该方法在模拟流体与多孔弹性介质之间的波传播时表现出高精度，数值结果与所有三种界面条件的解析解高度吻合。
• 浸入界面法成功以亚格子分辨率捕捉界面条件，实现了对非对齐界面的精确处理。
• 时空网格局部加密有效解析了局部化且具有扩散特性的慢纵波。
• 采用Strang分裂的四阶ADER格式在刚性系统（高粘度）下仍能保持最优CFL条件下的稳定性。
• 使用归一化变量改善了数值条件数，使得采用真实物理参数的稳定模拟成为可能。
• 初边值问题的适定性得到严格证明，确保了模型的数学一致性。