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

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

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2011년도 정기 학술발표대회
• /
• pp.44-44
• /
• 2011

저수지와 같은 갇혀진 수체는 상류에서 유입되는 오염물질 뿐만 아니라 성층현상에 의해서도 오염될 수 있다. 안정된 성층은 혐기성 조건을 제공하여 바닥에서의 과도한 조류성장을 유발하고 연직순환흐름을 방해한다. 갇혀진 수체에서의 연직순환은 환경적인 문제를 감소시킬 수 있도록 중요한 역할을 한다. 갇혀진 수체에서의 연직순환은 이러한 오염을 줄이는데 중요한 역할을 하는데, 연직순환을 일으키는 인자로는 빛의 입사, 바람, 물의 온도 및 열의 확산 등이 있으며, 그중에서도 가장 중요한 것은 바람의 영향이다. 그러므로 성층화된 흐름에서 바람에 의해 발생하는 연직순환에 대한 수치모형을 개발하고 적용하는 것이 필요하다. 본 연구에서는 수온성층흐름에서의 순환흐름을 해석하기 위하여 유속성분을 계산하고 자유수면 변위와 온도, 염도등과 같은 스칼라양의 해석을 위해 다음과 같은 3단계의 방법을 이용하였다. 첫 번째 단계(정수압 계산단계)에서는 운동량 방정식의 경사항을 음해적으로 해석하고, 두 번째 단계(자유수면 보정단계)에서는 자유수면의 변화를 계산하고 수평방향 유속성분을 계산한다. 예측-수정자 방법(predictor-corrector step method)을 이용하여 자유수면변위와 유속을 구하였다. 마지막으로 세 번째 단계(이송-확산단계)에서는 이송-확산 방정식을 이용하여 스칼라양을 계산하였다. 본 연구에서 제시한 모형의 정확도를 검증하기 위하여 정사각형수조에서 진동하는 자유수면의 해석해와 비교하였고, 성층화된 흐름에서 발생하는 연직순환에 대하여 수치모의를 실시하였다. 그 결과, 본 연구에서 개발된 수치모형이 흐름 내부의 현상을 잘 묘사함을 알 수 있었다.

• 한국해양환경ㆍ에너지학회지
• /
• 제10권1호
• /
• pp.53-58
• /
• 2007

하구에서 오염물은 salt finger가 발생할 수 있는 적당한 조건하에서 혼합과 수송이 일어날수 있다(Hwang and Rehmann, 2004). 선형이론을 salt finger가 일어날수 있는 적절한 조건하에 유체의 운동을 예측하는데 적용하였다. 모의 결과는 기존의 실험 결과와 거의 비슷한 결과를 도출하였다. 밀도율이 2보다 클 때, Turner(1967)가 열과 소금을 이용한 실험에서 발견한 것처럼 혼합율(the flux ratio)은 $0.55{\sim}0.57$를 보이며, 소금과 설탕을 이용한 Griffiths(1980)의 실험에서와 같이 0.87의 혼합율을 도출하였다. 두개의 매개 물질의 분자확산계수가 증가함으로써 이송속도율도 밀도율과 함께 증가하였고, 높은 밀도율에서 이송속도율이 정상상태가 되는데 걸리는 시간이 증가하였다.

• 한국해안해양공학회지
• /
• 제10권1호
• /
• pp.37-44
• /
• 1998

본 연구에서 제안하는 혼합 방법(hybrid method)은 흐름이 우세한 영역에서의 전송 문제를 정확하고 효과적으로 해결하기 위하여 개발된 것으로 오일러-라그란쥐적 방법과는 달리 전방추적에 의하여 이송 과정이 수행되므로 보간 기법이 불필요하고 무작위 행보에 의한 라그란쥐적 방법과 달리 유한 차분법에 의하여 확산 과정이 수행되므로 많은 입자가 요구되지도 않는다. 한 점에 순간적으로 부하되는 오염원과 연속적으로 부하되는 오염원에 대한 이론적인 해와 비교하여 확산 계수와 무관하게 상당히 만족할 만한 결과를 얻었다. 현 방법은 또한 2차원 상에서 주변 5격자로부터 보간하는 오일러-라그란쥐적 방법과 무작위 행보로 입자 추적하는 순수 라그란쥐적 방법과 비교하여 정확성은 물론 계산 시간에 있어서도 상당히 월등한 방법임이 입증되었다.

• 물과 미래
• /
• 제28권4호
• /
• pp.125-136
• /
• 1995

특정곡선을 고려한 ELM을 이송-확산방정식에 적용하여 그 결과를 Eulerian 기법(Stone-Brian, QUICKEST)과 비교하였다. 이송항의 계산을 위해서는 Lagrangian 보간법과 Cubic spline 보간법을 이용하였고 확산항의 계산에 있어서는 Crank-Nicholson 방법을 이용하였다. 수치모형의 적용결과는 다음과 같다. (1) Gaussian hill에의 적용:Lagrangian 보간법을 사용하여 계산한 경우가 가장 정확한 결과를 보였다. Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우에는 Peclet수가 50인 경우에 감쇠현상을 보였다. Stone-Brian방법은 Peclet수 10,50에서 위상오차가 발생하였다. (2) Advanced front에의 적용: 모든 방법이 Peclet수 1,4에서 정확한 결과를 얻었다. Peclet수가 50인 경우에 Lagrangian 보간법을 사용하여 계산한 경우와 Stone-Brian 방법은 증폭오차가 발생하였고 Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우는 수치진동 현상을 보였다.

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

