[guided]Assume that $f$ is convex on $U$, and fix two points $x,y \in U$. The desired inequality compares $f(y)$ with the affine approximation to $f$ at $x$, so we examine $f$ only along the line segment beginning at $x$ and ending at $y$. Since $U$ is convex, every point of this segment lies in $U$, so the map
\begin{align*}
\gamma_{x,y}: [0,1] \to U, \qquad t \mapsto x+t(y-x)
\end{align*}
is well-defined. Define
\begin{align*}
\varphi_{x,y}: [0,1] \to \mathbb{R}, \qquad t \mapsto f(\gamma_{x,y}(t)).
\end{align*}
The point $\gamma_{x,y}(t)$ is the convex combination $(1-t)x+ty$. Therefore convexity of $f$ gives, for each $t \in (0,1]$,
\begin{align*}
f(x+t(y-x)) \leq (1-t)f(x)+tf(y).
\end{align*}
In terms of $\varphi_{x,y}$, this is
\begin{align*}
\varphi_{x,y}(t) \leq (1-t)\varphi_{x,y}(0)+t\varphi_{x,y}(1).
\end{align*}
We now isolate a difference quotient at the point $t=0$. Subtracting $\varphi_{x,y}(0)$ from both sides gives
\begin{align*}
\varphi_{x,y}(t)-\varphi_{x,y}(0) \leq t(\varphi_{x,y}(1)-\varphi_{x,y}(0)).
\end{align*}
Since $t>0$, division by $t$ preserves the inequality:
\begin{align*}
\frac{\varphi_{x,y}(t)-\varphi_{x,y}(0)}{t} \leq \varphi_{x,y}(1)-\varphi_{x,y}(0).
\end{align*}
The left-hand side is the one-sided difference quotient of $f$ at $x$ in the direction $y-x$. Because $f$ is differentiable at $x$, its first-order expansion at $x$ implies
\begin{align*}
\lim_{t\downarrow 0}\frac{f(x+t(y-x))-f(x)}{t}=\nabla f(x)\cdot(y-x).
\end{align*}
The right-hand side is independent of $t$, so taking $t\downarrow 0$ gives
\begin{align*}
\nabla f(x)\cdot(y-x) \leq f(y)-f(x).
\end{align*}
Rearranging gives the supporting inequality
\begin{align*}
f(y) \geq f(x)+\nabla f(x)\cdot(y-x).
\end{align*}[/guided]