We have successfully developed a least-squares finite element method (LSFEM) for unsteady-state, three-dimensional convection dominated flows. The serial algorithm was used for simulations on a single CPU. In this work, our goal is to use a domain decomposition approach to parallelize the serial algorithm without much modifications of the serial code. We illustrate the efficiency of the parallel algorithm by flow simulations for more than 25 million unknowns.