CFD Online Discussion Forums - simpleFoam: strange error

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - simpleFoam: strange error (https://www.cfd-online.com/Forums/openfoam/107319-simplefoam-strange-error.html)

simpleFoam: strange error

Dear All,

when I try to launch my case, using the simpleFoam solve, I get this error:

Code:

pc@pc:/OPENFOAM/cases/moving_door_def_flow/test/steadyNoBordiNoEnergy$ rm 0/* ; cp 0/originalFiles/* 0/ && simpleFoam 

rm: cannot remove `0/originalFiles': Is a directory

/*---------------------------------------------------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \\    /   O peration     | Version:  2.1.0                                 |

|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

Build  : 2.1.0-0bc225064152

Exec   : simpleFoam

Date   : Sep 24 2012

Time   : 09:54:39

Host   : "slnxepmi05"

PID    : 2344

Case   : /OPENFOAM/cases/moving_door_def_flow/test/steadyNoBordiNoEnergy

nProcs : 1

sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

fileModificationChecking : Monitoring run-time modified files using timeStampMaster

allowSystemOperations : Disallowing user-supplied system call operations



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time



Create mesh for time = 0



Reading field p



Reading field U



Reading/calculating face flux field phi



Selecting incompressible transport model Newtonian

Selecting RAS turbulence model realizableKE

--> Upgrading k to employ run-time selectable wall functions

    Backup original k to k.old

    Writing updated k

--> Upgrading epsilon to employ run-time selectable wall functions

    Backup original epsilon to epsilon.old

    Writing updated epsilon

--> Creating nut to employ run-time selectable wall functions

    Writing new nut

realizableKECoeffs

{

    Cmu             0.09;

    A0              4;

    C2              1.9;

    sigmak          1;

    sigmaEps        1.2;

}



No field sources present





SIMPLE: convergence criteria

    field p_rgh         tolerance 1e-08

    field U         tolerance 1e-08

    field h         tolerance 1e-08

    field "(k|epsilon|omega)"         tolerance 1e-08





Starting time loop



faceSource faceObj1:

    total faces  = 188

    total area   = 0.024





Time = 1







--> FOAM FATAL IO ERROR: 

attempt to read beyond EOF



file: /OPENFOAM/cases/moving_door_def_flow/test/steadyNoBordiNoEnergy/system/fvSchemes::divSchemes::div((nuEff*dev(T(grad(U))))) at line 37.



    From function ITstream::read(token&)

    in file db/IOstreams/Tstreams/ITstream.C at line 83.



FOAM exiting



pc@pc:/OPENFOAM/cases/moving_door_def_flow/test/steadyNoBordiNoEnergy$

How can I solve it?
Thanks a lot,

Samuele

it seems you may lose a semicolon or a bracket in your fvScheme, its a typo error

That's good: thanks a lot for help.

Now it runs, but I do have another (bigger?!?!) problem: I get this error message

Code:

Time = 729



smoothSolver:  Solving for Ux, Initial residual = 0.34966013, Final residual = 0.014146042, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.42149661, Final residual = 0.018716496, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.36537674, Final residual = 0.014728068, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.47155456, Final residual = 0.022915631, No Iterations 2

time step continuity errors : sum local = 1.2745821e+92, global = 5.0729639e+90, cumulative = 1.7462592e+91

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#2   in "/lib/x86_64-linux-gnu/libc.so.6"

#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) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"

#7  Foam::fvMatrix<double>::solve() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"

#8  Foam::incompressible::RASModels::realizableKE::correct() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"

#9  

 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/simpleFoam"

#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#11  

 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/simpleFoam"

Floating point exception

pc@pc:/OPENFOAM/cases/moving_door_def_flow/test/steadyNoBordiNoEnergy$ rm -r 1* 2* 3* 4* 5* 6* 7* 8* 9* faceObj1/ log ; rm 0/* ; cp 0/originalFiles/* 0/ && simpleFoam

The simulation does not converge. The point is that if I run the same case with buoyancy (bouyantSimpleFoam instead of simpleFoam) it works.

What can I do?

Thanks a lot,
Samuele

Quote:

time step continuity errors : sum local = 1.2745821e+92, global = 5.0729639e+90, cumulative = 1.7462592e+91

this problem can be caused by several reasons,
maybe you need to increase iteration or change BC or refine your mesh
however it seems your numerical calculation is diverging

As far as BC are concerned, what do you suggest about pressure? I have 2 inlets and 2 outlets. I need to have the same mass flow entering from inlet1 goes out from outlet1 and the same with inlet2 and outlet2.

Also, the mesh should be fine, since it does not diverge if I use the buoyantSimpleFoam solver.

Also, where should I increase the numeber of iterations? I can't find the right file to edit.

Thanks a lot,
Samuele

any news on this topic?

have you solved your problem samiam?

I did.

I fixed the massflow at inlet1, inlet2 and outlet1 and I set zeroGradient as BC on velocity at outlet2.

This works well!