☀️ Aula 02 - Energia Solar Fotovoltaica, Física dos Semicondutores, MPPT e Inversores

Navegação Estrutural

📍 MOC: 06 Sistemas de Energia Alternativa e Renováveis | ⬅️ Aula Anterior: Aula 01 - Transição Energética e Recursos Renováveis | ➡️ Próxima Aula: Aula 03 - Marco Legal da GD e Lei 14.300

Objetivos Pedagógicos e Escopo Técnico

• Compreender a física dos semicondutores e a mecânica quântica do efeito fotovoltaico na junção p-n de Silício cristalino com banda proibida (bandgap $E_g=1,12\text{ eV}$).
• Modelar matematicamente a célula solar pelo circuito equivalente de 1 diodo e 5 parâmetros através da equação transcendental de Shockley.
• Analisar o comportamento das curvas características $I\times V$ e $P\times V$, determinando o Ponto de Máxima Potência ($P_{\text{MPP}}$), Fator de Preenchimento ($FF$) e a dependência térmica dos parâmetros.
• Deduzir e comparar os algoritmos clássicos de rastreamento do ponto de máxima potência: Perturbação e Observação (P&O) e Condutância Incremental (IncCond).
• Analisar a eletrônica de potência de inversores fotovoltaicos (estágio elevador CC-CC e ponte inversora CC-CA com modulação SPWM) e dimensionar arranjos fotovoltaicos (strings) compatibilizados com as janelas operativas do inversor.

---

1. Física dos Semicondutores e o Efeito Fotovoltaico

A conversão direta da luz solar em energia elétrica ocorre no interior de materiais semicondutores dopados, governada pela interação da radiação eletromagnética com a estrutura de bandas de energia dos átomos.

1.1 A Junção P-N de Silício e a Banda Proibida (Bandgap)

O Silício cristalino (número atômico $Z=14$) é um elemento tetravalente que forma ligações covalentes estáveis em rede cristalina cúbica tipo diamante.

• Dopagem Tipo P (Aceitadora): A introdução controlada de átomos trivalentes (como o Boro) cria uma deficiência de elétrons na rede, gerando uma densidade predominante de portadores de carga positiva livres denominados lacunas.
• Dopagem Tipo N (Do Inventors): A introdução de átomos pentavalentes (como o Fósforo) adiciona elétrons de valência excedentes que passam a orbitar livremente na banda de condução como portadores majoritários negativos.

Quando as regiões P e N são unidas metalurgicamente:

1. Ocorre a difusão espontânea de elétrons da região N para a região P e de lacunas da região P para a N.
2. Na interface da junção, forma-se a Zona de Depleção (ou Zona de Carga Espacial), desprovida de portadores livres e povoada por íons fixos ionizados (cátions de fósforo no lado N e ânions de boro no lado P).
3. Esses íons estabelecem um intenso campo elétrico interno dirigido da região N para a região P, criando uma barreira de potencial de contato de cerca de $V_{bi}\approx 0,6\text{ a }0,7\text{ V}$.

1.2 Mecanismo de Geração e Separação de Cargas

A luz solar é composta por um fluxo contínuo de fótons de diferentes comprimentos de onda ($\lambda$), cada qual transportando uma quantidade discreta de energia dada pela relação quântica de Planck-Einstein:

$E_{\text{foˊton}}=h\cdot \nu =\frac{h\cdot c}{\lambda }$

Onde $h=6,626\times 10^{-34}\text{ J}\cdot \text{s}$ e $c=3\times 10^8\text{ m/s}$.

Para o Silício cristalino à temperatura ambiente de $300\text{ K}$ ($27^{\circ }\text{C}$), a energia de banda proibida (bandgap) é:

$\boxed{E_g=1,12\text{ eV}\approx 1,794\times 10^{-19}\text{ J}}$

O comprimento de onda de corte óptico absorvível pelo silício é:

