海集公司fluent技术支持，1.21日之前在线答疑

lhp1137

楼主好人啊。我也顺便问个问题。我计算的是一个二维自维持振荡问题（好多文献都这样说），我采用层流算法也得到了类似的结果，k-e-rng也可以。而别的模型都不行了，一般都是最后得到一个稳定的解（和文献上说的不同）。因为雷诺数比较小，不能确定什么时候转变成紊流，所以想用一个能够计算过渡流动的模型。不知道用k-e-rng模型是不是就可以说是准确，因为没有具体的试验数据，是不是可以根据它的计算流场和试验流场相似就确定计算的合理性和准确性呢？多谢多谢

hi-key

正好碰到问题了，我模拟区域内气流及污染物扩散运动特征，入口条件给出气流速度，点或线污染源怎么给？应该是和domain中的air不同的气体。二维
如果是模拟室内洁净方面的污染物扩散，你可以用DPM设定点或者线的污染源（用injection)，设定污染物的体积力（如果是细菌可以用水蒸气代替），然后求解没有污染物时的流场，收敛后打开DPM，再求解。我想这样应该行得通，因为前两天刚做了个类似的。是室内的污染物。

hi-key

我也问一下,我用udf时用compile 可以通过并且可以计算,但是用解释时却会提示我第几行有问题,这是什么原因啊.还有就是我算的是非稳态动网格,每个时间步网格会变化,我想问的是残差的问题,在一个时间步里残差达到收敛的标准,可是在下个时间步开始的残差会很大,然后又收敛,所以整个残差图是个锯齿型的,在一个时间步的最后收敛,在下一个时间步的第一次迭代又很大,但随着迭代又降下去了,请问这样正常吗?呵呵,还有个问题,就是我有时在计算时会出现某几个单元的温度小于最小温度的极限1k,请问这主要是网格问题还是算法的选择问题.最后一个问题是我可以用分离求解器算可压的超音速流吗?效果怎么样?谢谢了啊

一共四个问题：
1、解释型udf要求函数中不能包含子函数，有可能使这个原因。
2、动网格求解中的锯齿形残差是正常的。
3、计算过程中，某几个单元温度越界没有关系，只要它能回来，别发散就行，这个原因较多，一般来说不是网格的问题，比如在求解超音速流动时，使用segregate算法时就容易出现。
4、分离求解器可以用来求解可压超音速流，但要求结构网格同时网格间距跳跃不大，马赫数小于2.5时，将松弛因子全部设为0.1，是可以计算的。马赫数高了就完蛋了。还有计算结果跟couple的一样。
这些是我的经验，也不一定十分准确，仅供参考

zhubin

楼主真是大好人!!!我现在算的是弹丸出炮口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每个时间步很难收敛,我就换成分离的,而且我亚松弛因子调的很大,压力有0.9,动量为0.8,其他都是1.我用的是三角形的非结构网格,气流速度有1500多,马赫数达到了4左右,可是在每个时间步上它收敛很快,一般16步左右就收敛了.想问下你,象我这样设那么大的亚松弛因子可以吗?还有就想问你,象我这样的非稳态怎么看每一个时间步是否收敛了?我现在就是看残差图.

wormxie

正好碰到问题了，我模拟区域内气流及污染物扩散运动特征，入口条件给出气流速度，点或线污染源怎么给？应该是和domain中的air不同的气体。二维
如果是模拟室内洁净方面的污染物扩散，你可以用DPM设定点或者线的污染源（用injection)，设定污染物的体积力（如果是细菌可以用水蒸气代替），然后求解没有污染物时的流场，收敛后打开DPM，再求解。我想这样应该行得通，因为前两天刚做了个类似的。是室内的污染物。
先谢谢楼主
我模拟的开放空间，大气中的污染物扩散，可以这样吗？
我以为增加一个inlet就好，结果发现没法加污染物

hi-key

