We are concerned with non-constant positive radial solutions of the system $$ \left\{ \begin{ aligned} S_k(D^2 u)&=|\nabla u|^{ m} v^{ p}&&\quad\mbox{ in }\Omega,\\ S_k(D^2 v)&=|\nabla u|^{ q} v^{ s} &&\quad\mbox{ in }\Omega, \end{ aligned} \right. $$ where $S_k(D^2u)$ is the $k$-Hessian operator of $u\in C^2(\Omega)$ ($1\leq k\leq N$) and $\Omega\subset\mathbb{ R}^N$ $(N\geq 2)$ is either a ball or the whole space. The exponents satisfy $q>0$, $m,s\geq 0$, $p\geq s\geq 0$ and $(k-m)(k-s)\neq pq$. In the case where $\Omega$ is a ball, we classify all the positive radial solutions according to their behavior at the boundary. Further, we consider the case $\Omega=\mathbb{ R}^N$ and find that the above system admits non-constant positive radial solutions if and only if $0\leq m<k$ and $pq < (k-m)(k-s)$. Using arguments from three component cooperative and irreducible dynamical systems we deduce the behavior at infinity of such solutions.
Categories: math.AP