A generalization of orthogonal arrays, namely cluster orthogonal arrays (CLOA), is introduced and some properties and construction methods are studied. The universal optimality of the fractional factorial designs represented by cluster orthogonal arrays is proved.