Edge and Face Meshing

Transcription

1 dge and Face Meshing 5-1

2 Meshing - General To reduce overall mesh size, confine smaller cells to areas where they are needed Locations of large flow field gradients. Locations of geometric details you wish to resolve. Controlling cell size distribution dges, faces, and volumes can be meshed directly. A uniform mesh is generated unless pre-meshing or size functions are used. Pre-meshing dge meshes can be graded (varying interval size on edge) A graded edge mesh can be used to control the cell size distribution of a face mesh. Controlling distribution of cell size on face mesh also controls the cell size distribution of the volume mesh. ize functions and boundary layers Allow direct control of cell size distribution on edges, faces and volumes directly for automatic meshing. 5-2

3 dge Meshing dge mesh distribution is controlled through the spacing and grading parameters on the Mesh dges form. Picking Temporary graphics Links, Directions Grading/pacing pecial characteristics Apply and Defaults Invert and Reverse Options 5-3

4 dge Meshing ense ense is used to show direction of grading very picked edge will show its sense direction using an arrow The sense can be reversed by a shift+middle-click on the last edge picked (this is in addition to the next functionality) or by clicking the Reverse button dge mesh preview When you pick an edge, the edge mesh is displayed using white nodes. This is a temporary mesh that has not been applied to the edge. Displayed edge mesh is based on current grading and spacing parameters If you modify the grading or spacing, the temporary mesh will be updated immediately. Meshing the edge The edge mesh is generated by clicking the Apply button. The nodes will then be displayed in blue. 5-4

5 Grading Controls mesh density distribution along an edge. Grading can produce single-sided or double-sided mesh Doubled-sided mesh can be symmetric or asymmetric. ymmetric schemes produce symmetric mesh about edge center. Asymmetric schemes can produce asymmetric mesh about edge center. ingle-sided grading: Uses a multiplicative constant, R, to describe the ratio of the length of two adjacent mesh elements: ingle-sided grading ymmetric grading Asymmetric grading R can be a user-specified value (uccessive Ratio) or calculated by GAMBIT. GAMBIT also uses edge length and spacing information to determine R. 5-5

6 Double ided Grading Double-sided grading can be generated by activating the double sided option in the Mesh dge form. Asymmetric grading is possible when the double-sided option is used with: uccessive Ratio, First Length, Last Length, First-Last Ratio, and Last-First Ratio The mesh is symmetric if R 1 = R 2 The mesh is asymmetric if R 1 R 2. dge center is determined automatically. ome schemes implicitly generate double sided grading that is symmetric. 5-6

7 oft Links Picking and soft links Pick with links By enabling this option, soft-linked edges can be selected in a single pick Linked edges share the same information and can be picked in a single pick Modifying soft links At any time, you can Form links Break links Maintain links By default, GAMBIT will form links between unmeshed edges that are picked together By default, GAMBIT will maintain links between meshed edges that are picked together 5-7

8 pacing In all meshing forms, the following spacing functions can be specified: Interval Count (recommended for edge mesh only) xample ntering a value of 5 will create 5 intervals along the selected edge(s) (6 nodes, including end nodes) Interval ize (default setting) Requires input of distance between nodes. dge is meshed with average interval size if grading is used. xample: An edge length of 10 and a value of 2 creates 5 intervals on the edge hortest dge % Meshes the selected edge according to a percentage of the length of the shortest edge in the model. xample hortest edge in model has length of 1. ntering a value of 20 will create a mesh with interval size

9 First dge ettings Use First dge ettings option If enabled: First edge selected in pick list updates all entries in the form. This mode is useful to copy settings from one meshed edge to other edge(s). If disabled: Use this setting any time you pick two or more meshed edges where there is a difference in type or spacing. The local Apply button for that option will be turned off This allows you to maintain pre-existing grading and/or spacing settings for each edge. nforce a change in grading and/or spacing by enabling Apply button. 5-9

10 Meshing Options Mesh This option is useful in cases where you want to impose a scheme without prescribing the number of intervals The higher level meshing scheme will decide (and match) the intervals Remove old mesh Deletes old mesh When selected, option to also delete lower geometry mesh appears. Ignore size function Toggle to either obey or ignore size functions ize function takes precedence when this option is disabled. 5-10

