본 발명은 지하수 내의 오염물 거동을 수치적으로 해석하기 위한 방법에 관한 것으로, 본 발명에 따르면, 높은 페클릿(Peclet) 숫자에 대하여만 해석이 가능했던 기존의 라그랑지안-율러리안 기법의 문제점을 해결하여, 총유량 경계조건에 인접한 요소의 경우는 율러리안 기법을 적용하고, 내부의 요소의 경우는 기존의 라그랑지안-율러리안 기법을 병합하여 새로운 수치 행렬을 재구성하여 행렬을 푸는 것에 의해, 총유량 경계조건하에서 높은 페클릿(Peclet) 숫자는 물론 특별히 낮은 페클릿(Peclet) 숫자와 같은 다양한 페클릿(Peclet) 숫자에 대한 오염물 거동 문제에 대하여도 정확한 수치해석이 가능하도록 구성되는 개선된 라그랑지안-율러리안 기법을 이용한 총유량 경계조건을 가지는 오염물 거동의 분석방법이 제공된다.