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

하천으로 유입된 오염물질의 거동을 정확하게 예측하는 것은 하천 시스템의 수질 유지관리에 매우 중요하다. 본 연구에서는 1차 감쇠율(decay rate)을 가진 비보존성 오염물의 비정상 이송해석 방정식의 해를 위해 유한체적기법이 개발되어졌다. 하천 흐름 해석을 위해 자연형 단면에서의 마름-젖음 해석이 가능한 기법이 도입되었다. 이 기법은 2차-정도의 정확성와 Courant 수가 1 까지 안정함을 보장한다. 도입된 기법은 Godnov 형의 유한체적기법을 이용하여 St. Venant 방정식들을 해석하였고 질량 및 운동량 플럭스는 Roe 형의 Riemann Solver 를 사용하여 연산하였다. 오염물의 이송 해석은 추가적인 이송-확산 방정식을 도입을 통해 기존의 St. Venant 방정식과 함께 풀려질 수 있다. 추가된 방정식과 St. Venant 식은 3×3 eigenstructure를 구성하였고 이는 2차원 흐름해석 기법과 유사하게 해석될 수 있었다. 본 연구 모형의 검증을 위해 오염물의 계속적 주입을 가정한 가상 및 실제 하천에 적용되었다. 연구된 기법은 모든 적용에서 합리적 정확도를 가지고 오염물질의 연속적인 특성을 잘 모의하고 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1742-1746
• /
• 2008

지금까지의 대부분의 2차원 수리해석 모형은 상류와 사류가 혼재된 불연속적인 천이류 흐름을 해석하기에 계산의 정확도 및 현실성에서 많은 문제를 보이고 있으며, 특히 계산과정에서 나타날 수 있는 마른하도의 처리에 있어서 많은 어려움을 겪고 있다. 본 연구의 목적은 유한체적기법을 사용하여 상류와 사류가 혼재하는 불연속적인 하천 천이류를 안정적으로 해석하기 위해 개발된 고정확도 수치모형의 자연하도 적용에 있으며, 또한 마른 하도로 전파되는 흐름 모의 및 계산과정에서 나타날 수 있는 마른하도 처리의 어려움을 해결함으로써 모형의 정확도와 안정성을 검증하여 실제 하천에서의 모형 적용성을 검토함에 있다. 이를 위해 본 연구에서는 흐름의 전파양상을 정확하게 반영할 수 있는 상류이송기법인 Godunov 기법과 관심격자의 좌우 격자 정보를 모두 사용하는 대표적 중앙차분기법인 Beam-Warming 기법의 장점을 모두 반영한 가중평균흐름률 (Weighted Average Flux) 기법을 사용하여 사각격자망의 구성을 통해 자연하도에 적용시킬 수 있는 2차원 유한체적모형을 개발하고자 하였고, 개발된 모형의 안정성, 정확도, 적용성을 검증하기 위해 직사각형 수로, 큰 사행비를 가진 만곡수로에 적용하고, 그 결과를 수리모형 실험결과와 비교, 분석하였다.

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

댐 제방 등의 붕괴로 인하여 발생하는 급격한 유량의 변화와 흐름영역의 변화로 인한 천이류 및 도수의 발생, 불규칙한 하천단면에서 갈수기 저수기의 흐름해석은 기존의 수치해법의 한계로 인하여 수리모형실험 및 경험식 또는 단면의 단순화 등에 의존하고 있는 실정이다. 본 연구에서는 자연하천에서 비선형 흐름율 계산에 불연속초기조건의 해석해인 Riemann 근사해법을 사용하여 수치적으로 안정되고 정확한 1차원 모형을 개발하고자 한다. 이를 위하여 유한체적법을 사용하였고, 수위와 유량의 계산을 위하여 요구되는 유한체적을 유출입하는 흐름율의 계산에 HLL Riemann 해법을 사용하였으며, MUSCL 기법으로 2차 정확도기법으로 확장하였다. Riemann 해법을 통하여 계산된 비선형의 흐름율과 보존 특성을 만족시켜줄 수 있는 하상 및 하폭변화로 인한 생성항을 처리하는 기법을 제안함으로서 새로운 1차원 수치해석모형을 개발하였다. 개발된 모형의 실제하천의 적용성을 확인하기 위하여 하상과 하폭이 변화하는 부정류 흐름에 적용하여 모형의 적용성 및 정확성을 검증하였다.

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

