Researchers MohammadReza Najarzadegan and Ramtin Moeini developed and validated an integrated simulation-optimization framework that couples an improved artificial bee colony (IABC) algorithm with EPANET, the widely used hydraulic and water-quality modeling software from the United States Environmental Protection Agency. Rather than introducing a brand-new optimization technique, the authors’ contribution lies in formulating and testing a complete inverse problem: simultaneously estimating the contamination source location, the injection time, and the injection concentration from residual-chlorine observations. These three variables are deeply intertwined through contaminant transport, network hydraulics, and chlorine decay, which makes the search space nonlinear and non-convex, meaning that small changes in one variable can dramatically alter the optimal values of the others.

The artificial bee colony algorithm draws its inspiration from the foraging behavior of honeybees, in which employed bees explore food sources, onlooker bees select promising sources based on shared information, and scout bees search for new regions. In the improved version used in this study, the relative contribution of employed and onlooker bees is adjusted to control how search effort is divided between broad exploration of the solution space and intensive refinement of promising candidates. This balance matters for contamination-source identification because the source location is a discrete variable, while injection time and concentration are continuous. Too much local search risks premature convergence on a wrong node; too much exploration prevents the algorithm from sharpening an accurate estimate of timing and concentration.

The mechanics of the framework are elegantly cyclical. In each iteration, the IABC algorithm generates a candidate contamination scenario consisting of a source node, an injection start time, and an injection concentration. These values are transferred from MATLAB to EPANET through the EPANET-MATLAB Toolkit, which triggers a full hydraulic and water-quality simulation of how the contaminant would propagate through the pipes. The simulated concentrations at the monitoring sensors are then extracted and compared with the observed values, and the objective function, defined as the sum of squared differences between observed and simulated concentrations across all sensors and time steps, is evaluated. The algorithm updates its population of candidate solutions based on these scores, and the cycle repeats until the stopping criterion is met.

Chlorine plays a central role as the diagnostic signal. Utilities maintain free chlorine between roughly 0.5 and 1 milligram per liter throughout their networks, and any contaminant entering the system reacts with and depletes that residual disinfectant. The researchers modeled this depletion using first-order reaction kinetics, in which chlorine and the reactive contaminant decay at rates proportional to their concentrations, with stoichiometric coefficients handled automatically by EPANET. Water age, the hydraulic residence time of water at each node, was also computed as supporting information, since older water tends to carry less residual chlorine and can confound the interpretation of sensor signals.

To test the framework, the team used two case studies of very different scale. The first was the third example network from EPANET, a benchmark containing 92 nodes, 117 conduits, 3 tanks, 2 pumps, and two reservoirs, analyzed over a 24-hour extended period simulation. Two contamination scenarios were examined, with injections at node 100 and node 209 respectively, each beginning at 11 a.m., coinciding with peak demand, at a nominal concentration of 50 milligrams per liter lasting 3 hours. The second case study was far more ambitious: the real water distribution network of Baharestan, Iran, comprising 3,147 consumer nodes, 3,402 pipes, 12 pressure-reducing valves, and 11 quality sensors placed according to actual operating conditions.

Each optimization problem was solved 30 independent times to account for the stochastic nature of the metaheuristic search. In the EPANET benchmark, IABC achieved mean objective-function values of 0.506 and 0.198 for the two scenarios, compared with 0.579 and 0.352 for the shuffled frog algorithm, reductions of 12.60 percent and 43.75 percent respectively. In the Baharestan network, IABC’s mean objective-function value of 6.720 beat the shuffled frog algorithm’s 7.820, a 14.07 percent reduction. When compared with genetic algorithm results from a previous benchmark study, IABC’s mean absolute concentration errors at selected observation points were 92.88 percent and 93.47 percent lower in the two EPANET scenarios. However, the authors are careful to note that Wilcoxon rank-sum tests found no statistically significant difference between IABC and the shuffled frog algorithm at the 5 percent level, so the reported advantages are descriptive rather than proof of superiority.

Robustness was probed by injecting synthetic measurement noise of plus or minus 5, 10, and 20 percent into the observations. As expected, higher uncertainty inflated the objective-function values and degraded identification accuracy, underscoring how much the framework depends on reliable, well-calibrated sensors. A sensitivity analysis over contaminant concentration and injection duration showed that the identification problem becomes harder as events grow larger or longer; for instance, raising the source concentration from 50 to 100 milligrams per liter increased the mean objective-function value in the Baharestan network from 6.930 to 7.510. Notably, in the large Baharestan network, the identified source node varied across the 30 runs, revealing that source localization was not unique under the adopted sensor configuration even when the objective-function values remained low.

