颗粒碰撞中的move-collide subCycles与哪些参数有关？

星星星星晴

@zhangxc0223 好像没有具体的硬球模型，我们都是自己写的code。
不知道我猜的对不对，如果说你要是将软球的相关参数设置成硬球是不是就可以了？

我们用的stochastic collision model，所以不需要追踪所有的parcel。相较而言就算的快很多。

zhangxc0223

@星星星星晴我理解的硬球是球的变形非常小，杨氏模量非常大，但是这种情况下，反而move-collide subcycles的数量变得非常大，计算非常慢。是不是我理解错硬球模型了？

星星星星晴

不好意思我没用过软球，你具体查一下文献吧~

星星星星晴

@zhangxc0223 我最夸张的建议就是，你吧所有计算这些东西的都不需要，你只要保留跟踪的部分，然后自己写一个两者在什么条件下碰撞，碰撞以后速度是怎么变化的，手动计算一下质量动能动量守恒，就是硬求模型

zhangxc0223

@星星星星晴哈哈这个不夸张，确实就是这个理，谢谢你的建议。

古奇古奇啦啦

@星星星星晴在颗粒碰撞中的move-collide subCycles与哪些参数有关？中说：

stochastic collision model

您好，想请问一下，constant/kinematicCloudProperties下面设置stochasticCollisionModel的参数有没有可以参考的网址或者资料呢？我看网上的资料颗粒间碰撞设置的基本上是

collisionModel pairCollision; 
pairCollisionCoeffs
{
}
stochasticCollisionModel none;

那是否可以理解为stochasticCollisionModel是关于碰撞后颗粒自身反弹和溅起颗粒的角度和速度的概率分布的设置呢？

星星星星晴

@古奇古奇啦啦碰撞模型分为deterministic collision model 和stochastic collision model 这两个完全是两个不同的方向。
建议你找一下论文看看这两者的区别～前者是根据trackingdata 后者是根据碰撞概率

具体参数设置要看一下code

古奇古奇啦啦

@星星星星晴谢谢您的回复～那我看看这两个方向的文章。

我找了一些修改 collisionModel pairCollision的设置（下面的代码），虽然体现出了颗粒间碰撞的结果，但是不符合想实现的效果，因为我想实现颗粒与床面上颗粒碰撞的概率模型，所以觉得还是用stochastic collision model根据碰撞概率设置颗粒碰撞好一些～

我现在困惑的除了deterministic collision model 和stochastic collision model的理论区别，还有，如果不用deterministic collision model ，那么关于stochastic collision model在OpenFOAM中的实现可以修改的参数都有什么？
下面的代码是我参考GOFUN2018关于OPENFOAM颗粒模拟的讲义写的，类似这样的，不知道stochastic collision model有没有相关的tutorial可供学习一下。

    collisionModel pairCollision;

    pairCollisionCoeffs
    {
        maxInteractionDistance  320e-6;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.12;
            b                   1.5;
            mu                  0.52;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };
        
        wallModel wallSpringSliderDashpot;

        wallSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            youngsModulus   1e10;
            poissonsRatio   0.23;
            alpha           0.12;
            b               1.5;
            mu              0.43;
            cohesionEnergyDensity 0;
        };
        U U.air;
    }
    stochasticCollisionModel none;

    radiation off;
}

星星星星晴

@古奇古奇啦啦 stochastic collision model 是根据cell中的颗粒浓度和相对速度计算每两个parcel的碰撞概率，这部分你最好查一下openfoam的原始code看一下，如何判断碰撞，碰撞后的结果，速度，组分，温度等变化，直接查一下code就好， OF中原生的kinematicParcel是没有stochastic collision 模型的，你要在sprayParcel中查。

https://github.com/OpenFOAM/OpenFOAM-8/blob/master/src/lagrangian/spray/submodels/StochasticCollision/ORourkeCollision/ORourkeCollision.C

星星星星晴

@古奇古奇啦啦

至少of8中是没有的，不知道为什么of之前的版本一直都将sprayparcel单独拿出来，很费解为什么不放到一起，好像of9已经改过来了，但是我没用过。。

不符合理想效果可能是你pairspring model的东西设置的不对，杨氏模量什么的，这部分我不是很了解，因为不同的parcel可能有的性质都不一样，GOFUN那个是他们那个case的设置，你具体用什么parcel，parcel的物理性质什么的是怎么设置的，有没有热相关的东西？这些都有可能改变结果，具体问题具体分析吧

CFD中文网

颗粒碰撞中的move-collide subCycles与哪些参数有关？