Calculating the traction vector field for a surface

You should be able to do something along the lines of: Iterate connected mesh entities in parallel