${\lambda }_{\text{corte}}=\frac{h\cdot c}{E_g}=\frac{(6,626\times 10^{-34})\cdot (3\times 10^8)}{1,794\times 10^{-19}}\approx 1,108\times 10^{-6}\text{ m}=1.108\text{ nm}(\text{Faixa do infravermelho proˊximo})$

• Caso 1 ($E_{\text{foˊton}}<E_g$): O fóton atravessa a estrutura semicondutora sem ser absorvido (o silício é transparente a fótons com $\lambda >1.108\text{ nm}$).
• Caso 2 ($E_{\text{foˊton}}\ge E_g$): O fóton é absorvido, transfere integralmente sua energia a um elétron da banda de valência, que salta para a banda de condução, criando um par elétron-lacuna fotogerado. A energia excedente ($E_{\text{foˊton}}-E_g$) é dissipada instantaneamente na forma de vibrações na rede cristalina (aquecimento térmico por termalização de portadores).
• Separação pelo Campo Interno: Os elétrons livres criados na vizinhança da zona de depleção são acelerados pelo campo elétrico interno em direção à região N, enquanto as lacunas são empurradas para a região P. Conectando-se um circuito elétrico externo com carga, estabelece-se uma corrente contínua e a entrega de trabalho elétrico útil.

---

2. Modelagem Matemática da Célula Solar: O Circuito Equivalente de 1 Diodo

Para fins de projeto elétrico e simulação computacional, uma célula solar fotovoltaica real é modelada através do circuito equivalente de 1 diodo e 5 parâmetros:

2.1 A Equação Transcendental de Shockley

Aplicando a Lei de Kirchhoff das Correntes no nó superior:

$I=I_{\text{ph}}-I_d-I_{\text{sh}}$

A corrente do diodo ideal que emula a recombinação dos portadores na junção segue a equação de Shockley:

$I_d=I_0\cdot \left[\exp \left(\frac{V+I\cdot R_s}{n\cdot V_t}\right)-1\right]$

A corrente de fuga através da resistência em paralelo (shunt) $R_{\text{sh}}$ é:

$I_{\text{sh}}=\frac{V+I\cdot R_s}{R_{\text{sh}}}$

Substituindo na expressão principal, obtém-se a equação transcendental mestra da célula fotovoltaica:

$\boxed{I=I_{\text{ph}}-I_0\cdot \left[\exp \left(\frac{V+I\cdot R_s}{n\cdot V_t}\right)-1\right]-\frac{V+I\cdot R_s}{R_{\text{sh}}}}$

Onde:

• $I_{\text{ph}}$: Fotocorrente gerada pela luz, estritamente proporcional à irradiância solar incidente $G$ (${\text{W/m}}^2$): $I_{\text{ph}}=\left[I_{\text{sc,STC}}+{\alpha }_{\text{Isc}}\cdot (T-T_{\text{STC}})\right]\cdot \frac{G}{G_{\text{STC}}}$
• $I_0$: Corrente de saturação reversa de difusão do diodo (fortemente governada pela temperatura absoluta): $I_0(T)=I_{0,\text{STC}}\cdot {\left(\frac{T}{T_{\text{STC}}}\right)}^3\cdot \exp \left[\frac{q\cdot E_g}{n\cdot k}\cdot \left(\frac{1}{T_{\text{STC}}}-\frac{1}{T}\right)\right]$
• $V_t$: Tensão térmica gerada pela agitação térmica eletrônica: $V_t=\frac{k\cdot T}{q}$ Onde $k=1,3806\times 10^{-23}\text{ J/K}$ (Constante de Boltzmann), $q=1,6022\times 10^{-19}\text{ C}$ (Carga elementar do elétron) e $T$ é a temperatura em Kelvin. A $T=25^{\circ }\text{C}$ ($298,15\text{ K}$), $V_t\approx 25,69\text{ mV}$.
• $n$: Fator de idealidade do diodo (tipicamente entre $1,0$ e $1,3$ para silício monocristalino).
• $R_s$: Resistência série interna (decorrente da resistividade intrínseca do silício, da resistência de contato dos dedos metálicos de prata na face frontal e barramentos). Idealmente $R_s\to 0$.
• $R_{\text{sh}}$: Resistência paralelo ou shunt (decorrente de caminhos de fuga microscópicos e defeitos cristalinos nas bordas da célula). Idealmente $R_{\text{sh}}\to \infty$.

