CFD中文网

@Jasper-0 建议更换org的openfoam11+版本即可，改动没有想象的那么大

@coolhhh 发现问题了，因为我的是两相流，需要设定field 为U.water

@Jasper-0 你现在能堆积，是因为把youngsModulus系数增大从而增大了碰撞力。这会引起无法堆积的，你可以看你的结果，应该会一直蹦，因为排斥力过大。这个系数因材料而定，不能随便给。比如沙子大概e+06好像。碰撞力的大小需要在计算时调节需要的参数，不能改变粒子的物理属性。

@Jasper-0 您好，请问这个问题您现在解决了吗？

@jasper-0 更新一下目前进度，我尝试在KinematicLookupTableInjection.C里将injectors_设置为MUST_READ_IF_MODIFIED,

template<class CloudType> Foam::KinematicLookupTableInjection<CloudType>::KinematicLookupTableInjection ( const dictionary& dict, CloudType& owner, const word& modelName ) : InjectionModel<CloudType>(dict, owner, modelName, typeName), inputFileName_(this->coeffDict().lookup("inputFile")), duration_(this->coeffDict().getScalar("duration")), parcelsPerSecond_(this->coeffDict().getScalar("parcelsPerSecond")), randomise_(this->coeffDict().getBool("randomise")), injectors_ ( IOobject ( inputFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE ) ),

但是在运行solver的时候提示不支持自动重新读取，有没有大佬知道这种情况咋解决呢？

GlobalIOList<kinematicParcelInjectionData> kinematicLookupTableInjection constructed with IOobject::MUST_READ_IF_MODIFIED but GlobalIOList<kinematicParcelInjectionData> does not support automatic rereading.

@jasper-0 我觉得不需要

@李东岳 感谢李老师回复！

@星星星星晴 感谢大佬指点！

@李东岳 李老师，用的是reactingTwoPhaseEulerFoam里的class method

切面上拉格朗日粒子的速度？

点完全在面上？通过时，那就抓拉格朗日粒子的坐标，找到符合坐标的粒子输出速度。。。。如果刚好一个时间步穿过，坐标范围怎么选？大了包含很多在面上蹭蹭不穿过的，小了，快的粒子刷就跑了，那穿过时的速度用穿前穿后的均值？

可以在那个表面位置设置一层内部面，捕捉粒子，撞上去速度归零，归零前记录粒子速度:quwan:

@李东岳 楼上这俩一样的啊老师，我就截图截的大小不一样了。但是一楼那个算例其实是我想模拟的，但网格太大，所以我先用楼上这个网格小的进行试验的

搜 openQBMM

@jasper-0 在 拉格朗日粒子浓度与网格大小的关系 中说：

各位大佬好，
本人最近在想是不是拉格朗日粒子的浓度是与网格大小相关的？因为在openfoam中浓度定义为：

theta[celli] += p.nParticle()*p.volume(); theta.primitiveFieldRef() /= mesh_.V();

当网格变小的时候，浓度不就自然而然变大了吗？如果是这样的话，该如何跟实验进行验证呢？

有关，但是个人理解是没法与实验验证，

volume fraction的作用的个人理解：

确定拉格朗日是用什么coupling 确定当前的模拟是否physically成立，因为of中lagrangian 都是point-mass。所以有的时候DPM会出现问题，这也算是EL法的固有缺陷吧，我记得上线好像是65%，具体忘了在哪看到的（有可能类似摩尔定律的故事）。 如果你的parcel很多，（不然为什么用EL不是VOF），达到一个统计学上的匹配就行，比如说velocity profile, velocity RMS/Std, local size distribution，等，你要相信，茫茫parcel海中，只要统计学上会发生，总有一个parcel会发生某些表现，只是概率问题而已。

@李东岳 大佬 气泡可以用拉格朗日，就是流场别太复杂，而且还需要引入bubble induced turbulence还有不同的drag模型，dispersion模型，还有 eccentricity, wall effect,etc～

把 return 挨个屏蔽，看看屏蔽了哪个报错。。。:yes:
就知道调用哪个了

我想知道它是什么意思呢？

好像是用 u 和 dim 初始化了一个新的 fvVectorMatrix 类型的变量

大家好，我想请教一下KinematicParcelTrackingDataI.H中的TrackCloudType& cloud里的cloud指的是什么？在如下代码中：

rhoInterp_ ( interpolation<scalar>::New ( cloud.solution().interpolationSchemes(), cloud.rho() ) ),

cloud直接调用了函数rho()，为了弄清楚cloud是哪个类，我改变了rho()为任意的字母，比如a()，但是重新编译依然可以编译成功，这是为什么呢?

@bestucan 好的，非常感谢大佬的指点！

@李东岳 感谢大佬！我好好研读一下 :142:

@bestucan 嗯嗯，我发现alpha是在twoPhaseSystem.C中计算的

:chitang:

@sxz0823 你好，请问有关concentration和parcelconcentration您是怎么计算的，求指导。