Accuracy estimation of wave loading on cylinders in a pile group under different impact scenarios is essential for both the structural safety and cost of coastal and offshore structures. Differing from the interaction of waves with a single cylinder, less attention has been paid to pile groups under different arrangements. Numerical simulations of interactions between plunging breaking waves and pile group in finite water depth are performed using the two-phase flow model in REEF3D, an open-source computational fluid dynamics program to investigate the wave loads and flow kinematics characteristics. The Reynolds-averaged Navier–Stokes equation with the two equation k − ω turbulence model is adopted to resolve the numerical wave tank. The model is validated by comparing the numerical wave forces and free surface elevation with measurements from experiments. The computational results show fairly good agreement with experimental data. Four cases are simulated with different relative distances, numbers of cylinders, and arrangements. Results show that the wave forces on cylinders in the pile group are effected by the relative distance between cylinders. The staggered arrangement has a significant influence on the wave forces on the first and second cylinder. The interaction inside a pile group mostly happens between the neighboring cylinders. These interactions are also effected by the relative distance and the numbers of the neighboring cylinders.