CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    OpenFOAM的插值能够用并行计算吗?

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

      在OpenFOAM的pisoFoam代码改写里面用到了的cellPoint插值方法,在并行计算的时候报如下错误:

      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      Create time
      
      Create mesh for time = 0
      
      
      PISO: Operating solver in PISO mode
      
      Reading field p
      
      Reading field U
      
      Reading field IBForce
      
      Reading/calculating face flux field phi
      
      Selecting incompressible transport model Newtonian
      Selecting turbulence model type RAS
      Selecting RAS turbulence model realizableKE
      RAS
      {
          RASModel        realizableKE;
          turbulence      on;
          printCoeffs     on;
          A0              4;
          C2              1.9;
          sigmak          1;
          sigmaEps        1.2;
      }
      
      No MRF models present
      
      No finite volume options present
      
      Starting time loop
      
      
      Create IB points
      
      Time = 0.005
      
      Courant Number mean: 9.97566e-05 max: 0.0448352
      DILUPBiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 2.34051e-09, No Iterations 1
      DILUPBiCGStab:  Solving for Uy, Initial residual = 1, Final residual = 5.12408e-08, No Iterations 1
      [0] #0  [1] #0  [2] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[3] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
       at ??:?
       at ??:?
       at ??:?
      [2] #1  Foam::sigSegv::sigHandler(int)[1] #1  Foam::sigSegv::sigHandler(int)[0] #1  Foam::sigSegv::sigHandler(int)[3] #1  Foam::sigSegv::sigHandler(int) at ??:?
      [0] #2  ? at ??:?
      [3] #2  ? at ??:?
      [2] #2  ? at ??:?
      [1] #2  ?[DESKTOP-898LF9E:02873] 3 more processes have sent help message help-btl-vader.txt / cma-permission-denied
      [DESKTOP-898LF9E:02873] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
       in /lib/x86_64-linux-gnu/libc.so.6
      [0] #3  Foam::polyMeshTetDecomposition::cellTetIndices(Foam::polyMesh const&, long) in /lib/x86_64-linux-gnu/libc.so.6
      [3] #3  Foam::polyMeshTetDecomposition::cellTetIndices(Foam::polyMesh const&, long) in /lib/x86_64-linux-gnu/libc.so.6
      [2] #3  Foam::polyMeshTetDecomposition::cellTetIndices(Foam::polyMesh const&, long) in /lib/x86_64-linux-gnu/libc.so.6
      [1] #3  Foam::polyMeshTetDecomposition::cellTetIndices(Foam::polyMesh const&, long) at ??:?
      [0] #4  Foam::cellPointWeight::findTetrahedron(Foam::polyMesh const&, Foam::Vector<double> const&, long) at ??:?
      [3] #4  Foam::cellPointWeight::findTetrahedron(Foam::polyMesh const&, Foam::Vector<double> const&, long) at ??:?
      [2] #4  Foam::cellPointWeight::findTetrahedron(Foam::polyMesh const&, Foam::Vector<double> const&, long) at ??:?
      [1] #4  Foam::cellPointWeight::findTetrahedron(Foam::polyMesh const&, Foam::Vector<double> const&, long) at ??:?
       at ??:?
      [0] #5  [3] #5   at ??:?
      [1] #5   at ??:?
      [2] #5  ???? in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      [3] #6  __libc_start_main in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      [1] #6  __libc_start_main in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      [2] #6  __libc_start_main in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      [0] #6  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
      [3] #7   in /lib/x86_64-linux-gnu/libc.so.6
      [1] #7   in /lib/x86_64-linux-gnu/libc.so.6
      [2] #7   in /lib/x86_64-linux-gnu/libc.so.6
      [0] #7  ???? in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      [DESKTOP-898LF9E:02881] *** Process received signal ***
      [DESKTOP-898LF9E:02881] Signal: Segmentation fault (11)
      [DESKTOP-898LF9E:02881] Signal code:  (-6)
      [DESKTOP-898LF9E:02881] Failing at address: 0x3e800000b41
      [DESKTOP-898LF9E:02881] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f438443ef20]
      [DESKTOP-898LF9E:02881] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f438443ee97]
      [DESKTOP-898LF9E:02881] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f438443ef20]
      [DESKTOP-898LF9E:02881] [ 3] /home/samuel/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam24polyMeshTetDecomposition14cellTetIndicesERKNS_8polyMeshEl+0x88)[0x7f43859d5e58]
      [DESKTOP-898LF9E:02881] [ 4] /home/samuel/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam15cellPointWeight15findTetrahedronERKNS_8polyMeshERKNS_6VectorIdEEl+0x65)[0x7f43880a5d75]
      [DESKTOP-898LF9E:02881] [ 5] pisoIBMFoamIPL(+0x38faa)[0x7f438a838faa]
      [DESKTOP-898LF9E:02881] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f4384421b97]
      [DESKTOP-898LF9E:02881] [ 7] pisoIBMFoamIPL(+0x3ebca)[0x7f438a83ebca]
      [DESKTOP-898LF9E:02881] *** End of error message ***
      [DESKTOP-898LF9E:02878] *** Process received signal ***
      [DESKTOP-898LF9E:02878] Signal: Segmentation fault (11)
      [DESKTOP-898LF9E:02878] Signal code:  (-6)
      [DESKTOP-898LF9E:02878] Failing at address: 0x3e800000b3e
      [DESKTOP-898LF9E:02878] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f586703ef20]
      [DESKTOP-898LF9E:02878] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f586703ee97]
      [DESKTOP-898LF9E:02878] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f586703ef20]
      [DESKTOP-898LF9E:02878] [ 3] /home/samuel/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam24polyMeshTetDecomposition14cellTetIndicesERKNS_8polyMeshEl+0x88)[0x7f58685d5e58]
      [DESKTOP-898LF9E:02878] [ 4] /home/samuel/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam15cellPointWeight15findTetrahedronERKNS_8polyMeshERKNS_6VectorIdEEl+0x65)[0x7f586aca5d75]
      [DESKTOP-898LF9E:02878] [ 5] pisoIBMFoamIPL(+0x38faa)[0x7f586d438faa]
      [DESKTOP-898LF9E:02878] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f5867021b97]
      [DESKTOP-898LF9E:02878] [ 7] pisoIBMFoamIPL(+0x3ebca)[0x7f586d43ebca]
      [DESKTOP-898LF9E:02878] *** End of error message ***
       in ~/OpenFOAM/samuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoamIPL
      

      难道说插值不能用并行计算??
      因为有一个命令findCell找到包含插值点的网格编号,我在想分区后,如果刚好这个插值点处于两个分解区域的交界面,这样岂不是就有两个网格都包含这个插值点了?这样不会有问题么。。

      K 1 条回复 最后回复 回复 引用
      • K
        keke @Samuel-Tu 最后由 编辑

        @samuel-tu
        你好这个问题解决了吗?我也是遇到了同样的问题。

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