Large eddy simulation (LES) with 1.12×108 meshes was performed on three-dimensional jet in cross-flow(JICF) using lattice Boltzmann method (LBM) and multiple Graphic Processing Units (Multi-GPUs). Reynolds number based on free-stream velocity and jet diameter is 4 000, streamwise inclination angle of jet is α=30°, and jet-to-cross-flow velocity ratio is set as 0.5. Validity and feasibility of LBM-LES on simulating jet in cross-flow were verified by reasonable qualitative and quantitative results. Fine coherent structures were captured by large-scaled simulation which benefits study on mixing mechanism of jet-in-cross-flow(JICF). Moreover, 6 K20M GPUs were adopted in simulation and it took 15 402 seconds for LBM-GPU solver to simulate 71 680 LBM steps resulting in calculated performance of 521. 24 MLUPS.