Some Theoretical Aspects of Tertiary Treatment of Water/Oil Emulsions by Adsorption and Coalescence Mechanisms: A Review

: The massive increase in the volumes of oily contaminated produced waters associated with various industrial sectors has initiated considerable technological and scientiﬁc efforts related to the development of new cleaning strategies. The petrochemical industry (oil and gas production and processing) contributes to those volumes by approximately 340 billion barrels per year. The removal of emulsiﬁed oily components is a matter of particular interest because the high emulsion stability necessitates sophisticated technological approaches as well as a deep theoretical understanding of key mechanisms of oil/water separation. This review deals with the theoretical aspects of the treatment of emulsiﬁed oil/water mixtures and is particularly focused on tertiary treatment, which means the reduction of the oil content from 70–100 ppm to below 10 ppm, depending on national regulations for water discharge. The review concerns the mechanisms of oil/water separation and it covers the (i) adsorption isotherms, (ii) kinetics of adsorption, (iii) interfacial interactions between oil/water mixtures and solid surfaces, and (iv) oil/water separation techniques based on the wettability of solid/oil/water interfaces. The advantages and drawbacks of commonly used as well as newly proposed kinetic and adsorption models are reviewed, and their applicability for the characterization of oil/water separation is discussed. The lack of suitable adsorption isotherms that can be correctly applied for a description of oil adsorption at external and internal solid surfaces of both nonporous and porous structures is pointed out. The direct using of common isotherms, which were originally developed for gas adsorption, often leads to the incorrect data description because the adsorption of oily components at solid surfaces does not ﬁt the assumptions from which these models were originally derived. Particularly, it results in problematic calculations of the thermodynamic parameters of sorption. The importance of nonlinear analysis of data is discussed, since recent studies have indicated that the error structure of experimental data is usually changed if the original nonlinear adsorption isotherms are transformed into their linearized forms. The comparison between the pseudo-ﬁrst-order and pseudo-second-order kinetic models was performed. It was shown that the correlation between data and models strongly depends on the selection of data, particularly on the frequency of collected data in time scale. The wettability of solid surfaces by oil in air and under water is discussed, regarding the surface morphology of surfaces. We demonstrate that the combination of surface chemistry and topology strongly inﬂuences the separation of oil/water emulsions.


Introduction
Oil and grease (O&G) are organic substances composed of hydrocarbons, fatty acids, soaps (surfactants), lipids and waxes. Petroleum wastewaters from various sources (crude oil, gas, shale gas extraction and oil refineries) represent the largest volumes of oily polluted waters, followed by metal processing wastewater and food processing wastewater [1]. As an illustration, the global volumes of produced water volumes were 202 billion barrels in 2014 and were estimated to be approximately 340 billion barrels in 2020 [2]. The estimated produced water volume from shale reservoirs over the first 5-10 years of production ranges from 1.7 to 14.3 million liters per well [3]. Emulsification of oil/water mixtures is a phenomenon frequently occurring in water treatment processes because emulsifiers and surfactants alter the chemical character of an oil droplet surface. In the petroleum industry, chemically enhanced oil recovery (EOR) leads to the production of a vast amount of oilemulsified wastewater because surfactants (together with polymers such as polyacrylamide and alkaline compounds) are key additives that are flooded into wells in order to enhance the recovery of bypass oil and oil trapped in porous reservoirs [4]. Both anionic and cationic surfactants are used, dependent on the character of natural reservoirs of oil and gas [5]. More details about various aspects of EOR and applied surfactants can be found in the review of Gbadamosi et al. [6]. The promotion of emulsion instabilities is an important feature that contributes to the separation process in deep-bed filtration. The most important phenomenon in the destabilization of oil-in-water emulsions is coalescence. More details about the mechanisms of coalescence, experimental techniques, and modeling can be found, for instance, in the review of Kamp et al. [7], and the theoretical aspects of emulsions' stability are described in the chapter by Danov [8]. On the other hand, conservation of a long-term stability of emulsions is a crucial point in the food and cosmetic industry. Various types of surfactants (such as small molecule surfactants, phospholipids, proteins, and polysaccharides [9]) and their mixtures have been studied recently, as summarized in review of McClements et al. [10], gelatin-based emulsifiers were summarized by Ting Zang et al. [11], and water-in-oil emulsions stabilized by surfactants, biopolymers and particles were reviewed by Zembyla et al. [12].
In the petroleum industry, chemically enhanced oil recovery leads to the production of a vast amount of oil-emulsified wastewater. The key factors that influence the efficiency of emulsion separation are the properties of the filtration medium, namely, the (i) surface free energy, (ii) roughness and pore size, (iii) size, and (iv) permeability; the properties of the dispersed and continuous phase, namely, the (i) droplet size, (ii) interfacial tension, (iii) densities of components, (iv) viscosity, and (v) presence of surfactants; and the operating conditions, namely, the (i) inflow velocity, (ii) pressure drop, (iii) emulsion concentration, (iv) temperature, and (v) filter and flow orientation [13].
Deep-bed filtration involves three fundamental mechanisms: (i) transport mechanism, (ii) attachment mechanism, and (iii) detachment mechanism. The transport mechanism includes straining, interception, inertia, sedimentation, diffusion, and hydrodynamic action [14]. The attachment mechanism is realized by London-Van der Waals forces, Born repulsive forces, hydration, mutual adsorption, hydrophobic effects, and steric interactions [15]. The detachment mechanism involves reverse flow flushing and coalescence. The transport mechanism is not discussed in this review, and much relevant information can be found in the comprehensive review of Jegatheesan and Vigneswaran [15]. In this review, the main attention is given to the attachment and detachment mechanisms, including adsorption isotherms, kinetic models, diffusion in adsorbents, wettability of surfaces, and induced coalescence on solid surfaces.

