[samiam1000]

Dear All,

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

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]

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

[samiam1000]

That's good: thanks a lot for help.

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

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]

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]

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]

Any news? Any idea?

[mikeP]

any news on this topic?

have you solved your problem samiam?

[samiam1000]

I did.

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

This works well!