A variational nodal formulation for multi-dimensional unstructured neutron diffusion problems