초록
본 연구에서는 수치해석의 시간차분에서 발생하는 절단오차를 극복하는 방안으로 고유치 기법을 도입하였다. 고유치기법은 모의를 할 때 공간만을 이산화하는 특징을 가지며, 공간적으로 이산화된 방정식을 대각화시킴으로써 선형동력학적 시스템을 분리시킨 후 시간적분을 이용한 계산이 임의의 위치에서 임의의 시간에 대해 개별적으로 또 연속적으로 수행된다. 이러한 고유치기법을 이용하여 오염물 이동을 모의하고 이를 해석해와 비교 검증하였고, 동일한 조건에서 유한요소법을 이용한 수치모형과 고유치 기법을 이용한 용질이동의 예측을 실시한 결과 고유치기법을 이용할 경우 계산시간과 저장용량이 수치모형에 비해 절약됨을 확인할 수 있었다. 고유치 기법을 이용하여 지하유류저장 공동주위의 불균일 유속장에서 용질의 이동을 분석하였다 이 방법이 모의발생에 오랜 시간이 걸리는 문제에 유용하게 사용될 수 있으므로, 공동에 인접한 오염원으로부터 공동의 안전성을 평가하기 위한 민감도 분석에 이 방법을 적용하였으며, 모의결과에 의하면, 종분산지수와 횡분산지수가 각각 50 m, 5 m일 때 공동에 도달하는 시간은 약 50년으로 추정되었다.
The eigenvalue technique is introduced to overcome the problem of truncation errors caused by temporal discretization of numerical analysis. The eigenvalue technique is different from simulation in that only the space is discretized. The spatially discretized equation is diagonized and the linear dynamic system is then decoupled. The time integration can be done independently and continuously for any nodal point at any time. The results of eigenvalue technique are compared with the exact solution and FEM numerical solution. The eigenvalue technique is more efficient than the FEM to the computation time and the computer storage in the same conditions. This technique is applied to the solute transport analysis in nonuniform flow fields around underground storage caverns. This method can be very useful for time consuming simulations. So, a sensitivity analysis is carried out by using this method to analyze the safety of caverns from nearly located contaminant sources. According to the simulations, the reaching time from source to the nearest cavern may takes 50 years with longitudinal dispersivity of 50 m and transversal dispersivity of 5 m, respectively.