11 Meshing Options xample 1 pecify interval size, no grading, apply without meshing 3 Generate face mesh. Face Mesh Generated Using Quad Pave cheme (Pave face meshing schemes require an even number of elements on edge meshes) 2 pecify grading only, apply without meshing Face Mesh Generated Using ubmap cheme 5-11

12 Face Meshing Mesh Faces form Upon picking a face: GAMBIT automatically chooses quad elements GAMBIT chooses the type based on the solver/face vertex types Available element/scheme type combinations Quadrilaterial: Map, ubmap, Tri-Primitive, Pave Triangular: Pave Quad/Tri (hybrid): Map, Pave, Wedge Quad-to-tri conversion utility. 5-12

13 Face Meshing - Quad xamples Quad: Map Quad: ubmap Quad: Tri-Primitive Quad: Pave 5-13

14 Face Meshing - Quad/Tri and Tri xamples Quad/Tri Map Quad/Tri Wedge Face must be split to generate more than one cell across Quad/Tri Pave Tri Pave Triangular cell Quad cells Triangular cell 5-14

15 Deleting Old Mesh xisting mesh must be removed before remeshing. Mesh can be deleted using delete mesh form. Lower topology mesh can also be deleted (default) Alternatively, existing mesh can be deleted by selecting the Remove Old Mesh option Remove old mesh alone will leave all lower topology mesh Remove old mesh + remove lower mesh will delete all lower topology mesh that is not shared with another entity Undo after any meshing operation also works. 5-15

16 Face Vertex Types All vertices that are connected to a face are assigned initial face vertex types based on the angle between the edges connected to the vertex. Vertices shared by multiple faces can have multiple types, depending on which face you are considering. φ R The combination of vertex types describes the topology of the face. Face vertex types are used automatically to determine all quad face meshing schemes (except quad pave and tri pave). 5-16

17 Face Vertex Types nd () o 0 < φ <120 zero internal grid lines o ide () o 120 < φ < 216 One internal grid line o Corner (C) o 216 < φ < 309 Two internal grid lines o C C C Reverse (R) Three internal grid lines o 309 < φ < 360 o R R 5-17

18 Modifying Face Vertex Types Face vertex types can be changed from their default settings: Automatically By enforcing certain meshing schemes in face and volume meshing. Can sometimes result in undesirable mesh. Manually By direct modification in the Face Vertex Type form. elect Face symbols appear in graphics window elect New Vertex Type elect Vertices to be affected Vertex Types can be applied to just Boundary Layers as option. A vertex can have multiple types; one for each associated face. For a given set of face vertex types, GAMBIT will choose which meshing scheme to use based on predefined formulae. 5-18

19 Formula for ubmap cheme A face can be made submappable by manually changing vertex types Consider which vertex should be changed to type (side) In the et Face Vertex Type form, change vertex type to by enforcing the submap scheme. In the Face Mesh form, change the scheme from default to submap and click Apply. GAMBIT will attempt to change the vertex types so that the scheme is honored. User has less control the resulting mesh may be undesirable! R ubmap ( 4 nd) + ide + [( 2 nd) + Reverse] Which vertex to change? R R 5-19

20 Map ( 4 nd) + ( n ide) Formula for Map cheme + Periodic Map n ide Project intervals can be specified for more control. 5-20

21 How to Make a Face Mappable nforce the Map scheme (most common method) In the Face Mesh form, change the scheme from default to Map and click Apply. GAMBIT will attempt to change the vertex types so that the chosen scheme is honored Map ( 4 nd) + ( n ide) C C Manually change the vertex types In et Face Vertex Type form, change vertices (default) to "ide. Open the Face Mesh form and pick the face. GAMBIT should automatically select the map scheme) Map ( 4 nd) + ( n ide) 5-21

22 Formula for ubmap chemes ( ) ( ) ( ) ( ) ubmap: 4 nd + m ide + n nd + Corner + p Reverse + 2nd (additional terms when interior loops exist) C C Periodic ubmap n ide + m nd + Corner where m > 2. (additional terms when interior loops exist) ( ) ( ) C C C C C C C C C C C 5-22

23 Formula for Tri-Primitive heme Tri-Primitive ( 3 nd) + ( n ide) To mesh a rectangular face with the tri-primitive scheme: Manually change one of the vertex types to "ide" in this example The Tri Primitive scheme can not be enforced 5-23

