|
[Sponsors] |
[snappyHexMesh] Sphere in a channel by snappyHexMesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 14, 2016, 17:01 |
Sphere in a channel by snappyHexMesh
|
#1 |
Member
Arsalan
Join Date: Jul 2014
Posts: 74
Rep Power: 12 |
Dear Foamers,
I want to mesh a sphere inside a channel in both regions of inside and outside of the sphere. For this purpose, I'm trying to mesh these regions separately using sHM and merge them using mergeMesh.( Do you have any better suggestion??!) For the region outside of the sphere, I'm using Tobias' tutorial (many thanks to him for great tutorials) but with my defined regions not the stl surfaces. ( I also used the stl surfaces, but with no success!) My problem is snapping not works properly when a finer grid is used as the background mesh by blockMesh, or when I use high levels for refinement region around the sphere. And as a result of bad snapping, layering procedure not works! My sHM dict is as follows: Code:
// Which of the steps to run castellatedMesh true; snap true; addLayers true; geometry { sphere { type searchableSphere; centre (0 0 0); radius 0.0025; } }; // Settings for the castellatedMesh generation. castellatedMeshControls { maxLocalCells 1000000; maxGlobalCells 2000000; minRefinementCells 0; maxLoadUnbalance 0.0; nCellsBetweenLevels 3; features ( ); refinementSurfaces { sphere { // Surface-wise min and max refinement level level (3 3); } } resolveFeatureAngle 60; refinementRegions { } locationInMesh (0.01 0.01 0.01); allowFreeStandingZoneFaces no; } snapControls { nSmoothPatch 3; tolerance 1.2; nSolveIter 185; nRelaxIter 4; nFeatureSnapIter 5; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; } addLayersControls { relativeSizes true; layers { sphere { nSurfaceLayers 6; } } // Expansion factor for layer mesh expansionRatio 1.1; finalLayerThickness 0.7; minThickness 0.0001; nGrow 0; featureAngle 85; slipFeatureAngle 25; nRelaxIter 5; nSmoothSurfaceNormals 4; nSmoothNormals 3; nSmoothThickness 10; maxFaceThicknessRatio 0.5; maxThicknessToMedialRatio 0.2; minMedianAxisAngle 90; nBufferCellsNoExtrude 0; nLayerIter 50; } meshQualityControls { #include "meshQualityDict" //- Number of error distribution iterations nSmoothScale 4; //- amount to scale back displacement at error points errorReduction 0.75; } writeFlags ( scalarLevels layerSets layerFields // write volScalarField for layer coverage ); mergeTolerance 1e-6; Code:
convertToMeters 0.0005; vertices ( (-30 -30 -30) ( 30 -30 -30) ( 30 30 -30) (-30 30 -30) (-30 -30 30) ( 30 -30 30) ( 30 30 30) (-30 30 30) ); blocks ( hex (0 1 2 3 4 5 6 7) (60 60 60) simpleGrading (1 1 1) ); edges ( ); boundary ( maxY { type wall; faces ( (3 7 6 2) ); } minX { type wall; faces ( (0 4 7 3) ); } maxX { type wall; faces ( (2 6 5 1) ); } minY { type wall; faces ( (1 5 4 0) ); } minZ { type wall; faces ( (0 3 2 1) ); } maxZ { type wall; faces ( (4 5 6 7) ); } ); mergePatchPairs ( ); I really need to high levels of refinement around the sphere and a finer grid of the whole domain together. Any help and comment will be appreciated. Regards, Arsalan. |
|
Tags |
layer, shm, snap, snapping, snappyhexmesh |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] Meshing a sphere with snappyHexMesh (for beginners) | Tobi | OpenFOAM Meshing & Mesh Conversion | 7 | November 20, 2016 18:39 |
transient motion of sphere in a square channel | sajeesh | FLUENT | 0 | March 4, 2013 11:15 |
[snappyHexMesh] Meshing a sphere with snappyHexMesh | Cyberholmes | OpenFOAM Meshing & Mesh Conversion | 2 | July 19, 2011 17:46 |
meshing F1 front wing | Steve | FLUENT | 0 | April 17, 2003 13:37 |
Motion of sphere freely suspended in channel flow | faithkim | Main CFD Forum | 1 | August 12, 1999 11:42 |