Fotovoltaïsche Zonnecel

Inleiding¶

De werking van een fotovoltaïsche zonnecel is gebaseerd op de interactie tussen fotonen (met een bepaalde energie $E_p$ ) en de vrije elektronen in een halfgeleidermateriaal. De elektronen kunnen door interactie met een invallend foton naar een hoger energieniveau worden gebracht. Als het elektron vervolgens terugvalt naar zijn oorspronkelijke energieniveau, kan de vrijgekomen energie worden gebruikt om extern elektrische stroom te genereren. In een beschouwing van een goede keuze voor de bandgap-energie (die afhangt van het gekozen materiaal), moeten we twee ‘tegenstrijdige’ situaties beschouwen:

$E_p > E_g$
Wanneer een invallend foton een energie heeft die groter is dan de zogenaamde bandgap-energie $E_g$ van het halfgeleidermateriaal ( $E_p > E_g$ ), dan zal de energie van een elektron in het halfgeleidermateriaal toenemen met een hoeveelheid energie $\Delta E=E_g$ . Het “overschot” aan foton-energie ( $E_p-E_g$ ) wordt omgezet in warmte, en gaat “verloren”. Het deel van de foton-energie dat in nuttige energie wordt omgezet neemt dus toe als de bandgap-energie groter wordt.
$E_p < E_g$
Wanneer het invallend foton een energie heeft die kleiner is dan de bandgap-energie $E_g$ van het halfgeleidermateriaal ( $E_p < E_g$ ), dan zal de energie van een elektron niet toenemen. De foton-energie $E_p$ wordt dan omgezet in warmte, en gaat “verloren”. Het aantal fotonen dat leidt tot nuttige energie neemt dus af als de bandgap-energie groter wordt.

Het resultaat van bovenstaande twee effecten is dat er een optimale bandgap-energie is waarvoor de zonnecel een maximale hoeveelheid foton-energie omzet in nuttige energie. In deze opgave gaan we uitwerken wat de optimale bandgap-energie is van een fotovoltaïsche zonnecel, voor invallende straling die overeenkomt met het fotonen-spectrum zoals uitgestraald door de zon. We beschouwen de zon als een zwarte straler met temperatuur $T_z=5800\,\mathrm{K}$ . Daarbij, de bandgap-energie van silicium zonnecellen is $1.10 \; \mathrm{eV} = 1.10 \times 1.602 \times 10^{-19} \; \mathrm{J} = 1.76 \times 10^{-19} \; \mathrm{J}$ .

Stralingswet van Planck – analytisch en numeriek¶

Volgens de stralingswet van Planck wordt de energie (in $\mathrm{W/m^2}$ ) uitgestraald door een zwart oppervlak met temperatuur $T$ binnen een golflengte-interval tussen $\lambda$ en $\lambda+\Delta \lambda$ gegeven door de Planck-kromme:

E_{\lambda,T} = \frac{2 \pi h c^2}{\lambda^5} \cdot \frac{1}{e^{\frac{hc}{\lambda k T}} - 1} \, \Delta \lambda

(1)

met:

$h$ = constante van Planck = $6.626 \cdot 10^{-34} \; \mathrm{J s}$
$c$ = lichtsnelheid = $2.998 \cdot 10^8 \; \mathrm{m/s}$
$\lambda$ = golflengte in $\mathrm{m}$
$T$ = absolute temperatuur in $\mathrm{K}$
$k$ = constante van Boltzmann = $1.381 \cdot 10^{-23} \; \mathrm{J/K}$

Onderstaande figuur laat zien hoe de verdeling van de stralingsenergie over de golflengtes afhangt van de temperatuur van de straler. We zien dat de zon (in goede benadering een zwarte straler met een temperatuur van $6000 \; \mathrm{K}$ ) de meeste stralingsenergie uitzendt bij zichtbaar licht golflengtes ( $400-800 \; \mathrm{nm}$ ).

Figure 1:De intensiteit van de uitgezonden straling als functie van de golflengte voor verschillende temperaturen van zwarte stralers. Figuur van wiki

In deze opgave gaan we analyseren welke implicaties de Planck kromme van de zon heeft voor het rendement van een fotovoltaïsche zonnecel.

Tijdens de analyse colleges hebben we afgeleid dat de totale energie die door een zwart oppervlak met temperatuur $T$ wordt gegeven door

E_T = \int_{0}^{\infty} E_{\lambda,T}\, d\lambda = 2\pi h c^{2} \int_{0}^{\infty} \frac{1}{\lambda^{5}\left( e^{\frac{hc}{\lambda k T}} - 1 \right)} \, d\lambda = \frac{2\pi^{5} k^{4}}{15 h^{3} c^{2}} T^{4} = \sigma T^{4}

(2)

met $\sigma = \text{constante van Stefan–Boltzmann} = 5.670 \cdot 10^{-8}\ \mathrm{W/(m^2K^4)}$

