Top Banner
General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Downloaded from orbit.dtu.dk on: Apr 01, 2019 Model identification for hindered-compression settling velocity Plósz, Benedek G.; Climent, Javier; Griffith, Christopher; Haecky, Pia; Blackburn, Nick; Chiva, Sergio; Valverde Pérez, Borja Publication date: 2018 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Plósz, B. G., Climent, J., Griffith, C., Haecky, P., Blackburn, N., Chiva, S., & Valverde Pérez, B. (2018). Model identification for hindered-compression settling velocity. Abstract from 6th IWA/WEF Water Resource Recovery Modelling Seminar (WRRmod 2018), Quebec, Canada.
11

Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Apr 02, 2019

Download

Documents

dinhtuong
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.

Users may download and print one copy of any publication from the public portal for the purpose of private study or research.

You may not further distribute the material or use it for any profit-making activity or commercial gain

You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

Downloaded from orbit.dtu.dk on: Apr 01, 2019

Model identification for hindered-compression settling velocity

Plósz, Benedek G.; Climent, Javier; Griffith, Christopher; Haecky, Pia; Blackburn, Nick; Chiva, Sergio;Valverde Pérez, Borja

Publication date:2018

Document VersionPublisher's PDF, also known as Version of record

Link back to DTU Orbit

Citation (APA):Plósz, B. G., Climent, J., Griffith, C., Haecky, P., Blackburn, N., Chiva, S., & Valverde Pérez, B. (2018). Modelidentification for hindered-compression settling velocity. Abstract from 6th IWA/WEF Water Resource RecoveryModelling Seminar (WRRmod 2018), Quebec, Canada.

Page 2: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Plósz et al.

1

Model identification for hindered-compression settling velocity

Benedek G. Plósz1, Javier Climent

2, Christopher T. Griffin

1, Pia Haecky

3, Nick Blackburn

3, Sergio Chiva

2,

and Borja Valverde-Pérez4

1Department of Chemical Engineering, University of Bath, Claverton Down, Bath BA2 7AY, UK (Email:

[email protected]; [email protected]; [email protected]) 2Universitat Jaume I, Department of Mechanical Engineering and Construction, Av. Vicent Sos Baynat, s/n 12071

Castellón (Spain), (Email: [email protected], [email protected]) 3Bioras, Hejreskovvej 18B, 3490 Kvistgaard, Denmark ([email protected], [email protected])

4Dept. of Environmental Engineering, Technical University of Denmark, Bygningstorvet, Building 115, 2800 Kgs.

Lyngby, Denmark (Email: [email protected]).

Abstract

Two of the key questions regarding secondary settling are (a) Does a process model exist for

which all hindered and compression settling velocity parameters can be estimated using

experimental data?; (b) What is the minimum data that need be inferred, from a settling sensor

setup to identify process models?” This international research effort aimed to address these

questions by carrying out a comprehensive practical identifiability assessment of constitutive

functions for hindered and compression settling velocity using laboratory-scale measurements and

one-dimensional (1-D) simulation models. For model validation, the triangulation technique was

used, including independent laboratory- and full-scale measurements as well as 1-D and

computational fluid dynamics (CFD) simulation models.

Keywords Activated sludge settling velocity; computational fluid dynamics (CFD); model identification.

INTRODUCTION

Parameter identifiability of activated sludge settling velocity models remains a challenge. The

increasing frequency of hydraulic shock events – as a result of climate change – necessitates more

effective operation and control of secondary settling tanks (SSTs) in wastewater treatment plants

(WWTPs) in the future (Ramin et al., 2014a). Theoretically, the maximum permissible SST loading

capacity determines the maximum permissible hydraulic WWTP load. However, the SST capacity

varies with sludge settleability, and thus process operation and control necessitates effective sensor

technology and identifiable simulation models (Jeppsson et al., 2013; Plósz et al., 2009). Settling

sensors should ideally provide experimental data for estimating settling velocity parameters; yet, up

to date, no simple and robust methods exist to calibrate hindered and compression settling

