En la relatividad general y en la geometría diferencial, el tensor de Bel-Robinson está definido en la notación de índices abstracta por:

${\displaystyle T_{abcd}=C_{aecf}C_{b}{}^{e}{}_{d}{}^{f}+{\frac {1}{4}}\epsilon _{ae}{}^{hi}\epsilon _{b}{}^{ej}{}_{k}C_{hicf}C_{j}{}^{k}{}_{d}{}^{f}}$

Alternativamente,

${\displaystyle T_{abcd}=C_{aecf}C_{b}{}^{e}{}_{d}{}^{f}-{\frac {3}{2}}g_{a[b}C_{jk]cf}C^{jk}{}_{d}{}^{f}}$

donde ${\displaystyle C_{abcd}}$ es el tensor de Weyl. Fue introducido por Lluís Bel en 1959.^[1]​^[2]​ El tensor de Bel-Robinson se construye a partir del tensor de Weyl de forma análoga a la forma en que el tensor de energía-impulso electromagnético se construye a partir del tensor de campo electromagnético. Al igual que el tensor electromagnético de tensión-energía, el tensor de Bel-Robinson es totalmente simétrico y carece de traza:

${\displaystyle {\begin{aligned}T_{abcd}&=T_{(abcd)}\\T^{a}{}_{acd}&=0\end{aligned}}}$

En la relatividad general, no existe una definición única de la energía local del campo gravitatorio. El tensor de Bel-Robinson es una posible definición de energía local, ya que se puede demostrar que siempre que el tensor de Ricci desaparece (es decir, en el vacío), el tensor de Bel-Robinson no tiene divergencia:

${\displaystyle \nabla ^{a}T_{abcd}=0}$