A Review of Control-Oriented Bioelectrochemical Mathematical Models of Microbial Fuel Cells

: A microbial fuel cell (MFC) is a potentially viable renewable energy option which promises effective and commercial harvesting of electrical power by bacterial movement and at the same time also treats wastewater. Microbial fuel cells are complicated devices and therefore research in this ﬁeld needs interdisciplinary knowledge and involves diverse areas such as biological, chemical, electrical, etc. In recent decades, rapid strides have taken place in fuel cell research and this technology has become more efﬁcient. For effective usage, such devices need advanced control techniques for maintaining a balance between substrate supply, mass, charge, and external load. Most of the research work in this area focuses on experimental work and have been described from the design perspective. Recently, the development in mathematical modeling of such cells has taken place which has provided a few mathematical models. Mathematical modeling provides a better understanding of the operations and the dynamics of MFCs, which will help to develop control and optimization strategies. Control-oriented bio-electrochemical models with mass and charge balance of MFCs facilitate the development of advanced nonlinear controllers. This work reviews the different mathematical models of such cells available in the literature and then presents suitable parametrization to develop control-oriented bio-electrochemical models of three different types of cells with their uncertain parameters.


Introduction
Electricity demands have seen a rapid increase due to industrial growth and rise in population. Fossil fuels (coal, natural gases), nuclear, and renewable sources are the avenues for providing this need [1,2]. Electricity generated from fossil fuels negatively impacts the environment by polluting the atmosphere and resulting in global warming. Continued innovation in renewable energy sources can over time result in a drop in fossil fuel-based energy dependence [3]. Of late, a significant proportion of electricity is generated from renewable energy sources like wind, geothermal, solar, biomass [4,5] fuel cell, etc. as a result of a lot of research work carried out in this sector. The power generated from such sources is clean (do not produce pollutant gases) and are more efficient [6]. A hydrogen-based fuel cell is also a promising renewable energy source which transforms chemical energy into electrical energy. However, conventional hydrogen-based fuel cells have some disadvantages like high production cost and mass generation, hydrogen storage, and transportation [7].
Wastewater containing fermentative, anodophilic, and methanogenic bacterias are drastically growing because of industrial growth and accelerated human population, thereby necessitating efficient wastewater treatment solutions. However, most of such technologies are inefficient, expensive, and not sustainable [8]. Microbial fuel cells (MFCs) provide a new technology that can act as pollutant removal devices by using microorganisms available in wastewater as catalysts to oxidize substrates and produce much needed electric power. MFC technologies have unique features and abilities to solve energy and environmental issues in remote, rural, and dry areas. MFCs are used as power supplies for underwater monitoring devices, remote biosensors, biological oxygen demand monitoring, and bio-hydrogen production along with wastewater treatment applications. The major advantages of MFCs include low operational temperature, less excess activated sludge, higher conversion efficiency, no requirement of gas treatment, direct electrical energy conversion from the substrate, and easy implementation in areas with limited electrical facilities [9][10][11].
In the anode chamber, a mixture of microorganisms culture and suitable substrates are placed, which generates protons and electrons by the oxidation process. Substrates are chosen based on types of microorganisms culture and performance requirements [12]. Microorganisms have an inherent capacity to transport developed electrons on the anode surface in two ways: direct or indirect electron transfer. In the direct transfer process, electrons are transported via nanowires or intracellular mediator, whereas additional external mediator is used in an indirect way [13]. Most of the MFC setups have two electrodes, namely anode and cathode, nourishment for bacterias (substrates), membranes, and microorganisms culture. The fundamental organization of the different sections in an MFC is shown in Figure 1. Now, electrons at the anode surface are transported to the cathode via an external electrical load [14], and the protons are transported to the cathode through cation exchange membrane (CEM). Finally, generated electricity is captured across the electrical load between electrodes and also filtered water is produced at cathode by a combination of electrons and protons that came from the anode compartment.
MFC setups are mainly classified by three major aspects such as design, operation, and biological aspects. Each aspect has some subdivisions to understand the functionality and detailed behavior of MFC. Several review articles have focused on performance improvement, electrode and membrane material advancement, electron transfer mechanism, and design and parameter effects [15][16][17][18][19][20][21][22][23][24][25]. Recent advancement in bacterio-algal MFCs with operational parameter and their effect on power density and its efficiency have been investigated [26]. Another recent review work has briefly overviewed MFC technology for wastewater treatment, the possible future applications, and primary challenges [27]. However, for the most part, these comprehensive reviews have mostly focused on the development undertaken on MFCs, their challenges, and future applications. MFCs are a challenging technology still confined to the laboratory level and awaiting significant commercialization. Model-based design and control approaches would be helpful for the commercialization of MFCs technologies under different operating conditions. The mathematical models relate between various parameters, system inputs and outputs, and also system characteristics. Few bio-electrochemical MFC models are formulated and validated with certain assumptions and circumstances. There are some review articles available that have focused on the mathematical modeling of MFCs [28][29][30][31][32]. However, these reviews do not delve into control-oriented models which parametrize the uncertain parameters of MFCs.
This work provides a review of the different MFCs and control-oriented models of such MFCs. In this work, a literature review of mathematical modeling of different MFCs is provided along with their controllable models with identified uncertain parameters. It is the first attempt to present an extensive description of the three models of MFCs such as single species anode based MFC, two species anode based MFC, and single species anode-cathode based MFC with uncertain parameters. The paper is organized as follows: Modeling strategies and description of mathematical models of MFCs are summarized in Section 2. The control-oriented models with their uncertain parameters are summarized in Section 3. Modeling of MFCs based on its applications is summarized in Section 4. Developments in control strategies of MFCs are summarized in Section 5. Section 6 presents the perspectives and challenges of MFCs, followed by the concluding section.

