Wrong Boundary Conditions (maximum number of iterations) in compressibleInterDyMFoam

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

 September 11, 2017, 07:05 Wrong Boundary Conditions (maximum number of iterations) in compressibleInterDyMFoam #1 New Member   Oliver K Join Date: May 2017 Posts: 15 Rep Power: 9 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::T(scalar f, scalar T0, scalar (thermo::*F)(const scalar) const, scalar (thermo::*dFdT)(const scalar) const, scalar (thermo::*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.