Nodal solutions for Neumann systems with gradient dependence
We consider the following convective Neumann systems:\begin{equation*}\left(\mathrm{S}\right)\qquad\left\{\begin{array}{ll}-Δ_{p_1}u_1+\frac{|\nabla u_1|^{p_1}}{u_1+δ_1}=f_1(x,u_1,u_2,\nabla u_1,\nabla u_2) & \text{in}\;Ω,\\ -Δ_{p_2}u_2+\frac{|\nabla u_2|^{p_2}}{u_2+δ_2}=f_2(x,u_1,u_2,\nabla u_1,\nabla u_2)&\text{in}\;Ω, \\ |\nabla u_1|^{p_1-2}\frac{\partial u_1}{\partial η}=0=|\nabla u_2|^{p_2-2}\frac{\partial u_2}{\partial η}&\text{on}\;\partialΩ,\end{array}\right.\end{equation*}where is a bounded domain in () with a smooth boundary , are small parameters, is the outward unit vector normal to are Carathéodory functions that satisfy certain growth conditions, and ( for ) are the -Laplace operators ,for every In order to prove the existence of solutions to such systems, we use a sub-supersolution method. We also obtain nodal solutions by constructing appropriate sub-solution and super-solution pairs. To the best of our knowledge, such systems have not been studied yet.