ch_jaehnel |
December 13, 2016 08:00 |
interFoam blows up, caused by p_rgh
1 Attachment(s)
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;
}
|