|
[Sponsors] |
August 30, 2009, 03:45 |
Flow Around a Cylinder
|
#1 |
Member
toto
Join Date: Jun 2009
Posts: 71
Rep Power: 17 |
Hi,
I am trying to model a flow around a cylinder 2D using simpleFoam.The Velocity is 3m/s and The turbulence model is k-epsilon Standard. The dimension of the cylinder is 0.1 m diameter. How could i set the Boundary Condition for Wall- Top and Bottom to get no Slip? What about the Cylinder-Wall? Please it ist very important for me and I would greatly appreciate it if you could send me an Answer or a case like that. |
|
August 30, 2009, 17:10 |
|
#2 |
Senior Member
Steve Hansel
Join Date: Jun 2009
Location: Colorado, USA
Posts: 112
Rep Power: 17 |
In the boundaryField section of your first U file you would have something like:
minY { type fixedValue; value uniform (0 0 0); } where minY is what ever name you gave your sides and / or the cylinder-wall. See the tutorials for more examples. |
|
August 31, 2009, 04:51 |
|
#3 |
Member
toto
Join Date: Jun 2009
Posts: 71
Rep Power: 17 |
Thanks Hansel!
My Boundary Condition and so on look like: FoamFile { version 2.0; format ascii; class volVectorField; object U; } // // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { Rohrwand { type fixedValue; value uniform (0 0 0); } Outlet_Pres { type zeroGradient; } Inlet_Vel { type fixedValue; value uniform (2 0 0); } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { Rohrwand { type zeroGradient; ; } Outlet_Pres { type fixedValue; value uniform 0; } Inlet_Vel { type zeroGradient; } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object k; } // // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.87; boundaryField { Rohrwand { type zeroGradient; } Outlet_Pres { type zeroGradient; } Inlet_Vel { type fixedValue; value uniform 0.87; } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object epsilon; } // // dimensions [0 2 -3 0 0 0 0]; internalField uniform 8579.621; boundaryField { Rohrwand { type zeroGradient; } Outlet_Pres { type zeroGradient; } Inlet_Vel { type fixedValue; value uniform 8579.621; } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { Rohrwand { type nutWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } Outlet_Pres { type calculated; value uniform 0; } Inlet_Vel { type calculated; value uniform 0; } Sym_Wand { type calculated; value uniform 0; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class volSymmTensorField; object R; } // // dimensions [0 2 -2 0 0 0 0]; internalField uniform (0 0 0 0 0 0); boundaryField { Rohrwand { type zeroGradient; } Outlet_Pres { type zeroGradient; } Inlet_Vel { type fixedValue; value uniform (0 0 0 0 0 0 0); } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class polyBoundaryMesh; location "constant/polyMesh"; object boundary; } // // 5( Rohrwand { type wall; physicalType wallFunctions; nFaces 100; startFace 182530; } Outlet_Pres { type patch; physicalType outlet; nFaces 300; startFace 182630; } Inlet_Vel { type patch; physicalType inlet; nFaces 100; startFace 182930; } Sym_Wand { type symmetryPlane; nFaces 700; startFace 183030; } frontAndBackPlanes { type empty; nFaces 183130; startFace 183730; } ) // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object RASProperties; } // // RASModel kEpsilon; turbulence on; printCoeffs on; FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object transportProperties; } // // transportModel Newtonian; nu nu [0 2 -1 0 0 0 0] 1.46e-05; CrossPowerLawCoeffs { nu0 nu0 [0 2 -1 0 0 0 0] 0; nuInf nuInf [0 2 -1 0 0 0 0] 0; m <> [0 0 0 0 0 0 0] 0; n <> [0 0 0 0 0 0 0] 0; } BirdCarreauCoeffs { nu0 nu0 [0 2 -1 0 0 0 0] 0; nuInf nuInf [0 2 -1 0 0 0 0] 0; k <> [0 0 0 0 0 0 0] 0; n <> [0 0 0 0 0 0 0] 0; } FoamFile { version 2.0; format ascii; class volScalarField; object nuTilda; } // * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { Rohrwand { type zeroGradient; } Outlet_Pres { type zeroGradient; } Inlet_Vel { type fixedValue; value uniform 0; } Sym_Wand { type symmetry; } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // turbulenceModel kEpsilon; turbulence on; laminarCoeffs {} kEpsilonCoeffs { Cmu Cmu [0 0 0 0 0 0 0] 0.09; C1 C1 [0 0 0 0 0 0 0] 1.44; C2 C2 [0 0 0 0 0 0 0] 1.92; alphaEps alphaEps [0 0 0 0 0 0 0] 0.76923; } I get this Error when i start the Simulation: Inconsistent Patch and patchField types for patch type symmetryPlane and patchField type symmety When i change SymmetryPlane to slip or Wall , my solution is absurdity. I need Help please! |
|
August 31, 2009, 10:29 |
symmetry should be symmetryPlane
|
#4 |
Senior Member
Richard Smith
Join Date: Mar 2009
Location: Enfield, NH, USA
Posts: 138
Blog Entries: 4
Rep Power: 17 |
Judging by your error message I think you need to replace:
Sym_Wand { type symmetry; } with: Sym_Wand { type symmetryPlane; } within your case/0/* field files.
__________________
Symscape, Computational Fluid Dynamics for all |
|
September 14, 2009, 12:25 |
|
#5 |
Senior Member
Steve Hansel
Join Date: Jun 2009
Location: Colorado, USA
Posts: 112
Rep Power: 17 |
I think the units on k are supposed to be [0 0 1 0 0 0 0]
and m [0 0 1 0 0 0 0] |
|
September 18, 2009, 09:13 |
|
#6 |
Member
toto
Join Date: Jun 2009
Posts: 71
Rep Power: 17 |
It is the units on omega! The units on k is [0 2 -2 0 0 0 0 ]
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
benchmark: flow over a circular cylinder | goodegg | Main CFD Forum | 12 | January 22, 2013 12:47 |
Flow over a flat plate & Flow over a cylinder | cfdxue | Main CFD Forum | 0 | November 27, 2007 00:26 |
Anybody knows smth abt the flow around a cylinder? | Wenxuan | Main CFD Forum | 3 | March 20, 2007 17:19 |
Flow over a cylinder | Anna | Main CFD Forum | 9 | March 24, 2006 15:32 |
Flow Past a Cylinder | Annie | Main CFD Forum | 2 | March 31, 2002 01:02 |