---

3. Curvas Características $I\times V$ e $P\times V$

A avaliação do desempenho de células e módulos solares é padronizada mundialmente nas Condições Padrão de Ensaio (STC - Standard Test Conditions):

• Irradiância Solar: $G=1.000{\text{ W/m}}^2$;
• Temperatura da Célula: $T=25^{\circ }\text{C}$;
• Massa de Ar Espectral: $AM=1,5$ (distribuição espectral da luz solar atravessando a atmosfera com ângulo de elevação de $41,8^{\circ }$).

Figura 2.2: Curvas características I-V e P-V do módulo fotovoltaico, ponto de máxima potência (MPP) e dependência com irradiância e temperatura.

🕹️ Simulador Dinâmico Interativo (Laboratório Virtual de Energia Solar e MPPT)

Ajuste em tempo real a Irradiância Solar ($G$), a Temperatura da Célula ($T$) e a Tensão do Inversor ($V$), visualizando a transição das curvas $Iext-V$ e $Pext-V$ e experimentando o algoritmo de rastreamento do ponto de máxima potência (MPPT Perturb & Observe):

👉 ▶️ Clique para Abrir o Simulador Solar em Tela Cheia no Navegador (Arquivo no cofre: [[simulador_solar_pv_iv_mppt.html]])

3.1 Parâmetros Fundamentais de Catálogo

1. Corrente de Curto-Circuito ($I_{\text{sc}}$): Corrente obtida curto-circuitando os terminais de saída da célula ($V=0$). Como $R_s\ll R_{\text{sh}}$, tem-se rigorosamente $I_{\text{sc}}\approx I_{\text{ph}}$.
2. Tensão de Circuito Aberto ($V_{\text{oc}}$): Tensão medida com os terminais em aberto ($I=0$). Desprezando $R_{\text{sh}}$: $0\approx I_{\text{ph}}-I_0\cdot \exp \left(\frac{V_{\text{oc}}}{n{V}_t}\right)⟹\boxed{V_{\text{oc}}=n\cdot V_t\cdot \ln \left(\frac{I_{\text{ph}}}{I_0}+1\right)}$ Para uma célula de silício individual, $V_{\text{oc}}\approx 0,60\text{ a }0,72\text{ V}$. Um módulo com 72 células em série atinge $V_{\text{oc}}\approx 45\text{ a }50\text{ V}$.
3. Ponto de Máxima Potência ($P_{\text{MPP}}$): Ocorre nas coordenadas $(V_{\text{MPP}},I_{\text{MPP}})$ onde o produto $V\cdot I$ atinge seu ápice absoluto: $P_{\text{MPP}}=V_{\text{MPP}}\cdot I_{\text{MPP}}[\text{W ou Wp}]$
4. Fator de Preenchimento (Fill Factor - FF): Mede a “quadratura” da curva $I\times V$, representando a relação entre a potência máxima real entregue e o produto da corrente de curto pela tensão de circuito aberto: $\boxed{FF=\frac{V_{\text{MPP}}\cdot I_{\text{MPP}}}{V_{\text{oc}}\cdot I_{\text{sc}}}=\frac{P_{\text{MPP}}}{V_{\text{oc}}\cdot I_{\text{sc}}}}$ Em módulos modernos de alta tecnologia (PERC, TOPCon e Heterojunção HJT), $FF$ situa-se entre $0,78$ e $0,83$.
5. Eficiência de Conversão Global ($\eta$): Razão entre a potência máxima de pico produzida e a potência radiante incidente sobre a área física da célula $A$ (${\text{m}}^2$): $\eta =\frac{P_{\text{MPP}}}{G\cdot A}=\frac{V_{\text{oc}}\cdot I_{\text{sc}}\cdot FF}{G\cdot A}$

