[论文解读] Scalable smoothing strategies for a geometric multigrid method for the immersed boundary equations

Amneet Pal Singh Bhalla, Matthew G. Knepley+3 more|arXiv (Cornell University)|Dec 7, 2016

Lattice Boltzmann Simulation Studies参考文献 34被引用 5

一句话总结

本文针对浸入边界（IB）方程在斯托克斯流和纳维-斯托克斯流条件下的几何多重网格求解器，提出了一套可扩展、可并行化的平滑策略。引入了一种加法型Vanka平滑器和一种基于Schur补的平滑器，在低至中等雷诺数及高材料刚度条件下表现出优异性能，并在复杂流固耦合问题中展现出良好的可扩展性。

ABSTRACT

The immersed boundary (IB) method is a widely used approach to simulating fluid-structure interaction (FSI). Although explicit versions of the IB method can suffer from severe time step size restrictions, these methods remain popular because of their simplicity and generality. In prior work (Guy et al., Adv Comput Math, 2015), some of us developed a geometric multigrid preconditioner for a stable semi-implicit IB method under Stokes flow conditions; however, this solver methodology used a Vanka-type smoother that presented limited opportunities for parallelization. This work extends this Stokes-IB solver methodology by developing smoothing techniques that are suitable for parallel implementation. Specifically, we demonstrate that an additive version of the Vanka smoother can yield an effective multigrid preconditioner for the Stokes-IB equations, and we introduce an efficient Schur complement-based smoother that is also shown to be effective for the Stokes-IB equations. We investigate the performance of these solvers for a broad range of material stiffnesses, both for Stokes flows and flows at nonzero Reynolds numbers, and for thick and thin structural models. We show here that linear solver performance degrades with increasing Reynolds number and material stiffness, especially for thin interface cases. Nonetheless, the proposed approaches promise to yield effective solution algorithms, especially at lower Reynolds numbers and at modest-to-high elastic stiffnesses.

研究动机与目标

• 解决现有Vanka型平滑器在浸入边界（IB）方法几何多重网格求解器中并行可扩展性有限的问题。
• 为斯托克斯和纳维-斯托克斯条件下的半隐式IB格式，开发适用于并行高性能计算的新平滑策略。
• 提升在材料刚度和雷诺数变化多端的流固耦合问题中的求解器效率与可扩展性。
• 实现半隐式时间离散化IB方程所生成线性系统的鲁棒且可扩展的求解。
• 研究新平滑器在不同类型物理条件下的性能表现，包括薄壁与厚壁结构模型。

提出的方法

• 采用Vanka平滑器的加法变体以增强并行性，同时保持对斯托克斯-IB系统的鲁棒收敛性。
• 提出一种基于Schur补的平滑器，利用线性系统中的块结构以提升收敛速度。
• 在局部加密的结构化笛卡尔网格上，采用标准Galerkin限制与延拓算子实现几何多重网格。
• 对欧拉变量和拉格朗日变量均采用二阶有限差分进行空间离散化。
• 采用Peskin的四点正则化δ函数以实现精确且稳定的拉格朗日-欧拉耦合。
• 将平滑器集成至几何多重网格预条件子中，用于Krylov子空间求解器，从而实现鞍点系统的可扩展求解。

实验结果

研究问题

• RQ1 在几何多重网格框架中，加法型Vanka平滑器是否能有效应用于浸入边界方程并实现高效的并行化？
• RQ2 与传统平滑器相比，基于Schur补的平滑器是否在斯托克斯-IB系统中展现出更优的收敛特性？
• RQ3 在浸入边界格式中，随着雷诺数和材料刚度的增加，求解器性能如何退化？
• RQ4 所提出的平滑器在不同结构厚度和物理参数下的可扩展性与鲁棒性如何？
• RQ5 新平滑器在流固耦合模拟中是否能同时保持低雷诺数和中高雷诺数流动的效率？

主要发现

• 加法型Vanka平滑器为传统点对点Vanka平滑器提供了有效且可并行化的替代方案，显著提升了在分布式内存架构上的可扩展性。
• 基于Schur补的平滑器在高材料刚度和复杂界面几何问题中表现出强劲的收敛性能。
• 随着雷诺数和材料刚度的增加，线性求解器性能下降，尤其在薄界面情况下更为明显；然而，新平滑器在低至中等雷诺数范围内仍保持高效。
• 所提出的求解器在广泛的物理参数范围内（包括斯托克斯与非斯托克斯流动）均实现了可扩展性能。
• 平滑器对厚壁与薄壁结构模型均保持鲁棒性与高效性，其收敛速率在不同弹性刚度下保持稳定。
• 性能分析表明，RAS与RMS平滑器实现良好优化，适用于大规模高性能计算平台上的仿真。