CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] How to use Explicit feature edge refinement

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 12, 2010, 09:01
Default How to use Explicit feature edge refinement
  #1
Senior Member
 
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 145
Rep Power: 17
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, 146 views)
File Type: jpg mesh-zoom.jpg (98.5 KB, 158 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: 283
Rep Power: 18
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: 145
Rep Power: 17
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: 16
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


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
how to set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 06:09
[snappyHexMesh] Meshing a very thin meatus with SnappyHexMesh enginpower OpenFOAM Meshing & Mesh Conversion 29 September 9, 2017 03:46
[snappyHexMesh] problems using snappyHexMesh 2.1.0 on a supercomputer Sunxing OpenFOAM Meshing & Mesh Conversion 9 September 20, 2014 09:30
[snappyHexMesh] snappyHexMesh refinement regions ignored guitarbren OpenFOAM Meshing & Mesh Conversion 2 April 9, 2013 03:59
[snappyHexMesh] Weird snapping outside a refinement region jasimpson89 OpenFOAM Meshing & Mesh Conversion 0 July 29, 2012 20:46


All times are GMT -4. The time now is 16:56.