求教一个问题，关于大涡模拟时间步长的问题

鹊桥仙客

我想问一下同行的高手们，有没有用大涡模拟计算二维NACA0012翼型小攻角状态下的流动状况：MA=0.799，RE=9.0*10^6，攻角为2.26；以及RAE2822，MA=0.75，RE=6.2*10^6，攻角为2.81这个状态的。另外，计算这个状态，用大涡模拟计算的话，能否用当地时间步长进行时间推进。我用的是三阶龙格库塔法求解。但查找了好多文献，基本都是采用定时间步长进行时间推进求解大分离流动。对于小攻角绕翼型绕流这类非定常特性不明显的问题，能否用当地时间步长进行大涡模拟求解以加速收敛。雷诺平均方法即是用当地时间步长加速收敛，其算出的是平均量，但是大涡模拟控制方程中变量是瞬时的，用当地时间步长貌似没有人这样算过。我用LILY模型，采用当地时间步长可以算的跟实验值吻合较好，但是WALE模型却不行，在激波后面出现了压力震荡和一系列的涡。希望高手们予以解答。谢谢啦。