  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\'e 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.
