A new plant-wide modelling methodology for WWTPs

A comprehensive evaluation of process kinetics: A plant-wide approach for nutrient removal and biogas production

The present study investigated the deviations of operational parameters of a large-scale wastewater treatment plant (WWTP) from design basis through combining dedicated batch experiments with full-scale dynamic modeling results. The long-term process performance of a full-scale biological nutrient removal (BNR) plant equipped with anaerobic sludge digestion system was monitored to evaluate the process kinetics of both carbon and nutrient removal and anaerobic sludge digestion. In this respect, plant-specific characterization; chemical oxygen demand (COD) fractionation, batch kinetic studies and sludge settling velocity tests were performed together with plant-wide SUMO model simulation. Results showed that nitrification and anaerobic hydrolysis were found to be 30% and 70% lower than literature values, respectively. The anaerobic digestion test coupled with plant-wide model calibration showed that anaerobic hydrolysis was the bottleneck in biogas production. Correspondingly, performance of the anaerobic digestion in the full-scale plant was poor as low biogas production yields were observed. In addition, the degradation rate via anaerobic hydrolysis of primary sludge was found to be higher (∼2-2.5) compared to anaerobic hydrolysis of biological sludge. The results of this study provide insight into model-based experimental characterization as well as plant-wide modeling approach. Coupling model-based batch experiments with full-scale modeling enabled to reduce the number of kinetic parameters to be fine-tuned. Moreover, the information gathered from kinetic batch tests to the simulation platform yielded a satisfying prediction of long-term performance of the plant operation.

A plant-wide model describing GHG emissions and nutrient recovery options for water resource recovery facilities

In this study, a plant-wide model describing the fate of C, N and P compounds, upgraded to account for (on-site/off-site) greenhouse gas (GHG) emissions, was implemented within the International Water Association (IWA) Benchmarking Simulation Model No. 2 (BSM2) framework. The proposed approach includes the main biological N₂O production pathways and mechanistically describes CO₂ (biogenic/non-biogenic) emissions in the activated sludge reactors as well as the biogas production (CO₂/CH₄) from the anaerobic digester. Indirect GHG emissions for power generation, chemical usage, effluent disposal and sludge storage and reuse are also included using static factors for CO₂, CH₄ and N₂O. Global and individual mass balances were quantified to investigate the fluxes of the different components. Novel strategies, such as the combination of different cascade controllers in the biological reactors and struvite precipitation in the sludge line, were proposed in order to obtain high plant performance as well as nutrient recovery and mitigation of the GHG emissions in a plant-wide context. The implemented control strategies led to an overall more sustainable and efficient plant performance in terms of better effluent quality, reduced operational cost and lower GHG emissions. The lowest N₂O and overall GHG emissions were achieved when ammonium and soluble nitrous oxide in the aerobic reactors were controlled and struvite was recovered in the reject water stream, achieving a reduction of 27% for N₂O and 9% for total GHG, compared to the open loop configuration.

Development and calibration of a new mathematical model for the description of an ion-exchange process for ammonia removal in the presence of competing ions

Ammonia ion removal and recovery via an ion-exchange process using zeolites is a promising alternative to traditional biological treatments. The analysis of its efficiency is not straightforward as it depends on various factors, such as the cation exchange capacity of the zeolite, amount of zeolite available, initial ammonia concentration, contact time, ammonia speciation depending on pH or the presence of competing ions. Mathematical modelling and simulation tools are very useful to analyse the effect of different operational conditions on the efficiency and optimal operation of the process. This paper experimentally analyses the effect that the presence of competing ions has on the efficiency of ammonia removal. This experimental work has shown a reduction of around 21% of ammonia removal efficiency in the presence of competing ions. The main contribution of this paper is the development new mathematical model able to describe the ion-exchange process in the presence of competing ions. The mathematical model developed is able to analyse the performance of the IEX process under different empty bed contact times, influent loads, pH and concentrations of competing ions. The capability of the model to reproduce real data has been proven comparing the experimental and simulation results. Finally, an exploration by simulation has been undertaken to show the potential of the mathematical model developed.

Design of Feedback Control Strategies in a Plant-Wide Wastewater Treatment Plant for Simultaneous Evaluation of Economics, Energy Usage, and Removal of Nutrients

