|
[Sponsors] |
August 6, 2014, 14:34 |
simpleFoam+T NaN`s problem
|
#1 |
Member
Join Date: May 2013
Posts: 51
Rep Power: 12 |
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 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 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 seav Last edited by seav; August 7, 2014 at 08:03. |
|
August 7, 2014, 03:26 |
|
#2 |
Member
Join Date: May 2013
Posts: 51
Rep Power: 12 |
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. |
|
August 10, 2014, 04:05 |
|
#3 |
Member
Join Date: May 2013
Posts: 51
Rep Power: 12 |
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 |
|
|
|
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 |