CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > ANSYS Meshing & Geometry

[ICEM] trouble with small angles

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 3, 2010, 23:59
Default trouble with small angles
  #1
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Hello,


Iím having some difficulty meshing a geometry that has some small angles between surfaces. The geometry represents the topology of a river channel. I have provided the channel bottom and inlet and outlet in fig1. The water surface sits on top. The flow is just over the banks so that a small angle occurs between the water surface and the channel bottom.

For the first run, Iím trying to mesh this geometry with the tetra/mixed octree. I have set the edge criterion fairly low (0.05) and defined a thin cut between the surfaces. The resulting mesh is shown in fig2. There a several gaps in the mesh near the edges (an example is shown in fig3). Some of these gaps were able to be fixed with the ďfind/Close HolesĒ tool, others could not be fixed though. I would prefer not to lower the minimum size to keep the element count reasonable, as the geometry shown here is only about a quarter of the entire geometry. I have experimented with other tetra/mixed schemes and meshing a single surface then generating a volume mesh. However these results are no better. Maybe Iím going about this in the totally wrong way. It certainly wouldnít be the first time. Any help/suggestions are greatly appreciated,

john
Attached Images
File Type: jpg fig1.jpg (39.8 KB, 95 views)
File Type: jpg fig2.jpg (69.8 KB, 75 views)
File Type: jpg fig3.jpg (83.4 KB, 67 views)
johnp is offline   Reply With Quote

Old   November 4, 2010, 05:13
Default
  #2
Senior Member
 
AB
Join Date: Sep 2009
Location: France
Posts: 323
Rep Power: 12
BrolY will become famous soon enough
It looks like some curves are missing. Have you used the "bluild topology" option ?
BrolY is offline   Reply With Quote

Old   November 4, 2010, 06:51
Default
  #3
Member
 
Ludvik Lanicek
Join Date: Apr 2009
Location: Czech Republic
Posts: 51
Rep Power: 8
Ludvik is on a distinguished road
Try this:
In "Mesh > Volume Mesh Params (Octree)" decrease "Edge Criterion" from 0.2 to 0.0002 or "Define Thin Cuts" (see help).
Ludvik is offline   Reply With Quote

Old   November 4, 2010, 11:23
Default
  #4
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
thanks for the suggestions. I did use the "build topology" function and it found no problems. I further deleted some curves and points that were not necessary. I have defined the thin cuts as between two surfaces, the water surface and the channel boundary. Also, I decreased the edge criterion to 0.0002 and the problem persists.

I have begun trying to manually edit each area on the surface mesh to close the gaps. This is very time consuming, so I'm hoping that there is an easier way. Also, I'm not sure what will happen if I mesh the volume from the fixed surface mesh. Any thoughts?
johnp is offline   Reply With Quote

Old   November 4, 2010, 12:03
Default
  #5
Senior Member
 
AB
Join Date: Sep 2009
Location: France
Posts: 323
Rep Power: 12
BrolY will become famous soon enough
Could you post an image of your model with the mesh and the curves you kept?

The edge criteron is useful, but reducing it to 0,0002 seems very low to me, maybe too much !
BrolY is offline   Reply With Quote

Old   November 4, 2010, 12:38
Default
  #6
Member
 
Ludvik Lanicek
Join Date: Apr 2009
Location: Czech Republic
Posts: 51
Rep Power: 8
Ludvik is on a distinguished road
Open help and see Figure.20 in "Mesh>Global Mesh Setup> Volume Meshing Parameters> tetra/Mixed>Robust/Octree". This value is resolution for grid gereration.
On the other side. Some small examples are:
http://rapidshare.com/files/428877963/thin_cut_0.tin
http://rapidshare.com/files/428878091/thin_cut_2.tin
http://rapidshare.com/files/428878049/thin_cut_1.tin

"0" is with default settings -> wrong mesh,
"1" with thin cut -> good mesh and hold edge,
"2" with Edge criterion -> strange mesh but hold edge
Ludvik is offline   Reply With Quote