Adsorption Isotherms
The most important factors affecting the adsorption process are listed in Figure 1. Adsorbents occur in various forms such as powder, small pellets, beads or granules, with a very porous structure having fine pores and pore volumes up to 50% of total particle volume. The most important feature of an adsorbent is the porous structure resulting in high surface area; a highly porous solid may be carbonaceous or inorganic in nature, synthetic or naturally occurring, and in certain circumstances may have true molecular sieving properties. Many adsorbent materials are amorphous and contain complex networks of interconnected micropores (D < 2 nm), mesopores (2 nm < D < 50 nm) and macropores (D > 50 nm). Adsorbing molecules must first pass through the fluid film, which is external to the adsorbent particle, then through the macroporous structure into the micropores where the molecules are adsorbed. Figure 2 illustrates the different types of pores. El-Naas [16]).
Adsorbents occur in various forms such as powder, small pellets, beads or granules, with a very porous structure having fine pores and pore volumes up to 50% of total particle volume. The most important feature of an adsorbent is the porous structure resulting in high surface area; a highly porous solid may be carbonaceous or inorganic in nature, synthetic or naturally occurring, and in certain circumstances may have true molecular sieving properties. Many adsorbent materials are amorphous and contain complex networks of interconnected micropores (D < 2 nm), mesopores (2 nm < D < 50 nm) and macropores (D > 50 nm). Adsorbing molecules must first pass through the fluid film, which is external to the adsorbent particle, then through the macroporous structure into the micropores where the molecules are adsorbed. Figure 2 illustrates the different types of pores. Porosity strongly depends on the processing conditions. For instance, polyethylene pellets, which always have smooth surfaces as a consequence of the route of preparation (extrusion) ( Figure 3A) and linear density polyethylene (LDPE) powders, depending on the route of the preparation (mostly grinding; partly precipitation from a solution) can have very porous structures, as demonstrated by SEM and profilometry analysis and shown in Figure 3A,B,C. It results in an increase of the specific surface area in 32,300% at the same average dimensions.  Adsorbents occur in various forms such as powder, small pellets, beads or granules, with a very porous structure having fine pores and pore volumes up to 50% of total particle volume. The most important feature of an adsorbent is the porous structure resulting in high surface area; a highly porous solid may be carbonaceous or inorganic in nature, synthetic or naturally occurring, and in certain circumstances may have true molecular sieving properties. Many adsorbent materials are amorphous and contain complex networks of interconnected micropores (D < 2 nm), mesopores (2 nm < D < 50 nm) and macropores (D > 50 nm). Adsorbing molecules must first pass through the fluid film, which is external to the adsorbent particle, then through the macroporous structure into the micropores where the molecules are adsorbed. Figure 2 illustrates the different types of pores. Porosity strongly depends on the processing conditions. For instance, polyethylene pellets, which always have smooth surfaces as a consequence of the route of preparation (extrusion) ( Figure 3A) and linear density polyethylene (LDPE) powders, depending on the route of the preparation (mostly grinding; partly precipitation from a solution) can have very porous structures, as demonstrated by SEM and profilometry analysis and shown in Figure 3A,B,C. It results in an increase of the specific surface area in 32,300% at the same average dimensions. Porosity strongly depends on the processing conditions. For instance, polyethylene pellets, which always have smooth surfaces as a consequence of the route of preparation (extrusion) ( Figure 3A) and linear density polyethylene (LDPE) powders, depending on the route of the preparation (mostly grinding; partly precipitation from a solution) can have very porous structures, as demonstrated by SEM and profilometry analysis and shown in Figure 3A-C. It results in an increase of the specific surface area in 32,300% at the same average dimensions.
The physical models, which describe the dependence of the amount of adsorbed species on the equilibrium concentration of those species in bulk solution at constant temperature, are called adsorption isotherms. The majority of common adsorption isotherms were originally developed for gas adsorption and then adapted for the adsorption of liquids or species dissolved or dispersed in liquids. There are some features that differentiate the applicability of adsorption isotherms for liquids in comparison to gases [17], as will be discussed below.

Langmuir Isotherm (LI)
The LI is the most commonly mentioned isotherm in the literature for both gases and liquids [18][19][20][21][22]; however, adsorption from liquids does not frequently fit the assumptions used for LI derivation. These assumptions are that (a) adsorbed species can form only a monolayer and that (b) all the reaction sites on the surface of the adsorbent are equivalent. This implies that the enthalpy of adsorption is independent of the surface coverage, that (c) a molecule (individual species) can be adsorbed on at a vacant site, and that (d) the probability of adsorption is independent of the occupation of neighboring sites (no interaction between adsorbed species). The LI is expressed in the form given by Equation (1) or in its linear form.
where c e (mg/L) is the equilibrium concentration, q e (mg/g) is the equilibrium adsorption capacity of the adsorbent, q m (mg/g) is the maximum adsorption capacity, and K L = k a /k d is the Langmuir adsorption (equilibrium) constant.
The LI is often expressed in linear form (Equation (2)), which simplifies the determination of unknown parameters (q m , K L ) using linear regression. The boundaries of linearity are also the boundaries of validity of the LI.
A useful parameter associated with the LI is called the separation factor R L (Equation (3)).
R L corresponds to the adsorption processes according to the following criteria [23]: Case I. R L > 1: The adsorption is unfavorable (an increase in Gibbs free energy of adsorption).
Case II. 1 > R L > 0: The adsorption is favorable (a decrease in Gibbs free energy). Case III. R L = 1: It characterizes a linear adsorption (unoccupied sites at the adsorbent are randomly occupied by adsorbate proportionally to their concentration, and only one reaction site is occupied by one species).
Case IV. R L = 0: The desorption process is irreversible.
The key disadvantage of the LI is that it fails to describe adsorption at rough surfaces of the adsorbent. Rough inhomogeneous surfaces have multiple site types available for adsorption, and some parameters, such as the heat of adsorption, vary from site to site.

Brunauer-Emmet-Teller (BET) Isotherm
The BET isotherm, expressed by Equation (4), is the most widely employed isotherm for describing the multilayer adsorption of gases, and it has been adapted for the adsorption of/from liquids [24][25][26]. There is no limit on how many layers can be deposited on the surface.
where c e (mg/L) is an equilibrium concentration, q e (mg/g) is the equilibrium adsorption capacity of the adsorbent, q m (mg/g) is the maximum adsorption capacity, c s is the saturation concentration of the solute in water, and C is a constant defined as: ∆H m S (J/mol) is the molar enthalpy of adsorption related to the first layer, ∆H m L is the molar enthalpy of condensation, and g 0 is the entropic factor. The BET isotherm consists of two adjustable parameters (q m , C). Because Equation (4) did not give good fit of the experimental data, some authors proposed keeping c s as an adjustable parameter [27,28]. In this case, there are three adjustable parameters in the BET isotherm, and the equation gives a better fit of the experimental data; however, large discrepancies between the calculated c S and the experimental values were found. This resulted in attempts to redefine the meaning of c s in Equation (4). For instance, Miller et al. [29] defined c s as the concentration at which the adsorbent is saturated by adsorbate.
Ebadi et al. [30] proposed a different approach, and instead of using an ad hoc parameter such as c s , they directly modified the BET isotherm derivation to the form expressed by Equation (6).
K S is the equilibrium adsorption constant for the first (adjacent) layer, and K L is the equilibrium adsorption constant for the upper layers. These two constants define the constant C; C = K s /K L . This model was successfully applied for various adsorption processes, including oil adsorption [31][32][33].

