초록
원자력 분야에서 중성자 수송계산을 위해 개발되어 널리 사용되는 방향차분법을 시간 종속 복사 열전달식의 해를 구하는데 적용하였다. 광자의 방향별 밀도를 자체수반형 2계 편미분방정식으로 나타내어 해의 안정성을 높였고 매질의 온도방정식의 비선형성은 다단계 선형화법을 사용하여 근사하였다. 본 연구에서 개발된 해법을 전형적인 Marshak wave 문제에 적용하였고 계산 결과를 기존의 Monte Carlo의 계산결과와 비교하여 그 우월성을 보였다.
In this study, the discrete ordinates method which has been widely used in the solution of neutron transport equation is applied to the solution of the time dependent radiative transfer equation. The self-adjoint form of the second order radiation intensity equation is used to enhance the stability of the solution, and a new multi-step linearization method is developed to avoid the nonlinearity in the material temperature equation. This new solution method is applied to the well known Marshak wave problem, and the numerical result is compared with that of the conventional Monte-Carlo method.