Modeling Strategies and Description of Mathematical Models
Mathematical model of any system transforms the compound system operations into facile mathematics representing dynamical equations. A mathematical model of the MFC system incorporates the impact of parameters such as design (i.e., membrane and electrode materials, type and volume of chambers, etc.), operational (i.e., number of populations and source of microorganisms, electron transfer mechanisms, etc.), and biological parameters (temperature, pH, resistances, flow rate, etc.) on the system output [33,34]. There are several ways to develop a mathematical model of MFC by considering different factors. For instance, based on the culture of the microorganism accounted for in the model development, either it can be a pure culture with one or two microorganisms or mixed culture with multiple types of microorganisms or based on the way of transportation of extracellular electron considered in the model, either it can be direct transfer or indirect transfer, or whether biofilm is considered or not [35,36]. Apart from these factors, mathematical modeling of the MFCs is done in either engineering or statistical approaches as shown in Figure 2.
Furthermore, engineering-based models can be done in two ways, either ordinary differential equations (ODE) or partial differential equations (PDE). ODE based models are accounting for only time dependence and no spatial dimension. ODE models have simple formulation and computation that is quite faster at a relatively low computational cost. PDE based models are accounting for both time and spatial dimensions, which makes a more informative and comprehensive system. However, the computational time and cost both are quite higher than the ODE models. One can choose as per the requirements based on performance and applications. Table 1 provides a brief overview of the important factors accounted for in some of the familiar MFC models. Zhang and Halme proposed one of the first simplest mathematical models of dual-chamber MFC, which is based on ODE. In this model, Monod kinetics are used to represent the substrate consumption by the microorganisms and redox reactions. Nernst's equation and Faraday's law are used for formulating the dynamics of emf and electrochemical current, respectively. The proposed model is proposed by taking several assumptions but provides the best initiative model for MFC study and used as a benchmark for the advanced models. This model is validated through simulation and experimental setups with different electrical load conditions [37].  Picioreanu et al. proposed a dual-chamber MFC model with batch mode configuration. In this model, acetate is used as a substrate and Geobacter sulfurreducens as a microorganism culture. It is quite a simple mathematical model based on ODE. It represents the accurate dynamics of anode chamber with suspended cells and electron transportation through a diffusible mediator. The dynamics of current density are presented using the Butler-Volmer equation and also the effect of different resistances including internal and external resistances is studied to analyze the output of MFC [48]. This model is more comprehensive compared to Zhang and Halme's model, but they have considered only one simple substrate (acetate) for experimental purposes.
Zeng et al. proposed an MFC model that accounted for the phenomena in an anode and cathode chamber with a pure culture of microorganism. The proposed model is an approximately similar model of a conventional chemical fuel cell. The Butler-Volmer equation is used to represent the electrochemical reaction happening at anode and cathode chamber. In addition, Monod kinetics equation is used to represent the mass balances during MFC operation. In this model, positive ions are not considered in the cathode reaction. The total MFC voltage is determined after considering the voltage losses due to over-potentials and resistances. The model is validated using two distinct substrate feed such as acetate and glucose-glutamic acid. Glucose-glutamic acid is used as a solution to make artificial wastewater. The only limitation of this model is that it does not account for the biofilm phenomena at the anode that has appeared to remarkably affect the overall performance of MFC. However, this model is very simple for performing experiments and also assisting as a good benchmark for the development of dual-chamber MFC models [41].
Pinto et al. have developed an ODE based dynamic model of single-chamber MFC. It only illustrates the dynamics of the anodic chamber considering two microorganisms' populations (anodophilic and methanogenic bacteria). Monod kinetics are used to represent the dynamics of mass balances of the substrate, both microorganisms' concentrations and intracellular mediator. Nernst's equation is used to represent the electrochemical dynamics. The electrodes' potentials and electrochemical current are calculated using Ohm's law. The parameters of the model are estimated by applying the Neder-Mead simplex algorithm. This model is comparatively easy to implement and enables fast numerical computation [39]. Shankar et al. proposed mathematical modeling of continuous mode, dual-chamber MFC. In this model, glucose glutamic acid is used as a substrate and the study is carried out to analyze the effect of temperature variations on the current density in both the chambers and also heat transfer through the membrane [44].
Oliveira et al. proposed an extensive dual-chamber MFC model by incorporating heat balance and biofilm formation. A combination of Monod kinetics and Tafel's equation is used to model the anode chamber's dynamics. However, only Tafel's equation is used to model the cathode chamber's dynamics. Fick model is utilized to illustrate the mass transfer mechanisms in the electrodes and biofilm, and also Fourier's law is used to describe heat balance. The simulation and experimental studies on this model indicate the impact of current density on the over-potentials, biofilm, and profiles of different concentrations and temperature. This model is considered as an extended version of Zeng's model [42]. This study doesn't include dynamic analysis of this model. Thus, it can be challenging to find the influence of operational parameters on MFCs overall performances. Recio-Garrido et al. proposed a combined bioelectrochemical-electrical MFC model with a single microorganism culture. This model enables process optimization as well as real-time estimation of physical and electrical parameters [40].
Esfandyari et al. developed a mathematical model for batch mode as well as continuous mode dual-chamber MFC with a pure culture of Shewanella [43,50]. In this study, three different kinetics models such as Monod kinetics, Blackman model, and Tessier model are investigated to calculate the microorganisms growth rate. Activation, ohmic and concentration losses are accounted for to calculate the final MFC output voltage. This study shows that the Monod kinetics model is far better and more accurate than the other two kinetic models. The process of MFC modeling and optimization of parameters can be found in [51,52] The major limitation of this model is uniformly distributed quantities, which is considered one of the influencing factors of biofilm's dynamic performance.

