## Abstract

Groundwater contamination due to contaminants like trichloroethylene (TCE), tetrachloroethylene, dichloroethylene, phenol, etc., is an alarming concern for most of the manufacturing areas. It is important to identify the type of pollutant, concentration, location, and direction of the contaminant plume for groundwater remediation. Bioremediation has been identified as one of the important remediation techniques for these types of contaminants. Bioremediation modeling comprises solutions to biodegradation equations and fixing the time of remediation and locating the oxygen injection wells. In this study, a simulation-optimization (S/O) model based on the coupled meshfree point collocation method (MFree-PCM) and particle swarm optimization (PSO) is proposed for *in-situ* bioremediation design. The *in-situ* bioremediation process of groundwater contamination is explored using the developed PCM-BIO-PSO multi-objective model with different strategies of minimization of cost, number of wells and time of remediation. The proposed model can be effectively used for the *in-situ* bioremediation design of contaminated sites.

## INTRODUCTION

The requirement of groundwater increases with the increase of population, industrialization, agriculture, and navigation. The groundwater quality has also become a major concern in terms of color, odor, taste, and hardness. The conservation of the quality and quantity of the groundwater requires supreme priority. Any groundwater contamination study involves identification of the source of contaminant, the movement of the pollutants in the groundwater environment once they are introduced, proper management of resource utilization, and preventive measures to ensure suitable development and remediation of polluted sites (Eldho 2001).

For the remediation of groundwater contaminated with pollutants such as trichloroethylene (TCE), dichloroethylene or phenol, conventional approaches such as pump and treat (PAT), pump and use are not effective. In such cases, *in-situ* bioremediation is one of the effective methods for remediation (Bedient *et al.* 1999).

In a typical *in-situ* bioremediation system, an aerobic environment is created by adding oxygenated fluid for the growth of microorganisms. Biological degradation is facilitated by oxygenated fluid. The process continues until the concentration of contaminants is reduced to the approved levels. The type of oxygenated fluid to be used is a further challenge for this method. The type of oxygenated fluid is preferably air, but pure oxygen, steam, water vapor, etc. are some other types used (Bedient *et al.* 1999).

The purpose of modeling bioremediation of groundwater is for the prediction of contaminant concentrations at receptor positions and the estimation of source control actions on remediation. Further, it also gives the point allocations of observation wells while the remedial study is carried out and evaluates the remedial functioning for different designs. Modeling of bioremediation of groundwater is an effective tool to predict the efficacy and clean-up time to stated levels (Bedient *et al.* 1999).

Numerical models are necessary for developing an efficient operative design and upkeep of an *in-situ* bioremediation arrangement. Simulation-optimization (S/O) models can give an efficient injection approach, suitable positions for remediation wells and time of remediation. The efficacy of S/O models for bioremediation is demonstrated by various studies (e.g., Minsker & Shoemaker 1998; Yoon & Shoemaker 2001; Prasad & Mathur 2008; Mategaonkar & Eldho 2012c).

Meshfree (MFree) methods have gained acceptance in many engineering applications owing to their mesh-free character. With the removal of grid, pre-processing time of simulation is reduced. With the ease in computing the distances in any number of spatial dimensions, working in higher dimensions does not increase the intricacy and computational cost of the method (Liu 2003, 2006; Liu & Gu 2005).

MFree methods have been applied successfully to groundwater flow and transport problems by many researchers (e.g., Li *et al.* 2003; Praveen Kumar & Dodagoudar 2008, 2010). Mategaonkar & Eldho (2012a; 2012b) developed a coupled flow and transport model for groundwater.

Simulation models along with optimization models for the effective remediation of groundwater pollution give the complete solution for the complex groundwater contamination problem. Shieh & Peralta (2005) presented a parallel recombinative simulated annealing (PRSA) model along with BIOPLUME II to optimize *in-situ* bioremediation system design with a two-stage management approach. The first-stage design goal is to minimize total system cost (pumping/treatment, well installation, and facility capital costs) and the second-stage goal is to minimize the cost of a time-varying pumping strategy using the optimal system chosen by the first-stage optimization.