Integralen kunnen we benaderen met een zogenaamde Riemann som. Dit is geïllustreerd in onderstaande figuur (voor een andere functie dan de Planck-kromme).

Figure 2:De integraal, of oppervlak onder een grafiek, kan benaderd worden met de Riemann som. Figuur van wiki

E_T = 2\pi h c^2 \int_{0}^{\infty} \frac{1}{\lambda^{5}\left( e^{\frac{hc}{\lambda k T}} - 1 \right)} \, d\lambda \approx 2\pi h c^2 \lim_{\Delta \lambda \to 0} \sum_{i=1}^{\infty} \frac{1}{\lambda_i^{5}\left( e^{\frac{hc}{\lambda_i k T}} - 1 \right)} \, \Delta \lambda

(3)

met $\lambda_i = i \cdot \Delta \lambda$

Deze Riemann-som kunnen we oplossen met een Python-script. Uiteraard kunnen we niet een oneindig aantal stapjes berekenen in een computer. Daarom benaderen we in het Python-programma de Riemann-som voor een eindig aantal stapjes. Dat kan omdat de functie exponentieel afneemt voor grote golflengtes:

E_T \approx 2\pi h c^2 \lim_{\Delta \lambda_0 \to 0} \sum_{i=1}^{\infty} \frac{1}{\lambda_i^{5}\left( e^{\frac{hc}{\lambda_i k T}} - 1 \right)} \, \Delta \lambda_0 \approx 2\pi h c^2 \sum_{i=1}^{N} \frac{1}{\lambda_i^{5}\left( e^{\frac{hc}{\lambda_i k T}} - 1 \right)} \, \Delta \lambda_0

(4)

waarbij we ‘slechts’ $N$ stappen nemen met $\lambda_i = i \cdot \Delta \lambda_0$ . {note} Links, midden, rechts De Riemann-som kun je links, midden (centered) en rechts nemen. Voor de Planck kromme kiezen we deze rechts (door te beginnen bij $i=1$ ). Zouden we $i=0$ kiezen, dan krijgen we $\frac{1}{0}=\infty$ ... dat is niet gewenst.

{exercise} Riemannsom :label: ex_Riemann Schrijf een Python programma dat de benaderde Riemann som, vergelijking (4), uitrekent voor

$\Delta \lambda_0= 10 \; \mathrm{nm}$
$N=300$
$T=5800 \; \mathrm{K}$

Herhaal, om te controleren of we $N$ groot genoeg hebben gekozen en of we $\Delta \lambda_0$ klein genoeg hebben gekozen, je berekening voor

$\Delta \lambda_0=3 \; \mathrm{nm}$
$N=...$
$T=5800 \; \mathrm{K}$

Bedenk welke waarde van N je moet nemen voor een eerlijk vergelijk.

Vergelijk beide resultaten met de verwachte uitkomst gegeven als vergelijking (2).

import numpy as np
import matplotlib.pyplot as plt

#constanten:
h = 6.626e-34 # Planck's constant (J·s)
c = 2.998e8 # speed of light (m/s)
k = 1.381e-23 # Boltzmann constant (J/K)
T = 5800
a = h * c / (k * T) # vereenvoudiging van de som
o = 5.670*10**-8 # constante van boltzmann

# Functie

def E_T(l0, N):
    return (2*np.pi*h*c**2)*np.sum((1)/((l0*i)**5*(np.exp(a/(l0*i))-1))*l0 for i in range(1,N+1))


E = E_T(10*10**-9,300)
E2 = E_T(3*10**-9,1000)
check = o * T**4
print('De totale energie bij een delta lambda van 10 nm en een N van 300 is',E,'J')
print('De totale energie bij een delta lambda van 3 nm en een N van 1000 is',E2,'J')
print('volgens de eerste formule moet uit de gegevens komen dat deze een totale energie heeft van',check,'J')

De totale energie bij een delta lambda van 10 nm en een N van 300 is 62890429.07000511 J
De totale energie bij een delta lambda van 3 nm en een N van 1000 is 62886247.89776597 J
volgens de eerste formule moet uit de gegevens komen dat deze een totale energie heeft van 64164532.32 J

C:\Users\korne\AppData\Local\Temp\ipykernel_14412\3013309854.py:15: DeprecationWarning: Calling np.sum(generator) is deprecated, and in the future will give a different result. Use np.sum(np.fromiter(generator)) or the python sum builtin instead.
  return (2*np.pi*h*c**2)*np.sum((1)/((l0*i)**5*(np.exp(a/(l0*i))-1))*l0 for i in range(1,N+1))
C:\Users\korne\AppData\Local\Temp\ipykernel_14412\3013309854.py:15: RuntimeWarning: overflow encountered in exp
  return (2*np.pi*h*c**2)*np.sum((1)/((l0*i)**5*(np.exp(a/(l0*i))-1))*l0 for i in range(1,N+1))

