Plant-Wide Models for Optimizing the Operation and Maintenance of BTEX-Contaminated Wastewater Treatment and Reuse

: Benzene, toluene, ethylbenzene and xylenes, collectively known as BTEX compounds, are signiﬁcant emerging contaminants in municipal wastewater. Stricter efﬂuent quality regulations necessitate their removal, especially with concerns about organic micropollutant concentrations


Introduction
The consumption of raw materials, as well as emissions, tends to get concentrated within specific areas as human communities evolve and grow.Due to the continuous human impact on the environment, special attention must be paid to the operation of public works that are crucial for a sustainable modern civilization.Certain substances that may be hazardous to human health must not be allowed to reach consumers or to adversely affect environmental conditions.The pollutants building up throughout municipal sewer systems appear at water resource recovery facilities.
The practice of wastewater use in agriculture has become increasingly important due to the growing demand for water resources and the need for sustainable and efficient agricultural practices.The main reason is that it could help reduce the pressure on limited freshwater sources and thus mitigate the impacts of climate change by conserving the water quantity and reusing the nutrients from the treated wastewater.Low-income countries pose substantial risks by the direct use of untreated wastewater or the indirect use of polluted water courses [1].Scheierling et al. (2011) gave guidelines for developing a strategic plan on treated wastewater and highlighted the necessity of a paradigm shift for decision makers in wastewater treatment [2].Identifying the risk of being exposed to pathogens and heavy metals and the salinization of soil can be easily seen as drawbacks.Additionally, contaminants of emerging concern were given special consideration, taking into account the known or perceived environmental and health risks associated with these contaminants [3].Pedrero et al. (2010) found that with proper salts management, treated municipal wastewater appears to be a viable water resource for irrigating citrus trees, and no other concerns are related to treated municipal wastewater [4], whereas Jaramillo and Restrepo (2017) stated that the absence of a quantitative assessment of the microbiological risk, specifically related to the concentration of helminths, is the essential component that is necessary for the appropriate adoption of agricultural reuse [5].However, the utilization of untreated wastewater is prevalent in particular areas due to the unavailability of an appropriate alternative water supply [6].
Among the volatile organic compounds (VOCs) that are released into the environment from both natural and manufactured sources, benzene, toluene, ethylbenzene and xylene isomers (BTEX) are considered to be among the most significant [7].Due to their volatile nature, attention should be given to the health of the personnel working at and people living in the vicinity of treatment plants.To measure the concentration of BTEX compounds in the atmosphere surrounding wastewater treatment plants, air quality has been monitored in numerous studies.In an effort to assess the extent of exposure of wastewater treatment plant (WWTP) workers to BTEX compounds, Dehghani et al. (2022) conducted a study.The study revealed that the workers were indeed exposed to BTEX compounds, particularly benzene, in certain cases exceeding the permissible reference values.Thus, it is imperative to recognize WWTPs as a source of BTEX emissions, and workers employed in these facilities may be at risk of exposure to these volatile compounds [8].Numerous studies have investigated the potential adverse impacts of BTEX compounds on human health, including their effects on organs such as the nervous system, liver, heart and kidneys [9].
During biological wastewater treatment, BTEX pollutants are removed through various pathways such as biodegradation, sorption and volatilization.However, these pathways only offer partial reduction in comparison to the initial influent concentrations.For a substantial reduction in BTEX pollutants, further purification of the biologically treated water is necessary, primarily through physico-chemical methods.Advanced Oxidation Processes (AOPs) use highly reactive chemical species, such as hydroxyl radicals, to degrade or transform contaminants into less harmful compounds.AOPs have been shown to be effective in removing BTEX compounds from wastewater due to their ability to break down organic compounds into smaller, less toxic molecules [10,11].
Biofilm-based treatment processes have been employed in many cases; the microbial community-that attaches to a solid surface-is established and grown on a support material.During the process, the BTEX compounds diffuse into the biofilm, where they are metabolized by the microorganisms.Mello et al. (2019) achieved 99% removal efficiency by applying attached biomass-supported activated carbon media [12].The same reduction rate was achieved in horizontal-flow anaerobic immobilized biomass reactors where a denitrifying consortium was used as an inoculum to form a biofilm [13,14].Although advanced techniques such as the use of membranes [15], microalgae [16] or carbon nanotubes [17] have been developed for BTEX removal, the study did not aim to provide a comprehensive review of the state-of-the-art techniques in this area; the focus was more on the kinetics of BTEX biodegradation.
Fundamentally, the presence of dissolved oxygen (usually introduced through mechanical aeration) and chemically bound oxygen (in anoxic conditions) or the complete absence of oxygen (in anaerobic conditions) define the environmental conditions that support the growth and survival of biomass.In each environment, the BTEX removal kinetics were investigated; high biodegradation could be expected in aerobic conditions [18], while the buildup of toxic metabolites led to the incomplete degradation of both xylene and benzene, whether they were present as individual components or as a mixture.The addition of another pollutant always resulted in competition and inhibition, with xylene being the most inhibitory [19], whilst the external carbon source could cause enrichment on BTEX degradation by bacterial cultures [20].Anaerobic conditions could be developed in packed bed reactors or deeper layers of the biofilm, which could also show significant biodegradation [21].Lee et al. (1998) made significant steps towards modelling a highly hydrophobic compound in wastewater, and the authors explained that the biodegradation of a compound that has been adsorbed to activated sludge directly enhances its biodegradation and diminishes its discharge with the waste activated sludge.Additionally, the evaporation of compounds from the primary and secondary systems' surfaces is significant for compounds that possess moderate to high volatilities [22].Pomiès et al. (2012) gathered 18 mathematical approaches to modelling micropollutant removal from wastewater, and they concluded that most authors had held the view that only dissolved micropollutants are biodegradable, although no experimental evidence has been found to support this assertion.Simplified approaches based on first-order kinetics cannot take into account advanced metabolisms and the connection between the various stages of treatment (such as between mainstream and sludge line unit operations) or respond to process dynamics during operation.Incorporating particulate phase data into models would be extremely valuable in carrying out a comprehensive mass balance calculation, particularly for hydrophobic substances, which would include their fate in sludge.Furthermore, the recent incorporation of co-metabolism into models has revealed an interesting correlation between the biodegradation of macropollutants and micropollutants.These modelling tools are beneficial in conducting research to comprehend the mechanistic behavior of micropollutant removal, but they are not yet ready to be used in process design [23].
The objective of this paper is to develop numerical methodology for simulating the treatment and reuse of BTEX-contaminated sewage throughout a facility using dynamic mathematical models that can adapt to diurnal and seasonal variations and changes in operational approaches.The novelty of the study lies in its application of a modeling approach to a full-plant configuration, which investigates the potential for operational cost minimization.This is achieved by adjusting the sludge retention time or air flow to enhance BTEX degradation and stripping, respectively.
Modelling wastewater treatment processes is complex and involves several submodels, such as a hydrodynamic model, the characterization of wastewater components [24], a phase separation model [25], a gas dissolution model, a process control model and a biokinetic model.Biokinetic modelling is a widely used approach to describing the biological processes that occur in wastewater treatment systems.Activated Sludge Models (ASMs) are a type of biokinetic model that simulate the behavior of microorganisms in the activated sludge process, one of the most common wastewater treatment technologies [26].Plant-wide models, on the other hand, provide simulation capabilities of more advancedsuch as granular and biofilm-based-treatment processes, feature more biological processes such as anaerobic digestion, enable pH, precipitation and gas transfer calculations and allow for more complex plant configurations that may include mechanical pretreatment, sidestream treatment or even effluent wastewater recycling.They can further provide oversight of different aspects of facilities such as energy management and operational costs, and additionally, they take into account equipment operational limitations regarding industrial safety [27].
The International Water Association (IWA) published guidelines for good modelling practice in the field of wastewater treatment.These guidelines provide a framework for developing and applying mathematical models to treatment systems, aiming to ensure that models are used in a rigorous and consistent manner.The guidelines cover topics such as model development, verification and validation, sensitivity analysis and uncertainty analysis [28].
The aim of this study is to develop and present a methodology for robust calculations regarding the fate of BTEX compounds in biological wastewater treatment, which-through the use of plant-wide modelling-can be coupled with the concepts of physico-chemical treatment for effluent polishing and water reuse.The following chapters describe model development efforts that consisted of setting up a biokinetic model matrix and the introduction of gas transfer equations, all to be integrated into (and used alongside) existing concepts for plant-wide simulation.

