• 한국수자원학회논문집
• /
• 제35권1호
• /
• pp.77-90
• /
• 2002

하상보호용 사석은 초기 거동조건의 범위를 넘어 한계조건에 도달하면 흐름에 의해서 각각의 자갈 또는 석재가 분리되면서 점진적으로 파괴된다. 이와 반대로 하상보호용 매트리스의 경우는 비록 초기 거동조건의 범위를 넘어 한계조건에 도달할지라도 매트리스 철망에 의해 각각의 자갈 또는 석재의 분리가 억제가 됨을 알 수 있었다. 비록 변형된 하상보호용 매트리스일지라도 새로운 형태의 평형을 얻을 수 있으며, 더 이상의 변형과 나쁜 조건에서도 저항할 수 있음을 확인하였다. 동일한 수리학적 거동조건에서 하상보효용 매트리스의 채움재로 사용된 자갈 또는 석재의 평균 크기가 하상보호용 사석의 1/2정도이며, 하상보호용 매트릭스에서 Shield 계수는 하상보호용 사석의 약 2배임을 알 수 있었다. 또한 동일한 크기의 자갈 또는 석재가 채움재로 사용될 때 허용유속에 대해서 하상보호용 매트리스의 경우 하상보호용 사석보다 2배이상 심지어 3∼4까지 견딜 수 있음을 알 수 있었다.

• 터널과지하공간
• /
• 제11권3호
• /
• pp.257-269
• /
• 2001

지하 500 m의 화강암반내 단층지역에 위치한 지하 방사성폐기물 처분장의 구조거동을 이해하기 위하여 수치 해석을 수행하였다. 해석에는 2차원 해석코드인 UDEC을 사용하였다. 해석모델은 화강암반, 처분공내의 압축 벤토나이트로 둘러싸인 PWR 사용후 핵연료 처분용기 및 처분동굴내에 채워진 혼합 벤토나이트를 포함한다. 한 개의 단층이 처분동굴의 지붕과 벽이 만나는 지점을 33, 45, 및 $58^{\circ}$의 각도로 관통하는 세가지 다른 경우게 대한 구조거동을 비교, 분석하였다. 그리고 $45^{\circ}$단층의 경우에 대해서는 수리역학적, 열역학적, 및 열수리역학적 상호작용 거동을 해석하고 비교, 분석하였다. PWR 사용후 핵연료내의 방사성 물질로부터 나오는 시간의존 방사성 붕괴열에 의한 영향을 해석하였다. 지하수위는 지표면 아래 10 m로 가정하였고, 지하수해석은 정류 알고리즘을 사용하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2008년도 학술발표회 논문집
• /
• pp.1874-1878
• /
• 2008

본 연구에서는 물고기 서식처의 수리조건을 조사하기 위한 현장조사와 대상지에서의 미시서식처의 분포를 파악하기위한 수치계산을 실시하였다. 현장조사는 미호천과 금강이 합류되는 지점으로부터 약 24km 상류에 위치한 동경127도26분-북위36도41분 충북 청원군 오창면의 미호천 일부구간을 대상으로 수리학적 특성 조사, 지형측량 및 서식처별 대표어종 조사를 수행하였으며, 수치실험은 하천에서의 흐름거동을 2차원적으로 해석하는데 많이 사용되고 있는 상용프로그램인 SMS (Surface - Water Modeling System)의 RMA2 모형에 현장조사로부터 얻은 결과들을 입력함으로서 진행하였다. 그 결과 유속과 수심에 따른 서식처 구분이 가능하였으며, 현지특성을 재현한 수치모형 결과는 대상하천에서의 서식지 분포를 정량적으로 나타내었다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2010년도 학술발표회
• /
• pp.724-728
• /
• 2010