Freundlich Isotherm (FI)
The FI was originally an entirely empirical function first proposed for a description of gas adsorption onto a solid surface but is also frequently used for adsorption from solutions. Unlike the LI, the FI is used for a description of multilayer adsorption on heterogeneous surfaces, and a non-uniform (exponential) distribution of adsorption heats is considered. Heterogeneous surfaces with a non-uniform distribution of adsorption heats can be differentiated from homogeneous surfaces by calorimetrical measurements of the heat of adsorption. The heat of adsorption per mass of adsorbed solute on a homogeneous surface is constant, whereas for a heterogeneous surface, it depends on the degree of surface coverage. The Freundlich isotherm can also be derived in different ways [34,35], for instance, by attributing the change in the equilibrium constant of the binding process to the heterogeneity of the surface and the variation in the heat of adsorption [14].
lnq e = lnK F + 1 n lnc e (8) where 1/n is the heterogeneity factor, n characterizes the intensity of the adsorption process and characterizes the relative distribution of the energy and the heterogeneity of the adsorbent reactive sites, and K F (L/mg) is the Freundlich adsorption constant.

Dubinin-Radushkevich (D-R) Isotherm
The D-R model in general takes into account the effect of the porous structure of the adsorbents [36]. The model is based on adsorption potential theory considering that the adsorption process is related to micropore volume filling [37], and Gaussian adsorption energy distribution is considered [17]. The D-R isotherm is expressed by Equation (9): q e = q m exp −βε 2 (9) where β (mol 2 kJ −2 ) is a constant related to the adsorption energy, and ε (kJ mol −1 ) is the adsorption potential (also called the Polanyi potential). The Polanyi potential for gas adsorption corresponds to the change in the Gibbs free energy of an adsorbent after adsorption of 1 mol of gas [38]. The D-R isotherm has also been adapted for the adsorption of liquids. In this case, the term ε is defined by Equation (25), where c s is the saturation concentration of the solute in water, similar to the BET isotherm.
β can be used for an estimation of the average energy of adsorption (E) by Equation (11): The D-R isotherm is often used to distinguish between physical and chemical adsorption [39,40]. It was postulated [41] that if E ∈ (8.0-16.0) kJ/mol, then the adsorption process occurs by chemisorption, and if E < 8.0 J/mol, then the adsorption has physical character.
It should be mentioned at this point that the Polanyi potential is frequently used in the wrong form [22,41]: that does not respect a dimensional consistency and leads to the wrong description of experimental data, as discussed by Hu and Zhang [42].

Error Analysis
In the past a linear regression analysis has been mostly applied for analysis of experimental data obtained from the adsorption process due to its numerical simplicity. It has been used to quantify the distribution of adsorbates on sorbents and to verify the consistency of adsorption models and the theoretical assumptions of adsorption models [43,44]. However, recent studies have indicated that the error structure of experimental data is usually changed if the original nonlinear adsorption isotherms are transformed into their linearized forms [45]. For this reason, nonlinear regression analysis is currently often preferred, since it provides a mathematically rigorous method for determining adsorption parameters using the original form of isotherm equations [46,47]. The nonlinear regression includes the minimization of error distribution between the experimental data and the predicted isotherm based on its convergence criteria [48]. The most common error functions used in nonlinear analysis are summarized in Table 1.

Thermodynamics of Adsorption
The thermodynamic parameters of adsorption, particularly the isosteric Gibbs energy of adsorption (∆G o ad ), the isosteric enthalpy of adsorption (∆H o ad ), and the isosteric entropy of adsorption (∆S o ad ), are calculated from adsorption (equilibrium) constants, particularly from the Langmuir adsorption constant, in order to characterize the thermodynamics of adsorption at various temperatures [17]. However, the use of those constants should be made with caution. The first problem is that the equilibrium constant in the original Langmuir isotherm is not dimensionless, whereas the equilibrium constant K 0 in Equation (13) has no dimension [17].
Equation (15) is the van 't Hoff equation. This equation serves to calculate the thermodynamic parameters of sorption from the linear regression of the experimental dependence lnK 0 (T) on 1/T. The thermodynamic parameters are expressed in J/mol. It is evident that both terms on the right side of the Equation (15) are dimensionless, and thus, K 0 must be dimensionless as well. This means that if K L values obtained from the LI and Equation (15) are used, then K L must first be transformed into a dimensionless parameter. The values of K L published in the literature have various units (L/g, L/mg, L/mol, L/mmol) [49][50][51], and thus, some consensus is needed on how to transform the dimensional K L into the non-dimensional K L . Recently, some approaches addressing this issue were proposed. The first approach is based on classical thermodynamics. The K L value involved in the LI for gases is also not dimensionless (the unit is Pa −1 or its equivalent); however, if it is used as the equilibrium constant K 0 in Equation (13), K L is multiplied by the standard reference pressure, which is commonly 1 bar. This does not change the value of K L but makes it dimensionless. A similar principle was applied for the K L constant in the LI for the adsorption of liquids. Let us consider solute A in liquid phase (in solution) A (liquid) and adsorbed at sorbent A (solid). In the equilibrium A(liquid)↔A(solid), the equilibrium constant can be expressed as K = a solid /a liquid , where a solid and a liquid are activities in the related phase at equilibrium. The equilibrium constant of the A(liquid)↔A(solid) reaction is given by Equation (16) [52,53]: (16) where c solid and c liquid are the concentrations in two phases, c 0 solid and c 0 liquid are the same concentrations in the selected standard states, and γ solid and γ liquid are activity coefficients. For simplicity, in dilute solutions of nonelectrolytes, γ solid , γ liquid ≈ 1. If the influence of the activity coefficients is neglected and the validity of Equation (16) is respected, the last remaining problem is the determination of values related to the standard concentrations. Unfortunately, there is no generally accepted consensus in the scientific community in the case of adsorption from liquid solutions. The simplest choice is to set standard concentrations of 1 mol/L or 1 mol/kg for both components, depending on how a concentration in equilibrium is expressed. More practically, the concentration can also be 1 mmol/L or 1 mmol/kg. In this case, c e must also be expressed in the same units. This approach can be applicable for true solutions; however, it is not physically correct to define a standard state for emulsions in this way. Other methods that consider different standard states or just nullified the dimension have been published [54][55][56][57][58] and critically discussed [59].
A different approach proposed by Azizian [59] is based on the reformulation of the LI in such a way that it directly leads to a dimensionless constant. The modified Langmuir isotherm based on this assumption is expressed by Equation (17): This approach gives an equilibrium adsorption constant K ML = k a /k d (Equation (18)) as the dimensionless parameter, unlike the common LI, where K L depends on the units used (it is reciprocal to the unit in which a concentration of solute in equilibrium is expressed). Equation (17) also predicts that q m is reached when c e = c s instead c e →∞ as a result of the original Langmuir isotherm (Equation (2)).
However, the introduction of parameter c s in general into the models (including BET and D-R) is questionable. It is evident that parameter c s , which may be significantly higher than c e , strongly influences the values obtained from the fitting of experimental data by the modified LI. Second, this model is not applicable for substances that do not form true solutions, such as emulsions, for which the saturation state cannot be unambiguously defined. This means that the dimensionless K ML would also be calculated from Equation (18) for the adsorption of emulsions because no standard state is required; however, experimental data q e = f(c e ) cannot be fitted by Equation (17) because the c S value is undefinable for emulsions. It may be kept as an adjustable parameter, however, without any physical meaning. The modified BET isotherm (Equation (6)), as proposed by Ebadi et al. [30], looks to be a good approach for evaluating the adsorption of liquids and emulsions because (i) it is a physically plausible model that lacks the involvement of any speculative parameters, (ii) it is applicable for multilayer adsorption, (iii) it does not require the use of the constant c s , which is an undefinable parameter for emulsions, and (iv) it enables calculations of equilibrium constants and adequate thermodynamic parameters.