천수방정식을 사용하는 초기 수치모형은 프로드수($F_4$)가 변화하는 흐름 즉, 상류방향과 하류방향으로 전파하는 홍수파를 동시에 해석하기 위해 중앙 차분기법이 필요한 상류(sub-critical flow)와 흐름방향에 따른 상류이송(upwinding)기법이 필요한 사류(super-critical flow)가 나타나는 흐름해석에서 어려움이 있었다. 하지만, 근사 Riemann 해법의 등장으로 흐름방향에 관계없이 특성선을 따라 정확한 상향가중기법의 적용이 가능하게 되어, 천수방정식을 지배방정식으로 하는 수치모형이 더욱 실용적으로 적용될 수 있도록 하였다. 따라서, 현재 근사 Riemann 해법은 Godunov 형 유한체적 기법, 불연속 Galerkin 혹은 Petrov-Galerkin 유한요소기법 그리고 Boussinesq 기법에도 적용되고 있으며, 특히 Godunov 형 유한체적기법과 결합한 근사 Riemann 해법은 댐 붕괴, 하천 범람 그리고 도시 및 해안지역 침수에 이르기까지 여러 가지 문제에 폭넓게 적용되고 있다. 지금까지 홍수 모델링에 적용된 Godunov형 유한체적모형은 정형 사각격자나 비정형 삼각격자 중에서 한가지의 격자 종류만을 적용한 연구가 주로 수행되었으며, 유한요소모형과 같이 이 두 가지 격자를 동시에 적용한 연구는 거의 이루어지지 않고 있다. 일반적으로, 삼각격자는 사각격자와 는 달리 연구유역의 경계나 지형이 복잡한 경우에도 큰 노력없이 격자의 생성이 가능하나, 격자와 노드의 수가 사각격자보다 많아 계산시간이 많이 소요되는 단점이 있다. 반면, 사각격자는 하천과 같이 선형으로 변하는 지형에 대해서는 표현하기가 용이하며 계산시간의 효율성도 뛰어나다. 본 연구에서는 하천, 도시 그리고 해안지역에서의 효율적이고 정확한 홍수 모델링을 위해 삼각 및 사각격자 그리고 이 두 격자를 동시에 고려한 하이브리드 격자의 적용이 가능한 Godunov형 2차원 유한체적 모형을 개발하였다. 그리고 개발모형을 정확해가 있는 댐 붕괴 문제, 실측치가 존재하는 실험하도 및 실제하도에 삼각, 사각 그리고 혼합격자를 생성하여 모의를 수행하고, 각 적용 격자에 따른 정확성과 효율성 및 장점과 단점을 연구하였다.

• 한국방재학회:학술대회논문집
• /
• 2009.02b
• /
• pp.152.1-152.1
• /
• 2009

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

국내외 발생하는 재난 중 70% 이상이 물과 관련된 재해로 분류되며, 집중호우와 태풍으로 인한 하천범람 및 내수침수 등으로 많은 피해를 발생시키고 있다. 특히 최근 발생하는 피해 양상은 과거에 발생하지 않았던 극한 강우로 인해 돌발적으로 발생하는 경우가 빈번하게 발생하고 있어 이에 따라 사전에 예측하여 미리 대비하는 선제적인 홍수 대비 시스템이 요구되고 있다. 선제적인 홍수 대비 시스템의 구축 여부는 정확한 하천 흐름 예측을 필요로 한다. 하지만 하천의 흐름은 댐붕괴, 제방붕괴, 하천 하상의 변동 등 다양한 상황에서 급격한 흐름의 변동이 발생하며, 이는 하천 흐름 예측에 장애물로 작용하여 정확도를 떨어뜨리는 요인이 된다. 특히 국내에는 산악지형과 수공구조물에 의한 불연속 단면이 다수 존재하고 있어 그 예측 결과에 대한 정확성에 대한 요구가 더욱 부각되고 있다. 그렇기 때문에 해당 문제를 해결하기 위한 다양한 기법들이 개발되어 실무에 적용되고 있으나 어떤 기법이 국내 하천특성에 적합한지 파악할 수 없으며, 그 정확성과 안정성에 측면에서 여전히 많은 문제점을 가지고 있는 실정이다. 본 연구에서는 불연속 흐름이 빈번하게 발생하는 국내 하천 특성에 적합한 수치 기법을 확인하고자 유량보존특성을 만족하는 유한체적기법 중 국내외적으로 다수 사용되었던 기법들을 비교 및 평가하였다. 불연속 흐름의 대표적인 예제로서 'Dam-break problem'과 충격파 해석 및 홍수기와 갈수기에 따른 하천 하상의 마름/젖음에 대한 평가를 할 수 있는 Toro's Riemann problems에 적용하여 비교하였으며 그 결과 값을 정량적인 수치로 나타내었다. 이를 통해 국내 하천 특성에 적합한 수치 기법을 선정하였으며. 향후 국내하천 환경을 만족할 뿐만 아니라 하천 종사자들의 요구에도 부합한 하천흐름해석 모형의 개발 시 많은 기여가 될 것이라 판단된다.

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