24 Meshing Faces with Hybrid Quad/Tri chemes Quad/Tri: Tri-Map 2 Triangle The face vertex types must be changed manually to Trielement (T) The Tri-Map scheme must be selected. T T T Quad/Tri: Pave All vertex types are ignored except Trielement (T) and Notrielement (N) Trielement (T) will force a triangular element. Notrielement (N) will avoid a triangular element. C Quad/Tri: Wedge Used for creating cylindrical/polar type meshes The Vertex marked (T) is where rectangular elements are collapsed into triangles N 5-24 T

25 Hard Links Mesh linked entities have identical mesh Created for periodic boundary conditions Applicable to edge, face, and volume entities Best to use soft links for edge meshing To link volume meshes, all faces must be hard linked first. Hard links for faces elect faces and reference vertices The sense of each edge appears. Reverse orientation on by default Periodic option should be used for periodic boundary conditions, which creates a matched mesh even if the edges are split differently. Meshing one of the faces either before or after hard linking will generate an identical mesh on the linked face. Multiple pairs of hard links can be created. 5-25

26 Mesh moothing moothing can increase mesh quality beyond that of the default meshing algorithms Most noticable in complicated geometry. May have little or no effect in simple geometries. Mesh smoothing algorithms adjust interior node locations to obtain marginal improvement in mesh quality. Boundary meshes are not altered. The mesh at the boundary is not altered. Face and volume meshes are smoothed using a default scheme. Different schemes can be selected and applied after meshing. Face mesh smoothing Length-weighted Laplacian: Uses the average edge length of the elements surrounding each node to adjust the nodes. Centroid Area: Adjust node locations to equalize areas of adjacent elements. Winslow (quad meshes only): Optimizes element shapes with respect to perpendicularity. Volume mesh smoothing Length-weighted Laplacian: same as for face mesh smoothing quipotential: Adjusts node locations to equalize the volumes of the mesh elements surrounding each node. 5-26

27 xamining the Mesh Display Type Plane/phere View mesh elements that fall in plane or sphere. Range View mesh elements within quality range. Histogram shows quality distribution. how worst element zooms the view to the worst element elect 2D/3D and element type elect Quality Type Display Mode Change cell display attributes how Worst lement Automatically zooms the display to the worst element (based on current settings). Update button Will update values reported in the panel when options are changed. 5-27

28 Assessing Mesh Quality GAMBIT has several methods for assessing mesh quality. Aspect Ratio Diagonal Ratio dge Ratio quiangle kew quiize kew MidAngle kew ize Change tretch Taper Volume The most important of these quality metrics are quiangle kew and ize Change. 5-28

29 Mesh Quality quiangle kew The most important mesh quality metric is quiangle kew (Q A ). QA = θ θ max 180 θe max e, θ e θ θ e min θ max θ max θ θ max min θ e = Largest angle in face or cell = mallest angle in face or cell = Angle for equiangular face or cell θ min θ min Range of quiangle kew values 0 (best) 1 (worst) Q A = 0 describes a perfectly orthogonal element Q A = 1 describes a degenerate element θ e Quad/Hex o = 90 θ e θ e Tri/Tet = 60 θ e o 5-29

30 Mesh Quality ize Change Another important mesh quality metric is ize Change (QC). [ r, r, ] Q = max 2, C 1 K r n r i = Area or Volume of element i Area or Volume of neighbor element j V j=3 V j=2 Vi V j= 1 This metric applies only to 3D elements. V j=4 By definition, Q C > 0. 3D xample Q C,i = V j=1 /V i since j=1 has largest volume ratio Q C = 0 describes an element whose neighbor elements have exactly the same volume as the element of interest (i.e. uniform mesh). 5-30

31 triving for Quality A poor quality grid can cause inaccurate solutions and/or slow convergence. Minimize quiangle kew: Hex, Tri, Quad: kewness for all/most cells should be less than Tetrahedral: kewness for all/most cells should be less than 0.9. All elements: ize Change for cells in regions of interest should be less than 2 Minimize local variations in cell size, such as large jumps in size between adjacent cells. If xamine Mesh shows such violations: Determine the reason(s) for the violations Differences in spacing and grading on adjacent edges Geometry with small features or other defects Geometric complexity and size Mesh that grows too rapidly Delete mesh completely or partially. Clean and/or decompose geometry, premesh edges and faces or adjust meshing parameters Remesh the domain. 5-31

32 Appendix 5-32

