CFD中文网

请问多重网格求解器可以对不规则几何，非结构网格进行求解吗？

@李子橙 谢谢提供线索。确实是这个问题。我原单位是m，用mm建模后已无此问题

需要生成蒸汽轮机叶片几何，目前有的数据是叶片不同高度的叶型离散点数据。我用的是ansys自带的scdm进行建模。
1） 导入scdm，生成12个横截面曲线。由于曲线不闭合，因此先把曲线fill成面。

2） 放样生成实体模型，发现有的面之间会生成错误，具体如下图所示。

3）原始数据：scdm.zip

问题是有的横截面之间放样不成功，有畸形曲面产生。
想请教大家：
1） 这种只有横截面数据的情况是不是只能通过放样获得三维几何模型？
2） 有什么工具或软件能实现此功能？

用过ssr，不过经常节点被封，还得花时间找节点，不划算

最近想趁双十一升级下电脑，有点想换amd，不过还在犹豫

@史浩 我不是大佬。最近在做一个VAR的项目，考虑温度场-电磁场-流场-组织耦合计算，需要自编程序，有程序验证的问题。对OpenFoam，我没做过这方面的案例，但我想考察2个方面：一是看一下OpenFoam与自编程序的结果对比，二是看一下OpenFoam与自编程序的计算速度对比。
ESR和VAR比较类似，有相通的地方，你如果做了OpenFoam算例的话可以交流一下。

@史浩 在 interFoam动网格计算问题计算的问题 中说：

简化为slip壁面，电极下表面简化为一个液态金属的入口，

你这个不考虑凝固和电场？
我最近也在做类似的项目，需要考虑凝固、电场、流场耦合

xps pression 5510 ，感觉速度还可以，能再用两年没问题

没人有类似的需求？

@队长别开枪 请教一个信息遍历问题：cellEdges输出信息如下：

第一行为单元个数，接下来是每个单元含有的edge数及edge的具体编号。
请教一下如何遍历所有单元，获得每个单元包含的edge？
我知道遍历patch获得其开始编号和大小可以按照下述代码：

  const polyMesh& patches = mesh.boundaryMesh();

  forAll(patches, patchi)
    {
      const polyPatch& pp = patches[patchi];

      Info << patchi << '\t' << pp.start() << '\t' << pp.size() << endl;

    }

我自己参照utility写了一个，运行出错：

  const polyMesh& mesh;
  const labelList& cellLabels;
  forAll(cellLabels, i)
    {
        label celli = cellLabels[i];

        const labelList& cEdges = mesh.cellEdges()[celli];

        forAll(cEdges, cEdgeI)
        {
            const edge& e = mesh.edges()[cEdges[cEdgeI]];

            Info<<e[0]<<e[1]<<nl;
        }
    }

请教如何遍历单元获得单元包含edge？

我输出cellPoints，发现六面体单元的节点编号时按照从小到大的顺序输出的：

其顺序与用户指南上的顺序不一致：

请教如何才能输出与上图对应的节点顺序的cellPoints文件 ？

@队长别开枪
也就是说并没有现成的函数实现查找当前面的相邻面？
如果没有的话只好自己编程实现了

@东岳 能抽空看一下本帖，并给出一个例子吗？谢谢！

primitiveMesh.h文件里面的checkFaceFaces 函数，
我看了一下不知道怎么用，请问能给一个例子吗？
比如我在一个六面体单元中，共有6个面，当前面有4个相邻面，有1个向对面，怎么用这个函数把4个相邻面找出来？谢谢！

计算完毕后，规定任意一条直线及点数，如何插值求出直线上任意一点的值？

特别是计算区域分为几个小的区域，每个区域内有不同的物性参数的情况如何计算？

我这个是一个圆柱，分成1/26，成为周期对称模型。中间是空气，不能省略。

我也是觉得尖角输入，无法更改。所以想把扭曲的四边形对半分成2个三角形，1个hexahedron变成2个prism。 卡在不会分成2个prism了！

我这个角度很小，只有13°。Y型网格对角度很大的划分网格质量不错，但是角度很小的话，质量很差

分析问题可以用周期网格，master、slave边界。我在计算区域角部切出一个block用Y型网格划分。


但是角部网格为菱形，变形严重。我想把菱形棱柱网格分为2个三棱柱网格，示意图如下：

这个如何操作？

@Wayne 好的，多谢，我试试。