CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   simpleFoam: strange error (http://www.cfd-online.com/Forums/openfoam/107319-simplefoam-strange-error.html)

samiam1000 September 24, 2012 04:14

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

nimasam September 24, 2012 11:20

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

samiam1000 September 25, 2012 06:05

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

nimasam September 25, 2012 07:13

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

samiam1000 September 25, 2012 08:30

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

samiam1000 September 27, 2012 06:54

Any news? Any idea?

mikeP December 11, 2012 04:47

any news on this topic?

have you solved your problem samiam?

samiam1000 December 11, 2012 04:52

I did.

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

This works well!


All times are GMT -4. The time now is 07:39.