Warning: mkdir(): Permission denied in /home/virtual/lib/view_data.php on line 81

Warning: fopen(upload/ip_log/ip_log_2024-04.txt): failed to open stream: No such file or directory in /home/virtual/lib/view_data.php on line 83

Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 84
설상사주 형성조건 하에 있는 3차원투과성잠제 주변에서 내부유속변동의 특성

설상사주 형성조건 하에 있는 3차원투과성잠제 주변에서 내부유속변동의 특성

Characteristics of Velocity Fields around 3-Dimensional Permeable Submerged Breakwaters under the Conditions of Salient Formation

Article information

J Korean Soc Coast Ocean Eng. 2017;29(6):399-409
Publication date (electronic) : 2017 December 31
doi : https://doi.org/10.9765/KSCOE.2017.29.6.399
*Dept. of Energy and Plant Eng., Catholic Kwandong University
**Dept. of Civil and Environmental Eng., Graduate School, Korea Maritime and Ocean University
***Dept. of Civil Eng., Korea Maritime and Ocean Univ., 727 Taejong-ro, Yeongdo-ku, Busan 49112, Korea
이광호*, 배주현**, 안성욱**, 김도삼,***
*가톨릭관동대학교 에너지플랜트공학과
**한국해양대학교 대학원 토목환경공학과
***한국해양대학교 건설공학과
Corresponding author; Do-Sam Kim, Dept. of Civil Eng., Korea Maritime and Ocean Univ., 727 Taejong-ro, Yeongdo-ku, Busan 49112, Korea, Tel: +82-51-410-4463, Fax: +82-51-403-0656, kimds@kmou.ac.kr
Received 2017 December 7; Revised 2017 December 22; Accepted 2017 December 26.

Abstract

본 연구는 투과성잠제 배후에서 형성되는 설상사주의 주요외력으로 작용하는 평균유속, 연안류 및 난류운동에너지 등을 포함한 유속장의 특성을 수치적으로 검토하였다. 또한, 연안류에 의한 수송유량으로부터 해안선의 지형변동도 예측하였다. 수치해석에는 오픈소스 CFD 코드인 OLAFOAM에 기초한 3차원수치파동수조에 설상사주의 형성조건 하에 있는 투과성잠제를 설치하여 수치실험을 수행하였다. 수치실험에서는 규칙파의 조건 하에서 잠제의개구폭 및 설치위치에 따른 잠제 주변에서 유속변동의 특성을 평가하였다. 수치실험 결과로부터 잠제의 개구폭이 넓을수록 강한 연안류가 형성됨과 동시에 수송유량이 증가하고, 개구폭이 좁을수록 연안류가 수렴되는 지점이 제간부에서 제두부로 이동함을 확인하였으며, 이로부터 잠제 배후에 형성되는 설상사주의 형성원인을 파악할 수 있었다. 또한, 잠제 사이의 개구폭 및 잠제와 해안선 사이의 거리에 따라 연안류는 해안선에서 난류운동에너지의 변화와 밀접한 관계를 가지는 것을 알 수 있었다.

Trans Abstract

This study numerically investigates the characteristics of the velocity field including the average flow velocity, longshore current and turbulent kinetic energy acting as the main external forces of the salient formed behind the permeable submerged breakwaters. Shoreline response is also predicted by the longshore-induced flux. In this paper, a three-dimensional numerical wave tank based on the OLAFOAM, CFD open source code, is utilized to simulate the velocity field around permeable submerged breakwaters under the formation condition of salient. The characteristics of the velocity field around permeable submerged breakwaters with respect to the gap width between breakwaters and the installing position away from the shoreline under a range of regular waves for different wave height are evaluated. The numerical results revealed that as the gap width between breakwaters increases, the longshore currents become stronger. Furthermore, as the gap width becomes narrower, the point where flow converges moves from the center of the breakwater to the head part. As a result, it is possible to understand the formation of the salient formed behind the submerged breakwaters. In addition, it was found that the longshore currents caused by the gap width between breakwaters and the installation position away from the shoreline are closely related to the turbulent kinetic energy.

1. 서 론

잠제는 대표적인 연안표사제어구조물로 잠제와 지형변동과의 관계를 규명하기 위한 다양한 연구들이 수행되어 왔다. 잠제 배후에서 지형변동과 관련된 대표적 연구로 Groenewoud et al.(1997)은 다기(多基)잠제(복수의 잠제가 해안선에 평행하게 배치된 경우)의 제두부 주변에서 측정한 유속으로부터 잠제 배후에서 세굴을 실험적으로 예측하였고, Ranasinghe et al.(2006)은 삼각형상의 다기능잠제(경관, 해안환경, 파랑제어, 표사제어 및 어초 등과 같은 많은 기능을 동시에 구비한 잠제) 배후에서 평균흐름의 패턴 및 실험과 수치해석으로부터 해안선에서 침 · 퇴적을 예측하였으며, Ranasinghe and Turner(2006)는 하나로 구성된 단일잠제에서 해안선에 대한 평행입사파와 경사입사파의 경우에 잠제 배후에서 해빈류 패턴을 비교 · 검토하여 해안선의 침 · 퇴적을 예측하였다. 또한, Ranasinghe et al.(2010)은 MIKE21 모델을 이용하여 해안선과 단일잠제와의 거리에 따른 해빈류의 패턴을 비교하여 해안선의 침 · 퇴적을 수치적으로 예측하였다. 한편, Black and Andrews(2001)는 잠제의 평면적인 배치에 있어서 하나의 잠제의 경우에 다음의 식 (1)과 같은 3개의 매개변수를 사용하여 육계사주(tombolo)와 설상사주(salient)의 형성조건을 제시하였다.

