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

How to use Explicit feature edge refinement

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 12, 2010, 10:01
Default How to use Explicit feature edge refinement
  #1
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 129
Rep Power: 8
ziemowitzima is on a distinguished road
Hi Everyone !
I am doing flow in 3D abdominal aorta. To create mesh I utilized snappyHexMesh. But obtained strange jagged edges (see pictures). To solve this problem I tried to use Explicit feature edge refinement:
features
(
{
file "someLine.eMesh";
level 4;
}
);

In the file someLine.eMesh I defined outflow edge (set of 132 points going around it circumference) and then I definite edges (file listing below).
But this didnt improved anything. I think that it is because I defined 131 edges from point to point. Similar like in igloo example.
Does anyone know how to use it to define one edge as spline or polyLine for example ?
To let snappyHM know that this is one edge ?

Best Regards
ZM

FoamFile
{
version 2.0;
format ascii;

root "..";
case "cavity";
instance ""constant"";
local "polyMesh";

class featureEdgeMesh;
object points;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Points
(
( 3.046 1.8233 0.056377 )
( 3.0385 1.8612 0.064124 )
( 3.0356 1.8752 0.066909 )
( 3.0469 1.82 0.055774 )
( 3.0619 1.7737 0.04887 )
( 3.0318 1.9929 0.099428 )
( 3.0345 2.0318 0.111908 )
( 3.0359 2.0403 0.11499 )
( 3.031 1.9834 0.096357 )
( 3.0353 1.8769 0.067289 )
( 3.0312 1.9276 0.080246 )
( 3.0307 1.9367 0.082673 )
( 3.031 1.982 0.09594 )
( 3.0441 2.0823 0.130756 )
( 3.0487 2.1024 0.138624 )
( 3.0567 2.1327 0.150868 )
( 3.0813 2.182 0.175856 )
( 3.0815 2.1826 0.176144 )
( 3.0823 2.1841 0.176902 )
( 3.1911 2.3309 0.266749 )
( 3.2177 2.3529 0.284656 )
( 3.1855 2.3263 0.262974 )
( 3.2392 2.3674 0.298228 )
( 3.1809 2.3218 0.259684 )
( 3.1566 2.2974 0.242049 )
( 3.1305 2.267 0.221891 )
( 3.1303 2.2668 0.221749 )
( 3.1302 2.2666 0.221631 )
( 3.1079 2.2323 0.202024 )
( 3.2483 2.3735 0.303924 )
( 3.2569 2.3783 0.309053 )
( 3.2768 2.389 0.320822 )
( 3.3035 2.3998 0.33553 )
( 3.3116 2.4032 0.340016 )
( 3.3243 2.4068 0.346603 )
( 3.3478 2.4133 0.358684 )
( 3.3641 2.4156 0.366407 )
( 3.516 2.4075 0.429967 )
( 3.5371 2.4032 0.437892 )
( 3.4808 2.4146 0.416779 )
( 3.5936 2.384 0.45679 )
( 3.4297 2.4195 0.395989 )
( 3.4284 2.4196 0.395476 )
( 3.4271 2.4196 0.394908 )
( 3.3867 2.4188 0.377173 )
( 3.759 2.291 0.501544 )
( 3.787 2.2648 0.506088 )
( 3.725 2.3156 0.493952 )
( 3.7166 2.3214 0.491996 )
( 3.6734 2.3486 0.481132 )
( 3.7086 2.3266 0.490031 )
( 3.6034 2.3807 0.460095 )
( 3.6058 2.3798 0.460848 )
( 3.6097 2.3779 0.462023 )
( 3.7937 2.2585 0.507138 )
( 3.7992 2.2524 0.507798 )
( 3.8293 2.2171 0.510589 )
( 3.8589 2.1763 0.511574 )
( 3.859 2.1761 0.511581 )
( 3.8592 2.1759 0.511583 )
( 3.9017 2.1058 0.50966 )
( 3.9179 2.0739 0.507449 )
( 3.9269 2.0548 0.505791 )
( 3.8869 2.1321 0.510904 )
( 3.9382 2.0288 0.503165 )
( 3.955 1.9921 0.499763 )
( 3.9687 1.9541 0.494676 )
( 3.976 1.9328 0.491638 )
( 3.9803 1.9132 0.48782 )
( 3.0646 1.766 0.047793 )
( 3.0794 1.7261 0.042625 )
( 3.0842 1.7152 0.04151 )
( 3.0976 1.6909 0.040288 )
( 3.1062 1.6744 0.039228 )
( 3.114 1.6626 0.03917 )
( 3.1333 1.6321 0.038717 )
( 3.1588 1.5963 0.039361 )
( 3.1629 1.5909 0.039577 )
( 3.2047 1.5443 0.044177 )
( 3.1765 1.5762 0.041184 )
( 3.9877 1.8797 0.481313 )
( 3.9914 1.8482 0.473759 )
( 3.9933 1.8283 0.468823 )
( 3.992 1.7634 0.449387 )
( 3.9881 1.7303 0.438081 )
( 3.9842 1.7119 0.431042 )
( 3.9936 1.7781 0.454386 )
( 3.9938 1.7839 0.456113 )
( 3.9786 1.6847 0.420711 )
( 3.9713 1.66 0.410388 )
( 3.9661 1.6412 0.402653 )
( 3.959 1.6234 0.394398 )
( 3.9492 1.5991 0.383113 )
( 3.937 1.5792 0.372019 )
( 3.8547 1.4798 0.30745 )
( 3.8332 1.4616 0.292871 )
( 3.8171 1.4482 0.281961 )
( 3.876 1.5018 0.323061 )
( 3.7917 1.4327 0.266452 )
( 3.7706 1.4203 0.253671 )
( 3.7576 1.4136 0.246069 )
( 3.7303 1.3999 0.230251 )
( 3.8092 1.4427 0.276957 )
( 3.9246 1.5595 0.360921 )
( 3.8914 1.5179 0.33444 )
( 3.9117 1.5436 0.350713 )
( 3.7102 1.3928 0.219475 )
( 3.6805 1.3864 0.204743 )
( 3.6639 1.3819 0.196213 )
( 3.6427 1.3785 0.186032 )
( 3.6275 1.3768 0.17895 )
( 3.6001 1.3762 0.166923 )
( 3.5914 1.3759 0.163038 )
( 3.5795 1.3764 0.158025 )
( 3.5519 1.3783 0.146603 )
( 3.5157 1.3834 0.132337 )
( 3.5118 1.3839 0.130818 )
( 3.4715 1.3931 0.116023 )
( 3.4518 1.3982 0.10895 )
( 3.428 1.4042 0.100346 )
( 3.5049 1.3855 0.128322 )
( 3.4146 1.4105 0.096366 )
( 3.3836 1.4249 0.087127 )
( 3.3561 1.4383 0.079073 )
( 3.3434 1.4447 0.075436 )
( 3.3119 1.4627 0.066964 )
( 3.307 1.4655 0.065673 )
( 3.2914 1.4758 0.061877 )
( 3.2747 1.4867 0.057809 )
( 3.269 1.491 0.056605 )
( 3.2424 1.5109 0.050821 )
( 3.2135 1.5362 0.045639 )
)
// Edges
(
( 0 1 )
( 1 2 )
( 2 3 )
( 3 4 )
( 4 5 )
( 5 6 )
( 6 7 )
( 7 8 )
( 8 9 )
( 9 10 )
( 10 11 )
( 11 12 )
( 12 13 )
( 13 14 )
( 14 15 )
( 15 16 )
( 16 17 )
( 17 18 )
( 18 19 )
( 19 20 )
( 20 21 )
( 21 22 )
( 22 23 )
( 23 24 )
( 24 25 )
( 25 26 )
( 26 27 )
( 27 28 )
( 28 29 )
( 29 30 )
( 30 31 )
( 31 32 )
( 32 33 )
( 33 34 )
( 34 35 )
( 35 36 )
( 36 37 )
( 37 38 )
( 38 39 )
( 39 40 )
( 40 41 )
( 41 42 )
( 42 43 )
( 43 44 )
( 44 45 )
( 45 46 )
( 46 47 )
( 47 48 )
( 48 49 )
( 49 50 )
( 50 51 )
( 51 52 )
( 52 53 )
( 53 54 )
( 54 55 )
( 55 56 )
( 56 57 )
( 57 58 )
( 58 59 )
( 59 60 )
( 60 61 )
( 61 62 )
( 62 63 )
( 63 64 )
( 64 65 )
( 65 66 )
( 66 67 )
( 67 68 )
( 68 69 )
( 69 70 )
( 70 71 )
( 71 72 )
( 72 73 )
( 73 74 )
( 74 75 )
( 75 76 )
( 76 77 )
( 77 78 )
( 78 79 )
( 79 80 )
( 80 81 )
( 81 82 )
( 82 83 )
( 83 84 )
( 84 85 )
( 85 86 )
( 86 87 )
( 87 88 )
( 88 89 )
( 89 90 )
( 90 91 )
( 91 92 )
( 92 93 )
( 93 94 )
( 94 95 )
( 95 96 )
( 96 97 )
( 97 98 )
( 98 99 )
( 99 100 )
( 100 101 )
( 101 102 )
( 102 103 )
( 103 104 )
( 104 105 )
( 105 106 )
( 106 107 )
( 107 108 )
( 108 109 )
( 109 110 )
( 110 111 )
( 111 112 )
( 112 113 )
( 113 114 )
( 114 115 )
( 115 116 )
( 116 117 )
( 117 118 )
( 118 119 )
( 119 120 )
( 120 121 )
( 121 122 )
( 122 123 )
( 123 124 )
( 124 125 )
( 125 126 )
( 126 127 )
( 127 128 )
( 128 129 )
( 129 130 )
( 130 131 )
( 131 0 )


)
Attached Images
File Type: jpg inlet-patch.jpg (45.3 KB, 62 views)
File Type: jpg mesh-zoom.jpg (98.5 KB, 70 views)
ziemowitzima is offline   Reply With Quote

Old   March 30, 2010, 08:58
Default
  #2
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 277
Rep Power: 9
vinz is on a distinguished road
Hi,

I have the same problem, not with 132 points but with 12 which doesn't change anything I guess.
Placing this file in the features part doesn't change anything while it is working for the igllo tutorial.
Is there something special to do in our case? did you find a solution?
regards,

Vincent
vinz is offline   Reply With Quote

Old   March 30, 2010, 09:44
Default
  #3
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 129
Rep Power: 8
ziemowitzima is on a distinguished road
Yes,
I imported mesh from ICEM .
ziemowitzima is offline   Reply With Quote

Old   April 22, 2010, 09:16
Default
  #4
Member
 
Moritz Wied
Join Date: Mar 2010
Location: suttgart, germany
Posts: 35
Rep Power: 7
Mo-ITB is on a distinguished road
Hi,
i also have some issues with the feature edge refinement. probably it is just not working the way we expect it. for me it looks like it is working like that:

1. looking for a cell with and edge in it
2. refining this cell
3. looking for the refined cell, which contains the edge now.
4. refining this cell
5. ...

to be continued depending on the level ..; you define in the snappyHexMeshDict.

In other words, it makes the same as refinementSurfaces, so if you put for example

features
{
file "...";
level 3;
}



it would be the same as




refinementSurfaces
{
patch
{
level (2 3); <- here the second number
}

}

with a low resolveFeatureAngle to resolve all edges.

Is this right? If yes, anybody already has a trick to solve sharp edges?

all the best,
moritz
Mo-ITB 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
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43
Actuator disk model audrich FLUENT 0 September 21, 2009 07:06
Where's the singularity/mesh flaw? audrich FLUENT 3 August 4, 2009 01:07
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Star mesh import problem chris1980 OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 20 May 8, 2006 01:07


All times are GMT -4. The time now is 21:43.