The interaction between solitary waves and underwater barriers is investigated using our in-house code, entitled VPM (volume-average/point-value multi-moment)–THINC/QQ (THINC method with quadratic surface representation and Gaussian quadrature)-coupled model. The stability and accuracy of the proposed model are validated by comparing the numerical results with those of the well established two-phase flow solver interFoam. All the results indicate that the presented coupled model has the advantage of high fidelity in simulating solitary wave propagation. Subsequently, solitary waves passing over a single underwater barrier are simulated by the present model. Numerical results are compared with experimental results in terms of the free surface elevation, velocity profile, vorticity field, and wave forces. Great agreements are obtained. In the end, the interactions between solitary waves and double underwater barriers are investigated numerically. The results reveal that the reflection coefficient increases first, and then decreases, with the increasing space between the two barriers. For cases with different wave heights, the transmission coefficient decreases monotonically, and the dissipation coefficient is opposed to the transmission coefficient.