1. 서 론
우수관거 내부 슬러리 흐름은 이송⋅침전 및 퇴적의 세가지 현상학적으로 구분할 수 있다. 이송과정의 명확한 정의는 고체로 된 입자 알갱이들이 우수관거를 타고 유체와 함께 원활히 흘러가는 것을 의미한다. 이송현상은 수자원공학의 수리⋅수문학적 해석에 있어 중요한 과정이고 실제 주변에서 일어나는 현상이다. 대게 퇴적물의 이송은 하천, 저수지와 관개배수를 비롯한 과정에서 일어나며 해양환경과 같은 분야에서도 주로 발생한다. 이송현상은 물리학적⋅화학적⋅생물학적으로 접근하는 복잡한 현상에 속한다. 현행 국내 우수관로 설계에 있어 관 내 퇴적물의 침전을 방지하기 위한 최소 설계유속을 규정하고 있다. 그러나 퇴적물의 거동에 초점을 두고 실제적인 거동을 규명한 연구는 미흡하다고 할 수 있으며, 유지관리 측면에 있어 반드시 필요하다. 관련 연구나 설계기준 마련을 위해 뒷받침하는 이론 및 분석적 근거가 필요한 현실이다.
많은 연구자들이 슬러리 이송에 대해 관심을 갖고 한계유속(critical values of velocity)을 이용하여 관계식을 제안하거나 임계값을 제시하여 이송⋅침전 현상에 대해 정량화 하는 연구 결과들을 제안하였다. 이러한 연구들은 퇴적물 이송현상 예측에 중요한 역할을 했다. 더불어 하상에서의 침식을 예측함으로써 수로 설계에 있어서도 퇴적물을 일정 고려함으로써 통수능 검토에도 많은 도움이 되었다. Einstein(1942), Meyer-Peter and Müller(1948), Brown(1950), Engelund and Hansen(1967), Ackers and White(1973), Van Rijn(1993) 등 많은 연구들이 하상에서의 퇴적물의 흐름을 정의하려고 노력했다. 일반적으로 기존에 제안된 모든 모델들은 일정한 유동에서 균일한 직경을 갖는 비 점착성 입자에 대해 상대적으로 유사한 결과를 제공하였다. 그러나 하천을 대상으로 한 연구결과들이 주를 이루며 우수관거 내부 흐름에 대한 연구를 통해 제시된 결과는 부족하다.
기존에 제안된 하천에서의 경험적 공식을 우수관로의 수리학적 특성에 맞게 수정하면 적용이 가능하다. 하지만 정작 중요한 것은 내부 흐름 특성이다. 우수관거 내부 흐름은 하천의 흐름과 전혀 다른 특성을 가지므로 기존의 모델을 활용할지라도 정확한 예측은 다소 무리가 있다. 더불어 대부분의 우수관거는 침전물이 발생하지 않는 상태를 유지하기 위해 유속이나 전단응력을 고려하여 설계한다. 하지만 설계의 기준은 대게 일관된 값을 고려하여 설계된다.
본 연구에서는 토사 슬러리를 고려한 2상유동(two-phase flow) 해석을 바탕으로 관 내부에서의 유동현상을 검토하였다. 관련 선행 연구로는 Song et al.(2017a), Song et al.(2017b)는 우수관로 내 토사 슬러리의 흐름에 대한 연구를 진행하였다. 다만, 본 연구에서는 이송과정 동안 침전이 발생한 이후에 나타나는 흐름의 경계층을 검토하고자 한다. 흐름의 경계층을 구분함에 있어 중요한 요소는 전단유속과 응력이다. 침전에 따른 퇴적은 결국 관 저면의 하상전단응력이 퇴적물 이동에 대한 한계전단응력 보다 작은 경우에 발생한다. 즉, 하상이나 벽면 근처에서는 매우 큰 전단응력이 발생하는데 이 힘을 견딜 수 있는 충분한 전단강도를 갖는 퇴적물만이 하상에 쌓여 침전층을 형성한다고 볼 수 있다.
본 연구는 기존의 Song et al.(2017b)이 제안한 연구 결과들을 활용하였다. 우수관로 내부로 유입 가능한 토사 슬러리를 고려하여 유입 조건을 변화시켜 이송과정에서 발생하는 침전현상을 검토하고 발생 가능한 흐름의 경계층을 검토하는 것에 목적이 있다. 더불어 토사 슬러리를 고려한 2상유동의 관측 및 해석의 어려움으로 인해 연구결과가 드물어 해당 부분이 본 연구에서 검토하고자 하는 연구목표이다. 그러한 이유로 본 논문의 2상유동 해석의 조건을 다양하게 설정하였다. 분석조건은 2장 3절에 명시하였다.
2. 계산방법 및 해석조건
2.1 수치해석 방법
본 연구에서는 ANSYS-FLUENT Workbench 13.0으로 수치해석을 진행하였다. 고체입자인 토사 슬러리와 유체로 구성된 2상유동 해석에는 Euler-Euler 모델을 적용하였다. Euler-Euler 모델은 수많은 입자들의 움직임을 표현할 때 입자를 개별적으로 추적하는 것이 아니라 분산된 흐름의 형태를 표현하므로 다상유동의 모델링에 주로 사용된다. 해석 과정에서 각 상들은 개별적으로 처리되며 각 단계마다 모든 상들을 결합한 하나의 보존 방정식이 도출된다. 각 상들 사이의 결합은 압력 및 상간의 운동량 교환 계수를 통해 이루어진다. 결국 서로 다른 상이라 하더라도 전체 체적농도의 합은 1로 일정함을 유지해야 한다. 결국 2개의 연속적인 상을 갖는 2상유체 모델로써 지배방정식은 다음과 같다.
여기서 ρ는 상의 밀도, U는 평균 속도의 장, P는 상들의 평균 압력이다. 아래 첨자 k는 k번째 연속 위상을 가리킨다. 분산된 각 상들의 방정식은 주어진 계산 Mesh의 셀 내부에서 평균되어 대표값을 얻게 된다. 분산된 각 상을 하나의 연속체로 표현하기 위해서는 높은 체적농도가 요구되며, 따라서 해당 방법은 밀도나 농도가 높은 상을 포함한 유동해석에 적합하다(Jeong, 2009).
2.2 분석모델의 검증
본 연구에서는 Matousek(1997)의 1.4∼2.0 mm의 토사 입자들에 대한 실험적 결과를 토대로 동일한 제원을 반영하여 본 연구에서의 수치해석 결과를 대조하였다. 다음의 Fig. 1은 우수관 내부에 1.4∼2.0 mm의 입자가 26%의 체적농도 및 2.5∼3.0 m/s의 유속조건에서의 이송과정에 대한 해석결과이다. 모델링을 통해 도출된 결과의 흐름이 상당부분 일치한다. 결과적으로 검증에 적용된 두 가지 유속 조건의 변화에 대해 체적농도에 따른 내부 분포가 동일한 경향성을 나타내므로 해석을 위한 난류모형의 검증이 이루어졌다 판단하였다.
2.3 수치해석 조건
본 연구에서는 직경 0.6m, 관 연장 10 m의 원형 우수관로를 140,000개의 격자망 구성을 통해 모사하였고, 동일한 유입조건에서 계산효율이 높은 2차원 모델로 해석을 수행하였다. 토사 슬러리 혼합물의 유입유속, 혼합물 중 토사의 체적농도 및 토사의 입자크기를 변수로 선정하여 각각의 변화 조건에 따른 전체 63개 조건에 대해 혼합물의 유속과 입자 분포를 비교⋅검토하였다. 다음의 Table 1은 해석 대상으로서 유입되는 유체와 토사 슬러리 혼합물의 물성값 및 해석 경계조건을 나타낸다.
3. 우수관거 내부 슬러리 혼합물의 흐름 특성
3.1 우수관거 내부 토사 슬러리 혼합물의 유동 특성
관 저면에서 발생하는 침전층의 형성은 이송과정에 있어 전단경계층 형성에 중요한 원인이며 시설물의 유지관리가 가능한 주요 참고사항이 된다. 하천이나 하구의 침식이나 바람, 공기의 흐름에서 입자의 이송에 대한 연구에 사용되는 주요 매개변수에는 전단유속을 활용한다.
경계층은 전단유속의 상대적인 값에 따라 분류되며 해당값을 초과하는 힘을 갖는 입자들에 의해 퇴적층을 형성하고 나머지의 경우 이송과 침전층을 구성하게 된다. 이때 침전된 토사 슬러리 이송과정에 영향을 미치는 중요한 요소는 유동 저항(flow resistance)이다. 퇴적물은 관 저면에서 침전되어 유동 저항을 증가 시키며 통수단면적을 감소시킨다(Fig. 3).
Liu(1998)와 Butler et al.(2004)은 우수관거 내부에서의 침전물의 이동은 3가지 기본 단계를 포함한다고 언급하였다. 부상(entrainment), 이송(transport) 그리고 퇴적(sedimentation)으로 정의하였다. 세 가지 분류에 따른 전단특성을 표현하면 다음의 Fig. 4와 같다. 이때 경계층의 두께는 일정하지 않으며 전체 흐름구간의 대부분이 난류영역에 속한다.
Raudkivi(1998)는 Fig. 4와 같은 흐름의 경계층을 수치적으로 표현하고자 유동과정 중 경계층의 구분이 관 저면에서 발생하는 전단유속(u * v s
Fig. 4와 Table 2를 연관지어 서술하면 입자간 점성력이 작용하여 관 저면을 벗어나지 않고 입자가 이동하는 층을 bed load layer로 정의한다. 퇴적층 상단부로 경계층에 해당하며 큰 난류응력의 작용으로 인해 퇴적된 입자들이 이송과 침전에 따른 이동이 활발히 발생하는 지역을 transition layer로 정의하였으며, 이후 대게 난류응력만이 작용하여 이송과정을 포함하는 층을 wash load layer 및 suspended load layer로 구분하였다.
유동에 대한 마찰저항이 작거나 없기 때문에 이송 과정에서 입자들의 침강속도 역시 작은 값을 갖는다. 따라서 완전한 난류흐름으로써 난류응력이 존재한다. 관거의 흐름방향을 기준으로 전단응력이 상단부를 향하여(유속이 0이 되는 구간) 선형적으로 감소하며, 유체와 고체입자의 격한 혼입이 발생하는 구간이다. 와류나 소용돌이가 발생하나 대게 해당 구간에서의 유속은 일정한 것이 특징이다. 이처럼 토사 슬러리의 이송과정에 많은 흐름 경계층이 존재한다. 대표적으로 구분하였지만 명확히 체계를 구분하는 것은 불가능하다.
전단유속과 하상의 전단응력은 실제로 침전물을 동반하고 이송하는 힘을 의미한다. 임계 값(임계 경계 전단 응력 또는 임계 침식 속도)은 침전물 이동이 불가능한 한계 조건으로 항복이 시작되는 전단응력(운동 한계점)의 값이다. 하상의 거칠기의 경우 바닥에 쌓인 퇴적물 입자의 크기에 따라 값이 산정되며, 유속분포의 변화와 퇴적물의 이송 능력에 영향을 미친다.
결과적으로 침강속도는 유체 역학적 힘이 입자에 가해짐으로써 입자의 이동에 영향을 미치고 이에 따른 경계층의 형성을 의미한다. 난류 흐름조건에서 유속이 증가하면 입자에 작용하는 힘의 크기가 임계값을 초과하면서 관 저면의 퇴적된 입자가 다시 일정한 속도로 동반된다. 임계 조건에서의 관내 전단 응력이 침전물 이동 속도를 지배함을 의미함으로써 유속 의존성이 매우 강하다는 것을 나타낸다.
3.2 우수관로 내 유동 경계층 검토
본 연구에서는 관수로 내 토사의 이동에 대해 유동흐름의 경계층을 검토하고자 한다. 혼탁된 유체의 이동과정에서 경계층을 결정할 수 있는 주요 특성은 흐름과정에서 토사 슬러리의 관 내부 수직 분포와 바닥 경계면에서의 전단유속이다. 전단유속은 관 저면으로 침전되는 토사 슬러리를 재부상과정에 의해 부유시키는 역할을 하며, 이때의 체적농도의 특성에 기인하여 원활한 이송 또는 재 침강과정을 유도한다. 따라서 침전층의 영향으로 인한 전단유속을 검토함으로써 관 저면에서의 전단유속의 변화량을 측정하였다. 결국 이러한 전단유속의 변화량은 운동에너지의 변화량을 의미한다.
전단유속의 경우 일정 단면에서 흐르는 유속분포를 가정할 때, 아래의 식으로 근사치 추정이 가능하다. 다만, 하상 근처에서 상대적으로 얇은 층(z/H<0.2)을 설명하는데 주로 쓰이는 벽의 법칙(Law of the Wall)을 토대로 산정이 가능하다.
여기서 u는 하상을 기준으로 상부 거리 z에 대한 평균유속, z 0 u *
결국 관 저면의 경계에서 일어나는 연속적인 속도의 증가에 따른 변화량을 반영함으로써 이송⋅침전과정 동안의 전단유속과 전단응력을 산정하게 된다. 다만 임의 산정지점에서 y가 불연속적으로 증가한다면 위 Eq. (5)에서 언급된 ∂u/∂y는 무한히 증가한다. 결국 유속의 변화가 잦아드는 지점에서의 전단응력이 한계전단응력으로 정의할 수 있다. 위의 식에서 비례상수 μ를 유체의 점성계수(viscosity)라고 부르는데, 전단력은 유체가 관 내부를 흘러갈 때 유체의 유속으로 인한 관 벽에 생긴 마찰저항력을 관 단면적에 대해 고려한 값이다. 결국 이 마찰저항력은 유속에 비례하고 관 단면적에 반비례하는데 이를 전단율 또는 전단응력의 개념으로 사용하고자 표현한 것이다.
다음의 Fig. 5의 결과는 1.0 mm의 토사 입자를 고려하였을 때 관 벽에서의 전단유속의 분포를 대표 나열하였다. 해석과정동안 관 저면에서 수렴된 변화 값을 갖는 것을 볼 수 있으며 관 벽으로부터 일정 거리 이격된 이상에서는 전단유속의 변화가 없어 결국 전단마찰이 관 저면에서 강하게 작용하지만 관 중심부까지는 미치지 않음을 가리킨다.
본 연구에서는 수치해석 결과를 활용하여 관 저면의 얇은 층에 대한 전단유속을 산정하였다. 전단유속(u * v s
여기서 ρ s - ρ / ρ ρ s d p C D ρ S - ρ / ρ 1 / 3 g d p / C D
산정된 결과를 살펴보면 기존의 Raudkivi(1998)가 제안한 결과의 경우 약 2.0 보다 큰 값일 때 대부분의 퇴적물이 관 저면에 침전된 상태를 나타내고 있다. 즉 이송현상이 발생하지 않고 저면에 고착되었음을 가리킨다. 하지만 본 연구에서 고려한 우수관거의 경우 평균적으로 1.0 미만의 값에서 퇴적의 경향성이 두드려지게 나타났다. 이후에 급격한 전단유속의 변화된 값이 0.2∼0.9를 나타낸 것으로 보아 해당 층이 입자의 도약에 따른 이송⋅침전이 반복되는 층으로 구성되었으며 부분적으로 부유 된 퇴적물 이동 현상이 나타났다. 0.2의 값 미만의 값들 사이에는 전단유속의 변화가 잦아들어 결국 현탁상태로 부유하여 흘러가는 것으로 나타났다. 결과를 종합하여 관 내부 전단유속을 활용하여 경계층을 형성에 대한 기준을 Table 6에 수식적으로 표현하였다. 결국 우수관로 내부 슬러리의 흐름특성에 대하여 전단유속(u * v s
본 연구에서 제안된 결과는 우수관거로 유입되는 다양한 유입 조건에 대한 해석결과를 토대로 추정하였다. 입자의 크기를 주요 매개변수로 하여 침강속도를 반영하였기 때문에 관 내 침전층에서의 침전물 이동 현상을 나열하였다. 본 연구에서 제안한 경계층 검토 연구의 주된 이점은 토사 슬러리가 낮은 체적농도 조건에서도 이송을 위한 경계층의 두께를 예측할 수 있을 뿐만 아니라 다량의 토사 슬러리가 이송되는 높은 체적농도의 조건에서도 이송과정에 따른 경계층을 예측할 수 있다는 것이다. 더불어 기존 제안식들이 갖는 하천에서의 값과 달리 우수관거에 적용할 수 있다는 특징을 갖는다. 이후에도 개선된 연구를 통해 본 연구에서 제안한 경계층 모델과 실제 우수관 내부에서 나타나는 현상을 계측을 통해 오차를 비교 할 경우 보다 합리적인 연구방법이 될 것으로 판단된다.
4. 결론 및 고찰
본 연구에서는 이송⋅침전 현상이 발생하는 관거 내부에서 난류구조를 파악하고 토사 슬러리의 이송⋅침전의 영향을 검토하기 위해 수치해석을 실시하였다. 난류모델을 고려하였으며 유속, 체적농도의 변화에 대한 전단유속을 검토하였다. 이를 토대로 우수관로 내부 슬러리 이송에 대한 난류 전단경계층을 구분하는 새로운 방법론을 제시하였다. 먼저, 퇴적특성을 갖는 관수로에서 유사의 체적농도에 의한 항력으로 유속분포가 전반적으로 약화되었으며 침전에 따른 관 내 유속경사가 커짐을 확인하였다. 둘째, 전단유속의 경우 하상 경계면 부근에서 값을 산정할 경우 큰 값을 형성하여 이를 전반적인 기울기의 형태를 반영하여 추정하였다. 둘째, 경계층의 구분은 관 내 유속장 검토를 바탕으로 전단유속과 입자크기별 침강속도간의 관계로 표현하였다. 셋째, 관 내 평균유속과 체적농도에 대한 조건별 결과가 유사한 경향을 따르는 점에 착안하여, 침강속도와 전단유속의 상관관계를 평균화하였다. 본 연구에서 결정한 전단경계층의 매개변수는 기존 Raudkivi(1998)의 연구에서 제시된 매개변수의 범위에 포함된다. 즉, 본 연구에서 제안한 기법이 신뢰할 수 있는 결과를 나타내는 것을 확인하였다. 그러나 본 연구에서 제안한 결과는 단지 강우 초반에 다량의 토사 슬러리가 유입되는 상황을 고려한 결과에 불과하다는 한계를 갖는다. 따라서 보다 다양한 호우사상에 대한 적용을 통해 제안된 방법론의 검증이 필요할 것으로 판단된다. 나아가 본 연구에서 제안한 결과들의 평가를 위해서는 실제 관로시스템의 자료와 강우-유출 모델과 비교하여 검증되어야 한다.
PDF Links
PubReader
ePub Link
Full text via DOI
Download Citation
Print
