Numerical Analysis of the Hydraulic Characteristics of a Boundary Layer Streaming over Surf-Zone Using LES and Dynamic Smagorinsky Turbulence Model

Natural shoreline repeats its re-treatment and advance in response to the endlessly varying sea-conditions, and once severely eroded under stormy weather conditions, natural beaches are gradually recovered via a boundary layer streaming when swells are prevailing after storms cease. Our understanding of the boundary layer streaming over surf-zone often falls short despite its great engineering value, and here it should be noted that the most sediments available along the shore are supplied over the surf-zone. In this rationale, numerical simulation was implemented to investigate the hydraulic characteristics of boundary layer streaming over the surf zone in this study. In doing so, comprehensive numerical models made of Spatially filtered Navier-Stokes Eq., LES (Large Eddy Simulation), Dynamic Smagorinsky turbulence closure were used, and the effects of turbulence closure such as Dynamic Smagorinsky in LES and k-ε on the numerically simulated flow field were also investigated. Numerical results show that due to the intrinsic limits of k-ε turbulence model, numerically simulated flow velocity near the bottom based on k-ε model and wall function are over-predicted than the one using Dynamic Smagorinsky in LES. It is also shown that flow velocities near the bottom are faster than the one above the bottom which are relatively free from the presence of the bottom, complying the typical boundary layer streaming by Longuet-Higgins (1957), the spatial scope where boundary layer streaming are occurring is extended well into the surf zone as incoming waves are getting longer. These tendencies are plausible considering that it is the bottom friction that triggers a boundary layer streaming, and longer waves start to feel the bottom much faster than shorter waves.

LES와 Dynamic Smagorinsky 난류모형을 이용한 쇄파역에서의 경계층 Streaming 수치해석

자연 해빈은 해양환경에 따라 침·퇴적을 반복하며 고 파랑에 의해 해빈이 대규모로 침식되더라도 폭풍이 잦아들고 다시 너울이 우세한 해양환경이 회복되는 경우 점진적으로 복원되며, 이러한 해빈 복원은 경계층 streaming을 통해 이루어진다. 이처럼 경계층 streaming은 그 공학적 가치에도 불구하고 해안에서 가용한 표사의 대부분이 공급되는 쇄파 역에서의 경계층 streaming에 대한 우리의 이해는 아직 상당히 부족하다. 이러한 인식에 기초하여 본 연구에서는 쇄파역 경계층 streaming 수리특성을 살펴보기 위해 단조 해안과 사주를 포함한 해안에서의 천수 과정을 수치모의하였다. 수치 모의는 Spatially filtered Navier-Stokes Eq., LES(Large Eddy Simulation), Dynamic Smagorinsky 난류모형으로 구성된 정교한 수치모형에 기초하여 수행되었으며, 이 과정에서 k-ε 난류모형과 LES Turbulence Closure가 모의결과에 미치는 영향도 함께 살펴보았다. 모의결과 해안공학계에 잘 알려진 k-ε 난류모형의 한계로 인해 wall function에 기반한 k-ε 난류모형의 경우 LES와 비교하면 저면 인근 유속이 다소 과다하게 모의 되었다. 또한, 바닥과 가까운 해역에서의 유속이 바닥의 영향으로부터 비교적 자유로운 상층부에서의 유속보다 우월한 Longuet-Higgins(1957)가 이야기하는 전형적인 경계층 streaming이 천수 초입부에서부터 쇄파 역 깊숙이까지 존재하는 것을 확인하였다. 또한, 주기가 상대적으로 긴 경우 경계층 streaming의 세기와 생성범위는 해안 방향으로 확대되며 이러한 경향은 경계층 streaming이 바닥 인근에서 진행되는 마찰로 인한 파랑에너지손실로 결과되며 주기가 긴 경우 천수 과정이 일찍 시작된다는 사실을 상기하면 충분히 수용 가능해 보이며, Longuet-Higgins(1957)의 해석 해에서도 같은 경향을 확인할 수 있다.