Taylor’s series method for solving the nonlinear point kinetics equations

Abstract Taylor’s series method for solving the point reactor kinetics equations with multi-group of delayed neutrons in the presence of Newtonian temperature feedback reactivity is applied and programmed by FORTRAN. This system is the couples of the stiff nonlinear ordinary differential equations. This numerical method is based on the different order derivatives of the neutron density, the precursor concentrations of i -group of delayed neutrons and the reactivity. The r th order of derivatives are derived. The stability of Taylor’s series method is discussed. Three sets of applications: step, ramp and temperature feedback reactivities are computed. Taylor’s series method is an accurate computational technique and stable for negative step, negative ramp and temperature feedback reactivities. This method is useful than the traditional methods for solving the nonlinear point kinetics equations.

[1]  A. Aboanber,et al.  Generalization of the analytical inversion method for the solution of the point kinetics equations , 2002 .

[2]  Y. Chao,et al.  A Resolution of the Stiffness Problem of Reactor Kinetics , 1985 .

[3]  A. Aboanber Analytical solution of the point kinetics equations by exponential mode analysis , 2003 .

[4]  W. Stacey Nuclear Reactor Physics , 2001 .

[5]  Abdallah A. Nahla,et al.  Generalization of the analytical exponential model to solve the point kinetics equations of Be-and D2O-moderated reactors , 2008 .

[6]  T. Sathiyasheela Power series solution method for solving point kinetics equations with lumped model temperature and feedback , 2009 .

[7]  Abdallah A. Nahla Analytical solution to solve the point reactor kinetics equations , 2010 .

[8]  D. Hetrick,et al.  Dynamics of nuclear reactors , 1972 .

[9]  A. Aboanber Stability of generalized Runge-Kutta methods for stiff kinetics coupled differential equations , 2006 .

[10]  A. Aboanber,et al.  Solution of the point kinetics equations in the presence of Newtonian temperature feedback by Padé approximations via the analytical inversion method , 2002 .

[11]  Haofeng Li,et al.  A new integral method for solving the point reactor neutron kinetics equations , 2009 .

[12]  Barbaro Quintero-Leyva CORE: A numerical algorithm to solve the point kinetics equations , 2008 .

[13]  J. Sánchez On the numerical solution of the point reactor kinetics equations by generalized Runge-Kutta methods , 1989 .