parameters. Derlon et al. (2017) present a cost-effective camera-based method to monitor sludge

blanket height (SBH). Ramin et al. (2014b) propose a sensor setup with a TSS sensor installed in

the bottom of a settling column, thus inferring SBH and the TSS concentration (XTSS,bottom) time-

series. Valverde-Pérez et al. (2017) demonstrate, however, that SBH and XTSS,bottom time-series do

not provide sufficient information for reliable model identification, and proposed a novel sensor

setup, additionally monitoring TSS concentration at different heights in the side of the column

(XTSS,side). Results obtained using state-of-the-art settling velocity models (Torfs et al., 2017; Ramin

et al., 2014) still suggest limitations in terms of practical identifiability of compression settling

velocity model parameters – in line with work by Li and Stenstrom (2016). As for the uncertainty

Page 3: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Session xxx WRRmod2018

2

sources associated with settling model identification, the design of settling column setups can

significantly influence measured data and thus the parameter estimates (Vanrolleghem et al., 1996;

Ekama et al., 1997). However, more research is still needed to understand better how the impact of

column size propagates to model parameters estimated. Additionally, this study addresses the

uncertainty source represented by the use of 1-D simulation models for estimating model

parameters, which are subsequently used to calibrate CFD simulation models. Triangulation is the

strategic use of multiple inquiries to address the same question, each depending on different set of

assumptions with their strengths and weaknesses (Lawlor et al., 2016). Results agreeing across

different inquiries are more likely to be replicated reliably.

The aims set in this study are (1) identifying a process model for hindered-compression settling

velocity for which all parameters can be estimated using the experimental data with both good

settling and filamentous bulking; (2) evaluating the feasibility of the sensor setup as a means to

infer experimental data on compressive solid stress; (3) assessing uncertainty sources associated

with the model identification method and the settling column design; and (4) evaluating and

validating the new settling velocity process model using the triangulation approach.

MATERIALS AND METHODS

Sampling and sensor setup. Activated sludge samples were collected in three WWTPs in Denmark

(Fredericia and Avedøre WWTPs) and one in Sweden (Ellinge) with well-settling characteristics

(Fredericia and Ellinge with SVI3.5≤90 ml/g) and filamentous bulking (Avedøre, SVI3.5~200 ml/g).

The three activated sludge processes differed in terms of operating conditions. Secondary biological

treatment in Avedøre WWTP (320 000 PE – mostly municipal sewage) and Fredericia WWTP (350

000 PE – mostly municipal sewage) were operated at solids retention time, SRT=10-15 days, and

used polymers and chlorination for bulking control, respectively. Ellinge WWTP (330 000 PE –

mostly food industrial wastewater) was operated as a high-rate system, SRT~2 days, without any

bulking control measure taken. Settling tests were carried out using the sensor prototype by

Valverde-Pérez et al. (2017), which consists of a column equipped with TSS SOLITAX (Hach,

USA) infrared sensors installed at 0.21m height in the side wall and in the bottom of the column.

a b Time elapsed, sec

0 1000 2000 3000 4000

So

lid

s c

on

ce

ntr

atio

n,

mg/L

2000

4000

6000

8000

10000

12000

14000

16000

X_init=2950 mg/L

Regression line

X_init=3440 mg/L

Regression line

X_init=4200 mg/L

Regression line

X_init=4670 mg/L

Regression line

Figure 1. The multi-probe sensor prototype developed equipped with two SOLITAX TSS sensors

installed in the bottom and the sidewall of the settling column as well as an image analysis-based

sensor with an immersed internal visible light source; (b) TSS values measured at the bottom of the

settling column (XTSS,Bottom) versus experimental time and regression lines used to estimate XTSS,Infi

values for the four settling experiments with Fredericia WWTP sludge;

Page 4: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Plósz et al.

3

Image analysis based sensor (camera) and the immersed light source are used for measuring SBH

and also to provide parallel TSSSide measurements (Fig. 1a). In the full-scale monitoring, a

