\`x^2+y_1+z_12^34\`
Advanced Search
Article Contents
Article Contents

An optimal control problem applied to a wastewater treatment plant

  • * Corresponding author: M. Teresa T. Monteiro

    * Corresponding author: M. Teresa T. Monteiro 
Abstract Full Text(HTML) Figure(6) / Table(4) Related Papers Cited by
  • This paper aims to present a mathematical model that describes the operation of an activated sludge system during one day. Such system is used in the majority of wastewater treatment plants and depends strongly on the dissolved oxygen, since it is a biological treatment. To guarantee the appropriate amount of dissolved oxygen, expensive aeration strategies are demanded, leading to high costs in terms of energy consumption. It was considered a typical domestic effluent as the wastewater to test the mathematical model and it was used the ASM1 to describe the activated sludge behaviour. An optimal control problem was formulated whose cost functional considers the trade-off between the minimization of the control variable herein considered (the dissolved oxygen) and the quality index that is the amount of pollution. The optimal control problem is treated as a nonlinear optimization problem after discretization by direct methods. The problem was then coded in the AMPL programming language in order to carry out numerical simulations using the NLP solver IPOPT from NEOS Server.

    Mathematics Subject Classification: Primary: 49J15, 65K05; Secondary: 97M40, 97M60.

    Citation:

    \begin{equation} \\ \end{equation}
  • 加载中
  • Figure 1.  Schematic view of the activated sludge system

    Figure 2.  Hourly results during one day with $ (\omega_1,\omega_2) = (1,0) $

    Figure 3.  Hourly results during one day with $ (\omega_1,\omega_2) = (0.5,0.5) $

    Figure 4.  Hourly results during one day with $ (\omega_1,\omega_2) = (0,1) $

    Figure 5.  Hourly results during one day with $ S_{inf} = 125 $

    Figure 6.  Hourly results for the law limits during one day with $ S_{inf} = 125 $

    Table 1.  Variables of the problem and initial values

    Variable Init.val. Variable Init.val. Variable Init.val.
    $ Q $ 4000 $ Q_w $ 100 $ Q_r $ 2000
    $ Q_{ef} $ 1900 $ X_I $ 727.3 $ X_{I_r} $ 950.571
    $ X_{I_{ef}} $ $ 10^{-5} $ $ S_{S_{in}} $ 50 $ S_S $ 10
    $ S_{O_{in}} $ 1 $ S_{NO_{in}} $ $ 10^{-6} $ $ S_{NO} $ $ 10^{-6} $
    $ X_{BH_{in}} $ 0.2 $ X_{BH} $ 350 $ X_{BH_{r}} $ 711.2
    $ X_{BH_{ef}} $ $ 10^{-5} $ $ X_{S_{in}} $ 3000 $ X_S $ 350
    $ X_{S_r} $ 806.714 $ X_{S_{ef}} $ $ 10^{-5} $ $ X_{BA_{in}} $ $ 10^{-5} $
    $ X_{BA} $ $ 10^{-6} $ $ X_{BA_r} $ 1.9$ \times 10^{-6} $ $ X_{BA_{ef}} $ $ 10^{-5} $
    $ S_{NH_{in}} $ 10 $ S_{NH} $ 7.5 $ X_{P_{in}} $ 1500
    $ X_{P} $ 90 $ X_{P_{r}} $ 174.52 $ X_{P_{ef}} $ $ 10^{-5} $
    $ S_{ND_{in}} $ 0.5 $ S_{ND} $ 0.5 $ X_{ND_{in}} $ 200
    $ X_{ND} $ 20 $ X_{ND_{r}} $ 20 $ X_{ND_{ef}} $ 0.5
    $ G_S $ 10000 $ SSI $ 1500 $ SSI_{ef} $ 0.1
    $ SSI_r $ 3500 $ HRT $ 3.5 $ r $ 1
    $ S_O $ 2
     | Show Table
    DownLoad: CSV

    Table 2.  Parameters

    Kinetic Operational Stoichiometric
    $ \mu_H $ 6 $ T $ 20 $ Y_A $ 0.24
    $ \mu_A $ 0.8 $ P_{O_2} $ 0.21 $ Y_H $ 0.666
    $ k_h $ 3 $ SRT $ 20 $ f_P $ 0.08
    $ k_a $ 0.08 $ \theta $ 1.024 $ i_{X_B} $ 0.086
    $ b_h $ 0.62 $ \alpha $ 0.8 $ i_{X_P} $ 0.06
    $ b_a $ 0.04 $ \eta $ 0.07
    $ \eta_g $ 0.8 $ \beta $ 0.95
    $ \eta_h $ 0.4
    $ K_S $ 20
    $ K_X $ 0.03
    $ K_{OH} $ 0.2
    $ K_{NO} $ 0.5
    $ K_{NH} $ 1
    $ K_{OA} $ 0.4
     | Show Table
    DownLoad: CSV

    Table 3.  Characteristics of the wastewater entering the system

    $ Q_{inf} $ 530
    $ S_{I_{inf}} $ 5.45, 12.5 and 25
    $ S_S $ 44.55,112.5 and 225
    $ S_{inf} $ 50,125 and 250
    $ X_{BH_{inf}} $ 0
    $ X_{BA_{inf}} $ 0
    $ X_{P_{inf}} $ 0
    $ S_{O_{inf}} $ 0
    $ S_{NO_{inf}} $ 0
    $ S_{alk_{inf}} $ 7
    $ X_{I_{inf}} $ 90
    $ X_{S_{inf}} $ 168.75
    $ S_{NH_{inf}} $ 11.7
    $ X_{S_{inf}} $ 168.75
    $ S_{NH_{inf}} $ 11.7
    $ S_{ND_{inf}} $ 0.63
    $ X_{ND_{inf}} $ 1.251
    $ X_{II} $ 18.3
     | Show Table
    DownLoad: CSV

    Table 4.  Results for three different amounts of pollution in the wastewater

    Scenario A Scenario B Scenario C
    $ (w_1,w_2)=(1,0) $ ($ w_1,w_2)=(0.5,0.5) $ $ (w_1,w_2)=(0,1) $
    $ S_{inf} $ 50 125 250 50 125 250 50 125 250
    $ S_O $ 60.6 66.4 80.2 71.5 73.8 77.9 72.0 73.4 80.2
    $ QI $ 2390.6 2849.9 813.1 692.3 795.4 857.2 750.2 750.4 813.1
     | Show Table
    DownLoad: CSV
  • [1] L. Amand and B. Carlsson, Optimal aeration control in a nitrifying activated sludge process, Water Research, 46 (2012), 2101-2110. 
    [2] M. S. Bazaraa and C. M. Shetty, Nonlinear Programming: Theory and Algorithms, John Wiley & Sons, New York-Chichester-Brisbane, 1979.
    [3] C. Binnie, M. Kimber and G. Smethurst, Basic Water Treatment, Royal Society of Chemistry, Cambridge, UK, 2002.
    [4] S. CristeaC. de PradaD. Sarabia and G. Gutierrez, Aeration control of a wastewater treatment plant using hybrid NMPC, Computers & Chemical Engineering, 35 (2011), 638-650. 
    [5] J. CzyzykM. P. Mesnier and J. J. Moré, The EOS server, IEEE Computational Science and Engineering, 5 (1998), 68-75.  doi: 10.1109/99.714603.
    [6] I. Espírito Santo, Desenho Óptimo de Estações de Águas Residuais Através da ModelaÇão dee FunÇões Custo, Ph.D thesis, Universidade do Minho, Braga, 2007.
    [7] M. FikarB. Chachuat and M. A. Latifi, Optimal operation of alternating activated sludge processes, Control Engineering Practice, 13 (2005), 853-861.  doi: 10.1016/j.conengprac.2004.10.003.
    [8] R. FourerD. M. Gay and B. Kernighan, AMPL: A mathematical programming language, Algorithms and Model Formulations in Mathematical Programming, 51 (1989), 150-151.  doi: 10.1007/978-3-642-83724-1_12.
    [9] G. FuD. Butler and S. T. Khu, Multiple objective optimal control of integrated urban wastewater systems, Environmental Modelling & Software, 23 (2008), 225-234. 
    [10] M. Henze, C. P. L. Grady Jr, W. Gujer, G. V. R. Marais and T. Matsuo, Activated Sludge Model No.1, Technical Report 1, IAWPRC Task Group on Mathematical Modelling for design and operation of biological wastewater treatment, 1986.
    [11] H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000), 599-653. 
    [12] B. HolendaE. DomokosÁ. Rédey and J. Fazakas, Aeration optimization of a wastewater treatment plant using genetic algorithm, Optimal Control Applications and Methods, 28 (2007), 191-208.  doi: 10.1002/oca.796.
    [13] R. HreizabN. RocheB. Benyahiac and M. A. Latifia, Multi-objective optimal control of small-size wastewater treatment plants, Chemical Engineering Research and Design, 102 (2015), 345-353. 
    [14] L. S. Pontryagin, The Mathematical Theory of Optimal Processes, Mathematics & Statistics for Engineers, 2017. doi: 10.1201/9780203749319.
    [15] J. F. QiaoY. C. BoW. Chai and H. G. Han, Adaptive optimal control for a wastewater treatment plant based on a data-driven method, Water Sci Techno, 67 (2013), 2314-2320.  doi: 10.2166/wst.2013.087.
    [16] H. S. Rodrigues, Optimal Control and Numerical Optimization Applied to Epidemiological Models, PhD thesis, University of Aveiro, 2012.
    [17] H. S. RodriguesM. T. T. Monteiro and D. F. M. Torres, Optimization of dengue epidemics: A test case with different discretization schemes, AIP Conference Proceedings, 1168 (2009), 1385-1388.  doi: 10.1063/1.3241345.
    [18] H. S. Rodrigues, M. T. T. Monteiro and D. F. M. Torres, Optimal Control and numerical software: An overview, Systems Theory: Perspectives, Applications and Developments, (2014), 93–110.
    [19] C. SweetappleG. Fu and D. Butler, Multi-objective optimisation of wastewater treatment plant control to reduce greenhouse gas emissions, Water Research, 55 (2014), 52-62.  doi: 10.1016/j.watres.2014.02.018.
    [20] A. Wätcher and L. T. Biegler, On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming, Math. Program., 106 (2006), 25-57.  doi: 10.1007/s10107-004-0559-y.
  • 加载中

Figures(6)

Tables(4)

SHARE

Article Metrics

HTML views(199) PDF downloads(249) Cited by(0)

Access History

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return