CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

simpleFoam: strange error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   September 24, 2012, 04:14
Default simpleFoam: strange error
  #1
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 453
Rep Power: 8
samiam1000 is on a distinguished road
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
samiam1000 is offline   Reply With Quote

Old   September 24, 2012, 11:20
Default
  #2
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,063
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
it seems you may lose a semicolon or a bracket in your fvScheme, its a typo error
nimasam is offline   Reply With Quote

Old   September 25, 2012, 06:05
Default
  #3
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 453
Rep Power: 8
samiam1000 is on a distinguished road
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
samiam1000 is offline   Reply With Quote

Old   September 25, 2012, 07:13
Default
  #4
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,063
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
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
nimasam is offline   Reply With Quote

Old   September 25, 2012, 08:30
Default
  #5
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 453
Rep Power: 8
samiam1000 is on a distinguished road
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 is offline   Reply With Quote

Old   September 27, 2012, 06:54
Default
  #6
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 453
Rep Power: 8
samiam1000 is on a distinguished road
Any news? Any idea?
samiam1000 is offline   Reply With Quote

Old   December 11, 2012, 03:47
Default
  #7
Member
 
Join Date: Jun 2011
Posts: 35
Rep Power: 5
mikeP is on a distinguished road
any news on this topic?

have you solved your problem samiam?
mikeP is offline   Reply With Quote

Old   December 11, 2012, 03:52
Default
  #8
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 453
Rep Power: 8
samiam1000 is on a distinguished road
I did.

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

This works well!
samiam1000 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
polynomial thermophysical properties II sebastian OpenFOAM Running, Solving & CFD 53 April 10, 2014 05:45
Compile calcMassFlowC aurore OpenFOAM Programming & Development 12 March 18, 2014 05:22
Installation of Netgen in SuSE Linux 92 edvardsenpriv Open Source Meshers: Gmsh, Netgen, CGNS, ... 23 January 16, 2009 06:12
How to get the max value of the whole field waynezw0618 OpenFOAM Running, Solving & CFD 4 June 17, 2008 05:07
user defined function cfduser CFX 0 April 29, 2006 10:58


All times are GMT -4. The time now is 01:38.