Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
@星星星星晴 在 不知道为什么会生成特别蛋疼的field 中说:
AUX_npIter
会不会这个值在边界处太太太小了,导致结果边界处太大了
能否看一下AUX_npIter和ROOTSMALL的量级比较
@东岳 npiter 是数 parcel 有多少的,所以再某些地方会出现0,因此加上ROOTVSMALL,但是我改了SMALL,结果一样,只是值小了,异常数据的位置还是相同的
@yfclark
这个是npiter的 cell值,ROOTSVMALL 一般都在 e-300左右?open foam自带的
@星星星星晴
IVT_UpMean = AUX_SumUp/(AUX_npIter + ROOTVSMALL);
这个公式是做什么?如果AUX_npIter(粒子数)是0,那IVT_UpMean必然是一个非常大的值。我觉得你要不然障眼法一下:
IVT_UpMean
IVT_UpMean = pos0(AUX_npIter)*AUX_SumUp/(AUX_npIter + ROOTVSMALL);
这个公式是计算平均速度,npiter这几步中,将parcel的速度加和,然后求平均值 IVT_UpMean = AUX_SumUp/(AUX_npIter + ROOTVSMALL);
但是情况是这样的,我的cell值都没啥事,但是point的值就很诡异了
@东岳 请麻烦再问一下 pos0是干嘛用的。。我刚才搜了一下 没看到
找到以前的case看了一下, 图1 是SumUp的cell值,靠近壁面比较小
图2 是nptier的cell值,可以看到靠近墙壁的值有一点,25左右
最后两个field 相除的的结果是这样
应该是除法的问题,壁面附近的粒子数为0,分母和分子都很小,但是分子的量级更小,才会出现局部大值,建议特别处理这种情况
我的cell值都没啥事,但是point的值就很诡异了
那可能是边界条件的问题,