interFoam blows up, caused by p_rgh

ch_jaehnel · December 13, 2016, 09:00

Hey,
I'm quite new here in this forum, but as I'm working on this problems since weeks and don't find a solution for my problem, I wanted to try it here.

I have a case of several bassins between an upper inletbassin (z=1.5m) to a lower outletbassin (z=0). It is an interFoam k-epsilon-modell. I tried my boundary conditions on a smaller, more simple modell, with two slots (delta z was 0.5m), which worked fine.

Now when I try to start my calculation it blows up. As I can see in the log-file, it has to be a problem with p_rgh, as the Number of Iterations rises rapidly in there, until it blows up.

Here is my logfile:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.0-665f1db4c1f1
Exec   : interFoam -parallel
Date   : Dec 13 2016
Time   : 13:50:03
Host   : "iwd-pool-linuxVB"
PID    : 9566
Case   : /home/student/Schreibtisch/Test3
nProcs : 4
Slaves : 
3
(
"iwd-pool-linuxVB.9567"
"iwd-pool-linuxVB.9568"
"iwd-pool-linuxVB.9569"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}


Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.322874e-11, No Iterations 537
time step continuity errors : sum local = 5.7149548e-16, global = -1.5843079e-17, cumulative = -1.5843079e-17
Courant Number mean: 0.00059570362 max: 0.031257202

Starting time loop

Courant Number mean: 0.00059570362 max: 0.031257202
Interface Courant Number mean: 0 max: 0
Time = 0.001

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077205142  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077202135  Min(alpha.water) = 0  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.036189843, No Iterations 6
time step continuity errors : sum local = 0.00021506936, global = -6.5651335e-09, cumulative = -6.5651335e-09
DICPCG:  Solving for p_rgh, Initial residual = 0.014877625, Final residual = 0.00070570922, No Iterations 49
time step continuity errors : sum local = 1.9685094e-05, global = 3.1296087e-06, cumulative = 3.1230436e-06
DICPCG:  Solving for p_rgh, Initial residual = 0.002576988, Final residual = 9.8808431e-08, No Iterations 263
time step continuity errors : sum local = 3.1365924e-09, global = 1.3125848e-10, cumulative = 3.1231749e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.0084387872, Final residual = 7.4673803e-09, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 1.850843e-09, No Iterations 7
ExecutionTime = 27.76 s  ClockTime = 28 s

Courant Number mean: 0.0010808219 max: 0.079187784
Interface Courant Number mean: 1.1182117e-07 max: 0.040518929
Time = 0.002

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077202142  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077202152  Min(alpha.water) = 0  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0049595596, Final residual = 0.0002463037, No Iterations 61
time step continuity errors : sum local = 9.9291397e-06, global = -1.987743e-06, cumulative = 1.1354318e-06
DICPCG:  Solving for p_rgh, Initial residual = 0.00058948804, Final residual = 2.8963494e-05, No Iterations 71
time step continuity errors : sum local = 1.1727486e-06, global = -1.0655696e-08, cumulative = 1.1247761e-06
DICPCG:  Solving for p_rgh, Initial residual = 0.00014898443, Final residual = 9.8067311e-08, No Iterations 186
time step continuity errors : sum local = 4.0036311e-09, global = -2.6948826e-10, cumulative = 1.1245066e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.00036108033, Final residual = 3.3546116e-09, No Iterations 4
bounding epsilon, min: -1.5371717e-05 max: 1.3778721 average: 0.003022388
smoothSolver:  Solving for k, Initial residual = 0.37398095, Final residual = 9.2001002e-09, No Iterations 10
ExecutionTime = 41.31 s  ClockTime = 41 s

Courant Number mean: 0.0011115087 max: 0.076620589
Interface Courant Number mean: 7.855838e-06 max: 0.069906882
Time = 0.003

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077202159  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.07720217  Min(alpha.water) = 0  Max(alpha.water) = 1.0000001
DICPCG:  Solving for p_rgh, Initial residual = 0.5625793, Final residual = 0.024324459, No Iterations 5
time step continuity errors : sum local = 0.0021578958, global = 0.00062811044, cumulative = 0.00062923495
DICPCG:  Solving for p_rgh, Initial residual = 0.99876805, Final residual = 0.045502165, No Iterations 7
time step continuity errors : sum local = 4.3216929, global = 2.5899022, cumulative = 2.5905314
DICPCG:  Solving for p_rgh, Initial residual = 0.99999669, Final residual = 9.5122869e-08, No Iterations 372
time step continuity errors : sum local = 5.6957605, global = 0.3613578, cumulative = 2.9518892
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 8.0703001e-09, No Iterations 21
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 7.6771196e-09, No Iterations 23
ExecutionTime = 57.83 s  ClockTime = 58 s

