卓越并没有注意到两人出去,他全部心神都集中在实验上。
看着电脑上的数据,他认真思索,许久后心道“应该采用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方程都解决了。
“万事俱备,只欠东风!”
说完他开始收拾东西,向四周看去,“咦,他们出去了啊!”
说完他将屋子里的机器关闭,扫视一圈没有任何遗漏后走出房间。
。