De waarden van de riemannsom komen overeen met elkaar, maar gaan beter kloppen met wat het hoort te zijn als de N veel hoger wordt genomen.

Bepaling optimale bandgap-energie van een fotovoltaïsche zonnecel¶

Een foton met golflengte $\lambda$ heeft een energie $E_p=\frac{hc}{\lambda}$ . Met vergelijking (1) vinden we nu dat het aantal fotonen dat wordt uitgestraald door een zwart oppervlak met temperatuur $T$ binnen een golflengte-interval tussen $\lambda$ en $\lambda+\Delta \lambda$ wordt gegeven door

N_{\lambda, T} = \frac{E_{\lambda,T}}{E_p} = \frac{2\pi h c^{2}}{\lambda^{5}} \cdot \frac{1}{e^{\frac{hc}{\lambda k T}} - 1} \cdot \frac{\lambda}{hc}\, \Delta \lambda = \frac{2\pi c}{\lambda^{4}} \cdot \frac{1}{e^{\frac{hc}{\lambda k T}} - 1}\, \Delta \lambda \quad \left[ \frac{\text{aantal fotonen}}{m^{2}\, s} \right]

(5)

Fotonen met een energie $E_p$ groter dan de bandgap-energie $E_g$ leveren een nuttige energie $E_g$ . Dit is het geval voor fotonen met een golflengte $\lambda < hc/E_g$ . Fotonen met een energie $E_p$ kleiner dan de bandgap-energie $E_g$ leveren een nuttige energie 0. Dit is het geval voor fotonen met een golflengte $\lambda > hc/E_g$ . De totale energie geleverd door de zonnecel bedraagt dus

E_{zc, E_g, T} = 2\pi c E_g \int_{0}^{hc/E_g} \frac{1}{\lambda^{4} \left( e^{\frac{hc}{\lambda k T}} - 1 \right)} \, d\lambda \approx 2\pi c E_g \sum_{i=1}^{hc/(E_g \Delta\lambda_0)} \frac{1}{\lambda_i^{4} \left( e^{ \frac{hc}{\lambda_i k T} } - 1 \right)} \, \Delta\lambda_0

(6)

met $\lambda_i = i \cdot \Delta \lambda_0$ .

{exercise}

Schrijf een Python programma dat bovenstaande Riemann som (4) uitrekent voor een gegeven $E_g$ en

$\Delta \lambda_0=10\mathrm{nm}$
$T=5800\mathrm{K}$

Bereken en plot als functie van $E_g$ voor $0.5 \mathrm{eV} \leq E_g \leq 2 \mathrm{eV}$ .
Wat is, volgens deze modelberekeningen, de optimale bandgap-energie (in eV) voor een zonnecel?

#constanten:
h = 6.626e-34 # Planck's constant (J·s)
c = 2.998e8 # speed of light (m/s)
k = 1.381e-23 # Boltzmann constant (J/K)
T = 5800
a = h * c / (k * T) # vereenvoudiging van de som

def E_T_2(l0,  Eg):
    lambda_g = h * c / Eg
    N = int(lambda_g / l0)
    return (2*np.pi*c*Eg)*np.sum(((1)/((l0*i)**4*(np.exp(a/(l0*i))-1)))*l0 for i in range(1,N+1))
Eg_ev = 1.1
Eg = 1.1*1.602*10**(-19)
b_1 = E_T_2(10**-9, Eg)


Eg_t = np.linspace(0.5,2,100)
Eg_test = Eg_t*1.602*10**(-19)

E_tot = [E_T_2(10**-9, Eg) for Eg in Eg_test]
plt.figure()
plt.plot(Eg_t,E_tot,"r-")
plt.plot(Eg_ev,b_1,"k.")
plt.xlabel("Bandgap E_g [eV]")
plt.ylabel("Totale nuttige energie")
plt.show()

C:\Users\korne\AppData\Local\Temp\ipykernel_14412\3085693558.py:11: DeprecationWarning: Calling np.sum(generator) is deprecated, and in the future will give a different result. Use np.sum(np.fromiter(generator)) or the python sum builtin instead.
  return (2*np.pi*c*Eg)*np.sum(((1)/((l0*i)**4*(np.exp(a/(l0*i))-1)))*l0 for i in range(1,N+1))
C:\Users\korne\AppData\Local\Temp\ipykernel_14412\3085693558.py:11: RuntimeWarning: overflow encountered in exp
  return (2*np.pi*c*Eg)*np.sum(((1)/((l0*i)**4*(np.exp(a/(l0*i))-1)))*l0 for i in range(1,N+1))

E_g = 1.602E-19 * np.linspace(0.5, 2, 1000)
E_g_T = np.array([])

#your code/answer