초록
본 연구는 오염물 거동에 대한 수치해석을 위해 보편적으로 사용되고 있는 수치 방법들의 장단점을 총괄적으로 나타내고, 효율적인 수치모델링 기법 개발을 위해 ELLAM과 LEZOOMPC를 비교분석하였다. 지하수 분야에서 가장 많이 사용되는 수치 방법은 Eulerian-Lagrangian 방식과 Eulerian 방식인데, Eulerian-Lagrangian 방식은 수치영역 내에서 일반적으로 질량을 보존하지 못하고, 경계조건을 체계적으로 처리하지 못하는 한계를 갖고 있다. 반면에 Eulerian 빙식은 시간 및 공간 절삭 오차로 인해서 시간 간격 및 격자 크기를 극히 줄여야 하는 제약을 갖고 있다. 최근 10 년간 지하수 분야에서 크게 대두되고 있는 수치기법인 ELLAM(Eulerian Lagrangian Localized Adjoint Method)은 Eulerian-Lagrangian 방식과 Eulerian 방식에서 나타나는 수치 제약점이나 한계점을 동시에 해결하는 수치기법으로 알려져 왔다. 그러나 본 연구에서는 ELLAM의 장단점을 파악하고 보완점을 제안한다. ELLAM의 단점을 파악하기 위해, mesh Peclet number가 다른 예제들을 설정하고, 그 예제들에 대한 ELLAM, LEZOOMPC(Lagrangian-Eulerian ZOOMing Peak and valley Capturing)와 visual MODFLOW의 수치결과들을 해석해와 비교하였다. Mesh Peclet number가 무한대일 때 ELLAM의 수치결과는 수치진동으로 인해 해석해와 일치하지 않았으나, LEZOOMPC의 수치 결과는 해석해와 일치했다. 위의 결과는 ELLAM의 수치오차가 LEZOOMPC의 특성을 이용하여 개선 및 보완될 수 있는 가능성을 시사해 준다. 따라서 ELLAM에 LEZOOMPC의 후향 입지추적, 전향 입지추적, 선택적 국부 격자 세립화 과정과 최고/최저 농도점 이동 추적 과정을 결합하면 ELLAM의 수치적 장점을 유지하면서 mesh Peclet number에 제약을 받지 않는 효율적인 수치모델링 기법을 개발할 수 있을 것으로 판단된다.
This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.