A Catalytic Effectiveness Factor for a Microbial Electrolysis Cell Biofilm Model

The aim of this work is to propose a methodology to obtain an effectiveness factor for biofilm in a microbial electrolysis cell (MEC) system and use it to reduce a partial differential equation (PDE) biofilm MEC model to an ordinary differential equation (ODE) MEC model. The biofilm mass balances of the different species are considered. In addition, it is considered that all the involved microorganisms are attached to the anodic biological film. Three effectiveness factors are obtained from partial differential equations describing the spatial distributions of potential and substrate in the biofilm. Then, a model reduction is carried out using the global mass balances of the different species in the system. The reduced model with three uncertain but bounded effectiveness factors is evaluated numerically and analyzed in the sense of stability and parametric sensibility to demonstrate its applicability. The reduced ODE model is compared with a validated model taken from the literature, and the results are in good agreement. The biofilm effectiveness factor in MEC systems can be extended to the reduction of PDE models to obtain ODE models that are commonly used in optimization and control problems.


Introduction
Bioelectrochemical systems (BES) have recently gained significant interest due to their potential for energy recovery from renewable soluble organic matter. In particular, microbial electrolysis cells (MEC) provide a new alternative to treat wastewater and produce hydrogen simultaneously [1].
Through reactions catalyzed in the anode electrode by electroactive microorganisms within an active biofilm, the removal of organic compounds from agroindustrial wastewaters is carried out during the oxidation of organic matter under anaerobic conditions [1,2]. Protons, electrons, and CO 2 are generated in this process [1]. The electrons circulate from the anode surface to the cathode by applying a small additional external electrical potential (0.3-1.0 V) between both electrodes [2]. Protons (hydrogen ions) are transferred to the cathode electrode through a membrane, after which the electrical current (electrons) is used to reduce the protons to biohydrogen [1,2].
Unfortunately, despite the recent attention and growing research interest in using MECs for wastewater treatment, their development has not yet moved beyond laboratoryand pilot-scale implementations, and their adoption by the industry remains limited [3][4][5]. This is partly due to the significant challenges posed by the step towards the industrial scale [6][7][8][9][10]. For instance, the impact of MEC operating conditions [11], the optimization of the process design for its integration into more complex systems, the establishment of critical performance ranges, and many other operational and empirical issues; as well as Table 1. Microbial electrolysis cell model classification based on ordinary differential algebraic equations structure. Substrate  2OD  OD  OD  OD  EM  OD  OD  OD  OD  OD  FM  OD  MM  3OD  2OD  2OD  OD  The significant issues regarding ODE models are the oversimplification or nonconsideration of spatial variation. Despite providing a spatial variation in substrate and potential at steady state, reported PDE models do not take into account the variation in biofilm thickness. To the best of the authors' knowledge, there is only one MEC contribution (3AE+4ODE+2PDE) that models the variation in biofilm thickness [26].

