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

최근 들어 기후 변화로 인한 이상 강우 현상과 기록적인 집중 강우로 탁수 발생의 빈도와 강도가 증가하고, 탁수의 장기간 하류 방류로 인해 경제적인 피해가 발생할 뿐만 아니라 환경적으로도 중요한 문제로 대두되고 있다. 저수지로 유입된 탁수의 시공간적 분포 예측 결과의 신뢰도는 탁수 밀도류 해석의 정확도에 가장 큰 영향을 받으므로 유입 하천수의 밀도 거동에 대한 정확한 해석이 필요하다. 따라서 탁수의 거동해석은 수리 및 수질의 연동 모의(Coupled modeling)가 필수적이다. 본 연구의 목적은 수리 및 수질의 연동 해석이 가능한 2차원 횡 방향 평균 모형(CE-QUAL-W2)을 소양호에 구축하고, 기록적인 강우 사상을 보인 2006년 수문조건에서의 모형의 적용성과 한계점을 평가하는데 있다. 구축한 소양호 탁수 예측 모델을 2006년 탁수사상에 적용한 결과, 홍수 사상으로 이중 첨두 유량이 발생한 후 저수지 내에서 중층 밀도류와 바닥밀도류의 두 가지 밀도류 형태가 혼재하여 진행되는 현상이 적절히 모의되었다. 하지만, 저수지내에서 측정한 탁수 수괴(plume)의 시공간적 수직분포(profiles)는 모의결과와 다소 편차를 보였으며, 특히 저수지 바닥 인근에서 모델은 실측값을 과대평가 하였다. 오차의 원인은 모델 입력 자료의 불확실성과 수치모델의 불확실성으로 판단된다. 따라서 향후 탁수 거동 모델링과정의 불확실성을 최소화하고 모의결과의 신뢰도를 개선하기 위해서 극한탁수 사상에서 발생하는 수리학적 현상과 SS 입자 동력학에 대한 이론적 연구와 입력 자료의 정확한 산정을 위한 현장 실험이 필요하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.394-394
• /
• 2020

최근 빈번히 발생하고 있는 슈퍼 태풍과 집중호우로 인해 크고 작은 침수피해가 발생하고 있다. 우리나라는 4대강 사업을 통해 주요 국가하천에 대한 정비를 마친 바 있으나 이후 지속적으로 변화하는 하천 환경에 대한 홍수예측 모형의 반영은 미비한 실정이었다. 따라서 본 연구에서는 낙동강 본류 및 지류의 최신 단면자료를 수집하여 수리학적 홍수예측모형에 반영하고자 하였다. 또한, 기존의 모형에 비해 정확도를 개선할 수 있는 방안을 모색하여 적용성을 검증하고 이를 반영한 수리학적 홍수예측 모형을 제시하고자 한다. 본 연구에서는 낙동강 본류 및 주요 지류에 대한 최신 횡단면 측량자료를 활용하여 1차원 수리학적 홍수예측 모형을 구축하고 2012년 태풍 산바 사상에 대한 검보정을 실시하였다. 대상구간은 안동조정지댐으로부터 낙동강 하구둑 하류 8km지점이며, 상류단 경계조건은 안동조정지댐 방류량을 입력하고 하류단 경계조건은 가덕도 조위관측소의 조위를 활용하였다. 또한, 반변천, 내성천, 위천, 감천, 금호강, 황강, 남강 등 7개 지류에 대한 하도를 하도추적이 가능한 네트워크 모형으로 구축함으로써 지류에 대한 홍수예보에 활용할 수 있도록 하였다. 낙동강 본류의 하도 길이는 340km, 824개의 단면으로 구성하였으며, 지류를 포함한 전체 하도 길이는 572km, 1,570개의 단면으로 구성하였다. 또한, 낙동강 본류에 위치한 8개의 다기능보와 지류에 위치한 횡단구조물의 반영을 위해 다기능보의 제원 및 하천기본계획을 참고하여 내부경계조건으로 활용하였다. 본 연구에서 구축된 낙동강 본류 및 지류의 수리학적 홍수예측 모형을 낙동강 유역에 대한 홍수주의보 및 홍수경보 등 홍수특보 발령 업무에 활용함으로써 정확한 홍수예보가 가능하도록 하였다.

• Journal of Korea Water Resources Association
• /
• v.35 no.5
• /
• pp.485-500
• /
• 2002

Hydraulic model was developed to analyze the complex flow due to channel structures, tide, and tributaries in the lower Han river and Imjin river. DWOPER-2K model which can automatically process the data transformation in the model was developed as the 1-D hydraulic routing model. Observed data in tidal zone and the recent channel geometry data were collected for hydraulic model. And the flow over the Jamsil and Singok submerged weir was analyzed properly and roughness coefficient was optimized to each regions and each discharges. By the results of verification of the model, the model developed in this study may contribute to improvement of the accuracy of flood forecasting and channel management because this model can efficiently and properly analyze the various kind of flow occurred in the region of the lower Han river and Imjin river.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.587-587
• /
• 2012

