09、openfoam案例之圆柱绕流

1、原视频地址

https://www.bilibili.com/video/BV1ME411A73k/?spm_id_from=333.1007.top_right_bar_window_custom_collection.content.click&vd_source=33b50a4dd201d7564e6e63d321809ce9
2、网格划分及导入

2.1 网格划分

本案例使用ICEM划分网格，并导入openfoam中

2.2 网格转换

目前通过在 3 维中定义网格来处理 2 维几何，其中前平面和后平面定义为空边界块类型。读取二维 Fluent 网格时，转换器会自动在第三方向拉伸网格并添加空面片，将其命名为 frontAndBackPlanes。

fluentMeshToFoam：读取fluent.msh网格文件。（指南的5.5章）
命令：

1. fluentMeshToFoam <meshFile>

复制代码

2.3 检查网格

1. checkMesh

复制代码

system/blockMeshDict文件可以删除

3、边界条件

3.1 网格boundary

这是转换网格后自动生成的。

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       polyBoundaryMesh;
12.     location    "constant/polyMesh";
13.     object      boundary;
14. }
15. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17. 5
18. (
19.     INLET
20.     {
21.         type            patch;
22.         nFaces          30;
23.         startFace       13060;
24.     }
25.     OUTLET
26.     {
27.         type            patch;
28.         nFaces          30;
29.         startFace       13090;
30.     }
31.     WALL
32.     {
33.         type            wall;
34.         inGroups        List<word> 1(wall);
35.         nFaces          100;
36.         startFace       13120;
37.     }
38.     CYLINDER
39.     {
40.         type            wall;
41.         inGroups        List<word> 1(wall);
42.         nFaces          120;
43.         startFace       13220;
44.     }
45.     frontAndBackPlanes // 空，表示2维
46.     {
47.         type            empty;
48.         inGroups        List<word> 1(empty);
49.         nFaces          13200;
50.         startFace       13340;
51.     }
52. )
54. // ************************************************************************* //

复制代码

3.2 0/P文件

• 在OpenFOAM中，zeroGradient是一种边界条件，用于指定在边界上的量的梯度为零，这意味着该边界的值不会改变。例如，对于速度边界条件，zeroGradient表示速度沿该边界的梯度为零，即速度沿该边界的值不会改变。通常，这种边界条件在流体边界和对称边界上使用。在代码中，它表示为：。
  
• 在 OpenFOAM 中，fixedValue 是一种边界条件类型，表示在特定的边界上施加固定的值。这个值通常是在 0、uniform、nonuniform 或 list 中指定的，具体取决于场的类型和边界的特征。例如，在一个流体模拟中，可以将固定的速度值施加在流场的边界上，这可以通过设置边界条件类型为 fixedValue 并指定速度值来实现。同样地，可以将固定的温度、压力或其他物理量值应用于相应的场变量上。在 OpenFOAM 中，边界条件的设置是在 case 目录下的 0 文件夹中的 boundary 文件中完成的。使用 fixedValue 边界条件时，用户需要指定固定的值并将其应用于适当的边界。
  
• 在OpenFOAM中，empty边界类型用于定义没有物理边界或者该边界不需要与周围区域相互作用的情况。例如，计算空气在室内流动的情况，室内墙壁可以定义为empty边界。在empty边界中，OpenFOAM不会对物理量进行任何修改，因此需要在边界条件的设置中进行额外的注意。在empty边界中，通常会使用fixedValue或者fixedGradient来设置边界条件。
  

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       volScalarField;
12.     object      p;
13. }
14. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16. dimensions      [0 2 -2 0 0 0 0];
18. internalField   uniform 0;
20. boundaryField
21. {
22.     INLET // 压力梯度为0
23.     {
24.         type            zeroGradient; // 压力梯度为0
25.     }
26.    
27.     WALL
28.     {
29.         type            zeroGradient;
30.     }
32.     CYLINDER
33.     {
34.         type            zeroGradient;
35.     }
36.    
37.     OUTLET // 出口压力为固定值0
38.     {
39.         type            fixedValue;
40.         value           uniform 0;
41.     }
43.     frontAndBack
44.     {
45.         type            empty;
46.     }
47. }
49. // ************************************************************************* //

