[step:Verify the rarefaction candidate when $u_L<u_R$]Assume $u_L<u_R$. Define the rarefaction candidate $u: \mathbb{R}\times(0,\infty) \to \mathbb{R}$ by
\begin{align*}
u(x,t):=U(x/t).
\end{align*}
Here $U:\mathbb{R}\to\mathbb{R}$ is the profile in the statement. Since $U$ is bounded and continuous, $u$ is bounded and locally integrable.
For each $(x,t)\in\mathbb{R}\times(0,\infty)$, define the self-similar variable $\xi\in\mathbb{R}$ by
\begin{align*}
\xi:=x/t.
\end{align*} On the regions $\xi<f'(u_L)$ and $\xi>f'(u_R)$, the function $u$ is constant, so $\partial_tu+\partial_xf(u)=0$. On the fan region $f'(u_L)<\xi<f'(u_R)$, we have $f'(U(\xi))=\xi$. Since $U$ is $C^1$ on this open interval, the chain rule gives
\begin{align*}
\partial_t u(x,t)=-\frac{x}{t^2}U'(x/t)
\end{align*}
and
\begin{align*}
\partial_x f(u(x,t))=f'(U(x/t))\frac{1}{t}U'(x/t)=\frac{x}{t^2}U'(x/t).
\end{align*}
Thus $\partial_tu+\partial_xf(u)=0$ pointwise away from the two rays $x=f'(u_L)t$ and $x=f'(u_R)t$. Across these rays the function $u$ is continuous, so no jump measure appears in the [distributional derivative](/page/Distributional%20Derivative). Therefore $u$ is a weak solution.
For the Kruzhkov entropy condition, fix $k\in\mathbb{R}$ and use the same maps $\eta_k:\mathbb{R}\to[0,\infty)$ and $q_k:\mathbb{R}\to\mathbb{R}$ as in the shock step. Define $A_k:\mathbb{R}\to[0,\infty)$ and $B_k:\mathbb{R}\to\mathbb{R}$ by $A_k:=\eta_k\circ U$ and $B_k:=q_k\circ U$. These maps are continuous and piecewise locally Lipschitz. At every point where $U(\xi)\neq k$ and $U$ is differentiable,
\begin{align*}
B_k'(\xi)=\operatorname{sgn}(U(\xi)-k)f'(U(\xi))U'(\xi)
\end{align*}
and
\begin{align*}
A_k'(\xi)=\operatorname{sgn}(U(\xi)-k)U'(\xi).
\end{align*}
Inside the fan, $f'(U(\xi))=\xi$, so $B_k'(\xi)-\xi A_k'(\xi)=0$ for almost every $\xi$ in the fan. Outside the fan both derivatives vanish. The possible point $U(\xi)=k$ creates no singular part because $A_k$ and $B_k$ are continuous and piecewise locally Lipschitz.
Let $\psi\in C_c^\infty(\mathbb{R}\times(0,\infty))$ be a test function. Using the substitution $x=t\xi$, under which $d\mathcal{L}^1(x)=t\,d\mathcal{L}^1(\xi)$ for each fixed $t>0$, we compute
\begin{align*}
\int_0^\infty\int_{\mathbb{R}} \bigl(A_k(x/t)\partial_t\psi(x,t)+B_k(x/t)\partial_x\psi(x,t)\bigr)\,d\mathcal{L}^1(x)\,d\mathcal{L}^1(t)
\end{align*}
as
\begin{align*}
\int_0^\infty\int_{\mathbb{R}} \bigl(tA_k(\xi)\partial_t\psi(t\xi,t)+B_k(\xi)\partial_\xi[\psi(t\xi,t)]\bigr)\,d\mathcal{L}^1(\xi)\,d\mathcal{L}^1(t).
\end{align*}
Integrating by parts in $\xi$ and using $B_k'=\xi A_k'$ in the distributional sense gives
\begin{align*}
\int_0^\infty\int_{\mathbb{R}} A_k(\xi)\partial_t\bigl(t\psi(t\xi,t)\bigr)\,d\mathcal{L}^1(\xi)\,d\mathcal{L}^1(t).
\end{align*}
This last integral is zero by [integration by parts](/theorems/2098) in $t$, because $\psi$ has compact support in $\mathbb{R}\times(0,\infty)$. Hence
\begin{align*}
\partial_t\eta_k(u)+\partial_xq_k(u)=0
\end{align*}
in the sense of distributions, and in particular the Kruzhkov entropy inequality holds.
Define the Riemann initial datum $u_0: \mathbb{R}\to\mathbb{R}$ by setting $u_0(x)=u_L$ when $x<0$ and $u_0(x)=u_R$ when $x>0$; the value at $x=0$ is irrelevant for $L^1_{\mathrm{loc}}(\mathbb{R})$ convergence. The initial trace is this Riemann datum. Indeed, for each $t>0$, the function $u(\cdot,t)$ can differ from $u_0$ only on the convex hull of the three points $0$, $f'(u_L)t$, and $f'(u_R)t$. This interval has $\mathcal{L}^1$-measure at most
\begin{align*}
\bigl(|f'(u_L)|+|f'(u_R)|\bigr)t.
\end{align*}
Since $u$ is bounded between $u_L$ and $u_R$, every compact interval $K\subset\mathbb{R}$ satisfies
\begin{align*}
\lim_{t\downarrow 0}\int_K |u(x,t)-u_0(x)|\,d\mathcal{L}^1(x)=0.
\end{align*}[/step]