SOLITAX and a SONATAX (Hach, USA) probes were used to measure the SBH and the TSS in

the bottom of the SST in the OBVA WWTP, Vila-Real, Spain. For measuring the SBH, the

threshold TSS concentration was set to 0.3 kg m-3

.

Identifiability analysis. The identifiability analysis and model calibration were done using the Latin-

Hypercube-Sampled-priors-for-Simplex (LHSS) global method (Wágner et al., 2015). In the LHSS,

the Janus coefficient (J) is used to assess the impact of parameter value variability – for cases with

covariance >0.6 – via relative predictive accuracy obtained using the upper and lower parameter

boundaries. If J~1, then we conclude we have identified parameters. 1-D simulation models were

implemented in Matlab (Matworks). The Akaike’s information criterion (AIC) was used for model

discrimination (Torfs et al. 2017).

Regression analysis. Values of the maximum solids concentration (XTSS,Infi, kg m-3

) are estimated

using the XTSS,Bottom data series obtained for each settling experiment using the regression equation

, Eq. 1

in SigmaPlot 13 with kX and fX, denoting regression parameters (Fig. 1b).

CFD simulations. The software ANSYS-CFX® (Academic Res. Release 17.2) was used to develop

the solver according to Ramin et al. (2014b). That is the solver employs an average Eulerian 2-

phase flow model. Turbulence is modelled using the k-modelMolecular viscosity of sludge is

predicted using the Herschel-Bulkley model. Additionally, the solver included the novel hindered-

compression settling model implementation. The initialization of the 2-day transient state was

explored by three different approaches: (1) defining intuitively a SBH with a constant TSS; (2)

converging a previous steady-state case with a constant influential flow; (3) using a transient state

(very costly in terms of computing time). The second choice of initialisation was eventually used.

For simulating the column, the wall-with-no slip and smooth roughness were used with fluid

velocity on the walls equalling zero.

Model validation by triangulation (MVT). The MVT addresses the question of reliabile prediction

of hindered and compression settling using the process model developed. MVT comprises two

independent approaches, i.e. (a) practical model identification using two independent sets of

laboratory-scale measurements (Ellinge and Avedøre) using the 1-D simulation model; and (b)

forward simulations of independent sets of dynamic full-scale measurement data (SBH and TSSRAS)

using a CFD solver developed. Key sources of bias for approaches a and b are the highly

degenerated simulation model structure in 1-D and the lack of estimation of parameter values other

than settling velocity parameters through the calibration of the CFD simulation model, respectively.

No specific direction of bias of these sources can be made explicit. Results from these two

approaches are compared through the CFD simulation of column tests for well-settling and

filamentous sludge (Fig. 6).

Assessment of two uncertainty sources. One of the sources of uncertainty assessed using CFD

simulations, involved the design boundary conditions of the settling column setup. The impact of

the column sensor design on the model parameters estimated was tested via forward CFD

simulations, whereby the CFD solver was calibrated with model parameters obtained for the

Fredericia sludge at 3.44 g l-1

as initial concentration (Fig. 2) and the Avedøre sludge at 3.86 g l-1

as

initial concentration (Fig. 4). The base case scenario (F=1) was that of the real setup (Fig. 1a), and

Page 5: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Session xxx WRRmod2018

4

factors (e.g., F=0.7 means 70%) were applied to resize the height and dimeter of the column,

maintaining the original proportions. Additionally, the approach of using 1-D simulation models for

estimating parameters – that are then used to calibrate CFD simulation models – was identified and

assessed as an additional uncertainty source. The predictive efficiencies of the 1-D and the 2-D

CFD simulation models were benchmarked using measured data obtained with the Fredericia

sludge at Xini=3.44 g/l as initial concentration.

RESULTS AND DISCUSSIONS

Model identification. Through an iterative approach, involving testing the practical identifiability of

parameters in a plethora of rate equations, including 2-parameter (2P) modified power, 3P

sigmoidal and 3P exponential, a 3P exponential term was identified to describe compressive solids

stress gradient, i.e.

