cannot find file /kinematiccloud/collisionrecordspairaccessed

星星星星晴

@zhe 哈哈说实话我也是第一次知道拉格朗日下的也能mapfield，不过上次尝试过使用mapfield，map出来的效果很差，所以后来就不用map了，那个时候还没有parcel呢。
我们这边有把lagrangian field 做一个时均，然后写到欧拉场里的。

pairRecord是什么？了解一下？deterministic collision model？还是软球模型的什么啊？

zhe

@星星星星晴也不是什么必须的文档。只是在读取朗格朗日数据时，因为我用的就是DPM的碰撞模型，所以需要读取这个文件。我看了一下collidingParcelO.C, 在readField里应该可以去掉吧，得试试，看看影不影响DPM。
对了，方不方便问一下，你那边怎么把拉格朗日做的时均啊？理论概念是什么？

星星星星晴

@zhe hi 昨天看到忘了回了，就是求一些统计上的平均速度，速度的RMS什么的。这样比较容易求，而且能直接在paraview上看。

zhe

@星星星星晴再请问一下哈。因为在用openfoam追踪粒子运动的过程中，发现漏洞很大。尤其是openfoam主要在意粒子的运动模式，比如力的作用都在速度变化量中体现，那么对于位置的更新有没有什么别的可实施的？比如保持粒子不动，是不是只能不断地减小作用合力在粒子身上产生的变化量？有时间了在回复，不急的，谢谢🙏

星星星星晴

@zhe 不好意思，我没有get到你的点
拉格朗日不就是建立在 $F=ma$ 么
那如果不通过速度的变化，怎么知道加速度，怎么算作用力呢。

保持粒子不动的话，完全可以把粒子做成wall，或者setzone呗
看你关注的啥啊，关注流场还是关注粒子

zhe

@星星星星晴没错，通过力的作用增加速度变化量，我的意思是：在高浓度粒子流中，openfoam中本身的力都作用在单个粒子上，而相互间只有碰撞力，那不能维持粒子由动到静，或由静到动的过程，因为它可是视为永远在动。把粒子作成wall是永远不能动的意思对吗？那就失去了粒子的感觉。setzone是啥呀？

星星星星晴

@zhe 那你得看你浓度有多高了，不然怎么会有 1way 2way 4way的区别呢。

不过我还是没有get到什么叫维持parcel由动到静，或者由静到动的过程是什么
parcel受力就会发生速度的变化，速度变化展示了收到了多大的力的作用，然后对流场发生了多大的影响。主要就是计算那个$UTrans$嘛。
有人也有在solver中植入某种力场的教程什么的。
还有人提出过parcel induced turbulence的模型，比如bubble induced turbulence什么的。
你指的是不仅仅只有momentum exchange？比如reaction啊 parcel induced turbulence什么的？

我以为你是想做类似圆柱绕流之颗粒绕流～不好意思理解错啦。

zhe

@星星星星晴是我没表示清楚。我的粒子浓度最大0.62。绕柱是我的一个case，但是感觉做出来不是那么理想。粒子在柱周围无法很好的堆积起来。那么从力的角度来说，有一种摩擦效果的力，存在在粒子间，但没有被表达出来。如果有足够的摩擦，是可以使得粒子静止堆积。当然也不会妨碍更大的水流造成的迁移（不是像wall一样永远不能动）。您所说的水流对粒子的影响，无论是动量交换还是turbulent，或是dispersion都是在流-固之间的，但咱高浓度粒子流中，固-固之间的很重要。

星星星星晴

@zhe 哈哈我没到你这么高的浓度，所以不是很了解，既然你的浓度这么大，你用的drag model 考虑alpha（theta）什么的了么？另外有可能是turbulence 没有解决好？有可能是dispersionmodel的问题，不好说，前两天我还看到一篇文章讲这个东西来着。。我回头找找。

pp之间的作用力没想过，正常能想到的pp之间的作用力，磁力？静电力？不然还能有啥除了碰撞意外以外？
相对论么？

暂时不知道想不到能有啥

星星星星晴

@zhe

LES and RANS calculations of particle dispersion behind a wall-mounted cubic obstacle
Marco Atzori a,b,∗, Sergio Chibbaro c, Christophe Duwig d, Luca Brandt
2022

闲着看一下这篇文章，他的解释是可能是涉及到湍流没有解好啊，Stokes 什么的。

而且正常情况下 spray中只考虑了drag， lift, 如果固体parcel 可能会有shear lift 什么的，如果bubble 会有什么basset force， historyforce什么的，可能你真是缺了某个模型，找找看
万一呢，这不就一片文章出来了

zhe

