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

본 연구에서는 ${\sigma}$-좌표계를 기반으로 하는 3차원의 이송확산 모형과 depth-integrated eddy simulation 모형을 결합한 효율적인 3차원 근역 (near-field) 해석모형을 제시하였다. 흐름 모형은 Boussinesq-type equations과 stochastic backscatter model을 기본으로 하고 있다. 이 흐름 모형은 수면의 변화와 바닥으로부터 발생하는 전단력과 파랑의 유동으로부터 발행하는 수심방향의 유속분포를 예측할 수 있다. 이와 같은 흐름 정보를 3차원 ${\sigma}$-좌표계의 이송확산모형에 제공하고 scalar의 이송과 확산에 대한 거동을 계산한다. 기본적인 이송과 이송-확산에 대한 검증 및 개수로에서 정량적 검증과 정성적 검증을 수행하였다. 전반적으로 타당한 결과가 도출되어 모형의 적합성이 있음을 확인하였다.

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

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

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

기후변화 및 해안 구조물의 증가 등 여러 원인이 연안침식 및 해안선 변화와 같은 연안의 지형변화를 가속하고 있다. 빠르게 변화하는 연안의 지형변화예측 및 대응책 강구를 위해서는 연안의 유사이송 현상에 대한 신속한 예측이 필요하다. 본 연구에서는 GPU 엔진 기반 파랑해석모형인 Celeris Advent를 활용하여 실시간으로 연안의 유사이송 모의가 가능한 수치모형을 개발하였다. Celeris Advent는 GPU의 병렬코어를 활용해 실시간 연산과 GUI를 통한 사용자와의 실시간 상호작용이 가능한 모형이다. 지배방정식은 확장형 Boussinesq 방정식에 유사이송방정식을 양방향 결합하여 구성하였고, 지배방정식에는 하이브리드 유한체적-유한차분 수치기법을 적용하여 이송항은 유한체적법(Kurganov & Petrova, 2007), 소스항은 유한차분법을 통해 이산화하여 해석한다. 유사이송방정식은 수심적분형 이송확산방정식에 침식 및 퇴적 플럭스를 반영하는 소스항을 결합하여, 이송항 및 확산항을 통해 유사의 이송/확산을 고려함과 동시에 소스항을 통해 하상과의 상호작용을 고려하였다.

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

본 연구에서는 GPU 가속화 기반의 Boussinesq 모형인 Celeris Advecnt에 수심 적분된 2차원 이송-확산방정식을 추가하여 인터렉티브 시스템 기반의 추적자 이동 모형을 개발하였다. Celeris Advent는 최초로 개발된 인터렉티브 시스템을 갖춘 Boussinesq 모형으로, 시뮬레이션 중에 사용자가 모형의 파라미터뿐 아니라 모델 도메인 내 수위 및 수심을 바꿀 수 있다. 이를 통해 사용자는 모의가 진행되는 도중에 모델의 안정성 및 효율성을 위해 시간 간격을 조정할 수 있을 뿐 아니라 방파제 설치 등과 같은 지형 변화를 고려하기 위해 도메인 내 격자별 수심을 조정할 수 있다. 본 연구에서는 연안에서의 추적자 이동 모의를 위해 Boussinesq 방정식과 더불어 이송-확산방정식을 풀이하는 추적자 이동 모형을 개발하였다. 추적자의 확산항의 경우 분자 자체의 확산과 더불어 쇄파에 따른 난류 확산을 고려하였다. 난류 확산계수는 슈미트 수를 1로 두어 와동점성계수와 동일하게 두었으며, 와동점성은 단순화된 형태의 쇄파모형을 고려하여 계산하였다. 쇄파모형의 고려로 인해 이송-확산방정식과 더불어 운동량 방정식에서도 쇄파에 따른 운동량 소산이 고려되었다. 마지막으로, 추적자 농도에 대한 인터렉티브 시스템을 추가하여, 모델 구동 중에도 사용자가 수심적분된 추적자 농도를 조정할 수 있도록 하였다. 기수행된 2개의 수리실험 조건과 관측값을 이용하여 벤치마크 테스트를 수행하였으며, 관측값과 대체로 일치하는 것을 확인하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 1987년도 제29회 수공학연구발표회논문초록집
• /
• pp.201-214
• /
• 1987

하천에서 종, 횡방향이 고려된 확산-이송 현상을 수치기법인 유한요소법을 이용하여 2차원으로 해석하였다. 유한요소법으로는 Galerkin의 가중잔차 방법을 수십에 대해 적분을 취한 연속, 운동량 및 확산-이송방정식에 적용하였고, 선형보간함수와 선형삼각형요소가 이용되었다. 모형의 타당성을 입증하기 위해 단순화된 1차원 수로에서 수차례 검정한 결과 정확해와 거의 일치하는 만족할만한 결과가 도출되었다. 개발된 모형의 실험이 2차원수로에서 행하여져 지류의 유입에 따른 확산-이송현상이 모의되었으며, 실험적용은 개발사업후의 한강본류 9km 구간에 적용되어 탄천과 중량천의 지천 영향을 받는 오염 농도가 2차원적으로 도시되었다.

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

