Abstract
This paper proposes a control strategy to stabilize a reverse osmosis desalination system against hydraulic shocks with enhancing productivity and sustainability. First, the effects of hydraulic transients on water quality have been reviewed. The transient waves are approximated by sinusoidal functions so that their effects are incorporated into the controlled system as external disturbances. Next, the active control is implemented based on the adaptive super-twisting (STW) sliding mode control (SMC) algorithms. Then, the robust performance is guaranteed whenever the sliding variables reach the sliding surfaces in finite time despite disturbances. The STW SMC scheme is to eliminate the chattering problems for protecting the valves and to improve the convergence precision for water production. The control gains are adaptable to enable formation of an effective controller for dealing with large disturbances such as water hammer during desalination process. The simulation results reveal the superior performances on controlling water product, while eliminating shock waves. Especially, the effect of hydraulic shocks has been dramatically attenuated, hence the plant components are protected to avoid fracture. Finally, the robust stability and performance of the desalination plants are guaranteed against large disturbances to ensure the population with quality water as well as system sustainability.
HIGHLIGHTS
Hydraulic transients and dynamical analysis have been provided for desalination plants.
Adaptive sliding mode control is applied for suppression of hydraulic shocks.
Active controller provides robust performance and stability ensuring water production efficiency and product quality.
Complete system is implemented with prolonging life expectancy of desalination equipment by sliding mode control.
INTRODUCTION
Generally, water desalination is the process of removing dissolved salts and mineral components from target substances in order to attain quality water for animal consumption, irrigation, and human use. Water management based on desalination technology has become the solution for world water scarcity for many decades. Especially, desalinated water is the main water source for meeting the water demands in water-scarce regions. Most of the modern interests in desalination have focused on developing effective ways of providing fresh water for human consumption in regions where the availability of fresh water is extremely limited. Desalination plants have been designed to consistently produce fresh water from salty feed to system output in sufficient quality as economically as possible. The major technologies utilized for desalination processes include reverse osmosis (RO) and multistage flash distillation (MSF). Recently, the RO has gained more dominance, particularly for brackish water treatments (Alatiqi et al. 1999). This desalination technique has also been used in diverse processing industries such as chemical, nuclear, biotechnology, and petroleum (Jamal et al. 2004). The fundamental RO principle is based on the transport mechanisms of solvent and solute through semi-permeable membranes for water treatments, using high pressure. Yet, this system is of considerable complexity in terms of its biological and physical aspects. A number of desalination researches have been conducted to realize a variety of water treatment mechanisms in existing technologies. Generally, the RO transport mechanisms can be divided into four main types (Ghernaout 2017): the wetted surface (Reid & Breton 1959), the solution–diffusion (Lonsdale et al. 1965), the sieve (Banks & Sharples 1966), and the preferential sorption–capillary flow mechanism (Sourirajan 1986). Consequently, many RO models have been presented in literature, including the models of Senthilmurugan et al. (2005), Oh et al. (2009), Sundaramoorthy et al. (2011), and Chaabene et al. (2018). Some transfer functions have been also proposed for RO control purposes such as those of Robertson et al. (1996), Alatiqi et al. (1999), Riverol & Pilipovik (2005), Bartman et al. (2009), and Chaaben et al. (2011). A better quantitative understanding in terms of process dynamics is thus required to enhance plant performance. Due to uncertainties and disturbances in the real-time operations, the design of robust monitoring systems of a large RO process is one of the most challenging jobs in water desalination plants. The plant uncertainties can be parametric uncertainty due to modeling mismatch, unmodeled dynamics, and simplification or parameter changes because of disruptions such as concentration polarization, feed concentration, fouling and scaling. The external disturbances are possibly feed-concentration change, leakage or more severe disturbances caused by, for example, hydraulic transients. Therefore, a robust dynamic controller is required to cope with a RO desalination plant confronting these challenges to increase productivity and prolong the life of plant components (Sobana & Panda 2011). RO literature shows that although there are more and more challenges in controlling RO systems, current controllers are mostly proportional-integral-derivative (PID) based (Alatiqi et al. 1989; Kim et al. 2009; Rathore et al. 2013; Chithra et al. 2015; Ncube & Inambao 2019; Patnana et al. 2020), and model predictive control (MPC) (Abbas 2006; Bartman et al. 2009; Sobana & Panda 2014; Rivas-Perez et al. 2019; Kargar & Mehrad 2020; Mehrad & Kargar 2020). Furthermore, there are very few studies considering robust control synthesis of desalination plants under parametric uncertainties and disturbances such as water hammering. PID, MPC, and sliding mode control (SMC) algorithms are some of the best for regulating perturbed dynamical systems with uncertainties and disturbances in many industrial applications (Edwards & Spurgeon 1998; Utkin et al. 1999). The chattering problem of conventional SMC has been solved by high order SMCs (Adjoudj et al. 2011; Benbouzid et al. 2014; Patnaik et al. 2016). In fact, the super-twisting (STW) algorithm is considered as a second-order SMC. The concept of this algorithm is to generate a continuous control function which ensures the convergence of the sliding surface and its derivative in finite time. The standard STW control algorithm requires knowledge of the disturbance bound which is not easy to be estimated in practice (Boubzizi et al. 2018). Despite the robustness and applicability of SMC, there are still very few works reported to regulate RO systems (Vrkalovic et al. 2018; Zebbar et al. 2019) and none of them account for the effects of hydraulics transients such as water hammer. An important consideration in the safe operation of water treatment plants is the suppression of undesirable hydraulic transients, which are also related to water hammer or pressure surge. In the late 1960s, transient flow analyses using method of characteristics (MOC) and digital computers gained favor over graphical methods. In reality, the desalination plants suffer from more frequent and severe environmental variations such as hydraulic and influent flow shock. Most commonly, hydraulic shock occurs when a column of moving non-compressible liquid is suddenly stopped in line pressure. The excessive pressure surges may cause collapse of pipeline or damage of hydraulic equipment in the water treatment system (Bird et al. 2002). Larger pressure waves can certainly cause serious safety problems for desalination plants. In many industrial applications, high-pressure waves can burst pipes, valves, and open joints, flanges, etc., and low pressures can collapse pipes or suck contaminated water into pipes (Chaudhry 2014; Liu & Simpson 2018). Hydraulic shock loads should be kept within the prescribed limits to protect the system components. Traditional passive devices against pressure pulsations and hydraulic shocks do not fully provide emergency protections against system failure due to pipe collapse or bursting (Schmitt et al. 2006). To mitigate severe hydraulic surge effects, the active controller can be utilized to minimize pressure fluctuations. Hydraulic transient and its effects has become a major area of concern for many researchers, but few have been involved in actively controlling water surges of desalination plants (Alidai & Pothof 2016; Phuc et al. 2017). The goals are to robustly regulate the product water flow and system pressure to track some set-points and, especially, to stabilize the water productivity under water shock loads. In this paper, an adaptive STW SMC approach (Shtessel et al. 2010) is introduced for the desalination system. This algorithm can adapt control parameters to guarantee the convergence of the sliding surface to zero with unknown uncertainties against hydraulic shocks. The STW SMC strategy is effective for controlling of RO systems in normal conditions. However, in case of hash conditions such as water hammer, the robust stability of the RO system controlled by the fixed-gains of STW SMC may not be achieved. Therefore, adaptable control gains are integrated to guarantee the robust performance of the control scheme. Another advantage of this SMC-based controller is that it is readily applicable for deploying into the RO control hardware. The simulation results show that the active control algorithm provides superior performance and is perfectly fitting to water monitoring systems coping with fluid shocks. Finally, the active controller can effectively monitor all the components of desalination plants, taking automatic compensation actions when hydraulic pressure surges occur, with resulting advantages in terms of effective automation and general plant management.
MATERIALS AND METHODS
Hydraulic transients in desalination plants
Basic components of a desalination system are briefly described considering hydraulic interactions among the system devices. In, the mass and momentum transfer under hydraulic shocks are not well understood for water treatment systems. The disadvantages of an RO system are the difficulties of obtaining rigorous mathematical models of desalination process, which accounts for several operating factors such as feed temperature, concentration polarization, and fouling under hydraulic surge waves. The RO controlled system is developed by Bartman et al. (2009), in which the overall mass balance and local energy balances are exploited around the valves and combined with additional terms. The schematic diagram of the RO model is depicted in Figure 1. In this simplified system, feed water enters the system and is pressurized by the high-pressure pump to pass through the membranes to the product water. The feed water from the water line valve is first pushed through the pre-filters which protect the semi-permeable RO membranes from physical and chemical damage. The RO desalination system incorporates a high-pressure pump to create the necessary level of force to counter the osmotic pressure (Δπ). The remaining water with a high concentration is rejected through the concentration valve. This concentrated flow is controlled by the concentration valve using valve resistance Rvc. The bypass valve is to control the bypass flow by the valve resistance Rvb. The product water flow Fp and the system pressure Ps are selected as controlled outputs in this study. The feed stream velocity is kept constant at 10 m/s. The product water and system pressure will change their values according to the opening of the control valves. A lack of proper plumbing in the feed water system with pump operations will result in mechanical damage to the RO membrane. Suitable methods for active control systems must be considered to avoid such severe unacceptable flow situation. During operations of dynamical fluid flow systems, such as the desalination process, the unsteady flow conditions will be inevitably encountered. Transient events in the desalination plants usually occur because of actions at pumps and control valves that are detrimental to vulnerable hydraulic components such as filtration and membranes. When hydraulic transients occur too quickly, they induce a rapid change in flow rate within desalination plants and cause potential pressure surges that could lead to detrimental operating conditions. High-pressure upsurges can lead to desalination system failure and excess leakage, where low-pressure down surges can create vacuum conditions, and system equipment and devices can collapse or burst. What started as a hydraulic surge has now evolved into a vapor condition, itself causing its own bubble collapse hydraulic shock. The primary objectives of transient analysis are to determine the values of transient pressures that can result from flow control operations and to establish the design criteria for system equipment and devices so as to provide a tolerable level of protection. Careful considerations are thus required in the system design stages to make sure that the unsteady fluid operations do not give rise to unacceptable flow and excessive pressure transient conditions.
Dynamical analysis
A hydraulic shock typically occurs in pipes after a valve is shut off suddenly. Specifically, the quick starting or stopping of pumps can create damaging pressure spikes in the valves known as water hammer. Since liquids are incompressible, when there is a sudden change in fluid velocity, the kinetic energy of the moving fluid is immediately converted into potential energy, causing waves of pressure and flow velocity back to the fluid source. A flow velocity greater than 1.5 m/s, a valve closing in less than 1.5 s, and a high operating pressure will cause, water hammer to occur and the transient pressure can be as high as five times the initial working pressure. Therefore, water hammer phenomenon in RO systems must be sufficiently considered for RO plant design and to avoid plant damages or failure.

