CFD中文网

非牛顿流体力拖曳力的计算应该不一样，不知道你用的哪个方法计算，如果直接套用牛顿流体的模型会导致拖曳力计算偏小吧

@东岳 最新版本里官网推荐的就OF 5.x, CFDEM 的versionInfo.H文件：
word CFDEMversion="cfdem-3.8.1";
word compatibleLIGGGHTSversion="3.8.0";
word OFversion="5.x-commit-538044ac05c4672b37c7df607dca1116fa88df88";
可能新版本代码还没更新，不兼容？

@东岳 我用的是CFDEM 来算的，调用的应该是OF 5.x，或者说CFDEM里面计算颗粒位置的代码用的是老版本的？

为什么reconstruct的时候会提示报错
--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 5 the punctuation token ')'
file: /processor0/0.01/lagrangian/particleCloud/positions at line 5.
而且切换其他OF版本就可以正常reconstruct
除此之外当颗粒数目较多时，压力会出现异常（非气体时出现明显负压）

如题，不用paraFoam处理，怎么求某个标量在整个计算区域的总量

@李东岳 MULES(FCT)格式是什么意思，用mules 的vanleer 和不用的对比吗

@余正东 你这发的看不见啊

在二相流模拟中，捕捉界面后怎么让界面更平滑来避免计算的不稳定呢，openfoam中经常通过mules来求解相方程，它起到的什么作用？

5.0版本的可以正常安装，但是低版本如2.4.0安装就会报错：

-bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: syntax error near u                              nexpected token `;&'
-bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: `        [ "$FOAM_V                              ERBOSE" -a "$PS1" ] && echo "Sourcing: $1" 1>&2'

想用二相欧拉模型做一个固液二相流问题，下层是固体（体积分数0.6），上层是水，入口处水给一个固定流量，入口的边界条件该怎么给呢？示意图：

@dzw05 比如说我的网格check是这样
Mesh non-orthogonality Max: 28.8083 average: 6.97813
Non-orthogonality check OK.
那么我的fvschemes和PIMPLE里面是不是都得非正交修正，只有一个修正或者修正次数对结果会有多大的影响呢？对非正交修正不是特别理解

那非正交网格需要修正的话是不是fvsolution 里面的

PIMPLE
{
    momentumPredictor no;
    nCorrectors     1;
    nNonOrthogonalCorrectors 0;
    correctAlpha    0;
    nAlphaCorr      1;
    pRefCell        0;
    pRefValue       0;
}

这种非正交修正次数要从0改成1之类的？为什么有人前面用的是高斯修正，后面非正交修正次数却是0呢

主要就是ddt 的Euler 和Euler implicit , laplacianSchems 的Gauss linear corrected 和uncorrected等

@李东岳 把压力梯度视为流体产生加速度的原因，压力梯度力=质量*加速度，我在二相欧拉法里面看过类似的做法，拉格朗日方程不太懂

@xiaofenger 你把$internalField 直接写成uniform (0 0 0)不就好咯

@李东岳 嗯，调分辨率就可以了。东岳老师我最近在用二相欧拉模型求解一个水沙运动问题，但是算到一定时间就会蹦掉，您有空帮我看下我的边界条件设置是不是有问题吗？

@asininno 用line source 的话

@李东岳 默认就是 point source 啊，我用这个连涡都显示不出来

貌似paraview 里面的流线图只能画经过一条直线的流线图，但是如果流场里有分布较为分散的涡旋时就无法全部表现出来，有什么好的办法可以解决吗？

我的计算区域是1m*0.2m,圆柱半径是0.02m,流体是水，圆柱横放在水中。

入口边界给的log速度分布，U=Uf/K*log(y/z0) 摩擦速度Uf=0.01m/s，k,omega 给zeroGradient；

上单是自由液面，底面和圆柱都采用wallFunction，其中k=1.5*(U*I)^2;I=5%;
omega=Cmu^(-1/4)k^(1/2)/L,L是湍流长度，L=0.07D，D=0.04m.

我可以这样设置边界条件吗？