오염원의 이송확산에 관한 많은 연구들이 수행되어 왔으나 특히 비포화 영역에서 오염원 이송확산을 측정하는 것은 매우 어려운 것으로 알려지고 있다. 비포화 토양에서의 오염원 이송확산은 매질의 함수량 변화에 영향을 받기 때문에 오염원 거동특성을 이해하려면 비포화 흐름 분석을 선행한 후 오염원의 이송확산 특성을 분석하여야 한다. 본 연구에서는 비포화 영역에서의 오염원 이송특성을 분석하기 위하여 TDR(Time Domain Reflectometry)을 이용하여 비포화 흐름 및 오염원 이송을 측정하였다. 이를 위하여 본 연구에서는 TDR을 이용하여 오염원 이송을 측정하는 방법을 개발하였으며, 이 방법을 이용하여 1차원의 토양기둥시료에서 비포화 흐름 및 오염원 이송확산에 관한 실험을 수행하고 수치모형을 적용함으로써 비포화 영역에서 오염원의 이송확산에 관한 거동특성을 규명하였다. 본 연구에서는 두 종류의 국내 토양시료(SUS, KUS)를 사용하였는데, 토양의 물리적 특성을 예비실험을 통하여 규명한 후 토양기둥시료를 이용한 본실험을 수행하였다. 비포화 천이흐름하의 오염원 이송확산 실험에서는 급격한 습윤전선의 전진에 따른 종형의 함수량변화를 관측할 수 있었고, 이때 오염원의 농도는 함수량의 천이구간의 중심점으로부터 전방영역의 농도분포가 습윤전선에서의 함수량 분포와 유사한 종형을 이루고 있음을 관측할 수 있었다. 비포화 정상흐름하의 오염원 이송확산 실험에서는 오염원이 이송하며 농도 천이구간이 확장되어지는 전형적인 형태를 보였다. 또한 예비실험에서 측정한 매개변수를 입력자료로하여 수행한 수치결과와 실험결과를 비교하였는데 비포화 흐름특성은 실험결과와 수치결과가 정량적으로 일치하는 경향을 보였으나, 오염원 이송확산 특성은 정량적으로 수치결과가 실험결과보다 더 많이 확산되는 경향을 보였다. 따라서 수치모형을 현장에 적용할 경우 확산지수 결정에 주의하여야 할 것으로 판단된다. 즉, 수치모형에 적용할 확산지수는 BTC 실험을 통하여 측정한 확산지수, 수치확산, 흡착계수, 적용영역의 크기 등을 고려하여 결정하여야 한다. 특히 본 논문에서는 TDR을 이용하여 최초로 천이상태의 함수량과 오염원 농도를 측정하였는데 이를 위하여 전기전도도와 함수량관계를 추정하는 식을 제안하였으며, 전기전도도와 토양수 농도, 전기전도도와 함수량의 관계를 이용한 천이상태의 오염원 농도 측정방법을 개발하였다. 특히 제안식에서는 한계함수량의 개념을 도입하여 전기전도도와 함수량관계를 추정하므로 추정식의 실험값 반영 정도를 증가시켰다. 본 연구에서 제안된 식을 이용하여 추정된 전기전도도와 함수량관계는 다른 제안식에 비하여 개선된 결과를 보여 주었고, 본 연구에서 개발한 오염원 농도 측정법을 이용하여 측정한 결과 함수량이 0.15이하에서는 측정오차가 크지만 함수량이 0.15이상일 경우 매우 좋은 결과를 보였는데 질량평형을 검토한 결과 약 5-10%의 오차율을 보였다. 따라서 본 논문에서 개발된 천이상태의 오염원 농도측정법은 용존 오염물질의 이송에 관한 정확한 실험을 제공할 것으로 판단된다.

• 한국수자원학회논문집
• /
• 제34권5호
• /
• pp.499-509
• /
• 2001

자연하천에서의 이송-확산 과정의 모의를 위하여 입자위치의 이산확률분포에 기초한 2차원 수송 모형을 개발하였다. 제안된 모형에서는 단위 시간간격동안 격자간의 질량이송을 예측하기 위하여 평균과 분산의 함수로 나타내어진 확률분포를 사용하였다. 개발된 모형은 유속, 확산계수, 단면적이 일정한 단순영역에 대하여 수치확산이 없는 해를 구하였고, 양의 확률을 만족시키는 안정조건이 성립한다면, 해석해와 다른 유한차분법과 비교하였을 때, 좋은 결과를 나타내었다. 본 모형의 현장적용성을 검토하기 위하여 캐나다에 위치한 Grand River를 대상을 얻은 수치실험 결과를 정상상태의 색소실험 결과와 비교하였다. 그 결과로서 본 모형은 자연하천에서의 2차원 이송-확산을 잘 모의할 수 있는 것으로 나타났다.

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

