CFD中文网

@李东岳 我这么写的，你这么一说，我觉得是不是要按照啥格式给点参数的值啊

我使用的是icoUncoupledKinematicParcelFoam求解器，我想多加个颗粒上的作用力，比如压力梯度力（PressureGradient），我在算例里添加了关键词，然后我还在creatFields里声明了p(我觉得可能用到这个变量所以我声明了)，0文件里也加了p文件，但是最终报错说这个力需要一个特定的字典，请问这是问啥呢？该咋样添加呢？

@李东岳 嗯，你说的没错，这个问题解决了，我是用“cat a.out | grep "Current number of parcels" > new.log”这个提取了计算输出的想要的部分，然后再编个小程序把中间不需要的时间步给删掉，只留下整数秒，然后就解决了。

@Samuel-Tu 我在用拉格朗日方法计算颗粒的时候也出现了相似的问题，单核能算，但是并行算几步就不行了，想问一下你当初怎么解决了，你楼上的帖子说“往scalarField里面加数据的时候超过了定义的大小”我没看懂。

如何输出出口escape的颗粒数量随时间变化的曲线？之前算单向流的时候我输出出口浓度随时间变化曲线是通过paraview那个plot selection over time, 但是现在这个颗粒好像不能那么整了....

@李东岳 好的好的，谢谢大家

拉格朗日求解器中注射模型里的massTotal，这个总质量指的是入射颗粒的总质量吗？如果是的话，我在设置参数的时候已经设置了颗粒密度，入射颗粒数量，颗粒半径，这个massTotal必须满足密度乘体积的数量关系吗？

@李东岳 有两个问题，第一，计算稳态流场的情况下buoyantBoussinesqSimpleFoam可以采用这种方式吗？第二，这种方法是不是没计算离散相方程。

我在连铸中间包（一种冶金反应器）入口射入颗粒，反应器内流场为给定流场，分别采用两种方式射入即脉冲加入和连续加入，颗粒在壁面我设置了rebound边界条件（其中e为1，mu为0）。结果如下图所示，为什么颗粒并没有在反应器内传输呢？而是在底面聚集，我在之前计算流场的时候速度给的是无滑移边界条件，也就是边界处速度为0。是这个原因吗？有啥解决方法吗？

图1 连续加入

图2 脉冲加入

图3 脉冲加入底面颗粒分布细节

@dxl 给定流场，然后"kinematicCloudProperties"中关掉耦合那种就叫one way呗，用的是DPMFoam吗？还是你用的其他方法或者求解器？

@dxl 你好，我想问一下，你这个"solidparticle做one-way"是啥意思？我用icoUncoupleKinematicParcelFoam做固定流场里入口射入tracer颗粒出了点问题，想了解一下你这个是咋做的？

自问自答，我粗略的整出来了，如图1，要注意画红圈地方的选择

图1

OpenFOAM-6.0自带的DPMFoam里面的那个Goldschmidt算完之后怎么查看颗粒图啊，就像下面的图一样（如图1）。我咋服务器算完之后，用foamToVTK转格式之后，把VTK文件下载下来之后（如图2），用paraview咋处理也出不来呢（如图3）？

图1

图2

图3

@李东岳 最近论坛的搜索功能是不是不好使了

我最终的计算结果已经得到的情况下，通过什么方式可以将边界处的值赋成相邻内部网格的结果值？
另外，一个小问题，OpenFOAM里绝对值符号是啥呀？mag(标量)，可以表示一个标量的绝对值吗？

@李东岳 我觉得不能算是很标准的固定梯度边界条件，因为那个f，就是扩散通量的数值是随着边界处标量浓度的变化而改变的，有一个表达式。

如下图所示，我想写一个边界条件，当标量浓度在壁面处以通量F吸附，现在壁面通量速度f我通过codedFixedValue边界条件已经写出来了，但是浓度C的边界条件不知道怎么设置？我开始是设置为固定值0，但是计算后边界附近的浓度都是负值，这是不对的。我想向大家问一下想实现我最开始说的边界条件，应该怎么处理？

@李东岳 这个式子里的delta是大涡里过滤的那个delta吧，与网格尺寸有关，能不能提示一下输出它的方法呀

@李东岳 那这样的所谓强制稳态结果有意义吗？我的意思是物理意义，比如说一个反应器我想得到流动稳定条件(这里的稳定是相对稳定，现实情况下可能不存在稳定)下的流场，这时候我用大涡算，那就要从0时刻开始算，算很多很多秒吗？所谓的强制稳态的流场可以表达这种流场吗？

我把buoyantBoussinesqSimpleFoam求解器里湍流模型改成LES Smagorinsky，一样可以算，那这样算出的结果是稳态结果，可以叫LES的稳态结果？