Skip to content

PIMPLE Solver is acting weird #131

Description

@MarineAiry

Dear Bruno,

Unfortunately InterFoam exits without any error message or anything, and the pc turns itself off.
The last step saved was 3844, so the simulation stopped just before its completion (endTime = 4000).
The output at its beginning is as follows:

`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 (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Using LTS
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 kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}


Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.77775e-06, No Iterations 46
GAMGPCG:  Solving for pcorr, Initial residual = 0.552676, Final residual = 7.78633e-06, No Iterations 16
time step continuity errors : sum local = 0.17814, global = 9.87668e-16, cumulative = 9.87668e-16

Starting time loop

forces forces:
    Not including porosity effects
swak4Foam: Setting default mesh
forces forceCoeffs:
    Not including porosity effects
forceCoeffs forceCoeffs:
    Not including porosity effects
Flow time scale min/max = 0.010909, 0.177895
Smoothed flow time scale min/max = 0.010909, 0.107174
Time = 1

PIMPLE: iteration 1
smoothSolver:  Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 1
Phase-1 volume fraction = 0  Min(alpha.water) = 0  Max(alpha.water) = 0
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0  Min(alpha.water) = 0  Max(alpha.water) = 0
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00383769, No Iterations 13
GAMG:  Solving for p_rgh, Initial residual = 0.554143, Final residual = 0.00526051, No Iterations 3
time step continuity errors : sum local = 0.00829409, global = 4.89339e-17, cumulative = 1.0366e-15
GAMG:  Solving for p_rgh, Initial residual = 0.264695, Final residual = 0.00151991, No Iterations 11
GAMG:  Solving for p_rgh, Initial residual = 0.421563, Final residual = 9.97681e-08, No Iterations 35
time step continuity errors : sum local = 1.50186e-07, global = -1.52441e-17, cumulative = 1.02136e-15
smoothSolver:  Solving for omega, Initial residual = 0.078478, Final residual = 0.00444975, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0530609, No Iterations 3
ExecutionTime = 12.11 s  ClockTime = 12 s

forces forces write:
    sum of forces:
        pressure : (-14.4665 -0.883086 -8.65782)
        viscous  : (-0.0112302 0.000474369 0.000372435)
        porous   : (0 0 0)
    sum of moments:
        pressure : (1.08301 -0.825871 211.234)
        viscous  : (-6.58758e-05 -0.0027993 -0.00285493)
        porous   : (0 0 0)

Expression compute_ptot : [0] swak4Foam: Allocating new repository for sampledGlobalVariables
 min=-112.824 max=106.272 average=1.9263 weightedAverage=1.89152
Expression compute_ptot :  min=-112.824 max=106.272 average=1.9263 weightedAverage=1.89152
forceCoeffs forceCoeffs write:
    Cm    = 0.000164297
    Cd    = 0.000151967
    Cl    = -9.08732e-05
    Cl(f) = 0.000118861
    Cl(r) = -0.000209734`

I copied only the initial part of the log file; if you may need the entire part, please let me know and I will enclose here the gz file.
I let this case running twice, with the same identical result.
I have enough free space on my drive, so I might rule out that the simulation aborted for that.
How to deal with this?

Thank you for your help!

Berin

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions