방파제 형식에 따른 반사율 변화가 해저지형에 미치는 영향을 살펴보기 위한 수치 모의를 수행하였다. 수치 모형은 OpenFoam 기반 tool box인 OlaFlow와 물리기반 지형모형[SeoulFoam]으로 구성하였으며, 이 과정에서 침·퇴적으로 인해 변형을 겪는 해저지형과 내습하는 파랑 간의 상호작용은 Dynamic Mesh를 활용하여 기술하였다. 다양한 반사 특성을 보이는 사석 경사제, 직립제, 곡면 슬릿 케이슨 방파제는 서로 다른 정상파를 결과하였으며, 이는 해저지형에 상당한 영향을 미치는 것을 확인하였다. 이러한 결과는 저면 유속이 상대적으로 큰 정상파 마디[node]에서 연행된 모래가 경계층 drift에 의해 배[antinode]로 이송된다는 Nielsen(1993)의 연구 결과와도 일치한다. 이렇게재배치되는 모래로 정상파의 배[antinode]에는 sand wave의 마루, 마디[antinode] 인근에는 sand wave의 곡이 형성되었으며, sand wave 진폭은 반사계수가 우월한 곡면 슬릿 케이슨에서 가장 크게 관측되었다. 이러한 현상은 반사계수가 큰 경우 마디에서의 저면 유속 증가로 상대적으로 많은 모래가 연행되어 발생하는 것으로 판단된다.
Numerical simulations were implemented to look into the modified seabed topography due to the presence of breakwaters of varying reflection characteristics. The numerical model was composed of OlaFlow, an OpenFoam-based tool box, and a physics-based morphology model [Seoul Foam]. In doing so, the interaction between the seabed, which undergoes deformation due to siltation and scouring, and the incoming waves was described using Dynamic Mesh. The rubble-mound, vertical, and curved slit caisson breakwaters with varying reflection characteristics resulted in standing waves that differ from each other, shown to have a significant influence on the seabed topography. These results are in line with Nielsen’s study (1993) that sands saltated under the surface nodes of standing waves, where the near-bed velocities are most substantial, convected toward the surface antinodes by boundary-layer drift. Moreover, the crest of sand waves was formed under the surface antinodes of standing waves, and the trough of sand waves was formed under the surface antinodes. In addition, sand wave amplitude reaches its peak in the curved slit caisson with a significant reflection coefficient, and the saltation of many grains of sand would cause this phenomenon due to the increased near-bed velocity under the nodes when the reflection coefficient is getting large.