Mategaonkar & Eldho (2012c) developed a S/O model for *in-situ* bioremediation using meshfree PCM and PSO (PCM-BIO-PSO). It is observed that the PCM-BIO-PSO model is an effective model for *in-situ* bioremediation of groundwater contamination. Mategaonkar & Eldho (2014) also developed a multi-objective model for PAT method for the remediation of total dissolved solids (TDS).

Luo *et al.* (2014) and Akbarnejad-Nesheli *et al.* (2016) developed a multi-objective S/O model to achieve the best *in-situ* bioremediation system design for groundwater with contaminated dissolved hydrocarbons and observed that the Pareto-optimal solution with low variability, high reliability is a potentially effective tool for optimizing multi-objective groundwater remediation problems under uncertainty.

Kumar *et al.* (2015) presented the use of a new hybrid algorithm to optimize a multi-objective function which includes the cost of remediation as the first objective and residual contaminant at the end of the remediation period as the second objective. The hybrid algorithm was formed by combining the methods of differential evolution (DE), genetic algorithms, and simulated annealing (SA). The results obtained from the hybrid algorithms were compared with DE, non-dominated sorting genetic algorithm (NSGA II) and SA. It is found that the proposed hybrid algorithm can provide the best solution. Ehsai *et al.* (2017) used a BIOPLUME III simulation model coupled with a NSGA-II-based model for optimal design of *in-situ* groundwater bioremediation system. The multi-objective optimization model tries to minimize: (1) cost; (2) sum of contaminant concentrations that violate standard; (3) contaminant plume fragmentation. The simulation period is divided into smaller time intervals for more efficient optimization aiming to minimize the contaminant plume fragmentation.

In this study, an S/O model is presented using MFree PCM for the simulation of coupled groundwater flow and transport and PSO for getting the optimal solution for efficient *in-situ* bioremediation. Two strategies of cost minimization with number of wells and cost minimization with time of remediation are considered with different population sizes. The proposed model is applied to a case study to demonstrate the applications of the developed model.

## GOVERNING EQUATIONS AND BOUNDARY CONDITIONS

*et al.*1999). The following scheme of equations (Freeze & Cherry 1979; Borden & Bedient 1986; Mategaonkar & Eldho 2012c) can be used in the simulation of biodegradation:

Here, *h* is the piezometric head; is the specific yield; t is time; and are the velocities in the *x* and *y* directions, respectively; and are the hydraulic conductivities in the *x* and *y* directions. respectively; *Q _{w}* is flow rate from the well and

*q*is the volume rate of steady uniform recharge per unit area per unit thickness of the aquifer. Actual velocity is obtained as and ; is the porosity;

_{s}*c*is contaminant concentration; is the total microbial concentration; is maximum contaminant utilization rate per unit mass of microorganisms; is the contaminant half saturation constant;

*O*is oxygen concentration; is oxygen half saturation constant; is the concentration of microbes; is the time interval being considered;

*F*is the ratio of oxygen to contaminant consumed; is the retardation coefficient for the contaminant; are dispersion coefficient [L

^{2}T

^{−1}] in

*x*and directions, respectively.

## MESHFREE FORMULATION FOR TWO-DIMENSIONAL TRANSPORT AND OXYGEN EQUATIONS

*n*is the number of nodes in the support domain and is the multi-quadric–radial basis (MQ-RBF) shape function (Liu & Gu 2005). The shape function is given as: where

*x*and

*y*are the co-ordinates of the point of interest in the support domain; are the co-ordinates of node in the support domain; ; is the shape parameter and is the nodal spacing in the support domain. The first and second derivatives of the shape function with respect to

*x*and

*y*are calculated as given in Mategaonkar & Eldho (2012a; 2012b). Forward finite difference scheme is adopted for time discretization. Therefore, from Equations (1), (3), and (4) we get (Mategaonkar & Eldho 2012c): where is global matrix of shape function; is global matrix of first derivative of shape functions with respect to

*x*; is global matrix of second derivative of shape functions with respect to

*x*; is global matrix of first derivative of shape functions with respect to

