Let $\Omega$ be a bounded domain in ${ℝ}^n$ with smooth boundary $\partial \Omega$. We consider the equation $d^2\Delta u-u+u^{\frac{n-k+2}{n-k-2}}=0\text{in}\Omega$, under zero Neumann boundary conditions, where $\Omega$ is open, smooth and bounded and $d$ is a small positive parameter. We assume that there is a $k$-dimensional closed, embedded minimal submanifold $K$ of $\partial \Omega$, which is non-degenerate, and certain weighted average of sectional curvatures of $\partial \Omega$ is positive along $K$. Then we prove the existence of a sequence $d=d_j\to 0$ and a positive solution $u_d$ such that $d^2{\left|\nabla u_d\right|}^2\rightharpoonup S{\delta }_K\mathrm{as}d\to 0$ in the sense of measures, where ${\delta }_K$...