如何将一些固体particles的体积整合成一个球形的有效parcel

[zhe]

@星星星星晴 您好,首先谢谢回复.我看了spray里的stochastic collision model, 对于里面的code不知道理解对不对,希望您指正.像图里看到的,对于P1和P2的定义,是分别落在不同的cell[i] 和cell [j]里吗?那么这两个粒子都是虚拟的吗?

[星星星星晴]

@zhe 不是，这个是一个cell中选取两个真parcel，我个人觉得你的想法和这个模型比较相像，有些地方你可以借用一下~

[zhe]

@Jacob-Z 不好意思,我没用过CFDEM,不好对它有什么说法.但是对于它们的比较的话,你可以直接在论坛上搜索"DEM和DPM"这个题目.很多有见地的总结在里面.希望能帮到你.

[zhe]

@星星星星晴 好的,谢谢啦.我再试试,有进展了会分享一下,谢谢您!

[Jacob Z]

@zhe 好的 谢谢！

[zhe]

@星星星星晴 您好,我这两天都在研究您提到的stochastic collision model, 帮助很大,首先谢谢您.我有一个小问题,如果在MPPIC的方法中,既然应用了inter-particle stress (packing model)再加上damping model, 那么是不是stochastic collision model就不应该并存了?或者说是不是二选一比较妥当.虽然在现在OpenFOAM自带的MPPIC的tutorials,在kinematicCloudProperties文件里,两个都显示了,但是stochastic collision model都是默认了none. 那么如果我想要在MPPIC中加入另外的碰撞的模型,还是需要关掉packing,damping,isotropic这3种 stresses? 希望您有时间了可以回复.

[星星星星晴]

@zhe hi 因为我也没用过你的模型，我们的collision 模型是in house的，我觉得你可以组合一下他们，我个人以为 stochastic collision model 更偏向于去寻找碰撞对，而你用的packing,damping,isotropic模型更倾向于碰撞后的处理，我没用过OF原来的模型，所以了解不深。。
或者你可以combine一下他们

[星星星星晴]

@zhe 毕竟 Orouke 的后处理其实很简单，更偏向于硬球模型。

[zhe]

@星星星星晴 谢谢您的及时回复.那我明白了,我也是比较倾向于碰撞处理.真心感谢这几天来您的回复,帮助极大.

[星星星星晴]

@zhe 客气啦 你在美洲？

[zhe]

@星星星星晴 没有,在英国.您也是在国外吧,看您回复的时间,都跟国内有差距.

[星星星星晴]

@zhe 咱不用客气，邻居啊，你们刚被踢出去。
话说顺便问你个问题，你做average 的时候，是根据Nparcel做average 还是Nparticle做平均。
比如你求一个cell中 粒子的平均质量的时候。

[zhe]

@星星星星晴 我现在用的都是以parcel为计算单位的.在code里也是p.nparticle.每一个parcel里包含的particel的数量是自己的一项输入,那么就不能再以每一个particle再做计算.并且,parcel的概念本身也是为了减少计算量,如果每次平均再回到particle,那就达不到原来的目的了.

[zhe]

@星星星星晴 所以把所有的parcel(p.nparticle)总和的话,就变成你说的Nparcel了.

[星星星星晴]

@zhe 但是这么想，如果一个cell中有n个不同size的parcel，

size(um)	nParcel	nParticle	Average Parcel size	Average Particle size
100	50	300	(100*50+175*15+500*35)/(50+15+35)	(100*50*300+175*15*200+500*35*100)/(50*300+15*200+35*100)
175	15	200	251.25	175.5813953
500	35	100

[星星星星晴]

@zhe 最终parcel的概念到底怎么搞的 我有的时候还是很迷糊。。。涉及到后处理，处理的我一头包

[zhe]

@星星星星晴 这里面是两个环节的.一开始是粒子的注入.需要在注入后计算出它在cell里的state.这个state会store下来以计算后的形式,我认为这里是parcel了.

// Define local properties at beginning of time step
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
const scalar np0 = nParticle_;
const scalar mass0 = mass();

接着是追踪parcel运动过程,在这里就只是调用之前calc出来的parcels的质量了.不同zise的parcel是一直被延续的,中间是以整体形式在计算域里被计算和追踪的.

// averaging sums
forAllConstIter(typename TrackCloudType, cloud, iter)
{
    const typename TrackCloudType::parcelType& p = iter();
    const tetIndices tetIs = p.currentTetIndices();

    const scalar m = p.nParticle()*p.mass();

这是我的理解,不知道是不是都对,你先看看,是不是符合你的情况.

[星星星星晴]

@zhe
这就存在问题了啊，在开始的时候 parcel的mass 就是 他的size *volume 不考虑particle.
但是比如计算这个parcel受力的时候也用的是parcel的mass

OpenFOAM-8/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H
inline Foam::scalar Foam::KinematicParcel<ParcelType>::mass() const
{
    return rho_*volume();
}

但是你后半部分的mass是在cloud中

OpenFOAM-8/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H

就是代表着，在cloud中，追踪的mass是parcel这个合集的mass，真的很头大

[zhe]

@星星星星晴 在KinematicCloudI.H中计算了system mass也是以parcel来计算的,这是cloud层面上的.在我上个回复里的后半段是是追踪粒子的文档里的,不是cloud怎么的. 在追踪过程里,还是以parcel为追踪项的.

template<class ParcelType>
template<class TrackCloudType>
inline void Foam::MPPICParcel<ParcelType>::trackingData::updateAverages
(
    const TrackCloudType& cloud
)
{
    // zero the sums
    volumeAverage_() = 0;
    radiusAverage_() = 0;
    rhoAverage_() = 0;
    uAverage_() = Zero;
    uSqrAverage_() = 0;
    frequencyAverage_() = 0;
    massAverage_() = 0;

    // temporary weights
    autoPtr<AveragingMethod<scalar>> weightAveragePtr
    (
        AveragingMethod<scalar>::New
        (
            IOobject
            (
                cloud.name() + ":weightAverage",
                cloud.db().time().timeName(),
                cloud.mesh()
            ),
            cloud.solution().dict(),
            cloud.mesh()
        )
    );
    AveragingMethod<scalar>& weightAverage = weightAveragePtr();

    // averaging sums
    forAllConstIter(typename TrackCloudType, cloud, iter)
    {
        const typename TrackCloudType::parcelType& p = iter();
        const tetIndices tetIs = p.currentTetIndices();

        const scalar m = p.nParticle()*p.mass();

        volumeAverage_->add(p.coordinates(), tetIs, p.nParticle()*p.volume());
        rhoAverage_->add(p.coordinates(), tetIs, m*p.rho());
        uAverage_->add(p.coordinates(), tetIs, m*p.U());
        massAverage_->add(p.coordinates(), tetIs, m);
    }

[星星星星晴]

@zhe 啊 我突然明白了 谢谢，这个定义总在拐弯，老了，脑子不行了