*y*; and is global matrix of second derivative of shape functions with respect to

*y*.

Here, is the area of support domain in which the pumping well or recharge well lies and is the global matrix of the entire source and sink terms. The basis function and its derivatives are calculated for each support domain following the Kronecker delta property and are assimilated in the global matrix for whole problem domain.

Based on the above formulation, a two-dimensional meshfree model for flow, transport, and bioremediation has been developed. The developed flow, transport, and bioremediation equations were verified using available analytical and numerical solutions (Mategaonkar & Eldho 2012a; 2012b; 2012c). Further, the flow models are coupled to get an efficient model with PCM-BIO model for bioremediation simulation. To simulate the transport of pollutant and oxygen in the subsurface in aerobic bioremediation, Equations (7b) and (7c) are solved simultaneously.

## PSO-BASED OPTIMIZATION MODEL

*pbest*and

*lbest*positions. Acceleration is biased by a random term that splits random numbers being generated for acceleration to

*pbest*and

*lbest*locations (Parsopoulos

*et al.*2001). Position of individual particles is updated as follows (Kennedy & Eberhart 1995; Parsopoulos

*et al.*2001): with velocity calculated as: where is the updated position; is the particle position; is the updated velocity; is the particle velocity; is the best ‘remembered’ individual particle position; is the best ‘remembered’ swarm position; are the cognitive and social parameters; are the random numbers between 0 and 1 and

*w*is the inertia weight.

## SIMULATION OPTIMIZATION (S/O) MODEL

*in-situ*bioremediation of groundwater, simulation model PCM-BIO is coupled with the PSO model (PCM-BIO-PSO-MO). The flow chart for PCM-PSO-BIO is shown in Figure 1. Optimal cost is attained by considering multi-objectives of optimal injection rate of oxygen, number of injection wells, and time of remediation. The objective function is given as (Minsker & Shoemaker 1998): where is the injection rate in m

^{3}/day; is the relative cost coefficient;

*N*is the number of injection wells; and

*i*is the node number. PSO parameters are initialized and the objective function is assessed. The equations are solved for contaminant concentration and oxygen with the constraints. If the termination criteria are encountered, the simulation results post-processed, otherwise the particles and swarm best values, velocities and positions of particles are restructured, and the objective function is recalculated.

The purpose of the present study is to get optimal injection rates for the wells for minimization of the total cost of remediation along with remediation time and number of wells. Concentration distribution is analyzed for different strategies for the entire remediation period.

Here, two scenarios of optimization are considered: cost minimization with number of wells and cost minimization with time of remediation.

### Scenario A

The aim of optimization for this scenario is to identify bioremediation designs that minimize cost of remediation and the number of injection wells by fulfilling the limitations on the specified contaminant concentration. The result causes the least number of wells for clean-up with the consideration of the plume position. Thus, the objective function contains the oxygen injection rate and number of wells to optimize.

### Scenario B

The purpose of optimization for this scenario, is to identify bioremediation designs which reduce cost and time of remediation on the selected pollutant concentration. The solution gives the minimum time of remediation for different population sizes and thereby the remediation cost is reduced.

The objective function comprises the injection rate of water containing oxygen, time of remediation, and the number of wells. Here, the cost of setting up of the injection well is considered as fixed and hence it is not considered in the objective function. Based on the above formulation, a multi-objective model PCM-BIO-PSO-MO is developed.

The oxygen content in injected water is fixed and hence cannot be considered as a decision variable. In this study, aerated oxygen containing 8 mg/L of oxygen concentration is considered (Hinchee *et al.* 1987; Minsker & Shoemaker 1998). Further, it is assumed by the model that the microorganisms are present in substantial amounts in the aquifer for bioremediation (Minsker & Shoemaker 1998). The rate of decay is assumed to be equal to the rate of biomass growth ensuring the indigenous nature of biomass.

## MODELING PROCEDURE

A step-by-step procedure for the PCM-BIO-PSO-MO model is given below (Figure 1):

- 1.
Input all physiological and hydrological parameters.

