local mesh refinement at channel geometries
 

Dear All Formers

I just saw you post for making local mesh refinement
in channel geometries

I want to refine grid only near wall , say, quarter of channel height at each wall.
It looks like cyclic bc prevent mesh generation. Have you any experience to this
problem?

I specify cellSet
name c0;

action new;

topoSetSources
(
boxToCell
{
box ( 0 0 0 ) ( 4 0.25 2 );
// box ( 0 1.75 0) (4 2 2);
}
);

and refineMeshDict

set c0;

coordinateSystem global;

globalCoeffs
{
tan1 ( 1 0 0 );
tan2 ( 0 1 0 );
tan3 ( 0 0 1 );
}

directions ( tan1 tan2 tan3);

useHexTopology yes;

geometricCut no;

writeMesh no;


But It show a error below. Do you have any idea for this remedy

Thanks in advance

Kyungrock

3D case; refining all directions

--> FOAM Warning :
From function syncTools<class T, class CombineOp>::syncEdgeList(const polyMesh&, UList<T>&, const CombineOp&, const T&, const bool)
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/syncToolsTemplates.C at line 1312
There are decomposed cyclics in this mesh with transformations.
This is not supported. The result will be incorrect
cyclicPolyPatch:rder : Writing half0 faces to OBJ file "sides1_half0_faces.obj"
cyclicPolyPatch:rder : Writing half1 faces to OBJ file "sides1_half1_faces.obj"
cyclicPolyPatch:rder : Dumping currently found cyclic match as lines between corresponding face centres to file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/sides1_faceCentres.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch:rder(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1547
Patch:sides1 : Cannot match vectors to faces on both sides of patch
Perhaps your faces do not match? The obj files written contain the current match.
Continuing with incorrect face ordering from now on!
cyclicPolyPatch::getGeometricHalves : Writing half0 faces to OBJ file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/sides2_half0_faces.obj"
cyclicPolyPatch::getGeometricHalves : Writing half1 faces to OBJ file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/sides2_half1_faces.obj"
cyclicPolyPatch::getGeometricHalves : Writing half0 face centres to OBJ file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/sides2_half0.obj"
cyclicPolyPatch::getGeometricHalves : Writing half1 face centres to OBJ file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/sides2_half1.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch::getGeometricHalves(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 438
Patch sides2 gets decomposed in two zones ofinequal size: 4015 and 3985
This means that the patch is either not two separate regions or one region where the angle between the different regions is not sufficiently sharp.
Please adapt the featureCos setting.
Continuing with incorrect face ordering from now on!
cyclicPolyPatch:order : Writing half0 faces to OBJ file "inout1_half0_faces.obj"
cyclicPolyPatch:order : Writing half1 faces to OBJ file "inout1_half1_faces.obj"
cyclicPolyPatch:order : Dumping currently found cyclic match as lines between corresponding face centres to file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/inout1_faceCentres.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch:rder(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1547
Patch:inout1 : Cannot match vectors to faces on both sides of patch
Perhaps your faces do not match? The obj files written contain the current match.
Continuing with incorrect face ordering from now on!
cyclicPolyPatch:order : Writing half0 faces to OBJ file "inout2_half0_faces.obj"
cyclicPolyPatch:order : Writing half1 faces to OBJ file "inout2_half1_faces.obj"
cyclicPolyPatch:order : Dumping currently found cyclic match as lines between corresponding face centres to file "/storage/kkim/OpenFOAM/OpenFOAM-1.6/tutorials/incompressible/channelFoam/channel395/inout2_faceCentres.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch:order(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1547
Patch:inout2 : Cannot match vectors to faces on both sides of patch
Perhaps your faces do not match? The obj files written contain the current match.
Continuing with incorrect face ordering from now on!
face 0 area does not match neighbour 4000 by 191.079% -- possible face ordering problem.
patch:sides1 my area:0.0510945 neighbour area:0.00116558 matching tolerance:0.001
Mesh face:1430800 vertices:4((0 0 0) (4 0.00480001 2) (0.05 0.00480001 0) (0.05 0 0))
Neighbour face:1434800 vertices:4((2.7 0.356633 0) (2.65 0.356633 0) (2.65 0.379945 0) (2.7 0.379945 0))
Rerun with cyclic debug flag set for more information.

From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 180.

FOAM exiting

Hi, Kyungrock

Just change the boundary type "cyclic" to other one, "wall", for example at first and refineMesh, and then change back to "cyclic" again, hoho.
Don't know whether solve your problem in this way, good luck!

Lijun

Hi,
finally, I tried the case channel395 by myself.
There is a warning, but no error appear.
attach it for u

Best wishes,
Lijun

Dear Jiang

I so much appreciate for your reply

Kyungrock

channel source files
 

Here is channel source files.

Thanks in advance

Kyungrock

Hi, Kyungrock

I tried your case this morning, and encounter the same problem. It 's a little strange.
If I change the domain refined to {(0.1 0 0) (3.9 0.25 2)}, It works and gives a good result.
Maybe the problem comes from the two cyclic boundary interface each other.
I am not sure about that.

I am sorry for that cann't solve your problem at present, but I will pay attention to this issue.

Best wishes,
Lijun

Hi Jiang

Dr. Eugene recommended me to use createPatch option after adjusting refineMesh with different boundary patch, not cyclic. I'll do that way. If you succeed in creating local mesh refinement, please let me know.

Thanks

Kyungrock