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.
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,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.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)
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.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 rThe 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.
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.
2011
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.
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:
? 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.

| State parameter | |
| Rapid Mix | Flocculation |
| SS 5 |
| 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) | |||||||||||||||||||||||||||||||||||
| 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 | ||
| 2011 |
| August |
| 29 |
| Volume XI Issue v v Version I V |
| (a ) |
| Journal of Researches in Engineering |
| Global |
| 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 | + | ||||||||||||||||||||||||||||||||||||
| EPA WaterTreatment Plant | 1 B Suggested Model Water Treatment Plant | |
| Total Annualized Cost($) | 1046570 | 955757 |
| Reduction(%) | 9.5 | |
| VIII. |
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.
Integral Water Treatment Plant Modeling: Improvements for Particle Processes. Environ. Sci. Technol 2005. (17) p. .
a), "water treatment system design for turbidity removal. II; optimization, j. Envir. Engrg., ASCE 1994. 120 (4) p. .
Capital cost functions of a surface water treatment plant and its components. water resour. Res 1977. 13 (6) p. .
Optimal water treatment plant design. Journal of Environmental Engineering division 1987. ASCE. 113 (3) p. .
Formulations for surrogate-based optimization under uncertainty. Proc., 9th AIAA/ISSMO Symp. on Multidisciplinary Analysis and Optimization, (9th AIAA/ISSMO Symp. on Multidisciplinary Analysis and OptimizationAtlanta
Off-Line Quality Control, Parameter Design, and the Taguchi method. journal of Quality Technology August 9. 1985. 17 (4) p. .
Process model application in potable water treatment. 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)
Dynamic programming optimization of water treatment plant design. J. Envir. Engrg 1993. 119 (6) p. .
Settling and coagulating behavior of fractal aggregates. Water Sci. Technol 42 (3) p. .