- 2.
Discretize the aquifer into equidistant nodes in x and y directions. Consider a rectangular support domain for each node.

- 3.
Apply PCM-BIO-PSO-MO model for three different cases:

- a.
cost optimization with number of wells;

- b.
cost optimization with time of remediation;

- c.
cost optimization with number of wells and time of remediation.

- a.
- 4.
The model is run to reach the contamination level to 5 ppb or less.

The developed PCM-BIO-PSO-MO has been applied to a hypothetical problem and verified as a single-objective model and found to be satisfactory (Mategaonkar & Eldho 2012c).

## CASE STUDY

The application of the PCM-BIO-PSO model is investigated by considering a case study which is like a field problem (Figure 2). The primary components of the contamination are identified as TCE, originating from a nearby laundry. From the pumping and flow rate, the aquifer parameters are determined, and the location of contaminant plume is identified. The average groundwater velocity in the area is about 0.1 m/day. The aquifer is spread over approximately 1,350 m × 1,000 m with a thickness of 100 m (Figure 2). The maximum contamination concentration at the hotspot is about 70 μg/L. The adopted method should remediate the plume to a drinking water standard of 5 ppb (Hazen & Fliermans 1995). The hydrogeological and physical parameters for this case study are given in Table 1 (Minsker & Shoemaker 1998).

Parameter . | Value . |
---|---|

Contaminant | Tri-chloro-ethylene (TCE) |

Size | 1,350 m × 1,000 m |

Thickness of aquifer | 100 m |

Hydraulic conductivity | 4.4 × 10^{−3} m/s |

Porosity | 0.25 |

Retardation factor | 1 |

Longitudinal dispersivity | 9 m |

Transverse dispersivity | 4.5 m |

Substrate half-velocity coefficient | 49.6 mg/L |

Oxygen half-velocity coefficient | 1 mg/L |

Maximum gross specific growth rate | 6.48/d |

Ratio of oxygen to substrate | 3 |

Maximum injection hydraulic head () | 100 m |

Minimum injection hydraulic head () | 80 m |

Water quality standard, | 5 ppb |

Parameter . | Value . |
---|---|

Contaminant | Tri-chloro-ethylene (TCE) |

Size | 1,350 m × 1,000 m |

Thickness of aquifer | 100 m |

Hydraulic conductivity | 4.4 × 10^{−3} m/s |

Porosity | 0.25 |

Retardation factor | 1 |

Longitudinal dispersivity | 9 m |

Transverse dispersivity | 4.5 m |

Substrate half-velocity coefficient | 49.6 mg/L |

Oxygen half-velocity coefficient | 1 mg/L |

Maximum gross specific growth rate | 6.48/d |

Ratio of oxygen to substrate | 3 |

Maximum injection hydraulic head () | 100 m |

Minimum injection hydraulic head () | 80 m |

Water quality standard, | 5 ppb |

A specific quantity of oxygen is infused into the aquifer to start bioremediation. Decision variables do not influence the oxygen concentration in the injection of water. It is determined by the user quantified values. For the aquifer considered, 9 m and 4.5 m are the dispersivity values in longitudinal and transverse directions. The southern and northern boundaries are not flux boundaries. Constant head at the western boundary is 92.5 m and at the eastern boundary is 90.00 m.

### Application of the PCM-BIO-PSO-MO model

The case study, as mentioned above, is analyzed with the model. For the PCM-BIO-PSO model, 588 nodes are considered (Figure 3) with value as 3 and . The value of is taken as 150, day and for every square support domain nine nodes are considered. The nodal assembly for the PCM-BIO-PSO-MO model and the position of injection wells are shown in Figure 3.

The injection wells are provided at nodes 144, 183, 198, 227, and 239. The injection rates vary from 200 to 2,000 m^{3}/d for the simulation. Further parameters applied are and **.** For the S/O model, the hydraulic heads of every node ought to be fewer than and more than (Table 1). is selected to guarantee that the rise in the hydraulic head is equivalent to the most possible hydraulic head.

