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

Problem with geometry - concentric cylinders

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 15, 2014, 10:22
Default Problem with geometry - concentric cylinders
  #1
New Member
 
Rhoddwen
Join Date: Dec 2014
Posts: 9
Rep Power: 11
Rhoddwen is on a distinguished road
Hello everyone,

I am trying to simulate the flow between two concentric cylinders. The inner one rotating (R = 25), the outer one fixed (R=30). Moreover, the inner cylinder should be connected at the upper wall which rotates with-it. My problem is that I have tried everything but I cannot build my geometry....I am not using anything else but open foam.I was wondering if enyone can help me figure this out.The dimensions of inner and outer cylinders have no relevance..the geometry shape is important.

So..my case is a multiphase case, I use interFoam solver.

My geometry is build-ed as axi-simmetric wedge. I don't know if it is correct, cause I cannot manage to obtain the ''in-between of cylinders'' flow domain. It builds every time the block relying the inner cylinder to the axis (the cone - see also attached piks).
I tried everything:
- defining all blocks and set "empty" boundaries on the 'hex' that normally does not belong to the flow domain; but than open foam build it anyway, and fills it with air;
- not defing the 'cone' relying my flow domain to the axis; but then 'blockMesh' and 'check mesh' failed...as wedge is not recognized, even if I used a small angle;

I am sending you all the data...the log, the mesh file, set fields, alpha.org...

Please answer me these questions..if you know.. :
1. How do I build a geometry consisting within 2 concentric cylinders in openF.

2. How do I change the coordinate system from Cartesian to Cylindrical;

3. How do I set a precise volume of "water" within my flow domain - the gap between 2 concentric cylinders.

I hope that anyone can help me...
Thank you.
---------------------------------------------------------------------------------------------------------------------------
LOG FILE
---------------------------------------------------------------------------------------------------------------------------
BLOCK MESH LOG
-----------------------------------------------------------------------------------------

Creating block mesh from
"/home/FOAM_RUN/DC/interFoam/laminar/pocheGR/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

Basic statistics
Number of internal faces : 5
Number of boundary faces : 20
Number of defined boundary faces : 20
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1
Block 0 cell size :
i : 0.0333333 .. 0.0333333
j : 0.416667 .. 0.416667
k : 0 .. 0

Block 1 cell size :
i : 0.0166667 .. 0.0166667
j : 0.416667 .. 0.416667
k : 0 .. 0

Block 2 cell size :
i : 0.0333333 .. 0.0333333
j : 0.0333333 .. 0.0333333
k : 2.49896 .. 2.49896

Block 3 cell size :
i : 0.0166667 .. 0.0166667
j : 0.0333333 .. 0.0333333
k : 2.49896 .. 2.49896

Block 4 cell size :
i : 0.0333333 .. 0.0333333
j : 0.0666667 .. 0.0666667
k : 2.69888 .. 2.69888


Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0 0 -1.54935) (3 30.9613 1.54935)
nPoints: 36481
nCells: 18000
nFaces: 72180
nInternalFaces: 35700
----------------
Patches
----------------
patch 0 (start: 35700 size: 120) name: leftWall
patch 1 (start: 35820 size: 120) name: rightWall
patch 2 (start: 35940 size: 120) name: sides
patch 3 (start: 36060 size: 0) name: axis
patch 4 (start: 36060 size: 60) name: atmosphere
patch 5 (start: 36120 size: 60) name: upperWall
patch 6 (start: 36180 size: 18000) name: back
patch 7 (start: 54180 size: 18000) name: front

End

-----------------------------------------------------------------------------------------
CHECK MESH LOG
-----------------------------------------------------------------------------------------


Create polyMesh for time = 0

Time = 0

Mesh stats
points: 36481
internal points: 0
faces: 72180
internal faces: 35700
cells: 18000
faces per cell: 5.99333
boundary patches: 8
point zones: 0
face zones: 0
cell zones: 0

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

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
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
leftWall 120 242 ok (non-closed singly connected)
rightWall 120 244 ok (non-closed singly connected)
sides 120 242 ok (non-closed singly connected)
axis 0 0 ok (empty)
atmosphere 60 122 ok (non-closed singly connected)
upperWall 60 122 ok (non-closed singly connected)
back 18000 18301 ok (non-closed singly connected)
front 18000 18301 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 0 -1.54935) (3 30.9613 1.54935)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (0 1 1)
Wedge back with angle 2.86479 degrees
Wedge front with angle 2.86479 degrees
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (1.66689e-17 -5.15476e-15 5.04149e-15) OK.
Max cell openness = 1.65719e-16 OK.
Max aspect ratio = 24.7607 OK.
Minimum face area = 0.000555555. Maximum face area = 1.03127. Face area magnitudes OK.
Min volume = 0.000144435. Max volume = 0.0343755. Total volume = 132.329. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.330003 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

-----------------------------------------------------------------------------------------
SET FIELDS LOG
-----------------------------------------------------------------------------------------

Create time

Create mesh for time = 0

Reading setFieldsDict

Setting field default values
Setting internal values of volScalarField alpha.water

Setting field region values
Adding cells with centre within cylinder, with p1 = (0.2 25 -1), p2 = (0.2 26 1) and radius = 26
Setting internal values of volScalarField alpha.water

End

