lagrangian类的一个悬而未决的小Bug

同心圆柱网格的二维画法
上回说到利用二维圆柱网格计算气固反应流，遇到一个小问题，颗粒竟然碰到二维网格的wedge边界导致报错

// $FOAM_SRC/lagrangian/basic/particle/particleTemplates.C

template
void Foam::particle::hitWedgePatch(TrackCloudType& cloud, trackingData& td)
{
    FatalErrorInFunction
        << "Hitting a wedge patch should not be possible."
        << abort(FatalError);

    hitSymmetryPatch(cloud, td);
}

这个报错实在想不通，因为在一个二维算例（xy平面）中，颗粒怎么会碰到z方向上的边界？

想不通

而且，根据前期测试，这个Error可能只出现在某些机器上......

解决这个问题最简单粗暴的方法就是在OF源代码中注释掉上述的FatalErrorInFunction，这样当颗粒碰到wedge边界时会反弹（相当于碰到Symmetry边界）当然这样会污染OF源代码，并且换机器跑代码的话都需要对OF源代码重新编译，很不优美。

第二种方法就是将$FOAM_SRC/lagrangian/basic中的代码移植到自定义的求解器中，将上述Error注释掉然后重新编译。这个当然没有什么难度，求解器编译成功，算例也都正常运行，就在Case运行完毕End之后，问题来了...

*** Error in `/home/**/OpenFOAM/***/platforms/linux64IccDPInt32Opt/bin/***': double free or corruption (!prev): 0x0000000000d89ce0 ***

在CDFOnline^[1]上有人报告过类似的问题，所以这个Error的原因应该还是依赖没有完全解决，导致在最后算例完成时重复释放了内存。

因为这个Error只出现在加载了我自定义的$FOAM_SRC/lagrangian/basic库的时候，所以问题应该还是在这个库及其依赖中。basic所依赖的库只有meshTools一个，如下

// $FOAM_SRC/lagrangian/basic/Make/options
EXE_INC = \
    -I$(LIB_SRC)/meshTools/lnInclude

LIB_LIBS = \
    -lmeshTools

但是看了一眼meshTools，里面茫茫多的代码...

// $FOAM_SRC/meshTools
├── AABBTree
├── algorithms
├── AMIInterpolation
├── cellClassification
├── cellDist
├── cellFeatures
├── cellQuality
├── coordinate
├── edgeFaceCirculator
├── edgeMesh
├── fields
├── indexedOctree
├── lnInclude
├── Make
├── mappedPatches
├── meshSearch
├── meshStructure
├── meshTools
├── momentOfInertia
├── output
├── PatchFunction1
├── polyMeshZipUpCells
├── primitiveMeshGeometry
├── processorLOD
├── regionModel
├── regionSplit
├── regionSplit2D
├── searchableSurfaces
├── sets
├── surfaceSets
├── tetOverlapVolume
├── triSurface
└── twoDPointCorrector

我选择妥协，采用方法一，虽然不优美，但是它方便啊

当然，如果有Foamer有更好的办法，欢迎下方留言讨论~

---

1. https://www.cfd-online.com/Forums/openfoam-programming-development/96886-ptrlist-lagrangian-clouds-double-free-corruption-crash.html ↩