同心圆柱网格的二维画法
上回说到利用二维圆柱网格计算气固反应流，遇到一个小问题，颗粒竟然碰到二维网格的wedge边界导致报错
// $FOAM_SRC/lagrangian/basic/particle/particleTemplates.C
template<class TrackCloudType>
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有更好的办法，欢迎下方留言讨论~
