• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.1955-1959
• /
• 2009

본 연구에서는 토양층에서 지하수 수리경사와 주 흐름 방향을 산정하여 변동특성을 분석하고, 그 원인에 대해 고찰하였다. 연구부지는 부경대학교 환경연구동 앞 잔디밭이며, 2008년 10월 18일 ${\sim}$ 11월 1일의 기간 중 12일(10월 21, 28일은 관측 못함) 동안 24회 관측되었다. 관측된 지하수공은 9개공이며, 4 m ${\times}$ 4 m의 정방형으로 설계되었다. 지하수위에 의한 수리경사와 주 흐름 방향의 산정은 총 9가지 형태의 평면으로 가정하여 산정하였다. 수리경사와 주 흐름 방향을 산정하기 위한 평면은 모두 9가지 형태로 가정하였다. 첫 번째는 9개공 모두를 이용한 정방형 평면(TW), 두 번째는 6개공을 이용한 4가지 경우의 삼각형 평면(T1, T2, T3 및 T4), 세 번째는 6개공을 이용한 4가지 경우의 직사각형 평면(R1, R2, R3 및 R4)으로 가정하여 각각의 평면에서 수리경사와 주 흐름 방향을 산정하였다. 이상에서 가정된 9가지 평면에서 산정된 평균 수리경사는 0.0094(R2) ${\sim}$ 0.0123(T2)의 범위로 나타났으며, 9개공 모두를 이용한 평면(TW)에서의 평균 수리경사는 0.0109로 산정되었다. 관측기간 동안 T1과 T3 평면의 평균 주 흐름 방향은 각각 -84.9818과 -86.2487로서, 다른 7개 평면의 평균 주 흐름 방향이 79.7045(T4) ${\sim}$ 85.8405(TW)의 범위로 나타난 것과는 반대의 흐름 방향을 보였다. 이러한 결과는, 본 연구에 적용된 수리경사와 주 흐름 방향의 산정기법(Devlin, 2003)이 삼각형 형태의 평면에는 부적합할 수 있음을 보여주는 것이다. 관측기간 중 10월 12일에 내린 소량의 강우(12 mm)에 의한 수리경사와 주 흐름 방향의 변동은 크게 나타나지 않았다. 본 연구에서는 정방형, 삼각형 및 직사각형의 평면 토양층에서 지하수의 수리경사와 주 흐름 방향을 산정하여 전체 및 국부적인 부지에서의 차이를 파악하고, 평면의 형태에 따른 수리경사와 주 흐름 방향의 차이 또한 확인할 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.15-15
• /
• 2022

우리나라 홍수 발생은 강수량이 집중되는 여름철에 집중되어 있어 홍수피해 방지에 주의가 필요하다. 보와 제방 등 이러한 홍수 피해를 대비하기 위한 구조물에 설계를 위해서는 하천의 유량 조사가 필수적으로 요구된다. 하지만 홍수기 직접적인 유량 조사는 안전상의 이유로 거의 이루어지지지 않고 있으며, 수위를 측정하여 수위-유량 관계를 만들어 유량을 측정하고 있다. 그러나 중소 규모의 하천의 경우 하도 경사가 급해 사류가 발생하거나 하도 단면이 급변하는 경우가 있어 수위-유량 관계를 그대로 적용하기 어려운 문제가 있다. 따라서 이러한 문제를 해결하기 위해 하천에 진입하지 않고 유속을 측정 할 수 있는 영상유속계와 같이 흐름 영상을 사용하여 유속을 측정하는 방법들이 개발되었다. 영상유속계의 측정 방법중 Spatio Temporal Image(시공간 영상)을 사용하는 방법은 일정시간의 시간평균 유속을 산정할 수 있고 한 장의 시공간 영상을 분석하기 때문에 유속 산정에 걸리는 시간이 작은 장점이 있지만 영상 내 흐름 방향을 정확히 설정하지 못하면 오차가 생길 수 있는 문제가 있어 주 흐름 방향을 정확히 탐색할 필요가 있다. 따라서 본 연구에서는 영상에서 시공간 체적을 만들어 주흐름 방향을 찾아내는 기법들의 장단점을 비교하고 안동 하천실험센터의 실규모 하천수로에 적용하여 결과를 비교하고 적용성을 평가하였다. 이를 위해 하천수로에 추적자를 살포하여 영상으로 녹화하였으며 녹화된 영상을 자기상관법과 시간적분법을 적용하였으며, 이를 통해 주 흐름 방향을 판별하였다. 또한 두 방법을 통해 결정된 주흐름 방향을 적용하여 시공간 영상을 제작하고 이를 이용하여 유속을 산정하여 비교하였으며, 주흐름 방향을 산정하는데 생기는 오차가 유속 계산에 얼마만큼의 영향을 끼치는지 분석하였다. 이러한 실험을 통해 하천에서 시공간 영상을 활용한 표면영상 유속계측 방법을 활용하는데 있어 도움이 될 것으로 기대된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.303-303
• /
• 2023

