CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

simpleFoam+T NaN`s problem

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By seav
  • 2 Post By seav

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 6, 2014, 14:34
Default simpleFoam+T NaN`s problem
  #1
Member
 
Join Date: May 2013
Posts: 51
Rep Power: 12
seav is on a distinguished road
Hello I got this kind of error :
Code:
--> FOAM FATAL IO ERROR: 
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/jils/Dokumenty/mgr/mgr_projekt_dol/system/data.solverPerformance.p at line 0.

    From function operator>>(Istream&, Scalar&)
    in file lnInclude/Scalar.C at line 93.

FOAM exiting
And the question is how to solve it ? I had similiar problem with epsilon and change of BC helped, but here I am out of ideas.

I model flow with v=0.5m/s in rectangular channel with obstacles inside. Channel is 140x450x5 mm wide/long/high. I made control volumes before inlet and behind outlet to prevent compute errors.. And here errors started.

Check mesh is fine :
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : checkMesh
Date   : Aug 06 2014
Time   : 20:19:25
Host   : "ubuntuspace"
PID    : 8754
Case   : /home/jils/Dokumenty/mgr/mgr_projekt_dol
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 polyMesh for time = 0

Time = 0

Mesh stats
    points:           4154652
    faces:            11667836
    internal faces:   10875874
    cells:            3757285
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       140

Overall number of cells of each type:
    hexahedra:     3757285
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    inlet_dol           2870     3168     ok (non-closed singly connected)  
    outlet_dol          2870     3168     ok (non-closed singly connected)  
    boki_przednie       9270     10208    ok (non-closed singly connected)  
    boki_tylne          9270     10208    ok (non-closed singly connected)  
    defaultFaces        767682   770148   ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-96.5 0.5 0) (546.5 5.5 139.998)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-2.62359e-17 -1.90826e-15 -2.4542e-17) OK.
    Max cell openness = 3.20902e-16 OK.
    Max aspect ratio = 3.06432 OK.
    Minimum face area = 0.0760913. Maximum face area = 0.562604.  Face area magnitudes OK.
    Min volume = 0.0217404. Max volume = 0.216776.  Total volume = 449644.  Cell volumes OK.
    Mesh non-orthogonality Max: 55.498 average: 15.1534
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.76418 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
I attach my 0 and system folders.

Code:
Time = 25

DILUPBiCG:  Solving for Ux, Initial residual = 0.699179, Final residual = 5.2484e-08, No Iterations 21
DILUPBiCG:  Solving for Uy, Initial residual = 0.699603, Final residual = 8.64538e-08, No Iterations 20
DILUPBiCG:  Solving for Uz, Initial residual = 0.699334, Final residual = 5.93537e-08, No Iterations 20
DICPCG:  Solving for p, Initial residual = 0.686836, Final residual = nan, No Iterations 1001
time step continuity errors : sum local = nan, global = -nan, cumulative = -nan
DILUPBiCG:  Solving for T, Initial residual = nan, Final residual = nan, No Iterations 1001
smoothSolver:  Solving for epsilon, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = nan, Final residual = nan, No Iterations 1000
ExecutionTime = 9236.17 s  ClockTime = 9255 s
Help me :)

seav
Attached Files
File Type: gz system.tar.gz (1.2 KB, 5 views)
File Type: gz 0.tar.gz (775 Bytes, 14 views)
Ramzy1990 likes this.

Last edited by seav; August 7, 2014 at 08:03.
seav is offline   Reply With Quote

Old   August 7, 2014, 03:26
Default
  #2
Member
 
Join Date: May 2013
Posts: 51
Rep Power: 12
seav is on a distinguished road
Found something interesting here : http://www.cfd-online.com/Forums/ope...tml#post202437

---------
Edit :

I think there is problem with BC on pressure can anyone give an example of that BC ?

Last edited by seav; August 7, 2014 at 05:01.
seav is offline   Reply With Quote

Old   August 10, 2014, 04:05
Default
  #3
Member
 
Join Date: May 2013
Posts: 51
Rep Power: 12
seav is on a distinguished road
Ok. I have found out what was bad.

Boundary Conditions for epsilon and k were set incorrectly. I have used those patterns to calculated them properly.

seav
wyldckat and Ramzy1990 like this.
seav is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
simpleFoam - Convergence problem - Simple rectangular prism domain HakikiCanakkaleli OpenFOAM Running, Solving & CFD 5 May 9, 2013 04:59
Problem in running simpleFoam tutorial Lana OpenFOAM Running, Solving & CFD 1 April 10, 2013 11:08
Covergence Problem with simpleFoam sfigato OpenFOAM Running, Solving & CFD 1 February 6, 2013 07:31
2D foil with simpleFoam: stability problem Nolwenn OpenFOAM Running, Solving & CFD 0 July 5, 2012 05:16
SimpleFOAM + SST-Model + problem with convergence A.Devesa OpenFOAM Running, Solving & CFD 0 November 9, 2010 04:43


All times are GMT -4. The time now is 02:50.