A new method is developed for computing transient temperature fields in cylindrical reactor channels when the relative fuel power density, the coolant speed, and the channel inlet tem¬perature are arbitrarily specified functions of time. Through the use of Hankel transforms, numerical instabilities inherent in finite difference approximations of spatial derivatives are avoided. Solutions of high order systems of ordinary differen¬tial equations in the time coordinate are efficiently obtained by a computational algorithm based upon the concept of Lie series.
展开▼