  We present long-time simulations of expanding ultracold neutral plasmas,
including a full treatment of the strongly coupled ion dynamics. Thereby, the
relaxation dynamics of the expanding laser-cooled plasma is studied, taking
into account elastic as well as inelastic collisions. It is demonstrated that,
depending on the initial conditions, the ionic component of the plasma may
exhibit short-range order or even a superimposed long-range order resulting in
concentric ion shells. In contrast to ionic plasmas confined in traps, the
shell structures are built up from the center of the plasma cloud rather than
from the periphery.
