Boundary conditions get applied based on the DOF number, not the geometric entity ID.
You may want to try something like
dof_center = dolfinx.fem.locate_dofs_geometrical(V, center)
bc_center = dolfinx.fem.dirichletbc(u_center, dof_center)
although I haven’t tested it with your example.