复制代码

3.3 0/U文件

• internalField是指某个场量在整个计算域内的初始值。在进行流体流动计算时，internalField可以表示压力、速度、温度等场量的初始值。internalField可以被定义为标量、向量、对称张量和非对称张量场量。常见的初始化方法包括常数、分布式初始化和从其他场量插值得到。例如，对于速度场量，可以通过设置其为一个初始速度场量或者通过从初始条件的压力梯度中计算得到。一旦internalField被定义，将无法在后续的计算过程中修改。因此，需要根据计算需求仔细地定义initialField。
  

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       volVectorField;
12.     object      U;
13. }
14. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
16. dimensions      [0 1 -1 0 0 0 0];
18. // 场内初始值
19. internalField   uniform (50 0 0);
21. boundaryField
22. {
23.     INLET // 速度为固定值ux=50m/s
24.     {
25.         type            fixedValue;
26.         value           uniform (50 0 0);
27.     }
28.    
29.     WALL
30.     {
31.         type            fixedValue;
32.         value           uniform (50 0 0);
33.     }
34.    
35.     OUTLET
36.     {
37.         type            zeroGradient;
38.     }
40.     CYLINDER
41.     {
42.         type            fixedValue;
43.         value           uniform (0 0 0);
44.     }
46.     frontAndBack
47.     {
48.         type            empty;
49.     }
50. }
52. // ************************************************************************* //

复制代码

oepnfoam中zeroGradient和uniform (0 0 0)的区别：
zeroGradient只表示在边界上梯度不变，并不是不随时间变化。
uniform (0 0 0)更强硬一些，就是为某个值。

4、计算参数

4.1 controlDict

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    "system";
13.     object      controlDict;
14. }
15. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17. application     icoFoam;
19. startFrom       startTime;
21. startTime       0;
23. stopAt          endTime;
25. endTime         0.05;
27. deltaT          0.00001; // Co = U * t / x
29. writeControl    timeStep;
31. writeInterval   100;
33. purgeWrite      0;
35. writeFormat     ascii;
37. writePrecision  6;
39. writeCompression off;
41. timeFormat      general;
43. timePrecision   6;
45. runTimeModifiable true;
48. // ************************************************************************* //

复制代码

4.2 fvSchemes

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    "system";
13.     object      fvSchemes;
14. }
15. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17. ddtSchemes
18. {
19.     default         Euler;
20. }
22. gradSchemes
23. {
24.     default         Gauss linear;
25.     grad(p)         Gauss linear;
26. }
28. divSchemes
29. {
30.     default         none;
31.     div(phi,U)      Gauss linear;
32. }
34. laplacianSchemes
35. {
36.     default         Gauss linear orthogonal;
37. }
39. interpolationSchemes
40. {
41.     default         linear;
42. }
44. snGradSchemes
45. {
46.     default         orthogonal;
47. }
50. // ************************************************************************* //

复制代码

4.3 fvSolution

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    "system";
13.     object      fvSolution;
14. }
15. // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
17. solvers
18. {
19.     p
20.     {
21.         solver          PCG;
22.         preconditioner  DIC;
23.         tolerance       1e-06;
24.         relTol          0.05;
25.     }
27.     pFinal
28.     {
29.         $p;
30.         relTol          0;
31.     }
33.     U
34.     {
35.         solver          smoothSolver;
36.         smoother        symGaussSeidel;
37.         tolerance       1e-05;
38.         relTol          0;
39.     }
40. }
42. PISO
43. {
44.     nCorrectors     2;
45.     nNonOrthogonalCorrectors 0;
46.     pRefCell        0;
47.     pRefValue       0;
48. }
51. // ************************************************************************* //

复制代码

5、计算

1. icoFoam > log

复制代码

6、后处理

生成临时后处理文件

1. paraFoam -builtin

复制代码

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