with , Eq. 2

where the effective solids stress () gradient is formulated with vC (m2 s

-2) and rC (-) parameters.

The maximum solids concentration (XTSS,Infi, kg m-3

) is used to normalise local biomass

concentration values XTSS,I. For hindered settling velocity (vH, m s-1

), the model includes a pseudo 2-

parameter exponential constitutive function with v0 (m s-1

) and rH (m3 kg

-1), denoting the hindered

settling velocity parameters. For hindered settling, the 3-parameter logistic function by Diehl

(Diehl, 2015; Torfs et al., 2017) was also tested, in combination with the new compression model

with parameters shown in (Fig. 2).

Figure 2. Measured and simulated data for solids collected in Fredericia WWTP, posterior

parameter probability distributions, covariance matrix; AIC assessed using the new hindered-

compression process model and the Diehl hindered settling model combined with the new

compression model.

Furthermore, in Eq. 2, S and f are the sludge and water density, respectively; g denotes the

gravity acceleration constant; z is the depth in the settling column. For simulating batch settling

tests and SST, the compressive threshold concentration (XTSS,C) is set to the initial solids

concentration and the influent TSS, respectively (Guyonvarch et al., 2015). Instead of letting

Page 6: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Plósz et al.

5

parameters independently vary, the ratio of v0/rH was identified with v0 set as constant (0.0025, m d-

1). The v0/rH ratio gives a good indication of settling properties and can be linked, notably, to the

degree of sludge bulking (Wágner et al., 2015) which makes it a good controlled parameter. The

novel process model requires only three parameters to estimate (v0/rH, vC, rC) – all practically

identifiable using the experimental data obtained using the sensor. That is posterior parameter

distributions (Fig. 2) show comparably narrow confidence intervals, and although, the covariance

matrices show values >0.6 for compression parameters in some cases, parameter variability does

not significantly influence simulation outputs, i.e. J~1. (data not shown).

Validation using independent batch settling data. Independent experimental settling data – obtained

using solids with well-settling and filamentous bulking characteristics – were used to test practical

identifiability of and to validate the simulation model structure (Fig. 3 and Fig. 4).

Figure 3. Measured and simulated data for solids collected in Ellinge WWTP, posterior parameter

probability distributions (obtained using 250 LHSS simulations), covariance matrix.

As for the Ellinge data (Fig. 3), results obtained show close agreement with the outcomes in the

Fredericia case (Fig. 2) in terms of predictive accuracy and parameter covariance. TSSside is

effectively predicted through all three experiments. Additionally, at TSS=3.76 g/l, prediction of the

SBH improves, which is not the case for the TSSbottom data series, thereby leading to ~1.5 g/l

overestimation of the measured data.

Figure 4. Measured and simulated data for solids collected in Avedøre WWTP, posterior parameter

probability distributions (obtained using 250 LHSS simulations), covariance matrix.

In contrast to the Fredericia and Ellinge datasets, validation using solids collected in Avedøre

WWTP extended the model identifiability boundaries to filamentous bulking conditions (Fig. 4).

Page 7: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Session xxx WRRmod2018

6

Again, the outcomes of the identifiability test closely agree with the Fredericia case. Compared to

Fredericia and Ellinge cases, improved prediction of the SBH and TSSbottom is obtained with bulking

sludge. Taken together, the independent results obtained with Ellinge and Avedøre solids suggest

the validity of the identifiability approach and the simulation model structure. The reliability of the

process model is further supported by the improved predictive efficiency – in terms of both SBH

and TSSBottom - under filamentous bulking conditions – an important aspect for future development

of model-based control design structures for WWTPs.

Parameter intervals for the new model. Fig. 5 summarises all parameter values with confidence

intervals obtained with the three solids. As for Fig. 5a., fixing v0 was found to allow the estimation

of v0/rH values in a narrow range for the different initial concentrations and independently from the

compression parameters according to the covariance matrices obtained. This was otherwise

impossible to achieve with any of the functions tested. Fig. 5a also supports the hypothesis of v0/rH

