  The understanding of dynamical evolutions of interacting photon pulses in
Rydberg atomic ensemble is the prerequisite for realizing quantum devices with
such system. We present an approach that efficiently simulates the dynamical
processes, using a set of local functions we construct to reflect the profiles
of narrowband pulses. For two counter-propagating photon pulses, our approach
predicts the distinct phenomena from the widely concerned Rydberg blockade to
the previously less noticed significant absorption in the anomalous dispersion
regime, which can occur by respectively setting the pulse frequency to the
appropriate values. Our numerical simulations also demonstrate how spatially
extending photon pulses become deformed under realistic non-uniform interaction
over their distributions.
