  We study theoretically and numerically the velocity dynamics of fully chaotic
time-dependent shape-preserving billiards. The average velocity of an ensemble
of initial conditions generally asymptotically follows the power law $v =
n^{\beta}$ with respect to the number of collisions n. If the shape of the
chaotic time-dependent billiard is not preserved it is well known that the
acceleration exponent is $\beta = 1/2$. We show, on the other hand, that if the
shape of the billiard is preserved then the only possible values of $\beta$ are
$1/4$, $1/6$ and in a special case $0$. We show that the value of $\beta$ is
determined solely by whether the angular momentum of the billiard is preserved
or not. Our theory is centered around the detailed study of the energy
fluctuations in the adiabatic limit. We show that three quantities, two scalars
and one tensor, completely determine the energy fluctuations of the billiard
for arbitrary time-dependent shape-preserving transformations. Finally we
provide several interesting numerical examples all in a perfect agreement with
the theory.
