In addition, I think there might be a typo in line 619 ngsPETSc/ngsPETSc/utils/fenicsx.py at main · NGSolve/ngsPETSc, where Numpy() should be NumPy().
For the above uniform refinement in 3D, that typo doesn’t affect the results. But if you want to apply the adaptive refinement and don’t want to refine faces, currently I think the following works:
geoModel.ngmesh.Elements2D().NumPy()['refine'] = 0
Thank you!