How to compute the solution gradient on the boundary

Hi there,
I need to compute what is generally referred to as du/dn. Using some math I can get away with it by computing grad(u) \dot n . However, i can’t find any operator/API computing either the former or the latter.
Am I missing something. Ideally, grad(u) should be readily available if one uses shape functions.
Am I missing something?