blockMesh-四分之一同心圆柱体网格

上一篇介绍了四分之一同心圆柱网格的画法，得到网格的分辨率大概是0.4mm，远远达不到DNS的标准，但是网格数量已经有21万，再加上详细的化学反应机理，一般的PC根本跑不动，不适合初期的算例调试，所以今天给大家分享一个适用于射流（主射流+伴流）的二维同心圆柱网格画法。

下面是具体的blockMeshDict文件

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

scale 0.001; //mm

vertices
(
  (0.0000000000000000   0.0000000000000000  0.0000000000000000) // vertex number 0
  (9.9965732497555700   0.0000000000000000  -0.2617694830787310)    // vertex number 1
  (10.9962305747311000  0.0000000000000000  -0.2879464313866050)    // vertex number 2
  (39.9862929990223000  0.0000000000000000  -1.0470779323149300)    // vertex number 3
  (39.9862929990223000  0.0000000000000000  1.0470779323149300) // vertex number 4
  (10.9962305747311000  0.0000000000000000  0.2879464313866050) // vertex number 5
  (9.9965732497555700   0.0000000000000000  0.2617694830787310) // vertex number 6
  (0.0000000000000000   80.0000000000000000 0.0000000000000000) // vertex number 7
  (9.9965732497555700   80.0000000000000000 -0.2617694830787310)    // vertex number 8
  (10.9962305747311000  80.0000000000000000 -0.2879464313866050)    // vertex number 9
  (39.9862929990223000  80.0000000000000000 -1.0470779323149300)    // vertex number 10
  (39.9862929990223000  80.0000000000000000 1.0470779323149300) // vertex number 11
  (10.9962305747311000  80.0000000000000000 0.2879464313866050) // vertex number 12
  (9.9965732497555700   80.0000000000000000 0.2617694830787310) // vertex number 13
);

edges
(
  arc 1 6 (10 0 0)
  arc 2 5 (11 0 0)
  arc 3 4 (40 0 0)
  arc 8 13 (10 80 0)
  arc 9 12 (11 80 0)
  arc 10 11 (40 80 0)

);

blocks
(
// mesh 0.2mm
  hex (0 1 8 7 0 6 13 7) (50 400 1) simpleGrading (1 1 1) // block 0
  hex (1 2 9 8 6 5 12 13) (5 400 1) simpleGrading (1 1 1) // block 1
  hex (2 3 10 9 5 4 11 12) (145 400 1) simpleGrading (1 1 1) // block 1
);


boundary
(
  PrimaryIinlet // 主射流入口
  {
    type patch;
    faces
    (
     (1 0 0 6)
    );
  }

 Front
  {
    type wedge;
    faces
    (
     (7 0 6 13)
     (13 6 5 12)
     (12 5 4 11)
    );
  }

  Back
  {
    type wedge;
    faces
    (
     (0 7 8 1)
     (1 8 9 2)
     (2 9 10 3)
    );
  }
   
  CoflowInlet // 伴流入口
  {
    type patch;
    faces
    (
     (2 3 4 5)
    );
  }

  Outlet1 // 侧面出口
  {
    type patch;
    faces
    (
     (4 3 10 11)
    );
  }

  Outlet2 // 顶部出口
  {
    type patch;
    faces
    (
     (8 7 7 13)
     (8 13 12 9)
     (9 12 11 10)
    );
  }
  Symmetry // 圆柱轴为对称边界
  {
    type empty;
    faces
    (
     (7 0 0 7)
    );
  }
  
  Wall // 主射流和伴流之间是一个厚度为1mm的壁面
  {
    type wall;
    faces
    (
     (1 2 5 6)
    );
  }
);

mergePatchPairs
(
);

// ************************************************************************* //

以下是成品图片（分辨率减小10倍之后）

Mesh

目前的网格分辨率是0.2mm，网格数量是8万，调试算例还是可以跑得动。

关于二维和三维圆柱网格计算结果的差异，如果有CFDer研究过，欢迎在下方留言讨论