Control-Oriented Mathematical Models
MFCs performances are significantly dependent on the design, operational, and biological parameters, which are briefly discussed in Section 2. Therefore, MFC operation should be carried out under a controlled environment to get a steady output with the best possible efficiency. Thus, advanced control strategies are required for the controlled operation of MFCs. Consequently, a control-oriented mathematical model is needed to develop advanced control strategies. In this section, control-oriented mathematical models of three different MFCs with their uncertain parameters are proposed, which can be useful for developing advanced control strategies.

Single Chamber MFC Model with a Single Microorganism
Abul et al. proposed a single-chamber MFC model with pure microorganism culture. The principal operation of MFC can be illuminated by Figure 3. It illustrates the chemical reaction process and also shows the types of substrate and microorganism used in the modeling process. In this model, some reasonable assumption has been taken [38]. The overall chemical reaction during the operation of MFC is given by where CH 3 COO − is the substrate. HCO − 3 is the byproduct of this operation. In this modeling process, Monod kinetics are used to represent the growth rate of microorganisms, which is given in Table 2.
The substrate utilization rate and microorganisms growth rate are associated with microorganism yield, Y as where q max is the maximum value of substrate utilization rate. The comprehensive modeling of single-chamber MFC with a single microorganism is presented in [38]. Table 2. Equations in physical and bio-electrochemical models of Microbial Fuel Cells [28].

