⚡ Aula 03: Fluxo de Potência (Load Flow) - Formulação, Gauss-Seidel, Newton-Raphson e Compensação Reativa

Integração Curricular & Contexto

• Disciplina: Planejamento e Análise de Sistemas de Energia Elétrica | Hub Principal
• Material Didático Oficial: Complemento analítico e aprofundamento da apostila Fluxo de Potência.
• Aula Anterior: 📐 Aula 02: Representação do SEP, Sistema por Unidade (pu) e Modelagem de LTs
• Próxima Aula: 💥 Aula 04: Análise de Curto-Circuito no SEP - Simétrico e Assimétrico
• Objetivo da Aula: Compreender a formulação matemática não linear do fluxo de carga em regime permanente, a classificação operacional das barras do SEP ($PQ$, $PV$ e Slack), os métodos numéricos iterativos de solução (Gauss-Seidel, Newton-Raphson e Desacoplado Rápido), a formulação linearizada (DC Power Flow) e as estratégias físicas de controle de tensão e compensação de reativos.

---

📌 Sumário da Aula

1. 1. O Problema Fundamental do Fluxo de Potência em Regime Permanente
   • Leis Físicas Subjacentes:
2. 2. Classificação Nodal e Variáveis de Estado do SEP
   • 2.1 Barras de Carga (PQ)
   • 2.2 Barras de Geração ou Tensão Controlada (PV)
   • 2.3 Barra de Folga, Referência ou Balanço (Slack / Swing)
   • 2.4 Contabilidade de Variáveis e Graus de Liberdade
3. 3. Diagrama Arquitetural Unificado do Fluxo de Potência
4. 4. Montagem Sistemática da Matriz de Admitância Nodal ($\mathbf{Y}_{bus}$)
   • 4.1 Regras de Inspeção Direta
   • 4.2 Equações Algébricas das Injeções Nodais de Potência
5. 5. Métodos Numéricos de Solução do Fluxo de Potência
   • 5.1 Método de Gauss-Seidel
   • 5.2 Método de Newton-Raphson
   • 5.3 Método Desacoplado e Desacoplado Rápido (FDLF)
   • 5.4 Fluxo de Potência Linearizado (DC Power Flow)
6. 6. Controle de Tensão e Compensação Reativa no SEP
   • 6.1 O Acoplamento Forte Q-V em Redes de Alta Tensão
   • 6.2 Bancos de Capacitores e Reatores Shunt
   • 6.3 Máquinas e Condensadores Síncronos
   • 6.4 Transformadores com Comutação de Tap sob Carga (LTC/OLTC) e FACTS
7. 7. 📝 Estudo de Caso Numérico Resolvido Passo a Passo
   • Formulação do Problema:
   • Passo 1: Montagem da Matriz de Susceptância Nodal ($\mathbf{B}_{bus}$)
   • Passo 2: Redução Matricial para a Barra de Referência
   • Passo 3: Solução Analítica dos Ângulos de Barra
   • Passo 4: Cálculo dos Fluxos de Potência Ativa nos Ramos
8. 8. 📚 Questões Resolvidas das Materiais Teóricos de Referência
   • Questão 01 — Variáveis e Tipos de Barra no SEP
   • Questão 02 — Formulação Básica e Leis Físicas do Fluxo
   • Questão 03 — Modelagem Algébrica e Graus de Liberdade
   • Questão 04 — Métodos de Solução e Linearização do SEP
   • Questão 05 — Solução Linear de Sistema de 3 Barras
   • Questão 06 — Modelagem Não Linear de Newton-Raphson
   • Questão 07 — Fatores Influenciadores do Perfil de Tensão
   • Questão 08 — Técnicas de Compensação Reativa e Máquinas Síncronas
9. 9. Resumo Conceitual da Aula e Conexões Curriculares

---

1. O Problema Fundamental do Fluxo de Potência em Regime Permanente

O estudo de Fluxo de Potência (ou Load Flow) constitui a principal ferramenta de simulação em regime permanente senoidal da Engenharia de Sistemas de Potência. Seu objetivo precípuo é calcular o estado elétrico estático de uma rede interligada sob uma dada condição de carregamento e geração.

Matematicamente, determinar o estado do sistema significa calcular os fasores de tensão nodais em todas as barras da rede: ${\mathbf{V}}_k=|V_k|\angle {\theta }_k,k=1,2,\dots ,N$

Uma vez conhecidos o módulo $|V_k|$ e o ângulo de fase ${\theta }_k$ de cada nó, torna-se possível derivar analiticamente todas as demais grandezas operacionais do SEP:

1. Fluxos de potência ativa e reativa em todas as linhas de transmissão e transformadores ($P_{km},Q_{km}$);
2. Perdas elétricas por efeito Joule e dispersão magnética em cada ramo do sistema ($P_{loss},Q_{loss}$);
3. Injeções de potência líquida nas usinas e subestações;
4. Correntes transitando em condutores e barramentos ($I_{km}=(V_k-V_m)y_{km}$).

Leis Físicas Subjacentes:

A rede elétrica obedece à Primeira Lei de Kirchhoff (Lei Nodal das Correntes), expressa em coordenadas matriciais pelo produto da matriz de admitância nodal ${\mathbf{Y}}_{bus}$ pelo vetor de tensões: ${\mathbf{I}}_{bus}={\mathbf{Y}}_{bus}{\mathbf{V}}_{bus}$

A injeção complexa de potência líquida na barra $k$ é dada pelo produto da tensão fasorial pelo conjugado da corrente líquida injetada: $S_k=P_k+j{Q}_k=V_k{I}_k^{\ast }=(P_{Gk}-P_{Ck})+j(Q_{Gk}-Q_{Ck})$

Substituindo a expressão de corrente $I_k={\sum }_{m=1}^N{Y}_{km}{V}_m$ na equação da potência: $S_k^{\ast }=P_k-j{Q}_k=V_k^{\ast }{\sum }_{m=1}^N{Y}_{km}{V}_m$

Essa equação estabelece um sistema algébrico de equações estritamente não lineares, visto que a potência é proporcional ao produto quadrático entre módulos de tensão e funções trigonométricas de diferenças angulares.

---

2. Classificação Nodal e Variáveis de Estado do SEP

Em cada barra $k$ de um sistema elétrico de $N$ barras, existem 4 grandezas elétricas fundamentais:

1. Potência ativa líquida injetada: $P_k=P_{Gk}-P_{Ck}$
2. Potência reativa líquida injetada: $Q_k=Q_{Gk}-Q_{Ck}$
3. Módulo da tensão nodal: $|V_k|$
4. Ângulo de fase da tensão nodal: ${\theta }_k$

Para que o sistema matemático seja possível e determinado, exatamente 2 dessas variáveis devem ser especificadas previamente (dados de entrada), enquanto as outras 2 variáveis tornam-se incógnitas a serem determinadas numericamente.

