Dynamic Programming and Taguchi Method Optimization of Water-Treatment-Plant Design

Table of contents

1. INTRODUCTION

he water Industry is one of the principal industries in every country. The total expenditure for water supply, including operation and maintenance, accounts to nearly 6 billion dollars per year in the U.S alone, of this approximately 15% is spent on actual water treatment (Wiesner, 1985). This massive investment requirement calls for more economical system design to allow for an efficient use of public funds. This objective can be achieved using "Optimization". Because in many cases pilot scale studies are conducted prior to full-scale plant construction, the use of an optimization model can greatly reduce the time and money spent on the pilot scale study.

A water treatment system is a combination of several unit processes. Design of such a system is a difficult task and the least cost design is most difficult of all. The principal unit process in a conventional watertreatment systemincludes coagulation-flocculation, sedimentation, filtration and disinfection. The performance of each treatment unit affects the efficacy of the subsequent units. Ideally, therefore, the design decisions should be made with regard to the interaction between various unit operations. So for very few works have been done on the optimization of a water treatment system. Letterman & Iyer (1977) analyzed the effects of selected process decision variables on overall system cost and performance. The authors used simplified process models (mostly empirical relationships). One comprehensive work is that of Wiesner et. al(1978)which is an important step towards the optimal design of water treatment system. The authors have provided a number of optimization models for integral analysis and economic optimization of the components of a water treatment plant. Mhiasalkar(1993) has developed a mathematical model incorporating the Performance relationships and cost functions for the component units of conventional water treatment system, and an algorithm using dynamic programming has described by them for functional and minimal cost design of the system. The last and one of the comprehensive models in this ground belongs to Dharmappa et. al (1994). The authors believe the incorporation of particle size distribution (PSD) is necessary for optimal process design and selection. In this context this works includes not only all three levels of system design but also process design and selection using PSD and the Algorithm they have used was a GP and NLP programming.

The major disadvantages of GP and NLP algorithms are that, global optima is not assured and they can not be directly used for system optimization because of the presence of discrete decision variables and requirement high computer storage. Using Dynamic Programming solely for large problems requires very high computer storage. Accordingly, this paper addresses it self to the development of a model for minimal cost and Energy design of a conventional water treatment plant using dynamic programming and Taguchi Design Of Experiments methods optimization model. The scope of this research was restricted to the economic optimization of conventional water treatment system for turbidity removal comprising four water treatment processes, which are listed in table 1.

The serial and interactive nature of the various unit processes favors application of Dynamic Programming for minimal cost design. Dynamic Programming is a mathematical technique often useful for making a sequence of interrelated decisions. It provides a systematic procedure for determining the combination of decisions that maximizes overall effectiveness [Hillier 1967]. Deterministic dynamic programming can be described diagrammatically, as shown in Fig. 1 Fig. 1: The basics structure for deterministic dynamic programming.

Thus at stage n the process will be in some state S n . Making policy decision X n then moves the process to some state S n+1 at stage (n+1). From that point onward the objective function value for the optimal policy has been previously calculated to be f* n+1 (S n+1 ). Conventional water Treatment systems can be considered as a multistage process, with the stages represented by various unit processes. The states represented by the levels of water quality parameters such as turbidity and suspended solids . The decision of unit.

variables are design parameters depending on the type Traditionally, measuring nonspecific parameters like turbidity units (Tu) suspended solids (SS); etc assesses the performance of water treatment plants. Though their parameters describe the Important aspects regarding the water quality of the effluent form the treatment processes (AWWA: water 1990). Insofar as turbidity is used as a measure of the efficiency of the removal of particulate matter throughout the purification process of water, low turbidity in the finished water is an indication of effectiveness of unit processes (Guidelines 1984). Because mathematical models of water treatment describe the behavior of systems in relation to the removal of mass concentration of suspended solids. And we can in some ways correlate turbidity measurements with suspended solids data, suspended solids has been used as a state parameter to indicate the effectiveness of the system at various stages of treatment. Fig. 2 Show the sequence of stages and the stage return for each input-decision combination.

2. II.

3. TAGUCHI METHOD

