05、【算例】openFoam盖驱动空腔流动

管网：https://doc.cfd.direct/openfoam/user-guide-v9/cavity
一、算例实现

文件结构

• 0：存放初场
  
• constant：存放网格信息
  
• system：存放网格划分、计算等工具
  

1、画网格

1. blockMesh

复制代码

2、求解

1. icoFoam

复制代码

3、保存文件

1. touch cavity.OpenFOAM

复制代码

4、后处理

1. paraview

复制代码

使用openFoamReader打开

显示所有边界

apply

二、网格加密

1、网格划分

在icoFoam文件夹内创建cavityFine文件夹，复制cavity文件夹内的constant和system文件夹至该文件夹。
在system/blockMeshDict文件改变网格数量，由之前的（20,20,1）改为（40,40,1），之后使用blockMesh生成网格信息

2、mapFields映射结果

将粗网格中的0.5结果，映射到细网格中。
修改controlDict文件

映射
mapFields ../cavity/cavity -consistent

3、计算并输出log文件

1. // 计算并输出log文件
2. icoFoam > log
3. // 查看log文件
4. cat log
6. // 可以合并运行
7. icoFoam > log && cat log

复制代码

4、calcType

可以对速度的某个量进行可视化处理并咬着计算域内的一条线来绘制一个二维的图

1. postProcess -func "components(U)"

复制代码

此时会在时间步长文件夹下出现3个新文件

5、后处理

使用paraview加载数据，就会显示刚才计算出的值。

三、网格非均匀分布

1、修改blockMeshDict

在icoFoam文件夹内创建asMeshCavity文件夹，复制cavity文件夹内的constant和system文件夹至该文件夹。
修改constant/blockMeshDict文件,如下

1. /*--------------------------------*- C++ -*----------------------------------*\
2. =========                 |
3.   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4.    \\    /   O peration     | Website:  https://openfoam.org
5.     \\  /    A nd           | Version:  9
6.      \\/     M anipulation  |
7. \*---------------------------------------------------------------------------*/
8. FoamFile
9. {
10.     format      ascii;
11.     class       dictionary;
12.     object      blockMeshDict;
13. }
14. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16. convertToMeters 0.1;
18. vertices
19. (
20.     (0 0 0)
21.     (0.5 0 0)
22.     (1 0 0)
23.     (0 0.5 0)
24.     (0.5 0.5 0)
25.     (1 0.5 0)
26.     (0 1 0)
27.     (0.5 1 0)
28.     (1 1 0)
29.     (0 0 0.1)
30.     (0.5 0 0.1)
31.     (1 0 0.1)
32.     (0 0.5 0.1)
33.     (0.5 0.5 0.1)
34.     (1 0.5 0.1)
35.     (0 1 0.1)
36.     (0.5 1 0.1)
37.     (1 1 0.1)
38. );
40. blocks
41. (
42.     hex (0 1 4 3 9 10 13 12) (10 10 1) simpleGrading (2 2 1)
43.     hex (1 2 5 4 10 11 14 13) (10 10 1) simpleGrading (0.5 2 1)
44.     hex (3 4 7 6 12 13 16 15) (10 10 1) simpleGrading (2 0.5 1)
45.     hex (4 5 8 7 13 14 17 16) (10 10 1) simpleGrading (0.5 0.5 1)
46. );
48. edges
49. (
50. );
52. boundary
53. (
54.     movingWall
55.     {
56.         type wall;
57.         faces
58.         (
59.             (6 15 16 7)
60.             (7 16 17 8)
61.         );
62.     }
63.     fixedWalls
64.     {
65.         type wall;
66.         faces
67.         (
68.             (3 12 15 6)
69.             (0 9 12 3)
70.             (0 1 10 9)
71.             (1 2 11 10)
72.             (2 5 14 11)
73.             (5 8 17 14)
74.         );
75.     }
76.     frontAndBack
77.     {
78.         type empty;
79.         faces
80.         (
81.             (0 3 4 1)
82.             (1 4 5 2)
83.             (3 6 7 4)
84.             (4 7 8 5)
85.             (9 10 13 12)
86.             (10 11 14 13)
87.             (12 13 16 15)
88.             (13 14 17 16)
89.         );
90.     }
91. );
93. mergePatchPairs
94. (
95. );
97. // ************************************************************************* //

复制代码

2、画网格

1. blockMesh

复制代码

3、生成asyMeshCavity.OpenFOAM文件

1. touch asyMeshCavity.OpenFOAM

复制代码

4、后处理

1. paraview

复制代码

5、修改时间步长



6、映射场

将cavityFine中的0.7文件映射过来

1. mapFields -consistent ../cavityFine

复制代码

7、计算输出

1. icoFoam > log

复制代码

四、大雷诺数

1、创建文件夹并拷贝

创建文件夹

1. mkdir cavityHighRe

复制代码

拷贝cavity文件夹至该文件夹

1. cp -r ../cavity/cavity/ .

复制代码

2、修改雷诺数

修改constant/transportProperties文件

即可提高雷诺数10倍
3、修改控制文件

从0.5开始计算到2s

4、计算

1. nohup nice -19 icoFoam > log && cat log

复制代码

• nohup：党用户退出登录时，程序依然执行；
  
• nice：调整程序优先级，-20对应优先级最高，19对应优先级最低进程。
  

可以看出1.4s时Ux结束迭代，No Iterations 0表示速度求解停止。

五、RAS模型

进入pisoFoam/RAS/cavity

可以看到在0文件夹下有许多有关湍流模型的参数

在momentumTransport文件里定义模型

1. /*--------------------------------*- C++ -*----------------------------------*\
2.   =========                 |
3.   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4.    \\    /   O peration     | Website:  https://openfoam.org
5.     \\  /    A nd           | Version:  9
6.      \\/     M anipulation  |
7. \*---------------------------------------------------------------------------*/
8. FoamFile
9. {
10.     format      ascii;
11.     class       dictionary;
12.     location    "constant";
13.     object      momentumTransport;
14. }
15. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17. simulationType  RAS;
19. RAS
20. {
21.     model           kEpsilon;
23.     turbulence      on;
25.     printCoeffs     on;
26. }
28. // ************************************************************************* //

复制代码

simulationType：laminar, RAS and LES
printCoeffs：为on时，这些参数会被输出到终端
六、CavityClipped

1、划分网格

1. blockMesh

复制代码

2、复制0文件夹内容至0.5

1. cp -r 0 0.5

复制代码

3、设置

controlDict文件开始时间设置为0.5，结束时间设置为0.6
system/mapFieldsDict文件
有两个参数

1. patchMap        (lid movingWall);
3. cuttingPatches  ();

复制代码

• 当用户向把原始场的patches条件投影到背投影场的patches的时候，我们使用patchMap。
  
• 当用户打算把原始场内的数据投影到被投影场的边界时，我们使用cuttingPatches列表
  

这里lib边界条件根cavity算例的movingWall边界条件是一致的。

4、计算

1. icoFoam > log

复制代码

5、后处理


修改样式


修改Legend

绘制矢量图



免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！