AbstractThis article presents several sufficient conditions for the existence of at least one weak solution and infinitely many weak solutions for the following Neumann problem, originated from a capillary phenomena,
$$\begin{equation*}
\left\{\begin{array}{ll}
-{\rm div}\bigg(\bigg(1+\frac{|\nabla u|^{p(x)}}{\sqrt{1+|\nabla u|^{2p(x)}}}\bigg)
|\nabla u|^{p(x)-2}\nabla u\bigg)+\alpha(x)|u|^{p(x)-2}u\\=\lambda f(x,u) \mbox{in}\,\,\Omega,\\
\frac{\partial u}{\partial \nu}=0\mbox{on}\,\,\partial
\Omega
\end{array}\right.
\end{equation*}$$where $\Omega \subset \mathbb{R}^N$$(N\geq 2)$ is a bounded domain with boundary of class $C^1,$$\nu$ is the outer unit normal to $\partial \Omega,$$\lambda>0$, $\alpha\in L^{\infty}(\Omega),$$f:\Omega\times\mathbb{R}\to\mathbb{R}$ is an $L^1$-Carathéodory function and $p\in C^0(\overline{\Omega})$. Our technical approach is based on variational methods and we use a more precise version of Ricceri’s Variational Principle due to Bonanno and Molica Bisci. Some recent results are extended and improved. Some examples are presented to illustrate the application of our main results.