CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

local mesh refinement at channel geometries

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2010, 17:10
Default local mesh refinement at channel geometries
  #1
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
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!
cyclicPolyPatchrder : Writing half0 faces to OBJ file "inout1_half0_faces.obj"
cyclicPolyPatchrder : Writing half1 faces to OBJ file "inout1_half1_faces.obj"
cyclicPolyPatchrder : 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!
cyclicPolyPatchrder : Writing half0 faces to OBJ file "inout2_half0_faces.obj"
cyclicPolyPatchrder : Writing half1 faces to OBJ file "inout2_half1_faces.obj"
cyclicPolyPatchrder : 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 cyclicPolyPatchrder(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
Kr_kim is offline   Reply With Quote

Old   February 8, 2010, 09:29
Default
  #2
New Member
 
Jiang Lijun
Join Date: Mar 2009
Posts: 14
Rep Power: 8
L.J.Jiang is on a distinguished road
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
L.J.Jiang is offline   Reply With Quote

Old   February 8, 2010, 12:18
Default
  #3
New Member
 
Jiang Lijun
Join Date: Mar 2009
Posts: 14
Rep Power: 8
L.J.Jiang is on a distinguished road
Hi,
finally, I tried the case channel395 by myself.
There is a warning, but no error appear.
attach it for u

Best wishes,
Lijun
Attached Files
File Type: gz channel375.tar.gz (3.5 KB, 22 views)
L.J.Jiang is offline   Reply With Quote

Old   February 8, 2010, 14:29
Default
  #4
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
Dear Jiang

I so much appreciate for your reply

Kyungrock
Kr_kim is offline   Reply With Quote

Old   February 8, 2010, 18:18
Default channel source files
  #5
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
Here is channel source files.

Thanks in advance

Kyungrock
Attached Files
File Type: gz channel375Test.tar.gz (20.3 KB, 18 views)
Kr_kim is offline   Reply With Quote

Old   February 9, 2010, 09:48
Default
  #6
New Member
 
Jiang Lijun
Join Date: Mar 2009
Posts: 14
Rep Power: 8
L.J.Jiang is on a distinguished road
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
L.J.Jiang is offline   Reply With Quote

Old   February 9, 2010, 15:25
Default
  #7
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
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
Kr_kim is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] Hexa mesh, curve mesh setup, bunching law Anorky ANSYS Meshing & Geometry 4 November 12, 2014 01:27
Auto Mesh refinement in CFX Kushagra CFX 4 June 1, 2008 21:00
Mesh Refinement Help...sort of Adam CFX 1 February 3, 2008 17:38
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
mesh refinement on top of existent mesh? jx FLUENT 10 January 11, 2004 05:32


All times are GMT -4. The time now is 14:12.