수환경으로 유출되는 유해화학물질은 독성을 가지고 직접 유출되거나 다양한 매체와 반응하여 화재 및 폭발 등의 사고가 발생한다. 실제로 낙동강 유역에서는 1991년 페놀 유출사고를 시작으로 2009년 구미공단 '1,4-다이옥산' 유출사고, 2014년 11월 경북 봉화군의 황산유출사고 등 크고 작은 사고가 빈번히 발생하고 있으며 작년 6월에는 대구와 부산의 수돗물에서 과불화화합물이 검출되기도 하였다. 이러한 대규모 사고를 방지하기 위해 신속한 오염물의 거동 예측이 가능한 추적모델이 필요하며, 본 연구에서는 수환경으로 유출된 유해화학물질의 추적을 위한 1차원 저장대 모형을 개발하였다. 일반적으로 저장대 모형은 복잡한 하천 구조를 하천의 주 흐름이 존재하는 본류대와 하천 흐름이 정체되는 저장대, 그리고 하상구조로 단순화 하여 나타낸다. 본류대에서는 하천흐름에 의한 이송 및 횡방향 유속차로 발생하는 전단류에 의한 확산이 일어나며, 저장대와의 물질교환으로 발생하는 저장효과와, 하상구조와의 흡착 및 탈착, 그리고 생물화학적 반응 및 휘발이 발생한다고 가정한다. 본류대와 저장대간의 질량교환은 난류유속변동과 농도차에 의해서만 발생한다고 가정하고 오염물질의 이송과 분산과정을 해석한다. 저장대에서는 이송 및 전단류에 의한 확산은 일어나지 않으며, 본류대와의 물질교환으로 발생하는 저장효과와 하상구조로의 흡착, 그리고 생물화학적 반응 및 휘발이 발생한다고 가정하며, 하상구조에서는 본류대 및 저장대와의 흡착 및 탈착만 발생한다고 가정한다. 저장대 모형의 해석을 위해서는 리치(Reach) 별로 본류대 분산계수($K_F$), 본류대 면적($A_F$), 저장대 면적($A_S$), 그리고 저장대 교환계수(${\alpha}$)의 네 가지 저장대 매개변수가 필요하며 본 연구에서 개발된 저장대 모형은 흡탈착, 생물화학적 반응 및 휘발 과정을 모두 고려하여 유해화학물질의 확산 거동을 모의한다. 최적의 리치길이, 흡탈착, 반응 및 휘발 계수를 산정하여 모형의 정확도를 향상시켰으며, 신속하고 정확하게 오염물의 거동을 예측할 수 있었다.

• 물과 미래
• /
• 제23권1호
• /
• pp.119-127
• /
• 1990

점착성 부유사 농도분포는 해수유동과 물질 확산에 의해서 결정되며 지배방정식으로는 2차원 수심적분된 Reynolds운동방정식, 연속방정식과 Fick의 확산법칙에 근거를 둔 대류-확산방정식이 사용되었다. 해수유동과 점성퇴적물 확산인 두개의 모형은 유한차분법을 이용하였고 유동모형은 양해법, 확산모형은 다증법을 사용하여 부유사 이동의 현상을 파악하였다. 해수유동방정식의 적용시 이송항의 포함여부에 대해서 조사하였으며 물질확산 방정식에 대해서는 한계전단응력값의 변화가 부유사농도에 영향을 주는가에 대해서 비교하였다.

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