(1) (mζ)t+(myHu)x+(mxHv)y=0

여기서, Ls는 잠제의 길이, S는 해빈변형 전의 원해안선과 잠제와의 거리, Soff는 잠제 설치로 인한 변형 후의 해안선과 잠제와의 거리를 나타내며, Black and Andrews(2001)는 무차원매개변수 Ls/SSoff/Ls의 관계를 제시하였다. 그들은 Ls/S > 0.6의 경우는 육계사주가 형성되고, Ls/S < 2.0의 경우는 설상사주가 형성되는 조건으로 나타내었다. 특히, Sulis et al.(2017)식 (1)과 섬(island)에 대한 제안식 Soff/Ls = 0.4(Ls/S)−1.52(Black and Andrews, 2001)과를 결합한 Soff/Ls = 0.48(Ls/S)−1.23를 이탈리아 Sardinia 섬 주변의 많은 해빈에 형성된 설상사주를 모니터링한 결과에 적용하여 Black and Andrews(2001)가 제시한 설상사주의 형성조건이 95% 일치한다는 사실을 규명하였다.

또한, Harris and Herbich(1986)Dally and Pope(1986)는 이안제에 있어서 육계사주와 설상사주의 형성조건으로 각각 Ls/S > 1.0~1.5 및 Ls/S = 0.5~1.0를 제시하였고, 개구폭의 영향을 추가하여 다기이안제(복수의 이안제가 해안선에 평행하게 배치된 경우)에서 설상사주의 형성조건으로 GS/Ls2>0.5 (G는 개구폭)를 제시하였다. 여기서, Pilarczyk and Zeidler(1996)Harris and Herbich(1986)Dally and Pope(1986)가 제시한 결과식을 기초로 하여 다기이안제에서 설상사주의 형성조건으로 전달율을 포함한 다음의 식 (2)를 제안하였다.

(2) (mHu)t+(myHuu)x+(mxHvu)y-(mf+vmyx-umxy)Hv=-myH(gζ)x+Qu

여기서, Kt는 이안제에서 월파에 의한 단면 2차원적인 전달율로 그의 범위는 0.1~0.3이며, 개구폭 GL ≤ G ≤ Ls(L은 파장)의 범위에서 값을 취한다. 다기이안제에 대한 이러한 결과는 이안제 설치에 따른 해빈변형에만 적용되며, 다기잠제와 지형변동과의 관계에 적용하기에는 무리가 있다. 이상에서 검토한 바와 같이 잠제와 지형변동와의 관계는 대부분 2차원잠제(단면적인 관점에서 잠제)에 집중되어 있는 반면, 잠제의 3차원(단면 및 평면적인 관점에서 잠제)적인 영향과 효과에 대해서는 충분한 연구가 수행되어 있지 않은 것이 현실정이며, 특히 수치적인 관점에서는 보다 미진한 것으로 판단된다. 이는 3차원모델의 비싼 계산비용, 긴 계산시간 및 복잡한 물리적 과정 때문이라고 판단된다. 하지만, 최근 들어 MPI(Message Passing Interface), SMP(Symmetrical Multi Processor) 등을 이용한 병렬계산기법의 발달로 계산효율이 대폭 증가됨에 따라 CFD 코드를 다차원으로 직접 적용한 수치계산 예가 증가하고 있다. 따라서, 본 연구에서는 Lee et al.(2016)이 적용한 오픈 CFD 코드인 OLAFOAM(Higuera et al., 2015)을 이용하여 3차원잠제 주변에서 평균유속(해빈류)분포 및 평균난류운동에너지분포 등과 같은 수리학적 특성을 잠제의 배치와 입사파랑 특성의 변화에 따라 면밀히 검토함과 동시에 설상사주 형성과의 연관성을 규명한다. 한편, 2차원잠제에서 OLAFOAM에 대한 수치모델의 구성은 Lee et al.(2016)을, 3차원잠제에서 OLAFOAM에 대한 타당성의 검증과 파고분포 및 평균수위변화 등은 전편의 논문인 Lee et al.(2017)을 각각 참조하기 바란다.

2. 수치해석이론 및 조건

2.1 수치해석이론

본 연구에서 적용한 수치해석모델인 OLAFOAM은 비혼합 및 비압축성을 갖는 혼상유체(액체 및 기체)의 거동을 해석하기 위하여 식 (3)의 연속방정식과 식 (4)의 RANS(Reynolds-Averaged Navier-Stokes) 방정식을 유한체적법에 기초하여 이산화한다.

(3) uixi=0
(4) ρuit+xj[1ϕρuiuj]=-p*fxi+ϕgjXjρxi+xj[μeffuixj]-[CT]

여기서, t는 시간, xi는 공간벡터좌표, ρ는 유체밀도, ui(i = x, y, z)는 유체속도의 Cartesian 성분, p*는 Pseudo-dynamic 압력, Xj는 위치벡터, φD50은 해석영역 내에 위치하는 투과성구조물의 공극률과 중앙입경을 각각 나타내며, μeffμ + νturbμ는 점성계수이고, νturb는 난류동점성계수이며, [CT]는 투과성구조물에 의한 저항 항으로 다음의 식 (5)(6)으로 주어진다.

