[Virtual-iCFD]

Hi foamer,

I experienced the pisoFOAM crashes and am wondering any workaround on this issue.

Thanks,

Cheers,
Pete


reate time

Create mesh for time = 0


PISO: Operating solver in PISO mode

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No MRF models present

No finite volume options present


Starting time loop

Time = 1e-05

Courant Number mean: 0 max: 0
smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 9.62803e-06, No Iterations 66
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0317079, No Iterations 1
time step continuity errors : sum local = 0.0574054, global = -0.0364258, cumulative = -0.0364258
GAMG: Solving for p, Initial residual = 0.0438414, Final residual = 9.43287e-07, No Iterations 43
time step continuity errors : sum local = 3.30835e-06, global = -1.73591e-07, cumulative = -0.036426
ExecutionTime = 47.94 s ClockTime = 48 s

Time = 2e-05

Courant Number mean: 0.490971 max: 1208.82
[17] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[17] #1 Foam::sigFpe::sigHandler(int) at ??:?
[17] #2 ? in "/lib64/libc.so.6"
[17] #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[17] #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[17] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[17] #6 ? at ??:?
[17] #7 ? at ??:?
[17] #8 ? at ??:?
[17] #9 ? at ??:?
[17] #10 __libc_start_main in "/lib64/libc.so.6"
[17] #11 ? at ??:?

[anon_q]

Your courant number is very high it shouldn't be greater than 1 for PISO algorithm.
Just decrease your time step.

You cannot use adjustTimeStep and maxCo with pisoFoam unless you modify the original files to add this features.

In other hand, pimpleFoam can be run essentially as pisoFoam by setting nOuterCorrectors to 1. you can add this to system/controlDict:

adjustTimeStep    on;
maxCo            0.5;

[RobertHB]

I'd guess something is wrong with your geometry.