[关闭]
@daidezhi 2016-07-01T07:50:52.000000Z 字数 4103 阅读 4098

经典算例：顶盖驱动流

CFD OpenFOAM

1. 问题描述

方腔内部等温不可压流动算例常用于测试不可压流动求解器的计算准确性，测试数据通常使用Ghia等人[1]的数据。如下图所示，方腔左表面、右表面和底部均为静止壁面，顶部为运动壁面，速度为。

cavity

这里使用OpenFOAM标准不可压求解器icoFoam对雷诺数的工况进行数值模拟，并比较了upwind、linearUpwind和vanLeer三种对流格式的精度。

2. 求解器设置

2.1 网格

这里使用的网格为适用于OpenFOAM二维模拟的多边形棱柱体网格，如下图所示

square_OpenFOAM_dual_mesh

关于多边形棱柱体网格的生成可参考教程OpenFOAM二维多面体网格准备。

2.2 字典文件设置

本教程算例文件结构如下：

$case 算例根目录
├──constant 网格和输运参数目录
│ ├──polyMesh 网格目录，网格格式详见CFD多面体网格数据结构—OpenFOAM
│ └──transportProperties 输运参数设置
├──0 初始和边界条件目录
│ ├──p 变量p初始和边界条件设置
│ └──U 速度U初始和边界条件设置
└──system 求解器设置目录
├──controlDict 计算控制参数设置，设置起始终止时间，时间步长，输出控制等
├──fvSchemes 微分算子离散格式设置
├──fvSolution 代数方程组求解和PISO算法设置
└──sampleDict 数据采样设置

下面按照上面文件出现的次序依次列出其具体内容（略去文件头部说明）。

2.2.1 transportProperties

transportProperties字典文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       dictionary;             //类型为字典文件
  7.     location    "constant";             //所在目录
  8.     object      transportProperties;    //对象名
  9. }
  11. //设置运动粘性系数nu的单位和值
  12. nu              nu [ 0 2 -1 0 0 0 0 ] 0.001;

2.2.2 p

p文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       volScalarField;         //类型为体心标量场
  7.     object      p;                      //对象名
  8. }
  10. dimensions      [0 2 -2 0 0 0 0];       //设置单位
  12. internalField   uniform 0;              //内场初始化为0
  14. //边界条件
  15. boundaryField
  16. {
  17.     leftWall
  18.     {
  19.         type            zeroGradient;
  20.     }
  21.     rightWall
  22.     {
  23.         type            zeroGradient;
  24.     }
  25.     lowerWall
  26.     {
  27.         type            zeroGradient;
  28.     }
  29.     topWall
  30.     {
  31.         type            zeroGradient;
  32.     }
  33.     defaultFaces
  34.     {
  35.         type            empty;
  36.     }
  37. }

2.2.3 U

U文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       volVectorField;         //类型为体心矢量场
  7.     object      U;                      //对象名
  8. }
  10. dimensions      [0 1 -1 0 0 0 0];       //设置单位
  12. internalField   uniform (0 0 0);        //内场初始化为0矢量
  14. //边界条件
  15. boundaryField
  16. {
  17.     leftWall
  18.     {
  19.         type            fixedValue;
  20.         value           uniform (0 0 0);
  21.     }
  22.     rightWall
  23.     {
  24.         type            fixedValue;
  25.         value           uniform (0 0 0);
  26.     }
  27.     lowerWall
  28.     {
  29.         type            fixedValue;
  30.         value           uniform (0 0 0);
  31.     }
  32.     topWall
  33.     {
  34.         type            fixedValue;
  35.         value           uniform (1 0 0);
  36.     }
  37.     defaultFaces
  38.     {
  39.         type            empty;
  40.     }
  41. }

2.2.4 controlDict

controlDict字典文件的内容如下

  1. FoamFile
  2. {
  3.     version     2.0;                //版本号
  4.     format      ascii;              //文本格式
  5.     class       dictionary;         //类型为字典文件
  6.     location    "system";           //所在目录
  7.     object      controlDict;        //对象名
  8. }
  10. application     icoFoam;
  12. startFrom       startTime;
  14. startTime       0;
  16. stopAt          endTime;
  18. endTime         50.0;
  20. deltaT          0.01;
  22. writeControl    runTime;
  24. writeInterval   1.0;
  26. purgeWrite      0;
  28. writeFormat     ascii;
  30. writePrecision  6;
  32. writeCompression off;
  34. timeFormat      general;
  36. timePrecision   6;
  38. runTimeModifiable true;
  40. //设置探针类型和位置,用于判断流动是否达到稳态
  41. functions
  42. {
  43.     UProbe
  44.     {
  45.         type                    probes;
  46.         functionObjectLibs      ("libsampling.so");
  47.         outputControl           outputTime;
  48.         probeLocations
  49.         (
  50.             (0.5 0.5 0.125)
  51.         );
  52.         fields
  53.         (
  54.             U
  55.         );
  56.     }
  58.     pProbe
  59.     {
  60.         type                    probes;
  61.         functionObjectLibs      ("libsampling.so");
  62.         outputControl           outputTime;
  63.         probeLocations
  64.         (
  65.             (0.5 0.5 0.125)
  66.         );
  67.         fields
  68.         (
  69.             p
  70.         );
  71.     }
  72. }