(5) [CT]=Aui+B|u|ui+Cuit
(6) {A=α(1-ϕ)3ϕ2μD502B=β(1+7.5KC)1-ϕϕ2ρD50C=0.34

여기서, αβ는 선형 및 비선형마찰계수로 본 연구에서는 Jensen et al.(2014)이 층류흐름, 층류와 난류가 혼합된 흐름 및 난류흐름에 대한 수치실험결과를 바탕으로 제시한 α = 500 과 β = 2.0을 각각 적용하였다. 또한, KC는 Keulegan-Carpenter 수로 TouM/D50φ를 나타내며, To은 진동주기이고, uM은 유체입자의 최대진동속도이다.

액체와 기체의 경계면은 계산격자 내에서 유체가 점유하고 있는 체적을 나타내는 VOF 함수 F를 이용하는 VOF 법을 적용하며, VOF 법에서는 액체셀의 경우 F = 1, 기체셀의 경우 F = 0, 그리고 액체와 기체를 모두 포함하고 있는 경계셀의 경우 0 < F < 1의 값을 갖는다. VOF 함수 F식 (7)의 이류방정식에 의해 계산된다.

(7) Ft+1ϕuiFxi+1ϕuicF(1-F)xi=0

위의 식의 마지막 항에 나타낸 uic는 인공압축항으로 |uic| =min[cα|ui|, max(|ui|)]와 같이 주어지며, 본 연구에서는 cα= 1을 적용하였다. 또한, 각각의 계산셀에 있어서의 유체밀도와 점성계수는 VOF 함수를 가중치로 적용하여 다음과 같이 산정할 수 있다.

(8) ρ=ρwaterF+ρair(1-F)
(9) μ=μwaterF+μair(1-F)

난류동점성계수 νturb는 LES 모델을 적용하여 산정하였으며, 그 밖의 OLAFOAM에 대한 상세한 수치해석기법은 Lee et al.(2016) 및 OLAFOAM 매뉴얼(https://gitgub.com/phicau/OLAFOAM/master/reference.zip)을 참조하기 바란다.

2.2 수치해석조건

본 연구의 수치시뮬레이션을 위해 Fig. 1에 나타낸 바와 같이 길이 12.5 m, 폭 10 m, 높이 0.4 m의 크기를 갖는 3차원 수치파동수조를 구성하고, 수심은 바닥의 경사를 고려하지 않는 일정수심 h = 25 cm를 적용하였으며, 수조의 끝단에는 1:20의 경사를 갖는 사빈을 설치하였다. 수치파동수조 내에 3차원투과성잠제를 설치하여 잠제 사이의 개구폭 및 잠제의 설치위치에 따른 수치시뮬레이션을 수행하였다.

Fig. 1.

Illustrative sketch of wave basin and submerged structure for numerical analysis.

이하에서는 Mizutani et al.(1998)과 같은 기존의 수리실험과 잠제 천단 상에서 강제쇄파의 유도 및 수치계산에서 소요되는 계산시간 등을 참조하여 잠제의 형상, 제체의 매질 및 사빈의 매질 등을 가정하였다. 먼저, Fig. 1에 보이는 바와 같이 해석대상의 잠제에서 천단고는 21 cm, 천단수심은 4 cm, 천단폭은 1.05 m의 크기를 각각 가지며, 사면경사는 2:1로 설정하였다. 그리고, 잠제의 제체는 공극율 φ = 0.33, 중앙입경 D50 = 3.0 cm인 투과성매질로 구성되며, 잠제의 배후에 설치되는 사빈은 공극률이 φ = 0.3, 중앙입경이 D50 = 0.2 cm인 모래로 각각 가정되었다. 잠제의 구성재료가 갖는 저항계수로 Jensen et al.(2014)이 제시한 α = 500, β = 2.0을 적용하였으며, 사빈의 경우는 모래 입경에 해당하는 D50 = 0.2 cm의 다공성매질에 대한 Billstein et al.(1999)의 실험으로부터 산정된 α = 0.0, β = 3.0을 각각 적용하였다. 또한, Black and 평균유속의 공간분포 Fig. 2는 입사파고 H = 5 cm인 경우 개구폭의 변화에 따른 Andrews(2001)가 제시한 설상사주가 형성되는 조건인 Ls/S < 2.0의 범위를 고려하여 잠제의 설치위치를 선정하고, 각각의 조건에 따른 잠제 사이의 개구폭을 설정하였다. Table 1에 수치실험에 이용된 잠제의 배치와 입사파랑에 대한 조건을 나타낸다.

Fig. 2.

Spatial distribution of mean velocities at three vertical layers of horizontal plane.

Incident wave condition and 3-dimensional layout of submerged breakwater

한편, 계산영역의 격자는 수평방향으로 일정격자를 적용하고, 연직방향으로는 가변격자를 적용한 Δx = 3 cm, Δy = 3 cm, Δz = 0.7~1.4 cm로 구성하였으며, 구성된 격자에 대한 Courant 상수는 0.3 이하가 되도록 시간간격을 조정하면서 수치실험을 수행하였다. 경계조건으로는 계산영역의 모든 벽면에서 Neumann 조건인 ∂f/∂xi = 0을 적용하였다.

3. 수치해석결과

3.1 개구폭의 변화에 따른 평균유속의 공간분포

Fig. 2는 입사파고 H = 5 cm인 경우 개구폭의 변화에 따른 CASES 1~3에 대한 잠제 주변에서 시간평균유속을 나타낸 것이다. 각 그림에서 (i)는 바닥 근방에서 평균유속분포로 저면에서 z = 0.02m인 지점에서의 평균유속분포를 나타내며, (ii)는 수심의 중간위치에서 평균유속을 나타내고, (iii)는 정수면 상에서 평균유속분포를 각각 나타낸다. 먼저, 바닥 근방에서 평균유속을 살펴보면 잠제의 천단이 위치하는 곳에서는 육측(onshore side)으로, 개구부에서는 해측(offshore side)으로의 흐름이 탁월하게 발생하며, 특히 개구부에서 해측흐름은 개구폭이 좁을수록 보다 탁월하게 발생한다. 제두부에서 이러한 육측과 해측으로의 흐름에 의해 잠제 전면의 개구부 근방의 양측에서 순환류의 셀이 형성되는 것을 확인할 수 있다. 한편, 해안선 근방에서는 순환셀이 형성되지 않지만 개구폭이 넓은 경우에는 연안류성분이 탁월해지고, 반면에 개구폭이 좁은 경우에는 이안류성분이 탁월해지는 것을 확인할 수 있다. 수심의 중간위치의 경우를 살펴보면 바닥 근방에서와 유사하게 잠제 천단 상에서는 육측흐름이, 개구부에서는 해측흐름이 탁월하게 발생하며, 또한 개구폭이 좁을수록 해측흐름이 보다 탁월하게 발생한다. 그리고, 잠제 전면 및 개구부 제두부 근방에서는 바닥 근방에서와 같이 순환셀이 형성되고, 개구폭이 넓은 경우에는 개구부의 양측 제두부에서도 순환셀이 형성되며, 여기서도 해안선 근방에서는 개구폭이 넓은 경우에 연안류성분이 탁월해지고, 반면에 개구폭이 좁은 경우에 이안류성분이 탁월해지는 것을 확인할 수 있다. 수면에서는 개구부를 제외하면 전체적으로 육측흐름이 지배적으로 나타나며, 개구폭이 넓을수록 해안선 부근에서 연안류가 강하게 나타나고, 잠제 천단 상으로 유입되는 강한 육측흐름과의 상호간섭에 의해 양측벽 근방에서 순환셀이 형성된다.

Fig. 2는 입사파고 H = 5 cm인 경우 개구폭의 변화에 따른 CASES 1~3에 대한 잠제 주변에서 시간평균유속을 나타낸 것이다. 각 그림에서 (i)는 바닥 근방에서 평균유속분포로 저면에서 z = 0.02m인 지점에서의 평균유속분포를 나타내며, (ii)는 수심의 중간위치에서 평균유속을 나타내고, (iii)는 정수면 상에서 평균유속분포를 각각 나타낸다. 먼저, 바닥 근방에서 평균유속을 살펴보면 잠제의 천단이 위치하는 곳에서는 육측(onshore side)으로, 개구부에서는 해측(offshore side)으로의 흐름이 탁월하게 발생하며, 특히 개구부에서 해측흐름은 개구폭이 좁을수록 보다 탁월하게 발생한다. 제두부에서 이러한 육측과 해측으로의 흐름에 의해 잠제 전면의 개구부 근방의 양측에서 순환류의 셀이 형성되는 것을 확인할 수 있다. 한편, 해안선 근방에서는 순환셀이 형성되지 않지만 개구폭이 넓은 경우에는 연안류성분이 탁월해지고, 반면에 개구폭이 좁은 경우에는 이안류성분이 탁월해지는 것을 확인할 수 있다. 수심의 중간위치의 경우를 살펴보면 바닥 근방에서와 유사하게 잠제 천단 상에서는 육측흐름이, 개구부에서는 해측흐름이 탁월하게 발생하며, 또한 개구폭이 좁을수록 해측흐름이 보다 탁월하게 발생한다. 그리고, 잠제 전면 및 개구부 제두부 근방에서는 바닥 근방에서와 같이 순환셀이 형성되고, 개구폭이 넓은 경우에는 개구부의 양측 제두부에서도 순환셀이 형성되며, 여기서도 해안선 근방에서는 개구폭이 넓은 경우에 연안류성분이 탁월해지고, 반면에 개구폭이 좁은 경우에 이안류성분이 탁월해지는 것을 확인할 수 있다. 수면에서는 개구부를 제외하면 전체적으로 육측흐름이 지배적으로 나타나며, 개구폭이 넓을수록 해안선 부근에서 연안류가 강하게 나타나고, 잠제 천단 상으로 유입되는 강한 육측흐름과의 상호간섭에 의해 양측벽 근방에서 순환셀이 형성된다.

Fig. 3Fig. 2와 동일한 조건인 입사파고 H = 5 cm인 경우 산정된 유속성분을 저면에서 수면까지 수심방향으로 적분하고, 동시에 시간평균한 평균유속(해빈류)의 공간분포를 나타낸다. 그림으로부터 잠제 천단 상에서는 육측으로의 흐름이, 반면에 개구부에서는 해측으로 흐름이 나타나며, 이는 일반적으로 잠제에서 잘 알려진 플럭스수송현상과 일치한다. 개구부 근방에서 잠제 천단 상에서 육측으로, 그리고 개구부에서 해측으로 급변하는 이러한 흐름 특성으로 인하여 개구폭에 관계없이 개구부 근방의 잠제 전면에서 순환셀이 형성된다. 이는 다기이안제와 다기잠제 배후에서 순환셀의 방향을 실험적으로 검토한 Loveless and MacLeod(1999)의 결과와도 일치한다. 즉, 이안제의 경우 제두부에서 육측흐름이 발생하여 제간부로 향하는 순환류가 발생하는 반면, 잠제의 경우 제간부에서 육측흐름이 발생하여 제두부로 향하는 순환류가 발생한다는 것을 실험적으로 규명한 Loveless and MacLeod(1999)의 결과와 일치하며, 단일잠제에 대한 Nobuoka et al.(1996)에서도 유사한 순환패턴이 나타나는 것을 알 수 있다.

Fig. 3.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

또한, 해안선 부근에서는 개구폭에 관계없이 순환셀이 형성되지만 개구폭이 가장 좁은 CASE 3보다 개구폭이 넓은 CASE 2에서 순환셀의 크기가 해안선과 평행한 방향으로 보다 확장되고, 개구부 배후의 해안선 부근에서 순환셀이 형성된다. 한편, 개구폭이 가장 넓은 CASE 1에서는 순환셀의 크기가 작아지는 반면, 연안류가 형성되는 것을 확인할 수 있다. 이것은 개구부 배후 및 제간부 배후의 각각 해안선 부근에서 나타나는 평균수위(Lee et al., 2017)의 차이에 기인하는 연안류로 판단된다.

Fig. 4는 입사파고 H = 7 cm인 경우에 Fig. 3과 같이 수심적분 및 시간평균된 평균유속의 공간분포를 나타낸다. 입사파고 H = 5 cm의 경우와 동일하게 개구부 근방의 잠제 전면에서 순환셀이 형성되고, 이는 제두부 근방의 제간부에서 육측흐름과 개구부를 통한 해측흐름과의 간섭에 의해 발생한다. 전술한 입사파고 H = 5 cm인 Fig. 2의 경우와는 상이하게 개구폭에 관계없이 연안류가 형성되지만 개구폭이 넓을수록 보다 강한 연안류가 형성되며, 사빈 경사면 상에서 보다 명확한 순환셀이 형성된다.

Fig. 4.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

3.2 잠제의 설치위치에 따른 평균유속의 공간분포

Fig. 5는 잠제 위치에 따른 수심적분 및 시간평균된 평균유속분포를 나타내며, 입사파고가 H = 7 cm, 개구폭이 G = 2.5 m인 경우에 잠제 위치 S = 6.5m를 중심(CASE 5)으로 해안선에서 1. 0m 더 가까운 S = 5.5m의 경우(CASE 7), 해안선에서 0.5 m 더 가까운 S = 6.0m의 경우(CASE 8)와 해안선에서 0.5 m 더 먼 S = 7.0m 경우(CASE 9)를 각각 제시한다. 결과를 살펴보면 개구부 근방의 잠제 전면에서는 모두 순환셀이 형성되며, 이는 Loveless and MacLeod(1999)에서도 지적되어 있는 바와 같이 제두부 부근에서 발생하는 제간부에서 육측흐름과 개구부에서 해측흐름과의 상호작용에 의해 발생된다. CASES 7과 8을 CASE 5와 비교하면 사빈 경사면 상에서 순환셀이 형성되는 현상은 동일하지만 CASE 5의 경우가 연안류의 형성이 보다 명확하다. 여기서, 잠제 위치가 해안선에 멀어질수록 CASE 5의 경우 사빈 경사면에서 종단방향으로 장축을 가진 타원형의 순환셀이 형성되고, 잠제 위치가 해안선에서 가장 멀리 떨어진 CASE 9의 경우 연안류성분이 강해지면서 순환셀의 형성은 불완전하게 되며, 최내측의 해안선 부근에서 순환셀이 형성되는 것을 알 수 있다.

Fig. 5.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

3.3 연안류분포

Ranasinghe and Turner(2006)의 결과에 의하면 해안선에 평행하게 설치된 단일잠제에 파가 직각방향으로 입사하는 경우 잠제 제간부 배후의 해안선 부근에서는 양측으로 발산하는 연안류가 발생하여 침식이 발생(Ranasinghe and Turner(2006)의 Fig. 11을 참조)하는 반면, 파가 해안선에 경사지게 입사하는 경우 파가 잠제를 통과한 후 해안선에서 연안류를 방해하여 연안류의 속도를 감소시키고, 이로 인하여 흐름이 정체되어 퇴적이 발생(Ranasinghe and Turner(2006)의 Fig. 12를 참조)한다. 또한, Ranasinghe et al.(2010)의 결과를 살펴보면 파가 해안선에 평행하게 입사하는 경우에도 잠제와 해안선과의 거리에 따라 해빈류 패턴이 달라지며, Ranasinghe and Turner(2006)와 동일하게 해안선 부근에서 양측으로 발산하는 연안류가 발생하는 경우에는 침식이 발생하고, 연안류가 수렴하는 지점에서 퇴적이 발생한다고 결론짓고 있다. 즉, 연안류는 해안선의 침 · 퇴적과 같은 지형변동에 큰 영향을 주며, 설상사주의 형성에 대한 예측에서 중요한 물리적 파라미터라고 할 수 있다.

Fig. 6Figs. 3~5에서 조파지점으로부터 10.5~11.4 m 이 격된 부근에 형성된 연안류를 x방향으로 평균하여 y방향에 대해 나타낸 것이며, (+)는 제두부에서 제간부로 향하는 연안류를, (−)는 제간부에서 제두부로 향으로 연안류를 각각 나타낸다. 개구폭에 따른 연안류의 변동을 살펴보면 입사파고가 H = 5 cm인 경우 Fig. 6(a)에서 CASES 2와 3은 y축을 따라 연안류가 진동하는 것을 알 수 있고, 잠제의 제두부 근방에서 감소하는 경향을, 제간부에서 증가하는 경향을 각각 나타내며, 이러한 결과는 순환셀이 형성되기 때문이다. 반면, CASE 1은 전술한 바와 같은 약간의 진동은 있지만 제두부 및 제간부에서 연안류가 상대적으로 안정된 큰 값을 나타낸다. 입사파고가 H = 7 cm인 경우의 Fig. 6(b)에서는 개구폭이 넓은 CASES 4와 5에서 연안류가 보다 빠르게 나타나고, 이는 제두부에서 최대치를 나타낸 후 제간부로 갈수록 감소되는 경향을 나타낸다.

Fig. 6.

Spatial distribution of longshore current velocities.

잠제 위치에 따른 연안류의 차이를 나타낸 Fig. 6(c)Fig.5에서 언급한 바와 같이 CASE 9의 경우 사빈 경사면 상에서의 순환셀이 연안류 형성에 영향을 미쳐 연안류의 속도가 상대적으로 빠르게 나타나며, CASE 8의 경우는 CASE 5보다 사빈 경사면 상에서 형성된 순환셀이 연안류의 흐름방향으로 확장되어 연안류 속도의 증가에 영향을 준 것으로 판단된다. CASE 7의 경우도 CASE 8과 같이 CASE 5보다 사빈 경사면 상에서의 순환셀이 연안류의 흐름방향으로 확장되어 있지만 잠제 배후의 공간이 좁기 때문에 연안류가 다소 약하게 형성된 것으로 판단된다.

본 연구는 다기로 구성된 잠제로 단일잠제에 대한 Ranasinghe et al.(2006)과는 상이하게 파가 해안선에 직각으로 입사하는 경우 제간부 배후 해안선 부근에서 양측으로 발산하는 연안류가 발생하지 않으며, 개구부 배후 해안선 부근에서 상승된 평균수위에 의해 제간부측으로 이동하는 연안류가 발생하는 것을 확인할 수 있다.

여기서, 연안류가 수렴되는 지점에서 퇴적이 발생된다는 Ranasinghe et al.(2010)의 결과에 따라 Fig. 6에서 제시된 연안류의 속도가 (+)에서 (−)로 변하는 변곡점을 본 연구에서는 연안류가 수렴되는 지점으로 판단하고, 그 지점에서 각 방향(제두부에서 제간부측으로 향하는 방향과 제간부에서 제두부측으로 향하는 방향)으로 연안류에 의해 수송되는 유량을 산정하여 Table 2에 그 값을 제시한다.

A longshore current’s inflection point on shoreline and total transport flow discharge due to longshore current

Table 2에 제시한 수송유량은 연안류의 흐름속도와 x-z 평면에서 연안류가 형성된 면적을 곱하고 누적한 값으로 추정되었다. 그리고, CASES 2와 3의 경우 Fig. 3에 나타난 바와 같이 순환셀에 의해 연안류가 완전히 형성되지 않아 연안류가 수렴되는 지점을 선정하기 어려웠기 때문에 CASE 1과 동일하게 잠제 제간부 배후의 지점을 연안류가 수렴되는 지점으로 선정하였다. 결과를 살펴보면 CASES 1~6의 경우는 개구폭이 넓을수록 연안류에 의한 수송유량이 증가하는 것을 알 수 있고, 입사파고가 H = 7 cm인 경우의 CASES 4~6은 개구폭이 좁을수록 연안류의 변곡점이 제간부 중앙에서 제두부측으로 이동하는 것을 확인할 수 있다. 이러한 결과로부터 퇴적이 발생되는 지점이 이동된다는 것을 알 수 있고, 또한 이로부터 설상사주가 형성되는 지점을 예측할 수 있다. 한편, 잠제 위치에 따른 연안류의 수송유량은 잠제 위치에 상관없이 CASE 9, CASE 8, CASE 5, CASE 7의 순으로 크게 나타났으며, 연안류의 변곡점도 잠제 위치에 상관없이 CASE 7, CASE 5, CASE 9, CASE 8의 순으로 제간부에 근접하여 발생되었다.

3.4 개구폭의 변화에 따른 평균난류운동에너지분포

Fig. 7Figs. 34에서와 같이 저면에서 수면까지 수심방향으로 적분한 유속을 적용하여 평균난류운동에너지를 나타낸 것이며, 이는 ux'2+uy'2+uz'2/2 으로 정의되는 난류에너지를 시간평균한 ux'2+uy'2+uz'2/2을 나타낸 것이다. 여기서, u'x, u'y, u'z는 각각 x, y, z방향으로 난류속도성분이다. 이러한 난류에너지는 난류성분의 크기와 그로 인한 파랑에너지의 소산 정도를 판정하는 중요한 요소로 전술한 바와 같이 파고의 변화 및 평균수위의 변화에 영향을 미친다. 여기서, Fig. 7(a)는 입사파고가 H = 5 cm인 경우를, Fig. 7(b)는 입사파고가 H = 7 cm인 경우를 각각 나타낸다. 그림으로부터 잠제 천단 상과 개구부 배후 해안선 부근에서 평균난류에너지가 전반적으로 집중 · 증가하고, 개구폭이 넓을수록 개구부 배후 해안선 부근에서 평균난류운동에너지가 상대적으로 크게 나타나며, 또한 종단방향으로 넓게 분포하는 것을 알 수 있다.

Fig. 7.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in horizontal plane.

Fig. 8은 해안선과 평행한 방향으로 시 · 공간평균난류운동에너지분포를 나타내며, x = 10.2m에서 해안선 위치 x = 11.47 m 사이의 평균난류운동에너지의 공간분포를 평균한 결과이다. 시 · 공간평균난류운동에너지분포는 파랑에너지의 소산과 연관되어 있으므로 파고분포와 관계를 가진다. 또한, 파고의 자승에 비례하는 잉여응력(radiation stress)은 쇄파 시에는 증가하고, 쇄파 후에는 감소하므로 해안선 근방에서 평균수위 상승이 발생하여 이에 균형을 취하는 것으로 알려져 있다. 결과를 살펴보면 전술한 Fig. 7에서와 같이 개구부 배후의 해안선 근방에서 난류운동에너지가 증가하고, 더불어 개구폭이 넓을수록 이러한 경향이 증가하는 것을 확인할 수 있으며, 이는 해안선 부근에서 평균수위 상승의 결과로 작용한다.

Fig. 8.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in longshore current region.

3.5 잠제의 설치위치에 따른 평균난류운동에너지분포

Fig. 9는 전술한 Fig. 7에서와 동일한 방법으로 잠제 위치에 따른 평균난류운동에너지분포를 나타낸 결과이다. 난류운동에너지의 분포 형상은 전술한 Fig. 7과 동일하게 잠제 천단 상 및 개구부 배후 해안선 부근에서 집중 · 증가되는 것을 알 수 있다. 반면, 해안선과 잠제 사이의 이격거리 S의 변화에 따른 뚜렷한 경향은 확인되지 않는다.

Fig. 9.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in horizontal plane.

Fig. 10Fig. 8과 동일한 방법으로 해안선과 평행한 방향으로 시 · 공간평균난류운동에너지분포를 나타낸 것이다. 그림으로부터 개구부 배후 해안선 부근에서 CASES 8과 9의 경우 난류운동에너지가 다소 큰 값을 나타내지만 전제적으로는 제두부에서 제간부로 갈수록 난류운동에너지가 작아지는 변동양상을 보이고, 각각에서 값의 차이는 그다지 크기 않다.따라서, 이러한 변화는 전술한 파고분포나 평균수위 상승에서의 변화와 매우 유사하다는 것을 알 수 있다.

Fig. 10.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in longshore current region.

이상으로부터 개구폭의 변화에 따른 개구부 배후의 평균수위변화 및 파고변화는 평균난류운동에너지와 크게 연관되어 있는 것을 확인할 수 있고, 평균난류운동에너지가 높을수록 개구부 배후에서 평균수위가 보다 상승하며, 이로 인하여 연안류의 속도에 간접적인 영향을 미치는 것으로 판단된다.

4. 결 론

본 연구에서는 오픈소스 CFD 코드인 OLAFOAM을 이용하여 다기투과성잠제 주변 및 배후에서 형성되는 평균유속, 연안류 및 난류운동에너지 등의 변동특성을 검토하였다. 본 연구로부터 도출된 중요한 결과를 다음에 요약 · 기술한다.

(1) 수치실험에 의해 산정된 평균유속의 공간분포로부터 개구폭 및 잠제 위치에 상관없이 제두부 전면에서는 순환셀이 대부분 형성되었고, 제두부 부근에서는 제간부로부터 유입되는 흐름으로 인하여 개구부로 순회하여 빠져나가는 순환셀이 형성됨을 확인하였다. 또한, 개구폭이 넓어질수록 해안선 부근에서 해안선과 평행한 방향으로 순환셀의 크기가 커지다가 연안류 형성에 관여한 후에 다시 크기가 작아지는 경향을 나타내었다. 잠제 위치가 해안선에서 멀어질수록 사빈 경사면에서 해안선에 직각방향으로 장축을 가진 타원형의 순환셀을 형성하며, 이는 연안류에 직접적인 영향을 미치는 것으로 판단되었다.

(2) 잠제의 배후에서의 연안류는 개구폭이 넓을수록 강하게 형성됨과 동시에 수송유량도 증가하였으며, 연안류가 수렴되는 지점은 개구폭이 좁을수록 제간부에서 제두부로 이동되는 것으로 추정되었다. 이러한 결과로부터 저질이 퇴적되는 지점을 예상할 수 있기 때문에 잠제 배후에 형성되는 설상사주와 같은 지형변동을 예측할 수 있을 것으로 판단된다.

(3) 평균난류에너지는 개구폭이 넓을수록 개구부 배후의 해안선 부근에서 시 · 공간누적평균난류에너지가 크게 나타났고, 이는 개구부 배후 해안선 부근에서 평균수위의 상승 정도와 직접적인 연관성을 갖는 것으로 판단되었으며, 더불어 연안류의 형성과 속도에도 직접적인 영향을 미치는 것으로 판단된다.

Acknowledgements

이 논문은 2017년도 가톨릭관동대학교 학술연구비에 의해 지원되었음(CKURF-201406810001)의 지원을 받아 수행된 사업임.

References

Black KP, Andrews CJ. 2001;Sandy shoreline response to offshore obstacles Part 1: Salient and tombolo geometry and shape. J. Coastal Research :82–93.
Dally WR, Pope J. 1986. Detached breakwaters for shore protection (No. CERC-TR-86-1) Coastal Engineering Research Center. Vicksburg MS:
Groenewoud MD, van de Graaff J, Claessen EW, van der Biezen SC. 1997;Effect of submerged breakwater on profile development. Proceedings of Coastal Engineering, ASCE :2428–2441.
Harris MM, Herbich JB. 1986;Effects of breakwater spacing on sand entrapment. J. Hydraulic Research 24(5):347–357.
Higuera P, Losada IJ, Lara JL. 2015;Three-dimensional numerical wave generation with moving boundaries. Coastal Engineering 101:35–47.
Jensen B, Jacobsen NG, Christensen ED. 2014;Investigations on the porous media equations and resistance coefficients for coastal structures. Coastal Engineering 84:56–72.
Lee KH, Bae JH, An SW, Kim DS, Bae KS. 2016;Numerical analysis on wave characteristics around submerged breakwater in wave and current coexisting field by OLAFOAM. Journal of Korean Society of Coastal and Ocean Engineers 28(6):332–349.
Lee KH, Bae JH, An SW, Kim DS. 2017;Characteristics of water surface variations around 3-dimensional permeable submerged breakwaters under the conditions of salient formation. Journal of Korean Society of Coastal and Ocean Engineers (in press).
Loveless JH, MacLoed B. 1999. The influence of set-up behind detached breakwaters. In : Proc. 26th International Conference on Coastal Engineering. ASCE. Long Island, USA. p. 2026–2041.
Mizutani N, Mostafa AM, Iwata K. 1998;Nonlinear regular wave, submerged breakwater and seabed dynamic interaction. Coastal Engineering 33(2):177–202.
Nobuoka H, Irie I, Kato H, Mimura N. 1997;Regulation of nearshore circulation by submerged breakwater for shore protection. Proceedings of Coastal Engineering, ASCE :2391–2403.
Pilarczyk KW, Zeidler RB. 1996. Offshore breakwaters and shore evolution control.
Ranasinghe R, Larson M, Savioli J. 2010;Shoreline response to a single shore-parallel submerged breakwater. Coastal Engineering 57(11):1006–1017.
Ranasinghe R, Turner IL. 2006;Shoreline response to submerged structures: a review. Coastal Engineering 53(1):65–79.
Ranasinghe R, Turner IL, Symonds G. 2006;Shoreline response to multi-functional artificial surfing reefs: A numerical and physical modelling study. Coastal Engineering 53(7):589–611.
Sulis A, Balzano A, Cabras C, Atzeni A. 2017;On the applicability of empirical formulas for natural salients to Sardinia (Italy) beaches. Geomorphology 286:1–13.

Article information Continued

Fig. 1.

Illustrative sketch of wave basin and submerged structure for numerical analysis.

Fig. 2.

Spatial distribution of mean velocities at three vertical layers of horizontal plane.

Fig. 3.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

Fig. 4.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

Fig. 5.

Spatial distribution of time- and depth-averaged mean velocities in horizontal plane.

Fig. 6.

Spatial distribution of longshore current velocities.

Fig. 7.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in horizontal plane.

Fig. 8.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in longshore current region.

Fig. 9.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in horizontal plane.

Fig. 10.

Spatial distribution of time- and depth-averaged turbulent kinetic energies in longshore current region.

Table 1.

Incident wave condition and 3-dimensional layout of submerged breakwater

CASE No. Wave maker theory Wave height H (cm) Wave period T (s) Gap width G (m) Submerged breakwater length Ls (m) Distance between shoreline and submerged breakwater S (m) Ls/S Wave breaking on crown
CASE 1 Stokes II 5 1.4 3.0 7.0 6.5 1.08 No
CASE 2 2.5 7.5 1.15
CASE 3 2.0 8.0 1.23
CASE 4 7 3.0 7.0 1.08 Yes
CASE 5 2.5 7.5 1.15
CASE 6 2.0 8.0 1.23
CASE 7 2.5 7.5 5.5 1.36
CASE 8 6.0 1.25
CASE 9 7.0 1.07

Table 2.

A longshore current’s inflection point on shoreline and total transport flow discharge due to longshore current

CASE No. Inflection point Transport flow discharge (cm3/s) (From head of breakwater to middle of breakwater) Transport flow discharge (cm3/s) (From middle of breakwater to head of breakwater) Total transport flow discharge (cm3/s)
CASE 1 y = 10.0 m 8,155.42 0 8,155.42
CASE 2 y = 10.0 m 6,227.68 0 6,227.68
CASE 3 y = 10.0 m 2,579.38 0 2,579.38
CASE 4 y = 10.0 m 25,631.62 0 25,631.62
CASE 5 y =9.0 m 22,857.21 219.85 23,077.06
CASE 6 y =8.11 m 8,490.11 2,073.05 10,563.16
CASE 7 y = 10.0 m 20,343.89 0 20,343.89
CASE 8 y = 8.38 m 23,161.30 1,713.69 24,874.99
CASE 9 y = 8.64 m 25,411.64 1,075.66 26,487.30