FC Coding is One step from 1st Principles to Solution!

Increased Productivity Example #7

Modeling a Nanostructured Solar Cell

Problem: How to develop solar cells with a new (higher) efficiency; grätzel cells?

There are many things said about what's most important for the solar cell. So what they need is a model to know what's the rate is limits for the whole system. By then they can choose what combination of parameters will give the best solar cell. The model shown is a one-dimensional non-steady state model; a start to compare it with the Laser experiments. The laser experiments are one of the things they use to predict the efficiency.

But without a model, does experiments really tell one anything? This model is only for one excitation from a laser beam and to analyze how the decay of all species are. There is a model done for steady state, but its not really working very good in practice. Simulating the non-steady state model for some time should converge to the steady state solution when there is equilibrium in the system. This means when the change of all species are zero over the film. This could be interesting to compare with other steady-state models.

When we are talking about efficiency, it should be for simulation of the whole system. Then we have to add certain things. There are continuous excitations of electrons which is the starting conditions in this model for the electrons and the excited dye. There are a few more reactions and we have to consider the other part of the solar cell which isn't contained by nano-structured TiO2.

The main thing about the efficiency is that we want as many electrons leaving the back contact which is at x=0. In the reality the electrons will go out in an outer circuit to make a full circuit. But in the Laser experiments this does not happen because the outer circuit is open. In reality we will get out a current dependent on the incident light. There are many ways to measure the efficiency. IPCE(\lambda) incident-photon-to-current efficiency says how much of the incident light was converted to external current.

Short review of the system

