Open access peer-reviewed chapter

Adsorption Cycle and Its Hybrid with Multi-Effect Desalination

Written By

Muhammad Wakil Shahzad, Kyaw Thu, Ang Li, Azhar Bin Ismail and Kim Choon Ng

Submitted: 30 October 2014 Reviewed: 12 February 2015 Published: 28 October 2015

DOI: 10.5772/60400

From the Edited Volume

Desalination Updates

Edited by Robert Y. Ning

Chapter metrics overview

2,488 Chapter Downloads

View Full Metrics


Adsorption (AD) cycle is recently pioneered for cooling and desalination applications. For water treatment, the cycle can be used to treat highly concentrated feed water, ranging from seawater, ground water, and chemically laden waste water. This chapter presents a review of the recent development of AD cycle and its hybridization with known conventional cycles such as the MED and MSF. We begin by looking at the basic sorption theory for different adsorbent–adsorbate pairs, namely the silica gel–water and the zeolite–water pairs. Under the IUPAC categorization, there are six types of isotherm behavior that capture almost all types of adsorbent–adsorbate behaviors and many isotherm correlations have been developed to described their uptake patterns, namely the Henry, Langmuir, Toth, etc. We have recently developed a correlation that can universally capture all six types of isotherms of IUPAC and it requires only four regression coefficients.


  • Thermal Desalination
  • Adsorption Desalination (AD)
  • Exergy Analysis
  • Multi-Effect Distillation (MED)
  • Hybrid Cycles

1. Introduction

Fresh water is a precious entity which is needed for economic development of every sector of a country such as agriculture, industrial, and domestic sectors. The persistent quest of economic development and exponential population growth in many countries has intensified the water demand which is projected to grow annually at a rate of 3–4% [13]. Even though 70% of the earth is covered by water, most of the part is in nonpotable state due to its salinity, either in the form of brackish, waste, and seawater [37]. About 20% of the world population is living below the acute water poverty level of 500 m3 per capita per year [8], risking their health due to poor water quality and substandard sanitation. The demand for fresh water in 2030 is expected to increase up to 6,900 billion cubic meters (Bm3), as compared to total available water resources of 4,500 Bm3 [9]. The fresh water supply and demand scenario is shown in Figure 1, and it is predicted that the sustainable natural water cycle of our earth cannot meet the projected future water demand.

Figure 1.

Fresh water supply demand gap: current and future estimates

In 2012, installed desalinated water capacity was 72 Mm3/day (mcmd) and it is estimated to increase to 98 mcmd by 2015 as reported in the literature [10]. Almost half of total desalination capacity is installed in the Gulf Cooperation Council (GCC) countries. In the 1960s and 1970s, many GCC countries relied heavily on the ground water supply for all sectors such as agriculture, industrial, and domestic sectors. However, the ground water resource in these countries has depleted rapidly over the decades due to excessive water extraction and insufficient aquifer recharge. Despite a higher desalination market share in GCC, the fresh water availability is dropping rapidly to below the acute water poverty level of 500 m3 per capital per year, caused by an exponential population growth and higher GDP thrust.

Figure 2.

Fresh water availability per capita in GCC countries

GCC countries population, fresh water availability, and desalination scenario are shown in Figure 2, spanning from the early decades in 1950 to the future years up to 2025 [1117]. Seawater desalination was started in GCC countries in the 1960s, but significantly desalinated water supply was observed in the past two decades, contributing to the overall water consumption share. Increasing trend of desalination capacities has resulted in higher energy consumption, more than 25% of total energy production in GCC countries. It is estimated that, with contracted desalination facilities, the total primary energy requirement will grow up to 255.45 GWh in 2025 with 55% increment as compared to 112.47 GWh in 2000 and CO2 emission will be doubled in 2025 as compared to 60 ktonne/year in 2000 as shown in Table 1. The income loss attributed to desalination was 3.6 billion USD in 2000 and this loss is expected to rise to 31 billion USD in 2025 as calculated by considering yearly fuel price in terms of US$/bbl as shown in Table 1.

Year Specific energy utilization for desalination (GWh_pe) CO2 production on the basis of fuel (0.527 tonne/MWh )
Displaced income due to desalination
(1 bbl = 1628kWh),
(million USD)
1970 6.32 3.3 517.2
1975 25.19 13.3 2066.0
1980 62.29 32.8 5107.6
1985 92.79 48.9 5103.7
1990 98.26 51.8 3930.6
1995 109.92 57.9 2748.0
2000 112.47 59.3 3598.9
2005 131.93 69.5 7256.2
2010 163.96 86.4 13116.6
2015 191.21 100.8 21989.5
2020 221.97 117.0 26636.7
2025 255.45 134.6 31164.5

Table 1.

Specific energy consumed and CO2 footprint in GCC countries due to desalination processes

Water production by desalination processes has significant effect on the energy requirement and environment. The intricate nexus between water, energy, and environment has encouraged scientists and engineers to innovate desalination methods with better energy efficiency and environment-friendly processes. Although RO processes are energy-efficient and dominantly used, they have certain limitations with respect to local conditions. For example, the frequent maintenance issues from high operating pressure; water quality problems in term of residuals of boron, chlorides, and bromides; and severe fluctuations in the seawater intake quality are some of the challenges faced by the RO membranes. In the GCC region, frequent occurrence of harmful algae blooms (HABs) that may contain neuroparalytic and diarrhetic toxins which can pass through the pores of membranes lead to health problems. Large fluctuations in the feed water quality have direct implication to the operation and maintenance costs of RO plants [1820]. Owing to the uncertainty of RO operation, thermal desalination methods are deemed as the dominant processes employed in desalination market in the Middle East countries, more than 65% of installed capacities.

