CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    MPPICInterFOAM 粒子大小对导致计算失败

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

      当我把manualinjection所使用的size distribution设置成 1e-3的时候就出现计算失误
      Cloud: kinematicCloud injector: model1
      Added 205 new parcels

      #0  Foam::error::printStack(Foam::Ostream&) at ??:?
      #1  Foam::sigFpe::sigHandler(int) at ??:?
      #2  ? in "/lib64/libc.so.6"
      #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
      #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
      #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
      #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
      #7  Foam::fvMatrix<double>::solve() at ??:?
      #8  Foam::PackingModels::Implicit<Foam::MPPICCloud<Foam::KinematicCloud<Foam::Cloud<Foam::MPPICParcel<Foam::KinematicParcel<Foam::particle> > > > > >::cacheFields(bool) at ??:?
      #9  ? at ??:?
      #10  ? at ??:?
      #11  ? at ??:?
      #12  ? at ??:?
      #13  __libc_start_main in "/lib64/libc.so.6"
      #14  ? at ??:?
      Floating point exception
      

      当我改成比如说8e-4的时候, 计算就没有问题, 能够正常进行. 我觉得是不是因为粒子之间有重叠的关系? 后来我使用scipy中的KDTree将有重叠的粒子删掉, 但是删掉以后发现还是出现错误, 甚至只有七个粒子的时候, 计算还是不能进行. 而且我肉眼观察, 粒子之间间隔已经比较远了, 而且距离边界也有一定的距离. 所以我想问一下, 有遇到这种情况的吗? 会是什么原因造成的呢?

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

        粒子太大了。要小于网格 。

        2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • Q
          qingdong_wang 最后由 编辑

          yes, actually, i also thought that. but is there such strict criterion? The mesh size of my case is 0.000166667, namely 1.6e-4 m. However, when i set the diameter of particle as 5e-4 even 8e-4 m, it also works?

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

            Double check if the lagrangian particle is two-way coupled with the fluid phase. If it is one-way coupled, the restriction can be abandoned.

            2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
            CFD高性能服务器 http://dyfluid.com/servers.html

            1 条回复 最后回复 回复 引用
            • First post
              Last post