  Primordial black holes (PBHs) may form before cosmological first-order phase
transitions, leading to inevitable collisions between PBHs and bubble walls. In
this Letter, we have simulated for the first time the co-evolution of an
expanding scalar wall passing through a black hole with full numerical
relativity. This black hole-bubble wall collision yields multiple far-reaching
phenomena including the PBH mass growth, gravitational wave radiations, and
momentum recoil that endows PBHs with additional velocities, approximately
doubling the formation rate for PBH binaries and hence strengthening the
observational constraints on the PBH abundances.