An innovative solution to strengthen the thermally driven and yet robust multi-effect desalination (MED) is its integration with the heat-driven adsorption desalination (AD) cycle. AD cycle is an emerging low-cost desalination system that requires only low temperature waste heat or solar energy to operate the cycle. The hybridization of both cycles, called MEDAD, extend the range of downstream (last stage) temperature of conventional MED system typically from 40oC to 5oC. The additional number of stages enhances the water production of the MED cycle by 2 - 3 fold at the same top-brine temperatures. In addition, AD integration cycle has the following advantages: (i) it increases inter-stage temperature differential of each MED stages due to the lowering of the bottom-brine temperature; (ii) it helps to scavenge the ambient energy in last part of the MED stages where the latent energy is further recycled; (iii) the AD cycle utilizes only low temperature waste heat; (iv) it has almost no major moving parts; (v) it reduces the chances of corrosion and fouling due to high concentration exposed to low temperature (5oC) in the last stages; (vi) it produces additional cooling effect from last stages of MED operating below ambient temperature; and (vii) significant increase in system performance. The basics of adsorption phenomenon, AD cycle, and its hybrids and desalination processes economics are presented in detail in this chapter.


2. Basic adsorption phenomena

The understanding of an AD cycle is predicated on the availability of basic equilibria–vapor uptake or isotherms of an adsorbent–adsorbate pair at assorted pressures and temperatures. From literature, all isotherms of adsorbent–adsorbate pairs can be categorized into six types (IUPAC) and they are described by many types of empirical and semi-empirical models. The simplest adsorption isotherm model is the classical Langmuir model [21] where it assumes a homogeneous surface with a monolayer vapor uptake where all adsorbent surfaces contain a uniform charged energy. Each pore vacant site is assumed to be filled by a single vapor molecule, forming a single sorption event. Invoking the rate of gas molecules filling the adsorption sites (dθdt), as given by Ward and co-workers [2225], the expression of the Langmuir isotherm model can be derived as follows:


where μ is the chemical potential in kJ/mol and subscripts g and a are the gaseous and adsorbed phases, T is the absolute temperature, and K′ is the dimensionless constant. The isotherm, θ, can be obtained by integrating the rate equation over the energy level from Ec to infinity and in this simple case, when 0 then θ 0, and as P is large, θ approaches 1.


A real solid surface of an adsorbent contains geometrical roughness that is formed during its formation or activation process and the energetic heterogeneity has to be accounted for. The gas molecules experienced varying potential at adsorption sites of uneven energy levels and the surfaces are subdivided into infinitesimal pieces. Thus, the total adsorption of a heterogeneous surface is the sum of the product of the sorption event or surface coverage, θ˜(εi) and its probability energy distribution, χ(εi), i.e.,


where the Langmuir model can be applied with a probability function over the entire surface being summed to unity, i.e.,


Applying a condensation approximation at moderate temperature, the local surface coverage can be simplified to a step-like dirac-delta function as shown in Figure 3, and expressed below;


where, εc relates to the adsorption energy in equilibrium conditions. Hence, Equation 3 can be simplified to


Figure 3.

Step-like profile yielded from Equation 5 at moderate temperatures

Figure 4.

An illustration of symmetrical or asymmetrical Gaussian function to represent the adsorption site energy distribution for classical isotherm models

A solution of χ(ε) to represent the site energy allocation for a real adsorbent surface is obtained by approximating it to a continuous probability distribution function. Depending on the adsorbent surface energy characteristics during adsorption interactions, symmetrical or asymmetrical Gaussian functions are usually assumed, as illustrated in Figure 4.

The details of EDFs, χ(ε), for the classic Langmuir–Freundlich [26], Dubinin–Astakhov [27], Dubinin–Raduskevich [28], and Tóth isotherm models are outlined in Table 2, and integrating EDFs from the cut-off energy εc to , the exact correlation of these isotherm models can be obtained.

Model Name Adsorption Site Energy Distribution χ(ε) Isotherm Equation
Langmuir–Freundlich χ(ε)=1cexp(εε0c)[1+exp(εε0c)]2 θ=[KPexp(ε0RT)]RTc1+[KPexp(ε0RT)]RTc
Dubinin–Astakhov χ(ε)=r(εε1)r1Erexp[(εε1E)r] θ=exp[(RTElnP0P)r]
Dubinin–Radushkevich χ(ε)=2(εε1)E2exp[(εε2E)2] θ=exp[(RTElnP0P)2]
Tóth χ(ε)=1RT[exp(εε3RT)]t{1+[exp(εε3RT)]t}t+1t θ=KPexp(ε3RT){1+[KPexp(ε3RT)]t}1t

Table 2.

Assorted forms of adsorption site energy distribution functions (EDF) for the Langmuir-Freundlich, Dubinin-Astakhov (DA), Dubinin-Raduskevich (DR) and Tóth isotherm models

2.1. Universal site-Energy probability Distribution Function (EDF)

In the recent development of adsorption isotherm theory, Li [29] proposed a universal model that was able to fit all types of isotherms, as specified in Equation 7. Type I to V patterns at various temperatures could be directly captured by the equation with four regression parameters.




