CFD中文网

你这是什么版本啊，为嘛那一通运算没有量纲检查。forAll这个只会遍历内场。

1.我猜是symmetry这个边界条件，试一下改成wall，然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查，外部软件的网格节点编号和blockMesh的逻辑不一样。
2.可能是前后的empty边界问题，你可以在starccm里生成3D网格后，重新在openfoam里extrude来生成前后empty边界
类似这样
constructFrom patch;
sourceCase "../w3-d_hc1-3";
sourcePatches (symFront);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName symBack;

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
extrudeModel linearNormal;

nLayers 1;

expansionRatio 1.0;

linearNormalCoeffs
{
thickness 1;
}

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces false; //true;

// Merge small edges. Fraction of bounding box.
mergeTol 0;

@XieXiaoyang 好文章，我去研究研究

最近多次使用foamDictinary感觉真是太用了，自己在文件里面加的keys也能识别出来并修改

我用fluent中dpm模型做壁面磨损，想请问一下颗粒的固相体积分数是用混合模型设定的还是需要自己计算，若需要自己计算应该怎么计算，还有里面的颗粒质量流率怎么计算，可以解答一下或者提供相关文献吗？

@benqing 会不会是因为固体域没有时间步，所以网格不发生变化？:xiezuoye:

@李东岳 好的谢谢您，我再试试

@浪迹天大 师兄可以加个微信吗？wallong95:laughing:

非常感谢您

@李东岳 请问在哪里下载？东岳流体？

今天在编译liggghts和cfdem的时候也出现如下报错：
185e1dd1-371e-43a4-9c79-1be84f630d3e-image.png

结果证明需要如下依赖包：

sudo apt-get install libvtk6-dev

另外好像OF8没法完成和liggghts的耦合，我用的是5.x的。

同理，下述网页中都需要处理$\sum_f p_f^{t}\bfS_f$为$\frac{1}{V_\rP}\sum_f p_f^{t}\bfS_f$，均已更正，谢谢

http://dyfluid.com/rhoPimpleFoam.html
http://dyfluid.com/rhoSimpleFoam.html
http://dyfluid.com/simplefoam.html

@李东岳 感谢李老师的回复，一下子就通透了

@Yan 好的，谢谢你，我已经学会从命令面板里设置怎么输出数据了。
再次感谢！

抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径D=0.01m的空气管流中，速度Uf=15m/s，雷诺数Re=10000。颗粒的斯托克斯数St=10，直径是大概dp=50×10^(-6)m，密度为ρp=1000kg/m^3，颗粒-空气的质量比是Φ=0.1，双向耦合，只考虑颗粒曳力。

在这种情况下，我查看OpenFOAM计算得到的kinamaticCloude:UTrans文件中，动量转移都是在10^(-14)量级（单位是kg×m/s），请问这合理吗，为什么只有这么小。。？

据我所知，这个UTrans文件表征的就是颗粒→流体相的动量转移，也就是颗粒相通过源项Sp，施加在流体相的曳力作用。不知道这个理解对不对？

那你先解释一下，为什么你认为这是个不可能事件

C1和C2差异很大，但Q3和w3几乎相同，请问如何解释

它违背了什么物理规律？和已有的什么结论相悖？

@walten 在 如何CFD自学成才 中说：

@闻久STU
乐观点， 六年:high: 这样才学得深入

不能提:135: 长期抗战是肯定的了

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。

密度就是单位体积的质量，速度的物质导数就是加速度。
所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？

力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。

物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。

要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。

另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。