Old   November 4, 2010, 17:06
Default more figures
  #7
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Here are figures showing the mesh with curves and just the curves alone after running "Build Diagnostic Topology". The view is looking downstream from above the channel, flow is approximately in the positive x-direction.

The model contains four surfaces: inlet, outlet, water surface, and channel boundary. Could the fact that the two surface defined as a thin cut interesct twice (i.e. the water surface and boundary)? I'll try to spilt the surfaces and see what happens.
Attached Images
File Type: jpg mesh-with-curves.jpg (89.3 KB, 47 views)
File Type: jpg curves.jpg (26.7 KB, 36 views)
johnp is offline   Reply With Quote

Old   November 9, 2010, 12:21
Default
  #8
New Member
 
June
Join Date: Dec 2009
Posts: 18
Rep Power: 6
examosty is on a distinguished road
Hello John,
did you solve your problem? I have met the same trouble. It would be grateful to share your experience.
Best Regards.
examosty is offline   Reply With Quote

Old   November 9, 2010, 17:58
Default
  #9
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Hi Examosty,

sorry to say I haven't resolved the issue. I found that lowering the edge requirement to a point helped. However, after it was reduced past some critical value, strange things happened, such as elements violating geometry surfaces. I can't say for sure that this is related to the edge criterion, except that it didn't happen for larger values.

Presently, I'm experimenting with a hexa-mesh on a simplified geometry. I'll let you know if come up with a solution,

john
johnp is offline   Reply With Quote

Old   November 10, 2010, 04:17
Default
  #10
Member
 
Fabian E.
Join Date: Nov 2009
Posts: 33
Rep Power: 7
galap is on a distinguished road
I would suggest a hexa blocking strategy as the geometry looks adequate and easy for this. Why you didn't try this? I suggest 2 blocks - a thin block right under your water surface could handle the edges.
galap is offline   Reply With Quote

Old   November 10, 2010, 09:28
Default
  #11
New Member
 
June
Join Date: Dec 2009
Posts: 18
Rep Power: 6
examosty is on a distinguished road
Hi John,

My problem is that, i have a very small surface (the blue surface in the attached picture). The angle between two edges is about 3į. I donīt know if itīs possible to get a gut quality. Because of the complicate geometry I can not use hexa. I hope you solve it soon.

Regards,
june
Attached Images
File Type: jpg small_surface.jpg (19.4 KB, 35 views)
examosty is offline   Reply With Quote

Old   November 10, 2010, 11:58
Default hexa
  #12
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Galap,

I have been working on a hexa mesh the last two days. The reason I initially avoided it is due to a section of the geometry that is not shown in the figure. There is a bend in the channel and the cross-section shape is more triangular than trapezoidal.

The initial work indicates that hexa may work. Right now i've just got some very poor quality elements (determinant 2x2x2 < 0) to deal with. We'll how it goes today...

john
johnp is offline   Reply With Quote

Old   November 11, 2010, 18:31
Default Tips...
  #13
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,650
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
The mesh is not fine enough to capture the shallow river bank...

Yes, adjusting the edge criterion may help because it causes refinement instead of moving the nodes when the mesh tries to fit to the geometry. I can explain that algorithm if you need it...

Yes, Turning on "Thin Cuts" may help, and here I will explain the algorithm.

"if a surface element has a node on Part "A", it may not have any nodes on part "B".

This prevents surface mesh from jumping between the parts. However, if any two parts specified for a "thin cuts" relationship "touch", it becomes impossible for the elements along that connection to fulfill the algorithm. Therefore the thin cut is violated and turned off. If the surface of the water touches the river bed, you must have a third part between them (such as a curve called "shoreline") that could separate the two thin cut parts.

Both of the above methods will work better if you give in and go with a smaller mesh size. Actually, at this coarse mesh size, you couldn't expect much from the CFD solver.

