In heavy ion collisions, the main goal is to create the quark-gluon plasma (QGP) and then study its properties in order to understand quantum chromodynamics at extreme conditions. Collective flow serves as an important probe to study the production and characterize the property of the QGP. In ultra high energy cosmic rays (UHECR), the collision energies can be easily achieved as the current collider experiment. It is naturally to believe the QGP to be created in UHECR collisions. In this work, collective flow are analytically studied within CORSIKA model, with EPOS-LHC for high energy hadronic interaction. The results suggest the collective flow still holds for all charged particles, further study indicates the collective flow of muons are much larger than inclusive charged particles. It seems promising that we can measure the collectivity of muons to search for QGP in UHECR. It can be tested at China’s large high altitude air shower observatory (LHAASO) experiments.