Bragg reflection of nonlinear waves is simulated by a 2D fully nonlinear numerical wave tank (NWT). The developed NWT was based on the Boundary Element Method (BEM) with potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme with Runge-Kutta 4th-order time integration. A spatial variation of wave elevations and their Fourier amplitudes of each component are compared to investigate the effect of sea bottom ripples and their relative heights. The incident waves over an undulated sea bottom are partially reflected and changed to partial standing waves due to Bragg reflection. The present results are verified with linear calculations and experimental data. It is found that the 1st-order wave component is mainly affected by Bragg reflection and its spatial modulation is significant in front of the bottom ripples.