CFD中文网

自带的ensightWriter就能整体输出结果

@oitocfd 如果只是交换边界两侧的值，试试syncTools::swapBoundaryCellList. 如果只是查找离边界两层的网格可以用下FaceCellWave这个类(这个类会迭代全场网格，你需要修改类的迭代行为）。如果是固定距离就trisurfaceSearch了。话说为嘛不直接移指下foam-Extend的ib呢？

@micro 这个会在每个核里存一个共享的偏移数组，根据当前核编号和核locasize累计。在decompose后，全局编号是当前核的局部编号+偏移，这个一般会跟未分解前不一致（比如blockmesh编号是按定义的x-y-z来的）。一般在需要reconstruct， distribute,聚合的地方用到吧，我也就用过一次。

@upc_ngh 试试在这两句lookup之前显式调用下massStick()呢？

你可能需要仿照MRFZone类写一个平动的类了，但是相对来说要更加容易些，主要工作是修正面通量，推荐用三角分解计算。还要注意下include的边界面上速度修正。

或者你应该在intU[1]=...之前加个if(cellID > -1).

@杨英狄 主程序里调用的构造参数不对呗，按匹配的传入就行

@oitocfd 如果粒子追踪的话可以尝试下lagrangian库，当前网格周围一圈网格可以直接用primitiveMesh.cellCells()

@oitocfd 我错了，meshSearch的查找效率很高，是primitiveMesh不行。你这是在写pic吗？

@oitocfd 可能是精度问题？手动写入的有截断误差？刚好这个点位于face上呢？还有不管是primitiveMesh和meshSearch的findcell效率都挺低的，都要去遍历所有的cell，不知道写什么需求要这么写？

@zz给我趴下 可以试下columeAveraged

@同学博 你这么操作map是对的，我开想的是查找key

@李东岳 这么写呢？

            List<DynamicList<label>> pL(U_.size());
            forAllConstIter(typename MomentumCloud<CloudType>, *this, iter)
            {                                                    
                     const parcelType& p = iter();                                              
                    pL[p.cell()].append(p.origId());                                           
            }            

或者

std::Multimap<label,label> Lp;
forAllConstIter(typename MomentumCloud<CloudType>, *this, iter)
{
     const parcelType& p = iter();  
     Lp.insert(std::pair<label,label>(p.cell(), p.origId()));
}

都是我云的

@同学博 map用迭代器访问效率并不高吧，而且最好不要直接访问value的吧，如果用value匹配，这样复杂度又上来了

用个List<dynamicList<label>> parcelsToCell, 然后遍历粒子，得到单元编号，然后将粒子塞到对应的dynamicList中。但是有空的。
或者用个std::Multimap<label>，cell编号就是key，value是粒子id。

1，不是，这个边界会构建一个虚拟的2D规则网格，然后利用读入的R和L计算出波动速度，然后加上U并map到真实边界网格上。
2，这些值是模拟计算的值，sample或者其他方式得到的。

@oitocfd 直接用MPPICInterFoam就可以了

像这样子嘛？

没得问题，只要在这里不使用具体类对象，而使用指针或引用，就不用预先知道类定义。然后模板类里一般都是用静态函数创建模板参数对象（静态函数属于类，不属于某一对象），有点类似全局函数。

@浪迹天大 不，我不相信，你都写了这么久了