강우강도의 증가와 식생대의 변화 등 기후변화로 인한 저수지내 대형탁수와 조류대발생으로 사회적 문제가 야기되고 있다. 이에 유역대책 평가와 저수지내 유체거동 및 영향평가를 위한 수단 확보가 절실한 상황이다. 본 연구에서는 국내 저수지에 대한 3차원 수리해석기술의 적용성을 검토하고자 하였다. 이를 위해 낙동강의 주요지류 중 하나인 남강에 위치한 남강댐에 대해 미 EPA의 오염총량지원도구인 EFDC 모델을 적용하였다. 저수지내 퇴사량 조사시 측정된 실측치와의 지형 재현성 검토, 실측수위와 모의결과의 비교를 그래프와 통계지표를 활용하여 실시하였고 염료 추적 모의를 통해 도달시간 등을 산정하였다.

• 터널과지하공간
• /
• 제28권5호
• /
• pp.400-425
• /
• 2018

본 논문에서는 국제공동연구인 DECOVALEX-2019 프로젝트 Task B의 연구결과와 현황을 소개하였다. Task B의 주제는 'Fault slip modelling'으로 유체의 주입으로 인해 발생하는 단층의 재활성(미끄러짐, 전단파괴)과 수리역학적 거동을 예측할 수 있는 해석기법을 개발하는 데에 그 목적이 있다. 1단계 연구는 참가팀들이 연구주제에 대해 숙지하고, 벤치마크 모델을 대상으로 단층의 투수특성과 역학적 거동의 상호작용을 모사할 수 있는 해석코드를 개발할 수 있도록 하는 준비 단계의 연구이다. 본 연구에서는 TOUGH-FLAC 연동해석 기법을 사용하여 물 주입으로 인한 단층의 수리역학적 연계거동을 모사하였다. TOUGH2 해석에서는 단층을 Darcy의 법칙과 삼승법칙을 따르는 연속체 요소로 모델링하였으며, FLAC3D 해석에서는 미끄러짐과 개폐가 허용되는 불연속 인터페이스 요소를 통해 모사하였다. 두 가지 수리간극모델에 대하여 수리역학적 커플링 관계식을 수치화하였으며, 연속체 요소(수리모델)와 인터페이스 요소(역학모델)의 거동을 연계할 수 있는 해석기법을 제시하였다. 또한, 단층의 역학적 변형(간극의 변화)으로 인한 수리물성 변화와 기하학적 변화(해석 메쉬의 변형)를 수리해석에 반영할 수 있는 해석기법을 개발하였다. 다양한 압력의 물을 단계적으로 주입하고 이로 인해 유도되는 단층의 탄성거동 및 전단파괴(미끄러짐)에 대해 살펴보았으며, 수리간극의 변화 양상과 원인, 압력 분포와 주입율의 관계 등을 면밀히 검토하였다. 해석 결과, 본 연구에서 개발한 해석기법이 물 주입으로 인한 단층의 미끄러짐 거동을 합리적인 수준에서 재현할 수 있는 것으로 판단할 수 있었다. 본 연구의 해석모델은 Task B에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 연구의 현장시험에 적용하여 타당성을 검증할 예정이다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2022년도 학술발표회
• /
• pp.237-237
• /
• 2022

