Solutions for a variable exponent Neumann boundary value problems with Hardy critical exponent (Q2340806)

In this paper, the authors consider the following nonlinear system involving variable exponents and Neumann boundary conditions with Hardy critical exponent \[ \begin{cases} -\operatorname{div}(|\nabla u|^{p(x)-2}\nabla u)+a(x)|u|^{p(x)-2}u = F_u(x,u,v) \, &\text{in } \;\; \Omega,\\ -\operatorname{div}(|\nabla v|^{q(x)-2}\nabla v)+b(x)|v|^{q(x)-2}v = F_v(x,u,v) \, &\text{in } \;\; \Omega,\\ \frac{\partial u}{\partial \nu}=\frac{\partial v}{\partial \nu}=0 \, &\text{on }\; \partial\Omega, \end{cases}\tag{1} \] where \(\Omega \subset \mathbb R^N\) is a bounded domain and \(\partial\Omega\) possesses the cone property. Under suitable conditions on \(a, b, p\) and \(q\), further if \(F\) satisfies some sub- or super-\((p(x),q(x))\) growth conditions, the authors prove the existence of solutions for problem (1).