Monod Kinetics
It illustrates microorganisms grwoth and substrate utilization rate. µ = µ max Cs Ks +Cs µ = specific growth rate, µ max = maximum value of µ, C s = substrate concentration, K s = Monod constant.

Nernst
It describes electrode potential based on an electrochemical reaction. The control-oriented mathematical modeling of single-chamber MFC with a pure culture of microorganism is given byẋ where x 1 , x 2 , and u are the state variables, and input of the system, respectively. X 1 and X 2 stand for substrate and microorganism's biomass concentrations. Dilution rate (ratio of feed flow rate to chamber volume) is considered as the input, u. S 0 and K d refer to the initial substrate concentration and microorganisms decay constant, respectively. The uncertain parameter is denoted by γ −1 , which is maximum microorganism growth rate, µ max . It significantly relies on the types of substrate and microorganisms used during the operation of MFC. The typical values and initial conditions of the parameters can be found in [53,54].

Single Chamber MFC Model with Two Microorganisms
Pinto et al. proposed a single-chamber MFC model with two microorganisms namely anodophilic and methanogenic. Several reasonable assumptions have been taken while modeling of MFC. The principal operation of MFC can be illuminated by Figure 1. The chemical reactions happening in the anode chamber are given by In this model, Monod kinetics is used to represent the growth rate of both the microorganisms. Moreover, the anodophilic microorganism's growth rate significantly relies on the substrate concentration and mediator, whereas the growth rate of methanogenic microorganism is only dependent on the substrate concentration. The comprehensive modeling of single-chamber MFC with two microorganisms is presented in [39].
The control-oriented mathematical model of single chamber MFC with two microorganisms is given by where x i , i = 1, 2 . . . 4 are the state variables, namely substrate concentration, anodophilic microorganism concentration, methanogenic microorganism concentration, and oxidized mediator fraction, respectively. u is the input, which is the dilution rate. Subscripts a and m refer to the anodophilic and methanogenic microorganisms. K M and V are the mediator constant and chamber volume, respectively. MFC current is denoted by I MFC . k 1 and k 2 are the inverse of the microorganisms yields. Other parameters' names can be found in Table 2. The uncertain parameters, θ * 1 and θ * 2 refer to the maximum growth rate of anodophilic and methanogenic microorganisms, respectively. The typical values of parameters and bounds of uncertain parameters can be found in [53].