현재 대부분 하천관리의 흐름해석에서 주가 되었던 1차원 분석은 하천 단면에 따른 횡적인 수면차, 유속분포를 분석할 수 없다는 단점을 가지고 있으며, 특히 유량 및 유속이 급속도로 늘어나는 홍수시에는 그 오차가 더욱 커질 수 있다. 반면에 2차원 모형의 흐름해석은 사행하천의 흐름 특성과 만곡부에서의 종 횡방향 수면경사 및 양안의 수면차와 합류지점의 횡방향 흐름 등의 영향을 고려할 수 있으며 1차원 해석과는 달리 전 단면에 걸쳐 유속 및 수위 분포를 나타낼 수 있어 실제흐름에 가까운 수리량을 얻을 수 있다. 본 연구를 위해 적용된 해석모형인 SMS는 미국 Brigham Young 대학의 환경모형연구실과 미공병단(USACE)의 수로실험국(WES) 등에서 개발한 프로그램으로서 RMA2, RMA4, SED2D 모형 등으로 구성되어 있다. 각각의 모형은 수리 동역학적 해석, 오염물 이송확산 해석, 유사의 이송 및 퇴적 해석이 가능하며 이 중 RMA2를 이용한 2차원 흐름해석을 통하여 보다 적합한 하천관리에 이용가능하도록 하고자 한다. 연구대상 지역은 제주도 한천 하류부로서 제주도 하천 특성상 평상시 건천의 상태를 이루고 있으나 태풍 및 집중호우시 홍수유출이 발생하여 수위가 급격하게 상승하는 양상을 보인다. 대표적인 예로 태풍 '나리'시 최대 일강우량 420mm로 인한 인근 지역에 0.5 ~ 1.5m의 침수흔적을 보이고 있다. 본 연구에서는 2011년 최고수위를 기록한 태풍 '무이파'를 대상으로 하였으며 대상지역의 Kalesto를 이용한 수위-유속 자료를 이용하여 산출된 유량을 경계조건으로 사용하였고 격자망 형성을 위한 지형 데이터는 지형도 및 측량자료를 이용하여 구축하였다. 사용된 대표적인 매개변수는 하상의 조도계수를 나타내는 Manning의 n값과 유체의 밀도, 속도구배, 구조 등 여러 가지의 유체조건에 따라 변하는 성질인 와점성계수(eddy viscosity)로 요약할 수 있으며 Manning의 n값은 하천설계기준에 따른 하천기본계획의 조도계수를 사용하였고 와점성계수는 적합한 흐름 분포를 결정하기 위해 흐름이 안정될 때까지 변화시켜 해석을 시행하였다. 해석결과 만곡부에서는 급한 흐름을 보이고 있으며 최대하폭 구간에서는 완만한 흐름이 나타나 사행하천의 흐름특성과 횡적인 하천단면에 따른 변화, 하상고 차이로 인한 유속분포를 확인할 수 있으며 이는 보다 유용한 하천관리에 이용가능할 것으로 사료된다.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.59-59
• /
• 2011

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

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

본 연구는 복단면구조의 수로에서 이차류 및 흐름특성을 비교, 관측하는 것이며 이러한 흐름특성을 연구하기 위해 폭 1.2 m, 길이 14 m, 높이 0.6 m인 실험수로를 이용하여 수심비 0.25 및 0.40 조건에서 0.05, 0.10 CMS의 유량으로 총 4가지 조건으로 실험을 진행하였고 이를 비교하기 위해 3차원 수치해석모형인 FLOW-3D를 이용하여 수치해석결과와 수리실험결과를 비교하였다. 주수로 및 홍수터의 수심과 유량의 조건이 각각 다른 4가지 조건 별 변화를 수리실험과 수치해석을 통해 분석한 결과, Shino & Knight(1991)에서의 언급한 내용과 부분적으로는 동일한 결과를 보였으나, 수치해석 부분에서의 결과와는 차이를 보였다. 일반적으로 복단면구조에서 발생하는 이차류는 수위 및 수심별 유속의 차이에서 발생하는 연직방향의 유속의 영향과 복단면과 주수로의 구조적 차이에서 발생하는 횡방향 유속의 영향과 같은 복합적인 작용으로인해 발생하게 되지만 수치해석의 결과에서는 수위차이에 의한 연직방향의 유속이 횡방향 유속보다 지배적인 결과를 보여 본 연구 중 고유속의 조건에서는 이차류의 형성이 수리실험 및 기존 연구결과들과 차이를 보였다.

• Journal of Korea Water Resources Association
• /
• v.42 no.4
• /
• pp.309-317
• /
• 2009