-------------------------------------------------------------------------------------------------------------------------
MESH FILE
----------------------------------------------------------------------------------------------------------------------
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;
//puncte aranjate pe verticala//
vertices
(
(0 0 0)//p0//
(2 0 0)//p1//
(3 0 0)//p2//
(0 24.96875650987 -1.24947923177)//p3//
(2 24.96875650987 -1.24947923177)//p4//
(3 24.96875650987 -1.24947923177)//p5//
(3 26.96625703066 -1.34943757031)//p6//
(2 26.96625703066 -1.34943757031)//p7//
(0 26.96625703066 -1.34943757031)//p8//
(0 30.96125807224 -1.54935424739)//p9//
(2 30.96125807224 -1.54935424739)//p10//
(0 24.96875650987 1.24947923177)//p11//
(2 24.96875650987 1.24947923177)//p12//
(3 24.96875650987 1.24947923177)//p13//
(3 26.96625703066 1.34943757031)//p14//
(2 26.96625703066 1.34943757031)//p15//
(0 26.96625703066 1.34943757031)//p16//
(0 30.96125807224 1.54935424739)//p17//
(2 30.96125807224 1.54935424739)//p18//
);

blocks
(
hex (0 1 4 3 0 1 12 11) (60 60 1) simpleGrading (1 1 1)
hex (1 2 5 4 1 2 13 12) (60 60 1) simpleGrading (1 1 1)
hex (3 4 7 8 11 12 15 16) (60 60 1) simpleGrading (1 1 1)
hex (4 5 6 7 12 13 14 15) (60 60 1) simpleGrading (1 1 1)
hex (8 7 10 9 16 15 18 17) (60 60 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
leftWall
{
type walls;
faces
(
(3 8 16 11)
(8 9 17 16)
);
}
rightWall
{
type wall;
faces
(
(5 6 14 13)
(7 10 18 15)
);
}
sides
{
type empty;
faces
(
(0 3 11 0)
(2 5 13 2)
);
}
axis
{
type empty;
faces
(
(0 1 1 0)
(1 2 2 1)
);
}
atmosphere
{
type wall;//patch;
faces
(
(9 10 18 17)
);
}
upperWall
{
type wall;
faces
(
(7 6 14 15) // (10 9 17 18)//
);
}
back
{
type wedge;
faces
(
(0 1 4 3)
(1 2 5 4)
(3 4 7 8)
(4 5 6 7)
(8 7 10 9)
);
}
front
{
type wedge;
faces
(
(0 1 12 11)
(1 2 13 12)
(11 12 15 16)
(12 13 14 15)
(16 15 18 17)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //

--------------------------------------------------------------------------------------------------------------------------
BOUNDARY CONDITIONS FILE
-------------------------------------------------------------------------------------------------------------------------
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
leftWall
{
type walls;
nFaces 120;
startFace 35700;
faces ( ( 3 8 16 11 ) ( 8 9 17 16 ) );
}
rightWall
{
type wall;
inGroups 1(wall);
nFaces 120;
startFace 35820;
}
sides
{
type empty;
inGroups 1(empty);
nFaces 120;
startFace 35940;
}
axis
{
type empty;
inGroups 1(empty);
nFaces 0;
startFace 36060;
}
atmosphere
{
type wall;
inGroups 1(wall);
nFaces 60;
startFace 36060;
}
upperWall
{
type wall;
inGroups 1(wall);
nFaces 60;
startFace 36120;
}
back
{
type wedge;
inGroups 1(wedge);
nFaces 18000;
startFace 36180;
}
front
{
type wedge;
inGroups 1(wedge);
nFaces 18000;
startFace 54180;
}
)
-----------------------------------------------------------------------------------------------------------------------
SET FIELDS FILE
------------------------------------------------------------------------------------------------------------------------
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object setFieldsDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

defaultFieldValues
(
volScalarFieldValue alpha.water 0
);

regions
(
cylinderToCell //boxToCell //box (3 0 -1.34947923177) (3 26.96625703066 1.34943757031);
{
p1 (0.2 25 -1); //Min
p2 (0.2 26 1); //Max
radius 26;
// radius 26;

fieldValues
(
volScalarFieldValue alpha.water 1
);
}
);

// ************************************************** *********************** //

-----------------------------------------------------------------------------------------------------------------------
ALPHA WATER FILE

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
leftWall
{
type zeroGradient;
}

rightWall
{
type zeroGradient;
}

axis
{
type empty;
}

upperWall
{
type zeroGradient;
}

front
{
type wedge;
}

back
{
type wedge;
}

atmosphere
{
type zeroGradient;//type inletOutlet;
//inletValue uniform 0;
//value uniform 0;
}

defaultFaces
{
type empty;
}
}

// ************************************************** *********************** //
Attached Files
File Type: gz Piks.tar.gz (80.5 KB, 1 views)
Rhoddwen is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] Volume Meshing Problem; YOUR GEOMETRY HAS A HOLE mariachi ANSYS Meshing & Geometry 5 May 30, 2018 03:52
Problem with building geometry - concentric cylinders Rhoddwen OpenFOAM Pre-Processing 9 December 17, 2014 11:02
problem in importing geometry R4RAHUL FLUENT 0 July 4, 2014 12:53
problem when imported geometry from 3D CAD to star ccm, TAREK GANAT STAR-CCM+ 1 May 21, 2013 23:15
vitual _ real deneb FLUENT 3 January 22, 2007 05:31


All times are GMT -4. The time now is 05:15.