Dans la première partie de cet article, nous établissons des estimées locales de gradient pour les fonctions $p$-harmoniques à l'intérieur et au bord, sur les variétés riemanniennes générales. Grâce à ces estimations et suivant une idée récente de R. Moser, nous obtenons un théorème d'existence de solutions faibles au sens de la formulation d'ensemble de niveau pour le flot $1/H$ (inverse de la courbure moyenne) des hypersurfaces dans les variétés ambiantes ayant la propriété de la croissance optimale du volume. Dans la deuxième partie, nous considérons deux types d'équations paraboliques pour les fonctions $p$-harmoniques et nous établissons une estimation optimale du type de Li-Yau pour les solutions positives pour ces équations sur les variétés à courbure de Ricci non-négative. Nous montrons aussi une formule de monotonie des entropies associées à ces équations. Cette formule généralise un résultat antérieur du deuxième auteur pour l'équation de la chaleur linéaire. Comme application, nous montrons que toute variété riemannienne complète à courbure de Ricci positive ou nulle et admettant une inégalité logarithmique $L^p$ optimale est isométrique à l'espace euclidien.