Simultaneous removal of nitrogen and phosphorous is a recommended practice while treating wastewater. In the present study, control strategies based on proportional-integral (PI), model predictive control (MPC), and fuzzy logic are developed and implemented on a plant-wide wastewater treatment plant. Four combinations of control frameworks are developed in order to reduce the operational cost and improve the effluent quality. As a working platform, a Benchmark simulation model (BSM2-P) is used. A default control framework with PI controllers is used to control nitrate and dissolved oxygen (DO) by manipulating the internal recycle and oxygen mass transfer coefficient (KLa). Hierarchical control topology is proposed in which a lower-level control framework with PI controllers is implemented to DO in the sixth reactor by regulating the KLa of the fifth, sixth, and seventh reactors, and fuzzy and MPC are used at the supervisory level. This supervisory level considers the ammonia in the last aerobic reactor as a feedback signal to alter the DO set-points. PI-fuzzy showed improved effluent quality by 21.1%, total phosphorus removal rate by 33.3% with an increase of operational cost, and a slight increase in the production rates of greenhouse gases. In all the control design frameworks, a trade-off is observed between operational cost and effluent quality.

Response of process performance and microbial community to ammonia stress in series batch experiments

To further clarify the key stage and microorganisms responsible for ammonia inhibition instability, three sequential batch experiments were conducted with various ammonia concentrations and different exposure modes. Acetate metabolism was most sensitive to ammonia, however, after continuous ammonia exposure, acetate metabolism was well restored by a shift in dominant microorganisms. In contrast, the metabolism of longer-chain volatile fatty acids (LCVFAs, C₃-C₅) was only inhibited under a high ammonia concentration (≥6000 mg/L), however, once inhibited, continuous exposure neither restored the abundance of functional microbes nor induced new microorganisms to perform metabolic functions. Therefore, LCVFA metabolism was the key stage responsible for process instability under ammonia stress. Moreover, the deterioration of LCVFA metabolism was caused by the inhibition of syntrophic acetogenic bacteria (SAB) induced by total ammonia nitrogen, rather than the feedback inhibition from methanogenesis. That is, SAB were the key microorganisms involved in process instability.

Plant-wide modelling in wastewater treatment: showcasing experiences using the Biological Nutrient Removal Model

Plant-wide modelling can be considered an appropriate approach to represent the current complexity in water resource recovery facilities, reproducing all known phenomena in the different process units. Nonetheless, novel processes and new treatment schemes are still being developed and need to be fully incorporated in these models. This work presents a short chronological overview of some of the most relevant plant-wide models for wastewater treatment, as well as the authors' experience in plant-wide modelling using the general model BNRM (Biological Nutrient Removal Model), illustrating the key role of general models (also known as supermodels) in the field of wastewater treatment, both for engineering and research.

A plant-wide modelling comparison between membrane bioreactors and conventional activated sludge

A comprehensive plant-wide mathematical modelling comparison between conventional activated sludge (CAS) and Membrane bioreactor (MBR) systems is presented. The main aim of this study is to highlight the key features of CAS and MBR in order to provide a guide for an effective plant operation. A scenario analysis was performed to investigate the influence on direct and indirect greenhouse gas (GHG) emissions and operating costs of (i) the composition of inflow wastewater (scenario 1), (ii) operating conditions (scenario 2) and (iii) oxygen transfer efficiency (scenario 3). Scenarios show higher indirect GHG emissions for MBR than CAS, which result is related to the higher energy consumption in MBR. The simultaneous variation of the investigated factors (scenario 4) exacerbates direct and indirect GHG emissions for both CAS and MBR. Indeed, during scenario 4 a maximum direct GHG emissions of 0.94 kgCO_2eq m^-3 and 1.56 kgCO_2eq m^-3 for CAS and MBR, respectively, was obtained.

Position paper - progress towards standards in integrated (aerobic) MBR modelling

Membrane bioreactor (MBR) models are useful tools for both design and management. The system complexity is high due to the involved number of processes which can be clustered in biological and physical ones. Literature studies are present and need to be harmonized in order to gain insights from the different studies and allow system optimization by applying a control. This position paper aims at defining the current state of the art of the main integrated MBR models reported in the literature. On the basis of a modelling review, a standardized terminology is proposed to facilitate the further development and comparison of integrated membrane fouling models for aerobic MBRs.

Advanced Supervisory Control System Implemented at Full-Scale WWTP—A Case Study of Optimization and Energy Balance Improvement

