Simulation Model A simulation model was constructed of a PBEC (photo-bioelectrochemical cell) using AnyLogic. It is to be used for illustrational purposes.
Click on the link below to run the simulation model in your browser (Note: Java is required) If your browser doesn't support Java, see the illustration of the model in action below:
"Black box" strategy
In order to reduce the complexity of the simulation model, a “black box” strategy was followed. A “black box model” requires specific information as input where after the system utilizes pre-programmed logic to return output to the user. Many “black box models” can be used which uses each other’s output (Investopedia, 2016). Making use of a “black box model” ensured that functional areas relevant to the proposed model were isolated without neglecting variability, as a result of the influence of other parameters.
Figure 1: A simple “black box” model
The first “black box” to be addressed relates to the influence of the material type used for the casing of the battery. Perspex and glass were considered as possibilities. Absorption properties of material influence transmitted light. This will in return affect the light reaching the attached photosystems. Perspex of 4.5 mm thick has a solar light transmittance of 37% and visible light transmittance of 33%. Glass has a solar transmittance of 90% and visible light transmittance of 91% (Altuglas International, 2016).
Figure 2: The first “black box”
The second "black box"
The second “black box” in the model represents the process of photosynthesis occurring in photosystems within the attached thylakoid membranes.
Figure 3: Photosystem II electron flux and various fluorescence yields as a function of light intensity
The estimate of electron flux
(open diamonds) is of importance in describing what happens in the second “black box” of this simulation model. A line of best fit was retrieved from the information in Figure 3 and transformed to a polynomial function plotted in Figure 3, described by the following equation:
Figure 4: Electron flux of photosystem II as a function of light intensity
The above figure describes the electron flux obtained from a single photosystem II. There is a single photosystem II inside each
thick thylakoid (Kirschner Lab, 2001).
Figure 5: The second "black box"
The third "black box"
Fick’s diffusion laws are used to describe the movement of electrons from the photosystems to the graphene electrode. Fick’s first law is given as:
Where
Equation [2] can be rearranged into Equation [3]:
Where
The concentration of electrons at point 1 is the changing electron flux produced by the thylakoid membranes and obtained from the second “black box”. It is assumed that the electron concentration at point 2, which is the surface area of the graphene electrode, is infinitely small. The exposed surface area (A) of the graphene electrode is where the thylakoid membranes are attached. The distance (x) resembles the space between the thylakoid membrane and the surface area of the graphene. Using aptamers, this distance amounts to
(Barendse, 2016). Without the use of aptamers an 50% increase in distance is assumed.
Where
Boltzmann’s constant of
is used. An absolute temperature of 273.16 K is assumed as the battery will be operating at room temperature. The electric mobility in water of
is used. The diffusion coefficient
is calculated as
Figure 6: The third "black box"
In the above Figure 6, electron flux
is used as input for the second black box, where after equation [3] and [4] describe the diffusion of electrons from the thylakoid membranes to the graphene electrode. An output of a number of electrons per second is obtained.
The Overall Model Logic
Combining all of the “black boxes” results in the model logic as illustrated in Figure 7. Changeable parameters include: light intensity
casing material (Perspex or glass), electrode size and whether or not to select thylakoid attachment by making use of aptamers. The main measurable is current (A).
Figure 7: The combination of "black boxes" used in the model
User Interface and Animation
Buttons and sliders are added to allow the user to change specific parameters. Graphs of both the current (A) measured and the influence of light intensity on the electron flux produced are drawn dynamically. A 3D illustration of the PBEC is also presented which is animated as the user changes parameters.
Model Results
In order to illustrate the results of the simulation model, 3D graphs were compiled of the output obtained from each “black box”. Following is a graph that displays the light intensity obtained after using Perspex and glass of both 0.2 and 0.5 mm thick at various input light intensities. As illustrated, the highest yield is obtained by using 0.2 mm thick glass at a light intensity of Figure 8: Model results of adapted light intensity related to different material types
The process of photosynthesis occurring in the thylakoid membranes is represented by the second “black box”. Following is a graph of the obtained electron flux as a result of both the electrode size (surface area) and light intensity. The highest yield is observed in the region using an electrode surface area of and light intensity of . Figure 9: Model results of electron flux as a result of electrode size and light intensity
The third “black box” represents electron transport from the thylakoid membranes to the graphene electrode via diffusion. Following is a graph of the adjusted electron flux after electron diffusion occurred. The highest output electron flux is obtained by using an input electron flux of μ mole electrons.m−2.s−1; with the smallest possible distance between the electrode and the thylakoid membranes. This shortest distance is obtained by making use of aptamers of length Figure 10: Model results of adjusted electron flux after diffusion occurred
After multiple runs, the highest yield resulting in a current of
was obtained with the following parameter settings: Light intensity of
using glass as casing material; keeping the electrode surface area at
making use of aptamers to attach thylakoids to the electrode.
AnyLogic is a simulation modelling tool that supports multiple simulation methods such as system dynamics, discrete event and agent based modelling. AnyLogic consists of a graphical modelling language working with Java code extensions.
A changed light intensity
is obtained as an output as seen in Figure 2 below.
Photosynthetic electron transfer’s relative proton to electron stoichiometries (H+/e- ratios) was obtained from living tobacco plant leaves under steady-state illumination. Electron transfer fluxes as well as turnover rates of linear electron transfer through the cytochrome (cyt) b6f complex were estimated (Colette A. Sacksteder, 2000). In Figure 3, fluorescence yields and electron flux of photosystem II is plotted. Yields from fluorescence during saturation pulses,
(closed squares) and yields from fluorescence in the steady-state,
>
(open circles), were obtained at varying light intensities. Photosystem II and associated light-harvesting complexes’
quantum yield and an estimate of the electron flux of photosystem II
were calculated, as described in the above mentioned report.
This was used to calculate the quantity of photosystems II possibly located on the surface area of the graphene electrode.
As seen in the above Figure 5, light intensity
is used as input for the first black box, where after equation [1] describes processes happening inside the box and gives an output of electron flux
The diffusion coefficient of electrons in solution is calculated using Einstein’s relation, also known as Einstein–Smoluchowski relation: