[emab]

I'm using Gmsh to generate a 2D Strutured Mesh with a boundary layer size field. When I try to run the simulation (on SU2) it returns "Zero-area CV face found for point 340."
I then went on to open the SU2 mesh file and searched for the location of point 340. Then, I found that if I view only the mesh edges in Gmsh (without geometry lines on top) I get some weird lines highlighted, one of them incluing point 340.
Screenshot 2024-11-23 221554.png
This probably has to do with the way I'm generating the boundary layer size field. I'm using the python API to paramatrically generate meshes, here's the section where I generate the field, curves and surfaces:

    f = gmsh.model.mesh.field.add('BoundaryLayer')
    gmsh.model.mesh.field.setNumbers(f, 'CurvesList', [c_tags[i] for i in 'tip nc_body rocket_body_long rocket_body_corner rocket_body_rad'.split()])
    gmsh.model.mesh.field.setNumber(f, 'Size', mesh.BL_first_cell)
    gmsh.model.mesh.field.setNumbers(f, 'NodesList', body_points)
    gmsh.model.mesh.field.setNumber(f, 'SizeFar', mesh.nc_tip)
    gmsh.model.mesh.field.setNumber(f, 'Quads', 0)
    gmsh.model.mesh.field.setNumber(f, 'Thickness', mesh.BL_thickness)
    gmsh.model.mesh.field.setAsBoundaryLayer(f)
    geom.synchronize()
    for l in curve_loops: geom.addCurveLoop(curve_loop_tags(curve_loops[l]), reorient=True)
    for i in range(1,len(curve_loops)+1): geom.addPlaneSurface([i])
    geom.mesh.setTransfiniteSurface(1, "Right", [p_tags[i] for i in 'tip_start rocket_body_0 offset_tip_start offset_point_0'.split()])
    geom.mesh.setTransfiniteSurface(2)
    geom.mesh.setTransfiniteSurface(3)
    geom.mesh.setTransfiniteSurface(4, "Right", [p_tags[i] for i in 'rocket_body_1 rocket_body_3 offset_point_1 offset_point_3'.split()])

I think it probably has to do with setting the NodesList to every point that the curves intersect, but thats the only way I have been able to make the boundary layer show up on top of the transfinite surfaces.
Any tips?

[bigfootedrockmidget]

Maybe start by generating a mesh with only 2 layers in the boundary so you can actually see what is going on there. It is now difficult to visually inspect the mesh correctness. Maybe the lines are so close together that they are overlapping? Difficult to say.