OpenFoam收缩扩张喷管（拉瓦尔喷管）边界条件的设置

简介

收缩扩张喷管（也成拉瓦尔喷管）广泛应用于火箭推进。将其流动特性定性描述如下：

（1）当入口流量较小时，不出现雍塞，流速先增大后减小，全程为亚声速。出口压力即为大气压，入口压力由出口压力决定。

（2）当入口流量较大时，出现雍塞，流速先增大到声速（喉部位置）然后继续加速到超声速。如果不出现气流分离现象，那么出口压力与大气压无关，而由入口压力决定。而入口压力由流量决定。

（3）特殊地，如果喷管的面积比（出口面积/喉部面积）较大，则有可能出现气流分离现象。此时射流核心区的压力由入口压力决定，回流区的压力由大气压力决定。

可见为了完成收缩扩张喷管的数值仿真，必须设置合理的边界条件。

接下来给出OpenFoam中边界条件的设置

设置方法

leftWall是喷管入口，rightWall是喷管出口，upWall是喷管壁面，downWall是中心对称轴。front和back是楔形面。使用可压缩气体求解器rhoPimpleFoam。

设置压力边界条件如下：

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 101325;

boundaryField
{
    leftWall
    {
        type            zeroGradient;
    }

    rightWall
    {
        type            waveTransmissive;
        field           p;
        psi             thermo:psi;
        gamma           1.4;
        fieldInf        101325;
        lInf            3;
        value           uniform 101325;
    }

    upWall
    {
        type            zeroGradient;
    }

    downWall
    {
        type            symmetryPlane;
    }

    front {type wedge;}
    back {type wedge;}
}

温度边界条件如下：

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    leftWall
    {
        type            fixedValue;
        value           uniform 3100;
    }

    rightWall
    {
        type            zeroGradient;
    }

    upWall
    {
        type            zeroGradient;
    }

    downWall
    {
        type            symmetryPlane;
    }

    front {type wedge;}
    back {type wedge;}
}

速度边界条件如下：

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    leftWall
    {
        type            fixedValue;
        value           uniform (6 0 0);
    }

    rightWall
    {
        type            zeroGradient;
    }

    upWall
    {
        type            noSlip;
    }

    downWall
    {
        type            symmetryPlane;
    }

    front {type wedge;}
    back {type wedge;}
}

结果

求解的具体过程不再展开，这里给出结果

可以看出，此工况下出现了气流分离。燃气由亚声速加速到超声速，经过激波之后减速到亚声速。出口处既有亚声速也有超声速气流，所以必须使用waveTransmissive压力边界条件，而不能使用固定压力出口。

waveTransmissive压力边界条件可以适用于各种工况下收缩扩张喷管的计算。