Método para la resolución de un sistema de ecuaciones diferenciales por transformación de subsistemas

Abstract

El análisis del comportamiento dinámico de una central nuclear implica resolver un conjunto do ecuaciones diferenciales que tienen en cuenta tanto la evolución cinética del reactor, como los procesos termohidráulicos de los sistemas vinculados con aquel. Cabe señalar que estos últimos sistemas presentan la particularidad de poseer una evolución temporal sensiblemente más lenta. Por lo tanto, al utilizar cualquiera de los métodos numéricos para ecuaciones diferenciales, se observaría que la solución de un sistema, con constantes de tiempo muy diferentes, resultaría correcta si se usa, únicamente, un paso de integración suficientemente pequeño. Pese a que el tamaño de dicho paso es realmente aceptable para las ecuaciones cinéticas, este es exageradamente chico para la parte térmica del modelo; con lo cual, el insumo de tiempo de cómputo se eleva notablemente. Lo que aquí se presenta es un nuevo algoritmo basado en el método de RUNG-KUTTA de 4° orden pero con un doble paso de integración.