飞行棋狂魔 发表于 2016-5-22 15:09:59

日常学习练手-10Ma二维锥体后加整流锥

一晚上算出的15Ma云图震楼,同一网格,以10马赫流场作为初值,迭代40000次。
有兴趣的可以对比10Ma的结果 看马赫角 激波和膨胀波 尾部回流区区别

飞行棋狂魔 发表于 2016-5-22 15:36:07

出于懒的考虑 整流锥跟前部锥体一样把
特征长度变成1.4m 第一层网格高度经cfd-onliney+计算器计算没变
网格增长值从1.2减少到1.1 共80层
几何模型和网格如下

飞行棋狂魔 发表于 2016-5-22 15:40:43

fluent设置均同上一篇 库朗数0.5
迭代开始(人可以休息 电脑不能╮(╯3╰)╭)

飞行棋狂魔 发表于 2016-5-22 19:23:38

一万三千次 还没收敛。。。
网格变密似乎库朗数可以加大,残差越小计算速度越慢(软件用得越多困惑越多23333333)

飞行棋狂魔 发表于 2016-5-22 20:22:03

终于收敛了。。。 电脑发烫 两万次
跟上次对比收敛过程: 1 网格密使库朗数可以继续增大
                                 2 精度有较大提高,可看到除能量曲线其他的都极低

飞行棋狂魔 发表于 2016-5-22 20:39:23

0.0 云图如下

飞行棋狂魔 发表于 2016-5-22 20:41:13

后缘温度依旧最高,我决定用matlab计算一下理论驻点静温应该是多少

飞行棋狂魔 发表于 2016-5-22 21:05:21

理论估值跟结果基本一致,结果偏大大概是因为壁面采用的是绝热壁面,绝热壁面估算前缘驻点基本都会略偏大

飞行棋狂魔 发表于 2016-5-22 21:09:10

因为温度高于2500度会发生电离,激发了空气中各分子的振动自由度,并会有复杂的化学反应。所以10马赫的这个结论仅供娱乐,但是后缘为什么依旧温度这么高 激波与边界层摩擦么?????233333

飞行棋狂魔 发表于 2016-5-22 21:30:18

加整流锥后的阻力是119.84 N,其中压力53.65N,粘性阻力65.69N
上一贴的       阻力是116.34 N,其中压力60.02N,粘性阻力56.32N
阻力竟然增大了。。。这不科学 不知道原因。
以下是瞎掰: 1 压力阻力减小,说明整流锥起作用了。
                     2 粘性阻力增大,因为增大较多,我猜是我把整流锥画太大了,导致粘性力增大明显。。。。。

                     3 压力阻力不知道是不是压差阻力
页: [1] 2
查看完整版本: 日常学习练手-10Ma二维锥体后加整流锥