하천, 호소 등에서 공사 시 준설작업으로 인해 고탁수가 발생되며 이러한 탁수의 이송 확산은 수환경에 변화를 초래한다. 탁수는 준설작업 시 커터헤드(cutter head)가 해저면에 닿는 순간부터 작동을 멈출 때까지 계속하여 발생하며, 이러한 과정들이 반복되면서 많은 양의 부유토사가 발생하게 되고 고탁수현상이 일어난다. 이렇게 발생한 탁수는 수체흐름에 따라 이송 및 확산된다. 탁수발생은 수중의 빛 투과를 감소시켜 일차 생산자인 부유성 및 부착조류의 생육을 저해하고, 이들의 생산성 감소와 군집구성의 변화는 수서생태계의 먹이사슬을 통해 이들을 먹이로 하는 저서무척추동물과 어류의 현존량 감소와 종 구성에 영향을 미치고, 고농도의 현탁 입자는 어류 아가미에 염증을 유발하거나 점막의 파괴와 감염을 유발하여 치사시킬 수도 있다. 또한 과도한 부유 입자는 하류로 침강되어 하천 바닥에 서식하는 부착조류, 무척추 동물 및 곤충의 생육에 피해를 주고, 이것은 어류의 먹이에 영향을 미쳐 어류 개체수를 감소시키거나 산란된 물고기 알을 매몰시키거나 질식시키는 등 여러 가지 방법으로 수서생물상에 영향을 미치게 된다.(낙동강수계관리위원회, 2005) 따라서 준설작업에 있어 탁수의 이송 확산범위를 사전에 예측하고 국내 실정과 환경여건에 알맞게 적용되고, 실용화될 수 있는 수치모델링에 대한 기반핵심 기술개발이 필요하다. 현재 낙동강에서 진행되고 있는 준설현장에서 발생하는 부유탁수의 이송 확산과정을 이차원 흐름해석모형인 RAM2 모형과 오염물 이송 확산해석모형인 RAM4 모형을 이용하여 수치해석을 하고 분석함으로써 수치해석에 의한 부유탁수의 이송 확산모의 결과가 환경영향 범위를 예측하는 데에 적용될 수 있는지를 알아보고자 한다.

• 물과 미래
• /
• 제27권2호
• /
• pp.155-166
• /
• 1994

Eulerian-Lagrangian 방법을 이용하여 1차원 종확산방정식의 수치모형을 비교·분석하였다. 본 연구에서서 비교·분석한 모형은 지배방정식을 연산자 분리방법에 의해서 이송만을 지배하는 이송방정식과 확산만을 지배하는 확산방정식으로 분리한다. 이송방정식은 특성곡선을 따라서 유체입자를 추적하는 특성곡선법을 사용하여 해를 구하고, 그 결과를 고정된 Eulerian 격자상에 보간하였고, 확산방정식은 상기 고정격자상에서 Crank-Nicholson 유한차분법을 사용하여 해를 구하였다. 이송방정식의 풀이에서 다양한 보간방법이 적용되었는데, 일반적으로 Hermite 보간다항식을 사용한 경우가 Lagrange 보간다항식을 사용한 경우보다 수치확산 및 수치진동 등의 오차를 최소화할 수 있어서 더욱 우수한 것으로 밝혀졌다.

• 물과 미래
• /
• 제29권4호
• /
• pp.187-198
• /
• 1996

난류전단 흐름에서의 비정상 수평 선오염원의 확산에 관한 수치모형을 개발하였다. 계산기법으로는 비정상 이송확산 방정식을 종방향 이송 및 연직방향 확산으로 분리하고, 이들 방정식을 방시간 간격에 대하여 번갈아 계산하는 단계분리 유한차분기법을 사용하였다. 종방향 이송방정식에 대해서는 Holly-Preissmann 기법을, 연직방향 확산방정식에 대해서는 Crank-Nicholson 기법을 각각 적용하였다. 개발된 모형을 난류전단 흐름에서의 정상 수평 선오염원의 확산 문제에 적용하여 계산결과를 반해석해와 비교함으로써 모형을 검증하였다. 또한 난류전단 흐름내로 순간적으로 방류된 면오염원의 확산문제에 계산모형을 적용하였다. 마찰계수에 대한 민감도 분석 결과, 동일한 무차원 시각에서의 혼합 정도는 마찰계수에 관계없이 거의 일정한 것으로 나타났다. 또한 동일한 정도의 혼합상태에 도달하는 데 소요되는 유하거리는 마찰계수의 제곱근에 반비례함을 알 수 있었다.