Original contaminant strength in the plume ranges from 0.5 to 70 μg/L with the peak concentration (>50 μg/L) at nodes 144, 183, 198, 227, and 239 (Figure 3). The primary biomass concentration at every node is assumed to be 0.005 mg/L. The plume is assumed to be enclosed in the oxygen concentration of 3 mg/L while at the center oxygen strength is assumed to be exhausted. The adopted method should remediate the plume to a drinking water standard of 5 ppb. The concentration levels in the 16 investigated observation wells are given in Table 2.

Observation well . | Contamination level (μg/L) . | Observation well . | Contamination level (μg/L) . |
---|---|---|---|

B1 | 6.3 | B9 | 48.9 |

B2 | 5.9 | B10 | 6 |

B3 | 44.1 | B11 | 10.5 |

B4 | 12.8 | B12 | 18.9 |

B5 | 0.9 | B13 | 40.6 |

B6 | 66.7 | B14 | 14.1 |

B7 | 22 | B15 | 8.2 |

B8 | 0.7 | B16 | 6.5 |

Observation well . | Contamination level (μg/L) . | Observation well . | Contamination level (μg/L) . |
---|---|---|---|

B1 | 6.3 | B9 | 48.9 |

B2 | 5.9 | B10 | 6 |

B3 | 44.1 | B11 | 10.5 |

B4 | 12.8 | B12 | 18.9 |

B5 | 0.9 | B13 | 40.6 |

B6 | 66.7 | B14 | 14.1 |

B7 | 22 | B15 | 8.2 |

B8 | 0.7 | B16 | 6.5 |

The PCM-BIO-PSO model is applied to remediate the groundwater by *in-situ* bioremediation. Three different cases are considered by applying the PCM-BIO-PSO with the cost and number of injection wells, cost and time of remediation, and cost and both number of wells and time of remediation. In all the cases, the three population sizes for PSO considered are 50, 75, and 100.

The purpose of the bioremediation is to decrease the contaminant plume up to 5 ppb; hydraulic heads at every node ought to be lower than and higher than . The assumed quantities of and are shown in Table 1. is selected in such a manner that it does not rise above the top of the aquifer whereas is selected to make sure that the rise in the hydraulic head is equivalent to the maximum possible injected hydraulic head (Minsker & Shoemaker 1998).

With the temporal interval of 1 day, the model is run for 100 days. The rate of aerated oxygen is taken as $0.56 (INR 28 approx.) and the rate of injection/extraction is taken as $14.41 (INR 720 approx.) per per month (Minsker & Shoemaker 1998).

#### Analysis 1

For the first case, three (at nodes 144, 183, and 198), four (at nodes 144, 183, 198, and 227), and five injection wells (at nodes 144, 183, 198, 227, and 239) are considered for remediation.

For this case, three population sizes of 50, 75, and 100 are considered. It is observed that the optimal solution for all the cases is achieved for a four wells' scenario. The cost variation with number of wells is shown in Table 3. The total cost ranges from Rs. 5.59 × 10^{7} to Rs. 9.24 × 10^{7} for a population size of 100 and 100 days.

Number of wells . | Population size . | Time of remediation . | ||
---|---|---|---|---|

70 days . | 100 days . | 200 days . | ||

Cost (Rs.) . | ||||

3 | 50 | 7.28 × 10^{7} | 9.15 × 10^{7} | 1.72 × 10^{8} |

75 | 6.82 × 10^{7} | 8.63 × 10^{7} | 1.70 ×10^{8} | |

100 | 6.87 × 10^{7} | 6.78 × 10^{7} | 1.57 × 10^{8} | |

4 | 50 | 7.61 × 10^{7} | 6.32 × 10^{7} | 1.31 × 10^{8} |

75 | 7.42 × 10^{7} | 6.25 × 10^{7} | 1.28 × 10^{8} | |

100 | 6.35 × 10^{7} | 5.59 × 10^{7} | 1.19 × 10^{8} | |

5 | 50 | 1.36 × 10^{8} | 9.30 × 10^{7} | 1.88 × 10^{8} |

75 | 1.35 × 10^{8} | 9.26 × 10^{7} | 1.86 × 10^{8} | |