In modern and cost-effective Wastewater Treatment Plants (WWTPs), processes such as aeration, chemical feeds and sludge pumping are usually controlled by an operating system integrated with online sensors. The proper verification of these data-driven measurements and the control of different unit operations at the same time has a strong influence on better understanding and accurately optimizing the biochemical processes at WWTP—especially energy-intensive biological parts (e.g., the nitrification zone/aeration system and denitrification zone/internal recirculation). In this study, by integrating a new powerful PreviSys with data driven from the Supervisory Control and Data Acquisition (SCADA) software and advanced algorithms such as Model Predictive Control (MPC) by using the WEST computer platform, it was possible to conduct different operation strategies for optimizing and improving the energy balance at a full-scale “Klimzowiec” WWTP located in Chorzow (Southern Poland). Moreover, the novel concept of double-checking online data-driven measurements (from installed DO, NO3, NH4 sensors, etc.) by mathematical modelling and computer simulation predictions was applied in order to check the data uncertainty and develop a support operator system (SOS)—an additional tool for the widely-used in-operation and control of modern and cost-effective WWTPs. The results showed that by using sophisticated PreviSys technology, a better understanding and accurate optimization of biochemical processes, as well as more sustainable WWTP operation, can be achieved.

Evaluation of anaerobic digestion post-treatment options using an integrated model-based approach

The objective of this paper is to present the main results of an engineering-research project dealing with model-based evaluation of waste streams treatment from a biotech company. This has been extensively done in domestic treatment systems, but is equally important, and with different challenges in industrial wastewater treatment. A new set of biological (activated sludge, anaerobic digestion), physicochemical (aqueous phase, precipitation, mass transfer) process models and model interfaces are required to describe removal of organics in an upflow anaerobic sludge blanket (UASB) reactor plus either traditional nitrification/denitrification (A₁) or partial nitritation (PN)/anammox (ANX) (A₂) processes. Model-based analysis shows that option A₁ requires a decrease in digestion energy recovery (E_recovery) in order to have enough organic substrate for subsequent post NO₃ reduction treatment (95 kWh.kg N^-1). In contrast, A₂ in an aerobic granular sludge reactor allows for higher UASB conversion since N removal is carried out autotrophically. The study also reveals that the addition of an aerated pre-treatment unit prior to the PN/ANX (A₂) reactor promotes COD and H₂S oxidation, CO₂ and CH₄ stripping, a pH increase (up to 8.5) and a reduction of the risk of intra-granular precipitation as well as sulfide inhibition. Simulations indicate clear differences regarding the microbial distribution/abundance within the biofilm in A₂ when comparing the two operational modes. Final results show the effects of different loading and operational conditions; dissolved oxygen (DO), Total Suspended Solids (TSS_op), energy recovery (E_recovery); on the overall process performance; N removal, aeration energy (E_aeration), net energy production (E_recovery); using response surfaces, highlighting the need of integrated approaches to avoid sub-optimal outcomes. The study shows the benefits of virtual plant simulation and demonstrates the potential of model-based evaluation when process engineers in industry have to decide between competing options.

Model-based opti mization of phosphorus management strategies in Sur WWTP, Madrid

Phosphorus has been considered as a pollutant to be removed from the wastewater. In the last years, however, it has been considered a valuable asset that needs to be recovered due to its shortage in nature. The study of optimum phosphorus management in wastewater treatment plants is not straightforward, due to the complexity of technologies and configurations that may be applied for phosphorus removal and recovery. In this context, plant-wide mathematical modelling and simulation tools are very useful for carrying out these studies. This paper introduces a study carried out at the Sur WWTP (Madrid) to assess optimum phosphorus management strategies based on the PWM. The mathematical model made it possible to describe the phosphorus flux and its characterization throughout the plant. Finally, an exploration by simulation with WEST™ was carried out to analyse different plant configurations and different operational strategies to optimize phosphorus management strategies in the Sur WWTP.

Microalgae and cyanobacteria modeling in water resource recovery facilities: A critical review

