使用OpenFOAM中的VOF模型仿真两相流

使用OpenFOAM中的VOF模型仿真两相流

OpenFOAM中的VOF(Volume of Fluid)方法是模拟两相流界面的常用方法，特别适用于不相溶流体的界面追踪。以下是使用OpenFOAM进行VOF两相流仿真的详细指南：

1. 选择求解器

OpenFOAM提供了几个基于VOF方法的求解器：

• interFoam：不可压缩、等温两相流
• multiphaseInterFoam：多相流扩展版本
• compressibleInterFoam：可压缩两相流
• interIsoFoam：使用等距表面重构的VOF方法

对于基本的两相流仿真，interFoam是最常用的选择。

2. 设置案例文件

典型的VOF案例文件结构如下：

case/
├── 0/
│   ├── alpha.water (或alpha1)  # 相分数场
│   ├── p_rgh       # 压力场
│   ├── U           # 速度场
│   └── ...
├── constant/
│   ├── transportProperties  # 传输属性
│   ├── turbulenceProperties # 湍流模型
│   └── ...
└── system/
    ├── controlDict  # 控制参数
    ├── fvSchemes    # 离散格式
    └── fvSolution   # 求解器设置

3. 关键文件配置

constant/transportProperties

phases (water air);  // 定义相的名称

water
{
    transportModel  Newtonian;
    nu              [0 2 -1 0 0 0 0] 1e-06;  // 水的运动粘度
    rho             [1 -3 0 0 0 0 0] 1000;   // 水的密度
}

air
{
    transportModel  Newtonian;
    nu              [0 2 -1 0 0 0 0] 1.48e-05; // 空气的运动粘度
    rho             [1 -3 0 0 0 0 0] 1;        // 空气的密度
}

sigma             [1 0 -2 0 0 0 0] 0.07;  // 表面张力系数

0/alpha.water (相分数场)

internalField   uniform 0;  // 0表示空气，1表示水

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1;  // 入口全为水
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    walls
    {
        type            zeroGradient;
    }
    ...
}

system/fvSchemes

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss vanLeer;  // 对alpha使用vanLeer格式很重要
    div(phi,k)      Gauss upwind;
    ...
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

system/fvSolution

solvers
{
    alpha
    {
        nAlphaCorr      2;
        nAlphaSubCycles 1;  // 可以增加以提高界面分辨率
        MULESCorr       yes;
        ...
    }
    
    p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-6;
        relTol          0.05;
    }
    
    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-6;
        relTol          0;
    }
}

PIMPLE
{
    nOuterCorrectors 3;
    nCorrectors     2;
    nNonOrthogonalCorrectors 1;
}

4. 运行仿真

使用以下命令运行仿真：

interFoam

对于并行计算：

decomposePar
mpirun -np 4 interFoam -parallel
reconstructPar

5. 后处理

可以使用ParaView查看结果，重点关注：

• alpha.water场（相分数）
• 速度场
• 压力场

在ParaView中，可以创建等值面(alpha=0.5)来可视化界面。

6. 提高VOF仿真质量的技巧

1. 网格分辨率：界面区域需要足够细的网格
2. 时间步长：使用自适应时间步长控制Courant数

   controlDict中设置：
   adjustTimeStep yes;
   maxCo 0.5;  // 通常0.2-0.5
   

3. 子循环：增加nAlphaSubCycles可以提高界面捕捉精度
4. 界面压缩：使用interfaceCompression参数控制界面锐度

   在constant/transportProperties中添加：
   interfaceCompression true;
   cAlpha 1;  // 压缩系数，通常1-2
   

7. 常见问题解决

• 界面过度扩散：尝试增加nAlphaSubCycles或使用MULES限制器
• 质量不守恒：检查边界条件，确保没有非物理的通量
• 数值振荡：尝试不同的对流格式或减小时间步长

通过以上设置和技巧，您可以在OpenFOAM中成功实现基于VOF方法的两相流仿真。