Courant Number mean: 1.3903213e+08 max: 6.1307199e+13
Interface Courant Number mean: 239.48244 max: 1.4235232e+08
Time = 0.004

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077202178  Min(alpha.water) = -4309546.6  Max(alpha.water) = 3625137
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.077066911  Min(alpha.water) = -1.5011323e+17  Max(alpha.water) = 4.5907267e+16
DICPCG:  Solving for p_rgh, Initial residual = 0.98769124, Final residual = 8587422, No Iterations 1001
time step continuity errors : sum local = 3.9060583e+14, global = 3456549.3, cumulative = 3456552.3
DICPCG:  Solving for p_rgh, Initial residual = 1.0884745e-11, Final residual = 1.0884745e-11, No Iterations 0
time step continuity errors : sum local = 1.9339287e+24, global = 72191674, cumulative = 75648227
DICPCG:  Solving for p_rgh, Initial residual = 7.0029875e-11, Final residual = 7.0029875e-11, No Iterations 0
time step continuity errors : sum local = 1.244244e+25, global = 4.01702e+24, cumulative = 4.01702e+24
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 9.9841113e-09, No Iterations 750
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 9.8686237e-09, No Iterations 27
bounding k, min: -58152392 max: 8.6966368e+48 average: 1.4095221e+43
ExecutionTime = 104.53 s  ClockTime = 106 s

Courant Number mean: 9.4159369e+24 max: 1.854573e+30
Interface Courant Number mean: 1.549583e+21 max: 4.6948108e+26
Time = 0.005

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = -5.3606365e+19  Min(alpha.water) = -1.4473739e+41  Max(alpha.water) = 4.3048548e+41
MULES: Solving for alpha.water
Phase-1 volume fraction = -7.3849356e+69  Min(alpha.water) = -6.7043022e+91  Max(alpha.water) = 7.7610196e+91
DICPCG:  Solving for p_rgh, Initial residual = 0.9999999, Final residual = 4.5334855e+23, No Iterations 1001
time step continuity errors : sum local = 4.0166775e+66, global = -7.963992e+60, cumulative = -7.963992e+60
DICPCG:  Solving for p_rgh, Initial residual = 6.5474258e-42, Final residual = 6.5474258e-42, No Iterations 0
time step continuity errors : sum local = 1.1375298e+120, global = 1.1282882e+120, cumulative = 1.1282882e+120
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 9568 on node iwd-pool-linuxVB exited on signal 8 (Floating point exception).

Does anybody has any experiences with such a case? I'm not quite shure, what h is in p_rgh? Could this be a problem of a reference pressure?
Here are my boundary conditions:

U:

Code:

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

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            variableHeightFlowRateInletVelocity;
        flowRate        0.3;
        alpha           alpha.water;
        value           uniform (0 0 0);
    }

    outlet.water
    {
        type            fixedValue;
        value            uniform (0 0.374158 0);
    }

    outlet.air
    {
        type            fixedValue;
        value            uniform (0 0 0);
    }

    top
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

    bassinright
    {
        type            fixedValue;
        value            uniform (0 0 0);
    }
    
    bassinleft
    {
        type            fixedValue;
        value            uniform (0 0 0);
    }
    
    bottom
    {
        type            fixedValue;
        value            uniform (0 0 0);
    }

p_rgh:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
     inlet
    {
        type           zeroGradient;
    }

    outlet.water
    {
       type           zeroGradient;
    }

    outlet.air
    {
        type           zeroGradient;
    }

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

    bassinright
    {
        type           zeroGradient;
    }
    
    bassinleft
    {
        type           zeroGradient;
    }
    
    bottom
    {
        type           zeroGradient;
    }

alpha.water:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            variableHeightFlowRate;
        lowerBound      0.0;
        upperBound      0.9;
        value           uniform 0;
    }

    outlet.water
    {
        type            zeroGradient;
    }

    outlet.air
    {
        type            zeroGradient;
    }

    top
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }

    bassinright
    {
        type            zeroGradient;
    }
    
    bassinleft
    {
        type            zeroGradient;
    }
    
    bottom
    {
        type            zeroGradient;
    }

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
k-e & GAMG interFoam Schemitisation Stability Issue	JFM	OpenFOAM Running, Solving & CFD	3	December 1, 2015 06:58
interFoam blows up with small channels and high pressure	NiFl	OpenFOAM	0	March 7, 2014 05:49
mpirun interFoam very stable --> then blows up	ghadab	OpenFOAM Running, Solving & CFD	3	October 27, 2013 11:34
Interfoam blows on parallel run	danvica	OpenFOAM Running, Solving & CFD	16	December 22, 2012 03:09
interFoam case blows up	jrrygg	OpenFOAM Running, Solving & CFD	8	November 14, 2012 04:16