@星星星星晴感谢回复啊。说的挺有见解的。因为目前所有的4-way基本上都是考虑弹簧碰撞，最多添加转动消散能量，但其实有些情况是不适合的。比如只有接触，并无overlap，那么一种类似摩擦的效果应该如何表达。这就对于高粒子流很重要，因为高粒子流中不会轻易受到流体影响，况且说的再细一些，流体在高粒子流中本就不该是连续的。而是像渗透在多孔介质。这只是我的想法，跟您交流交流，多听听别人的意见，我觉得会更能填补空缺。真心再次感谢您回复，以后有机会您多指教啊。

星星星星晴

@zhe 客气啥这玩意现在没一个通用的模型解释嘛，要不咱还咋毕业。。。
哈哈哈我这都是吹牛逼可以，实际操作就费了别客气

祝收敛～
祝稳定～

你是做流化床么？流化床应该和你说的概念差不多吧

星星星星晴

@zhe 我们这边考虑碰撞的时候回头一个impact efficiency 的考量，就是小parcel可能会因为大parcel周围的流场从而造成绕着大parcel 而不是与大parcel发生碰撞这种情况

sommerfeld and lain
influence of droplet collision….
2020

zhe

@星星星星晴不是流化床，我是做海岸的，哈哈。粒子启动和堆积很重要。所以研究粒子之间的比较多。

zhe

@星星星星晴那你们还是比较多的是考虑周围流体对粒子的影响。我已经不能再考虑那些了，否则又能做一个课题啦，哈哈哈

zhe

@星星星星晴有个问题请教一下，有时间帮忙看看。在计算粒子drag时，如果是ErgunWenYuDragForce。分了两个公式：
if (alphac < 0.8)
{
return forceSuSp
(
Zero,
(mass/p.rho())
(150.0(1.0 - alphac)/alphac + 1.75*Re)muc/(alphacsqr(p.d()))

    );
}
else
{
    return forceSuSp
    (
        Zero,
        (mass/p.rho())
       *0.75*CdRe(alphac*Re)*muc*pow(alphac, -2.65)/(alphac*sqr(p.d()))
    );
}

按道理将，这两个计算的都是beta值，也就是相对速度前面的那一部分。可是在第二个公式里，好像连两次的相对速度都已经包含在内了。还有，在KinematicParcel.C中，计算加速度时：const vector acp = (Fcp.Sp()*td.Uc() + Fcp.Su())/massEff;
这里乘以td.Uc()，感觉不对啊。应该在这里乘以相对速度才更合理一些阿。麻烦给看看，实在是推导出来跟公式对不上。谢谢。

zhe

@星星星星晴理解错了。第二个公式CdRe(Re)本就只有一个CdRe。但是连个公式都在分母处除以的alphac,这造成跟最后的公式不同。还有第二个问题，应该是乘以相对速度(Uf-Up)的，怎么乘以里td.Uc().请解惑。

星星星星晴

@zhe
在这个定义中 OpenFOAM-8/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/forceSuSp/forceSuSp.H

 Helper container for force Su and Sp terms.
        F = Sp(U - Up) + Su
    Explicit contribution, Su specified as a force
    Implicit coefficient, Sp specified as force/velocity

所以最后在return的时候 Sp和Su是完全不同的俩东西。 U就是td.Uc()

星星星星晴

@zhe
我用的of8 是这样的你用的那个版本啊？

template<class CloudType>
Foam::forceSuSp Foam::ErgunWenYuDragForce<CloudType>::calcCoupled
(
    const typename CloudType::parcelType& p,
    const typename CloudType::parcelType::trackingData& td,
    const scalar dt,
    const scalar mass,
    const scalar Re,
    const scalar muc
) const
{
    const scalar alphac =
        this->alphacInterp().interpolate
        (
            p.coordinates(),
            p.currentTetIndices()
        );

    if (alphac < 0.8)
    {
        return forceSuSp
        (
            Zero,
            mass/p.rho()*(150*(1 - alphac)/alphac + 1.75*Re)
           *muc/(alphac*sqr(p.d()))
        );
    }
    else
    {
        return WenYuDragForce<CloudType>::calcCoupled(p, td, dt, mass, Re, muc);
    }
}

星星星星晴

@zhe 我查了一下of5 好像和你看的wrgunwenyudragforce 这部分相同，我就简单的做了一下量纲分析，最后得到的都是 $\frac{kg}{s}$ 所以得到的是Sp，Sp乘以速度就是Su 就是力了$\frac{kgm}{s^2}$.

因为我也不做这种考虑$\alpha$的drag，仅仅是根据自己之前的理解和你强行解释了一波，这玩意属于of的基础定义了～

CFD中文网

cannot find file /kinematiccloud/collisionrecordspairaccessed