C'est bon j'ai trouvé, mais il a fallu que je remonte un peu loin dans mes souvenirs là. 
On utilise la méthode des caractéristiques, i.e. on va chercher des courbes x(t) telles que, la fonction h définie par h(t) = u(t,x(t)) soit constante.
En dérivant, on a
dh/dt = du/dt + dx/dt * du/dx.
On a envie de prendre dx/dt = 1-2u(t,x(t)), ce qui nous donnerait
dh/dt = du/dt + (1-2u)*du/dx = du/dt + d/dx (u(1-u)) = 0.
On se retrouve donc à vouloir résoudre
dx/dt = 1 - 2u(t,x(t))
avec la condition initiale
x(0) = x_0.
A priori on a l'air cons, vu que le x dépend de u, qui dépend de x,... mais en fait non. Il faut se rappeler qu'un tel h soit être constant (et sa valeur ne dépendra que du x_0), ce qui donne
dx/dt = 1 - 2u(0,x(0)) = 1 - 2u_0(x_0).
On peut donc intégrer :
x(t) = x_0 + t*(1 - 2u_0(x_0)) = x_0 + t*v(x_0) pour abréger.
On se retrouve donc avec une autre équation à résoudre :
u(t,x_0 + t*v(x_0)) = u_0(x_0), ou encore
u(t,x) = u_0(x - t*v(x)), où l'on a changé x_0 en x pour alléger les notations.
Finalement, la solution cherchée est donnée par la formule
u(t,x) = u_0(x - t*(1-2u_0(x))).
Voilà !
(Et donc, j'avais bien dit une connerie, Burgers 1D se résout entièrement explicitement avec la même méthode.
)