[silencebreak]

Hey there,

I've got a problem with running the case in compressible Mode. The case is running without any errors in incompressible, so I don't think the error could be in the mesh therefore only the BCs of the compressible solver could be chosen wrong.

The case:
The case consists an upper reservoir (kaverneWall, atmosphere) a pipe system with a valve (ACMI-patches, AMI, slide_interface) at the bottom which is closing. The water is set at a level in the reservoir which is flooding down the pipe. The case is going to be calculated in turbulence mode but for the start i tried the BCs in laminar mode to work.

The problem:
The pressure is rising and it occurs an error regarding the maximum limitation of iterations are exceeded.
The error looks like:

max(U) 799.31
min(p_rgh) -640853
PIMPLE: iteration 2
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = 0.148083  Min(alpha.water) = -8.56923  Min(alpha.air) = -11.9184
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = 0.148095  Min(alpha.water) = -1851.18  Min(alpha.air) = -50816.1
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = 0.255911  Min(alpha.water) = -7.80007e+06  Min(alpha.air) = -1.29547e+09
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = 46.493  Min(alpha.water) = -9.70012e+11  Min(alpha.air) = -6.11008e+11
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = 46.493  Min(alpha.water) = -9.70012e+11  Min(alpha.air) = -6.11008e+11
MULES: Solving for alpha.water
MULES: Solving for alpha.water
Liquid phase volume fraction = -3.36312e+21  Min(alpha.water) = -4.89659e+31  Min(alpha.air) = -3.17104e+24
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for T, Initial residual = 1, Final residual = 0.331259, No Iterations 1001


--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

alpha.water

atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    kaverneWall
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    pipeWall
    {
        type            zeroGradient;
    }
    AMI
    {
        type            zeroGradient;
    }
    slide_interface
    {
        type            zeroGradient;
    }
    ACMI1_couple
    {
        type            cyclicACMI;
        value           uniform 0;
    }
    ACMI1_blockage
    {
        type            zeroGradient;
    }
    ACMI2_couple
    {
        type            cyclicACMI;
        value           uniform 0;
    }
    ACMI2_blockage
    {
        type             zeroGradient;
    }

p

internalField   uniform 101325;

boundaryField
{
    atmosphere
    {
        type            calculated;
        value           $internalField;
    }
    kaverneWall
    {
        type            calculated;
        value           $internalField;
        //type            zeroGradient;
    }
    outlet
    {
        type            calculated;//was fixedValue
        value           $internalField;
    }
    pipeWall
    {
        type            calculated;
        value           $internalField;
        //type            zeroGradient;
    }
    AMI
    {
        type            calculated;
        value           $internalField;
        //type            zeroGradient;
    }
    slide_interface
    {
        type            calculated;
        value           $internalField;
        //type            zeroGradient;
    }
    ACMI1_couple
    {
        type            cyclicACMI;
        value           $internalField;
    }
    ACMI1_blockage
    {
        type            calculated;
        value           $internalField;
    }
    ACMI2_couple
    {
        type            cyclicACMI;
        value           $internalField;
    }
    ACMI2_blockage
    {
        type            calculated;
        value           $internalField;
    }
}

p_rgh

internalField   uniform 101325;
atmosphere
    { type            totalPressure;
        rho             rho;
        psi             none;
        gamma           1;
        p0              uniform 0;
        value           uniform 101325;
    }
    kaverneWall
    {
        type            fixedFluxPressure;//was zeroGradient latest fixedFluxPressure
        gradient        uniform 0;
        value           $internalField;
    }
outlet
    {
        type            zeroGradient;
        //type            fixedValue;
        //value           $internalField;
    }
    pipeWall
    {
        type            fixedFluxPressure;//was zeroGradient latest fixedFluxPressure
        gradient        uniform 0;
        value           $internalField;
    }
    AMI
    {
        type            fixedFluxPressure;//was zeroGradient latest fixedFluxPressure
        gradient        uniform 0;
        value           $internalField;
    }
    slide_interface
    {
        type            fixedFluxPressure;//was zeroGradient latest fixedFluxPressure
        gradient        uniform 0;
        value           $internalField;
    }
    ACMI1_couple
    {
        type            cyclicACMI;
        value           $internalField;//was uniform 0
    }
    ACMI1_blockage
    {
        type            fixedValue;
        value           $internalField;
        //type            zeroGradient;//was calculated
        //value           uniform 0;
    }
    ACMI2_couple
    {
        type            cyclicACMI;
        value           $internalField;//was uniform 0
    }
ACMI2_blockage
    {
        type            fixedFluxPressure;//was calculated
        gradient        uniform 0;
        value           $internalField;
        //value           uniform 0;
    }

T:

internalField   uniform 293.15;

boundaryField
{
    atmosphere
    {
        type            inletOutlet;//was zeroGradient
        inletValue      $internalField;
        value           $internalField;
    }
    kaverneWall
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
        //type            inletOutlet;//was zeroGradient
        //inletValue      uniform 293.15;
        //value           $internalField;

    }
pipeWall
    {
        type            zeroGradient;
    }
    AMI
    {
        type            zeroGradient;
    }
    slide_interface
    {
        type            zeroGradient;
    }
    ACMI1_couple
    {
        type            cyclicACMI;
        value           $internalField;
    }
    ACMI1_blockage
    {
        type            zeroGradient;//was zeroGradient latest inletOutlet
        //inletValue      uniform 293.15;
        //value           $internalField;
       // type            calculated;
       // value           uniform 0;
    }
    ACMI2_couple
    {
        type            cyclicACMI;
        value           $internalField;
    }
    ACMI2_blockage
    {
        type            zeroGradient;
        //type            calculated;
        //value           uniform 0;
    }
}

U:

internalField   uniform (0 0 0);

boundaryField
{
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
    kaverneWall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    pipeWall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    AMI
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    slide_interface
    {
        type            fixedValue;
        value           uniform 0;
    }
    ACMI1_couple
    {
        type            cyclicACMI;
        value           uniform (0 0 0);
    }
    ACMI1_blockage
    {
        type            zeroGradient;
        //type            fixedValue;//was calculated
        //value           uniform (0 0 0);
    }
 ACMI2_couple
    {
        type            cyclicACMI;
        value           uniform (0 0 0);
    }
    ACMI2_blockage
    {
        type            fixedValue;//was calculated
        value           uniform (0 0 0);
    }
}

Any suggestions where i went wrong with my boundary conditions or could the error be somewhere else?

I'm happy about any help!

Cheers,
Oli