An Algorithm to Optimize the Calculation of the Fourth Order Runge-Kutta Method Applied to the Numerical Integration of Kinetics Coupled Differential Equations  

作　　者：Sadao Isotani Walter Maigon Pontuschka Seiji Isotani 

机构地区：[1]Institute of Physics,University of Sao Paulo,Sao Paulo,Brazil [2]Institute of Mathematics and Computer Science,University of Sao Paulo,São Carlos,Brazil

出　　处：《Applied Mathematics》2012年第11期1583-1592,共10页应用数学（英文）

摘　　要：The kinetic electron trapping process in a shallow defect state and its subsequent thermal- or photo-stimulated promotion to a conduction band, followed by recombination in another defect, was described by Adirovitch using coupled rate differential equations. The solution for these equations has been frequently computed using the Runge-Kutta method. In this research, we empirically demonstrated that using the Runge-Kutta Fourth Order method may lead to incorrect and ramified results if the numbers of steps to achieve the solutions is not “large enough”. Taking into account these results, we conducted numerical analysis and experiments to develop an algorithm that determines the smallest non-critical number of steps in an automatic way to optimize the application of the Runge-Kutta Fourth Order method. This algorithm was implemented and tested in a variety of situations and the results have shown that our solution is robust in dealing with different equations and parameters.The kinetic electron trapping process in a shallow defect state and its subsequent thermal- or photo-stimulated promotion to a conduction band, followed by recombination in another defect, was described by Adirovitch using coupled rate differential equations. The solution for these equations has been frequently computed using the Runge-Kutta method. In this research, we empirically demonstrated that using the Runge-Kutta Fourth Order method may lead to incorrect and ramified results if the numbers of steps to achieve the solutions is not “large enough”. Taking into account these results, we conducted numerical analysis and experiments to develop an algorithm that determines the smallest non-critical number of steps in an automatic way to optimize the application of the Runge-Kutta Fourth Order method. This algorithm was implemented and tested in a variety of situations and the results have shown that our solution is robust in dealing with different equations and parameters.

关 键 词：Coupled Differential Equations RUNGE-KUTTA Kinetics Equations 

分 类 号：O1[理学—数学]