CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    并行边界上网格面心处数值通量(φ_f)的同步

    OpenFOAM
    并行边界数据同步
    1
    1
    157
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 新
      新喻庸 最后由 新喻庸 编辑

      最近在编写程序的时候遇到并行边界上网格面心处数值通量(φ_f)的同步问题,单核运行不存在并行边界,所以没有这个问题。以下图为例:

      Domain.png

      并行计算时,求解区域分为了4个区域,两个区域之间产生了并行边界面,对于内部网格面元,网格面心处的数值通量(φ_f)的计算跟单核运行时是一样的,但是对于并行边界处的网格面元,数值通量(φ_f)在两个相邻并行边界上应该是一样的,因为并行边界的存在不应该改变原来的解。因此,我用了如下的syncFaceList()函数来进行同步,alpha_f是定义在网格面心上的surfaceScalarField,这个函数在头文件"syncTools.H"中有定义,我用了以下三种同步方式,发现流体穿过并行边界时还是产生了误差,单核运行时误差(数值解与解析解之差)可以达到10E-12次方量级,但是并行运行时流体一穿过并行边界,误差就上升到了10E-4量级,我觉得这是因为并行边界上数据同步没有做好的原因,请问有小伙伴遇到过这个问题嘛?使用syncFaceList()函数的正确方式是什么?如果使用syncFaceList()函数进行同步网格面心处的数值通量(φ_f)还是有问题的话,打算自己写一下同步过程。

      // syncTools::syncFaceList(mesh(), alpha_f, maxEqOp<scalar>(), false);
      syncTools::syncFaceList(mesh(), alpha_f, plusEqOp<scalar>(), false);
      // syncTools::syncFaceList(mesh(), alpha_f, eqOp<scalar>(), false);
      
      1 条回复 最后回复 回复 引用
      • First post
        Last post