Tipo de Barra	Denominação Típica	Variáveis Pré-Fixadas (Conhecidas)	Variáveis a Calcular (Incógnitas)	Função Operacional no SEP
PQ	Barra de Carga	$P_k$, $Q_k$	$\Vert V_k\Vert$, ${\theta }_k$	Atendimento a cargas e subestações sem geração
PV	Barra de Geração (Tensão Controlada)	$P_k$, $\Vert V_k\Vert$	$Q_k$, ${\theta }_k$	Usinas geradoras com AVR e compensadores síncronos
Vθ / Slack	Barra de Folga, Referência ou Swing	$\Vert V_k\Vert$, ${\theta }_k=0^{\circ }$	$P_k$, $Q_k$	Fechamento do balanço de perdas e referência angular

graph LR
    subgraph Classificacao_Barras ["Classificação das Barras no SEP"]
        PQ["Barra PQ (Carga)<br/>Conhece: P, Q<br/>Calcula: |V|, θ"]
        PV["Barra PV (Geração)<br/>Conhece: P, |V|<br/>Calcula: Q, θ"]
        Slack["Barra Slack (Folga / Vθ)<br/>Conhece: |V|, θ = 0°<br/>Calcula: P, Q"]
    end
### 2.1 Barras de Carga (PQ)
Representam os nós aos quais estão conectadas cargas industriais, comerciais e residenciais, ou subestações puramente de manobra e distribuição sem equipamento de geração síncrona.
- Os consumos $P_{Ck}$ e $Q_{Ck}$ são conhecidos pela previsão de carga;
- Não há injeção de geração ($P_{Gk} = 0, Q_{Gk} = 0$), logo $P_k = -P_{Ck}$ e $Q_k = -Q_{Ck}$;
- O algoritmo de fluxo de carga deve determinar a magnitude da tensão $|V_k|$ e sua defasagem $\theta_k$.

### 2.2 Barras de Geração ou Tensão Controlada (PV)
Representam as usinas hidroelétricas, termelétricas e grandes parques de geração centralizada equipados com reguladores automáticos de tensão (AVR - *Automatic Voltage Regulator*).
- O operador programa a potência ativa despachada $P_k$ (definida pelo despacho econômico) e o valor de tensão desejado no barramento $|V_k|$;
- As incógnitas são o ângulo fasorial $\theta_k$ e a potência reativa $Q_k$ necessária para sustentar a tensão especificada;
- **Restrição Física de Reativos:** O gerador síncrono possui limites térmicos de enrolamento de campo e estator:
  $$Q_{k}^{min} \le Q_k \le Q_{k}^{max}$$
  Se durante o processo iterativo a potência reativa calculada violar um limite ($Q_k > Q_k^{max}$ ou $Q_k < Q_k^{min}$), o limite correspondente é fixado ($Q_k = Q_k^{lim}$) e a barra é **automaticamente convertida em uma barra PQ** na iteração subsequente, liberando a magnitude da tensão $|V_k|$ para flutuar.

### 2.3 Barra de Folga, Referência ou Balanço (Slack / Swing)
Em qualquer rede de transmissão, ocorrem perdas elétricas inerentes por efeito Joule ($\sum R I^2$). Contudo, essas perdas dependem do fluxo de corrente, que é função direta das diferenças de tensão e de fase entre as barras — as quais são justamente as variáveis desconhecidas antes da convergência do fluxo!
- Dessa forma, **é fisicamente impossível prever antecipadamente a potência ativa total que os geradores devem fornecer para cobrir a carga mais as perdas**.
- Define-se, obrigatoriamente, ao menos **uma barra de folga** (usualmente associada à maior usina hidroelétrica do sistema, com grande capacidade de regulação de velocidade).
- A barra slack tem seu ângulo de fase fixado como referência absoluta do sistema ($\theta_{slack} = 0^\circ$) e sua tensão terminal pré-ajustada (comumente $|V_{slack}| = 1{,}0$ ou $1{,}05\text{ pu}$).
- Após a convergência, a barra slack supre o déficit exato de potência ativa ($P_{slack}$) e reativa ($Q_{slack}$) necessário para fechar o balanço de potência e equilibrar as perdas totais do SIN.

### 2.4 Contabilidade de Variáveis e Graus de Liberdade
Para um sistema interligado de $N$ barras contendo $N_{PV}$ barras de geração e 1 única barra de folga:
- Quantidade total de barras de carga $PQ$: $N_{PQ} = N - 1 - N_{PV}$
- Vetor de variáveis de estado a determinar:
  - $N - 1$ ângulos de fase: $\theta_k$ para todas as barras exceto a Slack;
  - $N_{PQ}$ magnitudes de tensão: $|V_k|$ para todas as barras do tipo $PQ$.
- **Número Total de Equações Algébricas Não Lineares a Resolver:**
  $$\text{Dimensão do Sistema} = 2(N - 1) - N_{PV} = (N - 1) + N_{PQ}$$

---

## 3. Diagrama Arquitetural Unificado do Fluxo de Potência

O diagrama técnico vetorial abaixo sintetiza a montagem da matriz de admitância nodal, a classificação física das barras e a iteração matemática do método de Newton-Raphson:

![[sep_fluxo_potencia_newton_raphson_ybus.svg]]

---

## 4. Montagem Sistemática da Matriz de Admitância Nodal ($\mathbf{Y}_{bus}$)

A matriz $\mathbf{Y}_{bus} \in \mathbb{C}^{N \times N}$ encapsula integralmente a topologia da rede de transmissão, as admitâncias longitudinais dos ramos e os elementos capacitivos/reativos em derivação (*shunt*).

### 4.1 Regras de Inspeção Direta
Para cada ramo conectado entre as barras $k$ e $m$ com impedância série $z_{km} = r_{km} + jx_{km}$, a admitância série vale:
$$y_{km} = \frac{1}{z_{km}} = g_{km} + jb_{km} = \frac{r_{km}}{r_{km}^2 + x_{km}^2} - j\frac{x_{km}}{r_{km}^2 + x_{km}^2}$$

Se a linha possuir capacitância distribuída modelada em circuito $\pi$ com admitância *shunt* total $y_{sh,km} = jb_{sh,km}$, cada extremidade recebe metade da susceptância: $y_{km0} = j\frac{b_{sh,km}}{2}$.

As entradas da matriz $\mathbf{Y}_{bus}$ são estabelecidas diretamente por duas regras fundamentais:

