<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Chem. Eng.</journal-id>
<journal-title>Frontiers in Chemical Engineering</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Chem. Eng.</abbrev-journal-title>
<issn pub-type="epub">2673-2718</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1235125</article-id>
<article-id pub-id-type="doi">10.3389/fceng.2023.1235125</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Chemical Engineering</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Circular economy-based multi-objective decentralized controller for activated sludge wastewater treatment plant</article-title>
<alt-title alt-title-type="left-running-head">Subbian et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fceng.2023.1235125">10.3389/fceng.2023.1235125</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Subbian</surname>
<given-names>Sutha</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Natarajan</surname>
<given-names>Pappa</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Murugan</surname>
<given-names>Chitra</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2248067/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Bioprocess Laboratory</institution>, <institution>Department of Instrumentation Engineering</institution>, <institution>Anna University</institution>, <addr-line>Chennai</addr-line>, <country>India</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Rajalakshmi Institute of Technology</institution>, <addr-line>Chennai</addr-line>, <country>India</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/586331/overview">Santhana Krishnan</ext-link>, Prince of Songkla University, Thailand</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2364283/overview">Mullai Pandian</ext-link>, Annamalai University, India</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2371966/overview">Abirami Balakrishnan</ext-link>, Hindustan University, India</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Chitra Murugan, <email>chitramurugan.g@gmail.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>5</volume>
<elocation-id>1235125</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>06</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>08</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Subbian, Natarajan and Murugan.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Subbian, Natarajan and Murugan</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>
<bold>Introduction:</bold> Water scarcity and water pollution are two major issues in India. Circular economy-based wastewater treatment technology provides the most sustainable solutions for solving these issues. In this paper, a novel multi-objective decentralized controller (MODC) is proposed for benchmarking a multi-input multi-output (MIMO) activated sludge wastewater treatment plant (WWTP) to achieve maximum effluent quality with minimum cost. WWTPs with conventional control schemes consume more energy to achieve the desired effluent quality.</p>
<p>
<bold>Methods:</bold> In this study, a MIMO model is developed for the activated sludge process (ASP) from a physics-based model, and relative gain array (RGA) analysis are carried out to determine the interaction between the loops to identify a suitable control scheme for the MIMO process. In addition, a multi-objective decentralized control problem is formulated to achieve the conflicting multiple objectives of improving effluent quality and minimizing operational costs by efficient usage of energy.</p>
<p>
<bold>Results and discussion:</bold> The desired quality and cost reduction are verified by comparing the integral square error (ISE) and control effort (CE) values of a closed-loop WWTP. A multi-objective evolutionary algorithm (MOEA), namely, the non-dominated sorting genetic algorithm (NSGA)-II, successfully solves the multi-objective control problem. NSGA-II provides several optimal solutions in the Pareto front. In order to demonstrate the feasibility of the proposed controller, three optimal solutions are selected from the Pareto-optimal front, and their closed-loop performances are evaluated qualitatively and quantitatively for both servo and regulatory operations. Improving the quality of effluent enhances active sludge production, which in turn increases the methane production in the anaerobic digester.</p>
</abstract>
<kwd-group>
<kwd>circular economy</kwd>
<kwd>activated sludge process</kwd>
<kwd>MODC</kwd>
<kwd>WWTP</kwd>
<kwd>MIMO</kwd>
<kwd>MOEA</kwd>
<kwd>RGA</kwd>
<kwd>NSGA-II</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Sustainable Process Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Wastewater treatment plays a major role in developing countries like India due to water scarcity and environmental pollution. The population rise in major cities leads to more water consumption and wastewater discharge into natural water bodies. According to the Central Ground Water Authority, Government of India report, an average individual utilizes 150&#x2013;200&#xa0;L of water per day. Only 10&#xa0;L is consumed for cooking and drinking; 40&#xa0;L per day is spent on washing dishes, washing clothes, and floor cleaning; and the remaining 85&#xa0;L is utilized for toilet flushing and bathing. In general, the 80&#xa0;L of wastewater used for washing and bathing could be reused for gardening, toilet flushing, and other useful applications with proper treatment. Many cities in India have centralized wastewater treatment plants (WWTPs) with conventional treatment techniques. However, most of them are less effective because they require underground pipe lining, pumping stations, and extensive maintenance. In contrast, a WWTP with advanced technologies would improve the quality of water and reduce maintenance costs (<xref ref-type="bibr" rid="B6">Grobelak et al., 2019</xref>). These technologies could lead to the design of decentralized WWTPs that deploy sophisticated techniques closer to the water generation point to enable more effective reuse of water. Two types of treatment plants, namely, sewage treatment plants (STPs) and effluent treatment plants (ETPs), are used to treat sewage and industrial wastewater, respectively (<xref ref-type="bibr" rid="B15">Negi et al., 2022</xref>). In a WWTP, wastewater undergoes three stages: primary treatment, secondary treatment, and tertiary treatment.</p>
<p>Generally, untreated sewage has biological oxygen demand (BOD) and total solid substrate (TSS) values in the range of 250&#x2013;350&#xa0;mg/L and 150&#x2013;200&#xa0;mg/L, respectively. The characteristics of the effluent are predominantly influenced by its levels of BOD and TSS. A considerable proportion of sewage treatment plants (STPs) yield effluent containing BOD concentrations lower than 20&#xa0;mg/L and TSS concentrations lower than 10&#xa0;mg/L. In order to improve the quality of effluent, ETP requires conventional biological wastewater treatment techniques like the activated sludge process (ASP) to control dissolved oxygen concentration and substrate concentration in the effluent. The ASP is an effective secondary treatment process that uses an aeration tank and a sludge-settling tank. Organic matter in the sewage is digested by aerobic bacteria in the aeration tank. Aerobic bacteria within the aeration tank facilitate the digestion of organic substances present in sewage. Ensuring the growth of these aerobic bacteria necessitates an effective control system. According to <xref ref-type="bibr" rid="B12">Malaysia (2011)</xref>, there&#x2019;s a notable need for advanced control strategies to enhance or uphold the effluent quality of wastewater treatment plants (WWTPs). As a result, the processes of modeling, control variable selection, and designing control structures assume crucial importance in the pursuit of optimizing control objectives.</p>
<p>
<xref ref-type="bibr" rid="B7">Hodasz et al. (2016)</xref> developed a linear model to estimate the unmeasured process parameters of an activated sludge model (ASM1). They also suggested that good control of a WWTP leads to better water quality and an efficient use of energy. <xref ref-type="bibr" rid="B1">Castro-Amoedo et al. (2022)</xref> presented an overview of energy potentials present in a WWTP and demonstrated the concepts for local energy supply based on integrated perspectives. The control of dissolved oxygen levels in WWTPs has considerably reduced the energy required for aeration purposes in aerobic-activated sludge plants. <xref ref-type="bibr" rid="B3">de Araujo et al. (2013)</xref> discussed the application of sensitivity analysis for the selection of economic controlled variables for the optimal operation of a WWTP. <xref ref-type="bibr" rid="B19">Samuelsson et al. (2005)</xref> analyzed the cross coupling between loops in a bioreactor model using the relative gain array (RGA) and Hankel interaction index array (HIIA) methods. In addition, they compared the performances of a decentralized control scheme and a multivariable control scheme. <xref ref-type="bibr" rid="B23">Vilanova and Alfaro (2009)</xref> proposed a multi-loop decentralized control strategy for the control of a WWTP based on an ASP. <xref ref-type="bibr" rid="B24">Vinayak and Narayanan (2017)</xref> attempted a centralized nonlinear state feedback proportional integral (PI) controller for an ASP and compared its performance with a centralized linear state feedback controller. <xref ref-type="bibr" rid="B16">Nor&#x2019;Azlan et al. (2017)</xref> developed a centralized controller for a linear system using the Davison, Penttinen&#x2013;Koivo, and Maciejowski methods. They have also tuned the centralized controller parameters using single objective optimization techniques, such as a genetic algorithm (GA), particle swarm optimization (PSO), and the bat algorithm. <xref ref-type="bibr" rid="B8">Holenda et al. (2007</xref>, <xref ref-type="bibr" rid="B9">2008)</xref> applied a model predictive controller for an activated sludge wastewater treatment process to control the dissolved oxygen level in the reactor and optimize the effluent quality and aeration energy cost against the number of aeration cycles required per day using a GA. <xref ref-type="bibr" rid="B21">Sweetapple et al. (2014)</xref> demonstrated the potential of multi-objective optimization of WWTP control strategies for the reduction of greenhouse gas emissions. <xref ref-type="bibr" rid="B2">Chen et al. (2014)</xref> optimized benchmark simulation model 1 (BSM1) parameters to achieve the percentage of effluent violation (PEV), overall cost index (OCI), total volume (TV), and total solid substrate (TSS) using a non-dominated sorting genetic algorithm (NSGA)-II.</p>
<p>A multi-objective optimization is a valuable tool in control engineering for controller design and process optimization. <xref ref-type="bibr" rid="B10">Hu et al. (2014)</xref> proposed a multi-objective nonlinear PID controller parameter for greenhouse systems to minimize cost and maximize performance using NSGA-II. <xref ref-type="bibr" rid="B22">Tian et al. (2023)</xref> recently proposed that the installation position parameters of a submersible mixer can be optimized using multi-input multi-output (MIMO) for improved performance.</p>
<p>Compromising the performance and cost becomes a critical issue in WWTP. In recent years, applications of multi-objective optimization have played a major role in WWTP for control and process optimization. <xref ref-type="bibr" rid="B18">Reynoso-Meza and Carre&#xf1;o-Alvarado (2019)</xref> designed a multi-objective controller for a multivariable active sludge process with nitrification and denitrification. <xref ref-type="bibr" rid="B17">Qiao and Zhang (2018)</xref> developed a multi-objective control scheme for the Benchmark Simulation Model 1 (BSM1) MIMO process. <xref ref-type="bibr" rid="B5">Ghimire et al. (2021)</xref> reviewed two routes for developing more sustainable and circular economy-based wastewater treatment systems and challenges in those approaches. <xref ref-type="bibr" rid="B11">Kaur (2023)</xref> applied multi-objective optimization to improve the quality of WWTP effluent based on a data-driven model. The ASP utilizes suspended growth microorganisms to effectively remove organic matter and nutrients from wastewater, but it demands substantial energy and land resources and generates sludge. Trickling filters, in contrast, employ attached growth microorganisms and offer simpler operations with lower energy requirements but may not efficiently remove nutrients and face potential media clogging issues. Membrane bioreactors (MBRs) combine biological treatment with membrane filtration, producing high-quality effluent with reduced sludge, yet they come with higher capital and operational costs and are susceptible to membrane fouling. Sequential batch reactors (SBRs) feature flexible operation and good nutrient removal but may cause variations due to intermittent operation and require longer treatment time. Constructed wetlands provide an eco-friendly and aesthetically pleasing approach with low operating costs but necessitate larger land areas and may not be suitable for high-strength industrial wastewater. Advanced oxidation processes (AOPs) utilize chemical reactions to effectively remove recalcitrant pollutants without producing sludge, but their high cost and potential secondary waste generation remain drawbacks. Ultimately, the choice of the wastewater treatment model depends on the specific characteristics of the wastewater, treatment goals, available resources, and environmental considerations. Often, a combination of different treatment processes is adopted to optimize efficiency and achieve the desired effluent quality.</p>
<p>Much literature has been reported on circular-economy-based wastewater treatments (<xref ref-type="bibr" rid="B14">Neczaj and Grosser, 2018</xref>; <xref ref-type="bibr" rid="B20">Spreafico, 2022</xref>; <xref ref-type="bibr" rid="B25">Zhang and Liu, 2022</xref>). These studies have reported on circular economics in WWTPs and challenges, barriers that prevent the implementation of circular economics, and analysis and design for circular economics.</p>
<p>To our knowledge, using a multi-objective decentralized controller (MODC) for WWTP with ASM1 to control effluent substrate concentration and DO concentration by maximizing effluent quality and minimizing production cost has not yet been explored. MIMO systems allow for better control performance by considering multiple variables simultaneously. They can take into account cross-coupling effects between inputs and outputs, which often exist in interconnected processes. By controlling multiple variables concurrently, MIMO models can achieve tighter control and improved overall system stability. The main aim of this paper is to design a multi-objective decentralized PI controller for an ASP to achieve the conflicting multi-objectives of maximizing the performance of effluent quality and minimizing the production cost.</p>
<p>The main contributions of the proposed work are i) identifying inputs and outputs and developing an ASM1 model for a WWTP, ii) identification of a suitable controller (centralized/decentralized) for a MIMO process, iii) multi-objective control problem formulation for a decentralized controller, and iv) design of a circular economy-based MODC for WWTP using NSGA-II.</p>
<p>The paper is organized as follows: the modeling of benchmark ASP is described in <xref ref-type="sec" rid="s2">Section 2</xref>. RGA analysis and multi-objective problem formulation are elucidated in <xref ref-type="sec" rid="s3">Section 3</xref>. The design of the proposed MODC is discussed in <xref ref-type="sec" rid="s4">Section 4</xref>. <xref ref-type="sec" rid="s5">Section 5</xref> presents closed-loop simulation studies of ASP. The concluding remarks are presented at the end.</p>
</sec>
<sec id="s2">
<title>2 Modeling of activated sludge process</title>
<sec id="s2-1">
<title>2.1 Benchmark ASP</title>
<p>Secondary treatment plays a major role in treating wastewater in STPs. Secondary treatment consists of an aerobic bioreactor for aeration purposes and a settling tank, which is used to settle the sludge at the bottom while effluent is collected at the top of the tank. The schematic diagram of the ASP is given in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Schematic diagram of activated sludge process.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g001.tif"/>
</fig>
<p>The aerator is taken to be a well-stirred tank for complete mixing in which microorganisms start reacting with organic pollutants in the wastewater. Six major groups of microorganisms, namely, bacteria, protozoa, metazoa, filamentous bacteria, algae, and fungi, are generally found in the aeration basin of the ASP. Microorganisms take the organic pollutants as food material and start growing; they collide to form sludge. Then, the suspended microorganisms are given to the settling tank (gravitational settling tank). In order to maintain the microorganism level in the tank, a part of the concentrated biomass is recycled back to the bioreactor, and the remaining material is removed from the settling tank. This study assumes that there is no microorganism growth in the settling tank.</p>
</sec>
<sec id="s2-2">
<title>2.2 Modeling of ASP</title>
<p>Modeling and control of an ASP is more complex due to biological reactions and variations in the influent characteristics. The ASM1 model is developed from the real-time ASP under the following assumptions: i) complete mixing takes place in the aeration tank, ii) the substrate concentration of influent wastewater is constant, iii) all reactions take place only in the aeration basin, and iv) no reaction takes place in the settling tank. Hence, the concentration of substrate in the aerobic bioreactor is the same as the concentration of substrate in the settling tank. <xref ref-type="bibr" rid="B7">Hodasz et al. (2016)</xref> reported the dynamic model of the ASP, which is given in Eqs <xref ref-type="disp-formula" rid="e1">1</xref>&#x2013;<xref ref-type="disp-formula" rid="e5">5</xref>.</p>
<p>Mass balance equation of biomass in aeration tank:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>Mass balance equation of substrate present in aeration tank:<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#xb5;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mi>Y</mml:mi>
</mml:mfrac>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
</mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Mass balance equation of dissolved oxygen concentration:<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#xb5;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mi>Y</mml:mi>
</mml:mfrac>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Mass balance equation in the activated sludge in the settling tank:<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Biomass growth rate in aeration tank:<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:mi>&#xb5;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#xb5;</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>where <inline-formula id="inf1">
<mml:math id="m6">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> are the state variables. In the ASP, the dilution rate <inline-formula id="inf2">
<mml:math id="m7">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> and the air flow rate <inline-formula id="inf3">
<mml:math id="m8">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> are taken as manipulated variables to control the substrate <inline-formula id="inf4">
<mml:math id="m9">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> and dissolved oxygen <inline-formula id="inf5">
<mml:math id="m10">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, and the biomass production factor <inline-formula id="inf6">
<mml:math id="m11">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>Y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is taken to be the disturbance variable.</p>
<p>The transfer function plays a pivotal role in wastewater treatment plant and significantly impacts its performance. It represents the relationship between input variables, such as influent characteristics or control actions, and the corresponding output variables, usually effluent quality or process variables. Understanding the transfer function allows process engineers and operators to gain insight into the dynamic behavior of the WWTP and its response to changes in influent conditions or operational parameters. With this knowledge, they can design effective control strategies, optimize process conditions, and enhance treatment efficiency. Additionally, transfer functions are vital for process modeling and simulation, aiding in predicting system behavior under different scenarios and facilitating process optimization. They also play a crucial role in control system design and performance evaluation, enabling the WWTP to achieve stable operation, comply with regulatory standards, and consistently deliver high-quality effluent while minimizing environmental impact.</p>
<p>The transfer function model for MIMO is developed by linearizing the dynamic equations given in Eq. <xref ref-type="disp-formula" rid="e1">1</xref>&#x2013;<xref ref-type="disp-formula" rid="e4">4</xref> around the steady-state operating points tabulated in <xref ref-type="table" rid="T1">Table 1</xref> (<xref ref-type="bibr" rid="B7">Hodasz et al., 2016</xref>). The transfer function model is presented in Eq. <xref ref-type="disp-formula" rid="e6">6</xref>:<disp-formula id="e6">
<mml:math id="m12">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>;</mml:mo>
<mml:mi>W</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where<disp-formula id="equ1">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>688.8</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>4261</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>680.9</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.4246</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.651</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.277</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1.136</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.01056</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ2">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.01268</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.002405</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2.457</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:msup>
<mml:mn>10</mml:mn>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.651</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.277</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1.136</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.01056</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ3">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>11.1</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>278.7</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>52.54</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.5329</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.651</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.277</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1.136</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.01056</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ4">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>O</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.06151</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.06734</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.01078</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.0001004</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.651</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>7.277</mml:mn>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1.136</mml:mn>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.01056</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Steady-state operating points.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Variable</th>
<th align="left">Steady-state value</th>
<th align="left">Unit</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Biomass (<italic>X</italic>)</td>
<td align="left">565.5</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Substrate (<italic>s</italic>)</td>
<td align="left">42.35</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Dissolved oxygen (<italic>DO</italic>)</td>
<td align="left">6.442</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Recycled biomass (<italic>X</italic>
<sub>
<italic>r</italic>
</sub>)</td>
<td align="left">1244</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Dissolved oxygen concentration (<italic>DO</italic>
<sub>
<italic>i</italic>
</sub>)</td>
<td align="left">0.5</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Influent substrate (<italic>S</italic>
<sub>
<italic>in</italic>
</sub>)</td>
<td align="left">765</td>
<td align="left">[<italic>mg/L</italic>]</td>
</tr>
<tr>
<td align="left">Model constant (<italic>K</italic>
<sub>
<italic>0</italic>
</sub>)</td>
<td align="left">0.46</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Affinity constant (<italic>K</italic>
<sub>
<italic>s</italic>
</sub>)</td>
<td align="left">181</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Saturation constant (<italic>K</italic>
<sub>
<italic>DO</italic>
</sub>)</td>
<td align="left">0.2</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Ratio of waste flow to influent flow (<inline-formula id="inf7">
<mml:math id="m17">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="left">0.015</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Ratio of recycled flow to influent flow (<italic>r</italic>)</td>
<td align="left">0.8</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Oxygen mass transfer coefficient (<inline-formula id="inf8">
<mml:math id="m18">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="left">0.017</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Biomass production factor (<italic>&#x3b3;</italic>)</td>
<td align="left">0.54</td>
<td align="left">&#x2014;</td>
</tr>
<tr>
<td align="left">Air flow rate (<italic>W</italic>)</td>
<td align="left">340</td>
<td align="left">[<italic>m</italic>
<sup>3</sup>/&#x210e;]</td>
</tr>
<tr>
<td align="left">Dilution rate (<italic>D</italic>)</td>
<td align="left">0.0638</td>
<td align="left">[&#x210e;&#x2212;1]</td>
</tr>
<tr>
<td align="left">Biomass growth rate (<inline-formula id="inf9">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3bc;</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="left">0.229</td>
<td align="left">[&#x210e;&#x2212;1]</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The aforementioned transfer function model is used for RGA analysis to identify a suitable control scheme for the ASP in the subsequent section.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Identification of control strategy and multi-objective control problem formulation</title>
<p>In this section, RGA analysis is applied to the WWTP transfer function model, and the multi-objective control problem formulation is discussed<bold>.</bold>
</p>
<sec id="s3-1">
<title>3.1 RGA analysis</title>
<p>RGA is used to determine the best input (u)&#x2013;output (y) pairing for multivariable systems. It quantifies the change in steady-state gain between an input&#x2013;output pair that occurs when other control loops are closed. It provides a measure of steady-state interactions from gain. The relative gain between input j and output i is given by Eq. <xref ref-type="disp-formula" rid="e7">7</xref>:<disp-formula id="e7">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where<inline-formula id="inf10">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the open-loop gain with all other control loops open, and<inline-formula id="inf11">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="italic">&#x2202;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the closed-loop gain with all other control loops closed.</p>
<p>The ASP is a 2 &#xd7; 2 system; its RGA structure is presented in Eq. <xref ref-type="disp-formula" rid="e8">8</xref>:<disp-formula id="e8">
<mml:math id="m23">
<mml:mrow>
<mml:mo>&#x2227;</mml:mo>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>The steady-state gain matrix of the ASP is computed from Eq. <xref ref-type="disp-formula" rid="e6">6</xref> given in Eq. <xref ref-type="disp-formula" rid="e9">9</xref>:<disp-formula id="e9">
<mml:math id="m24">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>40.20833</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.00232670</mml:mn>
<mml:mo>;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>50.4640</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mn>0.0095075</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Using the gain matrix given in Eq. <xref ref-type="disp-formula" rid="e9">9</xref>, the RGA matrix is developed using Eq. <xref ref-type="disp-formula" rid="e10">10</xref>:<disp-formula id="e10">
<mml:math id="m25">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>11</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>12</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>21</mml:mn>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mn>22</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1.443293</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.443293</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.443293</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mn>1.443293</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>From the obtained RGA matrix, it is recommended that <inline-formula id="inf12">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> be paired with <inline-formula id="inf13">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf14">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> be paired with <inline-formula id="inf15">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. The obtained RGA matrix shows that there is a minimum interaction between the process loops. Hence, a decentralized control scheme is preferred for an ASP. <xref ref-type="bibr" rid="B13">Muntean et al. (2015)</xref> also carried out RGA analysis for the WWTP process with recycled sludge flow to influent flow and aeration rate from the bioreactor as manipulated variables and used the substrate and dissolved oxygen as controlled variables and reported that a decentralized control scheme is more suitable for the WWTP process with the identified model due to fewer interactions between the loops.</p>
</sec>
<sec id="s3-2">
<title>3.2 Decentralized PI controller</title>
<p>As stated earlier, the interaction between inputs and outputs is weak; hence, a decentralized control architecture with a PI controller is proposed. The decentralized control is given in Eq. <xref ref-type="disp-formula" rid="e11">11</xref>:<disp-formula id="e11">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>11</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>22</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>where <inline-formula id="inf16">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>11</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and<disp-formula id="equ5">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>22</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>To tune the PI controller parameters, a multi-objective control problem is formulated in the next subsection.</p>
</sec>
<sec id="s3-3">
<title>3.3 Multi-objective problem formulation (MOO)</title>
<p>The multi-objective control problem is formulated by identifying the PI control structure to improve the quality of effluent by minimizing the integral square error (ISE) and reducing production costs by minimizing control effort (CE). The conflicting objective functions of the multi-objective control problem are given in Equations <xref ref-type="disp-formula" rid="e12">(12)</xref> and <xref ref-type="disp-formula" rid="e13">(13)</xref>:<disp-formula id="e12">
<mml:math id="m33">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
</mml:mstyle>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m34">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
</mml:mstyle>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>The aforementioned two objective functions are optimized using eight decision variables. <inline-formula id="inf17">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf18">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are proportional gain and integral gain of the substrate (loop 1), respectively.<inline-formula id="inf19">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf20">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are proportional gain and integral gain of dissolved oxygen (loop 2), respectively.<inline-formula id="inf21">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf22">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are weighting factors for loop 1 and loop 2, respectively.<inline-formula id="inf23">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf24">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are weighting factors for loop 1 and loop 2, respectively.<inline-formula id="inf25">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf26">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are errors in loop 1 and loop 2, respectively.<inline-formula id="inf27">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf28">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are control efforts of loop 1 and loop 2, respectively.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Design of a multi-objective decentralized controller</title>
<p>Multi-objective optimization involves more than one objective function, and the functions must be optimized simultaneously. The multi-objective control problem is commonly solved by evolutionary algorithms (EAs). The primary reason for using EAs is their ability to find multiple Pareto-optimal solutions in one single simulation run. <xref ref-type="bibr" rid="B4">Deb et al. (2002)</xref> reported the performance comparisons of various multi-objective evolutionary algorithms (MOEAs) such as a niched Pareto genetic algorithm (NPGA), a Pareto-archived evolution strategy (PAES), a strength Pareto evolutionary algorithm (SPEA) and a non-dominant sort genetic algorithm-II (NSGA-II) for different test problems. Most EAs have been successfully applied for design optimization and controller optimization over the past two decades. However, NSGA-II outperforms other EAs for most multi-objective problems. It can provide a better spread of solutions and better convergence near the true Pareto-optimal front. In this paper, NSGA-II is selected to solve the multi-objective control problem of ASP.</p>
<sec id="s4-1">
<title>4.1 Non-dominant sort genetic algorithm-II</title>
<p>
<xref ref-type="bibr" rid="B4">Deb et al. (2002)</xref> proposed NSGA-II, which alleviates the difficulty of NSGA. The important features of NSGA-II are as follows:<list list-type="simple">
<list-item>
<p>&#x2022; NSGA-II is based upon a GA.</p>
</list-item>
<list-item>
<p>&#x2022; It is an extended version of NSGA.</p>
</list-item>
<list-item>
<p>&#x2022; It overcomes the difficulties faced in NSGA, which are complexity, lack of elitism, and choosing the optimal parameter for the sharing parameter.</p>
</list-item>
<list-item>
<p>&#x2022; The main aim of NSGA-II is to preserve elitism and diversity.</p>
</list-item>
<list-item>
<p>&#x2022; NSGA-II gives a Pareto-optimal front, which consists of a set of optimal solutions for conflicting objective functions.</p>
</list-item>
</list>
</p>
<p>The algorithm and detailed implementation procedure of NSGA-II are found in Kalyanmoy <xref ref-type="bibr" rid="B4">Deb (2002)</xref>. The flowchart of NSGA-II is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Flow chart of NSGA-II.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g002.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Implementation of NSGA-II</title>
<p>The parameters of NSGA-II used for execution are given in <xref ref-type="table" rid="T2">Table 2</xref>. The Pareto-optimal front provides the optimal solutions for the conflicting objective functions, and it is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Parameters of NSGA-II.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Point</th>
<th align="left">Process loop</th>
<th align="left">
<italic>K</italic>
<sub>
<italic>C</italic>
</sub>
</th>
<th align="left">
<italic>K</italic>
<sub>
<italic>i</italic>
</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="2" align="left">K1</td>
<td align="left">LOOP 1</td>
<td align="left">0.4833</td>
<td align="left">0.9666</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">120.8256</td>
<td align="left">241.6513</td>
</tr>
<tr>
<td rowspan="2" align="left">K2</td>
<td align="left">LOOP 1</td>
<td align="left">0.2534</td>
<td align="left">0.5068</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">63.3575</td>
<td align="left">126.7151</td>
</tr>
<tr>
<td rowspan="2" align="left">K3</td>
<td align="left">LOOP 1</td>
<td align="left">0.0796</td>
<td align="left">0.1593</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">19.9158</td>
<td align="left">39.8316</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Pareto optimal front of NSGA-II.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g003.tif"/>
</fig>
<p>To study the trade-off between the performance and cost of production, three controller parameters at points K1, K2, and K3 are selected from the Pareto-optimal front and are tabulated in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Controller parameters for the selected three points.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Point</th>
<th align="left">Process loop</th>
<th align="left">
<italic>K</italic>
<sub>
<italic>C</italic>
</sub>
</th>
<th align="left">
<italic>K</italic>
<sub>
<italic>i</italic>
</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="2" align="left">K1</td>
<td align="left">LOOP 1</td>
<td align="left">0.4833</td>
<td align="left">0.9666</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">120.8256</td>
<td align="left">241.6513</td>
</tr>
<tr>
<td rowspan="2" align="left">K2</td>
<td align="left">LOOP 1</td>
<td align="left">0.2534</td>
<td align="left">0.5068</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">63.3575</td>
<td align="left">126.7151</td>
</tr>
<tr>
<td rowspan="2" align="left">K3</td>
<td align="left">LOOP 1</td>
<td align="left">0.0796</td>
<td align="left">0.1593</td>
</tr>
<tr>
<td align="left">LOOP 2</td>
<td align="left">19.9158</td>
<td align="left">39.8316</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s5">
<title>5 Closed-loop simulation studies of ASP</title>
<sec id="s5-1">
<title>5.1 Servo response</title>
<p>The ASP is simulated using the set of differential equations given in the modeling section. Then, the closed-loop servo responses are obtained for the selected controllers and qualitatively compared in <xref ref-type="fig" rid="F4">Figure 4</xref>, and associated manipulated variables are shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. For quantitative comparison, the performance measures are given in <xref ref-type="table" rid="T4">Table 4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Servo response of Activated Sludge Process <bold>(A)</bold> Substrate <bold>(B)</bold> Dissolved oxygen.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Variations in manipulated variables under servo operation <bold>(A)</bold> Dilution rate <bold>(B)</bold> Airflow rate.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g005.tif"/>
</fig>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Performance measures of MODC.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="3" align="left">Point</th>
<th colspan="2" align="left">Servo</th>
<th colspan="2" align="left">Regulatory</th>
</tr>
<tr>
<th colspan="2" align="left">Operation</th>
<th colspan="2" align="left">Operation</th>
</tr>
<tr>
<th align="left">ISE</th>
<th align="left">CE</th>
<th align="left">ISE</th>
<th align="left">CE</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">K1</td>
<td align="left">0.0715</td>
<td align="left">1.978&#xd7;10<sup>6</sup>
</td>
<td align="left">0.0878</td>
<td align="left">1.19&#xd7;10<sup>6</sup>
</td>
</tr>
<tr>
<td align="left">K2</td>
<td align="left">0.0817</td>
<td align="left">9.947&#xd7;10<sup>5</sup>
</td>
<td align="left">0.1026</td>
<td align="left">5.979&#xd7;10<sup>5</sup>
</td>
</tr>
<tr>
<td align="left">K3</td>
<td align="left">0.09451</td>
<td align="left">2.558&#xd7;10<sup>5</sup>
</td>
<td align="left">0.1152</td>
<td align="left">1.56&#xd7;10<sup>5</sup>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The servo responses with controller gains at the K<sub>1</sub> point do not provide significant variation in the substrate concentration. However, the responses show the variation in the dissolved oxygen concentration in the effluent. From <xref ref-type="fig" rid="F4">Figures 4B</xref>,<xref ref-type="fig" rid="F5">5B</xref>, it is observed that the process with controller gains at the K<sub>1</sub> point takes less retention time and settles very quickly at the steady-state value due to large CE (refer <xref ref-type="table" rid="T4">Table 4</xref>). It is observed from <xref ref-type="table" rid="T4">Table 4</xref> that controller gains at the K<sub>1</sub> point give a much smaller integral square error. Hence, in order to improve the performance, the closed-loop system must compromise on CE.</p>
<p>The process with controller gains at the K<sub>2</sub> point gives the trade-off between the performance and CE (refer <xref ref-type="table" rid="T4">Table 4</xref>). The retention time of the process with controller gains at the K<sub>2</sub> point shown in <xref ref-type="fig" rid="F4">Figures 4B</xref>, <xref ref-type="fig" rid="F5">5B</xref> is larger than the controller gains at the K<sub>1</sub> point; meanwhile, it takes a smaller CE than a process with controller gains at the K<sub>1</sub> point.</p>
<p>The closed-loop servo response with controller gains at K<sub>3</sub> shown in <xref ref-type="fig" rid="F4">Figures 4B</xref>, <xref ref-type="fig" rid="F5">5B</xref> exhibited a longer retention time than other controllers. The controller gains at the K<sub>3</sub> point have less CE (refer <xref ref-type="table" rid="T4">Table 4</xref>). At the same time, it should compromise on the performance of the system.</p>
</sec>
<sec id="s5-2">
<title>5.2 Regulatory response</title>
<p>To study the regulatory operation of the ASP, biomass production factor (Y) is selected as a disturbance and is changed from 0.54 to 0.594&#xa0;at the tenth hour. The simulated regulatory responses of loop 1 and loop 2 are shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, and the corresponding manipulated variables are shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. The quantitative comparison of the regulatory operation with various controllers is tabulated in <xref ref-type="table" rid="T4">Table 4</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Regulatory response of activated sludge process <bold>(A)</bold> Substrate <bold>(B)</bold> Dissolved oxygen.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Variations in manipulated variables under regulatory operation <bold>(A)</bold> Dilution rate <bold>(B)</bold> Airflow rate.</p>
</caption>
<graphic xlink:href="fceng-05-1235125-g007.tif"/>
</fig>
<p>It is found from <xref ref-type="fig" rid="F6">Figures 6A, B</xref> that the effect of disturbance on substrate concentration and dilution rate is less than the effect on dissolved oxygen concentration and air flow rate. In quantitative comparison, it is noted that the CE required for regulatory operation is less than the CE required by the servo operation. However, the ISEs provided by the regulatory operation are larger than those of the servo operation. Hence, minimization of CE will affect the performance of the process. Meanwhile, minimization of CE reduces the production cost of freshwater.</p>
<p>The circular economy in wastewater treatment leads to freshwater, nutrient, and energy recovery from wastewater. Improvement in effluent water quality leads to the reuse of water for domestic applications, which will reduce water stress. Furthermore, nutrients recovered from active sludge can be used as fertilizers. Energy recovery can also be done by producing biogas from active sludge.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>In this work, a novel circular economy-based MODC design is proposed for the MIMO ASP. The transfer function model is developed from nonlinear differential equations to select an appropriate control scheme to control dissolved oxygen concentration and substrate concentration in the effluent of the ASP. Then, RGA analysis is performed to determine the interactions between the loops. From the interaction analysis, it is concluded that a decentralized controller is the best control scheme for the ASP due to less interaction between the loops. A multi-objective control problem is formulated to improve effluent quality and reduce costs. The multi-objective control problem is successfully solved by the MOEA NSGA-II. Some of the conventional wastewater treatment techniques are simple and have low production costs; however, they require a long treatment time and result in poor effluent quality. Whereas sophisticated treatment techniques are more expensive due to high production costs but can provide good-quality effluent. The advantage of the proposed multi-objective controller is that it can provide a number of optimal solutions in the Pareto-optimal front. The users can select their own solutions according to the applications from the Pareto-optimal front. The feasibility of the proposed controller is demonstrated by selecting three controller parameters from the Pareto-optimal front, and their closed-loop servo and regulatory performances are compared qualitatively as well as quantitatively. The proposed circular economy-based multi-objective controller will be very useful for the reuse of wastewater, nutrients, and energy.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>SS: writing&#x2014;original draft, conceptualization, methodology, and software. CM: writing&#x2014;review, editing, conceptualization, methodology, and software. PN: writing and conceptualization. All authors contributed to the article and approved the submitted version.</p>
</sec>
<ack>
<p>The authors gratefully acknowledge the Department of Instrumentation, Madras Institute of Technology, for their support with conducting the experiments.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Castro-Amoedo</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Morisod</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Granacher</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Mar&#xe9;chal</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>The role of biowaste: A multi-objective optimization platform for combined heat, power and fuel</article-title>. <source>Front. Energy Res.</source> <volume>9</volume>, <fpage>718310</fpage>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.718310</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Optimal design activated sludge process by means of multi-objective optimization: case study in benchmark simulation model 1 (BSM1)</article-title>. <source>Water Sci. Technol.</source> <volume>69</volume> (<issue>10</issue>), <fpage>2052</fpage>&#x2013;<lpage>2058</lpage>. <pub-id pub-id-type="doi">10.2166/wst.2014.119</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>de Araujo</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Gallani</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mulas</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Skogestad</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Sensitivity analysis of optimal operation of an activated sludge process model for economic controlled variable selection</article-title>. <source>Industrial Eng. Chem. Res.</source> <volume>52</volume> (<issue>29</issue>), <fpage>9908</fpage>&#x2013;<lpage>9921</lpage>. <pub-id pub-id-type="doi">10.1021/ie4006673</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deb</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Pratap</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Agarwal</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Meyarivan</surname>
<given-names>T. A. M. T.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>A fast and elitist multiobjective genetic algorithm: NSGA-II</article-title>. <source>IEEE Trans. Evol. Comput.</source> <volume>6</volume> (<issue>2</issue>), <fpage>182</fpage>&#x2013;<lpage>197</lpage>. <pub-id pub-id-type="doi">10.1109/4235.996017</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghimire</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Sarpong</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Gude</surname>
<given-names>V. G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Transitioning wastewater treatment plants toward circular economy and energy sustainability</article-title>. <source>ACS omega</source> <volume>6</volume> (<issue>18</issue>), <fpage>11794</fpage>&#x2013;<lpage>11803</lpage>. <pub-id pub-id-type="doi">10.1021/acsomega.0c05827</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grobelak</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Grosser</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kacprzak</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kamizela</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Sewage sludge processing and management in small and medium-sized municipal wastewater treatment plant-new technical solution</article-title>. <source>J. Environ. Manag.</source> <volume>234</volume>, <fpage>90</fpage>&#x2013;<lpage>96</lpage>. <pub-id pub-id-type="doi">10.1016/j.jenvman.2018.12.111</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Hodasz</surname>
<given-names>N. I.</given-names>
</name>
<name>
<surname>Br&#x103;dil&#x103;</surname>
<given-names>V. I.</given-names>
</name>
<name>
<surname>Nascu</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Lendek</surname>
<given-names>Z.</given-names>
</name>
</person-group> &#x201c;<article-title>Modeling and parameter estimation for an activated sludge wastewater treatment process</article-title>,&#x201d; in <conf-name>Proceedings of the 2016 IEEE International Conference on Automation, Quality and Testing, Robotics (AQTR)</conf-name>, <conf-loc>Cluj-Napoca, Romania</conf-loc>, <conf-date>2016, May</conf-date>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holenda</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Domokos</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>R&#xe9;dey</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Fazakas</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Aeration optimization of a wastewater treatment plant using genetic algorithm</article-title>. <source>Optim. control Appl. methods</source> <volume>28</volume> (<issue>3</issue>), <fpage>191</fpage>&#x2013;<lpage>208</lpage>. <pub-id pub-id-type="doi">10.1002/oca.796</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holenda</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Domokos</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Redey</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Fazakas</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Dissolved oxygen control of the activated sludge wastewater treatment process using model predictive control</article-title>. <source>Comput. Chem. Eng.</source> <volume>32</volume> (<issue>6</issue>), <fpage>1270</fpage>&#x2013;<lpage>1278</lpage>. <pub-id pub-id-type="doi">10.1016/j.compchemeng.2007.06.008</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Goodman</surname>
<given-names>E. D.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>NSGA-II-based nonlinear PID controller tuning of greenhouse climate for reducing costs and improving performances</article-title>. <source>Neural Comput. Appl.</source> <volume>24</volume>, <fpage>927</fpage>&#x2013;<lpage>936</lpage>. <pub-id pub-id-type="doi">10.1007/s00521-012-1312-8</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="thesis">
<person-group person-group-type="author">
<name>
<surname>Kaur</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2023</year>). &#x201c;<article-title>Modeling and multi-objective optimization of wastewater treatment process</article-title>,&#x201d;. <comment>Doctoral dissertation</comment> (<publisher-loc>Canada</publisher-loc>: <publisher-name>The University of Western Ontario</publisher-name>).</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Malaysia</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Control strategies of wastewater treatment plants</article-title>. <source>Aust. J. Basic Appl. Sci.</source> <volume>5</volume> (<issue>8</issue>), <fpage>446</fpage>&#x2013;<lpage>455</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Muntean</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Both</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Crisan</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Nascu</surname>
<given-names>I.</given-names>
</name>
</person-group> &#x201c;<article-title>RGA analysis and decentralized control for a wastewater treatment plant</article-title>,&#x201d; in <conf-name>Proceedings of the 2015 IEEE International Conference on Industrial Technology (ICIT)</conf-name>, <conf-loc>Seville, Spain</conf-loc>, <conf-date>March 2015</conf-date>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Neczaj</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Grosser</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Circular economy in wastewater treatment plant&#x2013;challenges and barriers</article-title>. <source>Proceedings</source> <volume>2</volume> (<issue>11</issue>), <fpage>614</fpage>. <pub-id pub-id-type="doi">10.3390/proceedings2110614</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Negi</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Verma</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Singh</surname>
<given-names>S. P.</given-names>
</name>
<name>
<surname>Mahapatra</surname>
<given-names>B. S.</given-names>
</name>
<name>
<surname>Jatav</surname>
<given-names>H. S.</given-names>
</name>
</person-group> (<year>2022</year>). &#x201c;<article-title>Global scenario of sewage-sludge management</article-title>,&#x201d; in <source>Sustainable management and utilization of sewage sludge</source> (<publisher-loc>Berlin, Germany</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>383</fpage>&#x2013;<lpage>401</lpage>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nor&#x2019;Azlan</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Selamat</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Yahya</surname>
<given-names>N. M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Multivariable PID controller design tuning using bat algorithm for activated sludge process IOP Conference Series: materials science and engineering</article-title>. <source>IOP Conf. Ser. Mater. Sci. Eng.</source> <volume>342</volume> (<issue>1</issue>), <fpage>012030</fpage>. <pub-id pub-id-type="doi">10.1088/1757-899x/342/1/012030</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Dynamic multi-objective optimization control for wastewater treatment process</article-title>. <source>Neural Comput. Appl.</source> <volume>29</volume>, <fpage>1261</fpage>&#x2013;<lpage>1271</lpage>. <pub-id pub-id-type="doi">10.1007/s00521-016-2642-8</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reynoso-Meza</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Carre&#xf1;o-Alvarado</surname>
<given-names>E. P.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>November). PI tuning of a multivariable activated sludge process with nitrification and denitrification with multi-objective optimization</article-title>. <source>Proceedings</source> <volume>48</volume> (<issue>1</issue>), <fpage>4</fpage>. <pub-id pub-id-type="doi">10.3390/ECWS-4-06434</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Samuelsson</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Halvarsson</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Carlsson</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Interaction analysis and control structure selection in a wastewater treatment plant model</article-title>. <source>IEEE Trans. Control Syst. Technol.</source> <volume>13</volume> (<issue>6</issue>), <fpage>955</fpage>&#x2013;<lpage>964</lpage>. <pub-id pub-id-type="doi">10.1109/tcst.2005.854322</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Spreafico</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>An analysis of design strategies for circular economy through life cycle assessment</article-title>. <source>Environ. Monit. Assess.</source> <volume>194</volume> (<issue>3</issue>), <fpage>180</fpage>. <pub-id pub-id-type="doi">10.1007/s10661-022-09803-1</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sweetapple</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Butler</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Multi-objective optimisation of wastewater treatment plant control to reduce greenhouse gas emissions</article-title>. <source>Water Res.</source> <volume>55</volume>, <fpage>52</fpage>&#x2013;<lpage>62</lpage>. <pub-id pub-id-type="doi">10.1016/j.watres.2014.02.018</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tian</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A study on the multi-objective optimization method and characteristic analysis of installation locations of submersible mixer for sewage</article-title>. <source>Front. Energy Res.</source> <volume>11</volume>, <fpage>1180607</fpage>. <pub-id pub-id-type="doi">10.3389/fenrg.2023.1180607</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Vilanova</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Alfaro</surname>
<given-names>V.</given-names>
</name>
</person-group> &#x201c;<article-title>Multi-loop PI-based control strategies for the activated sludge process</article-title>,&#x201d; in <conf-name>Proceedings of the 2009 IEEE Conference on Emerging Technologies &#x26; Factory Automation</conf-name>, <conf-loc>Palma de Mallorca, Spain</conf-loc>, <conf-date>2009, September</conf-date>, <fpage>1</fpage>&#x2013;<lpage>8</lpage>.</citation>
</ref>
<ref id="B24">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Vinayak</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Narayanan</surname>
<given-names>S.</given-names>
</name>
</person-group> &#x201c;<article-title>Nonlinear PI with state feedback control for activated sludge process</article-title>,&#x201d; in <conf-name>Proceedings of the 2017 International Conference on Information, Communication, Instrumentation and Control (ICICIC)</conf-name>, <conf-loc>Indore, India</conf-loc>, <conf-date>2017, August</conf-date>, <fpage>1</fpage>&#x2013;<lpage>7</lpage>.</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Circular economy is game-changing municipal wastewater treatment technology towards energy and carbon neutrality</article-title>. <source>Chem. Eng. J.</source> <volume>429</volume>, <fpage>132114</fpage>. <pub-id pub-id-type="doi">10.1016/j.cej.2021.132114</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>