A method-of-lines formulation for a model of reactive settling in tanks with varying cross-sectional area

Reactive settling denotes the process of sedimentation of small solid particles dispersed in a viscous fluid with simultaneous reactions between the components that constitute the solid and liquid phases. This process is of particular importance for the simulation and control of secondary settling tanks (SSTs) in water resource recovery facilities (WRRFs), formerly known as wastewater treatment plants. A spatially one-dimensional model of reactive settling in an SST is formulated by combining a mechanistic model of sedimentation with compression with a model of biokinetic reactions. In addition, the cross-sectional area of the tank is allowed to vary as a function of height. The final model is a system of strongly degenerate parabolic, nonlinear partial differential equations (PDEs)that include discontinuous coefficients to describe the feed, underflow and overflow mechanisms, as well as singular source terms that model the feed mechanism. A finite difference scheme for the final model is derived by first deriving a method-of-lines formulation (discrete in space, continuous in time), and then passing to a fully discrete scheme by a time discretization. The advantage of this formulation is its compatibility with common practice in development of software for WRRFs. The main mathematical result is an invariant-region property, which implies that physically relevant numerical solutions are produced. Simulations of denitrification in SSTs in wastewater treatment illustrate the model and its discretization.

[1]  Ingmar Nopens,et al.  A consistent modelling methodology for secondary settling tanks in wastewater treatment. , 2011, Water research.

[2]  Raimund Bürger,et al.  A difference scheme for a degenerating convection-diffusion-reaction system modelling continuous sedimentation , 2018 .

[3]  Raimund Bürger,et al.  A family of numerical schemes for kinematic flows with discontinuous flux , 2008 .

[4]  J Alex,et al.  A combined hydraulic and biological SBR model. , 2011, Water science and technology : a journal of the International Association on Water Pollution Research.

[5]  P A Vanrolleghem,et al.  A new one-dimensional clarifier model--verification using full-scale experimental data. , 2003, Water science and technology : a journal of the International Association on Water Pollution Research.

[6]  Ingmar Nopens,et al.  A consistent modelling methodology for secondary settling tanks: a reliable numerical method. , 2013, Water science and technology : a journal of the International Association on Water Pollution Research.

[7]  Michael K Stenstrom,et al.  Practical identifiability and uncertainty analysis of the one-dimensional hindered-compression continuous settling model. , 2016, Water research.

[8]  Ulf Jeppsson,et al.  Benchmarking biological nutrient removal in wastewater treatment plants: influence of mathematical model assumptions. , 2012, Water science and technology : a journal of the International Association on Water Pollution Research.

[9]  Ulf Jeppsson,et al.  A model of the settler coupled to the biological reactor , 1998 .

[10]  Ingmar Nopens,et al.  Extending and calibrating a mechanistic hindered and compression settling model for activated sludge using in-depth batch experiments. , 2008, Water research.

[11]  I Nopens,et al.  Impact on sludge inventory and control strategies using the benchmark simulation model no. 1 with the Bürger-Diehl settler model. , 2015, Water science and technology : a journal of the International Association on Water Pollution Research.

[12]  Stefan Diehl,et al.  A conservation Law with Point Source and Discontinuous Flux Function Modelling Continuous Sedimentation , 1996, SIAM J. Appl. Math..

[13]  Ingmar Nopens,et al.  Simulations of reactive settling of activated sludge with a reduced biokinetic model , 2016, Comput. Chem. Eng..

[14]  George Tchobanoglous,et al.  Wastewater Engineering: Treatment and Resource Recovery , 2013 .

[15]  I. Takács A dynamic model of the clarification-thickening process , 1991 .

[16]  Ingmar Nopens,et al.  On reliable and unreliable numerical methods for the simulation of secondary settling tanks in wastewater treatment , 2012, Comput. Chem. Eng..

[17]  Mogens Henze,et al.  Activated sludge models ASM1, ASM2, ASM2d and ASM3 , 2015 .

[18]  Ronald L. Droste,et al.  Theory And Practice Of Water And Wastewater Treatment , 1996 .

[19]  Raimund Bürger,et al.  Linearly Implicit IMEX Runge-Kutta Methods for a Class of Degenerate Convection-Diffusion Problems , 2015, SIAM J. Sci. Comput..

[20]  A. Friedman Partial Differential Equations of Parabolic Type , 1983 .

[21]  Gerasimos Lyberatos,et al.  Modeling and Pilot-Scale Experimental Verification for Predenitrification Process , 1992 .

[22]  R. LeVeque Numerical methods for conservation laws , 1990 .

[23]  J. Ph. Chancelier,et al.  Analysis of a Conservation PDE With Discontinuous Flux: A Model of Settler , 1994, SIAM J. Appl. Math..

[24]  Paul Serban Agachi,et al.  EVALUATION OF DIFFERENT CONTROL STRATEGIES OF THE WASTE WATER TREATMENT PLANT BASED ON A MODIFIED ACTIVATED SLUDGE MODEL NO. 3 , 2012 .

[25]  Stefan Diehl Dynamic and Steady-State Behavior of Continuous Sedimentation , 1997, SIAM J. Appl. Math..

[26]  Krist V Gernaey,et al.  Effect of nitrite, limited reactive settler and plant design configuration on the predicted performance of simultaneous C/N/P removal WWTPs. , 2013, Bioresource technology.

[27]  O. Ladyženskaja Linear and Quasilinear Equations of Parabolic Type , 1968 .

[28]  U Jeppsson,et al.  Impact of reactive settler models on simulated WWTP performance. , 2006, Water science and technology : a journal of the International Association on Water Pollution Research.

[29]  Raimund Bürger,et al.  A Model of Continuous Sedimentation of Flocculated Suspensions in Clarifier-Thickener Units , 2005, SIAM J. Appl. Math..

[30]  Raimund Bürger,et al.  A simulation model for settling tanks with varying cross-sectional area , 2017 .

[31]  Min Yang,et al.  Cost-performance analysis of nutrient removal in a full-scale oxidation ditch process based on kinetic modeling. , 2013, Journal of environmental sciences.