snappyHexMeshについて

[shin o-nishi]

OpenFoamユーザーの皆様

いつもお世話になっております

前回の質問から日を置かずに連続の質問申し訳ありません

２次元円柱周りの流れ解析を行っているのですが

メッシュの解像度を上げるとsnappyHexMeshで添付の画像のように

奇妙なメッシュが切れてしまいます

以前は円柱周りに細分化領域を作らずにメッシュを切ると

普通に切れるのですが

今回のように細分化領域を作ったところ

添付画像のような結果となりました

maxLocalCells maxGlobalCells

の値を大きくしたり

minTetQuality 

を-1e30にしたりなど

いろいろと試してみましたが

同じ結果となりました

ご存知の方どうかよろしくお願いします

注）blockmeshはx方向に6ｍを600分割

　　y方向に3ｍ300分割

　　z方向はextrudeMeshで一層を取り出し厚さは0.005です

以下はsnappyhexmeshdictの設定です

castellatedMesh true;

snap true;

addLayers true;

geometry

{

cylinder

{

type searchableCylinder;

point1 ( 0 0 1.0 );

point2 ( 0 0 -1.0 );

radius 0.05;

}

refinementBox

{

type searchableBox;

min (-0.3 -0.3 -1.0);

max ( 1.5 0.3 1.0);

}

};

castellatedMeshControls

{

maxLocalCells 10000000;

maxGlobalCells 200000000;

minRefinementCells 10;

maxLoadUnbalance 0.10;

nCellsBetweenLevels 3;

refinementSurfaces

{

cylinder

{

level (1 1);

}

}

refinementRegions

{

refinementBox

{

mode inside;

levels ((1E15 2));

}

}

locationInMesh (-0.8 0 0);

allowFreeStandingZoneFaces true;

}

snapControls

{

nSmoothPatch 6;

tolerance 4.0;

nSolveIter 0;

nRelaxIter 10;

//nFeatureSnapIter 10;

}

addLayersControls

{

layers

{

"cylinder.*"

{

nSurfaceLayers 10;

}

}

expansionRatio 1.0;

finalLayerThickness 1.0;

minThickness 0.5;

nGrow 0;

featureAngle 30;

nRelaxIter 3;

nSmoothSurfaceNormals 1;

nSmoothNormals 3;

nSmoothThickness 10;

maxFaceThicknessRatio 0.5;

maxThicknessToMedialRatio 0.3;

minMedianAxisAngle 90;

nBufferCellsNoExtrude 0;

nLayerIter 50;

}

meshQualityControls

{

maxNonOrtho 65;

maxBoundarySkewness 20;

maxInternalSkewness 4;

minVol 1e-13;

minTetQuality 1e-30;

minArea -1;

minTwist 0.02;

minDeterminant 0.001;

minFaceWeight 0.02;

minVolRatio 0.01;

minTriangleTwist -1;

nSmoothScale 8;

errorReduction 0.75;

}

debug 0;

mergeTolerance 1e-6;