effectively gauging sludge settling properties (Wágner et al., 2015) with v0/rH ~ 0.005, indicating

the boundary between well-settling and filamentous bulking solids.

a

v0/r

H [kg m

-2s

-1]

XIni

[kg m-3

]

2 3 4 5

0.002

0.004

0.006

0.008

Estimated - FredericiaEstimated - AvedoreEstimated - Ellinge

b

vC [m

2 s

-2]

XIni

[kg m-3

]

2 3 4 5

0.0

0.1

0.2

0.3

0.4Estimated - FredericiaEstimated - AvedoreEstimated - EllingeRegressionRegressionRegression

cr C

[-]

XIni

[kg m-3

]

2 3 4 5

0

1

2

3

4

5 Estimated - FredericiaEstimated - AvedoreEstimated - EllingeRegressionRegressionRegression

Figure 5. Posterior mean parameter values with confidence interval denoted with error bars for the

three WWTPs.

For vC – denoting the maximum compressive solid stress gradient parameter – significant

dependence on initial solids concentrations is obtained – an observation that cannot be fully

supported for rC (Fig. 5b and 5c). Notably, vC parameter values indicate different trends under well-

and filamentous-settling conditions. Under bulking conditions, vC values show an approximately

ten-fold increase towards low solids concentrations compared well-settling sludge. For the latter

case, the trend in vC values can be characterised with a minimum range at comparably low initial

TSS concentrations – close agreement between Fredericia and Ellinge data – and with a progressive

increase towards high initial TSS concentrations. Despite the considerable difference between the

three WWTPs in terms of operating conditions - in terms of SRT and bulking control measures –

settling parameters obtained show consistent and comparable trends.

Assessing sources of uncertainty using CFD simulations. CFD simulations of the settling column

setup (Fig. 6; at design factor F=1) indicate negligible uncertainties introduced by the 1-D

parameter estimation approach, and thus suggest the reliability of the parameter estimation

approach. We note that the improved predictive efficiency of CFD simulation model compared to

the 1-D case (Fig. 6a and 6b) – in terms of SBH – suggest that the overestimation of the

compressive SBH tail by the 1-D simulation model may be a bias caused by the degenerated 1-D

simulation model structure rather than the settling velocity process model structure. The latter was

the same in both the 1-D and the 3-D CFD model. Torfs et al., (2017) assessed the effect of

overestimation of the compressive SBH tail in more depth, suggesting that the 1-D simulation

Page 8: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Plósz et al.

7

model structure - in terms of hindered settling velocity formulation – as the potential cause of such

bias.

0.0

0.2

0.4

0.6

0.8

1.0

1.2

0 20 40 60

Experimental

F=1.50

F=1.00

F=0.90

F=0.70

F=0.50

SBH (m)

0.0

1.0

2.0

3.0

4.0

5.0

6.0

7.0

8.0

9.0

0 20 40 60

XSide, (kg/m3)

Figure 6. CFD simulation results obtained using solver calibrated according to parameter values

obtained with Fredericia sludge at Xini=3.44 g/l (Fig. 2) at different design similarity factors (F)

compared to the real setup (F=1; Fig. 1a). Fig. 6b shows an excerpt of 1-D results from Fig. 2.

Furthermore, to assess the variability of parameter values as a result of settling column design, CFD

simulations, carried out within a wide range of column design boundary conditions (Fig. 6), were

used to re-estimate settling velocity model parameters (Fig. 7). Results obtained suggest that, with

negligible wall effects, only values of rC can be expected to vary significantly in the wide design

boundary range studied for both well- and filamentous sludge settling.

Column sensor size factor0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

0.01

0.1

v0/rH

vC

rC

Reference

1

Para

mete

r valu

es

Column sensor size factor0.4 0.6 0.8 1.0 1.2 1.4 1.6

0.01

0.1

1v0/rH

vC

rC

Reference

ba

Figure 7. Settling model parameters estimated with different column design using CFD simulation

