有没有utility或者force，能使颗粒在RANS模拟中像LES那样发散开来

hotairballon1997

最近在做颗粒射流的RANS模拟，流体相可以实现与实验差不多符合。但是对于颗粒相，因为流体场都是雷诺时均的结果，没有湍流的涡结构来吹散颗粒。结果显然不准确，所以想问问，在OP中有没有相应的功能，能使得颗粒在RANS模拟中能像LES，DNS一样离散地发散开来？

李东岳

我们之前详细的研究过这个东西目前DPM这面的turbulent dispersion force结果都不太好你可以尝试一下

hotairballon1997

@东岳谢谢解答，刚刚看完您的贴子“拉格朗日中的湍流分散力模型”，也做了一些测试。
请问李老师，

想确认一下您指的turbulent dispersion force是不是指kinematicCloudProperties文件中的dispersionModel这一项，而不是像sphereDrag这样的某个particleForce？
在您帖子里的对拉格朗日模型的Stochastic tracking model，是不就是对应OpenFOAM中的kinematicCloudProperties→dispersionModel →stochasticDispersionRAS？
我试了一下用icoUncoupledKinematicParcelFoam，这个stochasticDispersionRAS可以正常运行，但是用DPMFoam会报以下错误，不知道应该从何修改？

...
Selecting RAS turbulence model kEpsilon
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No finite volume options present

Starting time loop

Courant Number mean: 0.000115464 max: 0.150963
deltaT = 0.000125
Time = 0.000125

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud


--> FOAM FATAL ERROR:
Turbulence model not found in mesh database
Database objects include:
36
(
U.air
alpha.air
alphaPhic
alphacf
boundary
cellZones
data
epsilon
faceZones
faces
fvOptions
fvSchemes
fvSolution
k
kinematicCloud
kinematicCloud:UCoeff
kinematicCloud:UTrans
kinematicCloudOutputProperties
kinematicCloudProperties
mu.air
neighbour
nu
nut
owner
p
phi.air
pointConstraints
pointMesh
pointZones
points
rho.air
tetBasePtIs
transportProperties
turbulenceProperties
volPointInterpolate(U.air)
volPointInterpolation
)


    From function Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::DispersionRASModel<CloudType>::kModel() const [with CloudType = Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > >]
    in file lnInclude/DispersionRASModel.C at line 51.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:?
#3  Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:?
#4  ? at ??:?
#5  ? at ??:?
#6  ? at ??:?
#7  ? at ??:?
#8  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9  ? at ??:?
Aborted (core dumped)

李东岳

是的，就是dispersionModel。
你这个问题是湍流模型的原因。感觉设置有问题

hotairballon1997

@东岳我查了下，好像说是DPMFoam中，流体相的湍流模型和分散模型定义在lagrangian/turbulence/submodels/Kinematic中，而颗粒相的湍流分散模型在$FOAM_SRC/lagrangian/turbulence/submodels/... 中，这样颗粒的dispersionModel找不到turbulence model是不是导致我这里报错的原因？

我对修改code的问题实在还没太上手，想问一下遇到这样的问题应该怎么解决？

李东岳

报错是什么？

hotairballon1997

抱歉没看到回复提醒，就是报错三楼贴出的问题，不过后来发现是我之前把DPMFoam里面的湍流模型调用修改了，导致出了问题。OpenFOAM里本来的DPMFoam是没有问题的~

hotairballon1997

@李东岳在有没有utility或者force，能使颗粒在RANS模拟中像LES那样发散开来中说：

我们之前详细的研究过这个东西目前DPM这面的turbulent dispersion force结果都不太好你可以尝试一下

想请教一下李老师，在这个stochasticDispersionRAS.C的代码中，dt代表的是什么，并没有看出它的出处在哪。。
从结果上来看，这个dt好像不是模拟中的timestep，它的值对于每个颗都不一样。

     // Member Functions
  
         //- Update (disperse particles)
         virtual vector update
         (
             const scalar dt,
             const label celli,
             const vector& U,
             const vector& Uc,
             vector& UTurb,
             scalar& tTurb
         );