Microalgal and cyanobacterial resource recovery systems could significantly advance nutrient recovery from wastewater by achieving effluent nitrogen (N) and phosphorus (P) levels below the current limit of technology. The successful implementation of phytoplankton, however, requires the formulation of process models that balance fidelity and simplicity to accurately simulate dynamic performance in response to environmental conditions. This work synthesizes the range of model structures that have been leveraged for algae and cyanobacteria modeling and core model features that are required to enable reliable process modeling in the context of water resource recovery facilities. Results from an extensive literature review of over 300 published phytoplankton models are presented, with particular attention to similarities with and differences from existing strategies to model chemotrophic wastewater treatment processes (e.g., via the Activated Sludge Models, ASMs). Building on published process models, the core requirements of a model structure for algal and cyanobacterial processes are presented, including detailed recommendations for the prediction of growth (under phototrophic, heterotrophic, and mixotrophic conditions), nutrient uptake, carbon uptake and storage, and respiration.

Validation of a multi-phase plant-wide model for the description of the aeration process in a WWTP

This paper introduces a new mathematical model built under the PC-PWM methodology to describe the aeration process in a full-scale WWTP. This methodology enables a systematic and rigorous incorporation of chemical and physico-chemical transformations into biochemical process models, particularly for the description of liquid-gas transfer to describe the aeration process. The mathematical model constructed is able to reproduce biological COD and nitrogen removal, liquid-gas transfer and chemical reactions. The capability of the model to describe the liquid-gas mass transfer has been tested by comparing simulated and experimental results in a full-scale WWTP. Finally, an exploration by simulation has been undertaken to show the potential of the mathematical model.

Modelling oxygen transfer using dynamic alpha factors

Due to the importance of wastewater aeration in meeting treatment requirements and due to its elevated energy intensity, it is important to describe the real nature of an aeration system to improve design and specification, performance prediction, energy consumption, and process sustainability. Because organic loadings drive aeration efficiency to its lowest value when the oxygen demand (energy) is the highest, the implications of considering their dynamic nature on energy costs are of utmost importance. A dynamic model aimed at identifying conservation opportunities is presented. The model developed describes the correlation between the COD concentration and the α factor in activated sludge. Using the proposed model, the aeration efficiency is calculated as a function of the organic loading (i.e. COD). This results in predictions of oxygen transfer values that are more realistic than the traditional method of assuming constant α values. The model was applied to two water resource recovery facilities, and was calibrated and validated with time-sensitive databases. Our improved aeration model structure increases the quality of prediction of field data through the recognition of the dynamic nature of the alpha factor (α) as a function of the applied oxygen demand. For the cases presented herein, the model prediction of airflow improved by 20-35% when dynamic α is used. The proposed model offers a quantitative tool for the prediction of energy demand and for minimizing aeration design uncertainty.

Cost minimization in a full-scale conventional wastewater treatment plant: associated costs of biological energy consumption versus sludge production

Energy consumption and sludge production minimization represent rising challenges for wastewater treatment plants (WWTPs). The goal of this study is to investigate how energy is consumed throughout the whole plant and how operating conditions affect this energy demand. A WWTP based on the activated sludge process was selected as a case study. Simulations were performed using a pre-compiled model implemented in GPS-X simulation software. Model validation was carried out by comparing experimental and modeling data of the dynamic behavior of the mixed liquor suspended solids (MLSS) concentration and nitrogen compounds concentration, energy consumption for aeration, mixing and sludge treatment and annual sludge production over a three year exercise. In this plant, the energy required for bioreactor aeration was calculated at approximately 44% of the total energy demand. A cost optimization strategy was applied by varying the MLSS concentrations (from 1 to 8 gTSS/L) while recording energy consumption, sludge production and effluent quality. An increase of MLSS led to an increase of the oxygen requirement for biomass aeration, but it also reduced total sludge production. Results permit identification of a key MLSS concentration allowing identification of the best compromise between levels of treatment required, biological energy demand and sludge production while minimizing the overall costs.

Quantitative assessment of energy and resource recovery in wastewater treatment plants based on plant-wide simulations