100 | 1.35 × 10^{8} | 9.24 × 10^{7} | 1.85 × 10^{8} |

Number of wells . | Population size . | Time of remediation . | ||
---|---|---|---|---|

70 days . | 100 days . | 200 days . | ||

Cost (Rs.) . | ||||

3 | 50 | 7.28 × 10^{7} | 9.15 × 10^{7} | 1.72 × 10^{8} |

75 | 6.82 × 10^{7} | 8.63 × 10^{7} | 1.70 ×10^{8} | |

100 | 6.87 × 10^{7} | 6.78 × 10^{7} | 1.57 × 10^{8} | |

4 | 50 | 7.61 × 10^{7} | 6.32 × 10^{7} | 1.31 × 10^{8} |

75 | 7.42 × 10^{7} | 6.25 × 10^{7} | 1.28 × 10^{8} | |

100 | 6.35 × 10^{7} | 5.59 × 10^{7} | 1.19 × 10^{8} | |

5 | 50 | 1.36 × 10^{8} | 9.30 × 10^{7} | 1.88 × 10^{8} |

75 | 1.35 × 10^{8} | 9.26 × 10^{7} | 1.86 × 10^{8} | |

100 | 1.35 × 10^{8} | 9.24 × 10^{7} | 1.85 × 10^{8} |

#### Analysis 2

For the second case, three different periods of remediation of 70, 100, and 200 days are considered. For this case, also three population sizes (50, 75, and 100) are considered. Cost variation with time for all the population sizes is shown in Figure 4. From Table 3 it is observed that the optimal solution is obtained for 100 days. The total cost for all scenarios range from Rs. 5.59 × 10^{7} to Rs. 1.35 × 10^{8} for all time periods.

#### Analysis 3

For the third case, both the number of wells and time of remediation were optimization functions along with cost. The results are shown in Table 4. It is observed that even if the number of injection wells is reduced, the cost is not reduced because of high injection rate, and cost is not reduced even if the time of remediation is reduced.

No. of wells . | Well locations (nodes) . | Time of remediation (days) . | Time of remediation (days) . | ||||
---|---|---|---|---|---|---|---|

70 . | 100 . | 200 . | 70 . | 100 . | 200 . | ||

Rate of injection (m^{3}/s)
. | Cost (Rs.) . | ||||||

3 | 144 | 1,582.44 | 750.98 | 624.91 | 6.87 × 10^{7} | 6.78 × 10^{7} | 1.57 × 10^{8} |

183 | 1,193.79 | 675.61 | 592.68 | ||||

198 | 887.86 | 516.69 | 320.99 | ||||

4 | 144 | 1,218.45 | 979.75 | 804.14 | 6.35 × 10^{7} | 5.59 × 10^{7} | 1.19 × 10^{8} |

183 | 981.76 | 750.34 | 657.34 | ||||

198 | 868.19 | 664.18 | 571.24 | ||||

227 | 753.44 | 577.78 | 497.01 | ||||

5 | 144 | 1,107.89 | 902.56 | 753.12 | 1.35 × 10^{8} | 9.24 × 10^{7} | 1.85 × 10^{7} |

183 | 958.78 | 813.22 | 639.26 | ||||

198 | 833.11 | 719.67 | 540.77 | ||||

227 | 720.19 | 657.76 | 489.23 | ||||

239 | 672.72 | 603.56 | 402.77 |

No. of wells . | Well locations (nodes) . | Time of remediation (days) . | Time of remediation (days) . | ||||
---|---|---|---|---|---|---|---|

70 . | 100 . | 200 . | 70 . | 100 . | 200 . | ||

Rate of injection (m^{3}/s)
. | Cost (Rs.) . | ||||||

3 | 144 | 1,582.44 | 750.98 | 624.91 | 6.87 × 10^{7} | 6.78 × 10^{7} | 1.57 × 10^{8} |

183 | 1,193.79 | 675.61 | 592.68 | ||||

198 | 887.86 | 516.69 | 320.99 | ||||