실용적인 표면영상유속계를 만들기 위해서는 적절한 하드웨어와 소프트웨어로 시스템을 구성해야 한다. 이를 위해서 본 연구에서는 하드웨어로 CCTV를 선택하고, 초음파 수위계를 이용하여 수위를 지속적으로 읽어들이도록 구성하였다. 한편, 소프트웨어적으로는 11변수 투영법을 적용하여 변화하는 수위에 따라 정확한 측정점을 재구성하도록 하고, 아울러 각 측정점에서 주흐름방향으로 정확한 시공간영상을 작성하고, 이 시공간영상(spatio-temporal image)들을 분석하였다. 그결과, 5분 간격으로 촬영된 1분 길이의 영상을 지속적으로 촬영하고 분석하여 유량을 산정하는 표면영상유속계측 시스템을 구축하였다. 본고에서는 이러한 소프트웨어 개선방향중 하나인 주흐름방향의 시공간영상 작성법을 소개한다. 먼저, 11변수 투영법을 이용하여 하천의 표면영상에 대한 좌표변환계수를 산정하였다. 그리고 이 좌표변환계수를 이용하여, 하천의 수위변화에 따라 표면영상내의 측정점이 적절히 수정될 수 있도록 하였다. 그 다음 이 측정점에서 측정횡단면과 수직이 되는 방향을 선정하고, 이 방향이 영상내에서 하천 측정횡단면과 수직인 방향, 즉 주흐름방향이 되도록 하였다. 촬영된 1분간의 동영상의 각 측정점 위치에서 잘라낸 시공간체적(spatio-temporal image volume)에서 주흐름방향의 시공간영상을 잘라내고 이를 상호상관법이나 고속푸리에변환을 이용하여 분석하였다. 이 때 만들어진 시공간영상은 주흐름방향과 정확하게 일치하여, 기존의 표면영상유속계의 문제이던, 일부 측정점의 유속벡터가 주흐름방향과 일치하지 않던 문제를 해결할 수 있었다. 개발된 방법으로 표면영상유속계를 제작하여 인수천에 시험 설치하고 호우 사상에 대해 검토한 결과 정확하고 신속하며 연속적인 유량측정이 가능하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.136-136
• /
• 2019

