# [blockMesh] Cavity with Internal Solid Cylinder

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 21, 2018, 16:15 Cavity with Internal Solid Cylinder #1 New Member   Join Date: Feb 2018 Posts: 5 Rep Power: 5 I'm modelling a water filled cavity with a heated central cylinder. I've created geometry as shown below I've done this by creating the outer cavity, and then an inner cylinder block. Code: ```/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 1; vertices ( (0 0 0) (1 0 0) (1 0.5 0) (0 0.5 0) (0 0 1) (1 0 1) (1 .5 1) (0 0.5 1) (0.40454 0.022 0.40454) (0.59546 0.022 0.40454) (0.59546 0.478 0.40454) (0.40454 0.478 0.40454) (0.40454 0.022 0.59546) (0.59546 0.022 0.59546) (0.59546 0.478 0.59546) (0.40454 0.478 0.59546) ); geometry { }; blocks ( hex (0 1 2 3 4 5 6 7) (50 50 50) simpleGrading ( 1 1 1 ) hex (8 9 10 11 12 13 14 15) (10 10 10) simpleGrading ( 1 1 1 ) ); edges ( arc 8 9 (0.5 0.022 0.365) arc 10 11 (0.5 0.478 0.365) arc 9 13 (0.635 0.022 0.5) arc 10 14 (0.635 0.478 0.5) arc 13 12 (0.5 0.022 0.635) arc 14 15 (0.5 0.478 0.635) arc 12 8 (0.365 0.022 0.5) arc 15 11 (0.365 0.478 0.5) ); boundary ( defaultFaces { type wall; faces ( (9 13 14 10) (13 14 15 12) (12 8 11 15) (8 9 10 11) (8 9 13 12) (11 10 14 15) ); } fixedWalls { type wall; faces ( (0 3 2 1) (4 5 6 7) ); } leftWall { type wall; faces ( (0 4 7 3) ); } rightWall { type wall; faces ( (1 2 6 5) ); } front { type cyclic; neighbourPatch back; faces ( (0 1 5 4) ); } back { type cyclic; neighbourPatch front; faces ( (2 3 7 6) ); } ); mergePatchPairs ( ); // ************************************************************************* //``` I've then specified the cylinder wall temperature as uniform. However, only the interior of the cylinder heats, and not the surrounding cavity. Code: ```/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 293.15; boundaryField { defaultFaces { type fixedValue; value uniform 353.15; } leftWall { type zeroGradient; } rightWall { type zeroGradient; } fixedWalls { type zeroGradient; } front { type zeroGradient; } back { type zeroGradient; } } // ************************************************************************* //``` How do I get the exterior of the cylinder to interact, or is there a better way of going about this?

 March 22, 2018, 00:34 #2 Senior Member     Uwe Pilz Join Date: Feb 2017 Location: Leipzig, Germany Posts: 588 Rep Power: 10 You need to define the boundary between the cylinder and the surroundings as wall or patch and give boundary conditions. __________________ Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)

March 22, 2018, 04:17
#3
New Member

Join Date: Feb 2018
Posts: 5
Rep Power: 5
Quote:
 Originally Posted by piu58 You need to define the boundary between the cylinder and the surroundings as wall or patch and give boundary conditions.
Is that not done? (Poor naming but defaultFaces=the cylinder)

Code:
```  defaultFaces
{
type wall;
faces
(
(9 13 14 10)
(13 14 15 12)
(12 8 11 15)
(8 9 10 11)
(8 9 13 12)
(11 10 14 15)
);
}```
Code:
```boundaryField
{
defaultFaces
{
type            fixedValue;
value           uniform 353.15;
}```

 March 22, 2018, 05:16 #4 Senior Member     Uwe Pilz Join Date: Feb 2017 Location: Leipzig, Germany Posts: 588 Rep Power: 10 I don't recommend to use defaultFaces. Give the boundary a name and set the condition. Then you know for sure everything is ok. __________________ Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)

 March 22, 2018, 08:34 #5 New Member   Join Date: Feb 2018 Posts: 5 Rep Power: 5 I've changed the face name to "Core" for all boundary conditions but the problem remains. The interior wall of the Cylinder is heating the internal volume (I don't care what happens in this space), while the exterior isn't heating the surrounding air filled cavity.

 March 22, 2018, 08:44 #6 Senior Member     Uwe Pilz Join Date: Feb 2017 Location: Leipzig, Germany Posts: 588 Rep Power: 10 Now I understand. A patch has only one direction, because the other side is thought to be the boundary of the domain to be computed. To define it for two adjacent regions you have to note the patches twice, with an alternativ counting of nodes. __________________ Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)

 March 22, 2018, 08:52 #7 New Member   Join Date: Feb 2018 Posts: 5 Rep Power: 5 So how do I reverse the direction so the cylinder becomes outward facing, and computation occurs between the cavity/cylinder. Thank you very much for your help so far!

 March 22, 2018, 11:24 #8 Senior Member     Uwe Pilz Join Date: Feb 2017 Location: Leipzig, Germany Posts: 588 Rep Power: 10 If you don't need the region inside the cylinder you should not mesh it but let it empty. __________________ Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)

 March 23, 2018, 05:10 #9 New Member   Join Date: Feb 2018 Posts: 5 Rep Power: 5 That result in a blockMesh error? Code: ```Creating topology blocks Creating topology patches Creating block mesh topology --> FOAM FATAL ERROR: face 0 in patch 0 does not have neighbour cell face: 4(9 13 14 10) From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID) in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 124.```

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post [ANSYS Meshing] Meshing solid cylinder Sunlight311 ANSYS Meshing & Geometry 0 September 30, 2018 13:32 Scabbard OpenFOAM Running, Solving & CFD 37 March 21, 2016 16:16 Andy Singh FLUENT 5 August 31, 2004 07:31 Mirek Kabacinski FLUENT 0 July 23, 2003 18:40 Sriram Popuri Main CFD Forum 6 July 13, 1999 18:09

All times are GMT -4. The time now is 02:46.

 Contact Us - CFD Online - Privacy Statement - Top