第62章 湍流模型校准实验(2 / 2)

卓越并没有注意到两人出去,他全部心神都集中在实验上。

看着电脑上的数据,他认真思索,许久后心道“应该采用jason中心有限体积+人工粘性方法,对守恒形式的n-s控制方程进行空间离散。”

“高reynolds数流动,边界层厚度很薄,粘性影响主要集中在靠近物面很近的区域中,且对于粘性计算所采用的物面王哥为大展弦比网格,粘性作用主要集中在垂直于物面方向上。”

“所以,在计算粘性作用时要采用薄层假设,只考虑了垂直于物面方向上的粘性作用。”

“另外,由于中心格式本身不具有耗散性,所以需要引入人工粘性想以有效抑制数值震荡,加速收敛到定常状态,得到光滑的压力分布。”

“所以,离散后的控制方程可以这样写。”

【v?,?,?dw?,?,?/dt+?,?,?-?,?,?-d?,?,?=0】

“jason的人工粘性是二阶和四阶差分的组合,即。”

【d?,?,?=d??2?+d??2?+d??2?-d????-d????-d????)w?,?,?

其中d??2?w?,?,?=??(λ?+?/?,?,?e?2??+?/?,?,?)△?w?,?,?

d????w?,?,?=??(λ?+?/?,?,?e?2??+?/?,?,?)△???△?w?,?,?

……】

卓越打开c型网格,网格点数为161x45。

远场取在12倍的弦长以外,从物面向外的第一层网格厚度为20x10(-5)倍弦长。

只见屏幕上出现一张密密麻麻的网格,网格所有线条成半圆形,中间有一条平均线将它们分为两半,中心还有一个如水滴一般的白色空隙。

“最终得出,流场计算状态。”

很快,电脑上就显示出经过公式计算的结果。

【∞=075,α=319?,re=63x10?】

他又是打开两个压强分布图,并让s-a模型与b-l模型的计算结果和实验结果进行对比。

只见屏幕上出现五个不同的等压线,线条上有许多点,显示着每一个点的具体数据。

“s-a模型对分离的模拟能力,应该选取较大迎角状态。”

【∞=0844,α=506?,re=117x10?】

看着屏幕上的曲折线,他心道“从这可以看出,机翼上表面翼根部的两道激波在靠近翼梢部分逐渐汇合,所以,我推导出的湍流模型方程是对的。”

“一次不行,最起码再计算三次!”

说完他又开始重复计算,连续计算了五次。

许久之后,他心道“经过五次计算,我可以肯定方程是对的。”

说完脸上露出兴奋的笑容。

边界层方程、湍流模型方程、湍流摩擦系数公式、雷诺数、紊动机理和n-s方程都解决了。

“万事俱备,只欠东风!”

说完他开始收拾东西,向四周看去,“咦,他们出去了啊!”

说完他将屋子里的机器关闭,扫视一圈没有任何遗漏后走出房间。