원형 실린더를 주변 흐름에 관한 연구는 오랜 기간 유체역학 전 영역에서 모형실험이나 수치모형으로 광범위하게 연구되었다. 이 흐름은 하천의 교각이나, 바다의 시추선과 같은 수공구조물 주변에서 관측된다. 난류와 와류가 공존하는 복잡한 특성 때문에, 이 흐름은 수공학에서 유사이송, 세굴, 오염물 확산 등에 영향을 준다. 본 연구는 실험실 수로에 설치된 원형 실린더(D=9cm) 후방의 근접 와류 구간에서(x/D<5) 유속을 ADV로 측정한 후, 난류 특성을 Power Spectral Estimation(PSE)와 Singular Spectral Analysis(SSA) 방법으로 연구하였다. PSE는 샘플 스펙트럼의 한계를 보완하고자 자료를 분할하고, window 함수를 적용하여 ensemble 평균을 구하는 경험적 방법이다. PSE를 이용하여 스펙트럼을 계산한 결과, 주 흐름 및 횡방향 흐름은 Inertial subrange에서 Kolmogorov의 가정과 일치하는 추세를 보였다. 그러나 수심방향 흐름의 스펙트럼은 -5/3보다 빠르게 감소하는 추세를 보였다. Inertial subrange 스펙트럼에서 난류 에너지 소산율은 원형 실린더에서 멀어짐에 따라 감소하는 추세를 보였고, 주 흐름방향과 횡방향 흐름은 비슷한 크기를 보였다. 난류 에너지 소산율과 동점성계수를 이용하여 Kolmogorov 길이, 유속, 시간 스케일을 계산했다. 난류의 운동에너지를 계산하기 위해 Triple decomposition 방법 중 하나인 SSA를 적용하였다. SSA는 유속행렬을 이용하여 고윳값과 고유벡터를 계산하고, 유속에서 기여도가 큰 부분을 추출하는 방법이다. SSA를 통해 실린더 후방 흐름에서 와류 성분과 난류 성분을 나누었다. 횡방향 흐름은 강한 와류로 큰 기여도를 갖는 고유벡터가 나타났지만, 주 흐름과 수심방향 흐름은 상대적으로 낮은 기여도를 갖는 고유벡터가 나타났다. 와류를 제외한 흐름에서 난류 운동에너지는 실린더와 멀어짐에 따라 감소하고, 흐름 중앙에서(y/D=0) 가장 큰 값을 보였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.434-438
• /
• 2005

본 연구에서는 개수로 흐름에서 전단속도의 분포를 조사하였다. 이를 위해 직사각형 개수로 흐름과 복단면 개수로 흐름에 대한 DNS 자료를 이용하였으며, 계산된 전단속도는 기존의 실험자료 및 수치해석결과와 비교$\cdot$분석되었다. 평균흐름장에서는 직사각형 개수로 흐름의 경우, 자유수면 부근에서 내측이차흐름이 발생하는 것을 볼 수 있었으며, 복단면 개수로 흐름에서는 주수로와 홍수터의 접합부 부근에서 쌍와가 생성되는 것으로 나타났다. 이차흐름의 최대크기는 체적평균된 주흐름방향 유속의 약 $5\%$로서, 직사각형 개수로 흐름에서 발생되는 것보다 더 강하다는 것을 알 수 있었다. 이러한 이차흐름의 영향으로 주흐름방향 평균유속의 등속선이 편향된 것으로 나타났다. 본 DNS에서 전단속도가 최대값을 갖는 지점은 직사각형 개수로 흐름의 경우 자유수면이며, 복단면 개수로 흐름의 경우 접합부인 것으로 나타났다. 이는 기존 LES 결과와 동일하지만, 실험연구와는 다른 결과이다. 접합부에서 이러한 차이를 보이는 이유는 실험에서 마찰속도를 산정하기 위해 사용하는 방법의 부적합성에 있는 것으로 나타났다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.883-888
• /
• 2007

일반적인 하천의 흐름방향으로 발생하는 주흐름(primary flow)에 중첩하여 주흐름 방향의 수직단면에 이차류(secondary flow)가 발생하게 되며 이러한 이차류의 발달은 투입된 오염물질의 횡혼합을 증대시킨다. 오염물질의 혼합은 이송(advection)과 확산(diffusion) 또는 분산(dispersion)의 과정으로 설명되며 본 연구에서는 수로전체의 혼합과정을 설명하기 위해서 이송 확산 방정식을 적용하였다. 본 연구에서는 실험수조를 $150^{\circ}$의 중심각을 갖는 S자 형태의 만곡수로를 제작하여 유량조건은 15, 30, $60\;{\ell}l/sec$의 세 가지 경우로, 수심은 15, 20, 30, 40 cm의 경우로 총 12 케이스의 실험을 수행하였다. 유속장의 측정은 Sontek사의 3차원 micro-ADV(Acoustic Doppler Velocimeter)를 이용하였다. 오염물질 확산실험은 소금물 용액에 주변수와의 밀도차를 없애기 위해서 메탄올 용액을 첨가하여 추적자로 이용하여 농도장의 분석을 일본 KENEK사의 전기전도도계(conductivity meter)와 Gartner사의 DAS(data acquisition system)를 이용하여 횡방향 유속장의 분포와 오염운의 거동을 비교하여 다음과 같은 결론을 얻게 되었다. 주 흐름은 직선구간에서는 중앙에서 최대 유속을 나타내며, 좌우대칭적인 유속분포의 모습을 보이고, 만곡부에서는 수로안쪽을 따라 최대유속이 발생하였다. 수로의 직선구간에서는 최대유속이 발생하는 즉, 중앙에서의 오염물질의 분산이 가장 활발하게 이뤄졌으며 농도의 퍼짐형상인 오염운 역시 만곡부에서는 수로만곡부의 안쪽을 따라 확산 이동함을 알 수 있었다. 만곡부 외측에서는 오염물질의 정체현상이 일시적으로 발생하며, 유속구조의 횡방향 비대칭구조로 인한 종 횡방향의 분리현상이 발생하고, 오염운의 중첩현상이 종방향으로 연속되게 나타난다. 향후 수심방향 거동을 포함한 3차원적 분석이 요구되며 이 연구결과는 2차원적 수치해석의 적용 및 분석 자료로써 이용이 가능하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.545-549
• /
• 2009