우리나라 하천 내에 설치된 횡단구조물은 보 35,000개소 이상으로 전국 하천의 0.6 km마다 수공구조물이 위치해 있으며 이로 인해 종적(longitudinal) 연결성에 영향을 미치는 것으로 나타났다. 농어촌연구원의 조사에 따르면 하천 횡단구조물이 설치된 구간에 어도 설치율은 14.9%로 나타났으며 이는 수생태 관점에서 연속성이 매우 열악한 상태임을 확인 할 수 있다. 이에 우리나라 각 부처 국토교통부, 환경부 등에서는 수생태 연속성을 확보하고자하는 노력이 지속되고 있으며 대표적으로 어류의 이동통로인 어도를 설치하거나 기존 어도의 효율을 향상시키기 위하여 개보수 작업을 지속적으로 실시하고 있다. 어도의 기능을 평가하기 위해서는 어도 내의 수리특성을 정확히 파악하는 것이 중요한데 연속적인 구조물로 구성된 어도 내 흐름은 매우 복잡하다. 특히 어도 내 구조물간 상호작용에 의하여 비정상 흐름이 발생하며 구조물 뒤에서는 사수역(dead zone)이 형성된다. 사수역에서 나타나는 와류의 거동은 구조물의 기하학적 특성에 따라 변화한다. 본 연구에서는 2차원 수심적분 수치모형을 활용하였으며 벽면 근처 점성저층의 유속분포을 재현하기 위하여와 점성항에 감쇠함수(damping function)를 고려하였다. 수치모형의 검증을 위해 실내 실험수로의 직선 개수로에서 PIV(particle image velocimetry)를 활용하여 연속적으로 배치된 구조물에 의한 유속자료를 활용하였다. 이 결과는 향후 새로운 어도설계 혹은 기 설치된 어도의 수리학적 기능을 평가하는데 활용이 가능할 것으로 보인다.

• 터널과지하공간
• /
• 제23권5호
• /
• pp.372-382
• /
• 2013

고준위폐기물처분장에 대한 열-수리-역학적 거동 모델링은 처분장의 성능 및 안전성 평가를 위해서 선행되어야 할 중요한 연구과제이다. 본 연구에서는 $\ddot{A}$sp$\ddot{o}$ 원형처분장의 열적 거동을 해석하고, 현장 실험데이터와의 비교를 통해 모델 계산결과의 타당성을 검증하였다. 모델 시뮬레이션에서는 처분공과 처분터널 및 그 주변 암반에 대한 온도분포를 분석하였다. 현장 실험데이터와의 비교는 처분공 DH-6를 대상으로 수행하였다. 그 결과 모델 계산치는 측정위치에 따라 실험치 보다 약 2-$5^{\circ}C$ 정도 높은 온도값을 보였으나, 온도변이 곡선은 비슷한 패턴을 보여 주었다. 실험치와 모델 계산치의 이러한 차이는 모델에 의한 열적거동 해석에서 암반을 제외한 완충재, 벤토나이트 펠렛, 뒷채움재 내의 수리학적 및 역학적 거동을 고려하지 않았기 때문으로 판단되었다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2015년도 학술발표회
• /
• 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값에 미치는 영향을 반영하여 수리학적으로 매끄러운 벽면에서 적용이 가능한 수정된 대수중복층 식을 제시하고자 한다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.423-427
• /
• 2007

현재까지 하천의 흐름, 유사이동, 수질해석을 위해서는 외국에서 개발된 소프트웨어를 주로 사용해 오고 있었다. 학술 분야에서는 국내의 모형들이 연구되어졌지만 그에 따른 GUI나 가시화 시스템에 대한 실용화는 거의 이루어지지 못하였다. 본 연구에서는 범용 2차원 하천 흐름, 수질, 유사이동 해석을 위한 GUI 및 가시화 시스템(이하 RAMS, River Analysis and Modeling System)을 개발하여, 하상변동 및 오염물 이송확산에 미치는 수리학적 영향을 규명할 수 있도록 하였다. RAMS는 크게 mesh generator, 해석 모형의 입력 GUI 모듈, 입출력 파일 생성 모듈, 그리고 모의 결과의 가시화 시스템 등으로 이루어져 있다. Mesh generator는 지형자료(이미지 또는 DXF 파일)를 백그라운드 이미지로 가져올 수 있으며, 삼각형 노드와 사각형 노드를 지원한다. 또한 thin triangle들을 제거하는 기능, 선택된 요소(elements)를 제거하는 기능, triangle들을 서로 병합하여 사각형 요소를 만드는 기능, mesh의 renumbering 기능 등을 구현하였다. 특히 사용자가 잘못 생성한 요소들을 바로 이전 상태로 환원하는 undo/redo 기능을 구현하여 능률적인 mesh 생성이 가능하다. 해석 모형의 입력 GUI 모듈에는 각 해석 모형(흐름, 수질, 유사이동)에 특화된 GUI를 설계하여 사용자는 보다 친숙한 환경에서 편리하게 자료를 입력할 수 있다. 입출력 파일 생성 모듈에서는 사용자가 GUI를 통해 입력한 자료를 파일로 변환하여 즉각적으로 모의를 수행하며, 그 출력 파일을 읽어 모의 결과를 자동적으로 가시화한다. 모의 결과의 가시화 시스템에서는 수많은 모의 결과를 체계화하여 등고선 및 화살표 등으로 표현하며, time step 별 결과를 바로 확인할 수 있다. 또한 애니메이션 기능을 구현하여 사용자가 원하는 구간의 time step에서의 모의 결과를 연속적으로 볼 수 있으며, 이 애니메이션을 AVI 파일로 변환하여 다른 동영상 프로그램에서도 재생할 수 있다. 본 연구에서 개발한 RAMS를 이용하여 하천 설계 시 그 공학적 근거를 제시하고, 국내 실정에 맞는 국산 소프트웨어를 제공함으로써 하천의 흐름, 수질, 유사이동 해석에 의한 하천의 수리학적 거동을 보다 편리하고 정확하게 모의할 수 있을 것으로 기대된다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2023년도 학술발표회
• /
• pp.127-127
• /
• 2023

