초록
양빈이 수행된 해빈의 침식률 산정을 위한 물리기반 해빈 지형모형이 제시되었다. 동수역학 모형은 OpenFOAM에 기반 한 tool box인 IHFOAM으로 구성되며, Morphology 모형은 부유사를 대상으로 한 이송확산방정식, 소류사 이송을 포함한 표사 수지 개념으로부터 유도된 Exner 식으로 구성하였다. 표사 이송여부, 부유사 농도저면 경계치, 소류사 이송률 산출과정에는 Shields Diagram, 삼차원 수치 모의된 유동계 정보로부터 직접 산출된 저면 전단응력이 활용된다. 본 논문에서 제시된 지형 모형을 검증하기 위해, 경사가 1/6인 단조 해안에서의 천수과정, 쇄파과정, 이에 따른 저면 변화를 수치 모의하였다. 모의결과 비선형 천수과정에서 예상되는 왜곡되고 왜도된 저면전단응력이 비교적 정확히 모의되었다. 또한 전빈에서 진행되는 쇄파로 인해 부유되고 침식된 표사가 up-rush에 의해 후빈으로 이동되어 형성되는 swash bar와 up-rush 정점에서 방향을 바꾸어 먼 바다 방향으로 진행되는 back-wash에 의해 쓸려간 표사가 수심증가에 따라 back-wash 흐름 강도가 약해지면서 퇴적하여 형성되는 breaker bar가 성공적으로 모의되는 것을 확인할 수 있었다.
In this study, a physics-based 3D morphology model for the estimation of an erosion rate of nourished beach is newly proposed. As a hydrodynamic module, IHFOAM toolbox having its roots on the OpenFoam is used. On the other hand, the morphology model comprised a transport equation for suspended sediment, and Exner type equation derived from the viewpoint of sediment budget with the bed load being taken to accounted. In doing so, the incipient motion of sediment is determined based on the Shields Diagram, while the bottom suspended sediment concentration, the bed load transport rate is figured out using the bottom shearing stress directly calculated from the numerically simulated flow field rather than the conventional quadratic law and frictional coefficient. In order to verify the proposed morphology model, we numerically simulate the nonlinear shoaling, breaking over the uniform beach of 1/m slope, and its ensuing morphology change. Numerical results show that the partially skewed, and asymmetric bottom shearing stresses can be successfully simulated. It was shown that sediments suspended and eroded at the foreshore by wave breaking are gradually drifted toward a shore and accumulated in the process of up-rush, which eventually leads to the formation of swash bar. It is also worth mentioning that the breaker bar formed by the sediments dragged by the back-wash flow which commences at the pinnacle of up-rush as the back-wash flow gets weakened due to the increased depth was successfully duplicated in the numerical simulation.