1. 서 론
15세기 중엽에 대항해 시대가 시작될 무렵 거친 파랑으로 부터 자유로운 정온한 수역을 얻으려는 인류의 노력은 시작되었다. 이러한 노력으로 자연에서 쉽게 얻을 수 있는 사석을 활용한 초기형태의 사석 방파제가 완성된다. 투석으로 조성된 경사식 사석방파제는 파랑에너지가 집중되는 평균 해수면과 이웃한 사면에서 진행되는 세굴로 변형되며 종국에는 평균 해수면 인근에 단이 형성된다. 이후 사석 경사제는 안정화 단계에 진입하여 상당 기간 추가적인 변형 없이 방파제로 기능한다, 이러한 변형과정은 파랑에 대한 내구성이 수월한 방파제 단면 형태를 우리에게 시사하며, 방파제 조성단계부터 평균 해수면 인근에 단을 배치하는 berm breakwater는 같은 맥락으로 이해할 수 있다.
방파제는 물동량의 증가와 조선 기술의 발전으로 선박이 대형화되고 파랑에 대한 우리의 이해가 깊어짐에 따라 다양한 형태로 진화한다. 현재 가용한 여러 가지의 형태의 방파제 중 내습하는 파랑에너지의 소산이란 관점에서 보면 여전히 사석 경사제가 가장 효율적이나 석재를 얻기 위한 채석과정에서 발생하는 환경훼손으로 인해 최근 들어 대형 케이슨을 활용한 직립제가 주류를 이룬다. 해상을 통한 물동량의 증가와 이로 인한 선박의 대형화로 현재 직립제는 수심 20~25m인 수역에 조성되며, 사석 경사제와 견줄만한 에너지 감쇄 효과를 얻기 위해 유수실, slit 등이 순차적으로 도입되며 진화한다. 그러나 전술한 유수실, slit등은 유수실에 포획된 공기층의 이완⋅압축으로 생성되는 방사파로 인해 목적한 에너지 감쇄 효과가 효과적으로 구현되지 못하는 경우도 상당하며 이 경우 상당한 반사가 야기되며, 이렇게 형성된 정상파는 방파제 전면해역의 쇄굴로 이어져 방파제 안전을 위협할 수 있다[Fig. 1 참조]. 현재 기부 침식 방지공은 하천 수리학에서 이야기하는 평형하상[Equilibrium regime] (Sumer and Fredsoe, 2002) 혹은 최대 세굴 깊이 (Sumer and Fredsoe, 2002)등의 개념에 준거해 돌을 다층 혹은 하나의 층으로 쌓아 보호하는 방법이 주류를 이루나, 이러한 방법은 잔류수압으로 인한 액상화현상으로 인해 방지공 밑에 있던 모래가 용출되는 등의 한계를 지니며 이는 방지공의 침하로 이어진다[Fig. 2, 3 참조]. 이와 더불어 해저지형은 내습하는 파랑에 대응하여 쉼 없이 그 모양이 바뀐다는 사실에 준거하면 쇄굴이 일정량 진행되면 더는 진행되지 않는다는 Equilibrium regime 혹은 최대 쇄굴심 등의 개념은 그 물리적 근거를 찾을 수 없다는 사실이 해안공학계에서 끊임없이 제기되고 있어, 전술한 기부 보호공은 한계를 지니는 것을 알 수 있다. (Cho and Kim, 2019; Cho and Chang, 2019; Cho et al., 2019).
본 연구에서는 이러한 인식에 기초하여 더욱 합리적인 기부침식 방지공 개발을 위한 기초연구로 사석 경사식 방파제 전면해역에서의 해저지형이 내습하는 파랑과 부분반사로 형성되는 정상 파동계에 따라 변화되는 양상을 살펴보고자 한다. 이 과정에서 openFoam 기반 물리기반 지형모형인 SeoulFoam이 활용되며 수치모의는 scale effect로 인한 오류를 최소화하기 위해 실해역을 대상으로 수행된다.
2. 수치모형
수치모의는 OpenFoam 기반 tool box인 OlaFoam과 물리 기반 지형모형인 SeoulFoam (Cho, 2020) 을 활용하여 수행하였으며, 이해를 돕기 위해 2.1절과 2.2절에는 OlaFoam과 물리 기반 지형모형인 SeoulFoam을 정리하였다.
2.1 OlaFoam
OlaFoam에서 파랑 모형은 RANS(Reynolds Averaged NavierStokes equation)와 질량 보존식으로 구성되며, 자유수면은 VOF [Volume Of Fraction] 법을 활용하여 추적된다. 난류 모형으로는 k-ε 모형을 사용하였다 (Pope, 2004).
OlaFoam (Higuera et al. 2013a; Higuera et al., 2013b)에서 사용되는 파랑 모형인 RANS, 연속 방정식을 기술하면 다음과 같다.
식 (1), (2)에서 U는 속도벡터, g는 중력가속도벡터, X는 위치벡터를 각각 나타낸다. μeff는 μ+ρvturb로 정의되는 동점성계수, p*는 유사 동압력을 각각 나타낸다.
이상 유체를 해석하기 위한 VOF식을 기술하면 다음과 같으며,
여기서 γ는 VOF법에서 자유수면을 추적하기 위해 도입된 지표로, 1에 가까울수록 액체 상태, 0에 가까울수록 기체 상태를 나타내며 본 논문에서는 자유수면의 위치를 γ = 0.5로 정의하였다.
2.2 Porous media에서의 파랑모형-VARANS equations
사석 경사식 방파제의 본체는 사석으로 구성되며, 이 경우 다공성 매질에서 발생하는 형상항력과 관성항력으로 인해 파랑 모형은 적절하게 수정되어야 한다. Darcy(1856)의 사질토 내부에서의 유체 흐름에 관한 연구 이후, 다공성 매질에서의 유체 해석을 형상항력으로 기술하려는 시도는 다공성 매질에서의 유동 해석에 큰 흐름을 형성한다. 층류의 경우 선형 형상항력만으로도 해석할 수 있으나, Forcheimer (1901)는 Reynolds 수가 큰 흐름에도 적용할 수 있도록 이차형상항력 항을 추가하였으며, Polubarinova와 Kochina(1962)는 1901년에 Forcheimer가 제시한 식을 비정상류로 확장하기 위해 관성력을 추가하였다.
본 논문에서는 del Jesus et al.(2012)에 의해 제시된 VARANS [Volume Averaged Reynolds Averaged Navier Stokes Eq.]를 사용하였으며 이를 기술하면 다음과 같다.
여기서 u는 평균유속 또는 Darcy속도, n은 공극율, A와 B는 형상항력 계수, α1은 VOF 계수, c는 관성력 형상항력계수로 물성치에 종속한다.
다공성 물질에서 비어있는 공간에 존재하는 유체의 속도를 의미하는 실제 속도 U는 평균유속 또는 Darcy속도를 나타내는 u와 다음과 같은 관계식을 추종하며
선형 형상항력계수 A, 비선형 형상항력 계수 B는 다음과 같이 정의되며[Fig. 4 참조]
여기서 D50은 다공성 물질의 명목직경, KC는 Keulegan-Carpenter수를 각각 나타내며 KC는 다음과 같이 정의된다.
식 (10)에서 uM은 최대 진동 속도, To는 진동주기를 나타내며, KC는 vortex induced vibration으로 인해 야기되는 마찰을 고려하기 위해 도입되었다. 마찰 매개 변수들과 연관되어 α와 β는 조정되며, 관성력 계수 c값은 0.34로 유지된다.
본 논문에서는 del Jesus et al. (2011)의 실험결과를 토대로, 방파제를 구성하는 사석의 명목직경 Dn50 = (M50/ρr)1/3은 0.015 m, 투수계수는 0.41, 선형 형상 항력계수 A = 10,000, 비선형 형상항력 계수 B= 6로 취하였으며, 여기서 M50는 누적 분포에서 중앙에 해당하는 사석 질량, ρr는 사석 밀도를 각각 나타낸다.
2.3 물리기반 지형모형 SeoulFoam
물리기반 지형모형 SeoulFoam은 부유사 해석을 위한 이송확산 방정식과 표사수지 개념으로부터 유도된 Exner식으로 구성된다 (Cho, 2020; Jacobsen et al., 2014). 전술한 Exner 식의 유도과정에는 소류사와 바닥으로 회귀하는 부유사를 고려하였으며, 표사 초기이송, 부유사 농도 저면 경계치, 소류사 이송률 산출과정에는 Shields Diagram, 삼차원 수치 모의된 유동계 정보로부터 직접 산출된 저면 전단응력이 활용된다(Cho, 2020; Jacobsen et al., 2014).
2.3.1 이송확산 방정식
이송확산 방정식은 다음과 같이 기술할 수 있으며,
여기서 c는 단위체적당 부유사 농도, wf와 u는 각각 표사침강속도와 유속 벡터를 나타내며, v와 vt는 각각 동점성계수와 와동점성계수를 나타낸다.
SeoulFoam에서는 식(11)로 정의되는 초기⋅경계치 문제를 완결하는 데 필요한 바닥에서의 순간 체적 표사 농도 cb는 다음과 같이 산출되며 (Cho, 2020; Jacobsen et al., 2014),
여기서 co는 최대 부유사 농도를 나타내며 본 논문에서는 0.35의 값을 취하였으며, λb는 바닥에서의 순간 선형 표사 농도를 나타내며 다음과 같이 기술될 수 있다 (Bagnold, 1954; Einstein, 1950).
식(13)에서 θ ′는 Shields parameter를 나타내며 이를 기술하면 다음과 같으며,
여기서 τb는 저면 전단응력, s는 모래 비중, ρ는 해수 밀도, d는 모래 지름을 각각 나타낸다. 식 (13)에서 pEF는 표사가 움직이는 확률을 나타내며 다음과 같이 기술할 수 있다.
식 (15)에서 μd는 동적 마찰계수를 나타낸다.
2.3.2 물리기반 지형모형
물리기반 지형모형 Exner식은 다음과 같이 기술할 수 있으며 (Cho, 2020; Jacobsen et al., 2014),
여기서 ed는 공극률을 나타내며, 침식률을 나타내는 E는 다음과 같이 산출될 수 있다.
식(17)에서 n은 저면에서의 단위 외향 법선벡터를 나타내며, 따라서 외향 법선벡터 N은 다음과 같이 정의되며,
여기서 |N|은 경계 인근 육면체로 구성된 셀에서 바닥과 접하는 면의 면적을 나타낸다. 식(16)에서 중력에 의해 바닥으로 회귀하는 모래를 나타내는 D는 침강속도 wf와 유속 u의 합을 저면법선 벡터 N에 투영하여 산출되며 이를 내적을 활용하는 경우 다음과 같이 기술될 수 있다.
표사에 작용하는 중력 w은 다음과 같이 기술될 수 있으며,
바닥과 나란한 방향으로 작용하는 접선방향 성분 wr은 외적을 활용하는 경우 다음과 같이 기술될 수 있다.
표사가 ub라는 속도로 움직이는 경우 생성되는 마찰력 ff은 운동마찰계수와 수직반력으로 기술하면 다음과 같다.
모래 입자에 작용하는 항력 fD 는 다음과 같이 기술될 수 있으며,
여기서 ur 은 상대속도를 나타내며 unw, τ - ub로 정의되며, cs는 항력과 물에 잠긴 모래 입자의 부력에 의한 중량감소를 설명해 주는 계수로 Luque (1974)의 실험결과에 준거하면 다음과 같이 기술될 수 있다.
3. 수치 모의
내습하는 파랑에 노출되는 시간에 따른 해저지형의 변화와 사석방파제 기부에서의 침식양상을 살펴보기 위해 수치 모의를 수행하였으며, 수치 모의는 OpenFoam 기반 tool box인 OlaFoam과 물리기반 지형모형인 SeoulFoam을 활용하여 수행하였다. 이해를 돕기 위해 Fig. 5에는 수치수조와 사석 방파제 제원을 정리하였으며, 도시된 계산영역은 1,152,000(3600 × 320 × 1) 개의 절점을 사용하여 이산화하였다. 수심은 해안선으로부터 60∼ 80 m 떨어진 지점까지 경사가 1/10 내외인 급경사의 전 빈이 포진하며 이후 수심이 10 m 내외인 간조 단구[low tide terrace]가 넓게 포진하는 동해안 지형 특성을 고려하여 8m, 9m, 11m로 취하였으며 (Cho, 2019; Cho and Chang, 2019, Cho and Kim, 2019, Cho et al., 2019), 파형은 간조 단구에서의 수심을 고려하여 유한 수심에서의 대표적 비선형 파랑모형으로 알려진 Cnoidal wave로 취하였다(Cho, 2019; Dean and Dalrymple, 2002). 모래 직경은 동해안에 분포하는 모래 특성을 고려하여 d50 = 0.2mm, 0.4mm, 0.8mm로 취하였다(Cho, 2019; Cho and Chang, 2019, Cho and Kim, 2019, Cho et al., 2019). Table 1에는 수치모의에 사용된 파랑 제원을 정리하였다.
4. 수치 결과
4.2 파랑
Fig. 8에는 RUN1에서 관측되는 자유수면 snapshot을 순차적으로 도시하였다. Fig. 9에는 수심이 처오름과 월파에 미치는 영향을 예시하고 검증자료로 활용하기 위해 RUN1, RUN13 [h=8m]에서 관측되는 월파 양상을 비교하였다. Fig. 10에는 Wave Gauge NO 1, 2, 3, 4에서 관측되는 자유수면 시계열자료를 도시하였다[RUN 4].
목표했던 파형이 정확하게 조파되어 정연하게 사석방파제를 향하여 진행하는 것을 확인할 수 있으며[Fig. 8 참조], 수심이 얕은 경우[RUN 13] 비선형 파랑의 대표적인 성정인 뾰족한[peaky] 파마루와 얕고 긴 파곡을 관측할 수 있다 [Fig. 9(b) 참조]. 이러한 파랑 거동으로 미루어 보면 수치모의는 상당한 정도로 수행된 것으로 판단된다. 사석방파제 사면에 거치된 Wave Gauge NO4의 경우 천수 효과로 인해 파고가 증가하는 것을 확인할 수 있으며 강해진 비선형성으로 파동계에 출현하는 고차 조화성분으로 인해 이중 파마루가 형성되는 것을 관측할 수 있다[Fig. 10(b) 참조].
Fig. 11에는 Wave Gauge NO1에서 관측되는 수평 방향 유속성분의 수직 방향 분포가 하나의 단위 주기에 걸쳐 변화하는 양상을 도시하였다[RUN 7]. 위상별 흐름 특성을 모두 담아내기 위해 하나의 단위 주기를 13개의 frame으로 나누어 자세하게 구성하였다. 전 수심대역에서 먼바다 방향의 저류가 우월한 것으로 드러나 사석방파제에서 상당한 반사가 진행되는 것을 확인할 수 있다. Fig. 12, 13에는 Wave Gauge NO2, 4에서 관측되는 [u, w]로 구성된 위상공간에서 속도벡터 u가 그리는 궤적을 도시하였다[RUN 7]. 자유수면 수면에 가까워질수록 궤적이 반시계 방향으로 틀어진 모양을 띠는 것이 이채롭다. 이러한 현상은 먼바다를 향하는 저류와 방파제를 향하는 Stokes drift가 만나서 생성되는 전단력에 기인하는 것으로 판단되며, 이러한 추론의 근거는 수평 방향 유속성분의 연직 방향 분포에서 찾을 수 있다 [Fig. 11 참조]. 또한 전 수심대역에서 먼바다를 향하는 흐름의 세기가 방파제를 향하는 흐름의 세기보다 우월하다. 방파제 사면에 거치된 Wave Gauge NO4의 경우 처오름의 세기가 처내림의 세기보다 우월하다.
4.4 해저 지형변화
Fig. 16에는 사석방파제 거치가 해저지형에 미치는 영향을 살펴보기 위해 사석방파제 거치 전후의 해저지형을 비교하였다. 사석방파제 사면에서 진행되는 부분반사와 이로 인해 전면해역에 형성되는 정상파는 해저지형에 상당한 영향을 미치는 것을 쉽게 확인할 수 있다. 연흔의 크기는 증가하며 전면으로 왜도된 모양을 띠는 것을 알 수 있다.
Fig. 17에는 내습하는 파랑에 노출되는 시간에 따라 변하는 해저지형을 정리하였다[RUN1]. 왜도된 모양의 연흔이 형성되는 것을 확인할 수 있으며 노출되는 시간이 증가함에 따라 연흔의 진폭도 같이 증가하며, 방파제 기부 인근 수역에서의 쇄굴심도 깊어지는 것을 알 수 있다.
Fig. 18에는 연흔 형상에 주기가 미치는 영향을 살펴보기 위해 94 s 경과 후 RUN9, RUN12에서 관측되는 해저지형을 비교하였다. 연흔의 진폭은 비슷하였으나, 연흔의 파장은 입사 파랑의 주기가 상대적으로 긴 RUN9에서 크게 관측되었으며, 연흔의 진폭은 입사 파랑의 파고에 의해 결정되는 것을 확인할 수 있다. Fig. 19(a)에는 모래 입경 차이가 결과하는 연흔 모양 변화를 살펴보기 위해 111.2 s 경과 후RUN1, 2, 3에서 관측되는 연흔 형상을, Fig. 19(b)에는 RUN4, 5, 6에서 관측되는 연흔 형상을 비교하여 수록하였다.
모래가 거칠수록 연흔 진폭과 방파제 기부에서의 세굴심은 증가하며, 연흔 마루가 형성되는 위치도 먼바다 쪽으로 이동하였다[Fig. 19(a) 참조]. 이러한 경향은 주기가 상대적으로 짧은 RUN4, 5, 6에서도 유사하게 관측된다. 이러한 차이를 만드는 기제를 다루기 위해서는 본 논문에서 활용한 물리기반 지형모형인 SeoulFoam에서의 저면 전단응력 산출과정을 이해할 필요가 있으며, 이를 정리하면 다음과 같다.
저면 전단응력의 정확한 모의는 물리 기반 지형모형의 정확성을 가름하는 중요한 인자로 보인다. 현재 가용한 대부분 지형모형에서는 예외 없이 마찰계수 f, quadratic friction Law, 수심 방향으로 평균된 유동계 정보를 활용하여 저면전단응력이 산출되며, quadratic friction Law를 기술하면 다음과 같다 (Nielsen, 1993).
여기서 저면 조도와 wave Reynolds NO.에 종속하는 마찰계수 f는 Lab scale에서 수행된 수리모형실험에서 측적된 자료에 준거하여 산출한다. 이 과정에서 활용되는 수리모형 실험의 상당 부분은 flat bottom, 정상 상태의 파동계를 대상으로 수행되었다. 그러나 왜도 혹은 왜곡된 파동계의 경우 연안 방향의 흐름 지속시간이 먼바다 방향의 흐름보다 짧은 지속시간으로 인해 경계층이 충분히 발달하지 못한다. 상기한 사실을 상기하면 기존 지형모형은 비선형 파동계에서 표사 이송을 기술하기에는 한계를 지니는 것을 알 수 있다.
SeoulFoam의 경우 저면 인근 경계층에서의 유속이 정밀하게 모의 되며, 이렇게 모의 된 유동계로부터 저면 전단응력은 다음과 같이 산출된다.
이는 기존에 사용되던 마찰계수 f로부터 자유로우며 왜도 혹은 왜곡된 저면 전단응력의 모의가 가능하다. 이러한 장점으로 인해 SeoulFoam의 경우 정상 파동계의 마디[node]에서 연행된 모래가 경계층 boundary layer streaming에 의해 배[antinode]로 이송되며 형성되는 사주[sand bar]를 정확하게 재현한 바 있다 (Nielsen, 1993)[Fig. 20 참조]. 따라서 상기한 SeoulFoam에서의 저면 전단응력 산출과정으로 가늠하면, 거친 모래로 구성된 저면 인근에서의 큰 전단응력으로 많은 모래가 연행되며, 이렇게 연행된 거친 모래는 상대적으로 큰 침강속도로 바닥으로 회귀하는 양도 증가하여 큰 진폭의 연흔이 결과되는 것으로 판단된다[Fig. 21 참조].
5. 결 론
현재 가용한 여러 가지 양식의 방파제 중 사석 방파제는 내습하는 파랑에너지의 소산이란 관점에서 보면 가장 효율적이나 사면에서의 부분반사로 전면 해역에 형성되는 정상파로 인해 기부에서 상당한 침식이 진행될 수 있다. 기부에서의 세굴은 피복층의 붕괴로 이어질 수 있으며, 이러한 피복층 붕괴는 사석 경사식 방파제의 주 파괴 양식으로 여겨진다. 현재 기부 침식 방지공은 Equilibrium regime 혹은 최대 세굴 깊이와 같은 개념에 준거해 돌을 다층 혹은 하나의 층으로 쌓아 보호하는 방법이 주류를 이루나, 해양환경에 따라 쉼 없이 그 모양을 바꾸는 자연 해빈을 고려하는 경우 전술한 Equilibrium regime 혹은 최대 세굴 깊이라는 개념의 논리적 근거는 상당히 취약해 보인다. 이와 더불어 사석 적층에 준거한 기부 침식 방지 공은 정상파에 노출되는 경우 잔류수압의 증가로 인한 액상화현상으로 인해 방지공 밑에 있던 모래가 용출되는 등의 한계를 지니며 이는 방지공의 침하로 이어진다.
본 논문에서는 전술한 가설을 확인하기 위해 사석 방파제 사면에서의 진행되는 부분반사와 이로 인해 전면해역에 형성되는 정상파, 정상파 형성에 따른 해저지형 변화, 사석 방파제 기부에서의 세굴 등을 수치 모의하였다. 이 과정에서 수치 모의는 OpenFoam 기반 tool box인 OlaFoam과 물리 기반 지형모형인 SeoulFoam (Cho, 2020)을 활용하여 수행하였다. 먼저 물리 기반 지형모형인 SeoulFoam의 정성적 검증을 위해 일정 수심부에서의 파랑 전파과정과 해저 지형변화를 수치 모의하였다. 모의 결과 연흔의 파장은 내습하는 파랑의 주기에 의해 특정되며, 연흔의 진폭은 내습하는 파랑의 파고에 의해 결정되는 것을 모의 되었으며, 유한 수심 해역에서 비선형 정도를 나타내는 수심 대비 진폭 비를 나타내는 a/h = 0.18가 상대적으로 작아 연흔은 정연한 모습으로 비정형성은 찾아볼 수 없었다. 또한 사석 방파제 거치로 인해 해저지형은 상당한 변화를 겪으며, 거치 전과 비교해 연흔의 크기는 증가하였으며 파랑 진행 방향으로 왜도 된 모양을 띠는 것으로 모의 되었다. 이와 더불어 정상파에 노출되는 시간이 증가함에 따라 연흔의 진폭도 같이 증가하며, 방파제 기부 인근 수역에서의 세굴도 깊어지는 것을 확인하였다. 그 밖에도 모래가 거칠수록 연흔 진폭과 방파제 기부에서의 세굴 깊이는 증가하였으며, 연흔 마루가 형성되는 위치도 먼바다 쪽으로 이동하였다. 이러한 경향은 주기가 상대적으로 짧은 경우에서도 유사하게 관측된다. 전술한 현상은 거친 모래로 구성된 저면 인근에서의 큰 전단응력으로 많은 모래가 연행되며, 이렇게 연행된 거친 모래는 상대적으로 큰 침강속도로 바닥으로 회귀하는 양도 증가하여 큰 진폭의 연흔이 결과되는 것으로 판단된다.
상기한 모의 결과는 사석을 다층 혹은 하나의 층으로 쌓는 방식의 기부 침식 방지 공의 이론적 근거로 활용되는 Equilibrium regime과는 각을 이루는 것으로 사석 적층에 준거한 현재의 기부 침식 방지 공을 재평가하는 작업이 필요해 보인다.