|
[Sponsors] |
February 21, 2021, 20:06 |
Pyramids generated when extruding cylinder
|
#1 |
New Member
Join Date: Jun 2017
Posts: 3
Rep Power: 8 |
I am fairly new at generating meshes using OpenFoam's blockMesh and/or snappyHexMesh. I am attempting to generate the meshes for a geometry that consists of several cylinders moving in relation to one another. I intend to use dynamic overset meshes. As a first step, I tried to reproduce the example in .../tutorials/incompressible/overPimpleDyMFoam/cylinder using my own .stl geometry (trying to reconstruct the steps so I can do so by myself). Nevertheless, when extruding the mesh along the axis of the cylinder, what I obtain is a series of pyramidal volumes instead of hexahedrals. Does anybody know where the problem may be?
I've attached the images of my cylinder, with the pyramidal volumes along the axis of the cylinder, and that of the example, which does extrude the mesh properly. The code of extrudeMeshDict is the following: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1812 | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object extrudeMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // What to extrude: // patch : from patch of another case ('sourceCase') // mesh : as above but with original case included // surface : from externally read surface //constructFrom mesh; //constructFrom patch; constructFrom surface; // If construct from patch/mesh: sourceCase "../cavity"; sourcePatches (movingWall); // If construct from patch: patch to use for back (can be same as sourcePatch) exposedPatchName movingWall; // If construct from surface: surface "constant/triSurface/c0s.stl"; // Flip surface normals before usage. Valid only for extrude from surface or // patch. flipNormals false; //- Linear extrusion in point-normal direction extrudeModel linearNormal; //- Single layer linear extrusion in point-normal direction // with empty patches on front and back //extrudeModel plane; //- Linear extrusion in specified direction //extrudeModel linearDirection; //- Sector extrusion //extrudeModel sector; //- Wedge extrusion of a single layer // with wedge patches on front and back //extrudeModel wedge; //- Extrudes into sphere around (0 0 0) //extrudeModel linearRadial; //- Extrudes into sphere around (0 0 0) with specified radii //extrudeModel radial; //- Extrudes into sphere with grading according to pressure (atmospherics) //extrudeModel sigmaRadial; //- Extrudes by interpolating along path inbetween two (topologically identical) // surfaces (e.g. one is an offsetted version of the other) //extrudeModel offsetSurface; nLayers 20; expansionRatio 0.95; sectorCoeffs { axisPt (0 0 0); axis (0 -1 0); angle 360; } linearNormalCoeffs { thickness 4; } planeCoeffs { // thickness 0.1; nLayers 4; } linearDirectionCoeffs { direction (0 0 1); thickness 1; } /* linearRadialCoeffs { R 0.1; // Optional inner radius Rsurface 0.01; } radialCoeffs { // Radii specified through interpolation table R table ((0 0.01)(3 0.03)(10 0.1)); } sigmaRadialCoeffs { RTbyg 1; pRef 1; pStrat 1; } */ offsetSurfaceCoeffs { // Surface that mesh has been meshed to baseSurface "<constant>/triSurface/DTC-scaled-inflated.obj"; // Surface to fill in to offsetSurface "<constant>/triSurface/DTC-scaled.obj"; } // Do front and back need to be merged? Usually only makes sense for 360 // degree wedges. mergeFaces false; // Merge small edges. Fraction of bounding box. mergeTol 0; // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // |
|
March 4, 2021, 04:34 |
|
#2 |
Senior Member
M
Join Date: Dec 2017
Posts: 642
Rep Power: 12 |
I have no experience with extrude mesh, but I suspect there is no problem at all. Have you checked (checkMesh) that these are actually not hexahedral cells? You seem to use Paraview and when you have "decompose polyhedra" activated, it will show some of these cells split. Make sure its not just a representation issue in PV.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Flow past a 2D cylinder - High Re (1E+05) - Cd too high | Pervispasco | OpenFOAM Running, Solving & CFD | 4 | March 14, 2022 02:19 |
Flow through small cylinder recirculated inside larger cylinder. | Keesie | FLUENT | 1 | July 28, 2016 04:30 |
Oscillating cylinder in still water - interDyMFoam error | kl1812 | OpenFOAM Running, Solving & CFD | 2 | May 9, 2016 05:10 |
[ANSYS Meshing] complicated motion of a rotating cylinder? | Vamsi | ANSYS Meshing & Geometry | 0 | April 9, 2016 00:32 |
rotating cylinder using sliding mesh | Tim Daly | FLUENT | 1 | November 9, 2008 23:02 |