The growing development of technologies and processes for resource treatment and recovery is offering endless possibilities for creating new plant-wide configurations or modifying existing ones. However, the configurations' complexity, the interrelation between technologies and the influent characteristics turn decision-making into a complex or unobvious process. In this frame, the Plant-Wide Modelling (PWM) library presented in this paper allows a thorough, comprehensive and refined analysis of different plant configurations that are basic aspects in decision-making from an energy and resource recovery perspective. In order to demonstrate the potential of the library and the need to run simulation analyses, this paper carries out a comparative analysis of WWTPs, from a techno-economic point of view. The selected layouts were (1) a conventional WWTP based on a modified version of the Benchmark Simulation Model No. 2, (2) an upgraded or retrofitted WWTP, and (3) a new Wastewater Resource Recovery Facilities (WRRF) concept denominated as C/N/P decoupling WWTP. The study was based on a preliminary analysis of the organic matter and nutrient energy use and recovery options, a comprehensive mass and energy flux distribution analysis in each configuration in order to compare and identify areas for improvement, and a cost analysis of each plant for different influent COD/TN/TP ratios. Analysing the plants from a standpoint of resources and energy utilization, a low utilization of the energy content of the components could be observed in all configurations. In the conventional plant, the COD used to produce biogas was around 29%, the upgraded plant was around 36%, and 34% in the C/N/P decoupling WWTP. With regard to the self-sufficiency of plants, achieving self-sufficiency was not possible in the conventional plant, in the upgraded plant it depended on the influent C/N ratio, and in the C/N/P decoupling WWTP layout self-sufficiency was feasible for almost all influents, especially at high COD concentrations. The plant layouts proposed in this paper are just a sample of the possibilities offered by current technologies. Even so, the library presented here is generic and can be used to construct any other plant layout, provided that a model is available.

Plant-wide modelling of phosphorus transformations in wastewater treatment systems: Impacts of control and operational strategies

The objective of this paper is to report the effects that control/operational strategies may have on plant-wide phosphorus (P) transformations in wastewater treatment plants (WWTP). The development of a new set of biological (activated sludge, anaerobic digestion), physico-chemical (aqueous phase, precipitation, mass transfer) process models and model interfaces (between water and sludge line) were required to describe the required tri-phasic (gas, liquid, solid) compound transformations and the close interlinks between the P and the sulfur (S) and iron (Fe) cycles. A modified version of the Benchmark Simulation Model No. 2 (BSM2) (open loop) is used as test platform upon which three different operational alternatives (A₁, A₂, A₃) are evaluated. Rigorous sensor and actuator models are also included in order to reproduce realistic control actions. Model-based analysis shows that the combination of an ammonium ( [Formula: see text] ) and total suspended solids (X_TSS) control strategy (A₁) better adapts the system to influent dynamics, improves phosphate [Formula: see text] accumulation by phosphorus accumulating organisms (X_PAO) (41%), increases nitrification/denitrification efficiency (18%) and reduces aeration energy (E_aeration) (21%). The addition of iron ( [Formula: see text] ) for chemical P removal (A₂) promotes the formation of ferric oxides (X_HFO-H, X_HFO-L), phosphate adsorption (X_HFO-H,P, X_HFO-L,P), co-precipitation (X_HFO-H,P,old, X_HFO-L,P,old) and consequently reduces the P levels in the effluent (from 2.8 to 0.9 g P.m^-3). This also has an impact on the sludge line, with hydrogen sulfide production ( [Formula: see text] ) reduced (36%) due to iron sulfide (X_FeS) precipitation. As a consequence, there is also a slightly higher energy production (E_production) from biogas. Lastly, the inclusion of a stripping and crystallization unit (A₃) for P recovery reduces the quantity of P in the anaerobic digester supernatant returning to the water line and allows potential struvite ( [Formula: see text] ) recovery ranging from 69 to 227 kg.day^-1 depending on: (1) airflow (Q_stripping); and, (2) magnesium ( [Formula: see text] ) addition. All the proposed alternatives are evaluated from an environmental and economical point of view using appropriate performance indices. Finally, some deficiencies and opportunities of the proposed approach when performing (plant-wide) wastewater treatment modelling/engineering projects are discussed.

Diagnosis and optimization of WWTPs using the PWM library: full-scale experiences

Given the shift in perception of wastewater treatment plants as water resource recovery facilities, conventional mathematical models need to be updated. The resource recovery perspective should be applied to new processes, technologies and plant layouts. The number and level of models proposed to date give an overview of the complexity of the new plant configurations and provides a wide range of possibilities and process combinations in order to construct plant layouts. This diversity makes the development of standard, modular and flexible tools and model libraries that allow the incorporation of new processes and components in a straightforward way a necessity. In this regard, the plant-wide modelling (PWM) library is a complete model library that includes conventional and advanced technologies and that allows economic and energetic analyses to be carried out in a holistic way. This paper shows the fundamentals of this PWM library that is built upon the above-mentioned premises and the application of the PWM library in three different full-scale case studies.