where the alphabet ϕ, C are constants, t is surface heterogeneity factor, and Ec denotes the characteristic energy of the adsorbent–adsorbate pair. These four parameters are required to calculate in the regression process. The rest of the letters have their usual means.

Using the unified adsorption isotherm framework, a universal adsorption site energy distribution function (EDF) was proposed, which relates directly to their isotherm types, and the proposed EDF fitted well with the statistical rate theory of adsorption. The EDF yielded a single peak asymmetrical distribution for Type I which was similar to that for the classical LF, DA, DR, and Tóth isotherm models. The EDF is given as below;


where the variable β* is a function of the adsorption site energy ε and it is expressed as


From the data available from literature and the proposed Eqs. 9 and 10, the assorted isotherms as categorized by the IUPAC can be successfully captured succinctly by these equations using only four coefficients of regression, and these energy distribution functions and isotherms are depicted in Table 3. For each type of isotherm, the corresponding energy distribution functions (EDFs) have been developed.

The adsorbent–adsorbate interaction is the key in the design of adsorption cycle, which can be functionalized to adopt to warmer ambient temperatures, particularly for operation in the summer period of semi-arid or desert regions. Figures 5(a) and (b) depict two major types of useful adsorbents in water uptake: The former is silica gel Type 3A, suitable for an AD cycle operating below 33οC ambient such as tropical weather conditions. The latter figure depicts the isotherms of zeolite (Z0-alumina phosphate oxide, AlXPhYO.nH2O). It has properties that can be tailored for adsorption/desorption at ambient temperatures up to 50οC (corresponding to desorption at 3.8–4.5 kPa (60οC isotherm) and adsorption at 2–2.8 kPa (40οC isotherm). The thermos-physical properties of both adsorbents are tabulated in Table 4, showing the BET surface-pore areas of 600–800 m2/g.

Types of Isotherms (IUPAC Categorization) Adsorbate–Adsorbent Pair / References
Type I Water–silica gel Type RD
Type I Water–silica gel Type A
Qiu J. Characterization of silica gel–water vapor adsorption, MEng Thesis, 2004, NUS
Type II Water–boehmite
Type II Δ Water–polyvinyl pyrrolidone
1. Wang, S-L, Johnston CT, David L, White JL, Stanley LH. Watervapor adsorption and surface area measurement of poorly crystalline boehmite. J Colloid Interface Sci 2003;260(1):26–35.
2. Zhang, J, Zografi G. The relationship between “BET” and “free volume”‐ derived parameters for water vapor absorption into amorphous solids. J Pharm Sci 2000;89(8):1063–72.
Type III Water-sctivated carbon
Type III Δ Water–carbon S-W nanotube
Kim P, Agnihotri S. Application of water-activated carbon isotherm models to water adsorption isotherms of single-walled carbon nanotubes. J Colloid Interface Scie 2008;325(1):64–73.
Type IV Water–boehmite
Type IV Δ Water–polyvinyl pyrrolidone
Bansal RC, Dhami TL. Surface characteristics and surface behaviour
of polymer carbons—II: adsorption of water vapor. Carbon 1978;16(5):389–95.
Type II Water–boehmite
Type II Δ Water–polyvinyl pyrrolidone
Lagorsse S, Campo MC, Magalhaes FD, Mendes A. Water adsorption on carbon molecular sieve membranes: experimental data and isotherm model. Carbon 2005;43(13):2769–79.

Table 3.

A summary of the energy distribution functions and isotherms as categorized by the IUPAC

Figure 5.

(a) shows the isotherms of water-silica gel Type 3A at four temperatures: 30oC to 45o C at increasing pressure up to 10 kPa, (b) depicts the isotherms of Zeolite (Z01- Alumina Phosphate Oxide) from 30o C to 60o C.

Properties Silica Gel Type 3A Zeolite FAM Z01
BET surface area [m2/g] 680 147.3
Porous volume [ml/g] 0.47 0.071
Apparent density [kg/m3] 770 600–700
Thermal conductivity [W/m.K] 0.174 0.113 (30oC)
Heat of adsorption (H2O) [kJ/kg of H2O] 2800 3110 (25oC)
Specific heat capacity [kJ/kg.K] 0.921 (30oC)

Table 4.

Thermo-physical properties of the silica gel type 3A and Zeolite (Z01)


3. Design of AD batch-operated cycle

There are five main components of AD system namely: (i) evaporator, (ii) adsorption and desorption reactor beds, (iii) condenser, (iv) pumps, and (v) pretreatment facility. The detailed process diagram is shown in Figure 6. For the batch-operated AD cycle, it involves two main processes.

3.1. Adsorption-assisted-evaporation

In which the vapors generated in AD evaporator are adsorbed on the pore surface area of adsorbent. The heat source is circulated through the tubes of evaporator and seawater is sprayed on the tube bundle. It is observed that the evaporation is initiated by heat source, but during adsorption process the high affinity of water vapor of adsorbent drops the evaporator pressure and contribute in evaporation. The AD evaporator operation temperature can be controlled by heat source temperature that is normally circulated in terms of chilled water. The AD evaporator can operate at a wide range of chilled water temperature varying from 5οC to 30οC to produce the cooling effect as well at low temperature operation. The vapor adsorption process continues until the adsorbent bed reaches a saturation state.

Figure 6.

Detailed schematics of an adsorption desalination cycle. The circle (filled) dots are the vapor valves and the pair of triangles refers to the liquid valves. The range of vapor pressures in evaporator and condenser are 1-2 kPa and 5 to 7 kPa

3.2. Desorption-activated-condensation

In which saturated adsorbent is regenerated using the low-grade industrial waste heat or renewable energy (desorption temperature varies from 55°C to 85°C) and desorbed vapors are condensed in a water-cooled condenser and collected as a distillate water.

It can be seen that two useful effects produced by AD cycle are the cooling effect by the first process “adsorption-assisted-evaporation” and fresh water production by converting the seawater by second process “desorption-activated-condensation”. Useful effects which are cooling and water production can be produced simultaneously by introducing the multi-bed technique.

In multi-bed AD system, the operation and switching technique is used. During operation, one or pair of adsorbent reactor beds undergo the adsorption process and at the same time one or pair of adsorbent reactor beds execute the desorption process. The time for adsorbent reactor beds operation, either adsorption or desorption, depends on the heat source temperature and silica gel quantity packed in a bed. Prior to changing the reactor duties, there is a short time interval called switching in which the adsorber bed is preheated whilst the desorber bed is precooled to enhance the performance of cycle. In AD cycles, the operation (adsorption and desorption) and switching processes are controlled by automated control scheme that can open and close the respective bed hot/cold water valves. During switching operation, all vapor valves are closed so that there is no adsorption or desorption taking place.


4. MED-AD hybrid cycle

MEDAD is a hybrid of two thermal systems namely multi-effect desalination system and adsorption cycle. The main components of this novel thermal hybrid system are (1) multi-effect distillation (MED) system, (2) adsorption desalination (AD) cycle, and (3) auxiliary equipments. In this hybridization system, the last stage of the MED is connected to adsorption beds of AD cycle for the direct vapor communication to adsorption beds. Figure 7 shows detailed flow schematic of MED plant combined with AD system.

Adsorption-based desalination is investigated by many researchers [3045] and reported that optimal specific daily water production (SDWP) for four bed scheme is about 4.7 kg/kg silica gel. The first adsorption desalination plant was installed in the National University of Singapore (NUS) which consists of four silica gel beds. Ng et al. investigated the processes using chilled water at assorted temperature and demonstrated that the specific water production of the system [46, 47]. They also introduced and patented a novel hybrid desalination method “MEDAD cycle” that is a combination of conventional MED and AD cycle [48, 49]. This novel desalination cycle can mitigate the limitations of conventional MED system to increase the system performance. This combination allows the MED last stage to operate below ambient temperature typically at 5οC as compared to traditional MED at 40οC. This not only reduces the corrosion chances but also increases the distillate production to almost 2 – 3 fold as compared to traditional MED systems.

Figure 7.

Detailed P& ID diagram for an experimental MED facility of nominal heating capacity of 12 kW. The first stage comprises a water-fired steam generator whilst the last stage (3rd stage) is open to a water-cooled cooling tower

4.1. MED-AD simulation

MEDAD cycle simulation have been conducted [50–53] and presented in Table 5. The simulation is based on a fully transient model and the predicted results are compared with conventional MED system. It is observed that at same input parameters such as a top-brine-temperature (TBT), water production can achieve up to two fold increase when the hybridized MEDAD is compared with the MED cycle alone.

Equation Remarks
Modeling equations for steam generator
[(Mhw.Cphw)]dThwdt=(mhw,hf,Thw,in)(mhwhf,Thw,out)hin,o.Ain,i(ThwTtube,i) Energy balance for the hot water flowing inside the tubes of brine heater.
[(MHX,i.Cphx,i)]dTtube,idt=hin,i.Ain,i(Thw,iTtube,i)hout,i.Aout,i(Ttube,iTv,i) Energy balance for metal tubes.
dMb,idt=mf,imb,imv,i Mass balance for the seawater inventory in the evaporator side of the brine heater.
Qin=ho,iAi(TtTv,i) [(Mb,i.Cpb,Tb)+(MHX,i.CpHX)]dTidt=(mf,ihf,Tf)(mb,ihf,Tb)(mv,ihg,Tv)+Qin,i Energy balance for the evaporator side of the steam generator.
Mb,idXb,idt=(mf,iXf,i)+(mb,iXb,i)(mv,iXv,i) Material/concentration balance
Nu=hin,idin,iKtube,i=0.023Rel0.80Prl0.40 Convective heat transfer coefficient equation
Rwall,i=ln(dout,idin,i)2.π.Ktube,i.Ltube,i Tube wall resistance
ho(ν2k3g)1/3=0.0007Re0.2Pr0.65q0.4 Falling film evaporation heat transfer coefficient, Han and Fletcher’s correlation
Overall heat transfer coefficient
Modeling equations for intermediate stages
[(Ml,i+1.Cpl,Tcond)]dTcond,i+1dt=[mv,hfg,Tv]i[hin.Ain(TcondTtube)]i+1 Energy balance for the condenser side of the i th effect
[(MHX,i+1.Cphx,i+1)]dTtube,i+1dt=hin,i+1.Ain,i+1(Tcond,i+1Ttube,i+1)hout,i+1.Aout,i+1(Ttube,i+1Tv,i+1) Energy balance for tube metal
dMb,i+1dt=mf,i+1mb,i+1mv,i+1 Brine inventory balance
Qin,i+1=hout,i+1Ai+1(Tt,i+1Tv,i+1) [(Mb,i+1.Cpb)+(MHX,i+1.CpHX,i+1)]dTi+1dt=(mf,i+!hf,Tf)(mb,i+1hf,Tb)(mv,i+1hg,Tv)+Qin,i+1 Energy balance for evaporator side
Mb,i+1dXb,i+1dt=(mf,i+1Xf,i+!)(mb,i+1Xb,i+1)(mv,i+!Xv,i+!) Material/ concentration balance
Nu=hin,i+1Li+1Ktube,i+1=0.728[ghfg,Tcondρl,Tcond(ρlρv)TcondKl,Tcond3μl,Tconddi(Tv,i+1Ttube,i+1)]1/4 Nusselt film condensation correlation for the calculation of the heat transfer coefficient inside the condenser tubes
UiAi=11hin,iAin,i+Rwall,i+1hout,iAout,i Overall heat transfer coefficient equation
MED last stage connected with AD beds
[(Mb,n.Cpb)+(MHX,n.CpHX,n)]dTndt=(mf,nhf,Tf)(mb,nhf,Tb)(Msghg,Tv)dqadsdt+Qin,n Qin,n=hout,nAn(Tt,nTv,n) Energy balance for the condenser side of the i th effect

Table 5.

MEDAD modeling equations

Figure 8 shows the transient temperature profiles of a MEDAD cycle. It can be seen that last stages of MED are operating below ambient temperature due to hybridization. It can also be noticed that MED last stages profiles are affected by cyclic AD operation.

Figure 8.

MEDAD cycle components temporal temperature profiles

Figure 9 mapped the performance parameters of hybrid MEDAD cycle (concentration, GOR, PR, and WPR). It can be observed that the batched-mode water vapor uptake by the coupled AD cycle dominates the performance of the cycle. The performance of the MEDAD cycle with different additional effects with the conventional eight effect MED cycle as baseline cycle is studied in terms of water production rate.

Figure 9.

MEDAD hybrid cycle performance (concentration, GOR, PR and WPR)

Figure 10 shows the quantum jump in the water production rate of the proposed MEDAD cycle. Another aspect of this hybridization is that the desorbed water vapor from the AD cycle can be recycled back into the MED system for further energy recovery.

Figure 10.

Percentage improvement in water production rate by the MEDAD cycle

4.2. MED-AD experimentation

A three-stage MED system is designed [54], fabricated, and installed in NUS as shown in Figure 11. In MED stages, vapor emanation from feed seawater is achieved by falling film-evaporation process. Evaporation energy is recovered by series of reutilization of vapor condensing energy in successive stages of those produced in preceding stages. Process of vapor production and energy recovery by condensation continues until the last stage of MED. The vapors from the last stage are then directed toward AD beds where they are adsorbed on the adsorbent surface. Adsorbent high affinity for water vapor drops the pressure and hence the saturation temperature of last stages falls below ambient, typically up to 5οC. It is observed that this drop in pressure and temperature of last stage also affected the operational parameters of the few preceding stages.

Figure 11.

MEDAD system installed in NUS

Experiments are conducted in two steps. In first part, system is operated as a conventional MED at assorted heat source temperature ranges from 38οC to 70οC. In second part, experiments are conducted as a hybrid MEDAD system at assorted heat source temperature ranges from 15οC to 70οC and results are compared with conventional MED system.

Figure 12.

MEDAD components temperature profiles at 38oC heat source temperature

Figure 12 shows the instantaneous temperatures of MED and MEDAD components at a heat source temperature of 38οC. It can be seen that steady-state events (minimum temperature fluctuations) occur after 1 hour from start-up and experiments for distillate collection are continued for 4 to 5 hours. It is noticed that the inter-stage temperature difference (ΔT) is more than twice per stage as compared to the conventional MED stages. This is attributed to the vapor uptake by the adsorbent of AD cycle, resulting in the increase of vapor production. The MEDAD cycle yields a stage ΔT from 3οC to 4οC as compared to 1οC or less in the case of MED alone. Table 6 shows the comparison of MEDAD and MED components steady-state temperature values.

Figure 13 shows the distillate production trace at heat source 38οC from MED stages, AD condenser and combined. The batch operated AD production can be seen clearly. At the start of desorption, the production is higher and it drops with time to zero during the switching period while MED stages production is quite stable. Small fluctuations in MED water production may be due to the fluctuations in the spray of the feed that affect the condensation rate. It can be seen that hybridization boost water production 2 – 3 fold as compared to conventional MED system at same TBT. Water production profiles are similar as explained in simulated results.

Heat source temperature
SG SG Stage-2 Stage-2 Stage-3 Stage-3
70 63.9 54.4 62.8 50.8 62.3 46.7 56.6 26.1
65 59.6 50.7 58.7 47.0 58.1 42.5 53.9 25.8
60 55.4 49.2 54.4 45.7 53.8 41.3 49.4 25.5
55 51.5 48.5 50.5 44.8 49.9 40.8 45.8 25.1
50 47.4 44.6 46.7 41.6 46.2 37.9 42.7 23.2
45 43.9 41.1 43.4 38.6 43.0 34.9 40.3 22.4
40 38.7 36.6 38.4 34.8 38.1 31.6 36.1 21.0
38 37.4 35.4 37.2 33.5 37.0 30.3 35.3 18.2
Operating limit of Conventional MED. The lower operational points are from MEDAD Hybrids
35 - 30.9 - 27.9 - 23.5 - 16.1
30 - 26.0 - 23.7 - 19.1 - 12.2
25 - 22.1 - 19.9 - 16.3 - 10.1
20 - 18.1 - 14.2 - 11.3 - 7.1
15 - 13.2 - 11.5 - 9.4 - 5.6

Table 6.

A comparison of conventional MED and hybrid MEDAD systems components temperatures at different heat source temperatures

Figure 13.

Conventional MED and hybrid MEDAD cycle water production profiles at 38oC heat source temperature

Figure 14.

MED and MEDAD steady state water production at different heat source temperatures

Figure 14 shows the comparison of water production of MED and hybrid MEDAD cycles at assorted heat source temperatures. Quantum increase in water production (two- to threefold) can be observed at all heat source temperatures. These results have good agreement with simulation results. It can also be seen that in conventional MED system last stage temperature is limited to 38οC due to condenser operating with cooling water from cooling tower. While in the case of hybrid MEDAD, the last stage temperature can be as low as 5οC because there is no condenser and last stage is connected to AD beds for vapor adsorption. This higher overall operational gap in proposed hybrid MEDAD cycle helps to insert more number of stages (up to 19 stages) as compared to conventional MED system (about 4–6 stages). More number of stages increases the vapor condensation heat recoveries and hence the water production at same top brine temperatures.


5. Exergy analysis for operational cost apportionment

A computation model is developed for the cogeneration plant where the properties of expanding steam, such as enthalpy (h) and entropy (s) at a given temperature and pressure, are computed for the key states in the schematic diagrams of “PP”, “PP+MED,” and “PP+MEDAD” cycles. The approach employed here is to calculate the total exergy changes or destruction across the inlet and exit sections of the equipment. For example, the exergy associated with the turbines is the sum of all contributions from (a) the HP-T unit at the same mass flow rate across it, (b) the LP-T unit until the extraction point, and (c) the exergy changes after the point of steam extraction, i.e., ET,12=E1a+Eb1'+E1'2. On the basis of total exergy available across whole system, proportion utilized by power plant and desalination system is calculated.

For this comparison study, an assorted range of bled-steam is extracted at low pressure but the mass flow rates are varied to cover the expected practical operational ranges, typically from 10% to 50% of the total flow. At 20% bled-steam from the LP-T, the ratios of power-to-water for exergy and energetic analyses are found to be 95.7%:4.3% (exergy) and 72.2%:27.8% (energetic) methods, respectively. This implies that if the 20% steam were to be continuously consumed by the low pressure turbines (LP-T), its work contribution from LP_T would be insignificant, i.e., a maximum of 4.3%. However, the energetic value of bled-steam accounts for a disproportionate share of 27.8% due primarily to the high latent heat content of the low pressure steam. This ratio of energy-to-exergy shares of the total working steam is found to be 4 – 7 fold higher, descending with the larger amount of bled-steam as shown in Figure 15. The higher effectiveness of working steam incurred at the MED cycle is attributed to the better thermodynamic matching of steam’s latent energy.

Figure 15.

Exergetic and energetic proportions at different percent of steam extraction

Figure 16.

A comparison of life-cycle unit water cost for various desalination methods

On the basis of the above analysis for primary fuel cost and with data from the published literature [55, 56], the life-cycle cost (LCC) of water production is compared for all capital expenditure (Capex) and operation expenditure (Opex), across all proven industrial processes, as shown in Figure 16. Exergy factor calculated above is utilized only for thermal and electricity cost calculations. Energy based analysis has good agreement with GWI [56] data. It can be seen that by LCC, unit water production cost is highest for PP+MSF which amounts to US$ 1.201/m3, whilst the lowest unit cost is the PP+MEDAD method which is only at US$ 0.485/m3 and this unit cost is even lower than the LCC of reverse osmosis (RO) plants.


6. Summary

Recent developments in adsorption theory, adsorption desalination (AD), and conventional MED desalination cycles have been reviewed in this chapter. We highlight the key role of AD cycles which can be hybridized with the proven cycles such as the MED cycle, exploiting the thermodynamic synergy between the thermally driven cycles that significantly improve the water production yields. Experiments were conducted in a lab-scale pilot MEDAD and confirmed the excellent synergetic effects that boosted the water production up to two- to threefold over the conventional MED. We believe that if the hybrid MEDAD cycles are well optimized and operated, it can achieve high GOR and the projected LCC of water production can be lowered to as low as US$ 0.485/m3.


7. Abbreviation

MED; Multi-effect desalination

AD; Adsorption desalination

RO; Reverse osmosis

SDWP; Specific daily water production

EDF; Energy distribution function

GOR ; Gain output ratio

SG; Steam generator

PR; Performance ratio

WPR; Water production ratio

TBT; Top brine temperature

PP; Power plant

HP-T; High pressure turbine

LP-T; Low pressure turbine

LBT; Lower brine temperature

LCC; Life cycle costing

PDF; Probability distribution function

P Pressure [kPa]
P0 Saturation pressure [kPa]
R Universal gas constant [kJ/kmol.K]
T Temperature [K]
t Surface heterogeneity for Tóth [-]
Greek Symbols
γ Dubinin–Astakhow power factor [-]
ε Adsorption site energy [kJ/kmol, or kJ/kg]
εc Equilibrium adsorption site energy [kJ/kg]
ε1 to ε3 Reference energy [kJ/kmol]
θ Adsorption uptake [kg/kg of adsorbent]
θ˜ Local adsorption uptake [kg/kg of adsorbent]
μ Chemical potential [kJ/kmol]
χ(ε) Energy distribution function [-]
α Adsorbed phase
g Gaseous phase



The authors wish to thank National Research Foundation (NRF), Singapore (grant WBS no. R-265-000-399-281), and King Abdullah University of Science & Technology (KAUST) (Project no. 7000000411) for financial support for MED plant at the National University of Singapore.


  1. 1. Sources of fresh water, UNESCO module-7 (
  2. 2.
  3. 3. Earth’s water supply (
  4. 4. National oceanic and atmospheric administration, United State department of commerce. (
  5. 5. The encyclopedia of earth (
  6. 6.
  7. 7. Science education resource center, Carleton College, (
  8. 8. Hoekstra AY, Mekonnen MM. Global water scarcity: the monthly blue water footprint compared to blue water availability for the world’s major river basins, Value of water research report series no. 53, UNESCO-IHE Institute for Water Education.
  9. 9. Thu K, Chakraborty A, Kim Y-D, Myat A, Saha BB, Ng KC. Numerical simulation and performance investigation of an advanced adsorption desalination cycle. Desalination 2013;308:209–18.
  10. 10. GWI, Desalination Markets 2007, A Global Industry Forecast (CD ROM), Global Water Intelligence, Media Analytics Ltd., Oxford, UK, 2007,
  11. 11. Dawoud MA, Al Mulla MM. Environmental impacts of seawater desalination: Arabian Gulf case study. Int J Environ Sustain 2012;1(3),22–37.
  12. 12. Dawoud MA. The role of desalination in augmentation of water supply in GCC countries, Desalination 2005;186:187–98.
  13. 13. Abdulrazzak MJ. Water supplies versus demand in countries of Arabian Peninsula. J Water Res PL ASCE 1995;121(3):227–34.
  14. 14. Fath H, Sadik A, Mezher T. Present and future trend in the production and energy consumption of desalinated water in GCC countries. Int J Thermal Environ Eng 2013;5(2):155–65.
  15. 15. Nezhad H. A report on world energy scenarios to 2050: issues and options, Decision Sciences Metropolitan State University Minneapolis, MN, Sep-2009.
  16. 16. Al-Zubari WK. The Water-Energy Nexus in the GCC Countries Evolution and Related Policies, Sixth "Zayed Seminar" on Green Economy: Success Stories from the GCC, 8–9 May 2013, UN virtual water learning center, the Arabian regional center, Arabian Gulf University, Bahrain.
  17. 17. Clean Water for a Healthy World, Secretariat of the Convention on Biological Diversity. 2009. Drinking Water, Biodiversity and Poverty Reduction: A Good Practice Guide, World water day March 22, 2010.
  18. 18. Caron DA, Garneau ME, Seubert E, Howard MDA, Darjany L, Schnetzer A, Cetinic I, Filteau G, Lauri P, Jones B, Trussell S. Harmful algae and their potential impacts on desalination operations off southern California. Water Res 2010;44:385–416.
  19. 19. Huttner KR. Overview of existing water and energy policies in the MENA region and potential policy approaches to overcome the existing barriers to desalination using renewable energies. Desalin Water Treat 2013;51:87–94.
  20. 20. Ghaffour N, Missimer TM, Amy GL. Combined desalination, water reuse, and aquifer storage and recovery to meet water supply demands in the GCC/MENA region. Desalin Water Treat 2013;51:38–43.
  21. 21. Langmuir I. The adsorption of gases on plane surfaces of glass, mica and platinum, J Am Chem Soc 1918;40:1361–403.
  22. 22. Ward C, Findlay R, Rizk M. Statistical rate theory of interfacial transport. I. Theoretical development. J Chem Phys 1982;76:5599–605.
  23. 23. Elliott JAW, Ward CA. Statistical rate theory description of beam-dosing adsorption kinetics, J Chem Phys 1997;106:5667–76.
  24. 24. Elliott JAW, Ward CA. Statistical rate theory and the material properties controlling adsorption kinetics, on well defined surfaces. Stud Surf Sci Catal 1997;104:285–333.
  25. 25. Elliott JAW, Ward CA. Temperature programmed desorption: a statistical rate theory approach, J Chem Phys 1997;106:5677–84.
  26. 26. Rudzinski W, Borowiecki T, Dominko A, Panczyk T. A new quantitative interpretation of temperature-programmed desorption spectra from heterogeneous solid surfaces, based on statistical rate theory of interfacial transport: the effects of simultaneous readsorption. Langmuir 1999;15:6386–394.
  27. 27. Rudzinski W, Lee S-L, Panczyk T, Yan C-CS. A fractal approach to adsorption on heterogeneous solids surfaces. 2. Thermodynamic analysis of experimental adsorption data. J Phys ChemB 2001;105:10857–66.
  28. 28. Rudzinski W, Lee S-L, Yan C-CS, Panczyk T. A fractal approach to adsorption on heterogeneous solid surfaces. 1. The relationship between geometric and energetic surface heterogeneities. J Phys Chem B 2001;105:10847–56.
  29. 29. Li A. Experimental and theoretical studies on the heat transfer enhancement of adsorbent coated heat exchangers, Doctor of Philosophy, Department of Mechanical Engineering, National University of Singapore, 2014.
  30. 30. Wang X, Chakarborty A, Ng KC, Saha BB. How heat and mass recovery strategies impact the performance of adsorption desalination plant: theory and experiments. Heat Transfer Eng 2007;28:147–53.
  31. 31. Saha BB, Choon NK, Chakarborty A, Thu K. Desalination System and Method, in: W.P. WO2011010205A1 (ed.), 2011.
  32. 32. Thu K. Adsorption desalination: theory & experiments, in: Doctoral Thesis, National University of Singapore, 2010.
  33. 33. Saha BB, Alam KCA, Akisawa A, Kashiwagi T, Ng KC, Chua HT. Two-stage non-regenerative silica gel-water adsorption refrigeration cycle, ASME Advanced Energy Systems Division (Publication) AES 2000;40:65–9.
  34. 34. Ng KC. Recent developments in heat-driven silica gel-water adsorption chillers. Heat Transfer Eng 2003;24:1–3.
  35. 35. Ng KC, Wang XL, Chakraborty A. in: Apparatus and Method for Desalination, 2006.
  36. 36. Ng KC, Saha BB, Chakraborty A, Koyama S. Adsorption desalination quenches global thirst. Heat Transfer Eng 2008;29:845–8.
  37. 37. Wu JW, Hu EJ, Biggs MJ. Thermodynamic cycles of adsorption desalination system. Appli Energy 2012;90:316–22.
  38. 38. Wu JW, Biggs MJ, Pendleton P, Badalyan A, Hu EJ. Experimental implementation and validation of thermodynamic cycles of adsorption-based desalination. Appl Energy 2012;98:190–7.
  39. 39. Boelman EC, Saha BB, Kashiwagi T. Experimental investigation of a silica gel-water adsorption refrigeration cycle – the influence of operating conditions on cooling output and COP. 1995;101:358–66.
  40. 40. Khan MZI, Alam KCA, Saha BB, Akisawa A, Kashiwagi T. Study on a re-heat two-stage adsorption chiller – the influence of thermal capacitance ratio, overall thermal conductance ratio and adsorbent mass on system performance. Appl Thermal Eng 2007;27:1677–85.
  41. 41. Khan MZI, Alam KCA, Saha BB, Akisawa A, Kashiwagi T. Performance evaluation of multi-stage, multi-bed adsorption chiller employing re-heat scheme. Renewable Energy 2008;33:88–98.
  42. 42. Ng KC, Thu K, Saha BB, Chakraborty A. Study on a waste heat-driven adsorption cooling cum desalination cycle. Int J Refrigeration 2012;35:685–93.
  43. 43. Thu K, Ng KC, Saha BB, Chakraborty A, Koyama S. Operational strategy of adsorption desalination systems. Int J Heat Mass Transfer 2009;52:1811–6.
  44. 44. EI-Sharkawy I, Thu K, Ng K, Saha BB, Chakraborty A, Koyama S. Performance improvement of adsorption desalination plant: experimental investigation. Int Rev Mech Eng 2007;1:25–31.
  45. 45. Ng KC, Wang X, Lim YS, Saha BB, Chakarborty A, Koyama S, Akisawa A, Kashiwagi T. Experimental study on performance improvement of a four-bed adsorption chiller by using heat and mass recovery. Int J Heat Mass Transfer 2006;49:3343–8.
  46. 46. Wang X, Ng KC. Experimental investigation of an adsorption desalination plant using low-temperature waste heat. Appl Thermal Eng 2005;25:2780–9.
  47. 47. Ng KC, Wang XL, Gao LZ, Chakraborty A, Saha BB, Koyama S. Apparatus and method for desalination, SG Patent application number 200503029-1 (2005) and WO Patent no. 121414A1 (2006).
  48. 48. Ng KC, Thu K, Amy G, Chunggaze M, Al-Ghasham TY. An Advanced ADMED Cycle For Low-Temperature Driven Desalination, US Provisional Application No. 61/450,165.
  49. 49. Ng KC, Thu K, Amy G, Chunggaze M, Ghasham TYA. A Regenerative Adsorption Distillation System, in: WO Patent WO2012121675A1, 2012.
  50. 50. Shahzad MW, Ng KC, Thu K, Saha BB, Chun WG. Multi effect desalination and adsorption desalination (MEDAD): a hybrid desalination method. Appl Thermal Eng 2014;72:289–97.
  51. 51. Thu K, Kim Y-D, Amy G, Chun WG, Ng KC. A synergetic hybridization of adsorption cycle with the multi-effect distillation (MED). Appl Thermal Eng 2014;62:245–55.
  52. 52. Ng KC, Thu K, Oh SJ, Ang L, Shahzad MW, Ismail AB. Recent developments in thermally-driven seawater desalination: Energy efficiency improvement by hybridization of the MED and AD cycles. Desalination 2015;356:255–70.
  53. 53. Ng KC, Thu K, Shahzad MW, Chun W. Progress of adsorption cycle and its hybrids with conventional multi-effect desalination processes. IDA J Desalin Water Reuse, DOI 10.1179/2051645214Y.0000000020.
  54. 54. Shahzad MW, Myat A, Chun WG, Ng KC. Bubble-assisted film evaporation correlation for saline water at sub-atmospheric pressures in horizontal-tube evaporator. Appl Thermal Eng 2013;50:670–6.
  55. 55. Lattemann S, Kennedy MD, Schippers JC, Amy G. Chapter 2; Global Desalination Situation.
  56. 56.

Written By

Muhammad Wakil Shahzad, Kyaw Thu, Ang Li, Azhar Bin Ismail and Kim Choon Ng

Submitted: 30 October 2014 Reviewed: 12 February 2015 Published: 28 October 2015