  We develop a numerical method to estimate the average speed of the free
boundary in a Hele-Shaw problem with periodic coefficients in both space and
time. We test the accuracy of the method and present a few examples. We show
numerical evidence of flat parts (facets) on the free boundary in the
homogenization limit.
