Regular Article
Experimental study of water droplet break up in water in oil dispersions using an apparatus that produces localized pressure drops
^{1}
PETROBRAS – Petróleo Brasileiro S. A., CENPES/PDDP/TPP, Av. Horácio Macedo 950, Cidade Universitária, Ilha do Fundão, 21941915
Rio de Janeiro, RJ, Brazil
^{2}
UFRJ – Universidade Federal do Rio de Janeiro, Escola de Química, Av. Athos da Silveira Ramos 149, Cidade Universitária, Ilha do Fundão, 21941909
Rio de Janeiro, RJ, Brazil
^{3}
UNIFEI – NUSEC – Núcleo de Separadores Compactos, Instituto de Engenharia Mecânica, Av. BPS 1303, 37500903
Itajubá, MG, Brazil
^{*} Corresponding author: fabriciosoares@petrobras.com.br
Received:
17
July
2018
Accepted:
22
October
2018
Oil production facilities have choke/control valves to control production and protect downstream equipment against over pressurization. This process is responsible for droplets break up and the formation of emulsions which are difficult to treat. An experimental study of water in oil dispersion droplets break up in localized pressure drop is presented. To accomplish that, an apparatus simulating a gate valve was constructed. Droplet Size Distribution (DSD) was measured by laser light scattering. Oil physical properties were controlled and three different break up models were compared with the experimental results. All experimental maximum diameters (d_{max}) were above Kolmogorov length scale. The results show that d_{max} decreases with increase of energy dissipation rate (ε) according to the relation d_{max} ∝ ε^{−0.42}. The Hinze (1955, AIChE J.1, 3, 289–295) model failed to predict the experimental results, although, it was able to adjust reasonably well those points when the original proportional constant was changed. It was observed that increasing the dispersed phase concentration increases d_{max} due to turbulence suppression and/or coalescence phenomenon. Turbulent viscous break up model gave fairly good prediction.
© F.S. Silva et al., published by IFP Energies nouvelles, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The primary oil processing in a production facility aims at promoting produced fluids separation (oil, gas, brine and eventually solid sediments), in a way that the different phases can be obtained with a suitable purity for subsequent handling (Arnold and Stewart, 1999; Manning and Thompson, 1995). In the beginning of field exploitation, the production is basically oil and gas. Over the years, the oil field starts to produce an increasing amount of salty water. This brine can be connate water or be due processes of enhanced oil recovery (injection of seawater or produced water in the reservoir).
When the produced mixture reaches an offshore facility, it encounters a device called choke valve. Its function is to regulate the flow coming from the well and control its flux to guarantee an efficient yield, avoiding premature water production (Janssen et al., 2001). Another function of this valve is the protection of downstream equipment, by applying a pressure drop to adjust the flow pressure from the well for the downstream first stage separator operational pressure (Arnold and Stewart, 1999). Hence, this localized pressure drop promotes droplets break up since the flow energy is dissipated in a turbulent way. The result of this turbulent energy dissipation is stable dispersions or emulsions, which are very difficult to treat.
The design of separation stages and separators focuses in maximizing production and minimizing compressor horsepower (Arnold and Stewart, 1999; Manning and Thompson, 1995). It shall consider field characteristics, fluids properties, well flow pressure, temperature, GasOil Ratio (GOR), particle diameter and other variables (Lockhart et al., 1987; Manning and Thompson, 1995). The higher number of stages, more light components will be stabilized in liquid phase (Arnold and Stewart, 1999; Lockhart et al., 1987). Subsequent stages operate in lower pressures, therefore, besides the pressure drop in the choke valve, there are pressure reductions between separation stages. Those pressure drops are fulfilled by control valves and, therefore, they are also drop break up generators.
There are innumerous articles that study the break up phenomenon in pipe flows (Angelli and Hewitt, 2000; Kubie and Gardner, 1977; Simmons and Azzopardi, 2001). Other works studied this phenomenon in static mixers or rotor stator mixers (Grace, 1982; Lemenand et al., 2003; Rueger and Calabrese, 2013a, b).
Some studies tried to better understand droplet break up in flows passing through localized restrictions (Fossen and Schumann, 2016; Galinat et al., 2005; Janssen et al., 2001; Mitre et al., 2014; Paolinelli and Yao, 2017; van der Zande, 2000; van der Zande et al., 1998a, b, c, 1999a, b).
Most of the studies in this area were performed with low concentration of the dispersed phase (normally dispersions of oil in water) and/or with localized pressure drops below those observed in a production field. The experimental study, here presented, is inserted in the investigation of droplet break up in water in oil dispersions, with dispersed phase concentration varying between 5% and 10%, and localized pressure drops varying from 3 to 9 bar.
2 Theory
A droplet in a dispersion is subjected to forces applied by the continuous phase that is surrounding the droplet (Hinze, 1955; Taylor, 1934).Thus, normal and tangential tensions deform the droplet surface (Lemenand et al., 2003). The force per area that opposes deformation, known as Laplace pressure (P _{L}), is proportional to the interfacial tension (σ) and inversely proportional to the curvature radii of the interface (R _{1} and R _{2}), as shown in equation (1).(1)
For spherical droplet, R _{1} = R _{2} = d/2. From equations (1) and (2) can be written, where d is the drop diameter:(2)
The third force per area acting in the drop deformation process and break up is related to the viscous tensions (σv) generated internally in the drop, by the external velocity gradient (continuous phase), which, according to Hinze (1955), can be described by the following expression, where μ _{d} is the droplet viscosity, τ is an external tension and ρ _{d} is the droplet density:(3)
The viscous tension and the interfacial tension of the droplet oppose the external tension in the break up process. Therefore, it is assumed that in order to have droplet break up, the external forces (imposed by continuous phase flow) must be greater than the sum of viscous and interfacial tensions (Hinze, 1955).
The dimensionless number of Weber (We – Eq. (4)) is the ratio between the deforming force per area applied by the continuous phase on the droplet and the interfacial tension that resists to this deformation. The higher the Weber number, higher will be the external tension compared to the internal tension. In a critical value of the Weber number (We_{crit.}), the droplet break up can occur (Hinze, 1955; Lemenand et al., 2003)(4)
2.1 Break up in a laminar flow
Lemenand et al. (2003) have done a review of Taylor’s (1932, 1934) study, which demonstrates that the droplets break up, in a laminar flow, is governed by viscous shear stress.
Following Walstra and Smulders (1998), when the external forces have origins in viscous stress, the Weber number assumes equation (5) form and can be named Capillary number, where G is the droplet surrounding velocity gradient and μ _{c} is the continuous phase viscosity:(5)
The critical Weber number is dependent of the flow regime and phases viscosity ratio: R _{visc} = μ _{d}/μ _{c}, in which, for viscosity ratios bigger than approximately four (R _{visc} ~ 4), the droplets acts like a solid sphere and no break up can be observed. On the other hand, whenever R _{visc} is close to zero (0.01 for example), a huge shear stress is required for droplet breaking up. To break a drop, a viscosity ratio between 0.1 and 1.0 is preferable (Grace, 1982; Lemenand et al., 2003; Walstra, 1993; Walstra and Smulders, 1998).
2.2 Break up in a turbulent flow
Kolmogorov and Hinze theories brought important knowledge to the study of droplet break up in a turbulent flow. This regime is characterized by the presence of eddies in various unpredictable directions, that constantly change by mutual interactions. It makes a turbulent local instantaneous velocity (u) variation, in a chaotic way, leading to local velocity fluctuations (u′). The temporal mean of these fluctuations (u′) is equal to zero (Walstra, 1993).
A turbulent flow presents a large spectrum of eddies length (λ). The bigger the vortices length, the higher are its u′. These bigger eddies transfer kinetic energy to smaller ones that, nonetheless, presents higher velocity gradients. The velocity fluctuation for an eddy (with determined length) can be correlated with the energy dissipation rate per unit of mass (ε), by equation (6) (Davies, 1985; Lemenand et al., 2003; Walstra, 1993; Walstra and Smulders, 1998):(6)
The diameter of the smallest eddy is called Kolmogorov (λ _{K}) lengthscale. Vortices smaller than λ _{K} are dissipated by the viscosity of the fluid, transforming the kinetic energy into heat. Eddies, bigger than Kolmogorov lengthscale, are responsible for droplet break up (Walstra, 1993). The maximum stable droplet diameter (d _{max}) in an emulsion is determined by the equilibrium between the dynamic pressure generated by turbulent fluctuations, that tend to deform the drop, and the interfacial tension, that tends to keep the drop cohesive. Therefore, the Weber number assumes equation (7) form, where <u′^{2}> is the mean value of the square of velocity,ρ _{c} is the continuous phase density and d _{max} is the maximum droplet diameter that resists to break up (Hinze, 1955; Lemenand et al., 2003).(7)
Hinze (1955) suggests that, in case of isotropic turbulence, the major contribution for the kinetic energy comes from fluctuations in the region of wave length where the Kolmogorov energy distribution is valid. Thus, <u′^{2}> can be expressed as equation (8), where C _{1} is a constant equivalent to 2, according to Batchelor (1951):(8)
Hinze turbulent inertial break up model for predicting d _{max} is obtained from equations (7) and (8) (Eq. (9)). This model assumes that the dispersed fluid viscosity is negligible (Angelli and Hewitt, 2000; Davies, 1985; Hinze, 1955; Lemenand et al., 2003; Sleicher, 1962):(9)
Hinze (1955) applied this model to experimental tests of Clay (1940a, b), considering d _{95} = d _{max}, where d _{95} is the diameter which 95% of the droplets lay below this value, and found that, for this case, the constant C = 0.725, with standard deviation of 0.315. This model was proposed for diluted dispersed systems. Clay’s experimental tests were performed with dispersions in which the dispersed phase concentration was in the range from 0.9% to 4.8% (Paolinelli and Yao, 2017).
Hinze (1955) suggests that, when dispersed phase viscosity cannot be neglected, a viscosity group can be added to the model (). Sleicher (1962) argued that the viscosity group suggested by Hinze is related to the natural frequency of vibration of droplets, and that this group is small and not relevant. Sleicher (1962), then, suggests that the important viscosity forces are the ones that oppose to the flow fields force, that elongate the droplet. The author indicates that this force per unit area should be proportional to , where U is the mean pipe velocity.
After Hinze’s work, innumerous papers investigating dispersions behavior, mainly of oil in water dispersions, were published. Davies (1985) correlated different types of devices with droplet diameters obtained from turbulent energy dissipation. In this study, a comparison of experimental tests and d _{max} calculated by the Hinze’s model (with typical energy dissipation for each device) showed good predictions. Differently that was defended by Sleicher (1962), Davies (1985) stated that, instead of mean pipe velocity, the turbulent fluctuation velocity (u′) that should be considered on viscosity term. By convenience, the author added the group to the Laplace pressure term on Hinze’s model. The new model proposed by Davies is presented in equation (10):(10)
Both models, previously described, are valid for inertial subrange and, therefore, consider only inertial forces dominating breakup phenomena (λ _{K} < d _{max} < 0.1 pipe diameter). It is expected that continuous phase viscosity do not have any influence. Walstra (1983, 1993) indicate some exceptions where the increase of continuous phase viscosity has influence in droplet size. In this case, some shear stress could be responsible to some of breakup occurring.
As stated before, Grace (1982) indicates that, for very small and very high viscosity ratios, droplet breakup due to shear flow is improbable. However, Stone (1994), affirms that the breakup under these conditions is possible if some vorticity is present in the flow. Moreover, Stone (1994) explained that abrupt changes of the flow type should also be responsible for droplet breakup due to droplet relaxation after deformation. This could occur when the apparatus geometry changes substantially.
Karbstein and Schubert (1995) also stated that is tough to produce emulsions of high viscosity ratios with simple shear flow. They affirm that in turbulent flow, the local fluctuating velocities are responsible droplet deformations, forming bulges and, therefore, causing droplet breakup. These local velocities fluctuations can also lead to droplet oscillations that could end up breaking the droplets.
Cristini et al. (2003) and Vankova et al. (2007) described a drop break up model in viscous turbulent flow (SubKolmogorov regime). This model (Eq. (11)) was estimated by comparing viscous tensions with the droplet capillary pressure (Laplace pressure, Eq. (2)). The boundary condition for the validity of the viscous turbulent break up model is that the generated drop (d _{max v }) must be lower than the Kolmogorov length scale in turbulent flow. In Vankova et al. (2007) study, experimental tests were done to validate the model. In the continuous phase (water), some glycerin was added to enhance viscosity. Different levels of viscosity (18 mPa s max.) were tested and the dispersed phase was a mineral oil (95 mPa s):(11)
The studies of Hinze (1955) and Davies (1985) were the precursors of innumerous other works in this area. Recently, some researchers started to evaluate droplet break up in restrictions (localized pressure drop). In a restriction, the flow passes through an intense turbulent region characterized by a high velocity zone with large velocity gradients, followed by a strong recirculation zone, where the energy dissipation is extremely high (Kwakernaak et al., 2007). The energy dissipation rate (Ė) in a control volume is equal to the permanent difference of pressure (ΔP _{perm}) times the volume flowrate (Q) (van der Zande, 2000). In a monophasic flow, the turbulent energy is dissipated into heat in a dissipation volume (V _{dis}). This volume can be expressed as an area (S) times a dissipation length (L), and the mass of fluid in dissipation zone is obtained multiplying this term by fluid density. Thus, the mean dissipation energy per unit mass can be calculated, as shown in equation (12) (Husveg et al., 2009; van der Zande, 2000):(12)
The energy dissipation calculus considers a dissipation volume or length. van der Zande et al. (1999a) defined that the dissipation zone extends as far as 2.5 times pipe diameter. In posterior work, van der Zande (2000) applying computational simulations observed that this dissipation zone could be higher, varying from 4.1 to 5.4 times pipe diameter. van der Zande et al. (1998c) argued about the possibility to have inhomogeneous energy dissipation rate downstream an orifice. In addition, considering recirculation zones as noneffective region for dissipation, the effective dissipation volume decreases.
In order to be valid, inertial turbulent breakup models based in Hinze theory indicate that the d _{max} obtained must be higher than the Kolmogorov length scale (λ _{K}), defined by equation (13) (where is the kinematic viscosity of continuous phase), and below the higher eddies that are present in the flow (λ _{max}), equivalent to 10% of the pipe diameter (Paolinelli and Yao, 2017):(13)
Percy and Sleicher (1983) proposed a model (Eq. (14)) in which the fluid acceleration and, therefore, a pressure variation in the entrance of an orifice, is responsible for droplets break up. Their experiments were carried out with a low concentration organic dispersed phase in water as continuous phase, passing through different orifices. The term d _{o} is the orifice diameter:(14)
The constant C _{PS} that best fitted Percy and Sleicher experimental points was 3.1. van der Zande and van den Broek (1998b) fitted his experimental points in Percy and Sleicher model and the best adjust was obtained with constant equal to 5.4.
3 Material and methods
The experimental study was conducted in the Phase Separation Laboratory from Compact Separation Nucleus at Universidade Federal de Itajubá (UNIFEI).
This work was focused on the investigation of break up of water droplets dispersed in oil, when the flow passes through a restriction, which produces a localized pressure drop. The goal was to better understand water in oil emulsification when crude oil with residual produced water is submitted to pressure drops in petroleum primary process plants.
3.1 Fluids characterization
The dispersions were obtained by mixing tap water with a commercial mineral oil (USP mineral oil). Water and oil physical properties were evaluated (Tab. 1). A nonionic surfactant with low HLB (sorbitan monooleate) was dosed for a better emulsion homogeneity, for interfacial tension control and to diminish droplet coalescence probability. An experimental planning was performed in two levels for each parameter (viscosity 12 and 26 mPa s, interfacial tension 15 and 25 mN/m, ΔP 3 and 9 bar, and dispersed phase concentrations 5% and 10%). The mass flowrate was kept constant (4 kg/min). Central points were added. The test matrix can be seen in Table 2 in Section 4. Each experimental point was tested in a minimum of three replicates.
Fluids physical properties at 20 °C.
Experimental results.
3.2 Experimental bench
A fluid preconditioning and pumping area was set up. It contained a water tank (200 L), oil tank (500 L) and progressive cavity pumps (all pumps installed with Variable Speed Drives – VSD). Downstream the pumps, Coriolis flowmeters and electrical heaters (aiming temperature, and consequently, viscosity control) were installed. This part of the experimental bench fed the test section where the water was injected in the oil stream in a Tjunction. This commingle point was made with small orifices (diameter of 0.6 mm or 0.5 mm) selected accordingly to water concentration in order to have high velocity in the injection point. The bench was constructed with carbon steel 1/2″ pipes, sch40.
Downstream the Tjunction, a mixing needle valve promotes an initial droplets break up generating the dispersion with the chosen average diameters (volumetric). This dispersion was pumped by a progressive cavity pump, to achieve the test pressure and flowrate, being ready to be sampled as inlet condition or go to the test section. After sampled the inlet condition, the flow was aligned to the test section, the experimental point was adjusted and, after stabilization, the outlet sample was collected. Figure 1 shows a simplified process and instrumentation diagram of the test bench.
Fig. 1. Test bench simplified process and instrumentation diagram. 
3.3 Test section singularity
In the test section, a singularity (mimetic valve) was built and installed (Fig. 2). This dispositive is a channel with squared transversal section, with 5 mm side and 200 mm length. At inlet and outlet, a conical part was made so that pressure drops and recirculation zones could be avoided.
Fig. 2. Gate valve configuration. 
Internally, in the central part, three gates permit that a variety of geometric configurations could be tested. The gates have a micrometric position indicator in order to know their exact position. The squared section was chosen for better discretization in numerical simulation studies.
By adjusting the position of the central gate, the flow through a gate valve could be simulated, with both superior gates kept fully opened. Other geometries can be simulated with this apparatus by changing the position of the other two gates position. Figure 3 shows details of the gate valvelike element, where A is the squared section side, L is the gate opening and 2.5 D _{h} is the dissipation zone, correlated to the hydraulic diameter (D _{h}), defined as , where P is the wetted perimeter. This dissipation zone was applied for energy dissipation calculation.
Fig. 3. Detail of gate valve configuration and dissipation zone. 
The Droplet Size Distributions (DSD) were measured offline using a Malvern Mastersizer 2000. The sample point was located the nearest as possible of the test section in order to avoid coalescence. It was isolated by a ball valve which, when opened, gave full passage, avoiding droplet breakup possibility. Samples were collected and immediately diluted into beakers containing mineral oil, in order to avoid coalescence. Samples were gentle mixed with magnetic mixer and injected in Mastersizer by gravity in a regulated flowrate. All these procedures had the intention to preserve the experimental droplets representative without any possibility of droplet breakup or coalescence.
4 Results and discussion
When the gates were completely opened, the flow regime was set to be laminar. In this condition, the highest Reynolds number calculated for the lowest continuous phase viscosity was less than 1200. This limitation of Reynold number had the objective of minimizing droplets break up due to pipe flow. Under these conditions, the inlet and outlet droplet size distributions were almost the same for the gates completely opened (Fig. 4).
Fig. 4. Inlet and outlet droplet size distributions for the gate fully open (Re_{t} ~ 1100). 
Unlike what happens in pipes, where turbulent flow occurs in Reynolds (Re_{p}) superior to 2100, the critical Reynolds number for orifice (Re_{o}) is smaller. Following Lakshmana Rao et al. (1977) and Janssen et al. (2001), the critical orifice Reynolds number in which the flow becomes turbulent, when passing through an orifice (quadrantedged) with small β = d _{o}/D, is constant and has values of 300–340, increasing with the increase of β, where d _{o} is orifice diameter and D is the pipe diameter.
In a previous study, Johansen (1930) stated that the critical orifice Reynolds number could be even lower. The author observed on photographs that the flow started to become turbulent in Reynolds number above 150, depending on the diameter ratio, β.
Assuming that Re_{p} = βRe_{o}, and considering that the singularity can be compared to an orifice plate, theoretically in all experiment points the flow was laminar until reaching the singularity, where it became turbulent.
Table 2 shows the means of d _{0.9} ~ d _{max} of the experimental points results with its standard deviation. The mean energy dissipation rate (e) and the Kolmogorov (λ _{K}) length scale were calculated using equations (12) and (13), respectively. Besides, the condition of each test point was stated. Table shows the viscosity (μ), interfacial tension (σ), pressure drop across the apparatus (ΔP) and Concentration of water droplets in oil (Conc).
The experimental results are discussed separately. First, low concentration points (5% of dispersed phase in oily continuous phase) together with the central point (concentration of 7.5%) are analyzed. It was considered that d _{0.9} ~ d _{max} and that the turbulent energy dissipation (Eq. (12)) was calculated with the dissipation volume equivalent to the volume above of the gate plus the volume corresponding to the dissipation zone (Fig. 3), following van der Zande et al. (1999a) and Fossen and Schumann (2016).
The abovementioned experimental points means (with three replicates at least), separated by viscosity and interfacial tension levels, are shown in Figure 5. According to this figure, the continuous phase physical properties have no influence on the results. The Kolmogorov length scales curves for low and high viscosities are also presented in this figure.
Fig. 5. Approximated maximum diameter (d _{0.9} ≅ d _{max}) for experimental points at low concentrations (5% and 7.5%). 
Assuming that the larger eddy scale (λ _{max}) is equivalent to 0.1 D (Paolinelli and Yao, 2017), it is possible to say that λ _{K} < d _{max} < λ _{max} and, therefore, inertial break up models could be applied to the experimental results.
As can be seen in Figure 6, Hinze’s model underestimates the experimental d _{max}. This was also observed in other studies with higher continuous phase viscosity and with dispersed phase in higher concentration (Janssen et al., 2001; and Paolinelli and Yao, 2017).
Fig. 6. A comparison between the experimental points for the approximated maximum diameter (d _{0.9} ≅ d _{max}) at low concentrations (5% and 7.5%) and different droplet break up models. 
Despite the fact that Hinze’s model was not adequate to predict droplet size downstream the gate, it is important to notice that, according to the adjusted curve, the experimental d _{max} ∝ which is similar to Hinze’s model where d _{max} ∝ (Eq. (9)).
Applying equation (9) to the experimental points at low concentrations (5% and 7.5%), a constant of 3.2 is found, which is higher than Hinze’s original constant of 0.725. The predicted curve using this new constant is also shown in Figure 6.
Considering the orifice droplet break up model proposed by Percy and Sleicher (1983) and applying the proportionality constant suggested by the authors of 3.1, this model also underestimates d _{max} and is only slightly better than Hinze’s model. The calculation of ΔP _{max} (therm of Percy and Sleicher model) was carried out following Paolinelli and Yao (2017), in which ΔP _{max} ≈ ΔP _{perm}/(1 − β ^{2}), where β is the ratio between hydraulic diameter of the gate opening and the hydraulic diameter of the channel. van der Zande and van den Broek (1998b) better fit their experimental points with C = 5.4. By applying this constant value, fairly better predictions of d _{max} were obtained.
Figure 7 shows the results for the high dispersed phase concentration (10%) adding the central point (7.5%) for comparison. A higher dispersion of d _{max} is evident. As can be seen, points with lower viscosity presented higher maximum droplet size. The lower viscosity allied to higher concentration of droplets can make the dispersion more prone to coalescence effects (higher dispersed droplets interaction).
Fig. 7. Approximated maximum diameter (d _{0.9} ≅ d _{max}) for experimental points at high concentration (10%) and the central point (7.5%). 
The increase of d _{max} due to higher dispersed phase concentration was also reported by Paolinelli and Yao (2017). They suggested that this increase in d _{max} was due to two phenomena. First, there is a turbulent tension effectiveness suppression, also called turbulence attenuation, due to an increase of emulsion apparent viscosity caused by the presence of higher content of dispersed droplet. The second phenomenon is due to coalescence that, at higher content, increases the probability of droplet collision (Walstra and Smulders, 1998).
It can be noted that at lower concentration the power law adjusted reasonably well the result for droplet break up at the central point (Fig. 6). At high concentration, it was not possible to fit a single power law equation for the experimental points (Fig. 7). This suggests that the effects of turbulence suppression and/or coalescence just occur in concentrations above 7.5% of water dispersed in oil.
Figure 8 shows a comparison between the experimental points for the approximated maximum diameter (d _{0.9} ≅ d _{max}) at high concentration (10%) and different droplet break up models. In accordance with what was done previously, with the low concentration experimental points, a new constant for Hinze’s model was tested. The best value for the constant was C = 4.1. The curve of prediction considering the new constant to Hinze model is also shown in Figure 8.
Fig. 8. A comparison between the experimental points for the approximated maximum diameter (d _{0.9} ≅ d _{max}) at high concentration (10%) and different droplet break up models. 
Although all values of d _{max} are above Kolmogorov length scale (Figs. 5 and 7), the proximity to these values justified an evaluation of turbulent viscous droplet break up model. Figures 9 and 10 compare the experimental results of d _{max} with curves of prediction by the model proposed by Vankova et al. (2007) (Eq. (11)) for 5% and 10% concentrations, respectively. Curves were calculated considering the lowest and highest levels for viscosity, interfacial tension and density.
Fig. 9. Comparison of experimental d _{max} and curves predicted by the model proposed by Vankova et al. (2007) for low concentration (5%). 
Fig. 10. Comparison of experimental d _{max} and curves predicted by the model proposed by Vankova et al. (2007) for high concentration (10%). 
From Figures 9 and 10, it can be observed that for low concentration the model proposed by Vankova et al. (2007) gives better results than for high concentration.
5 Conclusion
The experimental maximum droplet diameter decreased with increased turbulent energy dissipation rate according to a power law with a power of −0.42 for lower concentrations (less than 7.5%). This power is similar to that of Hinze’s model (−0.4).
Increases in dispersed phase concentration increases the maximum diameters. Turbulence suppression, due to an increase of emulsion apparent viscosity, and/or coalescence, due to higher droplet population, can possibly explain it. Although lots of careful procedures adopted on sampling and measuring, some coalescence could have happened.
Hinze’s inertial model was not capable to describe the experimental results. An increase in his proportionality constant tends to improve prediction. Percy and Sleicher model also was also not able to describe the results with the constant suggested by the authors. Applying the van der Zande and van den Broek (1998b) constant in Percy and Sleicher model gave slightly better results.
Although all experimental d _{max} were above Kolmogorov length scale, which would indicate that inertial forces are governing the phenomenon, turbulent viscous model described fairly well the experimental points. This may indicate that viscous forces may have some influence on the break up process.
New experiments will be performed to better understand which of the models (corrected turbulent inertial or turbulent viscous) is the most adequate to describe droplets break up phenomenon for a low viscous liquid dispersed in a more viscous continuous phase.
Acknowledgments
The authors would like to thank PETROBRAS and UNIFEI for their support for this study.
References
 Angelli P., Hewitt G.F. (2000) Drop size distributions in horizontal oilwater dispersed flows, Chem. Eng. Sci. 55, 16, 3133–3143. https://doi.org/10.1016/S00092509(99)005850. [Google Scholar]
 Arnold K., Stewart M. (1999) Arnold K., Stewart M., Surface production operations: design of oilhandling systems and facilities, Second Edition, Chapter 1, 2, Gulf Professional Publishing, ISBN 9780884158219. https://doi.org/10.1016/B9780884158219/500026. http://www.sciencedirect.com/science/article/pii/B9780884158219500026. [Google Scholar]
 Batchelor G.K. (1951) Pressure fluctuations in isotropic turbulence, Math. Proc. Camb. Philos. Soc. 47, 2, 359–374. https://doi.org/10.1017/S0305004100026712. [NASA ADS] [CrossRef] [Google Scholar]
 Clay P.H. (1940a) The mechanism of emulsion formation in turbulent flow. I. Experimental Part, Proc. Roy. Acad. Sci. (Amsterdam) 43, 852. [Google Scholar]
 Clay P.H. (1940b) The mechanism of emulsion formation in turbulent flow. II. Theoretical part and discussion, Proc. Roy. Acad. Sci. (Amsterdam) 43, 879. [Google Scholar]
 Cristini V., Blawzdziewicz J., Loewenberg M., Collins L. (2003) Breakup in stochastic stokes flows: subKolmogorov drops in isotropic turbulence, J. Fluid Mech. 492, 231–250. https://doi.org/10.1017/S0022112003005561. [Google Scholar]
 Davies J.T. (1985) Drop sizes of emulsions related to turbulent energy dissipation rates, Chem. Eng. Sci. 40, 839–842. https://doi.org/10.1016/00092509(85)850363. [Google Scholar]
 Fossen M., Schumann H. (2016) Experimental study of the relative effect of pressure drop and flow rate on the droplet size downstream a pipe restriction, J. Dispers. Sci. Technol. 38, 6, 826–831. https://doi.org/10.1080/01932691.2016.1207184. [Google Scholar]
 Galinat S., Masbernat O., Guiraud P., Dalmazzone C., Noïk C. (2005) Drop breakup in turbulent pipe flow downstream of a restriction, Chem. Eng. Sci. 60, 23, 6511–6528. https://doi.org/10.1016/j.ces.2005.05.012. [Google Scholar]
 Grace H.P. (1982) Dispersion phenomena in high viscosity immiscible fluid systems and application of static mixers as dispersion devices in such systems, Chem. Eng. Commun. 14, 225–277. https://doi.org/10.1080/00986448208911047. [Google Scholar]
 Hinze J.O. (1955) Fundamentals of the hydrodynamic mechanism of splitting in dispersion process, AIChE J. 1, 3, 289–295. https://doi.org/10.1002/aic.690010303. [Google Scholar]
 Husveg T., Bilstad T., Guinée P.G.A., Jernsletten J., Knudsen B., Nordbo H.T. (2009) A cyclonebased low shear valve for enhanced oilwater separation, Offshore Technology Conference, OTC 2009, https://doi.org/10.4043/20029MS. [Google Scholar]
 Janssen P.H., Noik C., Dalmazzone C. (2001) Emulsion formation in a model chokevalve, Society of Petroleum Engineers, SPE 71473. https://doi.org/10.2118/71473MS [Google Scholar]
 Johansen F.C. (1930) Flow through pipe orifices at low Reynolds numbers, Proc. R. Soc. Lond. A 126, 231–245. https://doi.org/10.1098/rspa.1930.0004. Published 1 January 1930. [CrossRef] [Google Scholar]
 Karbstein H., Schubert H. (1995) Developments in the continuous mechanical production of oilinwater macroemulsions, Chem. Eng. Process. Process. Intensification 34, 3, 205–211. https://doi.org/10.1016/02552701(94)040052. [CrossRef] [Google Scholar]
 Kubie J., Gardner G.C. (1977) Drop sizes and drop dispersion in straight horizontal tubes and in helical coils, Chem. Eng. Sci. 32, 195–202. https://doi.org/10.1016/00092509(77)80105X. [Google Scholar]
 Kwakernaak P.J., van den Broek W.M.G.T., Currie P.K. (2007) Reduction of oil droplet breakup in a choke, Society of Petroleum Engineers, SPE 106693. https://doi.org/10.2118/106693MS. [Google Scholar]
 Lakshmana Rao N.S., Sridharam K., Alvi S.H. (1977) Critical Reynolds number for orifice and nozzle flows in pipes, J. Hydraulic Res. 15, 2. https://doi.org/10.1080/00221687709499654. [CrossRef] [Google Scholar]
 Lemenand T., Della Valle D., Zellouf Y., Peerhossaini H. (2003) Droplets formation in turbulent mixing of two immiscible fluids in a new type of static mixer, Int. J. Multiph. Flow 29, 5, 813–840. https://doi.org/10.1016/S03019322(03)000326. [CrossRef] [Google Scholar]
 Lockhart F.J., Chilingarian G.V., Kumar S. (1987) Chapter 3 Separation of oil and Gas, in: G.V. Chilingarian, J.O. Robertson, S. Kumar (eds), Developments in Petroleum Science, vol. 19, Elsevier, pp. 59–108, ISSN 03767361, ISBN 9780444416254, https://doi.org/10.1016/S03767361(08)705326. http://www.sciencedirect.com/science/article/pii/S0376736108705326. [CrossRef] [Google Scholar]
 Manning F.S., Thompson R.E. (1995) Chapter 1, 2 and 6 – Oil field processing, Vol. 2, Pennwell Books, Crude Oil, ISBN 9780878143542. [Google Scholar]
 Mitre J.F., Lage P.L.C., Souza M.A., Silva E., Barca L.F., Moraes A.O.S., Coutinho R.C.C., Fonseca E.F. (2014) Droplet breakage and coalescence models for the flow of waterinoil emulsions through a valvelike element, Chem. Eng. Res. Des. J. 92, 2493–2508. [CrossRef] [Google Scholar]
 Paolinelli L.D., Yao A.R.J. (2017) Characterization of droplets sizes in large scale oilwater flow downstream from a globe valve, Int. J. Multiph. Flow 99, 132–150. https://doi.org/10.1016/j.ijmultiphaseflow.2017.09.014. [CrossRef] [Google Scholar]
 Percy J.S., Sleicher C.A. (1983) Drop breakup in the flow of immiscible liquids through an orifice in a pipe, AIChE J. 29, 1, 161–164. https://doi.org/10.1002/aic.690290125. [Google Scholar]
 Rueger P.E., Calabrese R.V. (2013a) Dispersion of water into oil in a RotorStator mixer. Part 1: drop breakup in dilute systems, Chem. Eng. Res. Des. 91, 2122–2133. https://doi.org/10.1016/j.cherd.2013.05.018. [Google Scholar]
 Rueger P.E., Calabrese R.V. (2013b) Dispersion of water into oil in a Rotorstator Mixer. Part 2: Effect of phase fraction, Chem. Eng. Res. Des. 91, 2134–2141. https://doi.org/10.1016/j.cherd.2013.06.010. [Google Scholar]
 Simmons M.J.H., Azzopardi B.J. (2001) Drop size distributions in dispersed liquidliquid pipe flow, Int. J. Multiph. Flow 27, 5, 843–859. https://doi.org/10.1016/S03019322(00)000550. [CrossRef] [Google Scholar]
 Sleicher C.A. (1962) Maximum stable drop size in turbulent flow, AIChE J. 8, 4, 471–477. https://doi.org/10.1002/aic.690080410. [Google Scholar]
 Stone H.A. (1994) Dynamics of drop deformation and breakup in viscous fluids, Annu. Rev. Fluid Mech. 26, 65–102. https://doi.org/10.1146/annurev.fl.26.010194.000433. [Google Scholar]
 Taylor G.I. (1932) The viscosity of a fluid containing small drops of another fluid, Proc. R. Soc. Lond. A 138, 41–48. Published 1 October 1932 https://doi.org/10.1098/rspa.1932.0169. [CrossRef] [Google Scholar]
 Taylor G.I. (1934) The formation of emulsions in definable fields of flow, Proc. R. Soc. Lond. A 146, 501–523. Published 1 October 1934, https://doi.org/10.1098/rspa.1934.0169. [CrossRef] [Google Scholar]
 van der Zande M.J. (2000) Droplet breakup in turbulent oilinwater flow through a restriction, PhD Thesis, ISBN 9090138706. [Google Scholar]
 van der Zande M.J., Muntinga J.H., van den Broek W.M.G.T. (1998a) Emulsification of production fluids in the choke valve, Society of Petroleum Engineers, SPE 49173. https://doi.org/10.2118/49173MS. [Google Scholar]
 van der Zande M.J., van den Broek W.M.G.T. (1998b) Breakup of oil droplets in the production system, Proc. of ASME Energy Sources Technology Conference and Exhibition, Houston, ETCE984744. [Google Scholar]
 van der Zande M.J., van den Broek W.M.G.T. (1998c) The effect of tubing and choke valve on oildroplet breakup, Proceedings of the 1st North American Conference on Multiphase Technology, June 10–11, Banff, Canada, pp. 89–100. [Google Scholar]
 van der Zande M.J., van Heuven K.R., Muntinga J.H., van den Broek W.M.G.T. (1999a) Effect of flow through a choke valve on emulsion stability, Society of Petroleum Engineers, SPE 56640. https://doi.org/10.2118/56640MS. [Google Scholar]
 van der Zande M.J., Muntinga J.H., van den Broek W.M.G.T. (1999b) The effects of production rate and choke size on emulsion stability, Paper EXPL6MZ, presented at the 3rd International Seminar in Practices of Oil and Gas Exploitation, INGEPET 99, October 27–29, Lima, Peru. [Google Scholar]
 Vankova N., Tcholakova S., Denkov N.D., Ivanov I.B., Vulchev V.D., Danner T. (2007) Emulsification in Turbulent flow 1. Mean and maximum drop diameters in inertial and viscous regimes, J. Colloid Interface Sci. 312, 363–380. https://doi.org/10.1016/j.jcis.2007.03.059. [Google Scholar]
 Walstra P. (1983) Encyclopedia of emulsion technology, Vol. 1, Chapter 2, in: Becher P. (ed), Marcel Dekker, New York. ISBN13: 9780824718763, ISBN10: 0824718763. [Google Scholar]
 Walstra P. (1993) Principles of emulsion formation, Chem. Eng. Sci. 48, 2, 333–349. https://doi.org/10.1016/00092509(93)80021H. [Google Scholar]
 Walstra P., Smulders P.E.A. (1998) Modern aspects of emulsion science, Chapter 2, 1st edition (October 9, 1998), in: Binks B.P. (ed), Royal Society of Chemistry, UK, pp. 56–99. http://dx.doi.org/10.1039/978184755147400056. [Google Scholar]
