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 temperature 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 differential equations in the time coordinate are efficiently obtained by a computational algorithm based upon the concept of Lie series. The method is particularly easy when the heat transfer coefficient is a weak function of the coolant speed, as for liquid metals.
ASJC Scopus subject areas
- Nuclear and High Energy Physics
- Nuclear Energy and Engineering
- Materials Science(all)
- Safety, Risk, Reliability and Quality
- Waste Management and Disposal
- Mechanical Engineering