But in the end, a Hexa mesh would be best quality for this geometry, if you can set it up... Just make sure that when you create the Ogrid, you Face the water surface, as well as the upstream and down stream cut offs. Other wise your ogrid will need to make an awkward and low quality turn at the shore line.
PSYMN is offline   Reply With Quote

Old   November 12, 2010, 09:42
Default finally...
  #14
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Thanks for all the suggestions. I've been working on a hexa mesh the last couple of days and, as many of you stated, it worked much better. I've attached the basic blocking structure which is a C-grid at each cross-section. To remove negative volume elements on the inner bank of the bend required some compromise in terms of element size transition. The last picture shows the location of elements with quality below 0.1. I know there are probably some improvements that can be made but I think it's a good start.
Attached Images
File Type: jpg blocking.jpg (63.4 KB, 31 views)
File Type: jpg bend.jpg (90.6 KB, 31 views)
File Type: jpg lowqual.jpg (88.6 KB, 31 views)
johnp is offline   Reply With Quote

Old   November 12, 2010, 10:23
Default Wedges?
  #15
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,650
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Great,

My guess is that the shoreline on the inner bank would be better represented with a wedge block than a warped hexa block...

You can get this by moving the shoreline vertex into the river to improve the corner and then adding a wedge block in for the shore line... (there are easier ways to get that topology if you were starting out with the plan, but this is how to get there from where you already are).

Get me an image of a scan plane with an index plane and verts with numbers, and I can tell you which verts to merge...

Best regards,

Simon
PSYMN is offline   Reply With Quote

Old   November 12, 2010, 13:44
Default
  #16
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Simon,


I had something like this in mind but I wasn’t able to get it to work. I have attached an image that I hope is what you requested (not sure what a scan or index plane are). The image shows the section that results in highly skewed elements. If you need something else, let me know. Thanks,


john
Attached Images
File Type: jpg blocks-bend.jpg (55.7 KB, 35 views)
johnp is offline   Reply With Quote

Old   November 14, 2010, 22:42
Default Step by step...
  #17
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,650
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
An index is a single plane of your blocking, such as the front plane in this image that contains the block verts 73-135-69-134. ICEM CFD hexa gets a lot easier if you know how to use your index control... (right click on the blocking branch of the tree to display the index, or hit the "i" hotkey while in the blocking tab. Index control appears in the lower left corner).

A scan plane is a good way to diagnose problems... Right click on "premesh" in the model tree to bring up the scan plane control (also in the bottom left). I usually use the option for "scan plane by selection"...

Anyway, this image is enough to talk about.

For the block I mentioned before, slide vertex 73 (and 105, etc.) under vertex 135 (and 138, etc.) to improve the angle of the corner... Of course, you will have to clean up your edge associations (the edge between 73 and 105 should now project to surface.)

But now, you are not properly representing that shore line... So we need to create some new blocks...

Go to the first icon of the blocking tab to initialize blocking... Then choose the second option on the DEZ to create blocks by verts. You have 4 verts, but are missing the last two... This makes the first block a little tricky... First change the pulldown option to degenerate (wedge). It now expects 6 verts... Give it the first two... Now 73 and 135, and then hit the middle mouse button... Normally this would mean you were done, but since you only picked 2 out of the 6 expected verts, it knows you are not done and switches to select location mode... Pick the geometry location where 73 used to be... Then pick near the locations where 105, 138 and 105 used to be... The tool was selecting by locations, so the verts near 105 and 138 are not actually connected... so you will need to merge those (Blocking => Merge => Merge Verts). From then on, you will have the base set and can pick your first 5 verts before you need to select the 6th geometry location... Follow the bed of the river and add these blocks one at a time or (skip a few indecies, split the block degenerate and merge it in.)

Note, there is a Y block option instead of Degenerate, but I don't think it will work very well (you will still have a hexa squeezed into the shore line)... You could also come back and change these wedge blocks into a stack of swept blocks if you prefer that sort of thing... (would avoid the stack of squished wedges at the shoreline)
PSYMN is offline   Reply With Quote

