CFD中文网

各位大神，我最近在写一个类的时候遇到一个问题，想请假一下。首先我在这个spanwiseAverage类的patchAddressing方法中需要用到某一patch的点坐标和面心坐标，我分别通过patch的localPoints()方法和facecentres()方法获取了两个pointField的引用。随后我需要通过点坐标来计算一个变量geoMin，问题就出在了此处。如果我在这个句代码前面通过Info输出一下patchCentre的大小或者是值是可以得到正确的结果的，后面的程序也可以正常的运行。但是如果不在计算geoMin变量之前对patchCentre进行操作，在计算完geoMin后再输出patchCentre的大小就会变为0。我在计算geoMin的时候并没有对patchCentre这个变量进行任何改变。想请问一下，有没有大神遇到过类似的问题或者这个bug产生的原因。

const Foam::labelList& Foam::spanwiseAverage::patchAddressing(const polyPatch& patch) const
{
    const edgeList& patchEdge = patch.edges();
    const pointField& patchPoint = patch.localPoints();
    const pointField& patchCentre = patch.faceCentres();

    Info << patchCentre.size() << endl;  // 若在此处输出patchCentre的size是可以获取正确的值的

    scalar geoMin = min(patchPoint & averagePatchsDir_);

    Info << patchCentre.size() << endl;  // 若注释掉上一条Info语句，在这里输出patchCentre的size得到的值变为了0

    ......
}

@东岳 是类似与MultiRegionFoam的方法吗？存在多套网格分别对应不同的region？这种方法我没试过。我用的求解器是interfoam，现在只是想的是从一套网格上生成一个subMesh，这样绝大部分的网格节点是对应的，求解标量C方程所需要的数据可以通过原始的网格map过来。

各位大神：
我需要在OF里面求解一个标量的对流扩散方程，这个方程有一个特点就是它的底边界和速度、压力等量的底边界不在同一个高度，就像下面的图一样，左侧为流体的网格，右侧为我所要求解的标量场C的网格。C的底边界和流体的底边界之间有一个高度差δ。在求解的时候采用的方法是在左侧的fvMesh的基础上，构造一个fvMeshSubset来用于求解标量场C，这样需要我设置一个cellSet来选中用于求解标量场C的网格区域，但是我设置cellSet的时候发现，对于最底层网格1来说，要么就是会被全部选中，要么就是无法选中，并不能从高度δ的位置来构建这个subMesh。因此想请教一些大神有什么好解决办法？谢谢！！