This study analyzed a hydrodynamic behavior using the EFDC model (Environmental Fluid Dynamics Code) in the downstream of the Nakdong River in the case of a storm surge and a localized torrential rainfall caused by a major typhoon, and the sea level rise caused by global warming. The study area is selected Gaduk island with the lower boundary and Jindong with the upper boundary, to investigate the total river hydrodynamic behavior including the estuary. In order to verify this numerical model, the calculated results was compared with the observed stage at each gauging point in case of the storm rainfall in 2003 and 2006. From the results, it was shown that the numerical model(EFDC) has high accuracy and is useful in simulating more various cases.

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

하천을 횡단하는 구조물은 흐름을 지체시켜 상 하류에 급격한 수위차를 발생시키며, 이에 동반되는 국부유속의 증가는 제방 및 하상을 세굴시킨다. 제방 및 하상 세굴은 하천횡단구조물의 구조적 안정성에 영향을 주어 불안정한 상태에 이르게 할 수 있다. 이와 같은 현상을 방지하기 위하여 하상유지시설 및 제방 보호공 등의 대책을 하천설계기준에 제시하고 있다. 하천설계기준에는 하천횡단구조물 보호시설의 종류, 설치길이 등이 명시되어 있으나, 보호공 제원 선정에 대한 구체적인 설계공식을 제시하지 않아 보완이 필요한 상황이다. 보호공 제원 선정 설계기준은 설계유속을 1차원 단면평균 유속으로 활용하는 국외의 설계공식을 활용할 수 있다. 그러나, 하천횡단구조물 주변의 흐름은 매우 복잡하므로 국외 설계공식을 그대로 적용하기에 앞서서 이에 대한 검토가 필요하다. 본 연구에서는 하천횡단구조물 주변의 흐름을 2차원 수치모형으로 재현하고 그 결과를 1차원 수치모의 결과 및 수리실험 관측결과와 비교하였다. 실험수로에 각기 다른 제원의 보호공을 설치하고 관측한 임계상태의 수리상황을 1, 2차원 수치모형으로 재현하였다. 각각의 모형에서 도출된 수리상황 결과를 비교 분석하여 보호공 제원 결정 공식을 검토하였다. 또한, 실제 하천구간에 임의의 보와 수제를 설치하고, 계획홍수량에서의 흐름특성을 비교하였다. 1, 2차원 모형에 의하여 도출된 관계를 분석함으로써, 하천횡단구조물 설계시 적용할 수 있는 보호공 제원 결정 설계공식을 제안하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.2
• /
• pp.166-177
• /
• 2006

The accuracy impact of using high-order Boussinesq-type model as compared to the typical order model is examined in this paper. The multi-layer model developed by Lynett and Liu(2004a) is used for simulating of wave breaking over a shelf region. The nonlinearity of the waves tested, ${k_0}{A_0}$, ranges from 0.029 to 0.180. The overall agreement between the two-layer model and the hydraulic experiments are quite good. The one-layer model overshoals the wave near the breakpoint, while the two-layer model shoals at a rate more consistent with the experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.541-545
• /
• 2012

이 연구에서는 SEEP2D 모형을 사용하여 하부에 투수지반을 갖는 제방 내부의 침투흐름에 대한 수치모의를 수행하였다. 제방의 재료는 4대강사업 낙동강 제2하구둑 건설 현장에서 채취한 시료를 사용하였으며, 입도분석과 들밀도시험을 통해 구한 제방의 건조단위중량과 다짐도 및 제방의 투수계수는 각각 $1.372g/cm^3$, 93%, 1.35 m/day이다. 투수지반층과 1:3 비탈경사를 가지는 제방모형을 제작하여 연직 이차원 정상류 침투해석 모형을 해석하는 SEEP2D 프로그램을 사용하여 나온 결과 값을 수리모형 실험 결과 값과 비교하였다. 모의조건은 0.45 m, 0.50 m, 0.55 m, 0.60 m의 4가지 수위(각각 Case1~4)조건에 대하여 수리모형실험과 동일한 조건을 수치모형에 적용하였다. 수치모형의 결과는 실험결과와 비교적 일치하였지만 제방의 비탈사면이 시작하는 부분에서는 수치모형의 결과 값보다 실험결과 값이 작은 것으로 확인되었다. 각각의 조건별로 수치해와 실험 결과값을 비교해본 결과 대체적으로 비슷한 양상을 보였지만 제방의 비탈사면이 시작하는 지점부터 수치모형의 결과 값보다 실험결과 값이 작은 것으로 확인되었다. Case3의 경우 4.35 m 지점부터 6.00 m 지점까지 유출이 발생하였으며 수치모형의 결과값과 실험결과 값의 차이가 가장 작게 나타났다. 사면유출 길이는 Case4에서 가장 길게 나타났으며 최대 4.10 cm로 발생하였다.