Dynamic Modeling of Anode Function in Enzyme-Based Biofuel Cells Using High Mediator Concentration

The working principle of enzyme-based biofuel cells (EBFCs) is the same as that of conventional fuel cells. In an EBFC system, the electricity-production process is very intricate. Analysis requires a mathematical model that can adequately describe the EBFC and predict its performance. This paper develops a dynamic model simulating the discharge performance of the anode for which supported glucose oxidase and mediator immobilize in the EBFC. The dynamic transport behavior of substrate, redox state (ROS) of enzyme, enzyme-substrate complex, and the mediator creates different potential changes inside the anode. The potential-step method illustrates the dynamic phenomena of substrate diffusion, ROS of enzyme, production of enzyme-substrate complex, and reduction of the mediator with different potential changes.


Notations
[ ] concentration in the system [

Introduction
The operating principle of biofuel cells (BFCs) (also called biological fuel cells) is similar to that of conventional fuel cells, regarded as potential substitutes for fossil fuels.BFCs can reduce greenhouse gas emissions and increase energy security.Enzyme-based biofuel cells (EBFCs) are one type of BFC that convert chemical energy into electrical energy.EBFCs are distinguishable from conventional fuel cells by the use of biomass (such as carbohydrates, cellobiose, ethanol, fructose, glucose, glycerol, hydrogen, lactose, methanol, pyruvate, sugars) and specific enzymes known as biocatalysts [1][2][3][4][5][6].Enzymes immobilized on bio-electrodes facilitate reuse of biocatalysts.In general, the substrate is the biomass oxidized at the anode to produce hydrogen protons and electrons.At the cathode, the oxidant (usually oxygen) reacts with electrons and hydrogen protons to generate water.Three critical challenges in developing direct EBFCs are short lifetime, lower current density, and poor power density.These are related to enzyme stability, electron transfer rate, enzyme loading, and inefficient electron conduction between biocatalysts and bio-electrodes [6,7].Also, some literature reports electron transfer from enzymes to electrodes as an important factor [2,[8][9][10][11][12][13][14][15], but mediators generally shuttle electrons from enzymes to electrodes [16,17].
Nanobiocatalysts overcome the disadvantages in EBFCs, improving electron transfer in the electrode.Many nano-structured materials, such as mesoporous media, nanoparticles, nanofibers, and nanotubes, have been proven to be efficient hosts for the immobilization of enzymes.It is evident that when nanostructures of conductive materials are used, the large surface area of these nanomaterials can increase enzyme loading and facilitate reaction kinetics, thus improving EBFC power density.Nanocarbon materials, such as carbon black, carbon nanotubes, carbon nanoballs and carbon aero-gel particles, form three-dimensional electrodes [14,[17][18][19][20][21][22][23][24][25].In these electrodes, mediators or redox polymers are not chemically immobilized to the nanocarbon materials [14].
Figure 1 shows the schematic structure of the simulated bio-electrode in this study.Small biomolecules such as flavin, mediator, and quinone compounds can transfer electrons directly to an electrode, allowing the study of the electrochemical compounds' redox behavior [25].As mediators, various organic, inorganic compounds and some redox proteins have been used [26].Benzoquinone, hydroquinone, and pyrroloquinoline quinone also may function as mediators for glucose oxidase [27].Ferri/ferrocyanide [Fe(CN) 6 3− ] is highly diffusible and can be easily reduced to its ferrous counterpart by the well-defined reversible reaction simultaneous to increase the redox potential of the solution [16,28], one of the most commonly used inorganic mediators.Large biomolecules, conducting polymers (CPs), proteins, and enzymes are usually sluggish in the electron transfer reaction at an ordinary electrode.Relatively small redox proteins, such as cytochrome c, can exchange electrons rapidly with a promoter-modified electrode or an ordinary electrode under appropriate conditions.The enzymatic reaction for the oxidation or reduction of the substrate (S) is linked to the electrochemical reduction or oxidation of the mediator (M red /M ox ) using the electrode as a final electron acceptor or a donor.Enzyme-catalyzed electrochemical oxidation or reduction of the substrate is called bioelectrocatalysis [29].Figure 2 shows a typical electron and proton kinetic in mediated bioelectrocatalysis of an EBFC system.This reaction is very useful for constructing BFCs, biosensors, bioreactors and for amplified detection of mediators as analytes and enzyme kinetic measurements [26].To determine kinetic parameters for the bioelectrocatalytic reaction in the film, the system was modeled as described in detail elsewhere [30,31].The total reactions occurring within the anode can be written as follows: The reaction mechanism is given as follows [2,15,17,18,26,[31][32][33][34][35][36][37][38][39]: where (E ox , E red , ES, M ox , and M red ) are the oxidized and reduced forms of enzymes, enzyme-substrate complexes; and oxidized and reduced forms of mediator molecules, and (S) and (P) are the substrates (reactants) and products of enzymatic reaction.k 1b and k 1f are backward and forward rate constants, respectively (Equation 2); k 2f is the product formation rate constant; and k 3b and k 3f are back and forward rate constants, respectively (Equation 3).The simplest description of steady-state enzyme kinetics is based on the reports of previous works.This treatment based on the assumption that the substrate forms a complex (the enzyme-substrate complex) (ES) with the enzyme in a reversible step, that is maintained at equilibrium between the enzyme (E), substrate (S), and enzyme-substrate complex (ES).Irreversible breakdown of the enzyme-substrate complex yields the product (P) [32].For an oxidase such as glucose oxidase, the most studied redox enzyme, (E ox ) is active towards the substrate (S) and produces the inactive form (E red ) which can be reactivated by oxidation with the cosubstrate (M ox ) to form the reduced mediator (M red ), which subsequently is reoxidized at the electrode surface.For a glucose oxidase, the active enzyme is the reduced form and the oxidized form of mediator which is recycled at the cathode.The general mechanism for a mediated enzyme electrode is displayed in Equations (2)(3)(4) where the oxidized (reduced) form of the enzyme (E ox ) reacts in response to the substrate (S) to form an enzyme-substrate complex (ES), with further decomposition that yields the reduced (oxidized) form of the enzyme (E red ) and the product (P).The oxidized (reduced) form of the mediator (M ox ) reoxidizes E red (E ox ) to the active form of the enzyme E ox (E red ) and the reduced (oxidized) mediator (M red ).At the electrode surface, (M red ) is converted to (M ox ) Equation ( 4), which diffuses out from the electrode and reacts with the enzyme [3,21,26,30,32,35].A summary of recent key electrode models in EFCs is presented in Table 1.They are discussed in detail in the following sections.More information can be obtained by sweeping the potential with different scan rates from 0.001 V/s to 1000 V/s [36].To investigate the discharge performance of biofuel cells, the rate of linear scan voltammogram (LSV) recording the steady state of the enzyme reaction was very slow (1-5 mV/s) [37,38].Bartlett et al. [37] also used cyclic voltammetry to extract the kinetics parameters with a scan rate of 5 mV/s by steady-state modeling for substrate-enzyme reaction.LSV was used to study the kinetic parameters of enzyme and mediator reactions with a high scan rate 20 mV/s [37].In this region, the discharge performance and the electrode kinetics are transient transport phenomena.However the transient transport phenomena of enzyme-mediator reactions in the biofuel cell have not received enough attention regarding their relation to mediator dispersion.The function of a mediator is to enhance the electronic properties of electrode and reduce the path of electronic migration.The mediator dispersion depends on the concentration of the mediator.The literature contains several modeling strategies for enzymatic electrodes [46][47][48][49] and enzymatic biosensors [32], but modeling of enzymatic electrodes has been researched more thoroughly.Mathematical models of enzymatic electrodes in the literature usually consider the enzymatic reaction and the material balance of species which are participating in the enzymatic reaction.Recent literature on bioelectroanalytical systems of kinetic models rarely discusses the diffusion [42] and potential [44] of solution in the system.Enhancing the redox reaction of mediator is most effective in upgrading the performance for EBFCs.It is very important to understand the effects of mediators on the performance by a model to evaluate various electrode conditions.The composite model has been demonstrated to realistic for cathode [31,35], and is applied to the mediated oxygen biocathode based on macrohomogeneous and thin-film viewpoint.However, reactions that dominate degrees of redox state (ROS) of enzyme and then determine the discharge of EBFC are unknown.Hence, this study investigates the process of operating anodes to evaluate the rate determining process, a dynamic model to diagnose EBFC systems based on diffusion, solid and liquid potential supported glucose substrate, glucose oxidase, potassium ferricyanide, and pH 7.0 phosphate buffer (0.1 M), (all anode interactions in EBFC systems).

Enzyme Electrode Preparation
The enzyme composite electrodes were prepared as previously report [25].Commercial carbon paper cut into 10 mm × 10 mm squares was used in both the cathode and the anode.The anodic solution and cathodic solution were made as follows: 20 U/mL GOx, 20 mM and 40 mM of [Fe(CN) 6 3− ], and 20 mM polyaniline.The above materials were mixed with 0.1 M phosphate buffer solution (PBS) pH 7 for the anode electrode.20 U/mL LAc, 20 mM and 40 mM of 2,2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid), and 20 mM polyaniline, were mixed with 0.1 M citric buffer solution (CBS) pH 5 for athode electrode.Two mixture solutions, anode and cathode, were pipetted onto the hydrophilic carbon paper 30 s in air at ambient temperature prior to immobilization of the carbon paper.

Electrochemical Measurements
The power output performance of the tested EBFCs was characterized by two-electrode electrochemical measurements containing 10 mM glucose, used 0.1 M PBS pH 7, and 0.1 M CBS pH 5 as the electrolyte solution at working volume of 15 mL, respectively.A Nafion 117 membrane was inserted to separate the two half cells.The resistance of circuit was 200 ohmic (Ω).The distance between the anode and cathode was 10 mm.The 24-bit differential analogy module received the voltage signal from the two-electrodes.The signal was send to a computer and analyzes using the Laboratory Virtual Instrument Engineering Workbench (LabVIEW) evaluation software.All electrochemical experiments were carried out at temperature 310 K [25].

Mathematical Formulation
Consider the EBFC in a square area of length L filled with enzyme-substrate-mediator complex solution (Figure 1), where coordinate x is the distance measured along the left vertical wall.The length L dimension is 15 × 10 −5 m.According to Equations ( 2) to (4), the model considers eight species: substrate (S), oxidized enzyme (E ox ), reduced enzyme (E red ), enzyme-substrate complex (ES), oxidized mediator (M ox ), reduced mediator (M red ), hydrogen ions (H + ) and product (P).This study also considers the effect of the solid conductivity of mediators on electrode performance.Table 1 shows the comparision of the model in the study and those in the literature.The mass balance equation for the substrate, hydrogen ions, and product is given by: where D j and R j denote the diffusion coefficient for species and the source term for reaction, and j denotes the chemical species j, as listed in Table 2.The initial condition (I.C.) and boundary condition (B.C.) for Equation 5 are as follows: where the values of BC1 j and BC2 j are as listed in Table 2.The rate expression for oxidized enzyme (E ox ), reduced enzyme (E red ), enzyme-substrate complex (ES), oxidized mediator (M ox ), or reduced mediator (M red ) is given as: The initial values for Equation 6 ([j] = [j] i ) are listed in Table 3.The local current density inside electrode (i t ) can be expressed as: in which a is the surface area of the electrode per electrode volume, i o denotes exchange current density, and β denotes the anode transfer coefficient.Applied potential can be written as follows: (10) in which l  and s  denote potential in the liquid phase and the solid phase, respectively.
The charge conservation can be expressed as: in which θ and k s are the porosity and conductivity of the solid phase, respectively.The boundary conditions are: where E A is the applied potential.The liquid potential related to the electronic current is given by: where k l is the conductivity of the liquid phase.The boundary conditions are: The total current density calculated can be expressed as follows: Table 2. Source terms for species and charge conservation equations.

Species j
Source terms a These compounds are immobilized on the electrode.

Results and Discussion
The application of the low and high mediators of anode and cathode bioelectrodes for the EBFC system has been demonstrated during GOx and LAc bioelectrode testing under a galvanostatic regime.Figure 3 shows a plot of power density vs. cell voltage using different concentration of mediators (20 mM and 40 mM).After 9 h, the high mediator analysis has shown that the cell electro-oxidation current of glucose appears at 0.09 V with a power density of 45.46 μW/cm 2 ; for the low mediator, results were 0.06 V with a power density of 20.18 μW/cm 2 .Experimental results can be obtained high concentration of mediator better than the low concentration of mediator.Therefore, a model was developed and evaluated through computational simulation series by comparing them to experimental data.The results of simulation demonstrate that a good fitting might be learned the relationship between mediator and EBFC.According to the enzymatic reaction in the BFC in Equations ( 2)-( 4), eight compounds (E ox , E red , ES, M ox , M red , H + , S, and P) must be considered to describe the enzymatic reaction.Figure 1 shows the reactive schematic in EBFC.This is a solid reaction, and the diffusion problem during the reaction should be considered.The mass balance and rate equation were expressed in Equations ( 5)- (13).Table 3 lists simulated parameters.The COMSOL Multiphysics ® package software, which is based on the finite-element method, calculated the simulated results.Solutions converge when the relative error tolerance is less than 1.0 × 10 −6 .

Current Density, Potential and Polarization Loss
When the applied potential deviates from the equilibrium potential (E e = −0.2V in this study), the electrode creates voltage losses.The voltage losses include activation loss, ohmic loss, and concentration loss.Figure 4 illustrates the relationship between current density and voltage (IV curves) using typical anode efficiency.The calculation was carried out with a slow linear sweep of potential (0.1 mV/s).A low potential sweep ensures a steady state result.If only one source of voltage loss is attributable to kinetics (activation loss), the relationship between current density and potential is shown in region (a).Activation voltage loss is related to the intrinsic activity of enzymes.If the applied potential increases to region (b), it experiences voltage losses owing to activation loss and ohmic loss.In this region, the main controlling factor is the conductivity of mediators and electrolytes.If the applied potential is set at a higher potential in region (c), activation loss, ohmic loss, and concentration loss occur, and mass transfer and ohmic resistance dominate electrochemical behavior.Relevant literature has not explored the reaction, diffusion, and ROS of enzymes, the generation of enzyme-substrate complexes, and the electrochemical behavior of mediator-associated micro-nature.Figure 5(a) shows the current density with time at low applied potential.The applied potential is set at 0.0 V; the current density rapidly rises, and with a slow rate when the time is approximately 50 s.The current density, which tends to be a constant value such as 0.0 V, is approximately 0.043 mA/cm 2 .When the applied potential is far from the equilibrium potential (−0.2 V), current density responses occur [Figure 5(b)].When the potential is 0.1 V, the initial current density is still in an upward trend, but the current density first increases and then decreases slowly.The response behavior of the current density declines at high values of approximately 0.5 mA/cm 2 to 0.36 mA/cm 2 , and after 70 s, it stabilizes at a constant value.When the applied potential increases to 0.15 V to 0.2 V, the trends are similar (first an increase and then a decrease behavior).When the applied potential is 0.25 V, the current density first increases rapidly and then dramatically decreases.The operation time for steady state is approximately 30 s.When the current density decreases from 1.1 mA/cm 2 to 0.48 mA/cm 2 , the maximum value decreases by 25%.When the applied potential increases to 0.25 V, current density jumps to 2.1 mA/cm 2 , then quickly drops to 0.7 mA/cm 2 , and the maximum value declines 35%.Researchers chose −0.15 V, 0 V, 0.05 V, 0.15 V as examples to study the variation of current density in EBFC.

Distribution of Substrate on the Electrode
In the EBFC reaction system, the substrate must first diffuse from the bulk solution to the active site, and then react with the enzyme.Figure 6(a-d) simulate the substrate distribution across the entire anode with different applied potentials, step by step.
Figure 6(a) illustrates substrate distribution after a potential step from equilibrium potential (−0.2 V) to −0.15 V at time points of 5 s, 30 s, 55 s, 80 s, and 105 s.Substrate concentration increases to about 1.9 mM at a time of 105 s.That is, it takes almost 105 s to establish a steady-state substrate distribution.The substrate is consumed through electrochemical oxidation, and its concentration drops minutely from 2.0 to 1.9 mM. Figure 6(b) shows the substrate concentration profile as the potential moves toward 0.0 V.The maximum substrate value is approximately 1.23 mM at 105 s.This value is about 60% of the substrate in the bulk solution.Substrate consumption increases with electrochemical oxidation when the applied potential is held at a higher value.Figure 6(c) shows the substrate concentration profile at various time levels for a potential set at 0.05 V.The substrate slowly and steadily increases, and reaches a steady state after 105 s.The substrate concentration drops quickly from 2.0 mM to about 0.0 mM at the surface of the current collector (x = 0) after 105 s. Figure 6(d) shows the substrate concentration profile when the potential is held at 0.15 V.The substrate reaches a steady condition in a very short time (about 10 s), and the variation of the concentration gradient of substrate diminishes.There is no substrate in the region of 0 m to 9.0 × 10 −5 m inside the anode.

Distribution of Oxidized Enzyme on the Electrode
Efficient capacity of immobilized enzymes relates to the E ox concentration, the degree of reaction, and operation of potential or current density.Figure 7 shows the simulation results of enzyme concentration at different operation times and different applied potentials.
The applied potentials used were −0.15 V, 0.0 V, 0.05 V, and 0.15 V, corresponding to points (A) to (D) in Figure 4.The initial concentration of the oxidized enzyme E ox used was 0.05 mM. Figure 7(a) shows that the enzyme concentration mainly distributes in the thickness range of 0 m to 11 × 10 −5 m with the applied potential of −0.15 V; time was 5 s.The E ox concentration decreased with increasing operation time, and the concentration profile shifted to the left.E ox concentration approximated to 0 at t = 50 s.This occurs because when the applied potential operated at the active polarization area, the E ox concentration's degree of ROS of enzyme was influenced by the reactivity of bimolecular M ox and E red in Equation 3. When the applied potential operated in the active polarization region, mediator M ox concentration was low, resulting in decreasing ROS of enzyme of E ox concentration.Figure 7(b) shows simulation results for E ox concentration when the applied potential was controlled in the coexistence of active polarization and ohmic polarization regions (i.e., applied potential at 0.0 V).When the time was 10 s, the E ox concentration was similar to that in Figure 7 (d) concentration approached 0. Figure 7(c,d) show the curve of E ox concentration shifting to a stable value.The applied potential operates at the ohmic polarization control region as the electrochemical oxidation of the supplied M ox reaches a sufficient amount.M ox and E red bimolecular reactions are reversible.This shows that M red concentration increases with decreasing oxidation potential, and that M ox concentration increases with increasing oxidation potential.For the EBFC with a mediator, the efficiency of EBFC depends on the diffusion of reactants and the oxidation rate of the mediator.

Distribution of Enzyme-Substrate Complex on the Electrode
ES level can affect battery discharge efficiency, the amount of product, and hydrogen ion migration.Figure 8 shows the plot of ES concentration on the x-coordinate at different operation times and applied voltages.The concentration profile of the ES is a peak-type distribution that relevant literature never discusses.In Figure 8

Distribution of Mediator on the Electrode
The oxidation degree of M red correlates directly with the applied potential in BFC.M red concentration changes with the oxidation potential of the electrode.The increment of oxidation potential increases the oxidation degree of M ox .In other words, M red concentration decreases with increasing applied potential.The increment of M ox concentration helps increase the ROS of E ox concentration to enhance the enzyme reaction.Figure 9(a,b) shows the relationship between M red concentration and the x-coordinate at different operation times and applied potentials.The initial M red concentration is 10 mM.In Figure 9  When time increase exceeds 50 s, M red concentration continues to decline to the lower limit of 0.58 mM at 50 s and 0.2 mM at 70 s, repectively.This decline is mainly caused by increased oxidation potential, which reduces the concentration of mediator M red .The enzymes must be recovered with the reaction of M ox and E red .At 0.15 V potential, E ox and M red concentrations change and decline in the same area; the major changes in the concentration of E ox are in the region of 8 × 10 −5 m-12 × 10 −5 m, and the concentration of M red changes in the region of 7 × 10 −5 m-11 × 10 −5 m, as shown in Figures 7(d) and 9(d).The applied potential is higher, which causes M red to react to M ox in a mass transport controlled situation; consequently, considering the electron-transfer kinetics in theoretical treatment is unnecessary.At the same time, no M red exists in the bulk solution.This problem is then a special case of the more general problem of coupled diffusion and reaction in an immobilized enzyme layer [50].

Conclusions
This work studies the concentration profile of enzymes, mediators, and substrates in EBCS for high mediator concentration.At low potentials (e.g., 0.0 V), the current density increases with increasing operation time and reaches a stable value.The reactant at low potentials (e.g., −0.15 V or 0.0 V) easily diffuses in the interior of the reaction region.Hence, the substrate concentration in the electrode film approached that in the exterior of the film.The oxidized enzyme concentration approaches zero in the electrode film when the operation time is greater than 50 s.The ES concentration profile is a peaktype.M red concentration always remains at approximately 10 mM, so the E ox generation rate is small.However, at higher potentials (e.g., 0.25 V), current density decreases dramatically with increasing operation time.At higher potentials, the reaction rate of electrochemistry is larger than the diffusion rate of the substrate.No substrate is evident in the region of 0 × 10 −5 m-9 × 10 −5 m.The enzyme E ox concentration can be maintained at approximately 15% of initial values in the region of 0 × 10 −5 m-13 × 10 −5 m.The ES concentration profile includes a peak-platform type.The ES concentration at the platform increased when the applied potential increased.M red concentration changed in the region of mainly 10 × 10 −5 m-15 × 10 −5 m, and the lower limit was at 3.5 mM when the applied potential was 0.05 V. M red concentration decreased with increasing applied potential and even declined to 0.15 mM.From the model's prediction, it is recommended that the applied potential should be controlled between 0.05 and 0.15 V when the equilibrium potential is −0.2 V by using high mediator concentration.

Figure 1 .
Figure 1.Schematic configuration and function of the anode electrode in contact with an EBFC System.S: substrate which bulk concentration is constant; P: product; e − : electron; X: biocatalysis thin-film [25].

Figure 2 .
Figure 2. Dynamic scheme configuration and function on bioelectrodes in EBFC.E ox and E red : oxidized and reduced forms of enzyme, M ox and M red : oxidized and reduced forms of mediator, E: enzyme, M: mediator, A: anode electrode, C: cathode electrode [25].

Figure 3 .
Figure 3. Plot of power density and cell voltage at different mediators on bioelectrode.

Figure 4 .
Figure 4. IV curve for typical anodes in EBFC at 310 K.

Figure 5 .
Figure 5. Current density response of EBFC: (a) High potential state; (b) Low potential state.
(a).Similarly, the curve of E ox concentration shifted to left.The E ox concentration range for 50 s was between 0 μm and 20 μm.After 70 s, the
(a), when the applied potential is at −0.15 V and the time is 10 s, the concentration of ES exists mainly between approximately 6 × 10 −5 m and 9 × 10 −5 m, and the maximum value is approximately 0.012 mM.As the operation time increases, ES concentration also shifts to the left.After 40 s, the main peak of ES appears between 0 m and 1.5 × 10 −5 m, and the maximum of ES concentration is 0.0068 mM.When time exceeds 50 s, the ES concentration is 0.0 mM.Figure8(b)shows the simulated results of ES concentration at the applied potential of 0.0 V. in Figure8(b) are similar to those in Figure8(a) except that the plate concentration at steady state is approximately 3.0 × 10 −4 mM.Figure8(c)shows the simulated results of ES concentration when the applied potential increases to 0.05 V.The ES concentration profile combines a peak shape and a platform of 7.7 × 10 −4 mM.When the applied potential increases to 0.15 V, the platform concentration of ES in Figure8(d) increases to 4.9 × 10 −3 mM.This value is far higher (approximately six-fold = 4.9 × 10 −3 /7.7 × 10 −4 ) than that for 0.05 V.
(a), when the applied potential is at −0.05 V and the operation time is 10 s, the concentration profile for M red resides between 6 × 10 −5 m and 15 × 10 −5 m.The mediator M red concentration drops from 10 mM to 7 mM, declining approximately 30%.When the time increaseConcentration of ES [mM], M red concentration is recovered to the initial value because M ox concentration is low.That is, the oxidation rate of M red is slow.The applied potential changes from the active polarization control region to the coexistence of active polarization and ohmic polarization regions are shown in Figure9(b).When the time increases from 10 s to 40 s, the variation of M red concentration moves from the initial region of 6 × 10 −5 m-9 × 10 −5 m to the region of 0 m-2.2 × 10 −5 m.When time has exceeded 50 s, and there is movement to the 0 m-1.6 × 10 −5 m region, the M red concentration drops to 6.3 mM.When the applied potential is set at 0.05 V, the simulation results of the M red concentration are obtained in Figure9(c).The applied potential is set in the ohmic polarization control region, and M red concentration decreases with increased operation time (e.g., 50 s), dropping from 10 mM to 4.7 mM.When the time was increased to 130 s, the change of the M red concentration gradient was in the region of 0 m-2 × 10 −5 m μm.M red concentration drops from 10 mM to 8.5 mM in Figure9(c), and Figure9(d) shows the simulated result of M red concentration at an applied potential of 0.15 V.When the time is 10 s, the M red concentration changes dramatically in the region of 9 × 10 −5 m-15 × 10 −5 m and decreased to 4.2 mM.This low value is different from those obtained for −0.15 V and 0.0 V (lower limit of 7 mM).When the time increases to 30 s, M red concentration continues to decline in the 6 × 10 −5 m to 15 × 10 −5 m region.The lower limit changes from 4.2 mM at t = 10 s to 1.5 mM, decreasing by 15% from the initial value.
mM] a surface area of the electrode per electrode volume [L/m] D H diffusion coefficient for free hydrogen [m 2 /s] D p diffusion coefficient for free product [m 2 /s] D s diffusion coefficient for free substrate [m 2 /s] l  potential in liquid phase S  potential in solid phase θ volume fraction of liquid in the bioelectrode

Table 1 .
Summary of enzyme electrode models.
S/E/M: Substrate/Enzyme/Mediator; M-M: Michaels-Menten kinetics; S-S: steady state; U-S: unsteady state; CME: composite mediated electrode; E-bound: enzyme is bound; E-free: enzyme is free; M-free: mediator is free; M-bound: mediator is bound; k l /k s : conductivity of liquid phase/ conductivity of solid phase.