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

[ICEM] ICEM mesh problems for sinusoidal leading edge wing

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   March 18, 2013, 21:57
Default ICEM mesh problems for sinusoidal leading edge wing
  #1
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Hi everyone.

Ive only just started using ICEM (like 2 days ago) and have been trying to mesh a 3d wing with leading edge tubercles in ICEM. Essentially what I did was follow the 2D wing tutorial that everyone generally gets pointed to, and made it 3D with multiple blocks in the spanwise direction (someone suggested to do it this one on someone elses thread)

My problem is, that while my mesh looks okay superficially (kinda), when I run error checks it comes up with lots of errors, and weird stuff keeps happening when I try to import it into CFX. For instance my boundary conditions all change and become weird, and there are random floating blocks in domain.

Does anyone have any idea what is going on, or what I have done wrong? I know the mesh is not particularly pretty too, but I have been really trying just to get it to work.

If anyone could help me out and point me in the right direction, that would be great.

Here are my files. (both zip files have the same mesh, tubercle_mesh_1b just has the hex mesh. Figured I would give both so you did have to dl such a big file if you have slow internet speeds)

https://www.dropbox.com/s/7zpgq4ink8...le_mesh_1a.zip

https://www.dropbox.com/s/e30u7pzklr...le_mesh_1b.zip


Ans some pics of my mesh, and what is happening when I import it into ansys.


If anyone could help me out and point me in the right direction, that would be great.
Attached Images
File Type: jpg ICEM_mesh_1.jpg (92.8 KB, 30 views)
File Type: jpg ICEM_mesh_2.jpg (91.6 KB, 35 views)
File Type: jpg ANSYS_1.jpg (33.7 KB, 30 views)

Last edited by ShowponyStuart; March 19, 2013 at 05:00.
ShowponyStuart is offline   Reply With Quote

Old   March 18, 2013, 23:04
Default
  #2
Senior Member
 
Stuart Buckingham
Join Date: May 2010
Location: Melbourne, Australia
Posts: 253
Rep Power: 15
stuart23 will become famous soon enough
Your spanwise edges near the leading edge (on both sides) are associated to the leading edge curve. Delete these associations and try again. If you want to associate to the leading edge, create another split.
__________________
http://bc247.eng.monash.edu.au
stuart23 is offline   Reply With Quote

Old   March 19, 2013, 04:43
Default
  #3
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Quote:
Originally Posted by stuart23 View Post
Your spanwise edges near the leading edge (on both sides) are associated to the leading edge curve. Delete these associations and try again. If you want to associate to the leading edge, create another split.
That seemed to do the trick more or less, but I am still getting a couple of bad elements that are being weird and seem to be connected to the inlet and outlets (see pic)

I will do a whole pile of smoothing later, but it seems that it isnt set up well enough to start with and smoothing will just mask the poor quality a bit.

So does anyone have any advice on how to refine the mesh on the wing to capture the curvature better with out increasing the node count dramatically?

https://www.dropbox.com/s/7zpgq4ink8...le_mesh_1a.zip
Attached Images
File Type: jpg full_mesh.jpg (96.7 KB, 13 views)
File Type: jpg poor_surface_mesh.jpg (82.1 KB, 16 views)
File Type: jpg weird bits 2.jpg (16.7 KB, 15 views)

Last edited by ShowponyStuart; March 19, 2013 at 04:58.
ShowponyStuart is offline   Reply With Quote

Old   March 19, 2013, 05:50
Default
  #4
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Im thinking that I have done something wrong in the end section.

The first pic is coloured by connectivity, and the second picture shows that the blocks seem uneven at the ends, but I cant figure out if it is actually the problem or how to fix it.

Any thoughts?

https://www.dropbox.com/s/7zpgq4ink8...le_mesh_1a.zip
Attached Images
File Type: jpg connectivity.jpg (37.9 KB, 14 views)
File Type: jpg end section.jpg (55.3 KB, 16 views)
ShowponyStuart is offline   Reply With Quote

Old   March 19, 2013, 06:39
Default
  #5
Senior Member
 
Matthias Voß
Join Date: Mar 2009
Location: Berlin, Germany
Posts: 405
Rep Power: 10
mvoss is on a distinguished road
did you assoc. all faces and surfaces on the wing? try to highlight them by blocking-->faces-->"right click"-projected .. now you should see all blockfaces by color and name if not... do the assoc. and try again.
mvoss is online now   Reply With Quote

Old   March 19, 2013, 08:30
Default
  #6
Senior Member
 
Stuart Buckingham
Join Date: May 2010
Location: Melbourne, Australia
Posts: 253
Rep Power: 15
stuart23 will become famous soon enough
Matthias, I would refrain from associating faces unless absolutely needed. Unless there are internal faces, or faces close to each other (and the surface mesh is snapping to the wrong face), you do not need to do face associations. I think by making unneeded associations, you can over complicate the model, and sometimes it causes strange behavior if you modify the topology downstream.

Interested to see what you think though?
__________________
http://bc247.eng.monash.edu.au
stuart23 is offline   Reply With Quote

Old   March 19, 2013, 08:54
Default
  #7
Senior Member
 
Matthias Voß
Join Date: Mar 2009
Location: Berlin, Germany
Posts: 405
Rep Power: 10
mvoss is on a distinguished road
there is a point..yes. i was pointing that out because of the fact that at the outermost face the wingsurface is probably crossing the blockfaces. so projecting the faces would be my first try. btw. there is a premesh blocking option "project to faces" instead of "project to edges". maybe this would help correctly mapping the wingtip.
mvoss is online now   Reply With Quote

Old   March 19, 2013, 10:28
Default
  #8
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Hey Stuart junior

Some changes in geometry and blocking as usual as the blocking shown in Simon's tutorial was for sharp trailing edge wing/airfoil. Domain is changed to rectangle. Few curves were constructed to ease the blocking/association. See attached files.


https://dl.dropbox.com/u/68746918/si...alwing_Far.zip
aerosaravanan likes this.
Far is offline   Reply With Quote

Old   March 19, 2013, 21:58
Default
  #9
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Quote:
Originally Posted by neewbie View Post
there is a point..yes. i was pointing that out because of the fact that at the outermost face the wingsurface is probably crossing the blockfaces. so projecting the faces would be my first try. btw. there is a premesh blocking option "project to faces" instead of "project to edges". maybe this would help correctly mapping the wingtip.
I did attempt to find these and give them a try, but I didnt have much success. I am still learning ICEM so my knowledge is less than extensive haha.

Quote:
Originally Posted by Far View Post
Hey Stuart junior
Some changes in geometry and blocking as usual as the blocking shown in Simon's tutorial was for sharp trailing edge wing/airfoil. Domain is changed to rectangle. Few curves were constructed to ease the blocking/association. See attached files.
Thanks for that Far, that mesh seems too work quite well. I am still learning ICEM so I was wondering if you could quickly step me through what you did for that. So do you not need all the blocks in the spanwise direction to capture the curvature?

Ultimately I want to make a script for it so I can put various wings in (with different leading edge profiles) at different angles of attack easily. And the way you meshed it seems to work quite well.
ShowponyStuart is offline   Reply With Quote

Old   March 20, 2013, 03:35
Default
  #10
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
I have tried to copy the method Far used closely as I could and I think I have gotton reasonably close considering (see pic1 and 2) but I am still having a few problems.

The mesh is still reasonably nice but I am still getting a few anomalies in ansys where I think there are stray elements or something.

1- In pic 3 you can see an arrow. This arrow seems to be associated with the inlet somehow and im not really sure why (same thing happens with outlet to but arrow points in the opposite direction.)

2 - I cant figure out how to get the boundary conditions to be smooth in Ansys rather than clumping up where the mesh is more dense like you can see in pic 4.

Any thoughts on these issues?

Here are the most recent files if anyone wants a look to see what I may have done wrong, or see how the mesh could be improved.
https://www.dropbox.com/s/619byclo0j...mesh_2_aoa.zip
Attached Images
File Type: jpg pic_0.jpg (80.2 KB, 10 views)
File Type: jpg pic1_overview.jpg (96.2 KB, 10 views)
File Type: jpg pic2_inlet_random_bit.jpg (30.8 KB, 10 views)
File Type: jpg pic3_inlet_clumping.jpg (22.5 KB, 8 views)

Last edited by ShowponyStuart; March 20, 2013 at 04:02.
ShowponyStuart is offline   Reply With Quote

Old   March 20, 2013, 04:03
Default
  #11
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
More arrows mean there are more no of nodes. Fluent / CFD display the velocity field at those points (either in centre or cell faces). If you want to make them smooth then dont copy mesh sizing on far field edges. Instead select them separately and specify uniform spacing.

Also increase the size of edges at these ends some thing like twice or more. In this way you will get uniform mesh at inlet and outlet and hence flow field.


But that all should not affect your solution much if aspect ratio is within range (1000 for single precision and 10,000 for double precision). But you need to check by running both meshes with same conditions.
Far is offline   Reply With Quote

Old   March 20, 2013, 04:08
Default
  #12
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by neewbie View Post
did you assoc. all faces and surfaces on the wing? try to highlight them by blocking-->faces-->"right click"-projected .. now you should see all blockfaces by color and name if not... do the assoc. and try again.
By default ICEM project faces to surfaces. This is done for blocks at boundaries either internal or outer. For internal boundaries such as wing in this case, when you delete block or change its material, faces are automatically projected to closest surface.

Explicit face to surface association is needed when you want to make wall inside the flow domain.
Far is offline   Reply With Quote

Old   March 20, 2013, 04:18
Default
  #13
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Quote:
Originally Posted by Far View Post
Explicit face to surface association is needed when you want to make wall inside the flow domain.
Could that be the problem that is giving me the random arrow at the wing when I define my inlet and outlets? Or is that a separate issue?
ShowponyStuart is offline   Reply With Quote

Old   March 20, 2013, 04:20
Default
  #14
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
About wing I am not sure what it is showing. Run a sample simulation and see what happens. I think it is not important if wing is defined as static wall.
Far is offline   Reply With Quote

Old   March 20, 2013, 06:59
Default
  #15
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Quote:
Originally Posted by Far View Post
About wing I am not sure what it is showing. Run a sample simulation and see what happens. I think it is not important if wing is defined as static wall.
I did run a sim on both the one you did and mine and yours gave much better results I think.

I am trying to get the script working so I can start making a few of them to test but I just can get it working. What are some common problems people have with scripting? Did you use the replay control when you did the mesh you posted a few posts up?
ShowponyStuart is offline   Reply With Quote

Old   March 21, 2013, 05:12
Default
  #16
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
I think I have a mesh that works okay now https://www.dropbox.com/s/ti1nlnf6js...OA_refined.zip
(fundamentally I mean. i.e. No random elements or weird problems when its imported to ansys)

But even though I have a pretty refined mesh, I am still getting convergence problems with my lift and drag (even at relatively low angles of attack that are below the stall point the values oscillates up and down), which is weird because my residuals are still converging to 10^03 to 10^-4 or better. Note: I know in the pics I took the screenshot a couple of iterations in, but the same thing happened on multiple other sims, it will just keep oscillating like that.

So atm im trying to figure out if it is a problem with my mesh quality thats causing it or if it is becomes of some transient behavior (which would surprise me greatly considering that other people have simulated things like this with RANS before)
Attached Images
File Type: jpg 1_overview.jpg (73.6 KB, 10 views)
File Type: jpg 2_wing.jpg (61.3 KB, 14 views)
File Type: jpg residuals.jpg (56.7 KB, 11 views)
File Type: jpg lift_and_drag.jpg (51.7 KB, 5 views)
ShowponyStuart is offline   Reply With Quote

Old   March 21, 2013, 05:17
Default
  #17
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
place inlet at twice the distance presently it is from leading edge. What are the boundary conditions you are giving it?
Far is offline   Reply With Quote

Old   March 21, 2013, 05:35
Default
  #18
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
Boundary Conds:
Inlet ~76m/s zero gradient air@25C (gives a Reynolds number of 3million)
Side walls are symmetry
top walls are free slip walls
wing is no-slip wall
outlet is average static pressure outlet.

SST with gamma-theta transitional turbulence (high resolution for numerics and advection) (which worked well for my baseline wing)

Is there any easy way to move the inlet without rebuilding the thing from scratch. I cant get the replay control to work properly, it will do some of it then fail or do something random. Its really frustrating having to start from scratch every time I want to use a different wing or change the angle of attack.
ShowponyStuart is offline   Reply With Quote

Old   March 21, 2013, 05:43
Default
  #19
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,874
Blog Entries: 6
Rep Power: 37
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Just change geometry and update blocking (less than 100 seconds work ). Output mesh. I believe you have yplus below 1 every where.

Flow is transitional at Re = 3 Million !
Far is offline   Reply With Quote

Old   March 21, 2013, 06:16
Default
  #20
Member
 
Stuart
Join Date: Jun 2012
Posts: 59
Rep Power: 4
ShowponyStuart is on a distinguished road
I think I got it to work by doing that way (im not very experienced) but I ended up with a whole pile of blocks around the edge but they didnt seem to be hurting anything so I left them (mesh went crazy if I deleted them)

Yeah, thats why im using the gamma-theta transitional turbulence model. Im hoping that it will play nice and give me what I want haha It was working pretty well with the standard wing. (unfortunately im still fighting with the mesh at the moment)

Will run the sim with the longer domain now to see if that helps.
ShowponyStuart 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] Problems with export 2D mesh from ICEM to FLUENT newcomer ANSYS Meshing & Geometry 53 August 18, 2014 05:17
[ICEM] Meshing sinusoidal leading edge wing in ICEM ryerson_guy ANSYS Meshing & Geometry 39 April 10, 2013 02:57
wing with sharp trailing edge mesh through extrusion Hybrid Pointwise & Gridgen 8 March 18, 2013 17:17
Unstructured Mesh ICEM on a cube jerome_ ANSYS 0 May 30, 2012 05:34
[ICEM] 3d mesh around a wing using ICEM artemis64s ANSYS Meshing & Geometry 6 September 23, 2011 07:59


All times are GMT -4. The time now is 08:07.