Abstract
A numerical model has been developed by employing a finite element method to simulate the depth-averaged 2-D dispersion of the heat pollutant, which is an important pollutant material in natural streams. Among the finite element methods, the Streamline Upwind/Petrov Galerkin (SUPG) method was applied. Also both linear and quadratic elements can be applied so that irregular river boundaries can be easily represented. To show the movement of heat pollutants, the reaction term describing heat transfer was represented as an equation in which sink/source term is proportional to the difference between the equilibrium temperature and water surface temperature. The equation was expressed so that the water surface temperature changes according to the temperature transfer coefficient and the equilibrium temperature. For the calibration of the model developed, analytic and numerical results from a case of rectangular channel with full width continuous injection have been compared in a steady state. The comparisons showed that the numerical results were in good agreement with analytical solutions. The application site was selected from the downstream of Paldang dam to Jamsil submerged weir, and overall length of this site is about 22.5 km. The change of water temperature caused by the discharge from the Guri sewage treatment plant has been simulated, and results were similar to the observed data. Overall it is concluded that the developed model can represent the water temperature changes due to heat transport accurately. But the verification using observed data will further enhance the validity of the model.
본 논문에서는 하천에 유입되는 열 오염물질의 혼합거동을 분석하기 위해 2차원 수심 평균된 이송-분산방정식에 유한요소법을 적용하여 수치 모형을 개발하였다. 유한요소법의 여러 수치기법 중 SUPG법을 적용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각 요소망의 혼용이 가능하도록 하였다. 열 오염물질의 거동을 표현하기 위해서 열 교환을 묘사하는 반응항을 평형온도와 수온과의 차이에 비례하는 식으로 나타내고, 열교환 계수 및 평형온도에 따라 수온의 변화가 적용되도록 방정식을 구성하였다. 모형의 검증을 위하여 직사각형 수로에 선원으로 연속주입하여 얻은 수치해와 1차원 정상상태의 해석해를 비교하였다. 비교결과 수치해와 해석해의 결과가 서로 일치하는 것으로 밝혀졌다. 모형의 현장적용을 위해 상수원 보호구역인 팔당댐 하류부터 잠실수중보까지 22.5 km 구간을 대상영역으로 하였다. 구리하수처리장 방류수에 의한 수온 변화를 모의한 결과 수질측정망 측정자료와 비교적 비슷한 경향이 나타났다. 본 연구에서 개발한 수치모형이 열 오염원 유입으로 인한 수온 변화를 잘 표현하는 것을 알 수 있었다.