The authors are candid about limitations. The validation was conducted under modeled conditions with fixed sensor numbers and locations, and the analysis did not explicitly account for sensor failures, calibration errors, communication delays, missing data, or temporal changes in demand and hydraulics that occur in real systems. The noise perturbations represent a controlled approximation of field uncertainty rather than its full complexity. Broader validation using field-monitoring data, alternative sensor configurations, and dynamically varying hydraulic conditions is recommended before the framework can be generalized to operational networks.

Even with those caveats, the study marks a meaningful step forward for water security. Most previous research on contamination in water distribution networks has concentrated on where to place sensors, how to assess vulnerability, or how to mitigate consequences after detection. Far fewer studies have attempted the full inverse problem of characterizing an unknown contamination event from water-quality observations alone. By demonstrating that a bee-inspired optimizer, coupled to a physically grounded simulation engine, can simultaneously recover source location, timing, and concentration across networks ranging from a 92-node benchmark to a city-scale system of more than 3,000 nodes, the research offers utilities a practical template for turning sparse sensor data into actionable forensic intelligence. In an era when drinking-water infrastructure faces growing threats from emerging contaminants such as microplastics, pharmaceuticals, and antibiotic-resistance pollutants, tools that can rapidly answer where, when, and how much may prove indispensable for protecting public health.

Subject of Research: Inverse identification of contamination source location, injection time, and concentration in water distribution networks using an improved artificial bee colony algorithm coupled with EPANET simulation

Article Title: Using an improved artificial bee colony algorithm to identify contamination sources in water distribution networks

Article References: Najarzadegan, M., & Moeini, R. (2026). Using an improved artificial bee colony algorithm to identify contamination sources in water distribution networks. Case Studies in Chemical and Environmental Engineering, 14, Article 101500. https://doi.org/10.1016/j.cscee.2026.101500

Image Credits: AI Generated

DOI: 10.1016/j.cscee.2026.101500

Keywords: water distribution networks, contamination source identification, artificial bee colony algorithm, EPANET, metaheuristic optimization, residual chlorine, water quality monitoring, inverse problem, shuffled frog algorithm, genetic algorithm, chlorine decay, water security

The research, carried out by Derya Deliktaş of Kütahya Dumlupınar University, Ezgi Aktar Demirtas of Eskişehir Osmangazi University, and Mert Göksüzoğlu of the SAMPA Automotive R&D Center in Samsun, Türkiye, tackles a deceptively simple question: what happens when you blend chloroprene rubber with cheaper natural rubber, and how do you find the exact formulation that delivers the properties an engineer actually needs? Blending the two elastomers is an obvious strategy for balancing cost and performance, but the chemistry is far from trivial. The proportions of natural rubber, along with the amounts of accelerators, retarders, and vulcanization agents added to the mix, all interact in complicated ways that determine how the final material behaves.

To map this complicated landscape, the team first turned to Response Surface Methodology, a statistical framework that has been a staple of industrial optimization since the mid-twentieth century. Rather than testing every possible combination of ingredients, which would require an impractical number of experiments, the researchers used a Central Composite Design to select a structured set of experimental runs. This design allows a quadratic response surface to be fitted over the formulation space, meaning that the relationship between each ingredient and the resulting material property can be captured with a relatively small number of trials. From these experiments, the team built regression models describing three key response variables: hardness, tensile strength, and the loss factor known as Tan δ.

Each of these three properties matters for a different reason. Hardness, measured in Shore A units, reflects how resistant the rubber is to indentation and deformation, which is critical for seals, mounts, and bushings. Tensile strength, expressed in megapascals, measures how much stretching force the material can withstand before breaking, a fundamental durability metric. Tan δ, the ratio of the loss modulus to the storage modulus, captures the damping behavior of the rubber, essentially how efficiently it converts vibration energy into heat. For automotive components such as engine mounts and suspension bushings, damping is what keeps noise and harshness out of the cabin, making Tan δ a crucial dynamic property.

Here lies the central conflict that makes rubber formulation so difficult: the objectives pull against each other. A formulation that maximizes tensile strength may deliver mediocre damping. One that minimizes Tan δ for superior vibration isolation may compromise strength. Improving hardness can come at the expense of flexibility. Because these trade-offs cannot be resolved by optimizing one property at a time, the researchers framed the problem as multi-objective optimization, and this is where the artificial bees enter the picture.