오염물질의 이송${\cdot}$확산 과정을 해석하기 위하여 계산효율이 높은 3차원 퍼프모형을 개발하였다. 본 연구에서 개발된 퍼프모형은 추적방법에 따라 전방추적모형과 후방추적모형으로 나눌 수 있으며, 이 두 가지 추적방법은 계산효율과 수치오차에 있어서 상이한 특성을 나타내었다. 전방추적 퍼프모형은 일정한 시간간격을 가지므로 정상상태의 연속오염원의 경우에 각각의 퍼프들이 동일한 질량을 갖는다. 그러므로 전방추적 퍼프 모형은 각 퍼프들간의 중첩정도가 일정하지 않다. 이에 관한 오차분석을 수행하기 위하여 본 연구에서는 무차원 퍼프중첩계수를 정의하였다. 퍼프중첩계수란 퍼프의 크기에 대하여 퍼프중심간의 거리가 떨어진 정도를 나타내는 무차원수로서 너무 작은 경우에는 정확도가 떨어지고 너무 큰 경우에는 계산효율이 감소한다. 전방 추적 퍼프모형의 경우, 중첩계수가 작은 초기구간에는 정확도가 떨어지며, 시간이 지남에 따라 중첩계수가 필요이상으로 증가하여 계산효율이 떨어지는 것으로 나타났다. 이에 비해 일정한 중첩정도를 갖는 후방추적 퍼프모형의 경우에는 전 영역에 걸쳐서 정확도와 계산효율이 높은 것으로 나타났다. 하지만 일정한 시간간격 마다 농도장을 계산하고자 할 때, 전방추적법은 단 한번의 전체계산을 통하여 수행가능하지만 후방추적법의 경우에는 매 출력시간마다 초기시점까지 반복해서 계산해야하는 단점이 있다. 경계처리에 있어서 입자추적모형과 상이한 방법을 사용하는 퍼프모형은 폐경계에서는 입자추적모형과 동일한 결과를 나타내지만 개경계에서는 다른 결과를 나타내었다. 또한 오염원이 임의의 공간적 분포를 갖는 경우, 퍼프모형은 입자추적모형보다는 적은 수의 퍼프를 사용할 수 있지만 이에 따른 경계면에서의 수치오차를 발생하였다. 본 연구에서는 개발된 퍼프모형을 검증하고 장${\cdot}$단점을 분석하기 위하여 흐름이 일정한 경우와 전단흐름의 경우에 대하여 이송${\cdot}$확산 수치모의를 수행하였으면, 이를 각각의 경우의 해석해 결과와 비교${\cdot}$분석하였다. 후방추적 퍼프모형은 전방추적 퍼프모형에 비하여 사용된 퍼프수와 관계없이 작은 오차를 발생하였으며, 전체적으로 퍼프 모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 향후, 보다 다양한 흐름영역에서 장${\cdot}$단점 분석 및 오차해석을 수행한 후에 각각의 Lagrangian 모형의 장점만을 갖는 모형결합 방법을 제시할 수 있을 것으로 판단된다.

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

본 연구에서는 수치모델을 이용하여 대하천서 발생되는 조류의 공간적 농도 분포를 예측하였고, 현장실험을 통해 모델을 검증하였다. 국내하천은 다수의 지류가 본류로 유입됨에 따라 오염물질의 생산과 공급이 지속적으로 발생하고, 하천의 유로연장과 하폭에 비해 수심이 낮은 지형학적 특성을 지닌다. 따라서 지류 유입 이후 발생되는 조류의 거동 특성을 분석하기 위해 수심 적분된 2차원 이송-확산 모델을 사용하였다. 광합성 성장을 이루는 조류의 성장속도 계산을 위해 영양염류, 수온, 일사량과 수심 등을 변수로 하는 성장속도 함수들을 위의 모델과 결합하였다. 본 연구의 대상구간은 낙동강과 금호강 합류부를 포함한 강정고령보 하류 약 9.2 km 구간으로 모델 검증을 위한 현장실험을 수행하였다. 2차원 이송-확산 모델의 입력 값인 유속 및 수심을 계산하는 수리동역학 모델 검증을 위해 미국 Sontek사의 M9을 이용하여 낙동강과 금호강 각각 32개, 12개 측선에 대하여 수리량을 측정하였다. 수리량 측정결과, 금호강과 낙동강의 평균 유량은 각각 $240m^3/s$, $60m^3/s$로 측정되었고 측정된 유량을 모델의 상류단 경계조건으로 사용하여 측정 유속 및 수심과 유사한 결과를 모델로부터 취득할 수 있었다. 조류 농도 측정을 위해 독일 bbe사의 AlgaeTorch 10을 사용하였으며, 수리량 측정과 동일한 측선서 총 조류 세포수(cells/ml)를 측정하였다. 농도 측정결과, 하류로 내려감에 따라 조류의 농도가 증가하는 경향이 나타났고 금호강 합류 후 최대농도는 측정구간 최하류 우안서 4,460 cells/ml로 나타났다. 주 흐름이 발생하는 하천 중앙부에 비해 유속이 느린 하안서 상대적으로 높은 농도가 측정되었으며, 이와 같은 경향은 하류로 내려감에 따라 강하게 나타났다. 측정된 조류 농도를 이용한 2차원 이송-확산 모델 검증결과, 합류부 최상류 측선서 MAPE = 10.5 %의 최대오차가 발생하였고 최하류 측선서 MAPE = 6.7 %의 최소오차가 발생하였다. 인과 질소와 같은 영양염류의 농도가 높고 횡 방향 수온 분포가 균일한 대상구간의 특성상 영양염류 함수와 수온 함수로부터 계산된 성장속도 가중치 범위는 각각 0.8~1.0, 0.91~1.09로 공간적 변동성이 크게 나타나지 않은 반면, 수심을 변수로 하는 일사량 함수의 성장속도 가중치 범위는 0.05~1.00으로 상대적으로 매우 높은 공간적 변동성이 나타났다. 수심이 4 m 이하인 하천 양안서 0.8 이상의 가중치가 나타났으며, 수심이 7 m 이상인 하천 중앙서 0.4 이하의 가중치가 나타났다. 본 연구의 수치모의 결과, 수리동역학 모델로부터 계산된 수심이 모델 결과 값에 큰 영향을 미치는 것으로 판단된다.