output data obtained using calibration parameter sets for well-settling sludge from Fredericia

WWTP (a) and sludge with filamentous bulking collected in Avedøre WWTP (b)

Full-scale measurements and CFD simulations. As part of the MVT approach, forward CFD

simulation results (Fig. 8) closely agree with full-scale SST measurement data collected during

more than 40 hours – in terms of SBH and XTSS,RAS.

Page 9: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Session xxx WRRmod2018

8

a

Time elapsed, h

0 10 20 30 40 50

Slu

dg

e b

lan

ket

heig

ht,

m

0.6

0.7

0.8

0.9

1.0

1.1

MeasuredCFD simulation

b

Time elapsed, h

0 10 20 30 40 50

TS

S in

slu

dg

e r

ec

yc

le, g

/m3

2000

2500

3000

3500

4000

4500

5000

MeasuredCFD simulation

Figure 8. Measured and simulated (a) SBH and (b) TSSRAS data for the full-scale SST monitored.

Taken together, both approaches involved in the MVT support the hypothesis that the novel

constitutive function for hindered and compression settling velocity can reliably predict the real

physical phenomena, thereby validating the process model developed.

Quantifying compressive solid stress using sensor. This study also addressed the question whether

the multi-probe sensor setup could be used to quantify the -gradient – a variable approximated

using the sensor data according to

, Eq. 3

where the density difference between water and sludge ( ) was assumed constant. Eq. 3 was

identified based on force balance analysis – assuming only the gravitational, buoyancy and solids

pressure acting on particles – and by assuming quasi steady-state conditions (Xu et al., 2017).

Simulation results obtained (Fig. 9) reasonably agree with the sensor-gradient values for sludge

with well-settling and filamentous bulking characteristics.

a X(i)/Xinfin

0.0 0.2 0.4 0.6 0.8

d

0.5

1.0

1.5

2.0

2.5

2.95

2.95 simulation

3.44

3.44simulation

4.22

4.22simulation

4.67

4.67simulation

b X(i)/Xinfin

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4

0.5

1.0

1.5

2.0

2.5

2.43

2.43 simulation

2.83

2.83simulation

3.31

3.31simulation

3.38

3.38simulation

Figure 3. Measured/simulated compressive solid stress using (a) well-settling sludge from

Fredericia WWTP and (b) sludge with filamentous bulking collected in Avedøre WWTP.

This result indicates the feasibility of the sensor approach to quantify solid stress. More research is

needed to assess the error introduced by assuming quasi steady-state in approximating compressive

Page 10: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Plósz et al.

9

solids stress using the sensor data. Additionally, the benefits of using -gradient sensor data for

settling model calibration will be evaluated.

CONCLUSIONS

The concluding remarks drawn in the study include

A pseudo 2P and a 3P exponential term were identified to describe hindered settling velocity

and the compressive solids stress gradient, respectively;

The ratio of v0/rH was estimated with v0 set as constant;

Three parameters are required to estimate using LHSS (v0/rH, vC, rC) – all practically

identifiable using the data obtained using the innovative multi-probe sensor setup;

Only vC shows significant dependence on initial solids concentration;

The process model developed was validated using the triangulation approach, including

independent laboratory- and full-scale measurement data and using 1-D and CFD simulation

models;

Negligible uncertainties – assessed by means of CFD simulations – introduced by the 1-D

parameter estimation approach were obtained, thus suggesting the reliability of the practical

identifiability assessment approach.

The multi-probe settling sensor setup developed can be used to quantify the -gradient, and

future research should assess the benefits of using -gradient sensor data for settling model

calibration.

REFERENCES

Derlon, N., Thürlimann, C., Dürrenmatt, D., Villez, K. (2017). Batch settling curve registration via image data

modelling. Water Res., 114, 327-337.

Diehl, S., (2015). Numerical identification of constitutive functions in scalar nonlinear convection - diffusion equations

with application to batch sedimentation. Appl. Numer. Math. 95, 154-172.