Dual Chamber MFC Model with Single Microorganism
Zeng et al. proposed a dual-chamber MFC model with a pure culture of microorganism [41]. The principal operation of dual-chambered MFC can be illuminated by Figure 4. The overall chemical reactions in both anode and cathode chambers are given by The control-oriented dynamics of anode chamber operation are given by where x i , i = 1, 2, . . . , 5 are state variables. x 1 is the substrate concentration, S. Other state variables are defined by C CO 2 , C H , and X are the concentration of carbon dioxide, hydrogen, and biomass, respectively. Superscript in refers to influent. The dilution rate is considered as system input u. Other parameters can be found in [53]. The uncertain parameter is defined by δ * 1 = β F R T . It significantly relies on the operational temperature.
The control-oriented dynamics of cathode chamber operation are given bẏ where x i , i = 6, 7, . . . , 9 are state variables. x 6 is the oxygen concentration, C O 2 . Other state variables are defined by where β 2 is a constant defined as 4F V c A m . Other parameters can be found in [53]. The uncertain parameter is defined by δ * 2 = (α−1)F R T . It largely relies on the operational temperature and input feed flow rate. These uncertain parameters can be estimated online through adaptive control mechanisms. The above control-oriented models have considered uncertainties in parameters, although one can also consider another uncertainties (i.e., measured, combined, and expanded uncertainties) while modeling MFCs [55].
MFC Models presented in Sections 3.1 and 3.2 are based on the dynamics of microorganisms growth and decay rates. These models come up with the influence of substrates on the microorganisms' growth rates. Therefore, control engineers can develop controllers to control the microorganisms' growth rate by varying the substrate concentration as per the desired MFC performances. In this modeling process, the effect of internal and external temperature is neglected, although it is an important factor for bacterial growth rate (certain bacteria can remain alive within a specific range of temperature). Maximum efficiency of MFCs can be achieved by maintaining appropriate internal and external temperature during the whole operation. Moreover, one can choose the MFC models based on several criteria such as the number of microorganisms' populations (i.e., one or two populations), controllers' simplicity, individual applications, and economical aspects. However, the models discussed in Section 3.3 are based on the dynamics of chemical reactions happening in both chambers of MFCs. The temperature of the operation is considered during the modeling process. The dynamics of this model are similar to conventional hydrogen-based fuel cells. There are several well-known advanced controllers already available for conventional fuel cells. Therefore, it could be easy to develop the same kind of controllers for MFCs.

Modeling of MFCs Based on Applications
The key aim of MFC operation is to obtain maximum power density just as other energy-providing devices. Models discussed in previous sections are based on the basic principle of operation, which have longer computational time and complex microorganisms operations, and also must require interdisciplinary knowledge for practical applications. Thus, application-based MFC models are developed to overcome the aforementioned hurdles to commercialize MFCs for specific predefined applications. There are three main approaches to develop application based models such as equivalent electrical circuits, based on intelligence and a system identification approach.

Equivalent Electrical Circuit Based MFC Models
In this approach, performance characteristics and dynamical behavior of MFC are described by electrical parameters such as voltage, resistance, current, and capacitance. Thus, the voltage and current of the electrical circuit model are referred to as the biological and electrochemical reactions in MFCs. A simple equivalent electrical circuit of MFC shown in Figure 5 is proposed by [56]. In this model, R 1 and R 2 are the internal resistances. C and R load are the internal capacitance and external load of MFC, respectively.  [57]. Research studies in the literature showed that the output power of MFCs is significantly influenced by the internal resistances [58,59]. He and Mansfeld have reported a tool, namely electrochemical impedance spectroscopy, to calculate the internal resistance of MFC and attempted to show the effect of it on MFC's overall performances [60]. The parameters of the equivalent electrical circuits such as resistances, capacitance, and phase elements are determined by several approaches (i.e., polarization curve tests) [61][62][63]. Moreover, the types of electrodes are also one of the key factors to determine the internal resistances [64]. In [65], they have proposed two equivalent electrical models for a two-chamber MFC by using different material of electrodes and studied the effects of electrodes' material. The effect of the biological capacitance of anodes in MFCs has been studied in [66]. Various equivalent electrical models of MFCs have proposed based on the applications in the literature. Furthermore, these models of MFCs may be used to develop advanced control strategies and MPPT algorithms to extract maximum power from MFCs.