Water hammer wave profiles: (a) transient pressure and (b) concentration flow velocity.
Water hammer wave profiles: (a) transient pressure and (b) concentration flow velocity.

From Figure 3, it can be seen that the directions of the flows through the valves are alternately inverted, causing strong disturbances to the dynamic behaviors of the RO system. With no surge protection, transient wave profiles are unacceptably high, and if the reduced pressure is below the vapor pressure, a steam bubble will form (a large cavitation). It is noted that fluid fluctuation effects may disturb overall operations of the hydraulic system. This potentially dangerous (or undesirable) phenomenon can be mitigated with suitable design practice or active process controller. The oscillation functions in Equation (8) will be incorporated into the numerical simulation as external disturbances to check the system performance. Transient behaviors were be eliminated and maximum pressures were reduced until they coincided with steady-state pressures. The identification and calculation of pressures, velocities, and other abnormal behaviors resulting from hydraulic transients make possible the effective use of various control strategies in order to minimize energy and fluid losses or to improve the system capacity and the water quality.
Adaptive SMC synthesis



In the reaching part of the controller, the action w is designed to deal with disturbances and uncertainties to guarantee the reachability of the sliding surfaces in Equation (15). The control synthesis of w is based on the approach proposed by Shtessel et al. (2010), where the following theorem provides sufficient condition to guarantee the stability and robustness of the RO system. Again, the external disturbances such as hydraulic transients can create serious consequences for water utilities if not properly recognized and addressed by proper analysis.