4 | 144 | 1,218.45 | 979.75 | 804.14 | 6.35 × 10^{7} | 5.59 × 10^{7} | 1.19 × 10^{8} |

183 | 981.76 | 750.34 | 657.34 | ||||

198 | 868.19 | 664.18 | 571.24 | ||||

227 | 753.44 | 577.78 | 497.01 | ||||

5 | 144 | 1,107.89 | 902.56 | 753.12 | 1.35 × 10^{8} | 9.24 × 10^{7} | 1.85 × 10^{7} |

183 | 958.78 | 813.22 | 639.26 | ||||

198 | 833.11 | 719.67 | 540.77 | ||||

227 | 720.19 | 657.76 | 489.23 | ||||

239 | 672.72 | 603.56 | 402.77 |

### Sensitivity analysis

*N*is the total number of nodes in the aquifer; is the contaminant concentration at the end of remediation period; is the optimum contaminant concentration to be achieved at the end of remediation period (5 ppb). The divergence is calculated and plotted with different parameters on semi-log plot by changing the parameters one by one keeping other parameters constant. The attempt is to reduce the divergence to its minimum to achieve the minimum concentration of 5 ppb (Minsker & Shoemaker 1998) after remediation for 100 days. The model is tested for the time steps of 1 day, 2 days, 5 days, and 10 days for 100 days. It is found that the model results improve with a lesser time step. The results are shown in Figure 5. Second, the model is tested for inter-nodal distances, , , and . It is observed that the divergence is reduced with the lesser and values. The results after two years of remediation are shown in Figure 6. Finally, the model is tested for different values of 25, 50, 75, 100, and 150. The error reduces with the increasing value of . It indicates that divergence reduces with the greater number of nodes in the support domain. The results after two years are shown in Figure 7. The results indicate that the model is sensitive for all the above-mentioned parameters. However, to achieve the optimum value of contaminant concentration to 5 ppb after bioremediation, the model gives satisfactory results for day, , and for 100 days of remediation with optimal computational efforts.

## CONCLUSIONS

The simulation-optimization models help to design optimal remediation of groundwater. In this study, a PCM-based meshfree model is developed for bioremediation simulation. Also, a PSO-based optimization model is developed for optimization. PCM is a very simple method to work with and apply. With the appropriate selection of shape parameters, the model provides efficient results. The PCM models are further coupled with PSO-based optimization techniques (PCM-BIO-PSO) for *in-situ* bioremediation of groundwater contamination to get an efficient and optimal solution to the bioremediation problem.

In this study, the PCM-BIO-PSO model is applied to a case study similar to a field problem for *in-situ* bioremediation with three scenarios, namely, cost and number of injection wells, cost and time of remediation, cost and number of wells and time of remediation together.

The total cost for three wells for 70 days, 100 days, and 200 days is Rs. 6.87 × 10^{7}, Rs. 6.78 × 10^{7}, and Rs. 1.57 × 10^{8}, respectively. The total cost for four wells for 70 days, 100 days, and 200 days is Rs. 6.35 × 10^{7}, Rs. 5.59 × 10^{7}, and Rs. 1.19 × 10^{8} respectively. The total cost for five wells for 70 days, 100 days, and 200 days is Rs. 1.35 × 10^{8}, Rs. 9.24 × 10^{7}, and Rs. 1.85 × 10^{8}, respectively. In all the cases, it is observed that the PCM-BIO-PSO-MO model gives the optimal cost of Rs. 5.59 × 10^{7} for four wells, 100 days, and 100 population size for four wells' scenario.

It is observed that as the number of injection wells increases the cost is increased, but at the same time, the remediation time is reduced. The injection rate increases for less time of remediation. Also, the problem solution is optimized with the increasing population size. The proposed methodology can be effectively used in the *in-situ* bioremediation of contaminated sites.

## ACKNOWLEDGEMENTS

The authors are grateful to Dr Ulf Mohrlok, Karlsruhe Institute of Technology for his support.

## REFERENCES

*.*

*.*

*.*

.Patent no. 5384048, Assignee: The United States of America.Journal of Irrigation and Drainage Engineering, Issue date: 24 January 1995.