For the structure of the thin electrical double layer (EDL) and the property related to the EDL capacitance, we analyze boundary layer solutions (corresponding to the electrostatic potential) of a non-local elliptic equation which is a steady-state Poisson–Nernst–Planck equation with a singular perturbation parameter related to the small Debye screening length. Theoretically, the boundary layer solutions describe that those ions exactly approach neutrality in the bulk, and the extra charges are accumulated near the charged surface. Hence, the non-neutral phenomenon merely occurs near the charged surface. To investigate such phenomena, we develop new analysis techniques to investigate thin boundary layer structures. A series of fine estimates combining the Pohožaev's identity, the inverse Hölder type estimates and some technical comparison arguments are developed in arbitrary bounded domains. Moreover, we focus on the physical domain being a ball with the simplest geometry and gain a clear picture on the effect of the curvature on the boundary layer solutions. The content involves three contributions. The first one focuses mainly on the boundary concentration phenomena. We show that the net charge density behaves exactly as Dirac delta measures concentrated at boundary points. The second one is devoted to pointwise descriptions with curvature effects for the thin boundary layer. An interesting outcome shows that the significant curvature effect merely occurs in the part of the boundary layer close to the boundary, and this part is extremely thinner than the whole boundary layer. The third contribution gives a connection to the EDL capacitance. We provide a theoretical way to support that the EDL has higher capacitance in a quite thin region near the charged surface, not in the whole EDL. In particular, for the cylindrical electrode, our result has a same analogous measurement as the specific capacitance of the well-known Helmholtz double layer.