The Artificial Bee Colony algorithm, first developed by Turkish computer scientist Derviş Karaboga and collaborators, is a metaheuristic inspired by the foraging behavior of honeybee swarms. In the algorithm, artificial bees explore the search space in three roles: employed bees exploit known food sources, which in this context correspond to promising formulation candidates; onlooker bees choose among those sources with a probability proportional to their quality; and scout bees abandon exhausted sources and search randomly for new ones. This division of labor allows the algorithm to balance exploitation of good solutions with exploration of uncharted regions, avoiding the premature convergence that plagues simpler optimization methods. The ABC algorithm has proven effective in fields ranging from structural truss design to vehicle routing, and the Turkish research team brought it to bear on rubber chemistry.

To combine the three competing objectives into a form the algorithm could optimize, the researchers applied scalarization methods, which convert multiple objectives into a single weighted function. They compared three approaches: the classical weighted sum method, the Tchebycheff method, and the Conic Scalarization Method, a technique developed by mathematician Refail Kasimbeyli that can characterize efficient solutions even in nonconvex problems. The comparison revealed a clear winner. The Conic Scalarization Method delivered the most robust optimization performance, offering flexible trade-offs between mechanical strength and dynamic properties that the other methods could not match. This finding is significant for practitioners because the choice of scalarization technique can determine whether an optimization finds genuinely useful compromise solutions or gets stuck in inferior corners of the design space.

The practical payoff of the approach is best illustrated by the numbers. When the team prioritized minimization of Tan δ by assigning it a weight of 0.8, the optimizer produced a formulation yielding a Tan δ of just 0.076, with a hardness of 49.53 Shore A and a tensile strength of 16.23 megapascals. That is a formulation tailored for components where vibration damping is paramount. When hardness was prioritized instead, the resulting blend achieved 57.15 Shore A, a tensile strength of 19.38 megapascals, and a Tan δ of 0.198, suited to parts that must resist deformation. Prioritizing tensile strength with the same weight produced the most impressive strength figure of all: 24.86 megapascals, accompanied by 53.60 Shore A hardness and a Tan δ of 0.238. In other words, a single optimization framework can hand engineers three very different rubber recipes, each tuned to a different design requirement.

What makes this work notable beyond the rubber lab is the way it stitches together two traditions that rarely collaborate so closely. Response Surface Methodology provides statistically validated models grounded in real experimental data, while the swarm-based optimizer searches those models efficiently across thousands of candidate formulations without a single additional laboratory test. The study builds on the team’s earlier 2024 work in Polymer Bulletin, which applied response surface methodology alone to the same chloroprene and natural rubber blend system, and it reflects a broader trend in materials science of pairing designed experiments with artificial intelligence and metaheuristic algorithms. Similar hybrid approaches have recently been applied to polymer nanocomposites, machining of composites, and biodiesel production, suggesting that the RSM-plus-metaheuristic pipeline is becoming a standard tool for modern materials engineering.

For the automotive industry, the implications are concrete. Rubber component manufacturers working with chloroprene blends can now adjust their formulations to hit targeted hardness, strength, or damping values with confidence, rather than relying on trial-and-error compounding. Because natural rubber is cheaper than chloroprene, every formulation that shifts the balance toward more natural rubber while preserving performance translates directly into cost savings at industrial scale. The researchers, who acknowledge the facilities provided by the SAMPA Automotive R&D Center, describe their CSM-based multi-objective ABC approach as a versatile tool for developing cost-effective, high-performance automotive rubber components. As material costs and performance demands continue to rise together, it seems fitting that the answer to designing better rubber may have come from mimicking the collective intelligence of bees.

Subject of Research: Multi-objective optimization of chloroprene/natural rubber blend mechanical and rheological properties using response surface methodology and the artificial bee colony algorithm

Article Title: Response surface-guided multi-objective optimization of CR/NR blend mechanical and rheological properties using artificial bee colony algorithm

Article References: Deliktaş, D., Demirtas, E. A., & Göksüzoğlu, M. (2026). Response surface-guided multi-objective optimization of CR/NR blend mechanical and rheological properties using artificial bee colony algorithm. Polymer Bulletin, 83(11), Article 618. https://doi.org/10.1007/s00289-026-06653-6

Image Credits: AI Generated

DOI: 10.1007/s00289-026-06653-6

Keywords: chloroprene rubber, natural rubber, rubber blends, vulcanization, response surface methodology, artificial bee colony algorithm, multi-objective optimization, conic scalarization method, tensile strength, hardness, Tan delta, automotive rubber components