본 연구에서는 비구조격자 좌표계상에서 천수방정식 해석을 위한 수치모형을 개발하였다. Fractional step method를 이용하여 이차원의 천수방정식을 두개의 일차원 문제로 분리하여 계산 효율이 우수한 수치연산을 수행하였다. 분리된 일차원 문제는 이차정확도의 TVD 기법을 이용하여 해석하였다. TVD 기법은 HLLC 기법을 이용하여 Riemann 해를 구한후 WAF 기법을 기반으로 이용하였으며 WAF 흐름율 제한자를 이용하여 이차정확도 문제에 수반