CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Wrong Boundary Conditions (maximum number of iterations) in compressibleInterDyMFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 11, 2017, 07:05
Default Wrong Boundary Conditions (maximum number of iterations) in compressibleInterDyMFoam
  #1
New Member
 
Oliver K
Join Date: May 2017
Posts: 15
Rep Power: 9
silencebreak is on a distinguished road
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:
Code:
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
Code:
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
Code:
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
Code:
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:
Code:
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:
Code:
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

Last edited by silencebreak; September 11, 2017 at 11:03.
silencebreak is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
courant number increases to rather large values 6863523 OpenFOAM Running, Solving & CFD 22 July 5, 2023 23:48
Decomposing meshes Tobi OpenFOAM Pre-Processing 22 February 24, 2023 09:23
Suppress twoPhaseEulerFoam energy AlmostSurelyRob OpenFOAM Running, Solving & CFD 33 September 25, 2018 17:45
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24


All times are GMT -4. The time now is 14:05.