Hello,<br><br>I'm trying to create a second order quadrangular mesh around a circular geometry, using O-Grids and Transfinite meshes. My FEM code requires that the intermediate nodes are placed at the midpoint of the element edges. Thus, I use the option Mesh.SecondOrderLinear = 1 in my .geo file. That avoids that GMSH move the extra nodes to the circle. Again, I want straight-edges elements. Until here, no issues.<br>

<br>After the mesh is created, I try to refine it using the "refine by splitting" option on the GUI. Even if I do that with a 1st order mesh, the added nodes are placed at the midpoint of the edges of the previous mesh. According to the log file, GMSH creates a 2nd order mesh as part of the refinement process. Thus, it uses the Mesh.SecondOrderLinear = 1 option, and that causes that the refined mesh is not adapted to the circular geometry. Is there any way to turn off this option off for the "refine by splitting" step, but leave it active for the generation of the actual 2nd order mesh?<br>

<br>For me, the expected behaviour or any refinement should be the increase of the mesh density, and improve the definition of the geometrical features of the domain. With the Mesh.SecondOrderLinear active, the second objective is not achieved.<br>

<br>I'm including the .geo and the logfile<br><br>Thanks in advance for any help with this issue,<br><br>=======testestest.geo==========<br>N = 5;<br>Point(1) = {0, 0, 0};<br>Point(2) = {1, 0, 0};<br>Point(3) = {2, 0, 0};<br>

Point(4) = {0, 1, 0};<br>Point(5) = {0, 2, 0};<br>Circle(1) = {2, 1, 4};<br>Transfinite Line {1} = N;<br>Line(2) = {2, 3};<br>Transfinite Line {2} = N;<br>Line(3) = {3, 5};<br>Transfinite Line {3} = N;<br>Line(4) = {5, 4};<br>

Transfinite Line {4} = N;<br><br>Line Loop(1) = {3, 4, -1, 2};<br>Plane Surface(1) = {1};<br>Transfinite Surface {1} = {};<br>Physical Surface(1) = {1};<br><br>Geometry.Tolerance = 1.0e-16;<br><br>Mesh.Smoothing = 0;<br>
Mesh.SmoothNormals = 1;<br>
Mesh.Algorithm = 6;<br>Mesh.Points = 1;<br>Mesh.PointSize = 8;<br>Mesh.SecondOrderLinear = 1;<br>Mesh.RecombineAll = 1;<br><br>Mesh.ElementOrder = 1;<br>Mesh.Format = 30;<br>Mesh.SaveElementTagType = 2;<br><br>=====testestest.log==========<br>

Info    : Reading 'testestest'...<br>Info    : Done reading 'testestest'<br>Info    : Meshing 1D...<br>Info    : Meshing curve 1 (Circle)<br>Info    : Meshing curve 2 (Line)<br>Info    : Meshing curve 3 (Line)<br>

Info    : Meshing curve 4 (Line)<br>Info    : Done meshing 1D (0.004 s)<br>Info    : Meshing 2D...<br>Info    : Meshing surface 1 (transfinite)<br>Info    : Done meshing 2D (0 s)<br>Info    : 26 vertices 37 elements<br>Info    : Refining mesh...<br>

Info    : Meshing order 2...<br>Info    : Meshing curves order 2...<br>Warning : We don't have a valid parameter on curve 3760-3764<br>Warning : We don't have a valid parameter on curve 3764-3765<br>Warning : We don't have a valid parameter on curve 3765-3766<br>

Warning : We don't have a valid parameter on curve 3766-3762<br>Warning : We don't have a valid parameter on curve 3760-3767<br>Warning : We don't have a valid parameter on curve 3767-3768<br>Warning : We don't have a valid parameter on curve 3768-3769<br>

Warning : We don't have a valid parameter on curve 3769-3761<br>Warning : We don't have a valid parameter on curve 3761-3770<br>Warning : We don't have a valid parameter on curve 3770-3771<br>Warning : We don't have a valid parameter on curve 3771-3772<br>

Warning : We don't have a valid parameter on curve 3772-3763<br>Warning : We don't have a valid parameter on curve 3763-3773<br>Warning : We don't have a valid parameter on curve 3773-3774<br>Warning : We don't have a valid parameter on curve 3774-3775<br>

Warning : We don't have a valid parameter on curve 3775-3762<br>Info    : Meshing surfaces order 2...<br>Info    : Meshing volumes order 2...<br>Info    : Done meshing order 2 (0.001 s)<br>Info    : Done refining mesh (0.002 s)<br>

Info    : Meshing order 2...<br>Info    : Meshing curves order 2...<br>Warning : We don't have a valid parameter on curve 3760-3785<br>Warning : We don't have a valid parameter on curve 3785-3764<br>Warning : We don't have a valid parameter on curve 3764-3786<br>

Warning : We don't have a valid parameter on curve 3786-3765<br>Warning : We don't have a valid parameter on curve 3765-3787<br>Warning : We don't have a valid parameter on curve 3787-3766<br>Warning : We don't have a valid parameter on curve 3766-3788<br>

Warning : We don't have a valid parameter on curve 3788-3762<br>Warning : We don't have a valid parameter on curve 3760-3789<br>Warning : We don't have a valid parameter on curve 3789-3767<br>Warning : We don't have a valid parameter on curve 3767-3790<br>

Warning : We don't have a valid parameter on curve 3790-3768<br>Warning : We don't have a valid parameter on curve 3768-3791<br>Warning : We don't have a valid parameter on curve 3791-3769<br>Warning : We don't have a valid parameter on curve 3769-3792<br>

Warning : We don't have a valid parameter on curve 3792-3761<br>Warning : We don't have a valid parameter on curve 3761-3793<br>Warning : We don't have a valid parameter on curve 3793-3770<br>Warning : We don't have a valid parameter on curve 3770-3794<br>

Warning : We don't have a valid parameter on curve 3794-3771<br>Warning : We don't have a valid parameter on curve 3771-3795<br>Warning : We don't have a valid parameter on curve 3795-3772<br>Warning : We don't have a valid parameter on curve 3772-3796<br>

Warning : We don't have a valid parameter on curve 3796-3763<br>Warning : We don't have a valid parameter on curve 3763-3797<br>Warning : We don't have a valid parameter on curve 3797-3773<br>Warning : We don't have a valid parameter on curve 3773-3798<br>

Warning : We don't have a valid parameter on curve 3798-3774<br>Warning : We don't have a valid parameter on curve 3774-3799<br>Warning : We don't have a valid parameter on curve 3799-3775<br>Warning : We don't have a valid parameter on curve 3775-3800<br>

Warning : We don't have a valid parameter on curve 3800-3762<br>Info    : Meshing surfaces order 2...<br>Info    : Meshing volumes order 2...<br>Info    : Done meshing order 2 (0.004 s)<br>Info    : Writing 'testestest.log'...<br>

<br> <br clear="all"><font face="'courier new', monospace">-------------------------<br><b>Félix Salazar</b><br><a href="mailto:felix.a.salazar@gmail.com" target="_blank">felix.a.salazar@gmail.com</a></font><br><div>

<span style="font-family:'courier new',monospace">-------------------------<br></span></div><br>