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



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



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



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

    1. 想确认一下您指的turbulent dispersion force是不是指kinematicCloudProperties文件中的dispersionModel这一项,而不是像sphereDrag这样的某个particleForce?
    2. 在您帖子里的对拉格朗日模型的Stochastic tracking model,是不就是对应OpenFOAM中的kinematicCloudProperties→dispersionModel →stochasticDispersionRAS?
    3. 我试了一下用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。
    你这个问题是湍流模型的原因。感觉设置有问题



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

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



  • 报错是什么?



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2