[菜鳥求救] K-omega SST 头几秒就崩溃了
-
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} $
你这个我感觉完全是边界条件设置的问题,因为不涉及到算法,你会会搞定的
可是我怎么算怎么试,结果还是发散。
有方法可以分辨出k和omega的初始值是太大或太小吗? -
@dzw05 在 [菜鳥求救] K-omega SST 头几秒就崩溃了 中说:
@jinjolee 你拿个最简单的模型导入fluent就可以了,只是借用一下fluent的初始化k和omega值。其实k和omega的值与速度的关系是有理论公式的,fluent也是用这个理论公式算出来的,看一下fluent的帮助手册。或者看这个网址湍流边界条件
我就是用你给的那个网址里的论坛的工具去算的。只是我摸拟的是外部流场,turbulence length scale比较难预测,只能用eddy viscosity ratio去算,而那论坛里好像也没有相关的资料(例如在什么情况下,eddy viscosity ratio大约是多少之类的)。
至于ansys fluent,可能是我还没弄清楚设置,试了一下还是发散了。看来得再研究一下。
话说在一个普通的室外环境里(零风速),eddy viscosity ratio大约是多少?有参考值吗? -
@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的值是怎得来的,老实说我忘了...
我之后从不同文献中看到湍流的k和epsilon初始值应该要以模拟的环境在现实中的特征来作基础(例如我在模拟单车在行驶,虽然在CFD中我给了个15m/s的风速,但turbulence intensity level的初始应该不是以这个15m /s的速度来算,而是以真实环境(0风速)来算,所以turbulence intensity level应该很低,大约0.2~0.02%左右)。然而我换了turbulence intensity level=0.2%后,结果就一直发散。
快烦死了