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

simpleFoam for complicated geometry, non-orthogonal mesh, omega wall function

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 27, 2012, 06:00
Default simpleFoam for complicated geometry, non-orthogonal mesh, omega wall function
  #1
Member
 
achinta
Join Date: May 2010
Location: Sydney
Posts: 66
Rep Power: 7
achinta is on a distinguished road
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 (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (9.64958e-17 2.34786e-15 -6.13025e-16) OK.
Max cell openness = 4.58304e-15 OK.
Max aspect ratio = 299.485 OK.
Minumum face area = 3.41806e-11. Maximum face area = 0.000249108. Face area magnitudes OK.
Min volume = 1.22828e-14. Max volume = 4.16394e-07. Total volume = 0.20005. Cell volumes OK.
Mesh non-orthogonality Max: 87.0638 average: 19.2898
*Number of severely non-orthogonal faces: 18931.
Non-orthogonality check OK.
<<Writing 18931 non-orthogonal 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 1e-6, 1e-2

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
achinta is offline   Reply With Quote

Old   October 20, 2013, 09:22
Default how did u determine k and omega for complicated geometry
  #2
Senior Member
 
Alhasan's Avatar
 
Hasan K.J.
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 181
Rep Power: 5
Alhasan is on a distinguished road
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.
Alhasan is offline   Reply With Quote

Old   October 25, 2013, 00:58
Default k and omega
  #3
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 6
skyinventorbt is on a distinguished road
Quote:
Originally Posted by Alhasan View Post
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.
Dear Hasan,
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
skyinventorbt 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
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42
error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Native Meshers: blockMesh 2 March 14, 2012 10:56
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
due to wall function, continue refining mesh until when zandi CFX 8 January 23, 2010 05:16
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 00:35


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