CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    2
    8
    1449
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • H
      hotairballon1997 最后由 编辑

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

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

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

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        H 2 条回复 最后回复 回复 引用
        • H
          hotairballon1997 @李东岳 最后由 编辑

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

          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)
          
          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 编辑

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

            线上CFD课程开始报名:http://www.dyfluid.com/class.html

            CFD高性能服务器 http://dyfluid.com/servers.html

            H 1 条回复 最后回复 回复 引用
            • H
              hotairballon1997 @李东岳 最后由 编辑

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

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

              1 条回复 最后回复 回复 引用
              • 李东岳
                李东岳 管理员 最后由 编辑

                报错是什么?

                线上CFD课程开始报名:http://www.dyfluid.com/class.html

                CFD高性能服务器 http://dyfluid.com/servers.html

                1 条回复 最后回复 回复 引用
                • H
                  hotairballon1997 最后由 编辑

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

                  1 条回复 最后回复 回复 引用
                  • H
                    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
                             );
                    
                    1 条回复 最后回复 回复 引用
                    • First post
                      Last post