请教大家一个问题

caohan

我想请教一下大家，这是我做的2D的翼型模拟，当我计算的24步的时候就开始分散了，
temperature limited to 1.000000e+000 in 2 cells on zone 2 in domain 1
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2 cells
  iter continuity x-velocity y-velocity     energy        nut     time/iter
    23 2.0765e+02 1.0028e-01 4.9721e-02 3.6678e-03 1.7781e-02  0:00:15   27
    24 2.8248e+03 9.3107e-02 9.2223e-02 2.7621e-03 2.4444e-02  0:00:12   26
Error: divergence detected in AMG solver
Error Object: ()
  iter continuity x-velocity y-velocity     energy        nut     time/iter
    24 2.8248e+03 9.3107e-02 9.2223e-02 2.7621e-03 2.4444e-02  0:00:22   50
absolute pressure limited to 1.000000e+000 in 5 cells on zone 2
    25 1.3491e+06 8.8138e-01 1.3820e-01 2.2625e-02 1.4079e-01  0:00:27   49
absolute pressure limited to 1.000000e+000 in 2103 cells on zone 2
Error: Floating point error: invalid number
Error Object: ()
这个是错误的信息
我先具体描述一下我的边界条件，我用的求解器是压力基，是压力远场
我计算了24步就不行了，但是我有不知道问题出在哪里，请大家帮帮忙，谢谢了
我顺便把网格的一些信息给出来
Grid Check
Domain Extents:
   x-coordinate: min (m) = -1.000000e+001, max (m) = 1.500000e+001
   y-coordinate: min (m) = -1.500000e+001, max (m) = 1.500000e+001
Volume statistics:
   minimum volume (m3): 2.305301e-006
   maximum volume (m3): 8.006855e-001
     total volume (m3): 6.531861e+002

Grid Quality:
Applying quality criteria for quadrilateral cells.
Maximum cell squish = 9.88732e-001
Maximum 'aspect_ratio' = 1.25705e+003

[ 本帖最后由 caohan 于 2010-2-9 23:18 编辑 ]

通流

我提几点吧。其实这也适用于很多情况。
（1）马赫数是什么？如果马赫数大于0.5，或者有局部马赫数大于0.8的话，pressure based 方法是不合适的。用density based方法会好的多。
（2）这网格也有点奇怪。叶片的上下独有一个很密的像后掠机翼的区域，这可有问题。发散很可能于这片区域有关。

onesupeng

哈哈，就网格我看着很不爽，好几个人了，前面号称很好网格的也是
我在想大家是如何对待网格生成的啊？

caohan

我模拟的雷诺数是在0.7，而且我想问一下就这翼型的网格怎样划分比较好，我是个新手，不太清楚怎样划分网格，高手们有没有什么书推荐的，谢谢大家乐

通流

我想你说的是马赫数为0.7吧。
如果是新手的话，一开始不要把网格搞得这么密。如果先算无粘的，那么有个几千个网格就够了。像这种跨音的流动，还是时间推进好。对网格的要求也比较低，很皮实。
算有粘的，主要是加密叶片附近的（边界层）的网格。有激波的话，那得先算出激波位置，然后再加密。如果有分离的话，那就很麻烦了。一般怎么算，都有些地方算不好。

caohan

不好意思，我说错了，我的马赫数是0.07，这个值应该不会产生激波吧。

通流

那就当不可压得算。基本思路还是可用。先用较粗的网格，再用密的网格。

caohan

我有个问题，我把网格调粗了，为什么我用simplic算就可以，用coupled算，算到第14步就分散了，算不了了，是我的boundary condition设的有问题还是我的网格不行，我的边界条件只设置了压力远场，另外我怕计算分散还把库朗数减小了，但是还是分散。 实在不知道哪里出的问题，我顺便把我的模型附上吧，哪位好心人帮我看看，小弟万分感谢，我的fluent是6.3版本，gambit是2.3的

通流

我试了一下，用density based， 显式格式，没问题。由于马赫数太小，收敛的有点慢。

caohan

可以用Coupled算吗？ 库朗数为多少比较好呢，另外我不太清楚pressure based与density based 哪个更适合做翼型绕流呢？

通流

不知道什么是pressure based coupled 算法。density based 就是时间推进法， 而pressure based应该是SIMPLE算法。我只知道时间推进法是对的，或者说时间推进法有完整的数学基础，也有很好的物理基础，所以很可靠。
我的观点是所有解欧拉方程，NS方程的人，都应该用时间推进法。

caohan

原帖由 通流 于 2010-2-11 19:36 发表

                               
登录/注册后可看大图

不知道什么是pressure based coupled 算法。density based 就是时间推进法， 而pressure based应该是SIMPLE算法。我只知道时间推进法是对的，或者说时间推进法有完整的数学基础，也有很好的物理基础，所以很可靠。
...

恩，好的，谢谢你

caohan

原帖由 通流 于 2010-2-11 19:36 发表

                               
登录/注册后可看大图

不知道什么是pressure based coupled 算法。density based 就是时间推进法， 而pressure based应该是SIMPLE算法。我只知道时间推进法是对的，或者说时间推进法有完整的数学基础，也有很好的物理基础，所以很可靠。
...

你好，我重算了一下，发现与风洞数据对比误差有点大， 风洞试验中7°攻角，马赫数为0.07状态下，cl=1.287, cd=0.01098, cm=-0.1314, 我算出来的为 cl=1, cd=0.024, cm=0.09, 这种误差算大吗？ 一般残差下降几个数量级可以看做为计算收敛呢？

梅羡仙

原帖由 通流 于 2010-2-11 11:36 发表

                               
登录/注册后可看大图

不知道什么是pressure based coupled 算法。density based 就是时间推进法， 而pressure based应该是SIMPLE算法。我只知道时间推进法是对的，或者说时间推进法有完整的数学基础，也有很好的物理基础，所以很可靠。
...

通流，你是不是认为这个世间的所有道理都只有符合了数学，才是科学的？
有一种神算法，比方说34乘以34，手指头一比划就能得出答案，不仅答案百试百灵，而且计算神速，它采用的方法完全不是现行的逻辑运算法则。你觉得这些又符合哪门子数学呢？

通流

这个与风洞对比的分析，你可以向白文先生请教。站长也很懂。

要是我的话，我会除了对上Re，M外，再看看风洞的尺寸。计算时把风洞的都放进去一起算。还有，这到底是二元的，还是三元的。你要知道，天下没有真正的二元流动。