Mathematical Model Development
The newly introduced sub-models were implemented according to a review of the scientific literature, equipment manufacturer specifications and experience gained throughout industrial modelling projects.Following a review of the related literature, it was decided that the most significant BTEX removal mechanisms-biochemical reactions and gas-liquid equilibrium-can be most efficiently coded in the Gujer matrix format [29], containing the stoichiometric coefficients and process rates, from which ordinary differential equations can be expanded and solved using computer algorithms.In this study, the rate expressions for BTEX modelling are added to the existing matrix of the commercially available MiniSumo model (see Tables A1-A13) They can be linked to any other similar biokinetic model; however, we note that the presented nomenclature regarding the state variables corresponds to the logic of MiniSumo.For the purpose of this study, this model was chosen as the basis for the sake of simulation speed and transparency, since it is a simple yet well calibrated plant-wide model that focuses on the most important and general aspects of water resource recovery-one-step nitrification and denitrification, anaerobic digestion, the oxygen requirement and sludge production [30].

Biodegradation
The biological reactions are interpreted using an ASM-type concept based mainly on the Monod-kinetics of saturation and inhibition, which, combined with a hydraulic sub-model (considering the amount of liquid entering individual compartments of unit processes), form volume-integrated scalar transport equations that define the mass balance throughout the plant configuration with simplified hydrodynamics.The component balance in Equation (1) describes the concentration change in a liquid phase component i (e.g., benzene) in time within a completely stirred tank unit with the liquid volume of V r .This ordinary differential equation (ODE) is generalized, and variables marked with the lower index must be expanded for each state variable in the model to be integrated over time by a solver algorithm.The liquid phase states generalized as L i may further be categorized based on the particle size, S i designating soluble components, C i applied for colloids and X i applied for particulates as follows: The mass flow or load to a tank is calculated knowing the influent concentration of the state variable-L i,in -entering the tank, by Equation (2).
In a similar manner, the component mass flow leaving the tank is interpreted by Equation (3).
The kinetic reaction rate of the state variable-interpreted in the dimension of concentration change-expressed as rate i in Equation (4), is derived from the kinetic matrix by multiplying a stoichiometric coefficient (v j,i , regarding the interaction of the given process j with the state variable i) by process rates (r j for the process j), summing them up for every case where a given process affects the state variable (based on reading the matrix) [29].
Equation ( 5) determines rateF i , the mass rate of the process variable, which was seen earlier (regarding the liquid phase) as contributing to the component balance from Equation (1).
The results of the research conducted by Kasi et al. (2012) were used for the concept of BTEX degradation kinetics [20], to be implemented in Gujer matrix format for this study.Appendix A contains the final format of the developed matrix (Tables A1-A10).It is assumed that one common group of heterotrophic biomass (present in all ASM-type models) metabolizes the BTEX contaminants, either oxidizing them aerobically using oxygen or removing them in an anoxic environment with only nitrite or nitrate available as the electron acceptor; otherwise, they are degraded through fermentation when neither of the previously mentioned electron acceptors are available.The reaction rate for aerobic pathways is higher than those of non-aerated conditions.
In this study, the process model MiniSumo was modified to include the presented BTEX-related biological reactions and gas exchange processes.MiniSumo is a full-plant model developed for simulating oxygen uptake rates and sludge production.It considers the aerobic, anoxic and anaerobic removal of soluble and particulate biodegradable organics, nitrification, denitrification, methanogenesis and phosphorus removal through chemical addition [30].Regarding the model extension, BTEX-related ordinary heterotrophic organism (OHO) growth processes have to be introduced into a model that already contains the state variables oxygen, nitrite and nitrate (either separately or combined), as well as volatile fatty acids.The latter variable, for instance, is not originally featured in MiniSumo; therefore, it was necessary to include VFA production and uptake (along with the related missing biochemical processes, such as the fermentation of the readily biodegradable substrate) from the matrix of the Sumo1 process model [31] in our example.
Eight new state variables are included in the model matrix to represent both dissolved and gaseous states of the BTEX compounds in g COD m −3 concentration units.These include dissolved benzene (S BENE ), dissolved toluene (S TENE ), dissolved ethylbenzene (S EBENE ), dissolved xylene (S XENE ), gaseous benzene (G BENE ), gaseous toluene (G TENE ), gaseous ethylbenzene (G EBENE ) and gaseous xylene (G XENE ).All the dissolved BTEX compounds undergo biodegradation under anaerobic, anoxic and aerobic conditions carried out by Ordinary Heterotrophic Organisms (X OHO ) biomass.All the biological reactions are parallel reactions that completely mineralize the BTEX compounds to carbon dioxide and biomass.A complete list of rates and the model structure are disclosed as supplementary information in Appendix A (Tables A1-A10).
The temperature-dependent growth rates of the added metabolic processes can be calculated using Arrhenius-type functions, with this example study applying the MiniSumo implementation for heterotrophs.The kinetic rates and stochiometric yields associated with the biological reactions are described in Appendix B (Table A11), relying on the experimental data of Kasi et al. (2012) [20].