We have a dye sensitizer attached to nano-structured Titanium dioxide (TiO2) film. The nano-structured particles are in a dye which transports the electrons from the electrode to the dye sensitizer. Incident light at a certain wavelength excites electrons in the dye sensitizer. So what happens to this electron after the excitation? A very fast process in nanosecond scale injects the electron to the TiO2 and its making a random walk (that's what most people think its doing) to the back contact. A new electron from the dye is put in the place of the injected electron. The electron's goes through the nano-structured film to a back contact to the outer circuit and we have a total circuit.

But there are other reactions involved in the process. The excited electron can travel other ways then to the back contact like reacting with the dye or dye sensitizer. These reactions are limiting the efficiency of the cell.

Thus I thought it would be a good idea for the model to set up rate constants for all these reactions. Make a discretization along x which is the distance to the back contact. And the step through time and see how the kinetics, diffusion and the electric field is changing the concentration of the species along x for different times.

A macroscopic model for the concentration of s ( the
dye sensitizer) could look something like this:

kinetics:
ds(x,t)/dt = -k_3*s(x,t)*e(x,t)-k_4*s(x,t)*i(x,t)
(k_3 and k_4 rate constants e = electron concentration, i = iodine
conc.)

electric field E(x,t):
ds(x,t) / dt = my*ds(x,t) /dt*dE(x,t)/dt
(my = mobility for the species)

The electric field we get from integrating
concentrations of all the charged species along x.

Explanations of each color:
= the dye which is the charge carrier, giving new
electrons to the dyes and get new one at the anode. It is a redox
couple of Iodine. It can also react with the excited electrons which
gives a less good efficiency. There are also other leakage's that
contribute to decline.

= Dye molecules, the electrons of those
are excited at incident light of certain wavelengths

= The nanostructured semi conductor,
most used is TiO2, the electrons diffuse in this medium towards the
back contact.

= back contact (x=0), where the
electrons go to get to outer circuit, anode.

= end of the nanostructured film, x=8*10^(-6)

= the 'entrance' for the electrons
from outer circuit, the cathode.

We start with the species in the
solar cell
S+ = excited dye S = dye
I- = Iodine I3- =
three iodine
I0 = iodine radical I02- = di iodine radical

There are some reactions between
the species during simulation with reaction rates k1..k6
{S+} + {e-} -> {S} k_1
{S+} + {I-} -> {I0} k_2 {I0} +
{I-} -> {I02-} k_3
2{I02-} -> {I3-} + {I-} k_4 {I02-} + {e-} ->
2{I-} k_5 {I3-} + 2{e-} -> 3{I-} k_6

The concentrations of each
species is defined as
s(x,t) = {s+} i(x,t) = {I-}
e(x,t) = {e-}
w(x,t) = {I02-} q(x,t) = {I3-}
z(x,t) = {I0}

The starting conditions are (
after a laser pulse there is excitation of the dyes s(x,0) and we look
at the relaxation of all species after that)
s(x,0) = 360*10^(-9)*0.34*10^6*0.1*exp(-0.34*10^6*x)
i(x,t) = 0.5 e(x,t) = s(x,0)
w(x,t) = 0.0 q(x,t) = 0.05
z(x,t) = 0.0

The differential equations with
electric field diffusion and reactions with the diffusion constants:
Di, De, Dw, Dq, Dz; and mobility constants: my_s, my_i, my_e, my_w,
my_q we set the constants.

x will be between 0 and 8*10^(-6) which is the thickness of the film where those reactions are.

At the boundaries we need to approximate the derivatives each time step for the diffusion and the electric field. There is no flow of particles out of these boundaries so we could set the concentrations change to zero at the boundaries; i.e. for the boundaries:

So far, this is a initial value problem with six coupled nonlinear differential equations.

Future:

Making a full scaled 3 dimensional model and optimize the parameters for a optimal solar cell. Most of the parameters are adjustable, they are all dependent of the materials used. There are many different things said about what the cell efficiency really depends on. Some even say that the important thing is to have as good a cathode as possible and others say its something completely different. A model is needed to guide the research for a better solar cell.

Increased Productivity Example #7
Source Code:

Problem SolarCel include 'SolarCel.inc' C Ok we start with the species in the solar cell: C ---------------- C S+ = excited dye S = dye I02- = di idodine radical C I- = Ioidine I3- = three iodine I0 = iodine radical C we have some reaktions between the species during simulation with C reaction rates k1...k6: C ---------------- C {S+} + {e-} -> {S} k1 C {S+} + {I-} -> {I0} k2 C {I0} + {I-} -> {I02-} k3 C 2{I02-} -> {I3-} + {I-} k4 C {I02-} + {e-} -> 2{I-} k5 C {I3-} + 2{e-} -> 3{I-} k6 C the concentrations of each species is defined as: C ---------------- C s(x,t) = {s+} i(x,t) = {I-} e(x,t) = {e-} C w(x,t) = {I02-} q(x,t) = {I3-} z(x,t) = {I0} C The differential equations with electric field diffusion and reactions C with the diffusion constants: Di, De, Dw, Dq, Dz and mobility C constants: mys, myi, mye, myw, myq we set the constants: C diffusion:
constDs = 1.5*1.e-9 : constDe = 200*1.e-9
constDi = constDs : constDq = constDs : constDw = constDs
print *,'Const.',constDs,constDi,constDe,constDw,constDq,constDz
C mobility:
mys = 1.5*1.e-9*1.602*1.e-19/(1.38*1.e-23*273)
myi = mys: myw = mys: myq = mys
mye = 200*1.e-9*1.602*1.e-19/(1.38*1.e-23*273)
print *, 'My.', mys, myi, mye, myw, myq, myz
C rate constants:
k1 = 1.0 * 1.e-6: k2 = 3 * k1: k3 = 4 * k1: k4 = 5 * k1
k5 = 6 * k1: k6 = 9 * k1
print *, 'Ks.', k1, k2, k3, k4, k5, k6
C x will be between 0 and 8*1.e-6 which is the thickness of the
C film where those reactions are:
xfinal = 8*1.e-6: xprint = xfinal / 100: dx = xprint / 10
tfinal = 1.e2: tprint = tfinal / 100: dt = tprint / 10
C At the boundaries I suppose we need to approximate the
C derivativeseach time step for the diffusion and the electric field.
C There is no flow of particles out of these boundaries
C so we could set the concentrations change to zero at
C the boundaries; i.e for the boundaries.
didx=0: dedx=0: dwdx=0: dqdx=0: dzdx=0
C the starting conditions are ( after a laser pulse there is
C excitation of the dyes s(x,0) and we look at the relaxation
C of all species after that):
initiate JANUS; for distance; equations
* dsdx/x, d2idx/didx, didx/x, d2edx/dedx, dedx/x, d2wdx/dwdx,
* dwdx/x, d2qdx/dqdx, dqdx/x, d2zdx/dzdx, dzdx/x, dEsumdx/Esum;
* of x; step dx; to xf;
print *,' TIME DSDT S DIDT I'
xf=xprint
do while (xf .le. xfinal)
integrate distance; by JANUS
print '(7(1pg13.5))', x, s, i, e, w, q, z
C @curves('plot')
xf=xf+xprint
end do
C @show('plot')
end
model distance
include 'SolarCel.inc'
s = 360*1.e-9*0.34*10**6*0.1*exp(-0.34*10**6*x)
i = 0.5 : e = s : w = 0. : q = 0.05 : z = 0.
C movement from diffusion:
dsdt = 0. ! stationary
didt = constDi * d2idx : dedt = constDe * d2edx
dwdt = constDw * d2wdx : dqdt = constDq * d2qdx
dzdt = constDz * d2zdx
initiate ATHENA; for ide; equations dsdt/t, didt/t,
* dedt/t, dwdt/t, dqdt/t, dzdt/t; of t; step dt; to tf;
print *,' X TIME DSDT S DIDT I'
tf=tp
do while (tf .le. tfinal)
integrate ide; by ATHENA
print '(6(1pg13.5))', x, t, dsdt, s, didt, i
tf=tf+tp
end do
Find dsdx, didx, dedx, dwdx, dqdx, dzdx, dEsumdx; * in eForce; by AJAX( cntrl1); * to match xs, xi, xe, xw, xq, xz, xEsum C @show('plot') end model ide ! Implicit Differential Equations include 'SolarCel.inc' Find dsdt, didt, dedt, dwdt, dqdt, dzdt; * in kinetics; by AJAX( cntrl1); * to match ts, ti, te, tw, tq, tz end model kinetics include 'SolarCel.inc' C kinetics:
ts = dsdt - (- k1 * s * e - k2 * s * I)
ti = didt - (- k1 * s * e - k3 * i * z)
te = dedt - (- k1 * s * e - k5 * w * e -k6 * q * e)
tw = dwdt - (- k5 * w * e + k3 * i * z - 2 * k4 * w**2)
tq = dqdt - ( k4 * w**2 - k6 * q * e)
tz = dzdt - (- k3 * i * z + k2 * s * i)
end
model eForce
include 'SolarCel.inc'
C movement from electric force:
xs = dsdt - (mys * s * dEsumdx + mys * Esum * dsdx)
xi = didt - (myi * i * dEsumdx + myi * Esum * didx)
xe = dedt - (mye * e * dEsumdx + mye * Esum * dedx)
xw = dwdt - (myw * w * dEsumdx + myw * Esum * dwdx)
xq = dqdt - (myq * q * dEsumdx + myq * Esum * dqdx)
xz = dzdt - 0 ! not charged
! next comes from Poisson's equation
xEsum = dEsumdx - (s + i + e + w + q) ! is Objective xEsum = 0 ?
C I guess this is a initial value problem with six coupled nonlinear
C differential equations. /Jarl
end
controller cntrl1( AJAX)
summary=0
end

<a
href="http://fortranCalculus.info/example/solar-cell.html"><img
style="float:left; width:100px"
src="http://fortranCalculus.info/image/solar-pic1.png"/>
<strong>Math Modeling & Simulating a Nanostructured Solar
Cell</strong> </a>, Energy Solar, Higher Efficiency Solar
Cells.