Validation of a plant-wide phosphorus modelling approach with minerals precipitation in a full-scale WWTP

The focus of modelling in wastewater treatment is shifting from single unit to plant-wide scale. Plant-wide modelling approaches provide opportunities to study the dynamics and interactions of different transformations in water and sludge streams. Towards developing more general and robust simulation tools applicable to a broad range of wastewater engineering problems, this paper evaluates a plant-wide model built with sub-models from the Benchmark Simulation Model No. 2-P (BSM2-P) with an improved/expanded physico-chemical framework (PCF). The PCF includes a simple and validated equilibrium approach describing ion speciation and ion pairing with kinetic multiple minerals precipitation. Model performance is evaluated against data sets from a full-scale wastewater treatment plant, assessing capability to describe water and sludge lines across the treatment process under steady-state operation. With default rate kinetic and stoichiometric parameters, a good general agreement is observed between the full-scale datasets and the simulated results under steady-state conditions. Simulation results show differences between measured and modelled phosphorus as little as 4-15% (relative) throughout the entire plant. Dynamic influent profiles were generated using a calibrated influent generator and were used to study the effect of long-term influent dynamics on plant performance. Model-based analysis shows that minerals precipitation strongly influences composition in the anaerobic digesters, but also impacts on nutrient loading across the entire plant. A forecasted implementation of nutrient recovery by struvite crystallization (model scenario only), reduced the phosphorus content in the treatment plant influent (via centrate recycling) considerably and thus decreased phosphorus in the treated outflow by up to 43%. Overall, the evaluated plant-wide model is able to jointly describe the physico-chemical and biological processes, and is advocated for future use as a tool for design, performance evaluation and optimization of whole wastewater treatment plants.

Greenhouse gases from wastewater treatment - A review of modelling tools

Nitrous oxide, carbon dioxide and methane are greenhouse gases (GHG) emitted from wastewater treatment that contribute to its carbon footprint. As a result of the increasing awareness of GHG emissions from wastewater treatment plants (WWTPs), new modelling, design, and operational tools have been developed to address and reduce GHG emissions at the plant-wide scale and beyond. This paper reviews the state-of-the-art and the recently developed tools used to understand and manage GHG emissions from WWTPs, and discusses open problems and research gaps. The literature review reveals that knowledge on the processes related to N2O formation, especially due to autotrophic biomass, is still incomplete. The literature review shows also that a plant-wide modelling approach that includes GHG is the best option for the understanding how to reduce the carbon footprint of WWTPs. Indeed, several studies have confirmed that a wide vision of the WWPTs has to be considered in order to make them more sustainable as possible. Mechanistic dynamic models were demonstrated as the most comprehensive and reliable tools for GHG assessment. Very few plant-wide GHG modelling studies have been applied to real WWTPs due to the huge difficulties related to data availability and the model complexity. For further improvement in GHG plant-wide modelling and to favour its use at large real scale, knowledge of the mechanisms involved in GHG formation and release, and data acquisition must be enhanced.

A plant-wide aqueous phase chemistry module describing pH variations and ion speciation/pairing in wastewater treatment process models

There is a growing interest within the Wastewater Treatment Plant (WWTP) modelling community to correctly describe physico-chemical processes after many years of mainly focusing on biokinetics. Indeed, future modelling needs, such as a plant-wide phosphorus (P) description, require a major, but unavoidable, additional degree of complexity when representing cationic/anionic behaviour in Activated Sludge (AS)/Anaerobic Digestion (AD) systems. In this paper, a plant-wide aqueous phase chemistry module describing pH variations plus ion speciation/pairing is presented and interfaced with industry standard models. The module accounts for extensive consideration of non-ideality, including ion activities instead of molar concentrations and complex ion pairing. The general equilibria are formulated as a set of Differential Algebraic Equations (DAEs) instead of Ordinary Differential Equations (ODEs) in order to reduce the overall stiffness of the system, thereby enhancing simulation speed. Additionally, a multi-dimensional version of the Newton-Raphson algorithm is applied to handle the existing multiple algebraic inter-dependencies. The latter is reinforced with the Simulated Annealing method to increase the robustness of the solver making the system not so dependent of the initial conditions. Simulation results show pH predictions when describing Biological Nutrient Removal (BNR) by the activated sludge models (ASM) 1, 2d and 3 comparing the performance of a nitrogen removal (WWTP1) and a combined nitrogen and phosphorus removal (WWTP2) treatment plant configuration under different anaerobic/anoxic/aerobic conditions. The same framework is implemented in the Benchmark Simulation Model No. 2 (BSM2) version of the Anaerobic Digestion Model No. 1 (ADM1) (WWTP3) as well, predicting pH values at different cationic/anionic loads. In this way, the general applicability/flexibility of the proposed approach is demonstrated, by implementing the aqueous phase chemistry module in some of the most frequently used WWTP process simulation models. Finally, it is shown how traditional wastewater modelling studies can be complemented with a rigorous description of aqueous phase and ion chemistry (pH, speciation, complexation).