我计算的是一个二维自维持振荡问题（好多文献都这样说），我采用层流算法也得到了类似的结果，k-e-rng也可以。而别的模型都不行了，一般都是最后得到一个稳定的解（和文献上说的不同）。因为雷诺数比较小，不能确定什么时候转变成紊流，所以想用一个能够计算过渡流动的模型。不知道用k-e-rng模型是不是就可以说是准确，因为没有具体的试验数据，是不是可以根据它的计算流场和试验流场相似就确定计算的合理性和准确性呢？多谢多谢
这种问题的要求比较高，类似的问题我只算过卡门涡阶的。不过当时是用层流算得。你这个例子里面如果跟湍流模型敏感，我建议你可以尝试以下方法：
选用其他湍流模型，然后在进出口边界处的湍流相设置时，使用湍流强度和粘性比然后将这两个数值全部给零，再计算时使用绝对压力计算。也许会有变化，也许没用，你可以试试，把结果告诉我。谢谢，另外在所有的湍流模型中k-e-rng是最适合计算低雷诺数湍流模型的，当然你也可以尝试真正的低雷诺数湍流模型，需要在用户界面中输入命令行激活，至于怎么激活我忘了，哈哈，不好意思，等我查到了给你哈。
另外判断结果是否正确只能靠实验或者查文献了，流态特征相似只能证明大体上没错，但是精度就不知道了。

hi-key

我现在算的是弹丸出炮口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每个时间步很难收敛,我就换成分离的,而且我亚松弛因子调的很大,压力有0.9,动量为0.8,其他都是1.我用的是三角形的非结构网格,气流速度有1500多,马赫数达到了4左右,可是在每个时间步上它收敛很快,一般16步左右就收敛了.想问下你,象我这样设那么大的亚松弛因子可以吗?还有就想问你,象我这样的非稳态怎么看每一个时间步是否收敛了?我现在就是看残差图
直觉不太对，不过如果你的时间步长很小也许就可以避过每个时间步长内发散的可能性，如果时间步长比较大就不对了哈，另外炮弹出口如果直线运动用layering，比remesh好多了，又准确又快。
PS，膛内外压力差多大？Mpa级？

hi-key

我模拟的开放空间，大气中的污染物扩散，可以这样吗？
我以为增加一个inlet就好，结果发现没法加污染物

开放空间中的污染物排放用DPM也可以，但是因为粒子运动轨迹长，需要很长的跟踪时间，因此用这个模型计算不经济，你可以用欧拉颗粒流来计算，可以在进口位置指定污染物的排放量。但是点源，线源都必须以边界条件的形式存在，当然DPM模型也要求如此。

wormxie

我模拟的开放空间，大气中的污染物扩散，可以这样吗？
我以为增加一个inlet就好，结果发现没法加污染物

开放空间中的污染物排放用DPM也可以，但是因为粒子运动轨迹长，需要很长的跟踪时间，因此用这个模型计算不经济，你可以用欧拉颗粒流来计算，可以在进口位置指定污染物的排放量。但是点源，线源都必须以边界条件的形式存在，当然DPM模型也要求如此。
我一直在刷：）
我是以边界条件的形式给出的，问题是如何定义污染物是与domain中空气不同的气体，inlet中好像没有关于物质成分的定义，污染物是气体，物理特性不重要，只要结果可以看出浓度分布就好

hi-key

问题是如何定义污染物是与domain中空气不同的气体，inlet中好像没有关于物质成分的定义，污染物是气体，物理特性不重要，只要结果可以看出浓度分布就好
拿欧拉模型讲：
打开多相流的欧拉模型，
define-materials  定义污染物的物性，如gas
define-phase，指定primary phase 为空气，secondary phase 为gas，以及颗粒直径
define-boundary con，选中inlet，在下面的mixture那里向下拉，选中gas，设定其进口速度和组分分量（一般为1）。初始化时记得全场gas的浓度为0。就好了。
try try

zhubin

是啊,有60多mpa.我时间步长为0.000001s

zhubin

还要请教楼主,想我这样的非稳态计算怎么样看每一步是否收敛了?看残差吗?

wormxie

着急，我模拟的是机动车排放(假设为点或线源)在城市内的扩散，二维，k-e turbulence model
希望楼主更明白，谢谢

linfanmr

我计算车室内流场，出风口很多，回风口（排风口）只有一个，采用hex core网格、
RNG k-e 模型，计算过程监视出口（回风口）流量，计算到300步的时候差不多成了一条直线，而残差图也差不多成了直线，接着迭代300步，残差图还是下不去，也就是不能降低到低于1e-3的水平，在report里面看net imbalance 在0.1%左右，请问这样能不能算收敛了？
压力我用的standard ，其他的是二阶迎风，有什么办法能使残差降到1e-3以下，谢谢?

hi-key

还要请教楼主,想我这样的非稳态计算怎么样看每一步是否收敛了?看残差
可以采用稳态算法中的办法，监控某点的某个变量值，要在流场变化区域内，如喷口位置。如果平了就好了。