Mixing resulting from Kelvin-Helmholtz instability is studied by direct numerical simulations for an intermediate Reynolds number in stratified flows where vortex pairing occurs. We investigate how the interaction of several billows changes the mixing properties. Simulations in a long domain show that only considering the interaction of two billows overestimates the amount of mixing and mixing ...