We impement and perform a large-scale LES simulation for aerodynamics on a road racing of bicycles. The mesh-refined lattice Boltzmann method with coherent-structure Smagorinsky model is employed to perform a multi-GPU computing. The validation on alone running or group running of 4 cyclists had a good agreement with previous experiments or CFD simulations. The aerodynamics simulation on the group of 72 cyclists is performed as the benchmark of the largest-scale simulation, and it took 4 days using 192 GPUs. The computational resources for the simulations were enough reasonable to run the application studies.