Método de Newton-Raphson

El método de Newton-Raphson es el algoritmo iterativo estándar de la industria para resolver las ecuaciones no lineales de flujo de potencia en redes eléctricas. Funciona linealizando las ecuaciones de desbalance en cada iteración mediante la matriz Jacobiana, resolviendo el sistema lineal resultante para obtener las correcciones de tensión, y actualizando las estimaciones de tensión hasta que los desbalances caen por debajo de una tolerancia especificada.

Aspectos clave del método de Newton-Raphson:

• Convergencia cuadrática: Cerca de la solución, el método de Newton-Raphson presenta convergencia cuadrática, lo que significa que el número de dígitos correctos se duplica aproximadamente con cada iteración. Sistemas de transmisión bien condicionados convergen típicamente en 3–5 iteraciones independientemente del tamaño del sistema, lo que supone una ventaja práctica importante sobre métodos más simples como Gauss-Seidel.
• Matriz Jacobiana: En cada iteración el método forma la Jacobiana, una matriz cuadrada de derivadas parciales de las ecuaciones de desbalance de potencia respecto a las magnitudes y ángulos de tensión. La Jacobiana captura la sensibilidad de las inyecciones de potencia ante cambios de tensión y debe recalcularse (o al menos actualizarse) en cada paso.
• Álgebra lineal dispersa: Dado que las matrices Y-bus y Jacobiana son altamente dispersas (la mayoría de elementos fuera de la diagonal son cero), los solvers de producción emplean técnicas de matrices dispersas (ordenamiento óptimo, factorización LU simbólica y numérica) para resolver eficientemente las ecuaciones de corrección lineal, haciendo el método escalable a sistemas con decenas de miles de barras.
• Condiciones iniciales: El método requiere un punto de partida razonable (el arranque plano con todas las tensiones a 1,0∠0° es habitual). Condiciones iniciales deficientes, o sistemas cercanos al límite de estabilidad de tensión, pueden provocar divergencia, por lo que las implementaciones robustas incluyen limitación del tamaño de paso, acotación de magnitudes de tensión y heurísticas de ajuste de ángulos.
• Variantes: Se han desarrollado varias variantes para necesidades específicas: el método Desacoplado Rápido usa Jacobianas constantes y simplificadas para mayor velocidad; el flujo de potencia DC linealiza completamente el problema; y el flujo de potencia de continuación traza la curva P-V parametrizando el nivel de carga, útil para el análisis de estabilidad de tensión cerca de puntos de colapso.