摘要: 对于城市大气环境模拟问题,高保真度往往意味着庞大的计算量,几何多尺度、高雷诺数等因素致使其求解非常具有挑战性,对计算机的硬件与算法都有极高的要求。寻求稳定、可扩展的高效并行数值算法以获取时效性强且可置信的数值结果是数值模拟的关键所在。基于国产大规模计算系统平台,本文通过对大规模非线性系统进行高效求解器和预处理技术的研究,提出一种可扩展并行Newton-Krylov-Schwarz算法。非线性方程采用非精确Newton方法进行求解,在每个Newton步,Jacobian系统通过基于区域分解方法的限制加性Schwarz 预
条件子处理,然后使用以GMRES方法为代表的Krylov子空间迭代法作为线性求解器进行求解。作为应用,我们对深圳地王大厦附近约一平方公里区域进行大气流场进行计算。数值结果显示,本文的算法在扩展至万核处理器平台时仍具有非常良好的可扩展并行性能,为未来对整个城市区域进行高精度(如网格分辨率达到米级,未知数个数达十亿量级)的大气环境模拟提供了一种可参考的算法。