The proof is completed for the state variable . Similarly, the same justification can be made to prove the stability of the state variable x2. The adaption laws for control gains utilize the current information on the sliding surface to adjust the control input in real time.
RESULTS AND DISCUSSION
Transient events are also significant for water quality as well as productivity in the desalination plants. Hydraulic shocks in particular can cause damages or fracture to the related components of desalination plants. Thus, the hydraulic transient studies with water hammering should be included in the design stage for RO plants. The numerical simulation is carried out at working steady-state (vr = 4, vb = 1, and vp = 5). The current product flow rate and the system pressure are 2,286 L/h and 35.69 bar, respectively. The set-points are 2,000 L/h and 40 bar, respectively. Since closing the valve suddenly could cause severe hydraulic shock effects, the active controller must be manipulated carefully to minimize flow fluctuations. The adaptive STW SMC controller has been successfully designed for the RO desalination system with the selected control parameters as follows:
,
,
,
, and
. The proposed SMC scheme with adaption is intended to minimize the control activities along with a finite-time convergence against water hammer. Figure 4 shows the step responses of the controlled variables. It can be seen that the product flow rate Fp and the system pressure Ps can reach new set-points with settling time within approximately 1 s, which is better than the other results reported in RO literature as summarized in Table 1. It is worth noting that previous studies do not provide clarity on water hammer. As depicted in Figure 4, there is a slight overshoot in the flow rate and no overshoot in the pressure channel. The fast responses demonstrate that the active controller can quickly drive the system to desired steady-states under some variations.
Comparison of control strategies performance
Control algorithm . | PID (Gambier 2011) . | MPC . | SMC (Vrkalovic et al. 2018) . | STW SMC (Zebbar et al. 2019) . | Proposed Adaptive STW SMC . |
---|---|---|---|---|---|
Set-point (L/h) | 35 | 350 | 227 | NA | 2,000 |
Settling time (s) | 1.5 | 4 | 1.5 | 1.5 | 1 |
Control algorithm . | PID (Gambier 2011) . | MPC . | SMC (Vrkalovic et al. 2018) . | STW SMC (Zebbar et al. 2019) . | Proposed Adaptive STW SMC . |
---|---|---|---|---|---|
Set-point (L/h) | 35 | 350 | 227 | NA | 2,000 |
Settling time (s) | 1.5 | 4 | 1.5 | 1.5 | 1 |
Figure 5 shows the tracking errors and the finite-time convergence to the sliding surfaces. It is noted that the set-point tracking in the RO plant is guaranteed to converge to zero with in a time limit of approximately 1 s. When the system states reach the sliding surfaces, they will converge to the set-points and be kept there until new desired values are set.
One of important issues in the SMC algorithm is the control actions. The classic SMC uses a discontinuous control action. If the control action is discontinuous in the state variables, it causes chattering problem which is undesirable in the practical implementations. In this work, the control signals are plotted in Figure 6 by implementing smooth control input. It can be observed that there is almost no chattering and the signals are not saturated due to the characteristics of the STW algorithm. This control scheme also guarantees the fast convergence of the sliding surfaces as illustrated in Figure 5. It helps protect against machine damages and prolongs the lifespan of the operating valves.
For the control synthesis, the most interesting concern is the ability of disturbance attenuation, especially large disturbances such as water hammer, and more generally, fluid hammer. It is noted that this fluid hammer in liquid line can develop peak pressures several times greater than normal working pressures. In this simulation, the water hammer occurring at the permeate side is considered as a disturbance to the RO system. There are always pressure-reducing valves (PRV) installed in RO systems but they are not enough to eliminate hydraulic shock effects. As discussed earlier, the hammer transient pressure has a damping wave profile and the magnitude of the first peak pressure wave is extremely high. Vacuum pressure also occurs, and both phenomena is very harmful to the RO system. Therefore, not only high positive pressure waves but also negative ones must be regulated to avoid machine or component fatigue failures. The advanced controllers can effectively handle this kind of sporadic disturbance to protect the desalination plants from damages or successfully eliminate erratic sounds coming out.
In this paper, the adaptive STW SMC is proposed to minimize the effect of this water hammer. The initial permeate pressure is 22.97 bar. The simulated condition is chosen as the worst case with the feed water concentration Cf = 50,000 mg/L. As illustrated in Figure 7, the peak pressure in the first wave is calculated at 82.57 bar without the controller. By realizing the proposed controller, this peak pressure is reduced to 24.3 bar, in which the disturbance has been significantly rejected by active controller, or eliminated to approximately 97%. This result demonstrates superior performance compared to 33.33%, which is the best result reported in the literature, where water hammering is attenuated by a PRV; or 58% disturbance attenuation as reported by Phuc et al. (2017). It is also noted that the controller prevents air being sucked into the RO system while avoiding every possible vacuum. As a result, the membrane is safely protected, and no vapor cavitation and column separation occur. The ability to deal with this high disturbance proves that the proposed controller can cope with slight perturbation such as fouling or scaling. Although the proposed controller has the ability to diminish the effect of water hammer, a surge tank should always be installed at a suitable location in RO systems for ultimate safety (Liu & Simpson 2018; Shi et al. 2019; Zhang et al. 2019).
Referring to Figure 7, even though the transient pressure has the sinusoidal wave profile with positive and negative values, the RO system is regulated effectively by the active controller. As illustrated in Figure 8, the control signals produced by the proposed controller provide the reduced chattering without saturation. In fact, the chattering is the main obstacle for its implementation. This control action is done automatically by changing the valves opening in a manner appropriating for the values of the transient pressure. The greatly reduced pressures are realized when surge protection is provided by an active control system. Transient behaviors have been immediately reduced until they coincided with steady-state values.
As mentioned earlier, the control parameters are adaptable so that the controller is more flexible to deal with uncertainties and disturbances. Figure 9 shows that the control gains are updating their values during water hammering. It can be observed that their values increase rapidly when the water hammer has just occurred. When the pressure is dramatically dampened, these values will be fixed or constant. An adaption scheme is developed to update the controller gains providing superior responses, eventually resulting in minimization of chattering effects under hydraulic shock waves. This study presents the basic concepts associated with transient flow, discusses the hydraulic surges, and introduces features of RO system design that should be considered during transient behaviors. The active controllers can eliminate hydraulic transients to protect the desalination plants while enhancing overall performance and reliability of system components in the water treatment applications.
CONCLUSIONS
Hydraulic shocks can have detrimental effects on pipes, pumps, membranes, and valves in the water desalination plants. This paper deals with comprehensive dynamical analysis and control synthesis for hydraulic transients of desalination plants. The robust monitoring of the water desalination plant is challenged by uncertainties and disturbances such as fouling, scaling, and environment changes. First, the dynamic analysis was conducted to investigate the system parameters that affect the water hammer effects. Next, a robust control system was proposed based on SMC scheme with adaptive law. The STW algorithm was implemented to alleviate chattering caused by discontinuous control structure, and to ensure robust performances of the controller in transient states. After realizing adaptive STW SMC algorithm for desalination plants, numerical simulations were carried out in order to validate the effectiveness of the proposed control scheme. The proposed control strategy provides the following novel features to the desalination system:
The robust control synthesis exploits the adaptive STW SMC strategy to provide the finite-time convergence of the sliding surface while maintaining robust performance and stability against uncertainties and disturbances.
The proposed controller ensures the tracking of product water flow and pressure set-points in a short time.
Transient pressure waves can be suppressed to prevent damages or fractures caused by excess pressure or vapor cavitation.
The presented approach provides a chattering-free SMC synthesis for the hydraulic shocks.
The presented active controllers can eliminate hydraulic shocks to protect the desalination system while enhancing overall performance and reliability. Finally, the validity and the superiority of the water production performances of the adaptive STW SMC scheme over other control algorithms have been clearly demonstrated through numerical simulations as safe and reliable to use in a desalination process.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.