基于图剖分的多块结构网格负载平衡方法

负载平衡是影响并行计算性能的重要因素。针对多块结构网格,给出了一种改进的多层次图剖分负载平衡方法。该方法设计了新的网格剖分算法,采用改进的子块分裂方法与图剖分算法的循环调用实现结构对接网格剖分,并通过建立不同物体重叠网格间的连接关系,实现了结构重叠网格的负载平衡。采用2个典型算例对方法进行了对比验证,数值结果表明,子块分裂方法对剖分结果具有重要影响,采用循环调用算法及改进的子块分裂方法能有效地实现计算负载均衡及通信量优化,同时显著减少了网格块数及因虚网格导致的内存需求,有利于提高并行效率。该负载平衡方法与网格拓扑无关,适用于多块结构对接网格及重叠网格,且整体型剖分方式对于多块结构重叠网格具有更好的剖分效果。     

改进k-ω-γ转捩模式对不同雷诺数下HIAD的转捩预测

高超声速充气式柔性减速器（HIAD）在气动力作用下会变形为波纹状，从而促进流动转捩为湍流，准确预测其转捩位置和壁面热流对热防护系统的设计至关重要。拓展了分离诱导转捩预测性能的改进k-ω-γ模式，同时具备对第1模态、第2模态、横流模态以及流动分离失稳的预测能力。本文将其应用于不同雷诺数下壁面波纹变形的HIAD边界层转捩预测，并与原始k-ω-γ模式的预测结果进行了对比，以评估和验证其对复杂转捩现象的预测性能。在此基础上，细致剖析了改进k-ω-γ模式的转捩预测机制。结果表明，改进k-ω-γ转捩模式可准确预测不同来流雷诺数下HIAD的转捩起始位置、转捩阵面形态和壁面热流分布。波纹壁面波峰处的转捩预测主要由构造的分离间歇因子猝发。而在波谷位置，第1模态、横流模态以及流动分离的贡献都很重要。以上研究显示了改进k-ω-γ模式在复杂外形中的应用潜力，可为多重不稳定耦合作用下的转捩预测方法发展提供参考。     

新型三阶TVD限制器性能分析

在计算流体力学(CFD)方法中,限制技术是影响计算精度和计算稳定性的重要因素,目前应用较广的经典二阶总变差衰减(TVD)限制器虽能较好地满足计算要求,但性能差异大且分辨率和耗散的性能间并未得到良好权衡。对一种新型的三阶TVD插值限制器(T-3限制器)进行了研究并将其与3种经典限制器进行对比。首先通过一维黎曼问题,得出T-3限制器兼顾较高间断分辨率和良好稳定性的特点;接着通过高超声速双锥绕流和X-33外形飞行器的数值实验,得到T-3限制器具有刻画复杂流动的能力以及较优的气动热计算性能。     

雷诺应力与涡黏性模型的分离流预测对比分析

现代飞行器设计对流动分离的准确预测需求愈发迫切，但使用广泛的涡黏性模型预测结果却不尽如人意。雷诺应力模型凭借其扎实的理论基础有可能获得可信度更高的结果，但其性能优势仍需进一步评估与发掘。选取了SST模型、Stress-BSL(Baseline)模型分别作为涡黏性和雷诺应力模型的代表，对二维驼峰、二维跨声速凸块、跨声速三维ONERA M6机翼等算例进行了数值仿真。结果表明，相较于实验值，2种模型的预测结果均出现流动分离点提前，再附点滞后的现象，但Stress-BSL模型的预测误差更小，表现出了强逆压梯度下预测分离流动的优势。通过分析发现，2种模型对雷诺应力的低估导致了分离区较大。具体表现为SST模型引入的Bradshaw假设限制了湍动能的生成，使得模型计算的涡黏性系数偏小，强逆压梯度下低估边界层雷诺应力，导致流动分离提前。而分离区上缘处的雷诺应力预测偏小则被认为是流动再附滞后的主要原因。对于雷诺应力模型，误差主要来源于雷诺应力输运方程再分配项的模化不准。最后，针对上述原因，对SST模型关键封闭参数进行了重新标定，并进行了初步验证，结果表明修改后的模型预测表现好于原模型结果。