请高手帮忙看看这个case怎么收敛

hyfair

我自己变成算出一道明显的激波打在地板上
可是在fluent里却不停震荡
请高手帮忙看看怎么才能让他收敛

cjjj2000

我看了看，觉得你的网格问题太大。
疏密不合理，在边界附近你可以做边界层，但是在远些的地方你从疏到密靠近地板是不对的
你采用全部一致的四方网格算得肯定比你现在好。
难道你自己编程采用的也是这种网格？

hyfair

果然是网格的问题，谢谢

hyfair

好像还是有点问题，当我把收敛判据减小两个数量级
残值仍然会出现震荡，这种震荡的出现好像和初始条件有关，在我解偶啦方程的时候通过一个耗散项再叠带过程中销去了，网格的疏密关系不大，
但是在fluent里我已经把网格取得很秘了，但是仍然无法收敛到10e-7数量级

cjjj2000

理论上讲，不要做边界层，且
你将所有网格都做成一样大小，网格足够密，是可以收敛到你指定的
任何数量级的

hyfair

加边界层是因为原来算亚音速问题的
这次的网格取得很密了，残值震荡看不到了，不过好像不收敛了？到了300步以后残值根本不下去

cjjj2000

加边界层是因为原来算亚音速问题的?
这个很经典的激波算例是无粘的！所以我建议你不要做边界层的
但你却使用了湍流模型！

hyfair

感谢cjjj2000对这个问题的关注，这个问题的提出本来是很随便的，因为我自己编的解瓯拉方程的程序能很快地得到结果
开始我是算一个台阶的亚音速流动的
后来把这个算例拿来看看fluent算的结果
原来以为会很容易就算处我要的激波，就用了原来的网格，忘了取消超音速问题的边界层，
至于粘性，那是因为我看有残值震荡，以为把它加上去会起到耗散项的作用
但是从结果来看，粘性对计算并没有什么影响
这次的模型改回去，残值震荡在300步以后就很明显了
我个人认为可能初始条件或者边界条件给得不好，但是不知道该怎么改
也许我是错的，现在对这个问题也没有任何办法
大家一起研究研究，毕竟发现问题就意味着进步
Version: 2d, segregated (2d, segregated, inviscid)
Release: 6.1.18
Title:
Models
------
   Model                        Settings   
   -------------------------------------
   Space                        2D         
   Time                         Steady     
   Viscous                      Inviscid   
   Heat Transfer                Enabled   
   Solidification and Melting   Disabled   
   Radiation                    None      
   Species Transport            Disabled   
   Coupled Dispersed Phase      Disabled   
   Pollutants                   Disabled   
   Soot                         Disabled