如何理解openfoam案例里面的blockMesh文件里面的simpleGrading
总结：
simpleGrading参数分为xyz三个方向。如果你想使得网格在某个方向上更密集，可以在simpleGrading中将该方向的渐变率设置为小于 1 .更稀疏则设置大于1.
一、案例
比如我这个爆炸案例：
对应的blockMeshDIct文件如下：
// 定义划分的网络数目
u1 20;//100
u2 4;//20
u3 16; //80
v1 12;//60
v2 20;//100
w 16; //16;//160 // 3D d*32
//w 1; // 2D// 定义划分的网络渐变率，分为xyz三个方向。如果你想使得网格在某个方向上更密集，可以在simpleGrading中将该方向的渐变率设置为小于 1 .
r1 0.25;
r2 4; //((0.5 0.5 $r1) (0.5 0.5 $r2)) 表示Z方向中间密集两遍稀疏，因为爆炸点在最中间//blocks中包含了五个长方体块，其中hex是一个构成长方体块的8个顶点，($u1 $v1 $w)代表xyz三个方向的网格，simpleGrading代表三个方向划分网格的渐变率。可以某个方向网格逐渐变密（值[0,1]）或者变稀(值[1,无穷])
blocks
(hex (0 1 5 4 12 13 17 16)($u1 $v1 $w)simpleGrading (0.4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2))) hex (2 3 7 6 14 15 19 18)($u3 $v1 $w)simpleGrading (4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2)))hex (4 5 9 8 16 17 21 20)($u1 $v2 $w)simpleGrading (0.4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2)))hex (5 6 10 9 17 18 22 21)($u2 $v2 $w)simpleGrading (1 3 ((0.5 0.5 $r1)(0.5 0.5 $r2)))hex (6 7 11 10 18 19 23 22)($u3 $v2 $w)simpleGrading (4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2)))
);
解析：
u,v,w：代表三个方向划分的网络数目
r1,r2:定义划分的网络渐变率。
blocks：中包含了五个长方体块，其中hex是一个构成长方体块的8个顶点，($u1 $v1 $w)代表xyz三个方向的网格，simpleGrading代表三个方向划分网格的渐变率。可以某个方向网格逐渐变密（值[0,1]）或者变稀(值[1,无穷])
simpleGrading详解：
例子：simpleGrading (0.4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2)))
simpleGrading参数分为xyz三个方向。如果你想使得网格在某个方向上更密集，可以在simpleGrading中将该方向的渐变率设置为小于 1 .更稀疏则设置大于1.
simpleGrading (0.4 3 ((0.5 0.5 $r1)(0.5 0.5 $r2))) 是 OpenFOAM 中定义网格渐变率的语法。它定义了每个方向上网格单元的变化比例。具体来说：
- 0.4：沿 x 方向的网格渐变率。
- 3：沿 y 方向的网格渐变率。
- ((0.5 0.5 $r1)(0.5 0.5 $r2))：沿 z 方向的分段渐变率。
详细解析
0.4沿 x 方向：网格单元从左到右逐渐变小，密度更高。3沿 y 方向：网格单元从下到上逐渐变大，密度更低。((0.5 0.5 $r1)(0.5 0.5 $r2)))沿 z 方向：- 第一段占长度的一半，均匀分布，初始和末端的渐变率为 0.25（单个网格逐渐变小）。
- 第二段占长度的一半，均匀分布，初始和末端的渐变率为 4（单个网格逐渐变大）。
- 总结就是Z方向上两边网格稀疏，中间网格密集。（因为爆炸点在Z方向中间，需要希望格捕捉细小变化）
可以看出符合预期，越靠近中间的弹性挡板和爆炸点，xyz三个方向的网格就越密。