> [!important] Regras de Formação da $\mathbf{Y}_{bus}$
> 
> 1. **Termos Diagonais ($Y_{kk}$ - Auto-Admitância da Barra $k$):**
>    É igual à soma de todas as admitâncias diretamente conectadas à barra $k$, incluindo as admitâncias longitudinais de todas as linhas incidentes, as admitâncias *shunt* das linhas conectadas ao nó e quaisquer elementos em derivação na própria barra ([[Bancos Automáticos de Capacitores com Indutores Dessintonizados Anti-Ressonância|bancos de capacitores]] ou reatores):
>    $$Y_{kk} = y_{k0}^{sh} + \sum_{m \in \Omega_k} \left( y_{km} + \frac{y_{sh,km}}{2} \right)$$
> 2. **Termos Fora da Diagonal ($Y_{km}$ - Admitância Mútua entre Barra $k$ e Barra $m$):**
>    É igual ao negativo da soma das admitâncias de linha diretamente conectadas entre a barra $k$ e a barra $m$:
>    $$Y_{km} = Y_{mk} = -y_{km}$$
>    Se não houver linha de transmissão conectando diretamente a barra $k$ à barra $m$, tem-se rigorosamente:
>    $$Y_{km} = 0$$

Em redes reais de centenas ou milhares de barras, cada nó conecta-se a poucas linhas adjacentes (em média $2$ a $4$ conexões). Portanto, a matriz $\mathbf{Y}_{bus}$ é **extremamente esparsa** (mais de $98\%$ dos seus elementos são nulos), propriedade explorada por algoritmos computacionais avançados de fatoração LU esparsa e ordenação nodal (algoritmos de Tinney).

### 4.2 Equações Algébricas das Injeções Nodais de Potência
Representando a admitância nodal na forma retangular $Y_{km} = G_{km} + jB_{km}$ e as tensões fasoriais em coordenadas polares $V_k = |V_k|e^{j\theta_k}$ e $V_m = |V_m|e^{j\theta_m}$:
$$S_k^* = P_k - jQ_k = |V_k|e^{-j\theta_k} \sum_{m=1}^N (G_{km} + jB_{km}) |V_m|e^{j\theta_m}$$
$$S_k^* = \sum_{m=1}^N |V_k||V_m| (G_{km} + jB_{km}) e^{-j(\theta_k - \theta_m)}$$

Aplicando a identidade de Euler $e^{-j\theta_{km}} = \cos\theta_{km} - j\sin\theta_{km}$ (onde $\theta_{km} \equiv \theta_k - \theta_m$) e separando a parte real e imaginária:

$$\boxed{P_k(\boldsymbol{\theta}, |\mathbf{V}|) = |V_k| \sum_{m=1}^N |V_m| \big( G_{km} \cos(\theta_k - \theta_m) + B_{km} \sin(\theta_k - \theta_m) \big)}$$

$$\boxed{Q_k(\boldsymbol{\theta}, |\mathbf{V}|) = |V_k| \sum_{m=1}^N |V_m| \big( G_{km} \sin(\theta_k - \theta_m) - B_{km} \cos(\theta_k - \theta_m) \big)}$$

Essas são as **equações fundamentais do fluxo de potência**, válidas para qualquer topologia e nível de carregamento.

---

## 5. Métodos Numéricos de Solução do Fluxo de Potência

### 5.1 Método de Gauss-Seidel
O método de Gauss-Seidel é um método iterativo de relaxação sucessiva (ponto fixo). A equação nodal $I_k = \frac{P_k - jQ_k}{V_k^*} = \sum_{m=1}^N Y_{km} V_m$ é isolada para a tensão $V_k$:

$$V_k^{(\nu+1)} = \frac{1}{Y_{kk}} \left[ \frac{P_k^{esp} - jQ_k^{esp}}{(V_k^{(\nu)})^*} - \sum_{m=1}^{k-1} Y_{km} V_m^{(\nu+1)} - \sum_{m=k+1}^N Y_{km} V_m^{(\nu)} \right]$$

- **Vantagens:** Extremamente simples de programar; consome memória desprezível por não armazenar nem inverter matrizes Jacobianas.
- **Desvantagens:** Convergência de ordem linear; exige elevado número de iterações ($50$ a $300$ iterações para redes médias); extremamente sensível à escolha da barra de folga e a sistemas altamente carregados; frequentemente diverge em redes com linhas de baixa impedância ou compensação série.

### 5.2 Método de Newton-Raphson
O método de Newton-Raphson é o padrão-ouro computacional para a solução do fluxo de potência na indústria elétrica mundial. Baseia-se na linearização das equações de injeção por meio da **Série de Taylor multivariada truncada no termo de primeira ordem**.

#### Vetores de Resíduo (Incompatibilidades ou Mismatches):
Para a iteração corrente $\nu$, avalia-se a discrepância entre as potências líquidas especificadas e as potências calculadas pelas equações não lineares:
$$\Delta P_k^{(\nu)} = P_k^{esp} - P_k(\boldsymbol{\theta}^{(\nu)}, |\mathbf{V}|^{(\nu)})$$
$$\Delta Q_k^{(\nu)} = Q_k^{esp} - Q_k(\boldsymbol{\theta}^{(\nu)}, |\mathbf{V}|^{(\nu)})$$

O sistema de correção linearizado é formulado matricialmente pela **Matriz Jacobiana $\mathbf{J}$**:

$$\begin{bmatrix} \Delta \mathbf{P} \\ \hline \Delta \mathbf{Q} \end{bmatrix} = \begin{bmatrix} \mathbf{H} & \mathbf{N} \\ \hline \mathbf{M} & \mathbf{L} \end{bmatrix} \begin{bmatrix} \Delta \boldsymbol{\theta} \\ \hline \frac{\Delta |\mathbf{V}|}{|\mathbf{V}|} \end{bmatrix}$$

Onde as submatrizes representam as derivadas parciais analíticas:
$$\mathbf{H} = \frac{\partial \mathbf{P}}{\partial \boldsymbol{\theta}}, \quad \mathbf{N} = |\mathbf{V}| \frac{\partial \mathbf{P}}{\partial |\mathbf{V}|}, \quad \mathbf{M} = \frac{\partial \mathbf{Q}}{\partial \boldsymbol{\theta}}, \quad \mathbf{L} = |\mathbf{V}| \frac{\partial \mathbf{Q}}{\partial |\mathbf{V}|}$$

As derivadas parciais possuem expressões explícitas elegantes:
- **Para termos fora da diagonal ($k \ne m$):**
  $$H_{km} = \frac{\partial P_k}{\partial \theta_m} = |V_k||V_m| \big( G_{km} \sin\theta_{km} - B_{km} \cos\theta_{km} \big)$$
  $$N_{km} = |V_m| \frac{\partial P_k}{\partial |V_m|} = |V_k||V_m| \big( G_{km} \cos\theta_{km} + B_{km} \sin\theta_{km} \big)$$
  $$M_{km} = \frac{\partial Q_k}{\partial \theta_m} = -N_{km}$$
  $$L_{km} = |V_m| \frac{\partial Q_k}{\partial |V_m|} = H_{km}$$