[21] [22] [23] [24] [25]
On the other hand, multidimensional PDE models tend to require long solution times and large computational resources, and, despite recent advances and improvements in multiprocessing capabilities and programming, the use of such approaches is minimal. Therefore, the current trend is to simplify 3D models [8]. Reaction-diffusion biofilm models accurately describe biofilm and could be recommended for rigorous performance analysis and design [7]. In this regard, the reduction of a PDE MEC model should include biofilm dynamics but reduce to an appropriate ODE MEC model suitable for real-time process monitoring and control applications.
Due to the incorporation of the hydrodynamical effect in the modeling of biomass detachment in anaerobic digestion processes [29], by means of the parameter α, several contributions including this phenomenon have been reported in ODE MEC models. The concentration of microorganisms in biofilm is mainly a function of growth kinetics. However, additional biofilm mass retention terms could be included, for instance: (i) constant parameter [22,24]; (ii) continuous AE [30]; and (iii) piecewise constant AE [30], where the maximum value is an upper-bounded constant related to the maximum mass retention (i.e., X max ).
Biofilm growth and retention are modeled with constants or AEs. The approach has the advantage of simplicity and is a reliable alternative for optimization and control applications [25,31]. However, scaling up is limited if hydrodynamics does not hold. On the other hand, the dynamic description of biofilm on PDE models provides information about concentration and potential in biofilm, which could be important in scaling up [15,28]. However, the use of robust PDE models is still limited in scale-up design and, in some cases, prohibitive in process optimization and control applications. Therefore, an ODE model with a biofilm variable state can be easier to implement than biofilm PDE models [8]. In this sense, an effectiveness factor (EF) for MEC biofilm could address two apparently conflicting goals: (i) the description of mass transfer diffusion and electrical properties in biofilm from the PDE MEC model and (ii) a simple description of biofilm in the ODE MEC model. PDE models are rarely used for design, optimization, or control purposes due to the difficulty of finding a computationally inexpensive solution [32]. To address this problem, an alternative approach is the use of the aforementioned EF concept [32][33][34]. Then, the reaction rate in a catalyst is expressed as function of reaction rates at surface/bulk conditions and the effectiveness factors [32]. Indeed, the EF is a useful tool for the study of heterogeneous reaction systems. Moreover, the EF has been reported in chemical and biochemical engineering literature for: catalyst particles [32], immobilized enzyme catalysts [35], immobilized growing cell systems [36], biofilm growth and maintenance [37][38][39][40], biocatalytic membrane reactors [41], and hollow fiber membrane bioreactors [42,43]. Microbial growth kinetics have been modeled with: (i) zero-, first-, and second-order kinetics [37][38][39]41]; (ii) Monod kinetics [37,39,41,44]; and (iii) Michaelis Menten kinetics [35,40,42,43]. To the best of the authors' knowledge, there is no contribution in bioelectrochemical systems concerning the effectiveness factor concept.
The aim of this work is to propose a methodology to obtain an effectiveness factor for biofilm in MEC systems and use it to reduce a PDE biofilm MEC model to an ODE MEC model. Therefore, the present work develops a methodology to describe the complex combined effect of local mass transfer and potential on biofilm composed of a growing microbial consortium. It aims to obtain the equation of the effectiveness factor, which in turn allows us to capture the dynamic behavior of the catalytic biofilm. The numerical values obtained from the PDE MEC model are used to determine a bounded range of uncertain parameters that represent the complex combined effect of the biofilm of nonhomogeneous substrate and potential profile space-time dependence. Such parameters can be used in the ODE MEC model. Hence, the proposed methodology looks for an uncertain ODE MEC model in which uncertain parameters also contain biofilm catalytic variability. The resulting model can be used for real-time process monitoring and robust control.