All Tables
All Figures
Fig. 1. Test bench simplified process and instrumentation diagram. 

In the text 
Fig. 2. Gate valve configuration. 

In the text 
Fig. 3. Detail of gate valve configuration and dissipation zone. 

In the text 
Fig. 4. Inlet and outlet droplet size distributions for the gate fully open (Re_{t} ~ 1100). 

In the text 
Fig. 5. Approximated maximum diameter (d _{0.9} ≅ d _{max}) for experimental points at low concentrations (5% and 7.5%). 

In the text 
Fig. 6. A comparison between the experimental points for the approximated maximum diameter (d _{0.9} ≅ d _{max}) at low concentrations (5% and 7.5%) and different droplet break up models. 

In the text 
Fig. 7. Approximated maximum diameter (d _{0.9} ≅ d _{max}) for experimental points at high concentration (10%) and the central point (7.5%). 

In the text 
Fig. 8. A comparison between the experimental points for the approximated maximum diameter (d _{0.9} ≅ d _{max}) at high concentration (10%) and different droplet break up models. 

In the text 
Fig. 9. Comparison of experimental d _{max} and curves predicted by the model proposed by Vankova et al. (2007) for low concentration (5%). 

In the text 
Fig. 10. Comparison of experimental d _{max} and curves predicted by the model proposed by Vankova et al. (2007) for high concentration (10%). 

In the text 