Kinetics of Adsorption
Kinetic models serve to estimate the duration of adsorption processes and thus to estimation the time needed for the effective treatment of liquids. In batch systems, the solute concentration in the treated liquid gradually decreases with time until it reaches equilibrium with the adsorbed species. In percolating fluid systems, the solute concentration in the treated liquid is less constant and changes after passing a filter medium. The required outlet concentration of the solute should maintain some limit values (per internal specifications).
The most common kinetic models cited in the literature are the pseudo-first-order kinetic model (PFOM) and pseudo-second-order kinetic model (PSOM) for batch systems and the Thomas model and bed depth service time (BDST) model for continuous filtration processes. For the majority of common sorbents, the adsorption process is accomplished with diffusion into the bulk medium; therefore, diffusion is a key rate-limiting step in most adsorbent/solution systems. Some diffusion models, such as Douven's model and the Weber-Morris intra-particle diffusion model, will be discussed in this regard.

Pseudo-First-Order Kinetic Model (PFOM)
The PFOM, derived by Lagergren [60], is based on the simplified precondition that the rate of adsorption is controlled by an adsorbed amount (concentration of occupied sites on the sorbent) rather than a concentration of adsorbing substances in solution. The first-order kinetics implies that one element interacts with one unoccupied reaction site on the sorbent, and this process is described by Equation (19).
where q is the amount of adsorbed species per mass of adsorbent (mg/g), k 1 (min −1 ) is the pseudo-first-order rate constant, q e is the amount of adsorbed species per mass of adsorbent in equilibrium (mg/g), and t is time (min). An integrated form of Equation (1) gives Equation (20): which enables the determination of parameters q e and k 1 from nonlinear fitting. On the other hand, the linear form given by Equation (21) is mostly applied and cited in the literature. ln q e − q = lnq e − k 1 t , where c 0 (mg/L) is the initial concentration, c e (mg/L) is the concentration in equilibrium, V (L) is the volume of the investigated liquid and m (g) is the mass of the sorbent. The units mentioned here are the most common units referred to in the literature, but different units can be used as well.

