Info : Blossom recombination algorithm completed (0.040003 s) Info : Meshing surface 1 (Plane, Frontal Quad) Info : 39 points found in points clouds (2 edges) Info : 5 points found in points clouds (0 edges) Info : Started on Fri Apr 4 13:41:36 2014Įrror : 'Simple_quad.geo', line 23 : Point 3 already existsĮrror : 'Simple_quad.geo', line 24 : Point 5 already existsĮrror : 'Simple_quad.geo', line 25 : Point 6 already existsĮrror : 'Simple_quad.geo', line 26 : Point 7 already existsĮrror : 'Simple_quad.geo', line 27 : Point 8 already existsĮrror : 'Simple_quad.geo', line 28 : Point 9 already existsĮrror : 'Simple_quad.geo', line 29 : Point 10 already existsĮrror : 'Simple_quad.geo', line 30 : Point 11 already exists Info : Running '/usr/bin/gmsh Simple_quad.geo -2' Next we write the Gmsh input file and run the program to generate the mesh: Simple_quad.writegeo( 'Simple_quad',mopt) The mesh will take into account the refinement around the wells but the wells willcorrespond to elements not to nodes. We have found that for quadrilateral it is better not to constraint the mesh on the features, therefore we will let the mesh to ignore the wells. Here we will set as general maximum element size 200 m, element type quad for quadrilateral and element order 2 for quadratic. We initialized the mesh option structure and change the fields to the desired options. Simple_quad=Simple_quad.readshapefile(wells) Simple_quad=Simple_quad.readshapefile(dummyLine) Simple_quad=Simple_quad.readshapefile(dom) Note that the we read first the variable that describes the aquifer outline. The nex step is to create an empty Constructive Solid Geometry object Simple_quad=CSGobj_v2(2,1,10,20,10) % Dim, Npoly, Nline, Npoints, usertolĪfter the object is initialized we read the shapefiles using the method. These features are all what we need to represent the geometry of the modeling domain. Next we define the wells and set the fields for refinement. To refine the mesh along the top right boundary we create a dummy feature dummyLine.Geometry= 'Line' Note that not all the fields that a shapefile has, are necessary (i.e. In the flow example this corresponds to a discharge boundary.įist we define the domain in a structure similar to a matlab shapefile dom.Geometry= 'Polygon' We want also to refine the top right boundary of the mesh. The domain is a rectangle with dimensions 1km x 1km. In addition we will use quadratic quadrilaterals where each element has 9 degrees of freedom. The goal of this example is to show how to create quadrilateral meshes. Which I had to do to make this model work on my version of GMSH.| main | Tutorials | Functions | website | I explicitly setup the meshing algorithm and recombination inside the GEO.Notice slight modifications to your GEO file: I am able to generate the following mesh: With this code: Mesh.Algorithm = 9 // Packing for Parallelograms Is set fixed packing of parallelograms regression in 4.1.1. Advice appreciated.įor this example, I would highly recommend the new (experimental) meshing option in GMSH: 9 – packing for parallelograms.Īccording to the release notes, this option appeared in GMSH at least in 4.1.1:Ĥ.1.2 (January 21, 2019): fixed full-quad subdivision if Mesh.SecondOrderLinear I suspect this maybe due to my relative inexperience with geo/gmsh. Info : Done recombining 2D mesh (0.012205 s) Info : Blossom recombination completed (0.012128 s): 230 quads, 0 triangles, 0 invalid quads, 0 quads with Q < 0.1, avg Q = 0.799983, min Q = 0.502415 Info : Meshing surface 1 (Surface, Frontal) The console log shows: Info : Meshing 1D. geo file into gmsh gui (v 4.3.0) and run mesh 1D then 2D (Frontal-Delaunay option) and finally 2D recombination (Blossom option) commands the resulting surface mesh is not that regular: geo file simplifies the type of intended use case: lc = 0.1 I want to be able to generate a quadrilateral surface mesh that is highly regular (each face has, as far as possible, the same area) and aligned with the surface boundary.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |