본 논문에서는 Eulerian 좌표축에서 온도 경계층과 혼합된 비정상 경계층 유 동을 효과적으로 다룰 수 있는 좌표변환 방법과 수치해석 방법을 도입하였다. 이전 의 수치적 방법들은 축차적이며, 뒷정체점 부근의 경계층 두께가 시간과 더불어 지수 적으로 증가하므로, 격자점의 수를 경계층 두께의 증가에 따라 상당히 늘려 주어야 된 다. 그러나 여기에서는, 격자점의 수를 무리하게 늘릴 필요가 없으며 일반적인 비선 형 차분방정식을 정확도를 떨어뜨리지 않고 선형화시킴으로써 비축차적인 수치해들이 허용된다. 이런 선형화 방법은 Beam & Warming에 의해 최초로 압축성 Navier-Stokes 방정식에 사용되었고, Orlandi & Ferziger, Kim & Chang에 의해 경계층 유동에 확장되 었다.또한 뒷정체점 부근에서 경계층 두께의 증가로 인한, 격자점 증가의 필요를 피하기 위하여, 몇 가지 서로 다른 종류의 변환변수들을 시간과 공간에 따라 선별적으 로 사용하여 수치적인 경계층 두께가 거의 일정하도록 만들었다. 이와 같은 변환변 수들은 또한 쌍-포물선형인 현 지배방정식의 초기조건들을 쉽게 구할 수 있도록 허용 해 준다.
A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).