Pseudo-Second-Order Kinetic Model (PSOM)
The PSOM represents the most common model for a quantitative description of the sorption rate of substances onto sorbent media [61,62]. The PSOM is expressed by the second-order differential equation (Equation (23)): where k 2 is the pseudo-second-order rate constant. Unlike k 1 , which has always the dimension reciprocal to time, the constant k 2 may have various dimensions (mg/g.min, g/g.min, mmol/g.min, etc.), depending on the definition of q [20,63]. Analytical solutions of Equation (23) can be expressed in various forms [64]; however, the most common solutions are given by Equations (24) and (25), which serve for a nonlinear (Equation (24)) and linear (Equation (25)) fitting of experimental data and a calculation of q e and k 2 .
One of the crucial assumptions for a derivation of the PSOM is that an initial concentration of bulk solution should not change significantly during sorption in order to evaluate the experimental data correctly [65]. On the other hand, this requirement is in direct contradiction to a common practice to determine q and q e from changes in the bulk concentration (c e ), particularly in batch systems, where the initial concentration changes significantly over the sorption experiment.
Azizian [43] proposed the general kinetic model, in which the adsorption rate is described by (i) a fraction of vacancies (available reaction sites) at the surface of the sorbent (1−θ), (ii) a concentration of solute in solution (c) in time t, and (iii) the desorption rate, which is given by a fraction of occupied sites at a sorbent surface (θ, where 0 ≤ θ ≤ 1) [66].
The general differential equation of that model is expressed in the form Equation (26): where c = c 0 − βθ, β = mq m M w V , or β = c 0 −c e θ e . c 0 is the initial molar concentration of the solute, q m is the maximum sorption capacity of the adsorbent, M w (g/mol) is the molar weight of the solute, V (L) is the volume of the solution, c e is the equilibrium molar concentration of the solute and θ e is the equilibrium coverage fraction. From the general model (Equation (26)), Azizian derived the PFOM (Equation (20)) and PSOM (Equation (25)) as special cases of the general model (Equation (26)).
Both models describe experimental data for both low and high solute concentrations, highlighting that the sorption process at high initial concentrations obeys PFO kinetics, while the PSO model better describes the sorption kinetics at lower initial concentrations [43,[67][68][69].
Various authors have analyzed a description of experimental data widely published in the literature by the PFOM and PSOM, and they identified the following problems, which indicate doubts about the general applicability of these models without considering additional phenomena [64,[66][67][68][69][70]: (i) unoccupied sites can be incapable of reacting with adsorbing entities independently, (ii) adsorption processes can take significantly longer than expected (electrostatic interactions between the surface and sorbents occur instantaneously and the diffusion layer at the solid/liquid interface is thin enough), and (iii) the initial concentration of substances in solution changes significantly during the experiment, particularly in the case of batch systems. More details about the limitations of the PSOM can be found in Hubbe's review [64].
Regarding the applicability of the PSO and PFO models to experimental data, Canzano [71] pointed out that the correlation between data and models characterized by R 2 values strongly depends on the selection of data. If the data set does not include enough values related to the short duration of the experiment but consists of many more data for long periods, the long-term data have a higher weight on the statistical parameters. Canzano [71] showed that the PSO model better fits data including measurements from systems that approach an equilibrium. Another problem where mathematical procedures can overlap the physical reality is associated with a selection of the equation used from a fitting. As mentioned above, an integrated form of Equation (5) can be expressed by four different equations. It was shown that application of those equations to the same set of data leads to significantly different R 2 values in the range from 0.862 to 1. On the basis of this finding, the authors proposed, first, to minimize the number of data points nearest to equilibrium and, second, to use nonlinear fitting by Equation (24)  Diffusion is a key rate-limiting step in most adsorbent/solution systems unless the surfaces are smooth and impervious; therefore, this phenomenon has to be taken into account in the majority of sorption processes [48].
Douven's model [72] represents a generalization of Crank's model [73] based on the analytical solution of 2nd Fick's law with some preconditions such as (i) the formation of a boundary layer at the solid/liquid interface, (ii) the reversible adsorption of solutes, (iii) the continuous occurrence of adsorption during the diffusion process, and (iv), an effective diffusion coefficient is considered instead a "true" diffusion coefficient, taking into account a decrease in the rate of diffusion due to solute immobilization on the adsorbent surface. Douven's model is generally expressed by Equation (27): where n(t) is the total number of adsorbed entities at time t, n e is the total number of adsorbed entities at equilibrium, D eff is the effective diffusion coefficient, and R p is the radius of the adsorbed entity. Douven also developed the procedure for data analysis based on the experimentally determined half-times (t 1/2 ) needed for an adsorption of half of the full adsorbent capacity and showed that the PSO model is equivalent to Equation (28):

Weber-Morris Intra-Particle Diffusion (IPD) Model
The IPD model [74] is the most frequently applied model for systems in which adsorption is accomplished by the diffusion process and is described by Equation (29): where k ipd (mg/g.min) is the rate constant for intra-particle diffusion, and C (mg/g) is a constant related to the boundary layer thickness. This model does not distinguish among various processes that can simultaneously occur during solute sorption and considers just a single rate-limiting process of diffusion.

Adsorption Kinetics in Flow-through Systems (Fixed Bed Adsorption)
When contaminated fluid moves through a fixed bed, the pollutant to be adsorbed transfers from the bulk fluid to the adsorbent bed. Several steps are involved in the overall adsorption process of a single molecule of pollutants:

1.
A mass transfer from the bulk of the fluid to the surface of the adsorbent through the boundary layer around the particle.

2.
An internal diffusion through the pores of the adsorbent.

3.
An adsorption onto the surface of the adsorbent.
In wastewater treatment applications, the overall adsorption process is dominated by mass transfer, especially by an intra-particle mass transfer (Figure 4), dependent on the morphology of sorbents.  During a continuous adsorption process in a fixed bed, the bed can be divided into three distinguished zones in (i) the saturated zone, (ii) the adsorption zone, and (iii) a clean zone in which the adsorbent contains little or no adsorbed pollutant ( Figure 5). The size and location of these three zones within the bed change with time. When the concentration of pollutant on the bottom part of the filtration column reaches some critical level of saturation, concentrations of pollutants in an effluent are higher than the required value. This point is called the breakthrough point and the corresponding curve of pollutant concentration in the effluent versus time is called the breakthrough curve. In this phase, the filtration process is stopped and the adsorption medium is cleaned. Otherwise, the concentration of pollutant in effluent rapidly increases until it reaches an exhaustion point when the column approaches saturation, and the effluent releasing the bed has the same concentration of pollutant as the influent stream.    During a continuous adsorption process in a fixed bed, the bed can be divided into three distinguished zones in (i) the saturated zone, (ii) the adsorption zone, and (iii) a clean zone in which the adsorbent contains little or no adsorbed pollutant ( Figure 5). The size and location of these three zones within the bed change with time. When the concentration of pollutant on the bottom part of the filtration column reaches some critical level of saturation, concentrations of pollutants in an effluent are higher than the required value. This point is called the breakthrough point and the corresponding curve of pollutant concentration in the effluent versus time is called the breakthrough curve. In this phase, the filtration process is stopped and the adsorption medium is cleaned. Otherwise, the concentration of pollutant in effluent rapidly increases until it reaches an exhaustion point when the column approaches saturation, and the effluent releasing the bed has the same concentration of pollutant as the influent stream.  The performance of column breakthrough can be analyzed based on the Henry C. Thomas model (Equation (30)) [75,76], which was originally developed for solute chromatography. This is one of the most utilized models for evaluating column performance and can be applied for a prediction of adsorption kinetics and the adsorption capacity in flow-through systems [77]. In this model, the intra-particle diffusion during the mass transfer process is considered to be negligible.
Equation (30) can be expressed in linear form as Equation (31): where C 0 is the feed concentration (mg/L), C e is the column effluent concentration (mg/L), K th is the kinetic rate constant (L/(mg.min)), q m is the adsorption capacity of the column (mg/g), M is the mass of media in the bed (g), V is the column effluent volume (L) and Q is the feed flow rate (L/min). The breakthrough data can be fitted by plotting ln C 0 C e − 1 versus V; q m and K th can be determined depending on the intercept and slope [78].
Pintor et al. [79] developed the model given by Equation (32)a,b,c for testing the adsorption of oil from oil/water emulsions.
where c b0 (mg/L) is the initial concentration of oil in the bulk solution, ξ m is the batch capacity factor, τ (s) is the time constant characterizing diffusion through a boundary layer, V (L) is the volume of the solution (a volume of a reactor), ε b is the bulk porosity, K f (cm/s) is the mass transfer coefficient of the boundary layer, and a f (cm 2 /cm 3 ) is the external surface area of sorbent per volume of sorbent given by the relation a f = 3/r p , where r p (cm) is the particle (sorbent) radius. The adjustable parameters τ f and K f are obtained from a fit of experimental data by Equation (32).

Bed Depth Service Time (BDST) Model
The BDST model predicts a dependence of the service time (t b (min)) at the breakthrough point in the depth (X, (cm)) of the packed column. In this model, intra-particle diffusion is neglected, and direct adsorption of solutes on the adsorbent surface is assumed. The original BDST model was proposed by Bohart et al. [80] in the following form (Equation (33)): Hutchins [81] proposed a linearized form of that model, expressing a direct relation between tb and X by Equation (34).
N 0 (mg/L) is the dynamic bed capacity, v (cm/h) is the linear flow rate, defined as the ratio of the volumetric flow rate Q vol (cm 3 /h) to the cross-sectional area of the bed, S c (cm 2 ); c 0 and c b (mg/L) are the initial and breakthrough concentrations of the solute, respectively, and k ads (L/mg.h) is the adsorption rate constant.
From a practical point of view, the breakthrough curves can be described using various commercial and non-commercial software. For example, freely available software FAST, Fixed-bed Adsorption Simulation Tool () is the program for the prediction of breakthrough curves of fixed-bed adsorption filters used in water treatment. Inputs include an empty bed contact time, mass of adsorbent, bed porosity, particle density, particle diameter, influent concentration, flow rate, bed volume. The outputs are Freundlich isotherm parameters, film diffusion coefficient (m/s), surface diffusion coefficient (m 2 /s), breakthrough curve, solute distribution parameter (Dg), Biot number (Bi), and Stanton number (St).
Some examples of the separations of o/w emulsions as well as mostly used models are shown in Table 2.

Interfacial interactions and Wettability of Surfaces
The physical parameter that describes energy/force-related phenomena that act on the outer surface of a material or at the interface between adjacent bodies is called the specific free surface energy (σ (J/m 2 )) or the surface tension (γ (N/m)). Except for intramolecular (covalent, ionic) bonds, intermolecular interactions, generally called Van der Waals interactions, contribute to both internal and interfacial energy. The energy of those interactions strongly decreases with the distance between interacting bodies (~1/r 6 ) and with temperature (1/T) [17]. The range of the forces is approximately from 0.15 to 1 nm, and the energy depends on the type of interactions: (i) dipole-dipole interactions~5 to 25 kJ/mol, hydrogen bonds~10 to 40 kJ/mol, (ii) ion-dipole interactions~40 to 600 kJ/mol, (iii) ion-induced dipole interactions~2 to 10 kJ/mol, and (iv) London (dispersion) interac-tions~0.5 to 40 kJ/mol. These interactions are time-independent and act instantaneously; therefore, the time dependence of an adsorption process is governed by additional, mostly diffusional phenomena [96].
The most common approach for characterizing interfacial tension is based on the sessile drop method, which is a measurement and evaluation of the contact angles of liquids on a solid surface at equilibrium at room temperature and atmospheric pressure. In a sense of water/oil separation, not only are solid/oil/air interfaces but also mostly solid/oil/water interfaces investigated, and related contact angles at different interfaces indicate the suitability of separation methods for different morphologies of oil/water systems (emulsions, free oils). The typical features of the wettability of smooth and rough surfaces by oil-in-water surroundings are briefly described below.

A Droplet Placed at a Smooth and Rough Surface in Air
The fundamental dependence between the equilibrium contact angle and related surface tensions is given by Young's equation (Equation (35)), which describes the equilibrium state of a small droplet of pure liquid deposited on a uniform, smooth and non-penetrating surface [97][98][99][100]: where Θ is the contact angle, γ SV is the interfacial tension at the solid-vapor interface, γ SL is the interfacial tension at the solid-liquid interface, and γ LV is the surface tension of the liquid, which is in equilibrium with a saturated vapor phase.
Other types of surfaces are those that have a fine roughness (surface is interrupted by cavities with sharp edges having micro/nano inter-cavity distances); therefore, those cavities are not wettable by liquids but are filled with air. This case is called the Cassie state and is characterized by the Cassie and Baxter equation (Equation (36)) [101]: where f is the fraction of the solid surface that is in contact with the liquid. The assumption involved in Equation (36) is that the contact angle of the droplet at the air layer is 180 • . The Cassie state corresponds to the very low adhesion between liquid and solid, and therefore is usually unstable. The transition between the Wenzel case and the Cassie-Baxter case depends on both the chemical character of the surface and the geometry of the roughness, particularly the scale of the roughness ( Figure 6). Microscale roughness leads to strong contact angle hysteresis effects. Contact angle hysteresis cannot be determined directly but can be characterized via the measurement of advancing and receding contact angles. The contact angle hysteresis is defined as the difference between the advancing contact angle (ACA) and the receding contact angle (RCA) of a droplet that moves on the tilted surface. The differences between static sliding angle (SSA) and dynamic sliding angle (DSA), similar to those between ACA and RCA, may range from a few degrees to a few tens of degrees, and SSA > DSA [102,103]. The relation between the contact angle hysteresis and sliding angle is described by the Furmidge equation (Equation (37)) [38]: where m is the weight of a droplet, α is a sliding angle, w is the width of the droplet, and Θ adv and Θ rec are advancing and receding contact angles, respectively. hysteresis is defined as the difference between the advancing contact angle (ACA) and the receding contact angle (RCA) of a droplet that moves on the tilted surface. The differences between static sliding angle (SSA) and dynamic sliding angle (DSA), similar to those between ACA and RCA, may range from a few degrees to a few tens of degrees, and SSA  DSA [102,103]. The relation between the contact angle hysteresis and sliding angle is described by the Furmidge equation (Equation (37)) [38]: where m is the weight of a droplet,  is a sliding angle, w is the width of the droplet, and adv and rec are advancing and receding contact angles, respectively.  [104][105][106] In this case, the wettability of surfaces by oil is significantly influenced by the contribution of the surface tension of water because, unlike air, which has no surface tension, water has one of the highest surface tension values (~72 mN/m). In this case, considering

A Droplet Placed at a Smooth and Rough Surface in Water
In this case, the wettability of surfaces by oil is significantly influenced by the contribution of the surface tension of water because, unlike air, which has no surface tension, water has one of the highest surface tension values (~72 mN/m) [104][105][106]. In this case, considering a three-phase solid substrate (S)/oil (O)/water (W) system, Young's equation can be rewritten into the forms given by Equations (38)-(40) [107]: where γ OA , γ WA , and γ OW are the oil/air, water/air, and oil/water interface tensions, respectively, and Θ O , Θ W , and Θ OW are the contact angles of oil in air, water in air, and oil in water, respectively. Equations (38)- (40) indicate that a hydrophilic surface in air is also oleophilic in air because γ OA << γ WA , and hydrophilic surfaces in air behave as oleophobic in water, as is evident from Equation (40) [107].
If the roughness of a surface is considered, then equations for the Wenzel state and Cassie state can be used to describe the surface wettability using Equations (41) and (42): where Θ OW* and Θ OW are the contact angles on the oil droplet on the rough surface and the smooth surface in the water surroundings, respectively. The Wenzel state corresponds to the situation where an oil droplet is attached to a rough, superhydrophobic surface, which has pores filled with air, and capillary forces are strong enough to suck adjacent oil into the pores [106]. In this case, the surface is fully wetted by oil. The Cassie state corresponds to the situation where valleys on the surface are filled with water, which suppresses the penetration of oil into those cavities, and oil droplets can only be in contact with the "pin" objects. This situation leads to an enhanced oleophobicity of such surfaces [107]. The capability or incapability of water to fill cavities at a given surface plays an important role in the final wettability of the surface by oil. This capability is given by the surface free energy and surface topology. Superhydrophobic surfaces are characterized by a static contact angle of water of over 150 • , and the rolling angle is below 6-10 • . Superhydrophobic or superoleophilic behavior is a consequence of a combination of the chemical structure of the material and controlled roughness [108][109][110][111].
The above considerations are summarized in Figure 7 [52].
where OW* and OW are the contact angles on the oil droplet on the rough surface and the smooth surface in the water surroundings, respectively. The Wenzel state corresponds to the situation where an oil droplet is attached to a rough, superhydrophobic surface, which has pores filled with air, and capillary forces are strong enough to suck adjacent oil into the pores [106]. In this case, the surface is fully wetted by oil. The Cassie state corresponds to the situation where valleys on the surface are filled with water, which suppresses the penetration of oil into those cavities, and oil droplets can only be in contact with the "pin" objects. This situation leads to an enhanced oleophobicity of such surfaces [107]. The capability or incapability of water to fill cavities at a given surface plays an important role in the final wettability of the surface by oil. This capability is given by the surface free energy and surface topology. Superhydrophobic surfaces are characterized by a static contact angle of water of over 150°, and the rolling angle is below 6-10°. Superhydrophobic or superoleophilic behavior is a consequence of a combination of the chemical structure of the material and controlled roughness [108][109][110][111].
The above considerations are summarized in Figure 7 [52].

Free (Stratified) Oil Separation from Water (Oil Spills): Superhydrophobic and Superoleophilic Surfaces
This treatment is applicable for the separation of free oil (oil spills) through sorption. For this reason, various types of particles, mats and bulk polymeric sorbents have been treated and tested [2,3,112]. Surface modification is mostly performed on traditional sorbents such as powders (silica, clays, activated carbon, ground cellulose-based bio-waste including kapok, wood dust, kenaf, etc.), synthetic polymer fabrics (polypropylene), and polymer sponges (polyurethane, melamine, synthetic rubbers, etc.). Surface treatment improves soil sorption and reduces water permeability into the materials [106].

Oil/Water Emulsions-Separation of Oil Droplets by Membrane Filtration. Superoleophobic and Superhydrophilic Surfaces
Ultra-and microfiltration technologies, which use suitable membranes, are commonly employed for oil/water emulsion separation. The drawback of membrane filtration is associated with intensive fouling of their surfaces, which leads to a decrease in the flow rate through the membranes and enhances the time and cost needed for cleaning. Designing membranes with superoleophobic and superhydrophilic surfaces enables suppression of that unfavorable behavior. From a chemical point of view, materials should have a surface free energy in the range from 20-30 mN/m (common values for oils) to~72 mN/m (surface tension of water) and characteristic dimensions ranging from several tens to hundreds of microns. The low surface tension of materials is given by their composition; among the materials with the lowest surface energy are polydimethylsiloxane (19.8 mN/m), polytetrafluoroethylene (22.6 mN/m) and polypropylene (29.6 mN/m), which means materials where the cohesion energy is given by London dispersive forces. The roughness can be introduced into membranes through various physical, chemical and processing approaches [38,[104][105][106][107]. The Cassie wetting state corresponds to this filtration mechanism.

Oil/Water Emulsions-Separation of Oil Droplets by Sorption (Preferable in Batch Systems). Superhydrophobic and Superoleophilic Surfaces
Strategies similar to those for the separation of free oil from water can also be applied for the separation of O/W emulsions in batch systems. The oil component adsorbs at the sorbent surface and, in some cases, penetrates into the bulk structure, where it is localized at the internal surface of sorbents. The separation efficiency is given by the surface area available for sorption. Once a surface area is fully covered, the separation process stops. The Wenzel wetting state corresponds to the situation where the oil droplet is attached to a rough, superhydrophobic surface, which has pores filled with air, and capillary forces are strong enough to suck adjacent oil into the pores. In this case, the surface is fully wetted by oil, and the surface and adjacent pores are fully utilized for oil sorption.

Oil/Water Emulsions-Separation of Oil Droplets by Demulsification-Coalescence (Deep-Bed Filtration). Balanced Hydrophobic and Oleophilic Surfaces
Two basic filtration media for deep-bed filtration are used in industry. The first type of medium includes fibrous, porous woven and nonwoven mats. If the dominant mechanism of oil separation is demulsification via coalescence, these media are called coalescers. The majority of commercial coalescers separate droplets with diameters less than 100 µm, and if they are formed by a combination of different fibrous mats that vary in composition, thickness, and pore size, they can separate droplets smaller than 10 µm [113][114][115][116][117].
The separation of droplets smaller than 10 µm is a challenging task because small droplets behave as solids due to the high internal Laplace pressure and have a small contact area with fibers within mats; therefore, coalescence occurs insufficiently. An improvement of fibrous coalescers is based on the design of sandwich structures formed by a combination of different mats varying in the type of material, fiber thickness, pore size, and specific surface treatment, leading to tailored surface porosity and roughness [113][114][115][116][117].
The second type of medium for deep-bed filtration includes various organic, synthetic and inorganic granules, possessing roughly spherical shapes and sizes on the order of millimeters. All these media adsorb oil to various extents, and if the surface chemistry and topology are favorable, they initiate coalescence of oil droplets at their surfaces and therefore significantly contribute to the separation efficiency [118]. If coalescence is promoted at the filter surface, the droplets, which are attached within a short mutual distance for some time, have favorable conditions to merge with each other. Agarwal et al. [119,120] tested the adhesion of droplets of iso-octane on surfaces of different polarities (the surface free energy), porosities, and roughness and searched for suitable conditions, which led to the efficient coalescence of such small droplets. The work of adhesion (Equation (43)) is considered the parameter characterizing the ability of droplets to be attached onto a surface and to remain there for some time. The dependence of the work on adhesion on the surface free energy of the filler is shown in Figure 8. W ad = γ Ly 1 + cos Θ Ly (43) where γ Ly is the interfacial tension between the dispersed liquid (y-oil) and the surroundings (water, or y-air). The dependence of the ratio of W ad (water/iso-octane) on the surface energy has a sigmoidal character, and the curve characterizes two distinguished regions: the region with low surface energy in which the adhesion of iso-octane droplets is stronger than the adhesion of water (below~25-27 mN/m) [119] and the region with high surface free energy (over 30 mN/m) where the adhesion of water is stronger than the adhesion of oil. Materials with high-energy surfaces are not suitable for oil adsorption because they cannot replace adsorbed water. For this reason, the surfaces of coalescing media should be hydrophobic and oleophilic, which enable good wetting and sufficient contact time, thereby promoting the coalescence of droplets [119]. The second important point is the morphology of the surfaces, characterized by roughness (quantified by a dimensionless roughness factor) and porosity (quantified by the pore size in nm/µm). As discussed above, the surface topology enhances the intrinsic wettability of flat surfaces. However, these effects may not be favorable for droplet coalescence. Agarwal et al. [119] demonstrated that nanoscale roughness improves oil droplet spreading and adhesion at hydrophobic surfaces (Wenzel wetting), while nanocoated hydrophilic surfaces lead to an enhanced contact angle and improved mobility, characterized by low static and dynamic sliding angles (Cassie state) in water surroundings ( Figure 9). This effect may also contribute to emulsion formation, however, outside the solid surface, which takes a long time, especially for diluted, stabilized emulsions, and the separation efficiency is low. The Wenzel state enhances the probability of mutual contacts between oil droplets at a surface with high coalescence and separation efficiency. A pore size of cavities larger than the drop size, as well as high roughness, is favorable for coalescing media. On the other hand, oil absorption within porous oleophilic media is not favorable for coalescence [119]. The same authors in another paper [120] demonstrated those findings on woven and nonwoven fabrics treated by silica nanoparticles, and the influence of the fabric thickness, fiber diameter, porosity, pore size, and microscale roughness were systematically investigated. The best configurations enabled the separation of emulsions consisting of oil droplets with a diameter of 5 μm, with a separation efficiency of 80% at a pressure drop of 160 mbar and a superficial velocity of approximately 1.0 m/min.
The physical and chemical principles discussed above are the same for both fibrous coalescers and granular media; however, there are some differences resulting from different media geometries. First, the interstices (pores) between fibers in mats are much smaller than the interstices (free volume) in the bed formed by granules. From simple geometry it can be estimated that the pore size has the same order or one order of magnitude less than the diameter of fibers/granules. Common fabrics have pore sizes on the order of 100-101 μm, whereas pores within granular beds have a size of 100 mm. The smaller space is more favorable for coalescence, which can occur at a surface and in a solution; however, on the other hand, the space is filled faster with residual droplets, which leads to the pressure drop. For the granular bed, because the space of interstices is relatively large, coalescence proceeds at the surface of granules only. The advantage is a lower influence on the pressure drop and a less intensive backwashing procedure [121,122].
Granular nutshells, particularly walnut shells, are the most common media in the petrochemical industry for oil impurity removal from produced water. These filters are considered the standard technology, mainly for onshore use. The separation mechanisms include oil adsorption on the surface and coalescence. It was shown that droplets act at first as discrete particles, which are trapped by straining, and then the process continues by coalescence of those larger droplets, leading to the formation of large pools of oil localized within the filter interstices. An oil emulsion (unrefined Bakken crude oil; 100 ppm) was filtered through the bed (a height of 23 inches and diameter of 4 inches) over six hours. It was found that 95% of the oil was trapped in the first 12 inches of the bed and 99% in the first 18 inches of the bed. Large pools of coalesced oil were detected within the medium interspace [118].
Coalescing agents [123] have been used in deep-bed filtration using sand [124] and sawdust [125] as the filter media. Cetyltrimethylammonium bromide as a coalescing agent has been added for the separation of hexadecane/water emulsions by using a quartz sand bed packed in a Plexiglas column with an inner diameter of 2.5 cm and a length of 10 cm [124]. An n-hexadecane/water emulsion at a concentration of 1000 ppm was pumped from the bottom at a superficial velocity of 0.04 cm/s using a peristaltic pump. It was concluded that the deposition of oil is found to be a maximum at a cetyltrimethylammonium bromide The same authors in another paper [120] demonstrated those findings on woven and nonwoven fabrics treated by silica nanoparticles, and the influence of the fabric thickness, fiber diameter, porosity, pore size, and microscale roughness were systematically investigated. The best configurations enabled the separation of emulsions consisting of oil droplets with a diameter of 5 µm, with a separation efficiency of 80% at a pressure drop of 160 mbar and a superficial velocity of approximately 1.0 m/min.
The physical and chemical principles discussed above are the same for both fibrous coalescers and granular media; however, there are some differences resulting from different media geometries. First, the interstices (pores) between fibers in mats are much smaller than the interstices (free volume) in the bed formed by granules. From simple geometry it can be estimated that the pore size has the same order or one order of magnitude less than the diameter of fibers/granules. Common fabrics have pore sizes on the order of 100-101 µm, whereas pores within granular beds have a size of 100 mm. The smaller space is more favorable for coalescence, which can occur at a surface and in a solution; however, on the other hand, the space is filled faster with residual droplets, which leads to the pressure drop. For the granular bed, because the space of interstices is relatively large, coalescence proceeds at the surface of granules only. The advantage is a lower influence on the pressure drop and a less intensive backwashing procedure [121,122].
Granular nutshells, particularly walnut shells, are the most common media in the petrochemical industry for oil impurity removal from produced water. These filters are considered the standard technology, mainly for onshore use. The separation mechanisms include oil adsorption on the surface and coalescence. It was shown that droplets act at first as discrete particles, which are trapped by straining, and then the process continues by coalescence of those larger droplets, leading to the formation of large pools of oil localized within the filter interstices. An oil emulsion (unrefined Bakken crude oil; 100 ppm) was filtered through the bed (a height of 23 inches and diameter of 4 inches) over six hours. It was found that 95% of the oil was trapped in the first 12 inches of the bed and 99% in the first 18 inches of the bed. Large pools of coalesced oil were detected within the medium interspace [118].
Coalescing agents [123] have been used in deep-bed filtration using sand [124] and sawdust [125] as the filter media. Cetyltrimethylammonium bromide as a coalescing agent has been added for the separation of hexadecane/water emulsions by using a quartz sand bed packed in a Plexiglas column with an inner diameter of 2.5 cm and a length of 10 cm [124]. An n-hexadecane/water emulsion at a concentration of 1000 ppm was pumped from the bottom at a superficial velocity of 0.04 cm/s using a peristaltic pump. It was concluded that the deposition of oil is found to be a maximum at a cetyltrimethylam-monium bromide concentration of 5 × 10 −6 M, at which the sand surface has an opposite charge to that of the emulsion droplets. Below this concentration, the deposition efficiency is reduced due to the electrostatic repulsion between the negatively charged droplets and the sand. Above this concentration, the deposition is again reduced due to the surfaces of sand and droplets being positively charged.
The coagulator calcium sulfate has been used to promote coalescence of oil/water emulsions at a concentration of 24 ppm in deep-bed filtration experiments. The filter media consist of sawdust with heights ranging from 30 to 55 mm packed in a 20-cm column, with the inner glass tube replaced by a 26-mm diameter methacrylate column. The dissolved salt ions compress the electric double layer around droplets, neutralize repulsive forces among them and between droplets and sawdust particles and hence enhance the coalescence of oil droplets. Moreover, the lack of electrostatic repulsion improves the physical adsorption of the oil droplets on the sawdust surface, originally negatively charged. Therefore, the oil droplet size is increased, and its charge is neutralized so that the disperse phase can be retained in the filter media interstices [125].

Conclusions
This review concerns the mechanisms of a separation of emulsified water/oil mixtures and it covers the (i) kinetics of adsorption, (ii) adsorption isotherms, (iii) interfacial interactions, and (iv) water/oil separation techniques based on the wettability of solid/oil/water interfaces. The advantages and drawbacks of commonly used as well as newly proposed kinetic and adsorption models were reviewed, and their applicability for the characterization of oil/water separation was discussed. It was pointed out that the direct applicability of models that can be applied to describe the sorption process from true solutions is not straightforward, and these models should be used with caution, especially if equilibrium constants of adsorption and related thermodynamic parameters are estimated. The wettability of solid sorbents by oil within a water environment was discussed for both smooth and rough surfaces. The main focus was on the surface energy and surface topology of solid media and their influence on the coalescence separation mechanism.
The separation of emulsions through deep-bed filtration is a complex technological process, and, therefore, not all aspects of this technology were covered in this review. Among those mechanisms were transport mechanisms include straining, interception, inertia, sedimentation, hydrodynamic action, filter and flow orientation, flow velocity, pressure drop, etc. Some relevant references concerning those areas were mentioned. For theoretical models that address the description of different phenomena related to the oil/water emulsion separation processes, let us mention just one aspect here that deserves attention in the future: the lack of suitable adsorption isotherms that can be correctly applied for a description of oil adsorption at external and internal solid surfaces of both nonporous and porous structures. A simple use of common isotherms is often misleading because the adsorption of oily components at solid surfaces does not fit the assumptions for which these models were originally derived. This bottleneck then results in problematic calculations of the thermodynamic parameters of sorption.

Conflicts of Interest:
The authors declare no conflict of interest.