昨天半夜在手机上看不到图,只看文字说明在箭体构形和湍流模型方面提了一些可能造成计算发散的原因。现在从网格图片看,计算发散的主要原因更可能在网格上:只做了多块网格处理,没看到边界层内的网格加细,而这部分处理对于计算的精度和稳定性至关重要。
对于定常、准定常的层流模型或RANS湍流模型,边界层厚度内需要10~100层网格刻画边界层流动的速度剖面,在激波区和边界层分离区还可能要进一步局部加细网格。如果对非定常问题(如跨音速湍流颤振)使用URANS(非定常RANS)或LES(大涡模拟)等模型,还要进一步细分网格达到相应的湍流旋涡尺度。
模拟圆柱形的箭体一般用壁面最密,从壁面向外随着直径增大逐渐变稀疏的网格形式,以兼顾边界层的网格密度和全流场的计算量。多块均匀网格一般只用在难以实现网格密度渐变的特殊情况下(如雷诺数不高的LBM)。
另外,升阻力系数Cd、Cl的单位是什么?典型的Cd数值应该在零点几,Cl更小。