网格体积

飞火流星jyj

@李东岳李老师可以了，3.0.x没问题，生成的数据是怎么排列的呢？

飞火流星jyj

@队长别开枪请问您测试过吗，我测试了在createField.H加入

Info << "\n-Class: cell----------" << endl;
    cell &cl(cells[100]);
Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl;

报错如下：

createFields.H:33:14: error: ‘cells’ was not declared in this scope
     cell &cl(cells[100]);
              ^
createFields.H:34:52: error: ‘points’ was not declared in this scope
 Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl;
                                                    ^
createFields.H:34:60: error: ‘faces’ was not declared in this scope
 Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl;
                                                            ^

怎么修改就好了

队长别开枪

@飞火流星jyj 我今晚或者明天把整个测试代码发上来吧:happy:

飞火流星jyj

@队长别开枪好的，十分感谢

程迪

@飞火流星jyj
别用cl.centre()

这个东西和mesh.C()似乎不一样，OF有两种中心，几何重心和FVM意义上的中心。二者在三维网格有warp的情况下有一些微小的差别。

飞火流星jyj

@程迪好的，一个是polymesh的一个是扩展的FvMesh，那我该用哪个表示出某网格的x，y，z，这该怎么写啊

程迪

@飞火流星jyj
GeometricField类型有一个component函数，返回的是分量，参数是direction，实际上类型是char/uint_8，你试试'x','y','z'或者1,2,3就好。

飞火流星jyj

@程迪请问怎么查看是0,1,2还是1,2,3这样的编号:confused:

赵一铭

@飞火流星jyj

例如你要输出U的x方向，就是U.x();
不同的版本好像不一样，例如有时候可能是U.component(0);

:big_mouth:

队长别开枪

@飞火流星jyj 不好意思，这两天太忙没顾得上，刚刚把代码发上来。链接地址：meshInfoFoam，它会将体单元的体积、中心和面单元的面积、中心、面积矢量、单位矢量全部写入初始文件夹0，具体结果可以参考链接meshInfoFoam output。代码里还测试了一些涉及单个网格元素操作的成员函数。

zhanghan

@队长别开枪
您好！我试过您的meshInfoFoam了，非常棒的求解器。
我有一个疑问，请问我怎样才能将导出的faceArea与我的cellCenter一一对应上呢？？？

程迪

@zhanghan
不可能对上，两者数目都不一样，一个是面场，一个是体心场。

队长别开枪

@zhanghan 你是指获取属于某一体单元的面单元的信息么？

飞火流星jyj

@队长别开枪请问如果在forAll(,celli)循环里，mesh.c()是指循环到的网格的中心嘛，还是所有的，如果不是怎样写循环到的celli的中心坐标:crying: 求大神指教

队长别开枪

@飞火流星jyj mesh.C()返回的是所有体单元的中心，要获取特定体单元中心的话可以这样写

forAll(mesh.cells(), cellI)
{
    Info << mesh.C()[cellI] << endl;
    // or
    Info << mesh.C().operator[](cellI) << endl;
}

飞火流星jyj

@队长别开枪十分感谢，对于
多个网格的U[celli]的保存，我想用定义一个速度单位的数组存起来，dimensionedvector 好像只是保存一个值吧，这样的数组我不知道该写好，请大神指教

鲸落

@队长别开枪老师您好，请问要是想输出动网格的体积等信息应该作何修改呢？

队长别开枪

@鲸落哪种动网格？overset，AMR？

鲸落

@队长别开枪抱歉回复晚了，动网格文件的代码是这样的

dynamicFvMesh           dynamicMotionSolverFvMesh;           
solver                  velocityLaplacian;                     
velocityLaplacianCoeffs 
{
  diffusivity          uniform;                              
}
v0                      v0 [ 0 1 -1 0 0 0 0 ]   (0 0 0);

桎梏

@队长别开枪

感谢大佬提供的工具。但是用了以后为何总是自动停止呢，0文件夹里什么也没有，

请问是因为我的网格数太多（800w），虚拟机的内存（4g）不够么？

CFD中文网

网格体积