Intelligence Based MFC Models
Intelligence-based models did not require in-depth knowledge of MFC and its operation. It requires only a huge number of experimental data to obtain the model through a learning approach. The accuracy of learning models significantly relies on the number of data and its sampling time. MFCs have slow dynamical characteristics, which takes a long time to get the experimental data. The modeling process using learning or intelligence approach is shown in Figure 6. Machine learning models like Artificial neural networks (ANN), Support vector machines (SVM), and Extreme learning machines (ELM) are effective in the forecasting of time-series data [67,68], and may be applicable in the operation of multi-variable MFCs after they are constructed. The Relevance Vector Machine (RVM), a Bayesian model for regression of identical functional form to the support vector machine, is proposed by Fang et al. to get MFC's mathematical model for optimizing the output power with multiple variables [69]. Furthermore, ANN is used to analyze the performance characteristics of the MFC model with four manipulated input variables. Moreover, the advanced modeling approach based on ANFIS is acquired to analyze the performance of MFC [70]. Garg et al. have compared the performance of the MFC model obtained from three different approaches such as ANN, MGGP, and SVR, and concluded that the MGGP is far better than the other two approaches, and also stated that ANN is better than the SVM approach [71].

System Identification Based MFC Models
The system identification approach is ideally used when the inputs and outputs of the system are available through experiments. The model identified shows the association between nonlinear characteristics and time function of input and output. The fundamental of modeling process using system identification approach is shown in Figure 7. One can obtain different types of models such as transfer function, state-space, the process model, nonlinear ARX model, and polynomial model. Boghani et al. have proposed an MFC model using a system identification technique. In this model, the electrical load of MFC is considered as a manipulated input variable and voltage of MFC as an output variable [72]. Later, Boghani et al. have developed a gain scheduling controller for the aforementioned MFC model to maintain the output voltage constant [73]. Ravi and Deb have developed dual-chamber MFC transfer function models for both anode and cathode chamber dynamics using a system identification method [74].

Developments in Control Strategies of MFCs
MFCs have several challenges for controller development such as gains of controllers, nonlinearity in the system, etc.. Selection of manipulated input and output variables can be a critical factor for the development of MFC's controller. The output power of MFC is significantly varied with the change in influent substrate flow rate and initial substrate concentration. Thus, one can choose one of these as a manipulated input variable. It can be chosen based on the type of control strategy and application of MFC [75]. Yan and Fan proposed a fuzzy proportional integral derivative (PID) controller for dual-chamber MFC to maintain a constant output voltage [76]. Later, Fan et al. developed an adaptive fuzzy control for MFC to maintain a constant voltage under various load conditions, and this adaptive controller provides superior steady-state performance and faster response compared to fuzzy control [77]. The model predictive control strategy is applied to dual-chamber MFC by Fan et al. [78]. The robustness of the controller is presented by considering disturbances in operational temperature and substrate concentration and concluded that the proposed controller provides fast response under disturbances.
Patel and Deb proposed a novel adaptive backstepping control strategy for a single chamber MFC with pure culture. In this study, an uncertain parameter (microorganisms growth rate) is estimated online through adaptive control law, and the influent substrate flow rate is considered as a manipulated input variable [53,54]. Later, Patel and Deb developed an adaptive control technique for single chamber MFC with two microorganisms. In this control technique, the initial substrate concentration is considered as a manipulated input variable and two uncertain parameters are estimated [79]. An LMI based robust control strategy with norm bounded uncertainty for single chamber MFC is proposed in [80]. In this study, a nonlinear MFC model is linearized at its equilibrium points and dilution rate is considered as a manipulated input variable. An exact linearization of dual-chamber MFC is proposed in [81]. Conventional Model Reference Adaptive Control(MRAC) and MRAC with MIT rule (a rule developed by Massachusetts Institute of Technology for MRAC) strategies for dual-chamber MFC have been proposed in [74]. In this study, the MFC model is obtained from system identification toolbox in MATLAB. Luo et al. proposed an MRAC scheme for MFC to track the output voltage [82].

