[Libmesh-users] Defining hanging nodes in the XDA format

[Salazar De T. M. <sal...@ll...>]

Hello,

I have the mesh attached below. It is the result of having three quadrangular elements in line, refine the left and right ones, then on their leftmost and rightmost children, do another refinement. After this, I flattened out the mesh and removed the leftmost and rightmost layer of elements. The issue is that I’ve got hanging nodes that libmesh does not recognize. Is there any way to find them or tell libMesh which ones they are?

Is there an alternative way to generate a similar mesh? I was not able to remove the leftmost and rightmost children without flattening the mesh first. It seems as if each parent element must have four elements.

Thanks
Miguel

libMesh-1.3.0
13           # number of elements
26           # number of nodes
n/a         # boundary condition specification file
.               # subdomain id specification file
n/a         # processor id specification file
n/a         # p-level specification file
8              # type size
0              # uid size
0              # pid size
8              # sid size
0              # p-level size
0              # eid size
0              # side size
0              # bid size
0              # subdomain id to name map
13           # n_elem at level 0, [ type sid (n0 ... nN-1) ]
5 0 0 1 2 3
5 0 4 0 5 6
5 0 6 5 3 7
5 0 1 8 9 10
5 0 10 9 11 2
5 0 12 4 13 14
5 0 14 13 6 15
5 0 15 6 16 17
5 0 17 16 7 18
5 0 8 19 20 21
5 0 21 20 22 9
5 0 9 22 23 24
5 0 24 23 25 11
5.0000000000000000e+01 0.0000000000000000e+00 0.0000000000000000e+00
1.0000000000000000e+02 0.0000000000000000e+00 0.0000000000000000e+00
1.0000000000000000e+02 5.0000000000000000e+01 0.0000000000000000e+00
5.0000000000000000e+01 5.0000000000000000e+01 0.0000000000000000e+00
2.5000000000000000e+01 0.0000000000000000e+00 0.0000000000000000e+00
5.0000000000000000e+01 2.5000000000000000e+01 0.0000000000000000e+00
2.5000000000000000e+01 2.5000000000000000e+01 0.0000000000000000e+00
2.5000000000000000e+01 5.0000000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+02 0.0000000000000000e+00 0.0000000000000000e+00
1.2500000000000000e+02 2.5000000000000000e+01 0.0000000000000000e+00
1.0000000000000000e+02 2.5000000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+02 5.0000000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+01 0.0000000000000000e+00 0.0000000000000000e+00
2.5000000000000000e+01 1.2500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+01 1.2500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+01 2.5000000000000000e+01 0.0000000000000000e+00
2.5000000000000000e+01 3.7500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+01 3.7500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+01 5.0000000000000000e+01 0.0000000000000000e+00
1.3750000000000000e+02 0.0000000000000000e+00 0.0000000000000000e+00
1.3750000000000000e+02 1.2500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+02 1.2500000000000000e+01 0.0000000000000000e+00
1.3750000000000000e+02 2.5000000000000000e+01 0.0000000000000000e+00
1.3750000000000000e+02 3.7500000000000000e+01 0.0000000000000000e+00
1.2500000000000000e+02 3.7500000000000000e+01 0.0000000000000000e+00
1.3750000000000000e+02 5.0000000000000000e+01 0.0000000000000000e+00
0              # presence of unique ids
0              # sideset id to name map
0              # number of side boundary conditions
0              # nodeset id to name map
0              # number of nodesets
0              # sideset id to name map
0              # number of edge boundary conditions
0              # sideset id to name map
0              # number of shellface boundary conditions


--