距壁面第二个网格中心速度代码求助



  • 壁面函数代码里面距壁面第一个网格中心速度是Up=turbModel.U().boundaryField()[patchi].patchInternalField()
    我还需要第二个网格中心的速度,请问大家代码应该怎么写?
    :zoule: :zoule: :zoule:



  • 下面这种情况红点是第一个网格中心,哪个绿点是第二个网格呢?

    捕获.PNG



  • 抱歉问题描述不严谨
    是结构六面体网格的中心点:wocao:
    无标题.png



  • 即使六面体,边边角角的第二层网格也不好判断是哪个,
    目前OpenFOAM没有现成的可以调用其他层网格的代码,虽然寻找某个网格周围的网格是可以实现的,
    捕获.PNG


  • 离散相副教授

    这个问题主要是找到边界面正对的面,然后通过neighbor-owner找到单元。可以通过边界面的点找到共点的edge,然后找到另一顶点,再来查找面。呃,好像有点繁琐。第二个可以通过particle类,将边界面法相取反作为位移方向,edge长度作为位移,trackToFace找到面。记得用labelPair保存相邻单元。


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2