- **Para termos diagonais ($k = m$):**
  $$H_{kk} = -Q_k - B_{kk}|V_k|^2$$
  $$N_{kk} = P_k + G_{kk}|V_k|^2$$
  $$M_{kk} = P_k - G_{kk}|V_k|^2$$
  $$L_{kk} = Q_k - B_{kk}|V_k|^2$$

#### Algoritmo de Execução:
1. Definir estimativa inicial (*flat start*): $\theta_k^{(0)} = 0^\circ$ e $|V_k^{(0)}| = 1{,}0\text{ pu}$ para todas as barras $PQ$;
2. Calcular $P_k^{(\nu)}$ e $Q_k^{(\nu)}$ pelas expressões não lineares;
3. Avaliar os resíduos $\Delta P_k^{(\nu)}$ e $\Delta Q_k^{(\nu)}$. Se $\max(|\Delta P|, |\Delta Q|) < \varepsilon$ (tolerância, tipicamente $10^{-4}\text{ pu}$ ou $0{,}01\text{ MW}$), **convergência atingida**;
4. Avaliar a matriz Jacobiana $\mathbf{J}^{(\nu)}$;
5. Resolver o sistema linear por fatoração LU: $\mathbf{J}^{(\nu)} \Delta \mathbf{x}^{(\nu)} = \Delta \mathbf{S}^{(\nu)}$;
6. Atualizar as variáveis de estado:
   $$\theta_k^{(\nu+1)} = \theta_k^{(\nu)} + \Delta \theta_k^{(\nu)}$$
   $$|V_k^{(\nu+1)}| = |V_k^{(\nu)}| + \Delta |V_k^{(\nu)}|$$
7. Incrementar iteração ($\nu \leftarrow \nu + 1$) e retornar ao passo 2.

> [!tip] Desempenho do Newton-Raphson
> 
> Apresenta **convergência quadrática** na vizinhança da solução ($\|e^{(\nu+1)}\| \propto \|e^{(\nu)}\|^2$). O número de iterações para atingir a tolerância rigorosa é pequeno e **praticamente invariante com o tamanho do sistema** (comumente $3$ a $5$ iterações para redes de $10$ ou $5000$ barras).

---

### 5.3 Método Desacoplado e Desacoplado Rápido (FDLF)
Em linhas de transmissão de Extra Alta Tensão (EAT) e Alta Tensão (AT), as seguintes propriedades físicas se manifestam:
1. Relação $X/R \gg 1 \implies r_{km} \ll x_{km} \implies G_{km} \approx 0$;
2. As diferenças angulares entre barras vizinhas sob operação normal são muito pequenas: $\theta_{km} \approx 0 \implies \cos\theta_{km} \approx 1$ e $\sin\theta_{km} \approx \theta_{km}$;
3. A potência ativa $P$ é fortemente acoplada ao ângulo de fase $\theta$ e fracamente acoplada à magnitude da tensão $|V|$ ($\mathbf{N} \approx \mathbf{0}$);
4. A potência reativa $Q$ é fortemente acoplada à magnitude da tensão $|V|$ e fracamente acoplada ao ângulo $\theta$ ($\mathbf{M} \approx \mathbf{0}$).

Desprezando as submatrizes fora da diagonal principal ($\mathbf{N} = \mathbf{0}, \mathbf{M} = \mathbf{0}$), o sistema se desacopla em dois blocos independentes:
$$\Delta \mathbf{P} = \mathbf{H} \Delta \boldsymbol{\theta}$$
$$\Delta \mathbf{Q} = \mathbf{L} \frac{\Delta |\mathbf{V}|}{|\mathbf{V}|}$$

No **Método Desacoplado Rápido (FDLF - Stott & Alsac)**, simplifica-se ainda mais as matrizes $\mathbf{H}$ e $\mathbf{L}$, convertendo-as em matrizes puramente topológicas $\mathbf{B}'$ e $\mathbf{B}''$, que são **invariantes no tempo**:
$$\frac{\Delta \mathbf{P}}{|\mathbf{V}|} = \mathbf{B}' \Delta \boldsymbol{\theta}$$
$$\frac{\Delta \mathbf{Q}}{|\mathbf{V}|} = \mathbf{B}'' \Delta |\mathbf{V}|$$

As matrizes $\mathbf{B}'$ e $\mathbf{B}''$ são fatoradas por eliminação de Gauss **uma única vez** no início da execução. Em cada iteração realizam-se apenas substituições diretas e retroativas, tornando a execução até 5 vezes mais rápida que o Newton-Raphson tradicional.

---

### 5.4 Fluxo de Potência Linearizado (DC Power Flow)
Para estudos de planejamento da expansão de longo prazo, despacho econômico e análise massiva de contingências ($N-1$), adota-se o modelo linearizado (conhecido popularmente como **Fluxo DC**).

#### Hipóteses Simplificadoras:
1. **Resistência desprezível:** $r_{km} \approx 0 \implies g_{km} = 0$ e $b_{km} = -\frac{1}{x_{km}}$ (linhas sem perdas);
2. **Tensões nominais em todas as barras:** $|V_k| \approx 1{,}0\text{ pu}$;
3. **Pequenas aberturas angulares:** $\sin(\theta_k - \theta_m) \approx \theta_k - \theta_m$;
4. **Desconsideração do fluxo reativo:** estuda-se apenas o fluxo de potência ativa $P$.

A equação do fluxo de potência na linha entre a barra $k$ e $m$ simplifica-se para:
$$P_{km} = \frac{\theta_k - \theta_m}{x_{km}} = -B_{km}(\theta_k - \theta_m)$$

A injeção de potência ativa líquida na barra $k$ torna-se uma equação linear direta:
$$P_k = \sum_{m \in \Omega_k} \frac{\theta_k - \theta_m}{x_{km}}$$

Na forma matricial global:
$$\mathbf{P} = \mathbf{B}_{bus} \boldsymbol{\theta}$$

Eliminando a linha e a coluna correspondentes à barra de referência (onde $\theta_{slack} = 0$), a matriz $\mathbf{B}_{bus}^{red}$ torna-se simétrica, não singular e estritamente inversível:
$$\boldsymbol{\theta} = (\mathbf{B}_{bus}^{red})^{-1} \mathbf{P}^{red}$$

A solução do fluxo DC é obtida em **uma única resolução linear direta**, sem iterações numéricas.

---

## 6. Controle de Tensão e Compensação Reativa no SEP

