CFD中文网

壁面的速度确实不能为0，壁面的surface是0. 从这个图来看，【猜测】是粒子直径太大了导致的。这个速度场不太正确。

@zhangxc0223
Hi
因为并没有做deterministic的模拟，所以code也基本没看过。。
stochastic collision 在计算collision probablity的时候是有限制的，具体你可以看一下O‘roucke的论文。具体是怎么定义的，怎么得到一个non-dimensional collision probablity的，还有如何决定这两个parcel发生了碰撞。

虽然碰撞在这个instantaneously 随机的，但是在统计学上，比如你有1M次碰撞的话，总的碰撞的结果就是相似的。
这也就是为什么stochastic碰撞会比deterministic相对来说’便宜'的原因，通过统计学的方法达到相似的碰撞结果，但是不需要追踪所有的parcel。

lagnrangian maxCo 存在的意义，个人认为，在一定程度上可以理解为再一次离散了每个euler time step，所以maxCo一定程度上决定了每个lagrangian time step的大小，和Euler 的Co 定义是相似的。

碰撞是不会被错过的，因为deterministic理论上来说是追踪所有的parcel，所以下一步这个parcel会在哪，走的路径应该也是计算的，所以不会存在错过的，stochastic collision 会计算两个parcel的碰撞概率，当概率很高的时候，也是不会错过的，而且就算是这次错过了，在茫茫parcel海中，总会碰到的。。统计学上是不会错过的

以上仅为个人理解

最近对两个概念非常混淆，即摩擦系数和切线回弹系数，这两个系数是否就是同一个概念？

比如在以下的颗粒pairCollision参数设置中，法向碰撞的回弹系数根据alpha设定。但是切向的回弹系数，是不是就是这个mu?

pairSpringSliderDashpotCoeffs { useEquivalentSize no; alpha 0.01; b 1.5; mu 0.1; cohesionEnergyDensity 0; collisionResolutionSteps 1; };

@dingcy 关掉耦合应该是oneway，也可以这么用，要是DPMFoam满足你使用要求的话就可以。

抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径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，施加在流体相的曳力作用。不知道这个理解对不对？

加了引号好像也不起作用呢，都是显示读取了这个场文件，但是就是没有文件输出出来。

Reading set description: cloudPlane001 Time = 0.082 Reading fields: volVectorField: U.air volScalarField::Internal: kinematicCloud:UCoeff volVectorField::Internal: kinematicCloud:UTrans Executing functionObjects End

@马乔 想借这里请问您一下，如果我只想输出颗粒的Uc的话（写入到每个时间步的/lagrangian/kinematicCloud中），应该怎么操作？

我代码只有入门水平，实在看不太明白您的一系列操作修改，望您有空能指教:zoule:

@zhangxc0223 我不清楚你怎么做的，不过我刚发现可以在paraview中可以显示粒子ID，随着时间步推移，好像每个粒子的ID号是不变的，这样就可以通过后处理输出某个固定ID的各个时间步的信息了，可以交流下

感谢您的回复

@古奇古奇啦啦

至少of8中是没有的，不知道为什么of之前的版本一直都将sprayparcel单独拿出来，很费解为什么不放到一起，好像of9已经改过来了，但是我没用过。。

不符合理想效果可能是你pairspring model的东西设置的不对，杨氏模量什么的，这部分我不是很了解，因为不同的parcel可能有的性质都不一样，GOFUN那个是他们那个case的设置，你具体用什么parcel，parcel的物理性质什么的是怎么设置的，有没有热相关的东西？这些都有可能改变结果，具体问题具体分析吧

@zhangxc0223 我正在用矩方法算这个耶.. 针对你的问题，感觉跟这个很像，https://www.jianguoyun.com/p/DSS2YxIQ9s3ZBhjX_6oC 看图3，二者速度无差异

现在明白了，谢谢东岳老师讲解

@浪迹天大 ControlDict里面的writeControl和writeInterval似乎会覆盖前面的？

application twoPhaseEulerFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 2; deltaT 1e-4; writeControl timeStep; writeInterval 100; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable on; adjustTimeStep no; maxCo 0.9; maxDeltaT 1e-04; functions { volFieldValue1 { type volFieldValue; libs ("libfieldFunctionObjects.so"); log true; //writeControl writeTime; writeControl timeStep; writeInterval 1; //"这里的writeInterval 1 似乎会覆盖前面的100导致所有的输出都是每个时间步都输出？" writeFields true; regionType all; operation volAverage; //weightField alpha1; fields ( alpha.particles U.air U.particles p p_rgh ); } }

@zhangxc0223 您好，我们没有特别关注过长度对循环边界LES的影响，Kara2012的Turbulence statistics in compound channels with deep and shallow overbank flows是循环边界做的，在stoesser的书里面有更加详细的图。

@hurricane007 脑残了，之前写错地方了。现在找到了

@东岳 固体密度是2650kg/m3 液体是1000kg/m3,没有轨迹交叉现象，颗粒的体积分数小于1%

这样的话，管流的速度不知道为何会持续增涨，然后在入口处变成一个非常大的值后计算崩溃。我之后改成了入口为totalPressure，才能实现稳定的速度，计算也能保持稳定。请问这两个边界条件分别在什么情况下使用呀？

应该在头文件中包含

#include "fvOptions.H"

不过拉格朗日求解器不能处理cyclic边界，因此你这么做 :135:

@zhangxc0223 是的

数值上，湍流的出现主要取决于精度。在FVM这种精度不高的情况下须添加随机扰动来制造湍流同时进行维持。我不确定你这种情况湍流会不会维持。你可以试试