In this paper, we study discrete harmonic functions on infinite penny graphs. For an infinite penny graph with bounded facial degree, we prove that the volume doubling property and the Poincaré inequality hold, which yields the Harnack inequality for positive harmonic functions. Moreover, we prove that the space of polynomial growth harmonic functions, or ancient solutions of the heat equation, with bounded growth rate has finite dimensional property.