Dynamical PDAE Biofilm Model
According to [26], the biofilm model is based on a set of two PDEs which are coupled to four ODEs. The model (2PDE+4ODE) is depicted in Figure 1. The following assumptions are considered for the two PDEs: (i) the biofilm is a continuum and it is homogeneous; (ii) acetic acid (C 2 H 4 O 2 ) is the only substrate; (iii) competition between methanogenic and exoelectrogenic microorganisms for the carbon source is the only one that occurs, while an inert fraction does not consume substrate; (iv) substrate and potential gradients only take place in the biofilm through the z axis, and thus, axial diffusion, as well as back diffusion, are negligible; (v) substrate is transferred by diffusion only (Fick's law); (vi) the distribution of the substrate in the biofilm rapidly reaches a steady state once the concentration of substrate in the liquid phase changes; (vii) the distribution of the potential in the biofilm rapidly reaches a steady state once the applied potential in the anode changes; and (viii) physical and transport parameters are constant (i.e., density, conductivity, and effective diffusion).
Under these assumptions, the first PDE describes the substrate dynamics in the biofilm as follows [26]: with boundary conditions: where S bio and E bio a are the local substrate concentration and the local potential through the biofilm, respectively, S is the substrate concentration in the liquid phase, ρ x is the biofilm density, D e is the effective diffusion of the substrate in the biofilm, φ m and φ e are the mass fractions of the exoelectrogenic and methanogenic microorganisms, respectively, and φ i is an inert fraction (φ m + φ e + φ i = 1). Notice that the first boundary condition in (1) implies no substrate diffusion across the anode electrode interface. The second boundary condition corresponds to no substrate concentration gradient at the interface between the liquid phase and the biofilm surface.
Again, under the former assumptions, the second PDE describes the potential variation E bio a through the biofilm [26]: with boundary conditions: where E a is the potential at the anode surface, κ bio is the biofilm conductivity that can be seen as the sum of the effects related to the electron transfer mechanisms, F is the Faraday constant, and γ is the conversion time factor. Notice that the first boundary condition in (2) corresponds to no potential losses at the interface between the anode and biofilm surfaces. In contrast, the second boundary condition implies that electrons conduct only on the biofilm matrix. Due to the mass and charge transfer, active microorganisms in the biofilm are exposed to different substrate concentrations and potentials. Therefore, their growth and oxidation rates depend on their positions inside the biofilm. The reaction rates (i.e., specific growth rates and respiration rate) in the biofilm for exoelectrogenic and methanogenic microorganisms are [26]: To find the average reaction rates, the previous expressions are integrated along the biofilm as follows [26]:

Global Balance
The time-dependent mass balance is also based on a previous work [26] and describes an MEC system consisting in an ideal continuous stirred-tank reactor with two chambers.
where S is the acetic acid concentration (as the only carbon source), L f is the biofilm thickness, F is the volumetric flow, V a is the anode chamber volume, S in is the inlet acetate concentration, A a is the anode active area, and k i for i = 1, 2, 3, 4 are yield coefficients.
Under the assumption that all the generated electrons go to the anode, the expected current I (mA) is [26]: where k 5 and k 6 are yield coefficients. Then, the biofilm model consists of two PDEs (1)-(2), local reaction rates (3)-(5), average reaction rates (6)- (8), the global balance (9)- (12), and the expected current (13). Table 3 shows the parameters. For details, the reader is invited to refer to [26].
To solve numerically the previously reported model in (1)-(13), the numerical method ode15s in MATLAB (MathWorks TM , Natick, MA, USA) was used. The discretization of Equations (1) and (2) was firstly proposed and due to their boundary conditions, a double boundary problem was solved with the numerical method bvp5c in MATLAB . Simultaneously, local reaction rates (3)-(5) were obtained. Then, average reaction rates were calculated as indicated in Equations (6)- (8). Finally, the global mass balance (9)- (12) was solved and the current (13) was calculated. The entire procedure was repeated for each time step.

Effectiveness Factor for an MEC Biofilm
The study of reaction-diffusion phenomena in a heterogeneous reaction system commonly involves the solution of PDEs for catalyst particles, where the concept of the effectiveness factor (EF), η, is used as an alternative for avoiding the inherent computationally expensive efforts. The main idea is to express the reaction rate in a catalyst particle as a function of η. For biofilm systems, η has been reported for zero-order and first-order kinetics [41] or Monod- [44] and Michaelis-Menten-type kinetics [40]. For more complex kinetics, there is no analytical solution.
The effectiveness factor has been defined as: the ratio of the actual reaction rate to that which would be observed if the total surface area throughout the catalyst interior were exposed to a fluid at the same conditions as those prevailing at the outside surface of the particle [33,34]. Effectiveness factors can be defined for methanogenic and exoelectrogenic microorganisms and for the respiration process as follows: where S bio is the local substrate concentration (1), S is the bulk substrate concentration (9), E bio a is the local potential (2), and E a is the anode applied potential (operating condition).

Numerical Implementation
To obtain the numerical values of the effectiveness factors (14)-(16), a set of operating conditions for initial substrate concentration (S(t = 0)), inlet substrate concentration (S in ), and anode potential (E a ) were established to simulate the dynamical model in (9)- (12). Notice that the selected operating conditions for S in and E a have been reported in MEC experimental implementations [49][50][51]. Table 4 summarizes the operating conditions used for the numerical implementation. Volume, inlet flow, and active anode area were considered constant. Biofilm thickness initial condition 10 µm Figure 2 shows the solution of the dynamical model in (1)- (12) and the current production for different values of operating conditions (see Table 4). From the comparison of substrate profiles (Figure 2a), it can be said that a more significant consumption of substrate is obtained when higher voltage is applied. In addition, bigger values of biofilm thickness   Table 4: (a) substrate (Equation (9)), (b) biofilm thickness (Equation (10)), (c) mass fractions (Equations (11) and (12)), (d) current (Equation (13)). Arrows indicate the direction of increased applied voltage.
The mass fraction of methanogenic microorganisms is zero at steady state (Figure 2c) in all numerical simulations. For the operating conditions (see Table 4), only exoelectrogenic microorganisms survive. The MEC model in (1)- (12) predicts the noncoexistence of microbial species in biofilm, although an inert fraction remains at steady state. The coexistence of MEC microbial species has been discussed in [24].
The numerical evaluation of the average reaction rates (6)- (8) and reaction rates are shown in Figure 3. It can be seen that the average reaction rate µ m (6) and reaction rate µ m (Figure 3a,b, respectively) exhibit similar dynamical behavior. That is mainly due to the average reaction rate of methanogenic microorganisms being affected by the diffusional resistance offered by the biofilm; thus, µ m < µ m , which is expected. However, this effect is almost negligible under saturated conditions (i.e., high concentrations).
Nevertheless, the average reaction rate µ e (7) and reaction rate µ e exhibit significantly different dynamical behavior. That is mainly due to the average reaction rate of exoelectrogenic microorganisms being affected by diffusional resistance in the biofilm and by the local potential. Notice that µ e presents similar dynamical behavior as µ m in the sense of lower magnitude in reaction rates being developed under lower bulk substrate concentration (S) and low applied voltage (V) as well, which is expected too. The average reaction rate r res (8) and reaction rate r res exhibit different dynamical behavior due to the dependence of the local potential, i.e., r res is only a function of the applied potential, which, for this numerical implementation, is constant (see Table 4).   Table 4. Average reaction rates for (a) exoelectrogenic microorganisms (Equation (4)), (c) methanogenic microorganisms (Equation (3)), and (e) respiration processes (Equation (5)). Reaction rates into the bulk for (b) exoelectrogenic microorganisms (Equation (7)), (d) methanogenic microorganisms (Equation (6)), and (f) respiration processes (Equation (8)). Arrows indicate the direction of increased applied voltage. Figure 4a,c clearly show similar results for effectiveness factors η e and η r . In addition, the behavior of η e and η r are in good agreement with those obtained for µ e and r res (Figure 3c,e), respectively. These results indicate that the dominant phenomenon is conductivity (which is related to the local potential) upon mass transfer. In support information, an additional numerical case is explored when the diffusional resistance offered by biofilm is almost negligible.
It is important to point out that, as mentioned before, µ m and µ m are not a function of the potential. Therefore, the effectiveness factor η m approximates unity to the extent that the saturating substrate operating conditions hold. It is essential to notice that from a practical point of view, it is convenient to operate in such conditions to remove as much substrate as possible, avoiding mass transfer limitations. On the other hand, despite mass transfer limitations not being shown for methanogenic microorganisms, for all numerical simulations of operation conditions (see Table 4), the fraction φ m is zero at steady state. Interestingly, the coexistence of microorganisms does not depend on mass transfer limitations.
From Figure 4 (at steady state), the lower and upper bounds η lb and η ub represent the permissible range for the effectiveness factors (14)-(16) for a set of parameters and operating conditions. In the present example (data in Tables 3 and 4 (15)), (b) methanogenic microorganisms (Equation (14)), and (c) respiration process (Equation (16)). Arrows indicate the direction of increased applied voltage.

Dynamical Reduced Model
The effectiveness factor for an MEC biofilm can be defined as an uncertain but bounded parameterη i ∈ [η lb,i , η ub,i ] with i = e, m, r. The bounded values are η lb,i = min(η * i ) and η ub,i = max(η * i ), where * stands for the value of η i obtained at steady state from the numerical evaluation of different operating conditions (see Table 4). Without loss of generality, a nominal effectiveness factor can be defined as the average of range-bounded values at steady state [52]. Therefore, the effectiveness factor for an MEC biofilm is redefined as a function of the nominal value as follows:η i =η nom,i (1 + |∆ i |/100), where ∆ i is the maximum percentage of variation between the nominal valueη nom,i and the upper or lower value at steady state. The values of η lb,i , η ub,i ,η ub,i , and |∆ i | are obtained from Figure 4 at steady state and are shown in Table 5. Then, a reduction of the MEC model in (1)- (12) is then defined by the following set of ODEs: Consequently, the dynamic MEC model in (1)-(12) with the biofilm model described by the set of PDEs (1)-(2) and the set of ODEs (9)-(12) is reduced through an effectiveness factor for an MEC biofilm. The reduced MEC model, in a compact form, is defined as follows: It consists of the set of ODEs (17)- (20); state variables Φ = [S, L f , φ e , φ m ] T ; the set of parameters Π = [µ max,e , µ max,m , K S,e , K S,m , K i , Ek a , ρ X , b det , b res , b in , f exp , D e , k bio , k 1 , k 2 , k 3 , k 4 , k 5 , k 6 , e cat , k H2 ]; and operational conditions [S in , E a , T, V a , F, A a ] with uncertain but bounded parametersη e ,η m , andη r , whereη e =η nom,e (1 + |∆ e |/100), η m =η nom,m (1 + |∆ m |/100), andη r =η nom,r (1 + |∆ r |/100).
The main assumptions of the reduced ODE MEC model are the following: (i) acetate is the only substrate in the feed wastewater; (ii) the anodic chamber operates as an ideal continuous stirred-tank reactor; (iii) all variables are considered spatially uniform; (iv) exoelectrogenic and (possibly) methanogenic microbial populations are mostly attached to the anodic biological biofilm; (v) in consequence, biomass growth in the anodic bulk phase is negligible; (vi) microbial populations compete for the same substrate; (vii) there is instant gas transfer from the liquid to the gas phase; and (viii) pH = 5.5 and temperature T = 25 • C are constant. Figure 5 shows the schematic representation of the reduced ODE MEC model in (21).
It is important to remark that even when suspended bacteria may affect biofilm growth, mass transfer, and electricity generation in BES [53,54], the former assumptions (physically plausible) imply the inhibition of methanogenic archaea growth (mainly because of the acid pH) and favor that the only exoelectrogenic bacteria that survive are eventually attached to the anode [55].
Experimental results in the literature reported that there are lower efficiency and electricity generation when scaling up BES [3,6]. The reactor and electrode size influence the efficiency of the MEC system and therefore the electricity generation [3,6]. Despite ideal suppositions on the reduced ODE MEC model, the proposed EF approach could be applied to a more complex PDE MEC system. In this sense, a detailed PDE MEC model for current, potential, and nonideal flow patterns can be solved numerically for complex geometries and then applied to the proposed EF approach on a representative ODE MEC model. Thus, reactor-electrode size side effects on mass transfer and electricity generation could be included in the alternative reduced ODE MEC model in a scale-up approach.

Equilibrium Points
In order to recognize the possible steady-state solutions for the model in (21), it is necessary to compute the equilibrium points and then classify them in the sense of stability criteria. Then, the following steps are performed: (i) to compute an equilibrium point under parameters and operating conditions; (ii) to find the Jacobian matrix of the nonlinear reduced model in (21); (iii) to compute the eigenvalues of the Jacobian matrix at an equilibrium point; and (iv) to corroborate the criteria of stability. The equilibrium point is locally asymptotically stable if all eigenvalues have negative real parts. Otherwise, if there is even one eigenvalue with a positive real part, the equilibrium point is unstable. As an example of this procedure, Table 6 shows the numerical evaluation of the equilibrium points under the set of parameters in Table 3 and operating conditions S(t = 0) = 2 [g L −1 ], S in = 2 [g L −1 ], and E a = 0.3 [V], and the remaining operating conditions are shown in Table 4. Table 6. Equilibrium points ψ of reduced model in (21). Notice that only the equilibrium point ψ 01 exhibits stable behavior and has a physical meaning. It is said that the MEC system is operating in a desirable condition of noncoexistence of microorganisms in which only exoelectrogenic microorganisms degrade the substrate [24]. Equilibrium point ψ 02 exhibits stable behavior, but this point does not have physical meaning because this would imply a physically impossible negative value for the substrate at steady state. Equilibrium point ψ 03 is particularly interesting because it shows that the MEC system exhibits washout and that the steady-state substrate concentration is given by its inlet composition. However, ψ 03 will not arise because it is an unstable equilibrium point. The analysis is similar for the other equilibrium points ψ 04 to ψ 13 , from which it can be deduced that the reduced MEC model in (21) has only one reachable locally stable equilibrium point ψ 01 under the set of parameters Π (Table 3) and the operating conditions under consideration for this example.

Parametric Sensitivity Analysis
Departing from the reduced MEC model in (21), a sensitivity parametric analysis [52] is performed to determine which parameters in Π affect the dynamical behavior of the field f (Φ; π). For a set of parameters π, the approximate solution for the sensitivity function is computed by the simultaneous solution of the MEC model in (21) and the linear time-varying sensitivity equationΘ f as follows [52]: The parameters of Table 3 and operation conditions in Table 4 are used to numerically solve Equation (22). In addition, in order to appreciate the parameter subset having the most significant effect of the MEC model in (21), the following parameter values are changed from Table 4 Table 5. The following initial conditions for the numerical solution of (22) are considered: Φ(t 0 ) = [1 g L −1 , 10 µ m, 0.3, 0.3] T and Θ 0 = [1, 1, 1, 1] T . The numerical method ode15s in MATLAB was used to solve numerically Equation (22). Figure 6 shows the numerical solution of (22). Notice that at stable steady-state conditions, the parameter subset having the most significant effect of the MEC model in (21) is: K S,e , b in , µ max,e , andη e . Interestingly, only the effectiveness factorη e arises as the important parameter in the parametric sensitivity analysis in comparison toη m andη r .

Open-Loop Dynamical Behavior
The performance of the reduced MEC model in (21) is evaluated and compared in an open-loop simulation with the following ODE model reported by Pinto et al. [30]: where S is the substrate (acetate) concentration; X e and X m are the exoelectrogenic and methanogenic microorganism concentrations, respectively; M ox is the oxidized mediator fraction per exoelectrogenic microorganism; q e and q m are the acetate consumption rates by exoelectrogenic and methanogenic microorganisms, respectively; µ e and µ m are the growth rates; α is the dimensionless biofilm retention constant; D is the dilution rate (D = F in V −1 reac ); Y M is the oxidized mediator yield; γ is the mediator molar mass; m is the number of electrons transferred per mol of mediator; and F is the Faraday constant. The following expressions are included: if X e + X m > X max 0 otherwise (24) where M total is the total mediator fraction per microorganism; M red is the reduced mediator fraction per exoelectrogenic microorganism; q max,e and q max,m are the maximum acetate consumption rates; µ max,e and µ max,m are the maximum growth rates; K S,e and K S,m are the half saturation constants; and X max is the maximum attainable biomass concentration. Table 7 summarizes the set of parameters used in the model in (23) and in the auxiliary functions (24). Notice that the biofilm in the model in (23) is limited by X max . Therefore, for the sake of normalization, the biomass fraction for the model in (23) is defined as follow: The operating conditions for both models are: S in = 2 g L −1 , E a = 0.5 V, V a = 50 mL, and F = 50 mL. Additional operating conditions used in the reduced model in (21) Figure 7 illustrates the dynamical behavior of the state variables of the reduced model in (21) and the model in (23). In the same figure, it can be seen that S in both models converges to a near-steady state. Notice that the values for S predicted by the model in (23) are smaller than the ones predicted by the reduced model in (21) for 0 < t < 20 [d] (see Figure 7a). Moreover, the model in (23) predicts a minimum substrate concentration for 0 < t < 5 [d]. This is because there is an increase in the mass fractions of methanogenic ( Figure 7c) and exoelectrogenic (Figure 7d) microorganisms. In addition, for t > 5 [d], the model in (23) predicts: (i) an increase in S (Figure 7a) and the mass fraction of exoelectrogenic microorganisms ( Figure 7d) and (ii) a decrease in the mass fraction of methanogenic microorganisms. On the other hand, all the states of the reduced model in (21) exhibit smooth dynamic behavior. Moreover, states S and Φ m decrease, and states L f and Φ e slowly increase and finally reach a steady state. It should be noted that the dynamic smoothness of the reduced model in (21) may be desirable in control or real-time applications.
It is important to remark that the aim of the work was to propose a methodology to obtain an effectiveness factor for biofilm in an MEC system and use it in the reduction of the PDE biofilm MEC model to the ODE MEC model, but the experimental validation stricto sensu is not provided here. Indeed, the numerical implementation includes operating conditions from the literature [49][50][51]. Moreover, the reduced ODE model was compared with the validated model taken also from [30] and the results were in good agreement. Nevertheless, the parameters of the reduced model should be adjusted to be successfully applied to lab-scale or pilot-scale implementations. This task is planned as future work.

Conclusions
In this work, an effectiveness factor for the microbial electrolysis cell model was proposed. This effectiveness factor was defined for the biofilm respiration process and methanogenic and exoelectrogenic microorganisms. Departing from a partial differential equation and numerical data obtained from a set of operating conditions, a set of uncertain but bounded ranges of effectiveness factor values was obtained. The stability and parametric sensibility of the reduced model were evaluated. Then, a reduced model based on ordinary differential equations was obtained. The reduced model was analyzed in the sense of stability and parametric sensibility. The reduced model was also tested in numerical simulation using a validated model taken from the literature.
The proposed effectiveness factor was included in a microbial electrolysis cell model described by a set of ordinary differential equations. The procedure is an alternative to address the tradeoff between (i) the local description of multipopulation microorganisms, local mass, and local potential in the biofilm and (ii) a simple and reliable alternative description of biofilm in an ordinary differential equation representing a continuous microbial electrolysis cell system.
The procedure could be extended to models comprising mass and charge transport in biofilm, simulated in two or three dimensions and then coupled with global mass balances for microbial electrolysis cell systems.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: Partial Differential Equation