[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?