본 연구에서는 주흐름방향으로 식생 영역이 교차적으로 존재하는 개수로 흐름에 대한 3차원 수치 모의를 수행하였다. 지배방정식의 난류 폐합을 위해 $k-{\varepsilon}$ 난류모형을 이용하였다. 먼저, 하상의 일부만 식재된 부분 식생 수로를 수치모의 하고 기존의 실험 결과와 비교하였다. 그 결과 본 모형이 평균유속 분포를 매우 잘 예측하는 것으로 나타났으나, 레이놀즈응력 분포는 실험 결과에 비해 비식생영역에서는 다소 과소 산정하고 식생영역에서는 과대 산정하는 것으로 나타났다. 이는 본 모형이 등방성 모형이기 때문에 식생 경계부에서 발생되는 난류의 비등방성 효과를 정확히 예측 할 수 없기 때문인 것으로 판단된다. 또한 주흐름방향으로 식생 영역이 교차적으로 존재하는 대응 식생 수로를 수치모의하고, 계산 결과를 기존의 실험 결과와 비교하였다. 그 결과 본 모형이 대응 식생 수로에서의 유속 분포를 매우 잘 예측하는 것으로 나타났다. 또한 식생 밀도가 증가함에 따라 식생이 흐름 방향을 변화시켜 점차 만곡수로와 유사한 형태의 흐름이 형성되는 것으로 나타났다.

• Journal of Korea Water Resources Association
• /
• v.54 no.11
• /
• pp.933-942
• /
• 2021

The surface image velocimetry was developed to measure river flow velocity safely and effectively in flood season. There are a couple of methods in the surface image velocimetry. Among them the spatio-temporal image velocimetry is in the spotlight, since it can estimate mean velocity for a period of time. For the spatio-temporal image velocimetry analyzes a series of images all at once, it can reduce analyzing time so much. It, however, has a little drawback to find out the main flow direction. If the direction of spatio-temporal image does not coincide to the main flow direction, it may cause singnificant error in velocity. The present study aims to propose a new method to find out the main flow direction by using a fast Fourier transform(FFT) to a spatio-temporal (image) volume, which were constructed by accumulating the river surface images along the time direction. The method consists of two steps; the first step for finding main flow direction in space image and the second step for calculating the velocity magnitude in main flow direction in spatio-temporal image. In the first step a time-accumulated image was made from the spatio-temporal volume along the time direction. We analyzed this time-accumulated image by using FFT and figured out the main flow direction from the transformed image. Then a spatio-temporal image in main flow direction was extracted from the spatio-temporal volume. Once again, the spatio-temporal image was analyzed by FFT and velocity magnitudes were calculated from the transformed image. The proposed method was applied to a series of artificial images for error analysis. It was shown that the proposed method could analyze two-dimensional flow field with fairly good accuracy.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.168-168
• /
• 2015

