[FLUENT] 高手请进，请吃饭！高超流场分析的问题

280852462

各位老大，小弟最近在做高超声速流场分析，用fluent迭代总出现如下错误
0.->1.->2.->3.->4.->5.<<<<<
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2508 cells
2 7.3503e-01 5.2409e-01 1.0000e+00 1.0000e+00 7.4436e-01 8.3665e-04 1.6680e-02
2:58:53
48
0.->1.->2.->3.->4.->5.<<<<<
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 5572 cells 3 6.1596e-01 3.2578e-01 1.0000e+00 1.0000e+00 6.2597e-01 1.7648e-03 1.6484e-02
2:31:52
47
0.
time step reduced in 229 cells->1.->2.->3.->4.->5.<<<<<
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 18160 cells
6 3.4213e-01 1.7965e-01 7.7446e-01 7.0366e-01 3.3066e-01 5.0451e-03 1.6611e-02
1:39:17
44
0.

time step reduced in 608 cells8 2.8194e-01 1.4804e-01 5.9028e-01 6.0275e-01 2.7358e-01 6.4990e-03 1.7032e-02
1:19:18
42
0.

time step reduced in 619 cells
absolute pressure limited to 1.000000e+000 in 1 cells on zone 2 9 2.6466e-01 1.3616e-01 5.2649e-01 5.5905e-01 2.5549e-01 7.0629e-03 1.7152e-02
1:12:02
41
0.
time step reduced in 614 cells
absolute pressure limited to 1.000000e+000 in 2 cells on zone 2 11 2.3156e-01 1.1894e-01 4.3127e-01 4.7472e-01 2.1349e-01 8.1036e-03 1.7272e-02
1:00:35
39
0.
time step reduced in 628 cells
absolute pressure limited to 1.000000e+000 in 10 cells on zone 2
temperature limited to 5.000000e+003 in 1 cells on zone 2 ->1.->2.->3.->4.->5.<<<<<

iter continuity x-velocity y-velocity z-velocity
energy
k
omega
time/iter

Error: Floating point error: invalid number
Error Object: ()这是怎么回事？另外，同样的模型和网格，我用0.8马赫算的时候就没这个问题，迭代1000多步没报错，把马赫数改成2或者5，其余条件不变，迭代几步就出现这个问题了，我试过增大计算域，改变粘性比，换湍流模型，都不行啊。
以前也有人出现过这个问题，零几年的帖子还有，好像都没怎么解决啊，这是个老问题了，希望真正的高手现现身吧，一堆人等着呢，好人有好报，先谢谢了。
小弟提几个问题，希望高手指点一下。
1，现在的模型长2M，建了个流场区域长120M，半径30M的圆柱，pressure-far-field边界条件，应该够大吧？
2，turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2508 cells
     time step reduced in 628 cells


absolute pressure limited to 1.000000e+000 in 10 cells on zone 2


temperature limited to 5.000000e+003 in 1 cells on zone 2
出现这种情况是什么原因？
有什么解决办法？
3，高马赫数的比热和热传导率不能再设成常数，应该怎么设置？
真心希望高手现身，以后如果有用得着小弟的地方，小弟一定全力以赴。

xrs333

用基于密度求解器。

280852462

我是用的基于密度求解器，还是有问题

xrs333

做turbulent viscosity 等值面，看看问题在那里。

280852462

谢谢前辈的耐心指导，现在画密了网格，可以往下算，但是提示这个，有影响吗？
99 3.7027e-02 3.3664e-02 4.9107e-02 5.3666e-02 3.5075e-02 4.1338e-05  0:01:02    1 0.
time step reduced in 3 cells
absolute pressure limited to 1.000000e+000 in 4 cells on zone 2

另外，turbulent viscosity 等值面怎么看问题在哪里？做了些迭代100步时的turbulent viscosity 等值面截图，麻烦前辈帮忙看看。

[ 本帖最后由 280852462 于 2011-4-24 17:51 编辑 ]

飞翔鸟

飞行器周围的网格太疏了，也没有边界层。

280852462

哦，谢谢前辈，另请教，我如果做6马赫飞行器，一般来说边界层第一层高度要多高？画几层比较好？
边界层单位和模型的单位一样吗？比如说我以mm为单位做了个模型，边界层为5e-5，导入fluent后scale，变成m单位，那5e-5是5e-5mm还是5e-5m呢？

飞翔鸟

我不知道你使用的什么网格工具，一般网格工具的长度是没有单位的，导入求解器后选择你的单位。
边界层厚度取决于你的湍流模型。高Re数模型在30-60之间结果比较准确，低Re数模型在1左右，并且一般要求粘性支层内10-15层网格。
很明显，你的网格非常粗，没有边界层。这种外流场的计算，你完全可以再飞行器周围布非常密的网格，往外逐渐放疏，既可以保证计算精度，也可以降低网格数。

[ 本帖最后由 飞翔鸟 于 2011-4-29 07:51 编辑 ]

onebyte

学习中……

41906957

我觉得这种不是很复杂的网格还是做结构的比较好，非结构可能精度有点低

长歌－废墟

absolute pressure limited to 1.000000e+000 in 4 cells on zone 2
这个楼主要注意一下，个人觉得，绝对压力不能这么低。

莫口

网格质量太差了，你这个结构比较简单，可以用结构网格，壁面处得画边界层啊！

浪子轻狂

网格太稀，你看你给出的云图，梯度变化非常难看。

return776

lz改进网格后 结果好了吗？