Ekama, G.A., Barnard, J.L., Günthert, F.W., Krebs, P., McConcordale, J.A., Parker, D.S., Wahlberg, E.J., (1997).

Secondary settling tank: theory, modelling, design and operation. Scientific and Technical Report No. 6. IAWQ,

London.

Guyonvarch, E., Ramin, E., Kulahci, M., Plósz, B.G. (2015). iCFD: Interpreted computational fluid dynamics -

Degeneration of CFD to one-dimensional advection-dispersion models using statistical experimental design -

The secondary clarifier. Water Res., 83, 396-411.

Jeppsson, U., Alex, J., Batstone, D.J., Benedetti, L., Comas, J., Copp, J.B., Corominas, L., Flores-Alsina, X., Gernaey,

K.V., Nopens, I., Pons, M.N., Rodriguez-Roda, I., Rosen, C., Steyer, J.P., Vanrolleghem, P.A., Volcke, E. I. P.,

Vrecko, D., 2013. Benchmark simulation models, quo vadis? Water Sci. Technol. 68, 1, 1-15.

Lawlor, D.A., Tilling, K., Smith, G.D., (2016). Triangulation in aetiological epidemiology. Int. J. Epidemiology, 45, 6,

1866-1886.

Li, B., Stenstrom, M.K., (2016). Practical identifiability and uncertainty analysis of the one-dimensional hindered-

compression continuous settling model. Water Res., 90, 235-246.

Plósz, B.G., Liltved, H., Ratnaweera, H., 2009. Climate change impacts on activated sludge wastewater treatment: A

case study from Norway. Water Sci. Technol. 60, 2, 533–541.

Ramin, E., Sin, G., Mikkelsen, P.S., and Plósz, B.G., (2014a). Significance of settling model structures and parameter

subsets in modelling WWTPs under wet-weather flow and filamentous bulking conditions. Water Res. 63, 209-

221.

Ramin, E., Wágner, D.S., Yde, L., Binning, P.J., Rasmussen, M.R., Mikkelsen, P.S., Plósz, B.G. (2014b). A new

settling velocity model to describe secondary sedimentation. Water Res., 66, 447-458.

Torfs, E., Balemans, S., Locatelli, F., Diehl, S., Bürger, R., Laurent, J., Francois, P., Nopens, I., (2017). On constitutive

Page 11: Model identification for hindered-compression settling ...orbit.dtu.dk/ws/files/145256914/WRRmod2018_Settling_model... · implemented in Matlab (Matworks). The Akaike’s information

Session xxx WRRmod2018

10

functions for hindered settling velocity in 1-D settler models: selection of appropriate model structure. Water

Res., 110, 38-47.

Valverde Pérez, B., Penkarski-Rodon, E., Zhang, X., Wágner, D.S., Haecky, P., Blackburn, N.D., Plósz, B.G. (2017).

Development and validation of a novel monitoring system for batch flocculant solids settling process. IWA-

ICA2017–12th

IWA Specialist Conf. on Instrumentation, Automation and Control. Quebec City, Canada.

Vanrolleghem, P.A., Van der Schueren, D., Krikilion, G., Grijspeerdt, K., Willems, P., Verstraete, W., 1996. On-line

quantification of settling properties with in-sensor-experiments in an automated settlometer. Water Sci. Technol.

33,1, 37–51.

Wágner, D.S., Valverde-Pérez, B., Sæbø, M., Bregua de la Sotilla, M., Van Wagenen, J., Smets, B.F., Plósz, B.G.,

(2016). Towards a consensus-based biokinetic model for green microalgae–the ASM-A. Water Res., 103, 485-

499.

Wágner, D.S., Ramin, E., Szabo, P., Dechesne, A., Plósz, B.G. (2015). Microthrix parvicella abundance associates with

activated sludge settling velocity and rheology–Quantifying and modelling filamentous bulking. Water Res., 78,

121-132.

Xu, G., Yin, F., Xu, Y., Yu, H.Q., (2017). A force-based mechanistic model for describing activated sludge settling

process. Water Res., 127, 118-126.