We study existence and approximation of non-negative solutions of partial differential equations of the type
$$\partial_t u - \div (A(\nabla (f(u))+u\nabla V )) = 0 \qquad \mbox{in } (0,+\infty )\times \mathbb{R}^n,\qquad\qquad (0.1)$$
where A is a symmetric matrix-valued function of the spatial variable satisfying a uniform ellipticity condition,
$f:[0,+\infty) \rightarrow[0,+\infty)$ is a suitable non decreasing function, $V:\mathbb{R}^n \rightarrow\mathbb{R}$ is a convex function.
Introducing the energy functional $\phi(u)=\int_{\mathbb{R}^n} F(u(x))\,{\rm d}x+\int_{\mathbb{R}^n}V(x)u(x)\,{\rm d}x$,
where F is a convex function linked to f by $f(u) = uF'(u)-F(u)$,
we show that u is the “gradient flow” of ϕ with respect to the
2-Wasserstein distance between probability measures on
the space $\mathbb{R}^n$, endowed with the Riemannian distance induced by $A^{-1}.$
In the case of uniform convexity of V, long time asymptotic behaviour and decay rate to the stationary state
for solutions of equation (0.1) are studied.
A contraction property in Wasserstein distance for solutions of equation (0.1)
is also studied in a particular case.