求解三维流固耦合问题的一种全隐全耦合区域分解并行算法
三维流固耦合问题的非结构网格数值算法在很多工程领域都有重要应用,目前现有的数值方法主要基于分区算法,即流体和固体区域分别进行求解,因此存在收敛速度较慢以及附加质量导致的稳定性问题,此外,该类算法的并行可扩展性不高,在大规模应用计算方面也受到一定限制.本文针对三维非定常流固耦合问题,提出一种基于区域分解的全隐全耦合可扩展并行算法.首先基于任意拉格朗日-欧拉框架建立流固耦合控制方程,然后时间方向采用二阶向后差分隐式格式、空间方向采用非结构稳定化有限元方法进行离散.对于大规模非线性离散系统,构造一种结合非精确Newton法、Krylov子空间迭代法与区域分解Schwarz预条件子的Newton-Krylov-Schwarz(NKS)并行求解算法,实现流体、固体和动网格方程的一次性整体求解.采用弹性障碍物绕流的标准测试算例对数值方法的准确性进行了验证,数值性能测试结果显示本文构造的全隐全耦合算法具有良好的稳定性,在不同的物理参数下具有良好的鲁棒性,在"天河二号"超级计算机上,当并行规模从192增加到3072个处理器核时获得了91％的并行效率.性能测试结果表明本文构造的NKS算法有望应用于复杂区域流固耦合问题的大规模数值模拟研究中.
流固耦合、有限元、全耦合算法、区域分解法、非结构网格
54
O357.1(流体力学)
国家自然科学基金;广东省自然科学基金资助项目;广东省自然科学基金资助项目
2023-01-13（万方平台首次上网日期，不代表论文的发表时间）
共11页
3513-3523