지표수와 하상 경계층에서 발생하는 흐름 교환은 하천, 호소, 연안, 해안 등 자연계에 존재하는 수환경시스템에서 일반적으로 나타나는 수리학적 특성으로서, 흐름 교환이 발생하는 경계층 아래 하상층 영역을 혼합대(hyporheic zone)라 부른다. 수질오염사고 등에 의해 외부의 오염물질이 하천 내 유입될 경우, 혼합대 흐름에 의해 하상층으로 침투되고 지표수 대비 유속이 느린 하상 내 공극 흐름에 의해 거동함에 따라 이들의 하천 내 체류시간이 증가하게 된다. 따라서, 본 연구에서는 지표수와 하상 흐름을 연계한 수치해석 방법을 적용하여 혼합대가 지표수 내 용질 체류시간에 미치는 영향을 분석하였다. 먼저 연직 2차원 Reynolds 평균 Navier-Stokes(RANS) 방정식과 Darcy 방정식을 연계하여 지표수와 하상 내 흐름을 해석하였다. 지표수 영역은 RANS 방정식을 이용하여 모의하였고, 지표수 흐름해석에서 얻어진 하상의 압력장을 경계조건으로 하여 Darcy 방정식과 함께 하상 내 흐름을 모의하였다. 여기서 하상의 형태는 자연계 하천에서 일반적으로 관찰되는 사련하상(Ripple bed)으로 모사하였다. 이후, 지표수-하상 연계모의를 통해 얻어진 흐름 결과를 바탕으로 지표수-하상 경계층에서 용질거동을 모의하였다. 흐름 모의결과를 과거 실험자료와 비교한 결과, 지표수 영역 내연직흐름 분포를 정확하게 재현하였고, 동시에 혼합대 흐름 구조에 큰 영향을 미치는 지표수-하상 경계층 압력 분포 역시 관측값과 유사하게 나타났다. 용질거동 해석을 통해 얻어진 용질의 체류시간을 분석한 결과, 혼합대 흐름이 고려된 경우(투수성 하상)와 고려되지 않은 경우(불투수성 하상)를 비교했을 때 전자에서 체류시간 분포의 감수곡선이 길어지고 첨두농도가 감소하는 것으로 나타났다. 아울러, 지표수 영역의 유입부 경계의 평균 유속이 증가함에 따라 최대 체류시간이 감소하는 것으로 나타났는데, 이는 지표수-하상 경계층에서의 압력 경사가 커져 혼합대 내 유속이 증가함에 기인하는 것으로 분석되었다.