2.2.5 fvSchemes

fvSchemes字典文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       dictionary;             //类型为字典文件
  7.     location    "system";               //所在目录
  8.     object      fvSchemes;              //对象名
  9. }
  11. ddtSchemes
  12. {
  13.     default         Euler;
  14. }
  16. gradSchemes
  17. {
  18.     default         Gauss linear;
  19.     grad(p)         Gauss linear;
  20. }
  22. divSchemes
  23. {
  24.     default         none;
  25.     div(phi,U)      Gauss upwind phi;   //之后将"upwind"替换为"linearUpwind"和"vanLeer"
  26. }
  28. laplacianSchemes
  29. {
  30.     default         Gauss linear corrected;
  31. }
  33. interpolationSchemes
  34. {
  35.     default         linear;
  36. }
  38. snGradSchemes
  39. {
  40.     default         corrected;
  41. }
  43. fluxRequired
  44. {
  45.     default         no;
  46.     p               ;
  47. }

2.2.6 fvSolution

fvSolution字典文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       dictionary;             //类型为字典文件
  7.     location    "system";               //所在目录
  8.     object      fvSolution;             //对象名
  9. }
  11. solvers
  12. {
  13.     p
  14.     {
  15.         solver          PCG;
  16.         preconditioner  DIC;
  17.         tolerance       1e-08;
  18.         relTol          0;
  19.     }
  21.     U
  22.     {
  23.         solver          smoothSolver;
  24.         smoother        symGaussSeidel;
  25.         tolerance       1e-08;
  26.         relTol          0;
  27.     }
  28. }
  30. PISO
  31. {
  32.     nCorrectors     5;
  33.     nNonOrthogonalCorrectors 5;
  34.     pRefCell        0;
  35.     pRefValue       0;
  36. }

2.2.7 sampleDict

sampleDict字典文件的内容如下

  1. //文件说明
  2. FoamFile
  3. {
  4.     version     2.0;                    //版本号
  5.     format      ascii;                  //文本格式
  6.     class       dictionary;             //类型为字典文件
  7.     location    "system";               //所在目录
  8.     object      sampleDict;             //对象名
  9. }
  11. interpolationScheme cellPoint;          //设置插值方式
  13. setFormat raw;                          //常规文本输出
  15. //采样集合
  16. sets
  17. (
  18.     //垂直中心线
  19.     verticalCenterLine
  20.     {
  21.         type    uniform;
  22.         axis    y;
  23.         start   (0.5 0.0 0.25);
  24.         end     (0.5 1.0 0.25);
  25.         nPoints 100;
  26.     }
  27.     //水平中心线
  28.     horizontalCenterLine
  29.     {
  30.         type    uniform;
  31.         axis    x;
  32.         start   (0.0 0.5 0.25);
  33.         end     (1.0 0.5 0.25);
  34.         nPoints 100;
  35.     }
  36. );
  38. //采样对象
  39. fields
  40. (
  41.     U
  42. );

3. 计算结果

3.1 稳态判断

由于icoFoam求解器属于非定常求解器，因此需要按时间推进模拟顶盖驱动流，到达一定时间后流动趋于稳定。用来对比验证的数据是稳态数据，因此需要确定何时流动能够趋于稳定。

下面三幅图分别是压力p，速度U的x和y分量随时间的变化曲线，从图中可以看出流动大概在后趋于稳定。

time_history_of_kinematic_pressure_at_centroid 探针位置p随时间变化曲线

time_history_of_x_velocity_at_centroid 探针位置U.x()随时间变化曲线

time_history_of_y_velocity_at_centroid 探针位置U.y()随时间变化曲线

3.2 upwind对流格式

使用upwind对流格式得到的计算结果与基准解的对比如下所示

x_velocity_along_vertical_central_line 垂直中线上U.x()变化曲线

y_velocity_along_horizontal_central_line 水平中线上U.y()变化曲线

3.3 linearUpwind对流格式

使用linearUpwind对流格式得到的计算结果与基准解的对比如下所示

x_velocity_along_vertical_central_line 垂直中线上U.x()变化曲线

y_velocity_along_horizontal_central_line 水平中线上U.y()变化曲线

3.4 vanLeer对流格式

使用vanLeer对流格式得到的计算结果与基准解的对比如下所示

x_velocity_along_vertical_central_line 垂直中线上U.x()变化曲线

y_velocity_along_horizontal_central_line 水平中线上U.y()变化曲线

不难看出，计算结果和基准解的吻合度linearUpwind vanLeer upwind。

感谢您的阅读，欢迎讨论和批评指正。

作者：戴得志

[1] Ghia U, Ghia K N, Shin C T. High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method[J]. Journal of computational physics, 1982, 48(3): 387-411. ↩
添加新批注
在作者公开此批注前，只有你和作者可见。
回复批注