We develop and analyse Neumann-Neumann methods for hp finite‐element approximations of scalar elliptic problems on geometrically refined boundary layer meshes in three dimensions. These are meshes that are highly anisotropic where the aspect ratio typically grows exponentially with the polynomial degree. The condition number of our preconditioners is shown to be independent of the aspect ratio of the mesh and of potentially large jumps of the coefficients. In addition, it only grows polylogarithmically with the polynomial degree, as in the case of p approximations on shape‐regular meshes. This work generalizes our previous one on two‐dimensional problems in Toselli & Vasseur (2003a, submitted to Numerische Mathematik, 2003c to appear in Comput. Methods Appl. Mech. Engng.) and the estimates derived here can be employed to prove condition number bounds for certain types of FETI method