33 Mesh Quality Aspect Ratio The Aspect Ratio metric (Q AR ) applies to tri, tet, quad, and hex elements and is defined differently for each element type. The definitions are as follows: Tri/Tet Inscribed circle Circumscribed circle Q AR = f is a scaling factor R and r are radii of circles (tri elements) or spheres (tet elements) that inscribe and circumscribe the mesh element. f = 1/2 for tri elements and f = 1/3 for tet elements. r f R r R 5-33 Quad/Hex Q AR max = min [ e1, e2, K, en ] [ e, e, K, e ] e i is the average length of edges in a coordinate system local to the element. N is the number of coordinate directions associated with the element d N = 2 for quad elements N = 3 for hex elements c a b 1 e e N a + c = 2 b + d = 2

34 Mesh Quality Diagonal Ratio The Diagonal Ratio metric (Q DR ) applies only to quad and hex elements and is defined as follows: Q DR max = min [ d1, d2, K, dn ] [ d, d, K, d ] The d i are the diagonals of the element. N is the total number of diagonals for a given element N = 2 for quad elements N = 4 for hex elements. 1 2 N d 3 d 1 d 2 d 1 d 4 d

35 Mesh Quality dge Ratio The dge Ratio quality metric (Q R ) is defined as follows: Q R max = min [ s1, s2, K, sn ] [ s, s, K, s ] The s i are the edge lengths of the element. N is the total number of edges for the element of interest. 1 2 N Tri N = 3 Quad N = 4 Tet N = 6 Pyramid N = 8 Wedge N = 9 Hex N =

36 Mesh Quality quiize kew The quiize kew metric (Q V ) applies only to quad and hex elements and is defined as follows: Q V = eq eq is the area (2D) or volume (3D) of the element of interest. eq is the maximum area (2D) or volume (3D) of an equilateral cell the circumscribing radius of which is identical to that of the mesh element. Actual element Area = 0 < Q V < 1 quilateral element Area = eq Q V =

37 Mesh Quality MidAngle kew The MidAngle kew (Q MA ) applies only to quadrilateral and hexahedral elements. Defined by the cosine of the minimum angle formed between the bisectors of the element edges (quad) or faces (hex). For quad elements: Q MA = cosθ For hex elements: Q MA = max [ cosθ,cosθ, θ ] 1 2 cos 3 θ Bisectors 5-37

38 Mesh Quality tretch The tretch quality metric (Q ) applies only to quadrilateral and hexahedral elements and is defined as follows: d i is the length of diagonal i s j is the length of the element edge j, Q [ s1, s2, K, sm ] [ d, d, K, d ] n and m are the total numbers of diagonals and edges, respectively. Quad elements: n = 2, m = 4, and K = 2; Hex elements: n = 4, m = 12, and K = 3. = 1 K min max 1 2 n By definition, 0 < Q < 1. Q = 0 describes an equilateral element Q = 1 describes a completely degenerate element. 5-38

39 Mesh Quality Taper The Taper quality metric (Q T ) applies only to quadrilateral and hexahedral mesh elements and is defined as follows: For any quadrilateral (or hexahedral) mesh element, it is possible to construct a parallelogram (or parallelepiped) such that the distance between any given corner of the parallelogram (or parallelepiped) and its nearest element corner node is a constant value. As a result, any vector, T, constructed from an element corner node to the nearest corner of the parallelogram (or parallelepiped) possesses a magnitude identical to that of all other such vectors. ach vector T can be resolved into components, T i, that are parallel to the bisectors of the mesh element. Quad elements: two components Hex elements: three components T 2 T T 1 Corner node lement edge Bisectors The Taper quality metric is defined as the normalized maximum of all such components for the element. By definition, 0 < Q T < 1. Q T = 0 describes an equilateral element Q T = 2 describes a degenerate element 5-39

40 Mesh Quality Volume and Warpage Volume The Volume quality metric (Q V ) applies only to 3D elements and represents quality in terms of element volume. Warpage The Warpage (Q W ) applies only to quad elements and is defined as follows: Z is the deviation from a best-fit plane that contains the element a and b are the lengths of the line segments that bisect the edges of the element. By definition, 0 < Q W < 1 Z Q W = min a, [ b] Q W = 0 describes an equilateral element Q W = 1 describes a degenerate element. lement edge Best-fit plane Bisectors 5-40