  In wet liquid foams, slow diffusion of gas through bubble walls changes
bubble pressure, volume and wall curvature. Large bubbles grow at the expenses
of smaller ones. The smaller the bubble, the faster it shrinks. As the number
of bubbles in a given volume decreases in time, the average bubble size
increases: i.e. the foam coarsens. During coarsening, bubbles also move
relative to each other, changing bubble topology and shape, while liquid moves
within the regions separating the bubbles. Analyzing the combined effects of
these mechanisms requires examining a volume with enough bubbles to provide
appropriate statistics throughout coarsening. Using a Cellular Potts model, we
simulate these mechanisms during the evolution of three-dimensional foams with
wetnesses of $\phi=0.00$, $0.05$ and $ 0.20$. We represent the liquid phase as
an ensemble of many small fluid particles, which allows us to monitor liquid
flow in the region between bubbles. The simulations begin with $2 \times 10^5$
bubbles for $\phi = 0.00$ and $1.25 \times 10^5$ bubbles for $\phi = 0.05$ and
$0.20$, allowing us to track the distribution functions for bubble size,
topology and growth rate over two and a half decades of volume change. All
simulations eventually reach a self-similar growth regime, with the
distribution functions time independent and the number of bubbles decreasing
with time as a power law whose exponent depends on the wetness.
