|
发表于 2013-7-16 12:58:35
|
显示全部楼层
粘性流计算有两大类:一种是传统的“无粘流计算+粘性修正”的方式,即先用无粘流计算程序(基于位势流或欧拉方程的求解器)进行计算,再在壁面附近加入粘性修正,即将边界层位移厚度加进去,然后二者迭代几次得到最终的解。第二种是现在最流行的基于RANS方程求解的方式,即用雷诺平均NS方程(RANS)作控制方程对流场进行求解,湍流粘度则由湍流模型方程确定。通常整个求解过程可以分为两大步,第一步是RANS方程的求解,第二步是湍流方程的求解。湍流方程求解后得到的湍流变量(比如k-e模型中的湍动能k和耗散率e),再用湍流变量计算得到湍流粘度。在边界层的处理上,一种是要求最底层网格能够分辨出湍流边界层的粘性底层,也就是要求y+的值小于等于1,还有一种是用壁面函数将粘性底层的特性进一步模式化,即用一个经验公式给出对底层的估算,从而可以放宽对底层网格高度的要求,比如可以放宽到y+大于30至200。大致情况就是这样。需要提醒的一点是,边界层是高雷诺数流场的一个特征,即在雷诺数很高时,流场中存在这样一个粘性集中区。在雷诺数不高时,边界层可能并不存在,粘性的影响是遍及整个流场的。所以,实际计算中还是要先把流场的特点研究一下,再决定用哪种方式进行计算。 |
|