Packed bed reactors are the most frequently used devices to perform heterogeneously catalyzed reactions on industrial scales. An industrial real-life heterogeneous catalysis is complex process that combines fully three-dimensional mass, momentum and energy transport on several scales. In the present work, we leverage our previously developed CFD solver for non-isothermal heterogeneously catalyzed reactive flow based on the finite volume method and couple it with our in-house DEM-based method for preparation of random packed beds. The resulting framework is verified in the simplified cases against available analytical solutions and correlations and is used to study an industrially-relevant case of ethylene oxychlorination performed in a tubular packed bed comprising CuCl₂-coated catalyst carrying particles. In particular, we compare properties of three different industrially used catalyst carrying particles: Raschig rings, Reformax, and Wagon wheels.