Perspectives and Challenges
Recycled wastewater is widely used in dry areas to meet water demand. However, wastewater treatment plants require appropriate water (wastewater with a particular mixture of microorganisms, which can be easily and economically treated) and more energy for treatment. MFCs have potential as renewable energy sources and may be capable of resolving energy and environmental issues. MFCs are more certainly sustainable compared to conventional wastewater treatment technologies. MFC technologies and applications are still at the laboratory level but are great potential industrial wastewater treatment applications and renewable energy sources in remote areas. MFCs have some limiting attributes such as lower power density, materials cost and efficiency. Therefore, more research has been conducted on the electrode, catalyst and membrane materials, types of substrate, transport mechanisms of electrons, etc.. The voltage generated from a single MFC is very low for electrical appliances. Several MFCs can be stacked together to obtain the desired voltage level. Stack MFCs require individual control systems as well as combined control actions for controlled voltage and flexible operation. There are certain challenges of mathematical modeling of stack MFCs and physical connections of electrodes. Further improvement should pay a lot of attention to the modeling of stack MFCs and selection of the materials therein.
Apart from these challenges, some researchers have attempted to scale up MFCs as a wastewater treatment plant. In [83], authors have analyzed MFCs performance by treating 200 liters of municipal wastewater. They have combined 96 MFC modules in series-parallel connections and continuously operated for around 300 days. The results of this study concluded that the performance of MFCs is not as good as the laboratory level in terms of energy recovery. Therefore, more intense research on MFCs could help to fill the gap between theoretical and real-world applications to recover more energy from MFCs. The data (i.e., energy consumption and generation, types of wastewater, etc.) of current operating wastewater treatments plants and small scale MFCs are available [84,85].
MFC technologies require significant knowledge of multidisciplinary fields such as microbiology, material science, chemical reactions, and control engineering. The performance of MFCs is enhanced by providing appropriate control and optimization techniques. These techniques help maintain strict reaction conditions such as pH value, operational temperature, substrate and biomass concentration, and rate of oxygen supply, and also guarantee the maximization of performance. Advanced control schemes such as adaptive control, model predictive control, sliding mode control, etc. are applicable in optimizing the operational parameters of MFCs described in previous sections. Standard mathematical models of MFC are not well documented for all conditions and applications until now. Some attempt has been made to develop models while considering different aspects like biofilm, chemical reactions, activity in microorganisms and substrate oxidization, and electrons transfer. It is necessary to make reasonable modification in the mathematical models for effective control. Control-oriented mathematical models of MFC are necessary to develop suitable advanced control techniques. Certain parameters like growth and consumption rates are sensitive towards and affect the performance of MFCs, and finding appropriate bounds of these parameters are quite challenging and a lot of experimental attention is required. Hence, control-oriented models can contribute to scale-up commercial MFCs with suitable control actions.

Conclusions
Different mathematical models with appropriate control-oriented forms for Microbial fuel cells are described for an attractive renewable technology which provides electricity and manages wastewater simultaneously. Such devices have found potential applications such as bio-sensors, wastewater treatment plant, underwater monitoring, and remote sensing. Model-based control and optimization techniques developed by identifying and estimating different parameters which vary within practical known bounds are required to evaluate the behavior and performance of such diverse systems under undesirable conditions. Development in different models has been reviewed and analyzed in terms of the microorganisms culture, compartment modeled, modeling perspective, etc.. This paper has reviewed three different control-oriented bio-electrochemical models of MFCs and provides control-oriented mathematical models with uncertain parameters. According to the discussions in previous sections, the following recommendations can be inferred:

•
Mathematical model of MFCs, and its control and optimization strategies are chosen based on the specific applications and operational requirements.

•
Bounds of uncertain parameters should be practical with a higher confidence level. It is advisable to do more experiments to get accurate practical bounds. • Several reasonable assumptions may be required in mathematical modeling while developing advanced control and optimization strategies to boost the overall performance of MFCs. • Selection of appropriate manipulated input variables for successful development of efficient control actions is required. • Development of control and optimization strategies for MFCs should be economical, environment friendly, and reliable.
Author Contributions: Writing-original draft, R.P.; Writing second draft and Editing, R.P. and D.D.; Review, D.D., V.E.B., and R.P. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

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

Abbreviations
The following abbreviations are used in this manuscript: