沈泳星
上海交通大学
多场耦合问题交替迭代法收敛慢的原因剖析与可能对策
在多物理场耦合问题的求解策略中,交替迭代法应用广泛,它特别适用于各物理场需借助不同求解器进行求解的情况。在单场求解时,该方法通常能展现出良好的数值稳定性。然而,由于整体收敛往往需要在各物理场之间进行多次迭代,导致计算资源的大量消耗。
本文将首先深入阐释影响交替法收敛行为的关键因素,并定义一个名为耦合度(degree of coupling, DOC)的参数。该参数依据各物理场对应的雅可比矩阵(或称切刚度矩阵)计算得出。通过数学推导可知,交替法具有一阶收敛特性,各物理场的残差与误差折减因子大致都等于DOC。这表明,当DOC接近0时,交替法的收敛速度极快;当DOC小于1但接近1时,收敛速度会变得很慢;而当DOC大于等于1时,交替法甚至无法保证收敛。
以位移加载下的断裂相场法为例,研究发现 DOC 值会随着位移加载呈现单调递增的趋势,最终在略低于临界值1的区间内震荡。这一结论从理论层面解释了在相场断裂分析中,交替法虽能实现收敛,但位移场与相场之间迭代次数有时会居高不下的内在原因。
基于DOC理论框架,本文提出了一种高效的求解策略,旨在加快交替法的收敛速度。在断裂相场的基准算例以及热力耦合问题的实践中,该方法显著减少了迭代次数,缩短了运行时间。与现有的常用方法(如交替法、L-BFGS方法)相比,本文所提方法所需的迭代次数最少,运行时间也最短,展现出了明显的优势。
本科毕业于清华大学,博士毕业于美国斯坦福大学,曾在西班牙加泰罗尼亚理工大学担任讲师,2014年入选国家重大人才计划青年组,加入上海交通大学浦江国际学院(原密西根学院),2023年起任副院长。近年主要研究方向为断裂相场的算法与应用,部分发表在JMPS, CMAME, IJNME等计算力学权威期刊上。担任《Computer Modeling in Engineering & Sciences》及《力学季刊》编委,《计算力学学报》青年编委,曾承担多项国家自然基金项目。