Taguchi approach to parameter design provides the design engineer with a systematic and efficient method for determining near optimum design parameters for performance and cost (Kackar, 1985(Kackar, : phadke, 1989;;Taguchi, 1986). In this research we have used orthogonal arrays from design of Experiments theory to study a large number of variables with a Small number of experiments. Whereupon the time of DP model analysis reduced considerably. 2011 Functional relationships for each unit process, linking input and output states parameters with a characteristic loading parameter, have been established based on the mechanics of unit process. These relationships are largely empirical models. Table 2 shows these relationships and key assumptions for each of the process models.

August R4=f(Q,T r , G r ) R3=f(Q, T s ,G s ) R2=f(Sor, U n , SS 3 , SS 2 ) R 1 =f(Fr,

4. III.

5. COST FUNCTIONS

Since the present study seeks to incorporate multiple design criteria in the optimization, there are three design criteria to be evaluated from the cost functions:

1. Annualized Capital Cost (ACC) 2. Operation & Maintenance cost (O & M) 3.

6. Energy requirement (ER)

There is a substantial literature on cost functions for water treatment plant component (Hinomoto, 1977;Clark, 1982). The annualized capital and operation and maintenance cost functions are most taken from Clark (1982). Some changes have been made by Dharmappa (1994) so as to make them more generals. Capital cost functions as presented by Clark (1982) were based on the annualization factor 0.102 (interest 8% and amortization 20 years). There were modified to incorporate a general annualization factor, of, which can be calculated by (assuming zero salvation value:

Where a f = annualization factor; i= interest rate, Fraction; and n= amortization period, years. The relationships for energy requirement are developed from Gumerman et al. (1979). Cost functions for unit processes (extracted and compiled from Clark, 1982, with Dharmappa (1994)

7. MODEL OPTIMIZATION

There are a number of designs that satisfy water quality standards. The objective, therefore, is to minimize the system cost satisfying all constrains. The bounds used in the optimizations are given in table

DHR PR A OM a CCI A ACC f f f = = 4 B Filter surface wash 4 Opt .? [R i ] i=1 V.

8. RECURSIVE EQUATIONS

For deriving Recursive equations it is necessary the Return functions for concrete, steel, filter sand and appurtenances are also necessary. Data on discount rate, and chemical costs are required for economic evaluation. The algorithm involves four stages in DP programming, and the optimal function of the previous stage is implemented in the next stage. In each stage in addition to the previous optimal function, the return function concerning to the current stage is also used to optimize the existing stage. This procedure will be followed until the last stage. The solution of this system is computed by backward form. With this parametric solution of stages, we can derive the optimal objective function in a parametric form. This parametric form includes all design parameters of the water treatment system. This objective function could be optimized by appropriate software based on the Taguchi Design of Experiments again could be calculated.

+ = Flocculation G r G r Q A Q A o T Q T Q R r

9. VII.

10. CASE STUDY

The application of the optimization algorithm is illustrated for a 6250-m3/hr conventional water treatment plant. The optimal results of this model have been compared with EPA traditional water Treatment plant. Design criteria and cost calculation of EPA conventional water treatment plant exists in EPA Documents Gumerman et al. (1979). The optimal design results in saving of 9.5% in capital and Annualized cost. The comparison is presented in table 8.

11. SENSITIVITY ANALYSIS

One particularly important feature of this model is the ease with which a sensitivity analysis can be carried out. This technique analyzes the sensitivity of the optimal solution to changes in various decision variables without re-solving the problem for each new value. It is highly unlikely that cost data of sufficient accuracy will ever be available (Walski 1991). Therefore, variation in total Annualized cost was also studied. The results of this analysis have been presented in table 9.

12. Water treatment Plant Items

13. Global ( a )

2011

14. August

The DP algorithm for functional and minimal cost design of water treatment system was developed in this research. The major inputs required are the design data on flow; raw-water suspended-solids concentration; alum dose; water temperature; and effective size, uniformity coefficient, and depth of filter sand. Data on unit costs of excavation, plain and reinforced cement VI.

15. OPTIMIZATION ALGORITHM

The sequential structure of serial systems can be exploited to transform the N-decision, one-state problems. This is accomplished by the procedure called dynamic programming, due to Richard Bellman (1957). Now by using there Return functions we can derive the Recursive equations for all of the unit processes. The general forms of recursive equations for unit processes are given in table 7.

Table 7 : Recursive equations of different unit processes . -rapid mix, flocculation, sedimentation, and rapid send filters). Mathematical models describing the performance relationships and cost functions for component units of a conventional water-treatment system have been formulated. By combining these models, an algorithm using Dynamic programming and Taguchi Design of Experiments methods software was developed for functional and minimal design of the system. The proposed approach is tested with a case study. The conclusions are listed as follows:

16. Recursive equations of Units Units

Q , A ] R .[ opt ) SS ( f f 1 2 * 1 = Filter s c 2 * 1 2 3 * 2 G , Q , A )] SS ( f R .[ opt ) SS ( f + = Sedimentation Q , T , G )] SS ( f R .[ opt ) SS ( f

? The computer times requirements for this model is very less than other similar model doing the same task.

? The optimal design results in savings of roughly 9.5% in capital and annualized cost compared to the conventional design.

? The sensitivity analysis results signifies the objective function is more sensitive to design flow rate (Q), alum dose (A) and gravity sand filter head loss (H), respectively. Suspension characteristics such as particle size and distribution are of significance in the overall process of water treatment. There are some shortcomings in major requirements for the successful incorporation of PSD in the process design and selection. (Dharmappa, 1994). Future research dynamic programming & Taguchi methods optimization of water treatment plant design may be directed to test the feasibility of considering PSD and molecular particle size distribution in water and further research in the field of software development with more levels for control factors.

Figure 1. Fig. 2 :
2Fig.2 : Unit sequences & Return of Conventional Water Treatment system
Figure 2. Table 1 :
1
State parameter
Rapid Mix Flocculation
SS 5
Figure 3. Table 2 :
2
Figure 4. Table 3 .
3
Figure 5. Table 3 :
3
Figure 6.
Dynamic Programming and Taguchi Method Optimization of Water-Treatment-Plant Design
Annual Construction Cost(ACC)
Unit processes Operation and Maintenance Cost ( O & M)
2011 Rapid mix Gr 1 ( . a 00251 ) f ) ( 181 ) 0007 DHR . 1 ( 988 . ) ( . 0 717 ) 0 . 0 CCI ( PR ) ( 79 . 0 799 . ) 0 Tr Tr × ) × Q ( Q 08507 ( 56947 . 0 . 1 = = ACC OM Gr
August Alum dry stock 10 A × ( 30366 . 3 0197 . 0 × = = ACC OM 4 Q ? ( ) 656 A . 0 × 849 ) CCI . 0 ) ( Q 994 ( . 0 PR ) . 0 1847 ( PPI ) 0 . 0259 ( DHR ) . 0 743
28 ACC = . 642 Q ( 48 × Ts ) 0 6916 . CCI ( ) 0 992 . 1 ( ) 00383 . Gs ( a f )
Volume XI Issue V Version I 1 B Flocculation 2 B Clarifier 0 B Filter structure Filter Media(single) PPI 756 ( 357 . 0 . 0 00047 . 1 ( ) ( ) n n U DHR PR ) 993 . 0 ( 181 785 . 0 . 0 ( ) ( ) a ) 1 0 006 399 . . ) ) f U CCI ) ( Ts 701 × 969 Q ( 040581 . . 0 . 0 ( ) ( 487 . 53 0 = ) ( 527 . 1872 c c PPI A A = ACC OM OM = 183 . 0 147 . 0 549 . 0 989 . 0 671 . 0 ) ( 1 61 0081 ) . . 0 ( ) ( ) ( ) ( 377 . 913 ) ( ) ( ) ( 11 . 13317 DHR PPI PR A OM a CCI A ACC f f f = = 0 = OM ) ( ) ( ) ( 3413 . 71 996 . 0 9336 . 0 = a CCI A ACC f f Gs
( a ) Researches in Engineering 3 B Filter backwash 59 65 . 0 . ( ( 801 ) 0 7146 A f ) A 1 . 6 1 . 6 ( f . 0 ) . 0 ( ) ( 4533 127 . 108 ( 398 . 747 = = ACC OM ( 732 . 204 . 568 966 ( . 0 543 ) ACC . 0 ) PPI ( a PR ( 526 . 0 ) ( ) 982 . 0 ) ( ) ) 219 . 0 315 . f ) 0 ( DHR ) . 0 137
Journal of
Global
Design variables
Where: R= Stage Returns
© 2011 Global Journals Inc. (US)
Figure 7. Table 5 :
5
Item Value
Power(PR) 0.04$/kWh
Labor Wage Rate(DHR) 11.00$/hr
Construction Cost Index/100(CCI) 3.25
Producer Price Index/100(PPI) 2.44
Alum 140$/tone
Polymer 400$/ tone
Polymer demand 6*10 -8 moles/m 2 of particle
Polymer molecular weight 5*10 4
each unit process be prepared with economic data. The assumed chemical costs are listed in
values of parameters for determining Return functions
that are used in the optimizations, as well as the
Figure 8. Table 5 .
5
Figure 9. Table 6 :
6
2011
August
29
Volume XI Issue v v Version I V
(a )
Journal of Researches in Engineering
Global
Figure 10. Table 4 :
4
1 R = . 8636 ( 59 f A . ) 671 + . 4688 ( 94 f A . ) 546 + . 46 ( 877 f A . ) 9336 +
. 293 847 ( A f ) . 59 + 2006 . ( 13 A f ) . 65 + . 363 ( 86 A f ) . 801 + 119 . ( 195 A f ) 7146 .
2 R = 1223 . 52 ( c A ) 701 . ( U n 1 ) 00047 . + 599 . 165 ( c A ) 469 . ( U n 1 ) 006 .
Alum feed & Rapid mix 54067 . . ( 12662 = 4 . ( ) . 656 . r 00160934 ) 10007 ( . 79 . 0 + ) + ( 849 3817 . . 0 ) . ( . ) 794 . . 1 ( ) 00251 +
Figure 11. Table 8 :
8
EPA WaterTreatment Plant 1 B Suggested Model Water Treatment Plant
Total Annualized Cost($) 1046570 955757
Reduction(%) 9.5
VIII.
1

Appendix A

Appendix A.1 CONCLUSION

So far very few works have been done on the optimization of a water treatment system. The scope of this research was restricted to the economic optimization of conventional water treatment systems for

The results signifies the objective function is more sensitive to Design flow rate (Q), alum dose (A) and gravity sand filter head loss (H) respectively. If the design flow is varied in the range of ±20%, the total annualized cost will vary +12.6% and 15.2 % A practical approach to analyzing either uncertain or variable situation is to test the sensitive of the optimal design to variation in key system parameters.

Appendix B

  1. , J. Envir. Engrg. Div ASCE. 108 (5) p. .
  2. water treatment plant design, Awwa Asce , Csse . 1990. New York: American water work association, Inc.
  3. Integral Water Treatment Plant Modeling: Improvements for Particle Processes. F Desmond , Jeffrey A Lawler , Nason . Environ. Sci. Technol 2005. (17) p. .
  4. , F S Hillier , Operations Research . 1974. Holden-Day, Inc. Second Edition.
  5. Introduction to Quality Engineering, G Taguchi . 1986. Dearborn, Mi: American Supplier Institute Inc.
  6. a), "water treatment system design for turbidity removal. H B Dharmappa , O Fujiwara , J Verink , S Vingneswaran . II; optimization, j. Envir. Engrg., ASCE 1994. 120 (4) p. .
  7. Capital cost functions of a surface water treatment plant and its components. H Hinomoto . water resour. Res 1977. 13 (6) p. .
  8. Optimal water treatment plant configuration: effect of raw water characteristics, M R Wiesner . 1985. Baltimore, Md. Johns Hopkins University (Ph.D. thesis)
  9. Optimal water treatment plant design. M R Wiesner . Journal of Environmental Engineering division 1987. ASCE. 113 (3) p. .
  10. Formulations for surrogate-based optimization under uncertainty. M S Eldred , A A Giunta , S F Wojtkiewicz , Jr , T G Trucano , _2002_ . Proc., 9th AIAA/ISSMO Symp. on Multidisciplinary Analysis and Optimization, (9th AIAA/ISSMO Symp. on Multidisciplinary Analysis and OptimizationAtlanta
    ) (Paper AIAA-2002-5585)
  11. Off-Line Quality Control, Parameter Design, and the Taguchi method. Raghu Kackar . journal of Quality Technology August 9. 1985. 17 (4) p. .
  12. Dynamic programming, R Bellman . 1957. Princeton, N.J: Princeton university press.
  13. Estimating water treatment costs, R C Gumerman , R L Clup , S P Hansen . 1979. Cincinnati, Ohio. 2. Municipal Environmental Research Laboratory, U.S. Environmental Protection Agency (cost curves applicable to 1 to 200 mgd treatment plants)
  14. Process model application in potable water treatment. R D Letterman , D R Iyer . 10.1021/es050089e#es050089eAF1#es050089eAF1. Proc. En grg. Foundation Conf. On theory, Practice, and process principles, (En grg. Foundation Conf. On theory, Practice, and process principles) 1977. p. . (For physical separation, pacific Grove, calif)
  15. Cost estimating for conventional water treatment, R M Clark . 1982.
  16. Quality Engineering Using Robust Design, S M Phadke . 1989. Englewood Cliffs, N.j: Prentice Hall.
  17. Discussion of preliminary optimal design model for wastewater quality and treatment, T M Walski . 1991. 1990. New York , N.Y.: Mc Graw-Hill Book Co. , Inc. p. 177. (A handbook of community water supplies. 4th Ed., F. W. Pontius, ed.)
  18. Dynamic programming optimization of water treatment plant design. V A Mhaisalkar , K Bassin , P Paramasivam , Khanna . J. Envir. Engrg 1993. 119 (6) p. .
  19. Settling and coagulating behavior of fractal aggregates. X.-Y Li , B E Logan , _2000_ . Water Sci. Technol 42 (3) p. .
Notes
1
© 2011 Global Journals Inc. (US)
Date: 2011-06-04