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

[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找到包含插值点的网格编号，我在想分区后，如果刚好这个插值点处于两个分解区域的交界面，这样岂不是就有两个网格都包含这个插值点了？这样不会有问题么。。

[keke]

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