Your PDE is quasilinear, meaning it migh not have a classic solution for all time $t$. That said, we know that the quasilinear equation
$$
a\big(x,t,u(x,t)\big) u_x(x,t) + b\big(x,t,u(x,t)\big)u_t(x,t) = c\big(x,t,u(x,t)\big)
$$
where $a,\,b,\,c \in C^1$ with data $\mathcal{C}(\xi) = \big(x(\xi), t(\xi), u(\xi)\big) \in C^1$ and with
$$
\begin{vmatrix} \frac{dx}{d\xi} & a \\ \frac{dt}{d\xi} & b\end{vmatrix} \neq 0
$$
has a unique solution near $\mathcal{C}$ given by
\begin{align}
\frac{d x}{d \eta} &= a & x\big|_{\eta = 0}&= x(\xi)\\
\frac{d t}{d \eta} &= b & t\big|_{\eta = 0}&= t(\xi)\\
\frac{d u}{d \eta} &= c & u\big|_{\eta = 0}&= u(\xi)\\
\end{align}
In your case, $\mathcal{C}(\xi) = \big(\xi,0,f(\xi)\big)$. Near $\eta \sim 0$
$$
\begin{vmatrix} \frac{dx}{d\xi} & a \\ \frac{dt}{d\xi} & b\end{vmatrix} = \begin{vmatrix} 1 & u \\ 0 & 1\end{vmatrix} = 1
$$
and the solution is unique near the initial condition. Now,
\begin{align}
\frac{d x}{d \eta} &= u & x\big|_{\eta = 0}&= \xi\\
\frac{d t}{d \eta} &= 1 & t\big|_{\eta = 0}&= 0\\
\frac{d u}{d \eta} &= x & u\big|_{\eta = 0}&= f(\xi)\\
\end{align}
has as solution
$$
t = \eta, \quad u = f(\xi)\cosh \eta + \xi \sinh \eta, \quad x =f(\xi) \sinh \eta + \xi \cosh \eta.
$$
and the characteristics are given by the equation
$$
t= \text{arctanh}\left(\frac{x - \xi}{f(\xi)}\right)
$$
(you have to be carefull while inverting $\tanh t$).
Depending on the behavior of $f$, if the characteristics can meet. If the do, there will be a shock and the classical solution will cease to exist. We can see this by studying the transformation
$$
(x,t) \longrightarrow (\xi,\eta).
$$
The change of variables will be invertible iif
$$
\begin{vmatrix} \partial_\xi x & \partial_\eta x \\ \partial_\xi t & \partial_\eta t \end{vmatrix} = f'(\xi) \sinh\eta + \cosh \eta \neq 0
$$
meaning there is no solution when
$$
f'(\xi) \sinh\eta + \cosh \eta = 0
$$
or, inverting the transformation, the solution will develop a shock at time
$$
t = -\text{arctanh}\left(f'(\xi)\right),
$$
and the profile of $f(x)$ is crucial in the shock development. This is why, if you want to understand how this works, you have to study the Burger's equation, as suggested by user88595.
In the region where no shock has developed (might be all the domain), the solution is given by
$$
u(x,t) = f\big(\xi(x,t)\big)\cosh t + \xi(x,t) \sinh t
$$
where $\xi(x,t)$ is determined by inverting
$$
x = f(\xi) \sinh t + \xi \cosh t,
$$
which can be written in implicit form as
$$
f\big(x \cosh t - u(x,t) \sinh t\big) = u(x,t) \cosh t - x \sinh t
$$