자연하천의 주요한 특징 중 하나인 하천의 사행은 직선수로에서 예측되는 유속분포를 왜곡시키며 매우 복잡한 흐름구조를 형성한다. 이는 하상 경계면에서 발생하는 전단응력 분포의 변화를 야기하는데 하상 경계면에서의 전단응력은 다양한 경험적 관계에 의존하는 유사이동의 한계 소류력 산정 및 오염물질 거동해석의 분산계수 산정에 많은 영향을 미치게 된다. 물리적인 관측을 통한 하상 경계면에서의 전단응력의 관측은 다소 제한적이며 많은 비용을 요구한다. 따라서 하상 경계면에서 발생하는 전단응력의 경우 수심의 20% 이하의 연직 유속분포를 벽법칙에 적용하여 추정하는 방법이 주로 이루어지고 있다. 벽법칙을 이용한 하상 경계면의 전단응력을 계산하는 경우 대수중복층의 유속 분포 $u/u^*=(1/{\kappa})ln(zu^*/{\nu})+B$에서 무차원상수 ${\kappa}$와 B의 적절한 추정이 요구되어 진다. 일반적으로 무차원상수 ${\kappa}$와 B는 수리학적으로 매끄러운 벽면에서 대략 ${\kappa}=0.4$, B=5.5로서 경험적으로 이용되고 있다. 본 연구에서는 직선수로 및 다양한 사행수로의 3차원 흐름장 모의를 수행하여 벽법칙의 대수 중복층을 따르는 주흐름 방향 유속의 연직분포를 비교하였다. 수치모의 소프트웨어로서 Linux 기반의 OpenFOAM이 사용되었으며 모델의 검증을 위해 Chang(1971)에 의해 수행 된 사행수로에서의 유속장 관측 결과와 비교하였고 수치모의 결과가 실험 관측치와 잘 일치하는 것으로 판단되었다. 수치모의에 적용 된 사행수로의 형상은 Hey(1976)에 의해 제안 된 사행하천의 지형학적 인자들 간에 관계를 이용하여 사행도 1.03에서 2.42까지 총 7개의 사행수로 지형을 생성하였다. 사행도의 변화에 따라 만곡부 정점에서 대수중복층 구간의 주흐름 방향 유속의 연직분포를 비교한 결과, 본 연구에서 생성 된 모든 사행수로에서 대수중복층 구간의 무차원 유속 $u^+$와 무차원 거리 $z^+$가 로그 분포를 따르는 것으로 나타났으나 경험적으로 사용되었던 무차원상수 B의 경우 사행도가 증가 할수록 대수적으로 감소하는 경향이 나타났다. 본 연구에서는 이러한 관계가 무차원 상수 B값에 미치는 영향을 반영하여 수리학적으로 매끄러운 벽면에서 적용이 가능한 수정된 대수중복층 식을 제시하고자 한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.168-168
• /
• 2011

본 연구는 일산대교의 교각이 한강 하류부의 지형변화에 미치는 영향을 분석하기 위한 준비단계로서, 이차원 흐름모형을 사용하여 일산대교 교각에 의한 흐름변화를 분석하였다. 일산대교 교각의 영향을 검토하기 위해 일산대교 상류 쪽 1.2 km부터 하류 쪽 1.7 km까지 약 2.9 km 구간을 대상으로 RMA-2 모형을 사용하여 흐름을 예측하였다. 이 구간의 하상고는 EL. -5.01 m ~ EL. 2.20 m이다. 일산대교의 교각이 있을 경우에 대한 절점 개수가 7,083개, 유한요소의 개수는 2,799 개이며, 교각이 없을 경우에 대해서는 절점 개수가 7,121개, 유한요소의 개수는 2,856개이다. 일산 대교에 의한 영향 검토 결과, 교각으로부터 상하류로 200 m 정도 떨어진 영역에서는 유속 벡터의 경향에서 차이가 거의 없는 것으로 나타났다. 그러나 일산대교 근처에서는 교각의 방향에 따라 흐름 방향이 유도되고 직상류에서 흐름의 방향이 급격히 변하는 현상이 나타났다. 반면, 일산대교 하류에서는 교각의 영향을 거의 받지 않는 것으로 나타났다.