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 are an order of magnitude higher than the current ion colliders. 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 collision energy dependence of collective flow will be also presented. These results will help the understanding of UHECR behavior and can be tested at China's large high altitude air shower observatory (LHAASO) experiments.