[Math] Hexagon packing in a circle

Question

Suppose I want to pack hexagons in a circle, as on the drawing below (red indicates "packed" hexagons). I am wondering what is known about this problem. Specifically, I am interested in an approximation to how many fit (given the radius of the circle and the length of the side of the hexagon) and the error in such approximation. Any ideas?

Answer 1

Let us try to model the grid:

We have six base directions $$ u_k = (x_k, y_k) = d (\cos k \pi/3, \sin k \pi/3) \quad (k\in \{0,\ldots,5\}) $$ where $d$ is the incircle diameter of a hexagon cell. Starting from the origin $(0,0)$ each cell can be reached by walking along one of the six directions, e.g. $$ u = u_2 + 2u_5 - u_1 + 7u_2 + 2 u3 + \cdots $$ This shrinks down to the integer multiples of two directions, because $u_3 =-u_0$ etc and $u_2 = u_1-u_0$: $$ u = c_0 u_0 + c_1 u_1 $$ so we can use $(c_0,c_1)$ as coordinates for each cell.

The vertices of the center cell are $$ v_k = D (\cos (2k+1) \pi/6, \sin (2k+1) \pi/6) \quad (k\in \{0,\ldots,5\}) $$ where $D$ is the diameter of the circumscribed circle. It relates to $d$ via: $$ D = \frac{2}{\sqrt{3}} d $$

We can describe every vertex via $(c_0,c_1,k)$, which is not injective.

Each vertex within the circle fulfills $$ \lVert u + v_k \rVert \le R $$

Here is an image where the red circles (except for the one around the center) have radius $r_c = (2k + 1) r$.

Here the hexagons were coloured the same, when they have the same circumscribed circle (these are the easy to calculate cases):

The number of black circles between red circles increases with increasing radius. This is because there are more discrete possibilities to position the hexagons.

So it is possibly a good idea to let a machine do the counting, it seems to get only more complicated, with the limit of continous radius possibility for large $r$.

Please take the above graphics with a grain of salt: The solid graph is $f(r) = \frac{A(r)}{A_H} = \frac{\pi r^2}{2 \sqrt{3}}$ which is the area of the circumscribed circle expressed in hexagon cells. Due to the gaps at the rim, the number of hexagons (dots) is less than this ideal number. However the dots with the hexagon count for that radius seem too low, I need to hunt the bug and provide a version with corrected numbers.