Gas-Liquid Transfer
It is important to discuss the fate of BTEX components emitted from the wastewater into the atmosphere.The equilibrium between the vapor and liquid phases for these volatile compounds is simulated using Fick's two-film theory, analogous to existing gas transfer processes in MiniSumo, serving as an example of model development for absorption and desorption.The impact of equilibrium processes on state variables is described within the process model's Gujer matrix.We note that, within the gas transfer equations, standard conditions regarding gases are interpreted as 101,325 Pa for pressure (p NTP ) and 293.15K for temperature (T NTP,K ).
The internal gas phase components are simulated in units of mass per liquid volume.This makes mass and component balancing substantially simpler (compared to using the gas volume as the basis).Equation ( 6) describes the component balance for gas i.
F Gi,air,inp denotes the mass flow of gas i from the air supply and is calculated based on the composition of the input gas and the air flow based on the ideal gas law, as shown by Equation (7).Appendix C elaborates on the stoichiometry behind the equivalent molar mass MM EQ,i of each gaseous BTEX component derived from the theoretical oxygen demand.
For gas phase mass balancing, the output gas mass flow F Gi,air,outp is calculated by Equation (8), where the concentrations expressed per liquid volume must be converted per volume of gas.
The volume of the gas phase in the reactors is estimated based on the gas hold-up fraction ε gas .The gas hold-up is an input parameter, for which this paper uses a value of 0.01 m 3 gas m −3 in case of aerated conditions and 0.001 m 3 gas m −3 regarding non-aerated conditions [32].The gas volumes at field conditions and standard conditions (V gas and V gas,NTP , respectively) are calculated by Equations ( 9) and (10).
V gas,NTP = The gas phase pressure p gas is determined by Equation (11) in accordance with EPA guidelines on dissolved oxygen saturation calculations [33].
The vapor pressure p v,T is calculated by Equation ( 12) based on the water temperature and Antoine-coefficients, where a multiplication factor is used for unit conversion from mmHg to Pa.
The air pressure p air is calculated by the barometric Equation ( 13), considering the elevation of the facility above sea level h sea (the molar mass of air MM air is converted from the unit g mol −1 to kg mol −1 ).The parameter L air for the temperature lapse rate applies a value of 0.0065 K m −1 [34].
The calculation of the effective saturation depth depends on whether the tank is aerated, in which case h sat,eff is calculated by Equation ( 14), and in case there is no aeration, the diffuser submergence h diff is replaced by the total side water depth (h r ).The saturation depth fraction f h,sat,eff may be assigned according to various approaches; the mid-depth concept provides reasonable accuracy in diffused aeration [35].
For the diffuser depth calculation, the diffuser mounting height is subtracted from the liquid height, according to Equation (15).
Unlike liquid phase hydraulic balancing, where the effluent flow is assumed to be equal to the influent, gas phase volumetric balancing needs to take into account the input air flow, as well as the gas transfer activity in the volumetric unit.Equation ( 16) determines the outgoing gas volumetric flow, Q gas,outp,NTP , which, physically, cannot be negative; this shall be reflected in the model code (e.g., using a maximum function).
The volumetric flow of the gas transfer Q gas,transfer,NTP is derived from the summation of the individual gas component transfer molar flows, calculated by Equation (17).
Two gas transfer routes are modelled: between the atmosphere and water surface and between the bulk-of-the liquid and bubbles.According to Fick's first law, the process rate involving the transfer of gas i between the atmosphere and water surface is described by Equation (18), with the driving force of either dissolution or stripping is determined by the difference in the saturation concentration and the liquid phase concentration.The transfer rate of exchange at the gas bubble interface is fundamentally analogous, described by Equation (19).These equations are listed for processes involving the individual gas phase BTEX states in the Gujer matrix, located within Appendix A.
Regarding gas transfer through bubbles, the volumetric liquid-side mass transfer coefficients k L a i,bub,st,cw -in clean water for a standard 20 • C temperature-and k L a i,bubin wastewater at the field temperature-are calculated by Equations ( 20) and (21).
The interfacial transfer area a bub is calculated by Equation ( 22) on a geometrical basis.Based on the literature review, in this study, the bubble Sauter mean diameter d bub is input as 0.003 m for aerated conditions and as 0.01 m under non-aerated conditions [32,36,37].
It shall be noted that in cases of fully covered tanks, atmospheric gas transfer is not modelled; the specific surface area a bub is extended with the specific contact area between the water surface and the headspace, as shown by Equation (23).
On the basis of the Higbie penetration theory, values of the liquid-side mass transfer coefficient k L,i,bub,st,cw are calculated using the respective diffusivities of volatile and soluble gas components by Equation (24) [38], relying on k L,O2,bub,st,cw as a model parameter with the input value of 0.54 m h −1 [39].BTEX-specific diffusion coefficients were averaged from values found in the scientific literature [40][41][42] and are disclosed in Appendix C (Table A12) The parameter for the fraction in the liquid side f kL,i equals 1 for all gases in this study, as they diffuse slowly through the liquid film, with the exception of 0.05 in the case of ammonia, due to its very high solubility [43].
Throughout aerated compartments, the α correction factor is modelled dynamically from process variables representative of wastewater characteristics and loads, relying on a method published for increased accuracy in oxygen transfer and air flow requirement predictions [44].The θ compensation factor for simulating Arrhenius-type temperature sensitivity in mass transfer is taken as 1.024, following standard procedures in wastewater treatment design [33].
For aerated conditions, regarding oxygen, k L a O2,bub is calculated by Equation ( 25), incorporating a further correction factor for diffuser fouling, an input parameter of aerated unit processes [45], which, in practice, is best adjusted knowing the diffuser age and the time of the last cleaning procedure [46].For non-aerated conditions (and regarding coarse bubbles), the correction factor for fouling is not interpreted (equal to one).
Regarding oxygen, under aerated conditions, k L a O2,bub,st,cw directly relates to the clean water performance of diffusers; thus, it is derived-according to Equation ( 26)-from the standard oxygen transfer rate involving air bubbles, SOTR bub .
SOTR bub is derived from the specific standard oxygen transfer efficiency SSOTE, which describes the diffused aerator characteristics in clean water, as Equation ( 27) shows.
The prediction of SSOTE is implemented according to a novel concept proposed in the form of Equation ( 28), as a function of the diffuser-specific air flux Q air,NTP,sp [47].
This proposed concept simulates the phenomena concerning fine bubble aeratorswhere SSOTE is known to drop significantly with an increased air flow rate [48]-and regarding coarse bubble diffusers-which show an increasing trend of SSOTE with an increasing air flow [49].The specific air flow per diffuser, Q air,NTP,sp is calculated by Equation (29).
SSOTE is corrected to account for the diffuser submergence and density.In the case of coarse bubble aeration, the correction for the diffuser density has no practical meaning, and thus, it is neglected; however, in fine bubble aeration, increasing the diffuser floor coverage will raise the SSOTE [50].The correction for diffuser depth is implemented to reflect the phenomenon that the absolute transfer efficiency (SOTE) does not increase linearly with depth, especially in tanks deeper than 8 m [51].
The correction term for diffuser density, corr d,diff , is calculated by Equation (30).
The correction term for diffuser submergence, corr h,diff , is calculated by Equation (31).
The diffuser floor density d diff is the ratio of the total diffuser area to the tank surface, as Equation ( 32) shows.
The method for clean water aeration efficiency estimation relies on an exponential function with asymptotes, instead of setting strict Boolean-type threshold levels for SSOTE, so that the air flux may safely be extrapolated outside of the model calibration range.This avoids potential numerical issues when utilized in a dynamic simulator, with the further advantage of not carrying over conceptual errors in SOTE, when combining it with factors describing the process water performance.The model parameters were fine-tuned based on data collected from numerous diffuser tests for various equipment types [33,52].Appendix C contains the parameter set used for this paper, characterizing ceramic discs.
At the atmospheric interface, calculating the volumetric liquid-side mass transfer coefficients involves the specific mass transfer coefficients (k L ) estimated for the liquid surface.The clean water k L a i,sur,st,cw at standard conditions and the process water k L a i,sur for field conditions are calculated by Equations ( 33) and (34), respectively.
The liquid-side mass transfer coefficient regarding the water surface k L,i,sur is determined by Equation ( 35) using the same principle as for the gas bubble interface, as previously shown by Equation (24).
The specific surface area of the liquid a sur is calculated by Equation ( 36), incorporating the multiplication factors f cover for reactor coverage and f wave for turbulence (waviness).The model inputs of 0.54 m h −1 for k L,O2,sur and 1.9 for f wave in this study are adjusted based on the typical measured values of k L a O2,sur found in the relevant literature [53].
The water surface area A r is defined by basin geometry, according to Equation (37).
According to Fick's law, as illustrated previously by Equations ( 18) and ( 19), the saturation concentration of a component determines whether the gas-liquid transfer process is driven towards dissolution or stripping.They are calculated by Henry's law based on the partial pressure of a component, with temperature dependency correction implemented using the van't Hoff equation [54].The saturation concentrations of gases attributed to the liquid-gas bubble interface S i,bub,sat,st,cw , standardized at 20 • C for clean water (thus requiring temperature conversion from 25 • C, used as the basis of Henry's law), are modelled based on Equation (38).The variable for process water at field conditions, S i,bub,sat,st , is expressed by Equation (39).The values of the Henry's law model parameters regarding BTEX contaminants were averaged in the process of the literature review [55-57] and are listed in Appendix C (Table A12).
The β correction factor for impurities in process water uses an input parameter value of 0.95, which is typical of municipal wastewater [33].
With regard to the bubble-water interface, partial pressures for standard conditions, p partial,I,bub,st , and for process conditions, p partial,i,bub , are calculated by Equations ( 40) and ( 41), respectively.
p partial,i,bub = G i /MM EQ,i n gas,bub • p gas (41) Equation ( 42) calculates the molar quantity of gas bubbles per unit liquid volume, n gas,bub .
As explained by Equation (43), the pressure at standard conditions, p st,h,sat,eff , is also compensated for the effective saturation depth, as diffuser testing in clean water involves the design submergence [58].
The off-gas composition in molar or volume percent units (interchangeable according to Avogadro's law) may be calculated using Equation (44), expressing G i,percent , derived from the individual gas component molar concentrations and their sum.
The saturation concentrations for the atmosphere-liquid interface, S i,sur,sat,st,cw in the case of clean water and S i,sur,sat regarding process water, are calculated by Equations ( 45) and ( 46). S i,sur,sat,st,cw = p partial,i,sur,st • Henry i,25 • e Henry i,dt •((1/T NTP,K )−(1/T SATP,K )) • MM EQ,i , ( 45) For the atmospheric saturation concentration calculations, the composition of the atmosphere around the water surface is required.The compositional figures G i,atm are defined as model constants, measured in volumetric percentages.Equations ( 47) and ( 48) quantify the partial pressure of gases in the atmosphere, p partial,i,sur,st for standardized conditions and p partial,i,sur for field conditions.
To sum up, the gas phase composition and the amount of gas dissolving or stripping are dependent on individual gas state variables that impact one another, according to the laws of mass balance.In accordance, the sufficient injection of air or nitrogen gas production will result in BTEX components stripping out of the liquid.

Adsorption on Granular Activated Carbon
Describing the simulated mechanism of BTEX contaminants binding to activated carbon granules is essential.The developed methodology in this study features a granular activated carbon sub-model with an engineer-oriented concept; the GAC media consumption is estimated according to the bed volume and TOC adsorption capacity specified.
With the goal of shortcut steady-state runs, the GAC unit operation process model may be set up using a continuous process flow approach, purely based on algebraic equations, and focusing on the plain instantaneous mass-balances of states, with the effluent quality estimated for operational phases with loading.For more robust dynamic modelling aimed at evaluating maintenance or control strategies, a cyclical process flow model setup is implemented.The latter features asymmetric breakthrough curves for quantifying the effluent quality based on the accumulated pollutant load and requires an ODE solver for the integrated mass balance calculation.
Water treatment using GAC adsorption is applicable to dissolved organic components.The adsorbable states in the extended MiniSumo biokinetic model presented in this chapter include VFA, the liquid phase BTEX variables, soluble unbiodegradable and readily biodegradable substrates as well as soluble biodegradable organic N. The unit process implementation does not include a reactive volume, with the assumption being that no biodegradation of adsorbed COD occurs.This approach was selected as appropriately sufficient for the main target of this study, which applies GAC as a tertiary treatment step, after readily biodegradable and slowly biodegradable substrates have mostly been removed biologically and mechanically, respectively.
It is assumed that, per unit mass of AC, a specified amount of soluble organic pollutants-expressed in unit mass of C-can be adsorbed, until the medium reaches its breakthrough capacity.Thus, in order to quantify the number of adsorbed components, their respective mass concentrations have to be converted to carbon equivalent units.This involves stoichiometric conversions from unit g COD into unit g C regarding the carbonaceous state variables and from unit g N into g C in case of organic nitrogen.The COD of organic substrates are interpreted using the concept of theoretical oxygen demand.Appendix D (Table A13) lists the ratios of COD (and N) to TOC.The continuous process flow model instantaneously calculates the carbon equivalent for all adsorbed components on the GAC bed volume based on the influent concentrations, using Equation (49).The removal efficiency Rem GAC,i was assigned to be 99% regarding BTEX components, 92% for a soluble unbiodegradable substrate and 90% in the case of VFA, as well as soluble biodegradable carbon and nitrogen [59,60].
As shown by Equation ( 50), from the summed carbon equivalent variable EQ C,ad,total and the influent flow, the algebraic model can calculate the average replacement cycle frequency, knowing the set bed volume V ac and the TOC breakthrough capacity BTC.
Using Equation ( 51), BTC is converted from the breakpoint-related TOC adsorption capacity BTC m in the mass fraction unit and the apparent density of the media ρ ac , both of which are design settings of the GAC model.
The required mass flow of the replaced GAC media can be quantified by Equation ( 52), according to the mass of activated carbon filled per cycle M ac,cycle .
The mass of GAC filled into the towers within each cycle corresponds to the bed volume, as Equation ( 53) explains.
The replacement frequency can also be fixed based on the lifespan of the AC media.The cyclical process flow model schedules replacement events either based on the lifespan or relying on breakthrough-based or volumetric load-based maintenance logic.In the case of the breakthrough-based approach, either the amount of the total adsorbed carbon per bed volume EQ C,ad,total is monitored and compared to the BTC or the effluent TOC is monitored and compared to a TOC concentration threshold to trigger the bed replacement procedure.Regarding the volumetric load-based logic, the influent loading expressed in the number of bed volumes is registered and initiates the replacement process when it reaches the target volume ratio of the load per bed.
In the dynamic GAC model, EQ C,ad,total is expressed directly as a sum of the accumulated C equivalent mass of pollutants per bed volume, following Equation (54).
L i,ad, , the adsorbed mass per bed volume of state variable i, is determined by differential equations that have to be solved by a numerical integrator.The derivative of L i,ad is expressed by Equation (55), describing GAC operational loading phases.The authors note that the same equation can also be used to account for potential impurities from backwash streams.
dL i,ad dt During bed replacement phases, the adsorbed state variable masses per bed volume are nullified, so the change in L i,ad is described simply using a zero-order rate expression, as seen in Equation (56).The accumulated mass per bed volume prior to the replacement L i,ad,pre shall be registered at the start of the replacement event, which is linearly decreased towards 0 throughout the duration of the bed replacement t repl .
The actual removal efficiency of organic constituents, Rem GAC,i,actual , depends on how saturated the granular medium is; based on EQ C,ad,total , the accumulated adsorbed components, a logistic saturation function expresses the drop compared to the defined initial removal ratio Rem GAC,i , as Equation ( 57) shows.
Applying the logistic formula would, in itself, generate a symmetric S-shaped curve; therefore, as a function of EQ C,ad,total , the midpoint concentration C mid is corrected for the asymmetry involving breakthrough curves, as pointed out by Equation (58).
C mid,symm , the midpoint concentration of the breakthrough curve without asymmetry correction, is derived from the BTC, the breakpoint fraction f break and the breakthrough curve slope sl break , using Equation (59).The drop in removal efficiency at the breakthrough point adsorption capacity is defined by f break .The slope and the two model parameters of the asymmetry correction term pow mid,asymm and magn mid,asymm were adjusted based on measured data involving breakthrough curves found in the relevant literature, and they can be re-estimated depending on the type of GAC media applied for the treatment [61,62].The values of the parameters are disclosed in Appendix D (Table A13).
The actual removal efficiency calculated using the asymmetric sigmoid expression also serves as the basis of estimating effluent state variable concentrations, as shown by Equation (60).The model with continuous process flow implementation does not consider the saturated status of the bed, so it uses the ideal removal efficiency ratio input instead.It is important to note that this equation only applies to soluble organic states; other dissolved compounds (such as ammonia and orthophosphate) do not have removal efficiencies specified and will pass through the GAC unit without treatment.
Although in practice, TSS loading to the GAC process should be minimal, activated carbon beds can also act as filtration media.The developed model includes an instantaneous mass balance separator approach for the mechanical separation of particulate state variables.This was implemented following the principle applied for the sand filter unit process in the simulation software used for demonstrations in this study, Sumo22 [30].Regarding the cyclical process flow model setup, the deposition of suspended matter on the bed volume is considered in the mass balance, with derivatives calculated analogously to Equation (55).Backwash events can also be scheduled by specifying the loading phase duration.Similar to the concept of wastage during replacement events quantified by Equation ( 56), the decrease in the captured particulate mass per volume due to backwashing is determined using simplified zero-order expressions, relating to the backwash duration.

Model Configuration
For the demonstration of BTEX modelling applications, a variant of the BSM1 test configuration [63] was set up in the Sumo22 simulation software and used to evaluate potential operational cost savings through controlling the sludge retention time or by air stripping and through comparing strategies of GAC tower maintenance.The chosen technological sequence is a simulation benchmark that defines a relatively simple standard water resource recovery facility with influent loads provided, aimed at evaluating system performance based on test procedures, according to a given set of evaluation criteria.The BSM1 plant comprises a five-compartment activated sludge reactor cascade, including two anoxic zones and three aerobic zones.This setup combines nitrification with predenitrification, a widely adopted arrangement for achieving biological nitrogen removal.Downstream of the reactor cascade, the configuration features a secondary clarifier for phase separation [63].Table 1 summarizes the main characteristics of the benchmarking facility, regarding average dry weather conditions.In addition to the BSM1 layout, Sumo22's empirical blower and pump models were set up to assess energy and associated costs regarding aeration, return activated sludge pumping, nitrified mixed liquor recycle pumping as well as backwash pumping.Solving the rate expression and component balance equations of variables in wastewater treatment technologies requires a simulation environment based on biokinetics.From an operational strategic point of view, process modelling can reveal cost-down opportunities by experimenting with control strategies and pointing out necessary interventions for complying with environmental regulations.Sumo22 is a multipurpose, open-source environment that has been developed specifically for dynamic environmental modelling, with particular emphasis on municipal and industrial wastewater treatment plants.Sumo22 can simulate a wide variety of BNR (biological nutrient removal) plant configurations [30].
According to the model layout prepared for the study, portrayed by Figure 1, as a fictive retrofit of the municipal plant with effluent purification, the scope of the standard BSM1 configuration was extended with a granular activated carbon process and a microfiltration step prior to GAC for the sake of mechanical protection against the remaining solids overflowing from the secondary clarifiers.In the granular activated carbon unit process model, the replacement frequency of the exhausted bed is determined based on the TOC adsorption capacity of the media and the load of soluble organics removed.It can be run in either a quasi steady-state mode, where the AC demand is quantified from an instantaneous carbon mass balance, or in dynamic mode, where it is calculated cyclically, based on the effluent quality predicted using an empirical asymmetric breakthrough curve.The total bed volume of the towers was sized as 50 m 3 ; with the effluent quality and quantity associated with the BSM1 use case, this ought to require bed replacement roughly every month, assuming that the bed density is 450 kg m −3 and that the soluble organic carbon adsorption capacity of the bed is 0.2 g C g AC −1 , estimated based on operational expertise.The MF and GAC units are backwashed in half-hour and daily cycles, respectively, for industrial safety reasons against solids depositing on the filter and the granular media.We note here that bed replacement can be calculated automatically in Sumo22 from the mass flow of the BTEX compounds removed, simulating how-in a full-scale treatment plant-the operators would remove exhausted granules and reinstall fresh granules periodically, whenever the breakpoint of the media is reached.After sedimentation, an additional re-aeration zone with the hydraulic retention time of 0.3 h-with coarse bubble aeration-is used to maintain sufficient dissolved oxygen levels required for the effluent quality.The influent quality regarding BTEX concentrations was characterized according to typical levels in municipal wastewater, containing 303 µg L −1 of benzene, 290 µg L −1 of toluene, 249 µg L −1 of ethylbenzene and 933 µg L −1 of xylene [64].Using the scenario analysis functionality of Sumo22, the targeted sludge retention time and the air flow to the re-aeration tanks were selected as the input parameters to be changed in between the simulation runs in order to efficiently evaluate the changes in the blower energy demand, the yearly requirement of GAC and the amount of BTEX chemicals removed.SRT control is handled according to the BSM1 scenario guidelines; Sumo22 automatically adjusts the wasted sludge flow rate, varying the activated sludge solids concentration to meet the SRT target.

Results and Discussion
The simulation and results interpretation involved preparing the model output setup for exporting variables linked to pollutant loads, composition profiles, effluent quality, biodegradation, stripping and adsorption rates as well as calculations on resources, regarding aeration requirement and activated carbon consumption.Sumo22's combined global and local solvers were used to handle the integrated calculations necessary for plant-wide modelling.Based on a steady-state simulation featuring a BSM1 baseline scenario, Figure 2 displays a profile of the summed BTEX component uptake and stripping rates throughout the biological treatment, expressed in the mass flow of COD per reactor cell volume.The BTEX biodegradation and stripping rate profiles represent the plug-flow characteristics of the treatment plant; degradation is most significant at the first anoxic cell because, here, the BTEX compounds serve as substrates for the biomass in the highest concentration.There is a decline in biological uptake in the second anoxic zone as substrate availability becomes more limited, but as the wastewater enters the first aerated compart- Using the scenario analysis functionality of Sumo22, the targeted sludge retention time and the air flow to the re-aeration tanks were selected as the input parameters to be changed in between the simulation runs in order to efficiently evaluate the changes in the blower energy demand, the yearly requirement of GAC and the amount of BTEX chemicals removed.SRT control is handled according to the BSM1 scenario guidelines; Sumo22 automatically adjusts the wasted sludge flow rate, varying the activated sludge solids concentration to meet the SRT target.

Results and Discussion
The simulation and results interpretation involved preparing the model output setup for exporting variables linked to pollutant loads, composition profiles, effluent quality, biodegradation, stripping and adsorption rates as well as calculations on resources, regarding aeration requirement and activated carbon consumption.Sumo22's combined global and local solvers were used to handle the integrated calculations necessary for plant-wide modelling.Based on a steady-state simulation featuring a BSM1 baseline scenario, Figure 2 displays a profile of the summed BTEX component uptake and stripping rates throughout the biological treatment, expressed in the mass flow of COD per reactor cell volume.
The BTEX biodegradation and stripping rate profiles represent the plug-flow characteristics of the treatment plant; degradation is most significant at the first anoxic cell because, here, the BTEX compounds serve as substrates for the biomass in the highest concentration.There is a decline in biological uptake in the second anoxic zone as substrate availability becomes more limited, but as the wastewater enters the first aerated compartment, the degradation rate is boosted due to the faster aerobic reactions.The elimination of BTEX compounds through gas exchange is mainly visible at aerobic zones due to the contribution of air stripping driving out greater mass flows of volatile contaminants; stripping from anoxic compartments is mostly limited to the interface of the liquid surface and the atmosphere.Throughout the aerated cells, the profiles of removal rates attributed to both biodegradation and stripping correspond well with the decreasing concentration gradient of BTEX pollutants during activated sludge treatment.
biodegradation, stripping and adsorption rates as well as calculations on resources, regarding aeration requirement and activated carbon consumption.Sumo22's combined global and local solvers were used to handle the integrated calculations necessary for plant-wide modelling.Based on a steady-state simulation featuring a BSM1 baseline scenario, Figure 2 displays a profile of the summed BTEX component uptake and stripping rates throughout the biological treatment, expressed in the mass flow of COD per reactor cell volume.The BTEX biodegradation and stripping rate profiles represent the plug-flow characteristics of the treatment plant; degradation is most significant at the first anoxic cell because, here, the BTEX compounds serve as substrates for the biomass in the highest concentration.There is a decline in biological uptake in the second anoxic zone as substrate availability becomes more limited, but as the wastewater enters the first aerated compartment, the degradation rate is boosted due to the faster aerobic reactions.The elimination of BTEX compounds through gas exchange is mainly visible at aerobic zones due to the In this section on demonstrated studies, the analysis is primarily centered around examining the model concept in detail based on SRT, followed by an evaluation of the effect of air stripping and concluding with a presentation of GAC operational strategy modelling.

SRT-Based Scenarios
In the first modelling experiment, the test configuration ran at a steady state in dry weather conditions; using Sumo22's combined global and local solvers, a series of runs were performed with varied SRT to obtain the blower power consumption and the required mass flow of the GAC bed to be replaced.From these variables and the raw influent characteristics regarding BSM1 specifications, the blower energy and activated carbon demand specific to the influent BOD were plotted, as shown in Figure 3, to provide more generic results irrespective of the benchmarking plant capacity.
Environments 2024, 11, x FOR PEER REVIEW 16 of 34 contribution of air stripping driving out greater mass flows of volatile contaminants; stripping from anoxic compartments is mostly limited to the interface of the liquid surface and the atmosphere.Throughout the aerated cells, the profiles of removal rates attributed to both biodegradation and stripping correspond well with the decreasing concentration gradient of BTEX pollutants during activated sludge treatment.In this section on demonstrated studies, the analysis is primarily centered around examining the model concept in detail based on SRT, followed by an evaluation of the effect of air stripping and concluding with a presentation of GAC operational strategy modelling.

SRT-Based Scenarios
In the first modelling experiment, the test configuration ran at a steady state in dry weather conditions; using Sumo22's combined global and local solvers, a series of runs were performed with varied SRT to obtain the blower power consumption and the required mass flow of the GAC bed to be replaced.From these variables and the raw influent characteristics regarding BSM1 specifications, the blower energy and activated carbon demand specific to the influent BOD were plotted, as shown in Figure 3, to provide more generic results irrespective of the benchmarking plant capacity.The results show that the reduction in demand for the GAC replacement is most significant within the range of 3 to 6 days of SRT.However, this range would not be typical for a water reclamation plant, as even at warm/moderate water temperatures, longer retention times are needed for proper nutrient removal.It typically takes between 10 and 15 days for a thorough nitrification-denitrification process.If the SRT falls below this range, the process remains incomplete.Conversely, if the SRT exceeds this range, the plant's energy efficiency is compromised [65].After 6 days, the further decrease in the activated carbon demand is not significant, given that, proportionally, a much higher performance is demanded from the blowers to meet the dissolved oxygen concentration target at higher biomass concentrations.The reason behind this is that most of the organic load treated by the towers is of the soluble unbiodegradable organic fraction, the decomposition of which cannot be boosted by raising the sludge residence time.
Apart from the operating cost, another important aspect of micropollutant removal is the effluent quality.The BTEX concentrations in the treated water were also plotted The results show that the reduction in demand for the GAC replacement is most significant within the range of 3 to 6 days of SRT.However, this range would not be typical for a water reclamation plant, as even at warm/moderate water temperatures, longer retention times are needed for proper nutrient removal.It typically takes between 10 and 15 days for a thorough nitrification-denitrification process.If the SRT falls below this range, the process remains incomplete.Conversely, if the SRT exceeds this range, the plant's energy efficiency is compromised [65].After 6 days, the further decrease in the activated carbon demand is not significant, given that, proportionally, a much higher performance is demanded from the blowers to meet the dissolved oxygen concentration target at higher biomass concentrations.The reason behind this is that most of the organic load treated by the towers is of the soluble unbiodegradable organic fraction, the decomposition of which cannot be boosted by raising the sludge residence time.
Apart from the operating cost, another important aspect of micropollutant removal is the effluent quality.The BTEX concentrations in the treated water were also plotted against the same range of SRT, as displayed in Figure 4.The difference in BTEX removal is most significant between the SRT range of 3 to 6 days, as a greater load of contaminants, which were not exposed to the biomass for a sufficient amount of time to be treated, bleed through the GAC bed.Consequently, it is important to note that the activated sludge system must always be operated responsibly, not allowing the SRT to drop below the necessary level for nutrient removal, as the activated carbon treatment step might not adsorb the surplus in the load sufficiently to meet regulations, depending on the use case (e.g., reuse for agriculture, drinking water supply).Benstoem et al. [61] noted a decrease in micropollutant removal efficiency due to the improper functioning of a secondary clarifier, resulting in an increase in the average TSS concentration from 6 mg/L to a peak of 30 mg/L [66].
Environments 2024, 11, x FOR PEER REVIEW 17 of 34 tem must always be operated responsibly, not allowing the SRT to drop below the necessary level for nutrient removal, as the activated carbon treatment step might not adsorb the surplus in the load sufficiently to meet regulations, depending on the use case (e.g., reuse for agriculture, drinking water supply).Benstoem et al. [61] noted a decrease in micropollutant removal efficiency due to the improper functioning of a secondary clarifier, resulting in an increase in the average TSS concentration from 6 mg/L to a peak of 30 mg/L [66].For interpreting the results of the modelling scenarios in more detail, the contributors to the BTEX metabolism throughout the treatment plant are shown in Figure 5, based on comparing the total mass flow of the compounds biodegraded, stripped and adsorbed to the influent BTEX load.The mass flows are interpreted in equivalent units of chemical oxygen demand.When discussing the COD percent removal of BTEX pollutants, according to the modelling study, biodegradation is predicted to be the most significant process-by raising the plant-wide SRT, a steep increase can be achieved in the extent of anoxic and anaerobic biodegradation.The increased biodegradation is in accordance with decreased stripping and the reduction in the adsorbed mass flows on the activated carbon beds since the remaining micropollutant load, to be removed by gas transfer or effluent polishing, drops.However, as deduced from the results in Figure 4, most of the load fed to GAC treat- For interpreting the results of the modelling scenarios in more detail, the contributors to the BTEX metabolism throughout the treatment plant are shown in Figure 5, based on comparing the total mass flow of the compounds biodegraded, stripped and adsorbed to the influent BTEX load.The mass flows are interpreted in equivalent units of chemical oxygen demand.When discussing the COD percent removal of BTEX pollutants, according to the modelling study, biodegradation is predicted to be the most significant process-by raising the plant-wide SRT, a steep increase can be achieved in the extent of anoxic and anaerobic biodegradation.The increased biodegradation is in accordance with decreased stripping and the reduction in the adsorbed mass flows on the activated carbon beds since the remaining micropollutant load, to be removed by gas transfer or effluent polishing, drops.
However, as deduced from the results in Figure 4, most of the load fed to GAC treatment is in the form of soluble unbiodegradable compounds; therefore, there is only limited cost down capability in minimizing the bed replacement frequency through intensifying the activated sludge process.It can be concluded that, the more complete the nutrient removal is, the better the effluent quality is, including lower BTEX concentrations.However, the increased demand for aeration is more significant than the amount of activated carbon saved; thus, an actual SRT target is practically determined by how stringent the water quality requirements are.However, ultra-short SRT processes (<4 days) aiming for carbon removal from wastewater demonstrate high efficiency in organic matter separation and in enhancing energy recovery rates in wastewater treatment plants [67].
oxygen demand.When discussing the COD percent removal of BTEX pollutants, according to the modelling study, biodegradation is predicted to be the most significant process-by raising the plant-wide SRT, a steep increase can be achieved in the extent of anoxic and anaerobic biodegradation.The increased biodegradation is in accordance with decreased stripping and the reduction in the adsorbed mass flows on the activated carbon beds since the remaining micropollutant load, to be removed by gas transfer or effluent polishing, drops.However, as deduced from the results in Figure 4, most of the load fed to GAC treatment is in the form of soluble unbiodegradable compounds; therefore, there is only lim-

The Effect of Aeration Intensity
In another demonstrated case study of the developed numerical methods, the scenario analysis was conducted with varied aeration intensities of the facility's post-aeration tanks for the different steady-state runs.The decision to modify air flow to this re-aeration process-rather than changing DO setpoints in the activated sludge trains-came down to the fact that oxygen transfer will be more effective in secondary clarifier effluent water than in mixed liquor due to the absence of biomass oxygen uptake and higher α factors.The volume to be aerated is also much smaller this way, with sufficient HRT for the reaction rate of stripping.Figure 6 presents these scenarios' results on BOD-specific aeration energy and activated carbon demand as a function of superficial gas velocity-the air flow expressed per re-aeration tank surface area.
Environments 2024, 11, x FOR PEER REVIEW 18 of 34 ited cost down capability in minimizing the bed replacement frequency through intensifying the activated sludge process.It can be concluded that, the more complete the nutrient removal is, the better the effluent quality is, including lower BTEX concentrations.However, the increased demand for aeration is more significant than the amount of activated carbon saved; thus, an actual SRT target is practically determined by how stringent the water quality requirements are.However, ultra-short SRT processes (<4 days) aiming for carbon removal from wastewater demonstrate high efficiency in organic matter separation and in enhancing energy recovery rates in wastewater treatment plants [67].

The Effect of Aeration Intensity
In another demonstrated case study of the developed numerical methods, the scenario analysis was conducted with varied aeration intensities of the facility's post-aeration tanks for the different steady-state runs.The decision to modify air flow to this re-aeration process-rather than changing DO setpoints in the activated sludge trains-came down to the fact that oxygen transfer will be more effective in secondary clarifier effluent water than in mixed liquor due to the absence of biomass oxygen uptake and higher α factors.The volume to be aerated is also much smaller this way, with sufficient HRT for the reaction rate of stripping.Figure 6 presents these scenarios' results on BOD-specific aeration energy and activated carbon demand as a function of superficial gas velocity-the air flow expressed per re-aeration tank surface area.By raising the air supply to the post-aeration tanks, it was revealed that despite the less significant energy requirement increase compared to the trend seen through raising the SRT, it only provides negligible cost-saving opportunities in activated carbon usage.
The overall percent removal of BTEX chemicals is already quite high (99%), even with the lowest superficial gas velocity, and given that the majority of the organic load to be polished by GAC is not volatile and cannot be removed through air stripping, there is not much room for improvement by removing some of the remaining BTEX components.It must be noted that the water samples investigated in this paper are in the context of municipal wastewater treatment plants, with the inlet concentration of BTEX compounds being marginal compared to that of other pollutants that the technology focuses on removing.The proposal of the enhanced stripping of volatile pollutants from much more concentrated industrial wastewater, especially in the petroleum industry [68], to reduce the frequency of GAC bed replacement could make much more sense.
The impact of stripping enhancement on effluent BTEX concentrations is presented in Figure 7.It may be a viable method for stabilizing the effluent water quality regarding By raising the air supply to the post-aeration tanks, it was revealed that despite the less significant energy requirement increase compared to the trend seen through raising the SRT, it only provides negligible cost-saving opportunities in activated carbon usage.
The overall percent removal of BTEX chemicals is already quite high (99%), even with the lowest superficial gas velocity, and given that the majority of the organic load to be polished by GAC is not volatile and cannot be removed through air stripping, there is not much room for improvement by removing some of the remaining BTEX components.It must be noted that the water samples investigated in this paper are in the context of municipal wastewater treatment plants, with the inlet concentration of BTEX compounds being marginal compared to that of other pollutants that the technology focuses on removing.The proposal of the enhanced stripping of volatile pollutants from much more concentrated industrial wastewater, especially in the petroleum industry [68], to reduce the frequency of GAC bed replacement could make much more sense.
The impact of stripping enhancement on effluent BTEX concentrations is presented in Figure 7.It may be a viable method for stabilizing the effluent water quality regarding micropollutants through a relatively fast intervention.Regarding environmental security, emergencies in wastewater plants-such as operational difficulties with biological treatment-may cause a sudden load increase in BTEX contaminants fed to GAC columns, in which case online monitoring could alert operators to make a swift change in the air supply to the tanks prior to the adsorption process, preventing a violation of the tertiary treatment effluent water quality limits.emergencies in wastewater plants-such as operational difficulties with biological treatment-may cause a sudden load increase in BTEX contaminants fed to GAC columns, in which case online monitoring could alert operators to make a swift change in the air supply to the tanks prior to the adsorption process, preventing a violation of the tertiary treatment effluent water quality limits.Figure 8 shows the removal percentage of the influent COD-based BTEX load through individual treatment mechanisms, as a function of the varied air flux.There is effectively no change in the contribution by biodegradation, since the re-aeration zone is located after the bioreactors, following the separation of biomass from the liquid phase by secondary clarification.Enhanced air stripping promotes the removal of the remaining BTEX constituents following biological treatment, whose concentrations are much smaller than those of the other remaining organic contaminants, therefore not being able to provide further minimization of the activated carbon requirement but potentially being useful for stabilizing water quality as a safe alternative to comply with effluent regulations.The fractional contribution of the adsorption process is inversely proportional to this trend.Figures 5 and 8 emphasize that the developed plant-wide modelling tools also provide insight into variables that cannot be directly measured, which, in conclusion, can give engineers and operators a better process understanding.This is a useful feature of modelling tools that are calibrated and virtually linked to actual facilities in real time (e.g., through digital twin implementations) [69].
The authors stress that this paper focuses on the general implementation of a methodology for plant-wide process simulation involving the fate of BTEX pollutants in biological wastewater treatment, along with options of physico-chemical effluent polishing, and the experimental cases presented are for demonstrating the feasibility of the model applications in municipal wastewater treatment.Therefore, these results shall not be generalized to other sources of municipal wastewater and certainly not to industrial treatment.The verification of robustness and accuracy shall be performed as a future effort in modelling studies, not on a benchmarking tool but against measured data from existing facilities.Figure 8 shows the removal percentage of the influent COD-based BTEX load through individual treatment mechanisms, as a function of the varied air flux.There is effectively no change in the contribution by biodegradation, since the re-aeration zone is located after the bioreactors, following the separation of biomass from the liquid phase by secondary clarification.Enhanced air stripping promotes the removal of the remaining BTEX constituents following biological treatment, whose concentrations are much smaller than those of the other remaining organic contaminants, therefore not being able to provide further minimization of the activated carbon requirement but potentially being useful for stabilizing water quality as a safe alternative to comply with effluent regulations.The fractional contribution of the adsorption process is inversely proportional to this trend.Figures 5 and 8 emphasize that the developed plant-wide modelling tools also provide insight into variables that cannot be directly measured, which, in conclusion, can give engineers and operators a better process understanding.This is a useful feature of modelling tools that are calibrated and virtually linked to actual facilities in real time (e.g., through digital twin implementations) [69].
The authors stress that this paper focuses on the general implementation of a methodology for plant-wide process simulation involving the fate of BTEX pollutants in biological wastewater treatment, along with options of physico-chemical effluent polishing, and the experimental cases presented are for demonstrating the feasibility of the model applications in municipal wastewater treatment.Therefore, these results shall not be generalized to other sources of municipal wastewater and certainly not to industrial treatment.The verification of robustness and accuracy shall be performed as a future effort in modelling studies, not on a benchmarking tool but against measured data from existing facilities.

GAC Operational Strategies
A third modelling case featured dynamic runs to compare two different GAC operational strategies, as illustrated by Figure 9. Scheduling bed replacement on a fixed time interval basis may be beneficial for sustaining stable effluent quality; however, this is not a cost-effective option, as some useful carbon adsorption capacity is lost.The use of effluent monitoring (e.g., applying online TOC probes) to synchronize bed replacement with the breakthrough capacity of the bed, i.e., when a certain effluent concentration limit is reached, may enable a more than 15% operational cost reduction involving the GAC process.[70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves [70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint

GAC Operational Strategies
A third modelling case featured dynamic runs to compare two different GAC operational strategies, as illustrated by Figure 9. Scheduling bed replacement on a fixed time interval basis may be beneficial for sustaining stable effluent quality; however, this is not a cost-effective option, as some useful carbon adsorption capacity is lost.The use of effluent monitoring (e.g., applying online TOC probes) to synchronize bed replacement with the breakthrough capacity of the bed, i.e., when a certain effluent concentration limit is reached, may enable a more than 15% operational cost reduction involving the GAC process.

GAC Operational Strategies
A third modelling case featured dynamic runs to compare two different GAC operational strategies, as illustrated by Figure 9. Scheduling bed replacement on a fixed time interval basis may be beneficial for sustaining stable effluent quality; however, this is not a cost-effective option, as some useful carbon adsorption capacity is lost.The use of effluent monitoring (e.g., applying online TOC probes) to synchronize bed replacement with the breakthrough capacity of the bed, i.e., when a certain effluent concentration limit is reached, may enable a more than 15% operational cost reduction involving the GAC process.2016) [70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves [70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint  2016) [70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves [70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint based on the adsorption capacity specification, thus aiding in decision-making about scheduling GAC bed replacement and the evaluation of operational and control strategies.

Conclusions
In this paper, the proposed kinetic concepts have been shown to provide a reasonable framework for simulating the fate of BTEX compounds in water resource recovery facilities.They can be merged with existing biokinetic models to be run in mathematical simulation software environments, as potential methods of effluent quality estimation and operational optimization in the areas of cost and resources.
The model-based benchmarking tests of this study demonstrated the applicability of the presented concepts for plant-wide modelling by investigating the impacts of operational parameters in the activated sludge process on the removal efficiency of BTEX chemicals, alongside the demand for granular activated carbon usage in effluent polishing.The first studied case revealed that the reduction in bed replacement frequency is most significant within the range of 3 to 6 days of SRT.Above 6 days, the further decrease in the AC requirement is less significant, and proportionally, a higher blower power is required to meet the DO setpoint at higher MLSS concentrations.Since BTEX constituents are only present at low concentrations in municipal settings, potential cost reductions in GAC bed replacement-through optimizing the biological treatment-may not be significant due to the higher fraction of inert and non-volatile material fed to GAC treatment, which cannot be removed through biodegradation and stripping, respectively.However, according to another studied example case, enhanced aeration for increased stripping prior to the GAC process may be a useful approach to preventing the deterioration of effluent quality regarding BTEX compounds in case of operational emergencies.The demonstration of dynamic simulation to compare operational strategies showed that synchronizing bed replacement with the activated carbon breakthrough capacity-based on the effluent quality-may result in more than 15% operational cost savings regarding the GAC process.
Goals for future studies include linking the model with further processes in water reuse.Efforts shall also be focused on model developments for off-gas treatment, considering the health risks-especially of benzene-in the atmosphere.Although the kinetic model parameters were set up based on values found in the scientific literature, future studies should aim for thorough model validation based on measured data from various facilities in municipal sewage treatment.As an outlook on degradation modelling with regard to granular activated carbon, future development efforts are dedicated to enhancing the GAC adsorption methodology with biomass activity included, enabling the simulation of BAC applications and lowered loading to an AC bed.The focus of potential further research could also be placed on applications in industrial wastewater treatment with a high BTEX content.

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

Appendix A. Gujer Matrix Development
The following tables show the MiniSumo biokinetic model structure, alongside BTEXrelated model adjustments and the affected cells of the Gujer matrix based on reference [29].All process rates are interpreted in g COD m −3 d −1 units, with the exception of the elimination of surfactants, which is interpreted in d −1 .

Environments 2024 , 34 Figure 1 .
Figure 1.Modified BSM1 test configuration for modelling examples of BTEX removal.Arrows reprentent the direction of the fluid flow.

Figure 1 .
Figure 1.Modified BSM1 test configuration for modelling examples of BTEX removal.Arrows represent the direction of the fluid flow.

Figure 3 .
Figure 3. Specific AC demand and blower energy demand for SRT-based scenarios.

Figure 3 .
Figure 3. Specific AC demand and blower energy demand for SRT-based scenarios.

Figure 5 .
Figure 5. Contribution of BTEX removal processes for SRT-based scenarios.

Figure 5 .
Figure 5. Contribution of BTEX removal processes for SRT-based scenarios.

Figure 5 .
Figure 5. Contribution of BTEX removal processes for SRT-based scenarios.

Figure 6 .
Figure 6.Specific AC demand and blower energy demand for air flux-based scenarios.

Figure 6 .
Figure 6.Specific AC demand and blower energy demand for air flux-based scenarios.

Figure 7 .
Figure 7. Effluent quality concerning BTEX compounds for air flux-based scenarios.

Figure 7 .
Figure 7. Effluent quality concerning BTEX compounds for air flux-based scenarios.

Figure 8 .
Figure 8. Contribution of BTEX removal processes for air flux-based scenarios.

Figure 9 .
Figure 9. Illustration of GAC tower operational strategies.Zietzschmann et al. (2016)[70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves[70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint

Figure 8 .
Figure 8. Contribution of BTEX removal processes for air flux-based scenarios.

Figure 9 .
Figure 9. Illustration of GAC tower operational strategies.Zietzschmann et al. (2016)[70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves[70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint

Figure 9 .
Figure 9. Illustration of GAC tower operational strategies.Zietzschmann et al. (2016)[70] conducted small-scale tests with GAC to assess the adsorption of organic micropollutants (OMP) from sewage effluent.They observed that as long as influent concentrations remained below specific thresholds, the relative breakthrough behavior remained unaffected.Consequently, the capacity for OMP removal was directly correlated with the influent OMP concentration.However, the authors noted that using the specific throughput of dissolved organic carbon alone was insufficient to generate superimposed breakthrough curves[70].In practice, when operating a full-scale facility, complete breakthrough curves shall never be able to be registered, as stable effluent quality must be maintained.We emphasize that the purpose of this research is not to generate accurate breakthrough curves but to provide a robust simulation of the breakpoint based on the adsorption capacity specification, thus aiding in decision-making about scheduling GAC bed replacement and the evaluation of operational and control strategies.

Author
Contributions: Conceptualization, D.B. and T.K.; methodology, D.B., T.W. and F.H.; software, D.B., T.W. and F.H.; data curation, D.B., T.W. and F.H.; writing-original draft preparation, D.B. and T.K.; writing-review and editing, D.B. and T.K.; visualization, D.B.; supervision, T.K.All authors have read and agreed to the published version of the manuscript.Funding: Project no.TKP2021-NVA-18 has been implemented with the support provided by the Ministry of Culture and Innovation of Hungary from the National Research, Development and Innovation Fund, financed under the TKP2021-NVA funding scheme.Data Availability Statement: The data are contained within the article.

Table 1 .
Main input parameters of the BSM1 test configuration.