3.2 O Impacto Severo da Temperatura Operacional

Enquanto a elevação da irradiância solar eleva a corrente de forma linear quase pura ($I\propto G$), a elevação da temperatura ambiente e da célula degrada fortemente a tensão:

• O aumento de temperatura estreita a banda proibida $E_g$ e eleva exponencialmente a corrente de recombinação reversa $I_0$.
• O coeficiente térmico de tensão ${\beta }_{V_{\text{oc}}}$ é negativo, situando-se em torno de $-0,27\%{/}^{\circ }\text{C}$ a $-0,30\%{/}^{\circ }\text{C}$.
• O coeficiente térmico da potência máxima ${\gamma }_{P_{\text{MPP}}}$ situa-se tipicamente entre $-0,30\%{/}^{\circ }\text{C}$ e $-0,38\%{/}^{\circ }\text{C}$.

O Efeito Térmico no Brasil

Em dias ensolarados de verão, módulos fotovoltaicos instalados em telhados podem atingir temperaturas de operação de $60^{\circ }\text{C}$ a $70^{\circ }\text{C}$ ($\Delta T=+40^{\circ }\text{C}$ acima de STC). Isso provoca uma perda líquida de potência de até: $\Delta P=40^{\circ }\text{C}\times (-0,35\%{/}^{\circ }\text{C})=-14\%\text{ de perda contıˊnua de gerac¸a˜o!}$

---

4. Algoritmos de Rastreamento do Ponto de Máxima Potência (MPPT)

Como a irradiação solar e a temperatura variam continuamente a cada segundo com a passagem de nuvens e o vento, a curva $P\times V$ desloca-se dinamicamente. O inversor deve possuir um algoritmo de MPPT (Maximum Power Point Tracking) em seu estágio conversor CC-CC para forçar o gerador a operar ininterruptamente no pico da curva.

4.1 Algoritmo de Perturbação e Observação (P&O)

É o método heurístico mais popular devido à sua extrema simplicidade de programação em microcontroladores:

• O algoritmo mede a tensão $V(k)$ e a corrente $I(k)$ no instante amostral $k$, calculando a potência $P(k)=V(k)\cdot I(k)$.
• Compara com a potência calculada no passo anterior: $\Delta P=P(k)-P(k-1)$ e $\Delta V=V(k)-V(k-1)$.
• Lógica de Decisão:
  1. Se $\Delta P>0$ e $\Delta V>0$: O ponto está subindo a encosta esquerda $\to$ incrementa a tensão ($V_{\text{ref}}=V_{\text{ref}}+\delta V$).
  2. Se $\Delta P>0$ e $\Delta V<0$: O ponto está subindo a encosta direita $\to$ decrementa a tensão ($V_{\text{ref}}=V_{\text{ref}}-\delta V$).
  3. Se $\Delta P<0$ e $\Delta V>0$: O ponto ultrapassou o pico $\to$ inverte a busca e decrementa a tensão.
  4. Se $\Delta P<0$ e $\Delta V<0$: O ponto caiu à esquerda $\to$ inverte a busca e incrementa a tensão.
• Limitações: Em regime permanente, o P&O oscila continuamente em torno do ponto de máxima potência (gerando pequenas perdas estáticas). Além disso, sob variações bruscas e ultravelozes de irradiância (nuvens cúmulos com ventos fortes), o algoritmo pode tomar decisões na direção errada.

---

4.2 Algoritmo da Condutância Incremental (IncCond)

O algoritmo da condutância incremental baseia-se na condição matemática exata de derivada nula da curva de potência no ponto de extremo:

$\frac{dP}{dV}=\frac{d(V\cdot I)}{dV}=I+V\cdot \frac{dI}{dV}$

Dividindo ambos os membros pela tensão $V$ ($V\ne 0$):

$\frac{1}{V}\cdot \frac{dP}{dV}=\frac{I}{V}+\frac{dI}{dV}$

Analisando os três estados possíveis da derivada:

undefined