Comprehensive Numerical Modeling of Greenhouse Gas Emissions from Water Resource Recovery Facilities

A numerical model was developed to comprehensively predict greenhouse gas (GHG) emissions from water resource recovery facilities. An existing activated sludge model was extended to include a nitrifier-denitrification process and carbon dioxide (CO₂) as a state variable. The bioreactor model was coupled to a process-based digester model and an empirical model of indirect CO₂emissions. Direct emissions were approximately 90% of total GHG emissions for a plantwide simulation using the Modified Ludzack-Ettinger process. Biogenic CO₂, nitrous oxide (N₂O), and methane (CH₄) represented 10, 43, and 34% of total emissions. Simulating a dissolved oxygen controlled closed-loop system reduced both sensitivity and uncertainty of GHG emissions. Nitrous oxide emissions were much more sensitive under different design and operating conditions compared to CH₄and CO₂, indicating a significant mitigation potential. An uncertainty analysis found that the uncertainty in GHGs emissions estimates could be significant. Nitrous oxide emissions dominated in both magnitude and uncertainty.

A new general methodology for incorporating physico-chemical transformations into multi-phase wastewater treatment process models

This paper introduces a new general methodology for incorporating physico-chemical and chemical transformations into multi-phase wastewater treatment process models in a systematic and rigorous way under a Plant-Wide modelling (PWM) framework. The methodology presented in this paper requires the selection of the relevant biochemical, chemical and physico-chemical transformations taking place and the definition of the mass transport for the co-existing phases. As an example a mathematical model has been constructed to describe a system for biological COD, nitrogen and phosphorus removal, liquid-gas transfer, precipitation processes, and chemical reactions. The capability of the model has been tested by comparing simulated and experimental results for a nutrient removal system with sludge digestion. Finally, a scenario analysis has been undertaken to show the potential of the obtained mathematical model to study phosphorus recovery.

New systematic methodology for incorporating dynamic heat transfer modelling in multi-phase biochemical reactors

This paper presents a new modelling methodology for dynamically predicting the heat produced or consumed in the transformations of any biological reactor using Hess's law. Starting from a complete description of model components stoichiometry and formation enthalpies, the proposed modelling methodology has integrated successfully the simultaneous calculation of both the conventional mass balances and the enthalpy change of reaction in an expandable multi-phase matrix structure, which facilitates a detailed prediction of the main heat fluxes in the biochemical reactors. The methodology has been implemented in a plant-wide modelling methodology in order to facilitate the dynamic description of mass and heat throughout the plant. After validation with literature data, as illustrative examples of the capability of the methodology, two case studies have been described. In the first one, a predenitrification-nitrification dynamic process has been analysed, with the aim of demonstrating the easy integration of the methodology in any system. In the second case study, the simulation of a thermal model for an ATAD has shown the potential of the proposed methodology for analysing the effect of ventilation and influent characterization.

Anaerobic digestion of seven different sewage sludges: a biodegradability and modelling study