Old   November 15, 2010, 18:36
Default
  #18
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
Simon,


Thanks for the instructions on generating the wedge block. I think the blocking is correct (figs 1&2), but the prism (2D tri) elements are not aligned on the bank edge (fig 3). Is there a way to reorient the distribution to place the prism at the outer edge?

Also, eventually Iíll have to merge the wedge blocks with the C-grids in the straighter sections. Figs 4&5 show the current blocking, which obviously wonít work. I tried to split vertices (vert 73) to match, but didnít work. Is there a straight-forward way to do this? Thanks,

john
Attached Images
File Type: jpg block1.jpg (27.4 KB, 24 views)
File Type: jpg block2.jpg (28.2 KB, 19 views)
File Type: jpg mesh.jpg (56.4 KB, 22 views)
File Type: jpg 2blocks.jpg (37.5 KB, 19 views)
File Type: jpg 2blocks-mesh.jpg (77.9 KB, 21 views)
johnp is offline   Reply With Quote

Old   November 16, 2010, 17:25
Default Order of Nodes and 7 noded blocks...
  #19
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,650
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Oops, my bad. The focus point of the degenerate block is the first point you pic... This means you must start each wedge by selecting the point on the short line... Then things turn out ok.

WedgeBlock_johnp_1.jpg

I will ask around if there is a way to change this once the block is created, but I don't think there is.

As for what to do in other areas of the short line... I would normally just put wedges along the entire river... It should be fine quality wise and would be the best for bringing the refined boundary layer right up to the shoreline...

If you really wanted to transition back to square, you would still need to move the other blocks back into alignment, but you could put a hexa block (instead of a wedge) along the shoreline (I created a point at an XYZ coordinate for easy selection) and merge the one end to make a 7 noded hexa. Under Settings => Hexa Mixed, you will see an option that turns the corner collapsed hexas into pairs of pyramids during conversion to unstructured mesh...
WedgeBlock_johnp_3.jpg
WedgeBlock_johnp_2.jpg
PSYMN is offline   Reply With Quote

Old   November 17, 2010, 13:32
Default getting there...
  #20
Member
 
john
Join Date: Aug 2010
Location: blacksburg, va, usa
Posts: 36
Rep Power: 6
johnp is on a distinguished road
i couldn't find a way to change the tri location, so i just starting blocking from the beginning. i created the wedge before any associations, so i have a wedge plus C-grid around the entire domain.

It looks very nice except i have some negative volumes at the inner bank edge where the wedge elements are located (see figs 1 & 2). i assume this is because all the elements are being squished here, for lack of more technical phrasing. Just in case there were gaps, i set the Projection Limit to 0.1.

One solution i tried is to make the central H-block very thin and contain only a few elements (see figs 3 & 4). This works, but i don't like it. i feel like i'm sacrificing element quality in the area of interest to account for an area i'm not interested in. i could improve the mesh by moving some vertices and edge matching, etc. But i would rather find a better solution.

Another thought i had was to coarsen the wedge block in the vertical direction using Refinement. Unfortunately, it seem that i can't set the Level below 1.0, therefore no coarsening. Is there a better way to handle this? Thanks,

john
Attached Images
File Type: jpg block.jpg (79.5 KB, 24 views)
File Type: jpg prob-elements.jpg (41.0 KB, 17 views)
File Type: jpg bad-soln-blocking.jpg (40.3 KB, 18 views)
File Type: jpg bad-soln.jpg (74.7 KB, 15 views)
johnp 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
Time step too small? derz CFX 2 September 8, 2010 01:03
Small speed simpleFoam trouble milos OpenFOAM Post-Processing 0 November 17, 2009 09:00
Flows at High Angles of Attack Sarah C. FLUENT 2 March 18, 2009 09:39
Airfoil simulation in High Angles of Attack Yasser Nabavi FLUENT 0 April 21, 2006 18:28
small size holes simulation brian ko Phoenics 5 October 22, 2001 13:53


All times are GMT -4. The time now is 05:36.