大家好，我用 icoUncoupledKinematicParcelFoam求解器模拟静电除尘器中颗粒荷电运动轨迹，颗粒考虑曳力、重力和电场力，计算电场力需要电荷量q（通过方程计算，是动态变化的值）和场强E(矢量场，给定）两个参数，我在submodels/Kinematic/particleforces路径下添加了一个电场力子模型，怎么能做到q的实时读取和电场力的实时更新呢，求大神指点。
0.png

@wangfei9088 邮箱联系您了！期待您的回复

第一列为能量，第二列电离系数为关于能量的函数，能量由能量方程求解，第二列需要根据场中的能量数值插值计算当地的电离系数。
e6d35e3b-505f-4342-ad41-cc3aa7905308-image.png
22ee645c-9f2d-4540-84ef-84a11e942e6f-image.png 替代文字

@李东岳 东岳老师您好，我使用通过结合 waves2foam 和 overInterDyMFoam 的用overWaveDyMFoam求解器计算波浪与旋转模型相互作用的时候,重叠部分区域的流体不随时间改变，重叠的部分已经设置成overset边界条件了，想请教老师是否知道这个是什么原因导致的。
113.png
背景部分如图
111.png
重叠部分如图
112.png
boundary文件设置如下：
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v1912 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
oversetPatch
{
type overset;
inGroups 1(overset);
nFaces 56;
startFace 222062;
}
inlet
{
type patch;
nFaces 98;
startFace 222118;
}
bottom
{
type wall;
inGroups 1(wall);
nFaces 1132;
startFace 222216;
}
outlet
{
type patch;
nFaces 98;
startFace 223348;
}
atmosphere
{
type patch;
nFaces 1132;
startFace 223446;
}
frontBack
{
type empty;
inGroups 1(empty);
nFaces 221480;
startFace 224578;
}
sides
{
type overset;
inGroups 1(overset);
nFaces 2000;
startFace 446058;
}
cylinder
{
type wall;
inGroups 1(wall);
nFaces 80;
startFace 448058;
}
)

// ************************************************************************* //

@李东岳 果然跑起来了！

是的，大约20秒之后

@zhe 谢谢，我按照方法排查尝试了一下，发现并不是粒子的加入影响了液体速度场，二而是MPPICInterFoam对于setFields对于输入速度场就有问题。下图一个是MPPICInterFoam计算的，一个是interFoam计算的，同样设置的5m/s液体下落速度，但结果截然不同。@李东岳 也请麻烦李老师能否帮忙看一下，谢谢！
2a82db74-6e61-43e1-a687-5a6ca08d9055-image.png
7751b300-7c97-4692-97a7-48f19f694107-image.png
96ee3830-8f39-4ffe-8986-3988a9d011df-image.png

@RQLI_CFD 一些涡量很杂乱的区域，导致Cd系数会在某些时刻突然蹦很高。

@coolhhh 应该是平均值

我是openFoam10，如果只是删了processor没有改动constant的话，在更改controlDict起始时间后可直接输入decomposePar然后进行求解器求解。
decomposePar命令会根据controlDict中的起始时间拷贝对应的U、p作为初始数据以及当前constant文件夹里的polymesh文件夹读取网格拓扑信息。

请问题主解决此问题了吗？

一眼望去，全是进步空间

哈哈，幽默，有意思

@Samuel-Tu 请问这个功能在paraview的哪里可以找到呢

@李东岳 f084b820-6ac4-4a71-aa27-4cc598955778-image.png 李老师，这是我文献中看到的，这种是咋实现的，没什么思路现在

@李东岳 好的，东岳老师，我试试按照您说得这样算一下。

另外我应该补充一下，100μm直径气泡的实验测得的浮升终速度大约是0.005~0.008m/s.

e584ada7-1a05-4bb6-8cd8-f1e1997c43cc-image.png

参考文献 Pawliszak, P., Ulaganathan, V., Bradshaw-Hajek, B.H., Manica, R., Beattie, D.A., Krasowska, M., 2019. Mobile or Immobile? Rise Velocity of Air Bubbles in High-Purity Water. J. Phys. Chem. C 123, 15131–15138.

@fubianhanshu 你好可以加个联系方式吗，想请教一些问题

@Prometheus10

期待期待

@李东岳 在 质量入口边界 中说：

你考虑的是多相流还是单相流

李老师您好，感谢您的回复我考虑的是多相流，但是入口处只有气体，也只有气体进入

我也不是很懂，从代码角度来看，貌似只有进入了if判断语句中才会进行变量的输出？是不是得把IO操作放在外层

@李东岳 好的，谢谢老师