论文部分内容阅读
计算流体力学(Computational Fluid Dynamics,CFD)是研究、预测流体运动规律,揭示流动现象的重要方法之一,现已广泛应用于航空航天、气象海洋、武器装备、船舶制造等众多领域的科学研究和工程实践中。由于非结构网格相比于结构网格可以实现自动快速生成,非结构有限体积方法是目前求解CFD问题的常用离散方法。但非结构网格不具备结构网格的规则性、高效性等特征,有限体积方法也很难达到有限差分方法的高阶精度。梯度重构时有限体积方法的重要步骤之一,本文将有限差分结构化的特点应用到非结构有限体积方法的梯度重构中,为非结构网格单元生成局部方向,按照一定方向进行模板选择,使得梯度重构的过程可以利用更多沿着壁面法向方向的信息,从而提高梯度重构的计算速度、收敛速度、负载均衡等性能。本文的主要工作和创新点如下:
1)提出了基于阵面推进的非结构网格局部方向搜索算法(LDSMAF)。对二维三角形、四边形和混合非结构网格的局部方向进行了重定义,在此基础上基于阵面推进方法对重定义的局部方向进行求解。LDSMAF不仅考虑了网格拓扑关系,而且考虑了壁面法向信息。相比于曲线梯度重构的局部方向,LDSMAF的局部方向在壁面附近可以减少47%与壁面法向方向的角度偏移,且在整个计算域中正交性能最多可提高14%。
2)提出了基于局部方向的模板选择算法(SSMLD)及结构化最小二乘梯度重构算法(struLSQR)。基于LDSMAF获得的非结构网格局部方向,设计了模板生成算法SSMLD。SSMLD采用递归选择函数进行模板选择,并从两个方面对递归选择函数进行了优化。在此基础上,将SSMLD获得的模板与最小二乘梯度重构算法结合,提出了struLSQR梯度重构算法,设计和实现了并行化算法。struLSQR可使控制体获得可控且个数为常量的模板。在二维圆柱绕流算例中取模板个数为8时,struLSQR相比于扩展的最小二乘梯度重构算法(extended Least Squares,extLSQR),可节约35%的计算量,且计算和通信负载均衡性能可分别获得最多41%和36%的提高。
3)基于代码验证的要点以及梯度重构算法的特点,设计了测试函数对struLSQR的重构梯度进行了代码验证。代码验证过程测试并分析了单个网格单元的收敛精度分布、整体梯度误差及收敛精度以及不同网格类型的误差和收敛精度。对于四类四边形和五类三角形网格,struLSQR重构梯度的收敛精度均达0.9992以上,且在三角形网格中沿y方向的收敛精度大于其它三种梯度重构算法对应的收敛精度。
4)实现了较丰富的struLSQR应用并进行了深入分析。将struLSQR应用到了无粘、层流和湍流算例中,并与另两种最小二乘梯度重构算法进行性能对比与分析。测试的算例覆盖了无粘、层流和湍流流动,性能指标包括串行和并行指标,其中串行测试指标涵盖误差、收敛精度、计算速度、收敛速度以及流场特性等,并行执行指标涵盖进程执行时间、并行可扩展性、通信开销、理论与实测的负载均衡情况等,全面地反映了struLSQR在非结构有限体积应用中的优势。struLSQR总计算量减少为extLSQR总计算量的65%,计算性能和收敛速度为extLSQR的1.21倍和2.6倍,并可在粘性流动中获得更精确的法向速度型。在并行测试中,struLSQR的计算性能是extLSQR的1.27倍,并行通信性能提高为extLSQR的1.64倍。
1)提出了基于阵面推进的非结构网格局部方向搜索算法(LDSMAF)。对二维三角形、四边形和混合非结构网格的局部方向进行了重定义,在此基础上基于阵面推进方法对重定义的局部方向进行求解。LDSMAF不仅考虑了网格拓扑关系,而且考虑了壁面法向信息。相比于曲线梯度重构的局部方向,LDSMAF的局部方向在壁面附近可以减少47%与壁面法向方向的角度偏移,且在整个计算域中正交性能最多可提高14%。
2)提出了基于局部方向的模板选择算法(SSMLD)及结构化最小二乘梯度重构算法(struLSQR)。基于LDSMAF获得的非结构网格局部方向,设计了模板生成算法SSMLD。SSMLD采用递归选择函数进行模板选择,并从两个方面对递归选择函数进行了优化。在此基础上,将SSMLD获得的模板与最小二乘梯度重构算法结合,提出了struLSQR梯度重构算法,设计和实现了并行化算法。struLSQR可使控制体获得可控且个数为常量的模板。在二维圆柱绕流算例中取模板个数为8时,struLSQR相比于扩展的最小二乘梯度重构算法(extended Least Squares,extLSQR),可节约35%的计算量,且计算和通信负载均衡性能可分别获得最多41%和36%的提高。
3)基于代码验证的要点以及梯度重构算法的特点,设计了测试函数对struLSQR的重构梯度进行了代码验证。代码验证过程测试并分析了单个网格单元的收敛精度分布、整体梯度误差及收敛精度以及不同网格类型的误差和收敛精度。对于四类四边形和五类三角形网格,struLSQR重构梯度的收敛精度均达0.9992以上,且在三角形网格中沿y方向的收敛精度大于其它三种梯度重构算法对应的收敛精度。
4)实现了较丰富的struLSQR应用并进行了深入分析。将struLSQR应用到了无粘、层流和湍流算例中,并与另两种最小二乘梯度重构算法进行性能对比与分析。测试的算例覆盖了无粘、层流和湍流流动,性能指标包括串行和并行指标,其中串行测试指标涵盖误差、收敛精度、计算速度、收敛速度以及流场特性等,并行执行指标涵盖进程执行时间、并行可扩展性、通信开销、理论与实测的负载均衡情况等,全面地反映了struLSQR在非结构有限体积应用中的优势。struLSQR总计算量减少为extLSQR总计算量的65%,计算性能和收敛速度为extLSQR的1.21倍和2.6倍,并可在粘性流动中获得更精确的法向速度型。在并行测试中,struLSQR的计算性能是extLSQR的1.27倍,并行通信性能提高为extLSQR的1.64倍。