
[Sponsors] 
simpleFoam for complicated geometry, nonorthogonal mesh, omega wall function 

LinkBack  Thread Tools  Search this Thread  Display Modes 
July 27, 2012, 06:00 
simpleFoam for complicated geometry, nonorthogonal mesh, omega wall function

#1 
Member
achinta
Join Date: May 2010
Location: Sydney
Posts: 66
Rep Power: 13 
Hello everyone,
I am using simpleFoam on complicated geometry (non orthogonal mesh). Here is the checkMesh report:  Mesh stats points: 6487723 faces: 19200050 internal faces: 18938278 cells: 6356388 boundary patches: 16 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 6356388 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). . . Mesh (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (9.64958e17 2.34786e15 6.13025e16) OK. Max cell openness = 4.58304e15 OK. Max aspect ratio = 299.485 OK. Minumum face area = 3.41806e11. Maximum face area = 0.000249108. Face area magnitudes OK. Min volume = 1.22828e14. Max volume = 4.16394e07. Total volume = 0.20005. Cell volumes OK. Mesh nonorthogonality Max: 87.0638 average: 19.2898 *Number of severely nonorthogonal faces: 18931. Nonorthogonality check OK. <<Writing 18931 nonorthogonal faces to set nonOrthoFaces Face pyramids OK. Max skewness = 2.88371 OK. Coupled point location match (average 0) OK. Mesh OK.  I am using SST turbulence model. I am using 'flowRateInletVelocity' for few inlets and it's mentioned like this: INLET { type flowRateInletVelocity; flowRate constant XXXX; // volume flow rate value uniform (0 0 0); } The boundary conditions for k and omega are:  Feature k, omega Turbulent inlet 0.4, 4000 Other inlets 1e6, 1e2 Walls type kqRWallFunction; type omegaWallFunction; value uniform 0.4; value uniform 4000; Outlet zeroGradient, zeroGradient internal field 0.4, 4000  I am using DICPCG for pressure, smoothSolver for omega and DILUPBiCG for others. My fvSchemes were as follows:  ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwindV grad(U); div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; p; }  The simulation diverges after 150 steps. I tried Gauss linear limited 0.333 for laplacian schemes which didn't work. gradSchemes was changed to 'cellLimited Gauss linear 1' and it was not good either. Could someone give me tips to improve the simulation? Note: When i use epsilonWallFunction(as wall boundary condition) for omega, it converges to some decent results (not accurate though) with the 1st set of fvSchemes. It's strange!! Kind regards, Achinta Last edited by achinta; July 27, 2012 at 06:07. Reason: typos 

October 20, 2013, 09:22 
how did u determine k and omega for complicated geometry

#2 
Senior Member
Hasan K.J.
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 200
Rep Power: 12 
Hey achinta,
I am stuck with similar problem but I'm a one step before u How did u determine the K and Omega for the complicated geometry, and did u over come the problem ? i am having the same problem Regards, Hasan. 

October 25, 2013, 00:58 
k and omega

#3  
Member
Dr. B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 55
Rep Power: 11 
Quote:
Try to calculate k and epsilon from the intensity and velocity, which you know for inlet. If not, get it from experiments(available in literature). Then calculate omega by omega=epsilon/k  KANNAN 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
3D Hybrid Mesh Errors  DarrenC  ANSYS Meshing & Geometry  11  August 5, 2013 06:42 
[blockMesh] error message with modeling a cube with a hold at the center  hsingtzu  OpenFOAM Meshing & Mesh Conversion  2  March 14, 2012 09:56 
ParaView for OF1.6ext  Chrisi1984  OpenFOAM Installation  0  December 31, 2010 06:42 
due to wall function, continue refining mesh until when  zandi  CFX  8  January 23, 2010 04:16 
Error with Wmake  skabilan  OpenFOAM Installation  3  July 28, 2009 00:35 