Locating the peaks of least-energy solutions to a quasilinear elliptic Neumann problem (Q2381944)

The authors deal with the shape of certain solutions to the following quasilinear elliptic Neumann problem \[ \begin{aligned} \varepsilon^m\Delta_m u- u^{m-1}+ f(u)= 0 &\quad\text{in }\Omega,\\ u> 0 &\quad\text{in }\Omega,\\ {\partial u\over\partial\vec n}= 0 &\quad\text{on }\partial\Omega, \end{aligned} \] where \(m\) \((2\leq m< d)\) and \(\varepsilon> 0\) are constants and \(\Omega\subseteq \mathbb{R}^d\) \((d\geq 3)\) is a smooth bounded domain. The operator \(\Delta_m u= \text{div}(|\nabla u|^{m-2}\nabla u)\) is the \(m\)-Laplacian operator, and \(\vec n\) is the unit outer normal to \(\partial\Omega\). Using variation method, the authors show that as \(\varepsilon\to 0^+\), the global maximum point \(Q^\varepsilon\) of least-energy solutions goes to a point on the boundary \(\partial\Omega\) at the rate of \(O(\varepsilon)\) and this point on the boundary approaches to a point where the mean curvature of \(\partial\Omega\) achieves its maximum. Moreover, they give a complete proof of exponential decay of least-energy solutions.