Para una función dada $g\in L^2(\Omega)$ se define $\bar g:=\int_\Omega g\ dx$. Al probar la ecuación con la función constante $1$, realizando integración por partes en el espacio se obtiene la ecuación diferencial ordinaria $$ (\bar u)_t = \bar f, \ \bar u(0)=0. $$ Entonces $\|\bar u\|_{L^\infty(0,T)} \le \|\bar f\|_{L^1(0,T)}$.
Al probar la ecuación con $u-\bar u $, integrando por partes en el espacio, e integrando desde $0$ hasta $t$ en el tiempo se llega a $$ \|u(t)-\bar u\|_{L^2}^2 - \| u(0)\|_{L^2}^2 + \|\nabla u\|_{L^2(0,t;L^2)}^2 = \int_0^t \int_\Omega f (u-\bar u) dt. $$ Ahora $$ \int_0^t\int_\Omega f (u-\bar u) dt =\int_0^t\int_\Omega (f-\bar f) (u-\bar u) dt \le \|f-\bar f\|_{L^2(0,t;L^2)} \|u-\bar u\|_{L^2(0,t;L^2)} , $$ el último factor es menor que $c\|\nabla u\|_{L^2(0,t;L^2)} $ por la desigualdad de Poincaré ($c$ depende solo de $\Omega$). Y esta expresión puede ser compensada por el lado izquierdo.
De manera similar, al probar la ecuación con $-\Delta u$, integrando por partes en el espacio en el primer término, e integrando desde $0$ hasta $t$ en el tiempo se llega a $$ \|\nabla u(t)\|_{L^2}^2 - \|\nabla u(0)\|_{L^2}^2 + \|\Delta u\|_{L^2(0,t;L^2)}^2 = \int_0^t f (-\Delta u) dt. $$ Combinando todas estas estimaciones, obtenemos: $$ \|\bar u\|_{L^\infty(0,T)}^2 + \|u-\bar u\|_{L^\infty(0,T;L^2)}^2 + \|\nabla u\|_{L^2(0,T;L^2)}^2 + \|\Delta u\|_{L^2(0,T;L^2)}^2 \le c\left( \|\bar f\|_{L^1(0,T)}^2 + \|f-\bar f\|_{L^2(0,T;L^2)}^2 + \|f\|_{L^2(0,T;L^2)}^2 \right), $$ donde $c$ es independiente de $T$. Al estimar estas normas de $L^1$ y $L^\infty$ contra normas de $L^2$, obtendríamos factores adicionales dependiendo de $T$. Sin embargo, al no haber utilizado la desigualdad de Gronwall, las constantes no son exponenciales en $T$.