Simulation of first order confinement transition in relativistic heavy ion collision

We have carried out numerical simulation of first order phase transition in 2+1 dimensions to study the formation and evolution of Z(3) domain walls in relativistic Heavy Ion Collision using the effective potential proposed by Pisarski for QCD where Polyakov loop is the order parameter of the weak first order phase transition. Bubbles of the QGP phase are randomly nucleated on the lattice, which grow and coalesce. The spontaneous breaking of Z(3) symmetry in QGP phase gives rise to domain walls and topological strings. We discuss P _T enhancement due to reflection of quarks from the collapsing domain walls. We also discuss enhancement of doubly strange and triply strange hadrons due to larger concentration of s quarks inside collapsing wall. The decay of the domain walls when temperature drops below T _c results in the fluctuations of energy density.