代数涡粘模型在RANS%2fLES混合计算应用.pdfVIP

胡宁，徐晶磊等：代数涡粘模型在RANS／LES混合计算中的应用 代数涡粘模型在RANS／LES,昆合计算中的应用 胡宁1’2徐晶磊3 (1．中国航天空气动力技术研究院，北京市7201信箱56分箱，100074 2．北京大学湍流与复杂系统国家重点实验室，北京100871 3．北京航空航天大学，北京100083) 摘要：提出一种可以适用于RANS／LES混合计算的代数涡粘模型构造方案，在近壁处采用Prandtll约混合长模 型，在离壁面较远处恢复至lJSmagorinsky模型。新的模型在很粗的网格上得到了比DES更好的结果，消除了 对数区错位的问题。研究发现代数涡粘模型得到了更好的小尺度湍流结构，避免了结构退化。因此，新模 型从RANS向LES的过渡发生较早，得到的可解尺度雷诺应力较大，更好地保存了非定常湍流统计特性。 关键词：RANS／LES混合计算，槽道流 了足够多近壁涡的样本时，可以假设内层流动是由 0引言 RANS方程控制。距壁面足够远处，矛的值限制为 c脚△，而不是像原SA模型中那样线性增长。这 RANS／LES混合方法是针对有壁面存在的情况 r存-i2 而发展的湍流计算方法。在壁面附近，湍流涡结构 样，涡粘毁灭项cwif。l引的值增大，得到的涡粘系 L“J 的尺度很小，使用大涡模拟(LES)方法要求网格 数减小，使湍流脉动得以发展。 尺度以粘性尺度度量，导致计算量大大增加。根据 在存在大分离的流动中，由于分离的剪切层极 Piomelli等fl】的估算，LES需要的网格数随雷诺数 不稳定，任何扰动都会很快放大。因此DES方法 成Rel‘8的关系。由于目前的雷诺平均(RANS)湍 可以迅速产生小尺度涡结构，维持流动的湍流状态 流模型对附着流能够作出较好的预测，RANS／[，ES 及湍流统计量的正确预测。但在槽道湍流这样的内 混合方法在附着区利用较成熟的RANS模型降低流流动中，由于没有明显的附着流／分离流边界，导 计算成本，在分离区用LES方法提高计算精度。 致小尺度旋涡发展很慢【3】；另一方面，涡粘系数在 目前常用的RANS／LES混合方法多从原有的 距壁面达到一定距离后停止增长甚至开始下降，使 PANS模型发展而来。例如分离涡模拟(DES)【2】， 模型不能提供足够的雷诺应力，因此导致平均流场 其原型是一方程的Spalart．Allmarus(SA)模型： 的预测出现很大误差。 堡Dt=吉【v．“y+tT)vt7)‰(V盯】 严格地说，这种无分离或弱分离的流动模拟属 。’ 垧讥“阱 于带壁面模型的大涡模拟(WM。LES)，并不在DES 模型最初设计时所考虑的范畴之内【41。但是，DES 在原SA模型中，特征长度尺度0=d，其中d 在WM-LES中的失败严重影响了其作为一种湍流 为到壁面的距