Constructing coastal structures against wave action is one of the common and necessary works in port construction and shoreline stabilization. Construction cost and environmental impact are the important factors in the design and construction of the shoreline protective structures. The idea of Geosynthetic Cellular Systems (GCS) has been made considering the mentioned factors. Basic parts of the GCS are Geosynthetics (usually Geotextile) as container, frame for modification of the shape of the con-tainer, and fill material which may be the dredged materials. The stability of this type of breakwater is re-lated to the weight of filler materials. On the other hand, because of the low density of the dredged mate-rials and also saturation pf fill materials below water level, the probability of liquefaction under cyclic wave loading is possible. To assess the probability of liquefaction, physical models of this type of system with different geometrical characteristics (cube, rectangular cube and pyramid) with a variety of filler ma-terials (sand and silt) have been done. In this study, the comparison between the results of the physical model described above and the numerical modeling results of finite difference method are presented. The results show that there is a good agreement between physical and numerical modeling which indicates the ability of the numerical methods in prediction of liquefaction potential of these type of coastal protection structures. Then the effect of the soil parameters, the geometry of the structure and wave characterizations is investigated by numerical parametric studies.