[菜鳥求救] K-omega SST 头几秒就崩溃了



  • @东岳[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    CFD发散并没有一阵见血的解决办法。目前很难判断哦,但我觉得你一定会搞定的 :ok2:

    老师你太高估我了 我只是个新手渣渣 光凭我自己一个始终有点吃力:mihu:
    其实我在怀疑这会不会是k和omega的初始值所造成的问题,毕竟这都是估算的,没有实验根据,而且听说k omega对初始值很敏感。请问老师有没有方法可以验证初始值是否合适?



  • @jinjolee k和omega可以用fluent里面初始化的值做参考


  • CORE 网格教授 OpenFOAM教授 管理员

    @jinjolee[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    听说k omega对初始值很敏感

    是的,kOmega模型是这样的。不过Omega可以从epsilon计算出来:$\Omega=\frac{\varepsilon}{0.09 k} $

    你这个我感觉完全是边界条件设置的问题,因为不涉及到算法,你会会搞定的 :duang:



  • simpleFoam是稳态求解器,你这个问题是不是用瞬态比较好?如果你模拟的流场达不到稳态可能会这样。建议用pimpleFoam试试。
    而且,k-omega SST稳定性很好,对初值没有k-omega模型这么敏感吧。



  • @Foamer24[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    @jinjolee k和omega可以用fluent里面初始化的值做参考

    我的模型放进Ansys后貌似有点问题,无法顺利生成enclosure,要把其中一个部件去除才行,而且每次生成boundary later都会报错,无法生成。因此ansys fluent对我来说可能不是最佳方案
    ~~~~~~~~~~
    @tidedrinker[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    simpleFoam是稳态求解器,你这个问题是不是用瞬态比较好?如果你模拟的流场达不到稳态可能会这样。建议用pimpleFoam试试。
    而且,k-omega SST稳定性很好,对初值没有k-omega模型这么敏感吧。

    我试过用pimpleFoam,但貌似设置有点问题,发散问题更严重。所以想先用simpleFoam试试看。
    ~~~~~~~~~~
    @东岳[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    @jinjolee[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    听说k omega对初始值很敏感

    是的,kOmega模型是这样的。不过Omega可以从epsilon计算出来:$\Omega=\frac{\varepsilon}{0.09 k} $

    你这个我感觉完全是边界条件设置的问题,因为不涉及到算法,你会会搞定的 :duang:

    可是我怎么算怎么试,结果还是发散。:zoule:
    有方法可以分辨出k和omega的初始值是太大或太小吗?



  • 我指的是fluent选择k omega初始化的时候,就可以看到fluent里面的初始值。可以拿来和你现在的初始值对照下。



  • @Foamer24[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    我指的是fluent选择k omega初始化的时候,就可以看到fluent里面的初始值。可以拿来和你现在的初始值对照下。
    能请教一下怎么查看吗? 我没怎用ANSYS,不太清楚。
    是不是要先汇入模型,生成网格,然后进了setup才可以看到omega的初始值?



  • @Foamer24[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    我指的是fluent选择k omega初始化的时候,就可以看到fluent里面的初始值。可以拿来和你现在的初始值对照下。

    我上面的留言格式弄错了, 我再发一次。

    能请教一下怎么查看吗? 我没怎用ANSYS,不太清楚。
    是不是要先汇入模型,生成网格,然后进了setup才可以看到omega的初始值?


  • OpenFOAM副教授

    @jinjolee 你拿个最简单的模型导入fluent就可以了,只是借用一下fluent的初始化k和omega值。其实k和omega的值与速度的关系是有理论公式的,fluent也是用这个理论公式算出来的,看一下fluent的帮助手册。或者看这个网址湍流边界条件



  • @dzw05[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    @jinjolee 你拿个最简单的模型导入fluent就可以了,只是借用一下fluent的初始化k和omega值。其实k和omega的值与速度的关系是有理论公式的,fluent也是用这个理论公式算出来的,看一下fluent的帮助手册。或者看这个网址湍流边界条件

    我就是用你给的那个网址里的论坛的工具去算的。只是我摸拟的是外部流场,turbulence length scale比较难预测,只能用eddy viscosity ratio去算,而那论坛里好像也没有相关的资料(例如在什么情况下,eddy viscosity ratio大约是多少之类的)。
    至于ansys fluent,可能是我还没弄清楚设置,试了一下还是发散了。看来得再研究一下。:wanan:
    话说在一个普通的室外环境里(零风速),eddy viscosity ratio大约是多少?有参考值吗?


  • OpenFOAM副教授

    @jinjolee eddy viscosity ratio一般设大于1小于10,给个10应该就行。另外,你这个不一定是湍流模型导致的发散,直接用层流先试试,逐个排除原因。


  • CORE 网格教授 OpenFOAM教授 管理员

    感谢大家畅所欲言,看起来正在往正确的方向走,祝收敛!
    :xiexie:



  • @dzw05[菜鳥求救] K-omega SST 头几秒就崩溃了 中说:

    @jinjolee eddy viscosity ratio一般设大于1小于10,给个10应该就行。另外,你这个不一定是湍流模型导致的发散,直接用层流先试试,逐个排除原因。

    可是奇怪的是我之前用k-epsilon模型来算,当时用了:

    turbulentKE=0.17557;
    turbulentEpsilon=0.05677;
    

    这样来算,反而结果有收敛,另外reynold number也得出为1638000,所以我认为基本上得用湍流模型。不过这个CdA比正常的低(我的结果只有0.26左右,而从其他的实验结果得来是0.46左右)。我当时用的是:

    flowVelocity=12.5;
    

    ,而k是用turbulence intensity level = 2.73%。至于这个2.73%和epsilon的值是怎得来的,老实说我忘了...:chigua:

    我之后从不同文献中看到湍流的k和epsilon初始值应该要以模拟的环境在现实中的特征来作基础(例如我在模拟单车在行驶,虽然在CFD中我给了个15m/s的风速,但turbulence intensity level的初始应该不是以这个15m /s的速度来算,而是以真实环境(0风速)来算,所以turbulence intensity level应该很低,大约0.2~0.02%左右)。然而我换了turbulence intensity level=0.2%后,结果就一直发散。

    快烦死了:mihu: :mihu:


  • OpenFOAM副教授

    @jinjolee openFOam里的湍流模型有些时候确实不好用。(1)用你现在的条件计算出k、omega和nut的值,保证初始场initialValue也等于计算出的值;(2)检查壁面函数是否合适;(3)检查空间离散格式fvSchemes。


 

Forest
Mountains