댐 및 제방 등의 수공구조물 붕괴에 의한 극한홍수 전파를 해석하기 위한 2차원 홍수 해석기법에는 현재까지 다양한 기법들이 개발되어 왔고 다양한 모형들이 상용화 또는 범용화 되고 있다. 그 중 흐름의 전파양상을 정확하게 반영할 수 있는 상류이송기법인 Godunov형태의 유한체적기법은 충격파와 같은 불연속적인 해를 가지는 문제를 정확히 해석할 수 있고, 비구조적 격자 사용의 용이성 등의 장점 때문에 2차원 홍수파 전파 해석에 있어서 최근 십수년간 가장 활발하게 연구되어왔다. 하지만 이러한 기법은 양해법을 근간으로 하는 해석 기법으로써, 계산거리의 간격이나 계산시간의 간격, 격자망의 구성 등 엄격한 제한이 필요하다. 특히 방대한 계산시간을 요구하는 기법의 약점은 홍수예 경보 등을 위한 실시간 모형의 구동에 있어서 큰 제약이 되어 왔다. CPU의 성능이 지속적으로 발전하면서 이러한 문제는 점차 극복되어 왔으나, 발열 등의 문제와 이를 극복하기 위한 멀티코어의 등장으로 인해 단일 코어의 성능개발은 매우 더딘 것이 사실이다. 현재까지 연구되고 개발되어 온 모형들은 특별한 처리 없이는 단일 코어만을 사용하여 계산할 수 밖에 없기 때문에 멀티코어의 장점을 전혀 이용할 수 없다. 이러한 점을 극복하기 위해 프로그램을 병렬화하여 단일 문제에 대해 멀티코어를 사용할 수 있다면 계산시간 단축에 큰 효과를 거둘 수 있을 것이다. 현재까지 IT분야에서 다양한 병렬프로그래밍 기법들이 개발되고 소개되어 왔다. 본 연구에서는 병렬프로그래밍 기법중 가장 널리 사용되고 있는 MPI(Message Passing Interface)기법을 적용함으로써 기 개발된 고정확도 유한체적모형을 병렬화 하여 계산시간을 단축하고자 하였다. 개발된 모형을 장애물이 존재하는 실험하도의 댐 붕괴 문제와 2002년 태풍 루사 시 큰 피해를 입은 강릉시 일원의 섬석천 유역에 위치한 장현저수지와 동막저수지의 붕괴사례에 대해 적용하였다. 모형을 코어 개수별로 다양하게 모의함으로써 기존모형과의 결과에 대한 일치성을 확인하였고, 기존 모형 대비 계산시간 단축의 효과를 입증할 수 있었다. 개발된 본 모형을 실시간 홍수범람해석을 위한 시스템으로 구축할 수 있다면, 실시간 홍수예 경보에 있어 주요지점에서의 수위해석 뿐만이 아닌 제내지 범람 예보 분야까지 확대 적용할 수 있을 것으로 사료된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.198-201
• /
• 2010

지형이 불규칙한 자연하천에 대해 2차원 격자를 구성할 경우, 사각형 격자만을 사용한다면 지류와 본류의 합류부분에서 격자의 처리가 어려운 문제가 발생할 수 있으며, 삼각형 격자만을 사용하여 지형을 처리한다면 격자수가 많아져 계산시간이 다소 많이 소요되는 어려움이 존재할 수 있다. 혼합격자의 적용이 가능하다면 이러한 어려움은 어느정도 극복할 수 있다. 본 연구에서는 1차정확도 기법인 HLLC 기법을 적용하고, 지형이 복잡한 자연하천에 대한 격자처리의 유연성을 위해 삼각형 및 사각형 격자 그리고 이 두 격자가 혼용된 혼합격자의 적용이 가능한 2차원 유한체적모형을 개발하였다. 그리고 개발모형을 수리모형 실험을 통해 얻어진 실험자료가 존재하는 실험하도 및 실제 자연하천에서의 댐 붕괴에 대해 적용하여 결과를 비교하였다.

• Journal of the Korea Concrete Institute
• /
• v.19 no.3
• /
• pp.341-350
• /
• 2007

The volume control method which utilize a pressure node added into a finite shell element can overcome the drawbacks of conventional load control method and displacement control method. In this study, an improved volume control method is introduced for effective analysis of path-dependent behaviors of RC columns subjected to lateral cyclic loading or reversed cyclic loading along with compressive loading. RC shell structures and RC hollow columns are analyzed by discretizing the structures with layered shell elements and by applying in-plane two dimensional constitutive equations for concrete layers and reinforcement layers of the shell elements. The so-called path dependent volume control method as a finite element analysis technique is verified by comparing analysis results with other data including experimental results. The validity and applicability of the modeling technique is also confirmed by the comparison.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.12 no.3
• /
• pp.128-142
• /
• 2009

A lot of efforts have been made to resolve balancing problems between flux and source term and most finite volume models under development have been applied to simple channels such as imaginary and experimental channels. However, a number of numerical problems which can not be found in simple channels occur in the application to natural rivers and the problems should be resolved to apply finite volume models to natural rivers. In this study, 2D finite volume model which is applicable to natural rivers was developed and the accuracy of the developed model was validated through the application of partial dam break In addition, a simple and efficient 2D mesh generation method was suggested and the method can be accurately reflected to 2D mesh converted from surveyed cross sections in Han-river using GIS. The accuracy and applicability of the developed model on natural rivers were verified by performing simulation on Han-river using the generated mesh and comparing computed water elevation with measured water elevation.