A manutenção dos níveis de tensão dentro dos padrões regulamentares da ANEEL (PRODIST Módulo 8: faixa adequada de $0{,}95\text{ pu} \le |V| \le 1{,}05\text{ pu}$ em regime normal) é vital para a estabilidade de tensão e segurança operativa do SIN.mermaid
graph TD
    subgraph Equipamentos_Controle_Reativo ["Equipamentos de Suporte de Tensão e Reativos"]
        Cap["Bancos de Capacitores Shunt<br/>(Injeta Q proporcional a V²)"]
        Reat["Reatores Shunt<br/>(Absorve Q indutivo / Efeito Ferranti)"]
        Sync["Compensadores Síncronos<br/>(Controle contínuo: sub/sobre-excitação)"]
        LTC["Transformadores com Tap Sob Carga LTC<br/>(Ajusta relação de espiras a:1)"]
        FACTS["Sistemas FACTS (SVC, STATCOM)<br/>(Eletrônica de potência ultrarrápida)"]
    end
### 6.1 O Acoplamento Forte Q-V em Redes de Alta Tensão
Pelo modelo simplificado de dois nós conectados por reatância $X$:
$$|V_1| - |V_2| \approx \frac{R P_{12} + X Q_{12}}{|V_1|}$$
Como $X \gg R$, tem-se que a queda de magnitude de tensão ao longo do ramo é primariamente governada pelo **fluxo de potência reativa $Q_{12}$**:
$$\Delta |V| \approx \frac{X Q}{|V_1|}$$

### 6.2 Bancos de Capacitores e Reatores Shunt
- **Bancos de Capacitores em Paralelo (Shunt):**
  - Fornecem potência reativa capacitiva localmente ($Q_C = \omega C |V|^2$), liberando capacidade de transmissão nas LTs e elevando o perfil de tensão nas barras de carga durante horários de pico;
  - *Limitação Crítica:* Como $Q_C \propto |V|^2$, quando a tensão da barra sofre afundamento severo, a capacidade de injeção de reativos do capacitor cai quadraticamente, acelerando o colapso de tensão.
- **Bancos de Reatores em Paralelo (Shunt Reactors):**
  - Absorvem potência reativa indutiva excessiva ($Q_L = \frac{|V|^2}{\omega L}$);
  - São energizados nos períodos de carga leve (madrugada e finais de semana) para compensar a susceptância capacitiva shunt das longas linhas de transmissão de $500\text{ kV}$ e $800\text{ kV}$ e mitigar a sobretensão do **efeito Ferranti**.

### 6.3 Máquinas e Condensadores Síncronos
Os geradores síncronos e condensadores síncronos (máquinas síncronas operando sem carga mecânica no eixo) oferecem suporte contínuo de reativos mediante a variação da corrente de campo no enrolamento de excitação ($I_f$):
1. **Regime Sobre-excitado ($E_a > |V|$):** A máquina opera com [[Aula 02 - Análise de Potência em Regime CA - Potência Ativa, Reativa, Aparente e Fator de Potência|fator de potência]] indutivo visto de seus terminais, comportando-se fisicamente como um gerador de reativos (fornece $Q$ capacitivo ao sistema, elevando a tensão);
2. **Regime Subexcitado ($E_a < |V|$):** A máquina opera com corrente adiantada em relação à tensão terminal, comportando-se fisicamente como uma carga reativa indutiva (absorve o excedente de reativos capacitivos do sistema, reduzindo a tensão).

### 6.4 Transformadores com Comutação de Tap sob Carga (LTC/OLTC) e FACTS
- **LTC (Load Tap Changer):** Altera a relação de transformação $a = \frac{N_1}{N_2}$ em degraus de $\pm 10\%$ a $\pm 15\%$ sob plena carga para compensar flutuações sazonais de tensão na média e alta tensão;
- **SVC (Static Var Compensator):** Utiliza [[Tiristor Retificador Controlado de Silício (SCR)|tiristores]] (TCR e TSC) para prover injeção ou absorção de reativos com resposta na ordem de ciclos ($20 - 40\text{ ms}$);
- **[[Compensadores Estáticos de Reativos e Geradores Estáticos de Var (STATCOM e SVG)|STATCOM]] (Static Synchronous Compensator):** Baseado em inversores fonte de tensão (VSC) com transistores [[Transistor Bipolar de Gate Isolado (IGBT)|IGBT]], capaz de injetar corrente capacitiva nominal mesmo sob tensões de barra extremamente degradadas ($I_C$ constante), superior ao capacitor chaveado.

---

## 7. 📝 Estudo de Caso Numérico Resolvido Passo a Passo

### Formulação do Problema:
Considere um sistema elétrico de potência didático de **3 barras** interligado em anel:
- **Barra 1:** Barra de Folga (Slack), $V_1 = 1{,}0 \angle 0^\circ\text{ pu}$;
- **Barra 2:** Barra de Carga ($PQ$), demanda de potência ativa $P_2 = -5{,}0\text{ pu}$ e reativa $Q_2 = -2{,}0\text{ pu}$;
- **Barra 3:** Barra de Carga ($PQ$), demanda de potência ativa $P_3 = -10{,}0\text{ pu}$ e reativa $Q_3 = -4{,}0\text{ pu}$;
- Reatâncias série das linhas de transmissão (com resistências e capacitâncias shunt desprezíveis):
  - Linha 1-2: $x_{12} = 0{,}03333\text{ pu} \implies b_{12} = \frac{1}{0{,}03333} = 30\text{ pu}$
  - Linha 1-3: $x_{13} = 0{,}05000\text{ pu} \implies b_{13} = \frac{1}{0{,}05000} = 20\text{ pu}$
  - Linha 2-3: $x_{23} = 0{,}05000\text{ pu} \implies b_{23} = \frac{1}{0{,}05000} = 20\text{ pu}$mermaid
graph LR
    B1["Barra 1 (Slack)<br/>V1 = 1.0 ∠ 0°"] ---|x12 = 0.0333 pu| B2["Barra 2 (PQ)<br/>P2 = -5.0 pu"]
    B1 ---|x13 = 0.0500 pu| B3["Barra 3 (PQ)<br/>P3 = -10.0 pu"]
    B2 ---|x23 = 0.0500 pu| B3

Passo 1: Montagem da Matriz de Susceptância Nodal (${\mathbf{B}}_{bus}$)

Aplicando as regras de inspeção para a rede reativa:

• $B_{11}=b_{12}+b_{13}=30+20=50\text{ pu}$
• $B_{22}=b_{12}+b_{23}=30+20=50\text{ pu}$
• $B_{33}=b_{13}+b_{23}=20+20=40\text{ pu}$
• $B_{12}=B_{21}=-b_{12}=-30\text{ pu}$
• $B_{13}=B_{31}=-b_{13}=-20\text{ pu}$
• $B_{23}=B_{32}=-b_{23}=-20\text{ pu}$

A matriz ${\mathbf{B}}_{bus}$ completa é: ${\mathbf{B}}_{bus}=\left[\begin{array}{ccc}50& -30& -20\\ -30& 50& -20\\ -20& -20& 40\end{array}\right]$