Seven mixed sewage sludges from different wastewater treatment plants, which have an anaerobic digester in operation, were evaluated in order to clarify the literature uncertainty with regard to the sewage sludge characterisation and biodegradability. Moreover, a methodology is provided to determine the Anaerobic Digestion Model No. 1 parameters, coefficients and initial state variables as well as a discussion about the accuracy of the first order solubilisation constant, which was obtained through biomethane potential test. The results of the biomethane potential tests showed ultimate methane potentials from 188 to 214 mL CH4 g(-1) CODfed, COD removals between 58 and 65% and two homogeneous groups for the first order solubilisation constant: (i) the lowest rate group from 0.23 to 0.35 day(-1) and (ii) the highest rate group from 0.27 to 0.43 day(-1). However, no statistically significant relationship between the ultimate methane potential or the disintegration constant and the sewage sludge characterisation was found. Next, a methodology based on the sludge characterisation before and after the biomethane potential test was developed to calculate the biodegradable fraction, the composite concentration and stoichiometric coefficients and the soluble COD of the sewage sludge; required parameters for the implementation of the Anaerobic Digestion Model No. 1. The comparison of the experimental and the simulation results proved the consistency of the developed methodology. Nevertheless, an underestimation of the first order solubilisation constant was detected when the experimental results were simulated with the solubilisation constant obtained from the linear regression experimental data fitting. The latter phenomenon could be related to the accumulation of intermediary compounds during the biomethane potential assay.

Model selection, identification and validation in anaerobic digestion: a review

Anaerobic digestion enables waste (water) treatment and energy production in the form of biogas. The successful implementation of this process has lead to an increasing interest worldwide. However, anaerobic digestion is a complex biological process, where hundreds of microbial populations are involved, and whose start-up and operation are delicate issues. In order to better understand the process dynamics and to optimize the operating conditions, the availability of dynamic models is of paramount importance. Such models have to be inferred from prior knowledge and experimental data collected from real plants. Modeling and parameter identification are vast subjects, offering a realm of approaches and methods, which can be difficult to fully understand by scientists and engineers dedicated to the plant operation and improvements. This review article discusses existing modeling frameworks and methodologies for parameter estimation and model validation in the field of anaerobic digestion processes. The point of view is pragmatic, intentionally focusing on simple but efficient methods.

Assessment of the influence of thermal pre-treatment time on the macromolecular composition and anaerobic biodegradability of sewage sludge

Laboratory and pilot-scale experiments were carried out in order to evaluate the influence of thermal pre-treatment time on waste-activated sludge properties and anaerobic biodegradability. Six experimental conditions were analyzed from 0 to 30 min of hydrolysis time. Solubilization of macromolecular compounds, changes in the main sludge properties and anaerobic biodegradability of the sewage sludge were evaluated. A similar carbohydrate solubilization degree was achieved, from 53% to 70% and 59% to 75% for lab- and pilot-scale experiments, respectively. In the case of proteins, the values of solubilization were lower in the pilot-scale experiment than in the laboratory, with 31-45% and 47-70%, respectively. Ammonia and volatile fatty acid did not undergo important changes; however the sludge dewaterability enhanced at increased pre-treatment times. All the pre-treatment conditions had a positive effect with regard to anaerobic biodegradability and by fitting experimental data with a simplified mathematical model, it was concluded that the maximum biogas production rate is more influenced by the pre-treatment time than the total biogas production.

New generic mathematical model for WWTP sludge digesters operating under aerobic and anaerobic conditions: Model building and experimental verification

This paper presents a new mathematical model developed to reproduce the performance of a generic sludge digester working either under aerobic or anaerobic operational conditions. The digester has been modelled as two completely mixed tanks associated with gaseous and liquid volumes. The conversion model has been developed based on a plant wide modelling methodology (PWM) and comprises biochemical transformations, physicochemical reactions and thermodynamic considerations. The model predicts the reactor temperature and the temporary evolution of an extensive vector of model components which are completely defined in terms of elemental mass fractions (C, H, O, N and P) and charge density. Thus, the comprehensive definition of the model components guarantees the continuity of elemental mass and charge in all the model transformations and between any two systems defined by the model. The aim of the generic digester model is to overcome the problems that arise when trying to connect aerobic and anaerobic digestion processes working in series or to connect water and sludge lines in a WWTP. The modelling methodology used has allowed the systematic construction of the biochemical model which acts as an initial illustrative example of an application that has been experimentally verified. The variation of the temperature is also predicted based on a thermal dynamic model. Real data from four different facilities and a straightforward calibration have been used to successfully verify the model predictions in the cases of mesophilic and thermophilic anaerobic digestion as well as autothermal thermophilic aerobic digestion (ATAD). The large amount of data from the full scale ATAD and the anaerobic digestion pilot plants, all of them working under different conditions, has allowed the validation of the model for that case study.