CompressibleInterIsoFoam边界条件设定

Rachel0096

各位老师，晚上好。
我根据一个interIsoFoam的case，想把他按照原来的工况在compressibleInterIsoFoam中重新跑一下，主要是模拟横射流的，示意图如下图所示，但是出现以下问题：

其中，p_rgh、U和alpha.water都是之前interIsoFoam的设置
0/p:

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

internalField   uniform 0;

boundaryField
{
    INLET_AIR
    {
        type        fixedFluxPressure;
        value		uniform 0;
    }

    INLET_LIQUID
    {
        type        fixedFluxPressure;
        value		uniform 0;
    }

    WALL
    {
        type        fixedFluxPressure;
        value		uniform 0;
    }

    SYMMETRY
    {
        type        symmetry;
    }
  
    OUTLET
    {
        type            totalPressure;
        p0              uniform 0;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        value           uniform 0;
    }
}

0/alpha.water

dimensions      [0 0 0 0 0 0 0];


internalField   uniform 0;

boundaryField
{
    INLET_AIR
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    INLET_LIQUID
    {
        type            fixedValue;
        value           uniform 1;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    SYMMETRY
    {
        type            symmetry;
    }
    WALL
    {
        type            zeroGradient;
    }
}

0/U

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
    INLET_AIR
    {
        type            fixedValue;
        value           uniform (10 0 0);
    }
                  
    INLET_LIQUID
    {
        type            fixedValue;
        value           uniform (0 0 1);
    }
                  
    WALL
    {
        type            noSlip;
    }
    
    SYMMETRY
    {
        type            symmetry;
    }
    
    OUTLET
    {

        type            inletOutlet;
        inletValue	    uniform (0 0 0);
        value           uniform (0 0 0);
    }
}


// ************************************************************************* //

我自己主要添加了p、T、thermophysicalProperties.air和thermophysicalProperties.water文件，主要设置如下
p

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

internalField   uniform 101325;

boundaryField
{
    INLET_AIR
    {
        type            calculated;
        value           $internalField;
    }

    INLET_LIQUID
    {
        type            calculated;
        value           $internalField;
    }

    WALL
    {
        type            calculated;
        value           $internalField;
    }

    SYMMETRY
    {
        type        symmetry;
    }
  
    OUTLET
    {
        type            inletOutlet;
        inletValue	    uniform 0;
        value           uniform 0;
    }
}

T

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    INLET_AIR
    {
        type        fixedValue;
        value	     uniform 300;
    }

    INLET_LIQUID
    {
        type        fixedValue;
        value	     uniform 300;
    }

    WALL
    {
       type            zeroGradient;
    }

    SYMMETRY
    {
        type        symmetry;
    }
  
    OUTLET
    {
        type            inletOutlet;
        inletValue	    uniform 300;
        value           uniform 300;
    }
}

thermophysicalProperties.water

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;//认为其不可压缩，或许采用perfectfluid也可以？
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        molWeight   18.0;
    }
    equationOfState
    {
        rho         1000;
    }
    thermodynamics
    {
        Cp          4195;
        Hf          0;
    }
    transport
    {
        mu          1e-3;
        Pr          7;
    }
}

thermophysicalProperties.air

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1013;
        Hf          0;
    }
    transport
    {
        mu          16.337e-06;
        Pr          0.7;
    }
}

想问一下我的这些设置是否存在不合理的地方，当出现报错的时候，有没有什么方法可以找到出错的原因（我每次总是百度+论坛），怎么能够查看自己的边界条件是否合理呢，以及在thermophysicalProperties中的Pmin应该根据什么来进行设置？
问题有点多，麻烦各位老师啦！

李东岳

老铁你这啥都算啊，前几天还记得你弄IBM。为啥不用不可压缩呢。变量越多越容易出问题。这个错误就是发散了。

Rachel0096

@李东岳李老师，我上个周刚刚被我导师换课题了，现在做横射流雾化因为后期想要做到燃烧，所以需要改到可压缩的求解器里面

Rachel0096

@Rachel0096 省流版就是 IBM他等不及了

李东岳

射流雾化这个用VOF算，网格量不小。不太好debug。你这个问题，只能看出来是发散了。我看你最好先用不可压缩整收敛，然后再上弱可压缩。

Rachel0096

好的好的，我去试试，谢谢李老师

CFD中文网

CompressibleInterIsoFoam边界条件设定