高雷诺数湍流风场大涡模拟的并行直接求解方法1)

摘要： 在环境流体力学中,风场是风沙流、风雪流等自然环境特性问题研究的动力源和基础. 通常采用壁湍流模型进行风场大涡模拟(large eddy simulation, LES)计算,但受到计算规模的限制使得 高雷诺数风场的模拟计算难以实现. 并行计算技术是解决大规模高雷诺数风场大涡模拟的关键技术之一. 在不可压湍流风场的LES模拟中,压力泊松方程的并行计算技术是进行规模并行计算的困难点. 根据风场流动模拟计算的特点,采用水平网格等距而垂直于地面网格非等距,在解决规模并行计算中求解压力泊松方程的难点问题时,利用FFT解耦三维泊松方程使其变为垂向的一维三对角方程, 并利用可并行的三对角方程PDD求解技术,可建立三维泊松方程的直接并行求解技术. 结合其它容易并行的动量方程计算,本文建立风场LES模拟的并行直接求解方法(parallel direct method-LES, PDM-LES). 在超级计算机上对新方法进行并行效率测试,并行计算效率达到90${\%}$. 新的方法可用于进行湍流风场大涡模拟的大规模并行计算. 计算结果表明,湍流风场瞬时速度分布近壁面存在条带状的拟序结构,平均场的速度分布符合速度对数律特性,风场湍流特性基本合理