Passo 2: Redução Matricial para a Barra de Referência

Como a Barra 1 é a referência angular (${\theta }_1=0$), elimina-se a linha 1 e a coluna 1 da formulação: ${\mathbf{B}}_{bus}^{red}=\left[\begin{array}{cc}{B}_{22}& B_{23}\\ B_{32}& B_{33}\end{array}\right]=\left[\begin{array}{cc}50& -20\\ -20& 40\end{array}\right]$

O vetor de injeções líquidas de potência ativa nas barras 2 e 3 é: ${\mathbf{P}}^{red}=\left[\begin{array}{c}{P}_2\\ P_3\end{array}\right]=\left[\begin{array}{c}-5,0\\ -10,0\end{array}\right]$

O sistema linear do Fluxo de Potência DC é dado por: $\left[\begin{array}{cc}50& -20\\ -20& 40\end{array}\right]\left[\begin{array}{c}{\theta }_2\\ {\theta }_3\end{array}\right]=\left[\begin{array}{c}-5,0\\ -10,0\end{array}\right]$

Passo 3: Solução Analítica dos Ângulos de Barra

Calculando o determinante da matriz ${\mathbf{B}}_{bus}^{red}$: $\det ({\mathbf{B}}_{bus}^{red})=(50)(40)-(-20)(-20)=2000-400=1600$

A matriz inversa é: $({\mathbf{B}}_{bus}^{red}{)}^{-1}=\frac{1}{1600}\left[\begin{array}{cc}40& 20\\ 20& 50\end{array}\right]$

Resolvendo para o vetor angular: ${\theta }_2=\frac{40(-5,0)+20(-10,0)}{1600}=\frac{-200-200}{1600}=\frac{-400}{1600}=-0,250\text{ rad}$ ${\theta }_3=\frac{20(-5,0)+50(-10,0)}{1600}=\frac{-100-500}{1600}=\frac{-600}{1600}=-0,375\text{ rad}$

Convertendo para graus sexagesimais: ${\theta }_2=-0,250\times \frac{180^{\circ }}{\pi }\approx -14,32^{\circ }$ ${\theta }_3=-0,375\times \frac{180^{\circ }}{\pi }\approx -21,49^{\circ }$

Passo 4: Cálculo dos Fluxos de Potência Ativa nos Ramos

$P_{12}=\frac{{\theta }_1-{\theta }_2}{x_{12}}=\frac{0-(-0,25)}{0,03333}=\frac{0,25}{0,03333}=+7,50\text{ pu}$ $P_{13}=\frac{{\theta }_1-{\theta }_3}{x_{13}}=\frac{0-(-0,375)}{0,05000}=\frac{0,375}{0,05000}=+7,50\text{ pu}$ $P_{23}=\frac{{\theta }_2-{\theta }_3}{x_{23}}=\frac{-0,250-(-0,375)}{0,05000}=\frac{0,125}{0,05000}=+2,50\text{ pu}$

Verificação dos Balanços Nodais:

• Barra 1 (Slack): $P_1=P_{12}+P_{13}=7,50+7,50=+15,00\text{ pu}$ (geração fornecida pela barra slack);
• Barra 2: $P_{12}-P_{23}=7,50-2,50=5,00\text{ pu}=-P_2$ (atende exatamente à carga da Barra 2);
• Barra 3: $P_{13}+P_{23}=7,50+2,50=10,00\text{ pu}=-P_3$ (atende exatamente à carga da Barra 3). Balanço fechado com consistência analítica total de $100\%$.

---

8. 📚 Questões Resolvidas das Materiais Teóricos de Referência

Esta seção reúne e resolve passo a passo todas as 8 questões conceituais e numéricas presentes na Material Teórico de Referência Fluxo de Potência.

---

Questão 01 — Variáveis e Tipos de Barra no SEP

Enunciado Oficial

Na formulação básica do problema de fluxo de potência é preciso conhecer quatro variáveis para cada barra do sistema: módulo da tensão, fase da tensão, potência ativa e potência reativa. A partir dessas variáveis define-se três tipos de barra. Sobre os tipos de barra, é correto afirmar que:

• (A) As barras de geração fecham o balanço de potência do sistema.
• (B) A barra PQ refere-se às barras que fornecem potência ativa e reativa ao sistema.
• (C) O sistema pode possuir quantas barras de referência forem necessárias.
• (D) A barra PV tem tensão constante.
• (E) Na barra de referência é conhecida a potência ativa e reativa fornecida.

Resolução & Gabarito Comentado

Alternativa Correta: (D)

Análise Técnica:

1. As barras PV (também denominadas barras de geração ou de tensão controlada) possuem como variáveis especificadas (conhecidas) a potência ativa líquida injetada $P_k$ e a magnitude da tensão $|V_k|$. Logo, o módulo da sua tensão permanece constante durante a resolução do algoritmo de fluxo, dentro dos limites físicos de suporte de reativos da máquina ($Q_{min}\le Q\le Q_{max}$).
2. Análise dos distratores:
   • (A) Falso: quem fecha o balanço de perdas e potência ativa do sistema é a barra de folga (Slack / Swing), não as barras PV comuns;
   • (B) Falso: as barras $PQ$ são barras de carga, que consomem potência ativa e reativa do sistema;
   • (C) Falso: um sistema elétrico síncrono necessita de exatamente uma barra de referência angular única ($\theta =0^{\circ }$) para fixar o referencial fasorial global;
   • (E) Falso: na barra de referência, as potências $P$ e $Q$ são incógnitas a serem determinadas após o cálculo dos fluxos de carga.

---

Questão 02 — Formulação Básica e Leis Físicas do Fluxo

Enunciado Oficial

A respeito da formulação básica do problema de fluxo de potência, é possível afirmar que:

• (A) A modelagem das perdas é feita na própria equação da linha de transmissão.
• (B) O somatório das injeções de potência é igual ao somatório de potências que fluem pelas linhas de transmissão e transformadores.
• (C) O estudo de fluxo de potência é válido para o sistema operando em qualquer condição.
• (D) As equações que representam o fluxo de potência no sistema são do tipo diferenciais.
• (E) A matriz de admitância representa apenas os elementos série do sistema.

Resolução & Gabarito Comentado

Alternativa Correta: (B)

Análise Técnica:

1. A formulação matemática do fluxo de potência ancora-se na aplicação rigorosa da Primeira Lei de Kirchhoff (Lei Nodal das Correntes e conservação de potência complexa em cada nó). O balanço de potência em cada barra estipula que a potência líquida injetada no nó ($S_k=S_{Gk}-S_{Ck}$) deve ser idêntica à soma de todas as potências complexas que partem daquele nó através dos ramos de interconexão e elementos em derivação incidentes: $P_k={\sum }_{m\in {\Omega }_k}{P}_{km},Q_k={\sum }_{m\in {\Omega }_k}{Q}_{km}$
2. Análise dos distratores:
   • (A) Falso: as perdas surgem da integração global dos fluxos de tensão e corrente nos resistores dos circuitos equivalentes;
   • (C) Falso: o fluxo de potência clássico modela estritamente o regime permanente senoidal; não é aplicável a regimes transitórios eletromecânicos ou eletromagnéticos rápidos (estudados por equações diferenciais no domínio do tempo);
   • (D) Falso: as equações de fluxo de potência em regime permanente são estritamente algébricas e não lineares;
   • (E) Falso: a matriz de admitância nodal ${\mathbf{Y}}_{bus}$ engloba tanto as impedâncias longitudinais série quanto as admitâncias capacitivas e reativas transversais (shunt).

---

Questão 03 — Modelagem Algébrica e Graus de Liberdade

Enunciado Oficial

Uma vez que o sistema elétrico esteja operando em condições normais, o problema de fluxo de potência é modelado matematicamente por equações algébricas. O equacionamento para solução do problema considera os tipos de barra e variáveis envolvidas. A respeito da modelagem das equações e suas soluções é correto afirmar que:

• (A) Para cada barra do sistema são previstas 4 equações e 4 variáveis, sendo o sistema possível e determinado.
• (B) As equações para solução do problema de fluxo de potência são lineares.
• (C) O tipo de barra não influencia na quantidade de variáveis a ser determinada.
• (D) Os métodos de solução aplicados ao fluxo de potência determinam a tensão e frequência em cada barra do sistema.
• (E) Das variáveis de estado nas barras, 2 são previamente conhecidas e 2 são calculadas por métodos diversos de solução de equações.

Resolução & Gabarito Comentado

Alternativa Correta: (E)

Análise Técnica:

1. Em cada nó $k$ do SEP, existem exatamente 4 grandezas físicas associadas: $|V_k|,{\theta }_k,P_k,Q_k$. Por outro lado, a teoria de circuitos fornece apenas 2 equações nodais fundamentais independentes por barra (a equação do balanço de potência ativa $P_k$ e a equação de potência reativa $Q_k$).
2. Logo, para que o sistema algébrico seja matematicamente determinado (número de equações igual ao número de incógnitas), é imperativo que exatamente 2 variáveis sejam pré-fixadas e conhecidas, permitindo calcular as 2 variáveis restantes através dos algoritmos de solução.

---

Questão 04 — Métodos de Solução e Linearização do SEP

Enunciado Oficial

A respeito dos métodos de solução aplicados ao problema de fluxo de potência, é correto afirmar que:

• (A) A linearização do modelo de fluxo de potência permite boas aproximações das variáveis de estado em sistemas de extra alta tensão.
• (B) A solução de um fluxo de potência somente é possível em sistemas computacionais.
• (C) Os métodos lineares mais utilizados são o de Gauss-Jacobi e Gauss-Seidel.
• (D) O fluxo de potência linear é um modelo simplificado que modela apenas o fluxo de potência reativa no sistema.
• (E) O método de Newton para solução de equações não lineares é pouco utilizado na solução de fluxo de potência em razão de sua difícil convergência.

Resolução & Gabarito Comentado

Alternativa Correta: (A)

Análise Técnica:

1. Em sistemas de Extra Alta Tensão (EAT, com tensões de $230\text{ kV}$, $500\text{ kV}$ e $800\text{ kV}$), as linhas de transmissão operam com condutores múltiplos por fase e espaçamento amplo, resultando em reatâncias indutivas dezenas de vezes maiores que as resistências de condutor ($X/R$ na faixa de $10$ a $30$).
2. Além disso, os reguladores de tensão mantêm os perfis de tensão muito próximos de $1,0\text{ pu}$, e as defasagens angulares entre barras vizinhas raramente ultrapassam $10^{\circ }$ ($\sin \theta \approx \theta$, $\cos \theta \approx 1$). Sob essas premissas físicas legítimas, o modelo linearizado (Fluxo DC) reproduz o comportamento da potência ativa com exatidão superior a $95\%$, sendo amplamente utilizado para estudos operativos e de planejamento.
3. Análise dos distratores:
   • (C) Falso: Gauss-Jacobi e Gauss-Seidel são métodos iterativos de ponto fixo para resolver equações não lineares;
   • (D) Falso: o modelo linearizado (DC) modela exclusivamente o fluxo de potência ativa, desprezando o fluxo de reativos;
   • (E) Falso: o método de Newton-Raphson é o mais utilizado mundialmente devido à sua robusta convergência quadrática.

---

Questão 05 — Solução Linear de Sistema de 3 Barras

Enunciado Oficial

Considere o sistema elétrico composto por 3 barras com as seguintes matrizes reduzidas de injeção e admitância (adotando a Barra 1 como referência angular ${\theta }_1=0$): $\left[\begin{array}{c}-5,0\\ -10,0\end{array}\right]=\left[\begin{array}{cc}50& -20\\ -20& 40\end{array}\right]\left[\begin{array}{c}{\theta }_2\\ {\theta }_3\end{array}\right]$ Com base nos dados fornecidos, os ângulos das barras 2 e 3 em radianos serão:

• (A) ${\theta }_2=-0,35\text{ rad}$ e ${\theta }_3=-0,275\text{ rad}$
• (B) ${\theta }_2=-0,55\text{ rad}$ e ${\theta }_3=-0,175\text{ rad}$
• (C) ${\theta }_2=-0,25\text{ rad}$ e ${\theta }_3=-0,375\text{ rad}$
• (D) ${\theta }_2=-0,75\text{ rad}$ e ${\theta }_3=-0,355\text{ rad}$
• (E) ${\theta }_2=-0,45\text{ rad}$ e ${\theta }_3=-0,975\text{ rad}$

Resolução & Gabarito Comentado

Alternativa Correta: (C)

Dedução Matemática Completa:

