We show a locally uniform bound for global nonnegative solutions of the system $u_t=\Delta u+uv-bu$, $v_t=\Delta v+au$ in $(0,+\infty )\times \Omega$, $u=v=0$ on $(0,+\infty )\times \partial \Omega$, where $a>0$, $b\ge 0$ and $\Omega$ is a bounded domain in ${ℝ}^n$, $n\le 2$. In particular, the trajectories starting on the boundary of the domain of attraction of the zero solution are global and bounded.

Using as a main tool the time-regularizing convolution operator introduced by R. Landes, we obtain regularity results for entropy solutions of a class of parabolic equations with irregular data. The results are obtained in a very general setting and include known previous results.

We study regularity results for solutions $u\in H{W}^{1,p}\left(\Omega \right)$ to the obstacle problem $${\int }_{\Omega }𝒜(x,{\nabla }_{ℍ}u){\nabla }_{ℍ}(v-u)\mathrm{d}x\ge 0\forall v\in {𝒦}_{\psi ,u}\left(\Omega \right)$$ such that $u\ge \psi$ a.e. in $\Omega$, where ${𝒦}_{\psi ,u}\left(\Omega \right)=\{v\in H{W}^{1,p}\left(\Omega \right):v-u\in H{W}_0^{1,p}\left(\Omega \right)v\ge \psi \text{a.e.}\text{in}\Omega \}$, in Heisenberg groups ${ℍ}^n$. In particular, we obtain weak differentiability in the $T$-direction and horizontal estimates of Calderon-Zygmund type, i.e. $$\begin{array}{ccc}\hfill dT\psi \in H{W}_{\mathrm{loc}}^{1,p}\left(\Omega \right)& \Rightarrow Tu\in L_{\mathrm{loc}}^p\left(\Omega \right),{\left|{\nabla }_{ℍ}\psi \right|}^p\in L_{\mathrm{loc}}^q\left(\Omega \right)\hfill & \hfill \Rightarrow |{\nabla }_{ℍ}{u|}^p\in L_{\mathrm{loc}}^q\left(\Omega \right),\end{array}d$$ where $2<p<4$, $q>1$.

Let $\Omega$ be a bounded open subset of ${ℝ}^n$, let $X=(x,t)$ be a point of ${ℝ}^n\times {ℝ}^N$. In the cylinder $Q=\Omega \times (-T,0)$, $T>0$, we deduce the local differentiability result $$u\in L^2(-a,0,H^2(B\left(\sigma \right),{ℝ}^N))\cap H^1(-a,0,L^2(B\left(\sigma \right),{ℝ}^N))$$ for the solutions $u$ of the class $L^q(-T,0,H^{1,q}(\Omega ,{ℝ}^N))\cap C^{0,\lambda }(\overline{Q},{ℝ}^N)$ ($0<\lambda <1$, $N$ integer $\ge 1$) of the nonlinear parabolic system $$-\sum _{i=1}^n{D}_i{a}^i(X,u,Du)+{\displaystyle \frac{\partial u}{\partial t}}=B^0(X,u,Du)$$ with quadratic growth and nonlinearity $q\ge 2$. This result had been obtained making use of the interpolation theory and an imbedding theorem of Gagliardo-Nirenberg type for functions $u$ belonging to $W^{1,q}\cap C^{0,\lambda }$.