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

interFoam blows up, caused by p_rgh

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 13, 2016, 08:00
Default interFoam blows up, caused by p_rgh
  #1
New Member
 
Christian Jähnel
Join Date: Nov 2016
Posts: 11
Rep Power: 9
ch_jaehnel is on a distinguished road
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;
    }
Attached Images
File Type: jpg initialconditon.JPG (26.0 KB, 54 views)
ch_jaehnel 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
k-e & GAMG interFoam Schemitisation Stability Issue JFM OpenFOAM Running, Solving & CFD 3 December 1, 2015 05:58
interFoam blows up with small channels and high pressure NiFl OpenFOAM 0 March 7, 2014 04:49
mpirun interFoam very stable --> then blows up ghadab OpenFOAM Running, Solving & CFD 3 October 27, 2013 10:34
Interfoam blows on parallel run danvica OpenFOAM Running, Solving & CFD 16 December 22, 2012 02:09
interFoam case blows up jrrygg OpenFOAM Running, Solving & CFD 8 November 14, 2012 03:16


All times are GMT -4. The time now is 22:29.