1. Escrevendo o sistema de equações lineares: $\{\begin{array}{l}50{\theta }_2-20{\theta }_3=-5,0\text{(Equac¸a˜o 1)}\\ -20{\theta }_2+40{\theta }_3=-10,0\text{(Equac¸a˜o 2)}\end{array}$
2. Da Equação 2, dividindo todos os termos por $20$: $-{\theta }_2+2{\theta }_3=-0,5⟹{\theta }_2=2{\theta }_3+0,5$
3. Substituindo na Equação 1: $50(2{\theta }_3+0,5)-20{\theta }_3=-5,0$ $100{\theta }_3+25-20{\theta }_3=-5,0$ $80{\theta }_3=-30,0⟹{\theta }_3=-\frac{30}{80}=-0,375\text{ rad}$
4. Substituindo o valor de ${\theta }_3$ na expressão de ${\theta }_2$: ${\theta }_2=2(-0,375)+0,5=-0,750+0,500=-0,250\text{ rad}$
5. Portanto: $\boxed{{\theta }_2=-0,25\text{ rad},{\theta }_3=-0,375\text{ rad}}$

---

Questão 06 — Modelagem Não Linear de Newton-Raphson

Enunciado Oficial

O método de Newton será utilizado para a solução do ponto de operação de um sistema elétrico que contém 3 barras, interligadas por linhas de transmissão entre si. Com base na modelagem desse sistema, é correto afirmar que:

• (A) As equações do problema são do tipo lineares.
• (B) A matriz de admitância tem ordem 4x4.
• (C) A matriz jacobiana independe do tipo de barra no sistema.
• (D) As variáveis a calcular são o módulo e a fase da tensão nas barras.
• (E) O método de Newton não é eficiente para solução de fluxo de potência, pois necessita de muitas iterações para convergência.

Resolução & Gabarito Comentado

Alternativa Correta: (D)

Análise Técnica:

1. No fluxo de potência formulado pelo método de Newton-Raphson, o estado elétrico da rede é governado pelas equações não lineares de balanço de potência ativa e reativa.
2. As variáveis de estado fundamentais a determinar iterativamente a partir do resíduo das equações são o módulo $|V_k|$ e o ângulo de fase ${\theta }_k$ da tensão fasorial nos nós do sistema.
3. Análise dos distratores:
   • (A) Falso: as equações do fluxo de potência CA são inerentemente não lineares;
   • (B) Falso: para um sistema de 3 barras, a matriz de admitância ${\mathbf{Y}}_{bus}$ possui rigorosamente dimensão $3\times 3$;
   • (C) Falso: a estrutura e dimensão da matriz Jacobiana dependem explicitamente do tipo de cada barra (uma barra $PV$ elimina uma linha e coluna correspondente à correção de $|V|$);
   • (E) Falso: o método de Newton é extremamente eficiente por convergir em apenas $3$ a $5$ iterações.

---

Questão 07 — Fatores Influenciadores do Perfil de Tensão

Enunciado Oficial

Os sistemas elétricos de potência operam em condições de estabilidade quando o perfil de tensão nas barras está dentro de limites preestabelecidos. No entanto, diversos fatores podem modificar o ponto de operação de um sistema a ponto de ocorrerem mudanças bruscas nas tensões. Podemos citar como fator altamente relacionado à tensão de um sistema:

• (A) A variação de frequência.
• (B) A geração de energia ativa.
• (C) A variação súbita de carga.
• (D) A presença de componentes harmônicas de corrente.
• (E) O consumo de reativos por equipamentos indutivos.

Resolução & Gabarito Comentado

Alternativa Correta: (E)

Análise Técnica:

1. Conforme a teoria de transmissão em linhas altamente reativas ($X/R\gg 1$), a queda de tensão longitudinal ao longo dos ramos é ditada pela transferência de potência reativa: $\Delta |V|\approx \frac{RP+XQ}{|V|}\approx \frac{XQ}{|V|}$
2. Equipamentos fortemente indutivos — como motores de indução trifásicos operando com carga mecânica, transformadores operando a vazio ou com alta corrente de magnetização e reatores de linha — drenam elevados montantes de potência reativa do sistema. Essa circulação massiva de corrente indutiva reativa pelas reatâncias de dispersão das linhas e transformadores provoca quedas substanciais nos módulos de tensão nas barras adjacentes.

---

Questão 08 — Técnicas de Compensação Reativa e Máquinas Síncronas

Enunciado Oficial

As técnicas de compensação reativa são importantes não apenas para se evitar custos adicionais de energia excedente, mas também para garantir um melhor funcionamento dos equipamentos em uma instalação elétrica. A respeito das principais técnicas de compensação reativa é correto afirmar que:

• (A) Motores síncronos operando subexcitados podem ser utilizados para a correção de fator de potência em instalações predominantemente capacitivas.
• (B) Bancos de capacitores são pouco utilizados em função de seu custo de aquisição.
• (C) O aumento de potência ativa não contribui de forma eficiente no controle de fator de potência.
• (D) Máquinas síncronas necessitam de uma carga conectada ao seu eixo para produzir reativos para o sistema.
• (E) A instalação de bancos de capacitores é uma técnica de compensação de energia reativa para instalações predominantemente capacitivas.

Resolução & Gabarito Comentado

Alternativa Correta: (A)

Análise Técnica:

1. Um motor síncrono (ou condensador síncrono) permite controle contínuo e bidirecional de potência reativa através do ajuste da sua corrente de excitação de campo $I_f$:
   • Ao operar sobre-excitado, ele injeta potência reativa capacitiva no sistema (comportando-se como banco de capacitores);
   • Ao operar subexcitado ($E_a<|V|$), o fasor de corrente terminal fica atrasado em relação à tensão, fazendo com que a máquina absorva potência reativa indutiva da rede.
2. Portanto, em instalações industriais ou subestações que apresentem características predominantemente capacitivas (por exemplo, alimentadas por longas linhas de cabos subterrâneos ou com bancos de capacitores superdimensionados durante o período noturno de baixa carga), o motor síncrono subexcitado consome o excedente capacitivo, restabelecendo o fator de potência para níveis próximos da unidade e impedindo sobretensões inaceitáveis.

---

9. Resumo Conceitual da Aula e Conexões Curriculares

graph TD
    LF["Fluxo de Potência (Regime Permanente)"] --> Modelagem["Formulação Algébrica Não Linear (Kirchhoff + Injeções Nodais)"]
    Modelagem --> Ybus["Matriz de Admitância Nodal Ybus (G + jB)"]
    Modelagem --> Barras["Classificação Nodal: PQ, PV e Slack"]
    
    Ybus --> Metodos["Métodos de Solução Numérica"]
    Metodos --> GS["Gauss-Seidel (Linear, Lento)"]
    Metodos --> NR["Newton-Raphson (Quadrático, Jacobiana H, N, M, L)"]
    Metodos --> FDLF["Desacoplado Rápido (B', B'' Constantes)"]
    Metodos --> DC["Fluxo DC Linearizado (P = B · θ)"]
    
    Barras --> QV["Controle de Tensão e Reativos"]
    QV --> Comp["Compensação Shunt (Capacitores / Reatores), LTC e Síncronos"]

Conexão com o Módulo 04: Curto-Circuito

A matriz de admitância nodal ${\mathbf{Y}}_{bus}$ desenvolvida nesta aula para o cálculo de fluxo de carga em regime normal constitui a base analítica direta para a montagem da matriz de impedância nodal ${\mathbf{Z}}_{bus}={\mathbf{Y}}_{bus}^{-1}$, utilizada no cálculo de correntes de falta e dimensionamento de disjuntores na Aula 04.