CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Native Meshers: snappyHexMesh and Others (http://www.cfd-online.com/Forums/openfoam-meshing-snappyhexmesh/)
-   -   How to use Explicit feature edge refinement (http://www.cfd-online.com/Forums/openfoam-meshing-snappyhexmesh/72637-how-use-explicit-feature-edge-refinement.html)

ziemowitzima February 12, 2010 10:01

How to use Explicit feature edge refinement
 
2 Attachment(s)
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 )


)

vinz March 30, 2010 08:58

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

ziemowitzima March 30, 2010 09:44

Yes,
I imported mesh from ICEM :).

Mo-ITB April 22, 2010 09:16

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


All times are GMT -4. The time now is 04:58.