Next Article in Journal
Ultracapacitors for Port Crane Applications: Sizing and Techno-Economic Analysis
Next Article in Special Issue
Efficient History Matching of Thermally Induced Fractures Using Coupled Geomechanics and Reservoir Simulation
Previous Article in Journal
Numerical and Experimental Analysis of the Thermal Performances of SiC/Water and Al2O3/Water Nanofluid Inside a Circular Tube with Constant-Increased-PR Twisted Tape
Previous Article in Special Issue
Heavy Oil Laminar Flow in Corrugated Ducts: A Numerical Study Using the Galerkin-Based Integral Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Isothermal Treatment of Oily Waters Using Ceramic Membrane: A Numerical Investigation

by
Hortência L. F. Magalhães
1,*,
Gicelia Moreira
1,
Ricardo S. Gomez
2,
Túlio R. N. Porto
2,
Balbina R. B. Correia
2,
Anderson M. V. Silva
2,
Severino R. Farias Neto
1 and
Antonio G. B. Lima
2
1
Department of Chemical Engineering, Federal University of Campina Grande, Campina Grande 58429-900, Brazil
2
Department of Mechanical Engineering, Federal University of Campina Grande, Campina Grande 58429-900, Brazil
*
Author to whom correspondence should be addressed.
Energies 2020, 13(8), 2092; https://doi.org/10.3390/en13082092
Submission received: 21 February 2020 / Revised: 17 April 2020 / Accepted: 20 April 2020 / Published: 22 April 2020

Abstract

:
Currently, the oil industry deals with the challenge of produced-water proper disposal, and the membrane-separation technology appears as an important tool on the treatment of these waters. In this sense, this work developed a mathematical model for simulating the oil/water separation by a ceramic membrane. The aim was to investigate the thermal aspects of the separation process via computational fluid dynamic, using the Ansys CFX® 15 software (15, Ansys, Inc., Canonsburg, PA, USA). Oil concentration, pressure, and velocity distributions, as well as permeation velocity, are presented and analyzed. It was verified that the mathematical model was capable of accurately representing the studied phenomena and that temperature strongly influences the flow behavior.

Graphical Abstract

1. Introduction

In the oil fields, significant amounts of oil-contaminated water are generated which are commonly known as produced water. These waters are provided from the reservoir and water injection, during the secondary recovery process. The oily water is an aqueous solution that is rich in salts, which may contain dispersed or emulsified oil droplets.
Operations at refineries generate a large amount of difficult-to-treat wastewater, rich in organic pollutants. Effluent treatment plants are generally used to treat these waters. The purpose is to change the effluent physical and/or chemical compositions, to enable the water’s proper disposal into the environment, according to the actual legislation. During the treatment, dispersed or dissolved hydrocarbons, suspended solids, toxic organic compounds (BTEX—benzene, toluene, ethylbenzene, and xylenes), phenol, and inorganic compounds—such as cyanides, sulfides, metals, phosphates, and nitrogen compounds—are removed [1,2].
Currently, many processes are used to remove contaminants from produced water, such as electrochemical treatment, membrane filtration, biological treatment, flotation, and photocatalysis, among others.
The electrochemical treatment uses an electrical potential to oxidize or reduce substances present in contaminated water; this treatment is considered to be a highly efficient, versatile method with easy automation and control [3]. However, the electrodes are consumed over time and need to be replaced regularly.
Photocatalysis, in turn, is a process commonly used to decompose compounds by the action of the free radicals, such as hydroxyl, obtained from a catalyst that absorbs photons of ultraviolet light. This method has strong oxidative power and the ability to fully mineralize pollutants; however, it requires high energy consumption [4].
Biological treatment is another alternative; it is capable of degrading organic matter by the action of microorganisms in an aerobic (with the presence of oxygen) or anaerobic (without the presence of oxygen) way. It demands low energy cost, but requires an extensive area for implantation of the treatment system.
Flotation consists of a separation process that uses tiny air bubbles added to oily water, to carry the oil from the mixture to the surface, using the principle of the difference in density between the phases. Compared to the biological treatment method, it requires a smaller installation area. It also presents considerable contaminant removal [5].
The produced water treatment using membrane has been applied due to characteristics such as permeate uniformity, long working life, and operational simplicity. Acting as a semipermeable selective barrier, by pressure and concentration gradients, the membrane, usually constituted by ceramic or polymeric materials, promotes the physical separation (ultra and microfiltration) between water and oil [2,6,7,8,9]. The disadvantages of the method are the inefficiency in the treatment of ionic compounds and the need to operate with low contaminant concentrations. However, it is considered an efficient and economical treatment when compared to other conventional techniques, such as decantation, centrifugal separation, and flotation.
Studies were carried out by several researchers, aimed at improving membrane-separation technology [10]. Zsirai et al. [11] and Ding et al. [12] used different materials for the membrane manufacturing, making use of structural modifications, aiming at improving membrane characteristics, and consequently increasing the separation efficiency.
Within theoretical studies, Frederic et al. [13] studying the modeling of a multichannel ceramic membrane applied for produced water filtration, based on Darcy’s equation and series resistance model. The authors observed that the existence of channels can considerably affect the permeate flux distribution. Furthermore, it was verified that the initial operating conditions can contribute to attenuate incrustation formation.
Mota et al. [14] have evaluated the produced water treatment by ceramic membranes under the regime of turbulent flow induced at a tangential inlet, using a multiphase mathematical model, and RNG k-ɛ turbulence model. From the obtained results, it was verified that the geometric aspect of the feed tube (circular or rectangular) with the same transverse area does not affect the permeate flow.
Motin et al. [15] numerically investigated the performance of a rotating tubular membrane operating in crossflow for the treatment of produced waters, using a discrete phase model. The idea is to estimate the trajectory of the oil droplet. The authors have identified that the increase in shear stress, caused by angular velocity and cross-flow, decreases the oil concentration near the membrane surface, increasing the separation efficiency.
Several theoretical works are found in the literature, using CFD (computational fluid dynamics) to evaluate the oily-water treatment process, using membranes, as in the works by Serra et al. [16] and Serra et al. [17], using circular membranes; Geraldes et al. [18] and Darcovich et al. [19], using flat membranes; and Souza et al. [20], Cunha et al. [21], and Alves et al. [22], using modules equipped with a membrane. Most of the theoretical works adopt a two-dimensional or three-dimensional domain, assuming a fully developed velocity profile and laminar flow, thus solving the conservation and mass transport equations. Among the works that use this methodology are Damak et al. [23,24,25], Pak et al. [26], and Minnikanti et al. [27]. In general, the main purposes of these works were to evaluate the thickness of the concentration polarization and determine the pressure drop in the membrane.
Despite the importance, few works include the effect of the temperature in the produced-water separation process by membrane filtration, especially using CFD analysis and the proposed geometry. Thus, the present work seeks to evaluate the influence of the thermal effects on the water/oil separation process via CFD, as the thermal effects have not been well explored in the literature, presenting a new approach to treat oily water and providing data that enable the optimization of the currently used processes. Herein, it is realized a comparison between the two- and three-dimensional approaches of the process, aiming to contribute to the understanding of temperature and geometry influence in the concentration distribution and permeation velocity.

2. Methodology

2.1. The Physical Problem and Geometries

2.1.1. 2D Analysis

For the studying of the oil–water separation process, in a two-dimensional approach, we used a tubular membrane with a tube diameter of 0.03 m and a length of 3 m [23], subject to oily water flow (Figure 1). Due to the membrane angular symmetry, only one longitudinal section was used (YZ plane).
For the 2D numerical simulation, a mesh was generated, using ICEM CFD 15 software (15, Ansys, Inc., Canonsburg, PA, USA), having a total of 205,056 elements and 180,000 nodes, with inflations near the membrane surface, as illustrated in Figure 2.

2.1.2. Three-Dimensional Analysis

In this case, we used the water/oil separation module known as microfiltration, which consists of a tube shell, with diameter De, enclosing a tubular membrane with diameter di. The module has a rectangular input and output (Figure 3). The main dimensions of the module are shown in Table 1.
For the 3D mesh generation, we used the blocking strategy, with the creation of hexahedral elements within the domain, aiming to facilitate the refinement. The resulting final mesh has a total of 217,352 elements and 180,780 nodes, as illustrated in Figure 4.

2.2. Mathematical Modeling

The water/oil separation process through the tubular membrane was described by using mathematical modeling based on the conservation equations (mass, momentum, and energy) and mass transport. The aim was to evaluate the thermal effect on the fluid-flow behavior.

2.2.1. The Model

For the study of the microfiltration process using a 2D tubular membrane and a shell and tube membrane-separation module, the following considerations were used for the mathematical modeling:
(a)
Steady-state;
(b)
Laminar and incompressible flow;
(c)
Fluids density and viscosity are dependent on the temperature;
(d)
The local permeation rate was established by the series resistances theory;
(e)
Concentration layer is considered homogeneous, and the Carman–Kozeny equation is valid;
(f)
The axial velocity profile is fully developed in the membrane inlet;
(g)
The obstruction of the porous medium by the oil is neglected;
(h)
The mass diffusivity of the dispersed phase (oil) was assumed to be constant;
(i)
The gravitational effect was neglected;
(j)
The fouling mechanisms were neglected;
(k)
The temperature range from 35 to 95 °C was used for the analyzes.
From the considerations made, the following equations can be written:
(a)
Mass conservation equation
( ρ U ) = 0
where ρ is the density, and U is the velocity vector.
(b)
Linear momentum equation
( ρ U U ) ( μ U ) = p + [ μ ( U ) T ]
where μ is the solution viscosity, and p is the pressure.
(c)
Energy conservation equation
( ρ U H ) ( Γ e T ) = 0
where H is the enthalpy, and Γ e is the effective thermal diffusivity.
(d)
Mass transport equation
U C = D A B 2 C
D A B is the mass diffusivity (constant for each Schmidt number) (Equation (5)), and C is the oil concentration.
D A B = μ S c ρ

2.2.2. Boundary Conditions

(a) Inlet condition (2D and 3D analyses)
At the equipment inlet is considered fully developed flow, at an initial temperature (T0), with zero radial velocity and axial velocity expressed by Equation (6):
U = Re   μ ρ ( D e d i )
where Re is the axial Reynolds number.
Oil initial concentration (Equation (7)) and the inlet fluid temperature (Equation (8)) are expressed as follows:
C = C 0
T = T 0
(b) Output condition (2D and 3D analyses)
At the system outlet was assumed a pressure condition equal to atmospheric (1 atm) and parabolic condition for concentration (Equations (9) and (10)), as follows:
C z = 0
C y = 0
(c) Symmetry condition (only for the 2D analysis)
The following symmetry conditions were assumed on the tube axis:
U z y = 0
C y = 0
T y = 0
The following conditions were assumed for the transverse planes:
U z x = 0
C x = 0
T x = 0
(d) Porous wall condition (2D and 3D analyses)
A no-slip condition (zero axial velocity) (Equation (17)) was assumed in the membrane surface, as well as an equation responsible for convective heat transport (Equation (18)).
U z = 0
q = h ¯ c ( T s T )
where T is the fluid temperature (air), T s is the surface membrane temperature, and h ¯ c is the convective heat transfer coefficient (W/m2 K), calculated according to Equation (19).
h ¯ c = N u k L c
where k is the thermal conductivity, Nu is the Nusselt number, and Lc is the characteristic length.
Characteristic length, for flow over horizontal cylinders, is defined as the tube diameter (Equation (20)):
L c = d i
To calculate the Nusselt number, for natural convection on a horizontal cylinder, Equation (21) was used [28].
N u = { 0.60 + 0.387 R a 1 / 6 [ 1 + ( 0.559 / Pr ) 9 / 16 ] 8 / 27 } 2
where Pr is Prandtl number, Ra is the Rayleigh number (Equation (22)), and Gr is the Grashof number (Equation (23)).
R a = G r Pr
G r = ρ 2 g β d Δ T L C 3 μ 2
where μ is the dynamic viscosity, and β is the thermal expansion coefficient. Through the membrane, there is a radial velocity, U y , expressed by the permeation velocity, U w , as follows:
U y = U w = Δ P μ ( R m + R p )
where Rp is the resistance due to concentration polarization, and Rm is the membrane resistance.
Mass transport (Equation (25)) [27] was added, to be modeled as a source term. In this case, the following equation was used:
U w C R r = D A B C y
where Rr is the oil intrinsic retention by the membrane, considered constant and equal to one. This consideration means 100% of oil rejection by the membrane.
In Equation (24), the transmembrane pressure, ΔP, can be express as follows:
Δ P = P ¯ p P e x
where P ¯ p is the average permeate pressure, and Pex is the external pressure to the membrane.
The terms Rm and Rp in Equation (24) are expressed by Equations (27) and (28), respectively.
R m = e K
R p = r p δ p
where e is the membrane thickness, K is the porous media permeability, δ p is the concentration layer thickness, and r p is the resistance (when the concentration layer is considered homogeneous). The parameter r p is calculated by using the Carman–Kozeny equation, as follows:
r p = 180 ( 1 ε p ) 2 d p 2 ε p 3
where ε p is the porosity of the concentration polarization, and dp is the mean particle diameter (oil).
Determination of the local variation of the boundary layer thickness by polarization, δ p , can be performed by using Equation (30) [24], valid for Re = 300~1000.
δ p d i = 2 ( z d i ) 0.33 ( Re S c ) 0.33 Re w 0.3 { [ 1 0.4377 ( S c 0.0018 Re w 0.1551 ) ] }
where z is the axial coordinate along the membrane, Sc is the Schmidt number, di is the internal diameter, Re is the axial Reynolds number, and Re w is the permeation Reynolds number. According to the authors [24], this correlation is especially useful, as it allows for the evaluation of the boundary-layer thickness by polarization with a 10% average error, as a function of the operational parameters.

2.3. Physical Properties of the Fluid

Thermal conductivity, specific heat, density, and viscosity of the water were calculated by using expressions as a function of the temperature, as reported in the literature.
(a) Thermal conductivity (k)
The thermal conductivity of water as a function of temperature is given by Equation (31) [29]:
k = 0.6065 [ 1.48445 + 4.12292 ( T 298.15 ) 1.63866 ( T 298.15 ) 2 ]
where T is the temperature in Kelvin.
(b) Specific heat (cp)
Specific heat of the water was determined by Equation (32), as follows [30]:
c p = ( A + BT + CT 2 + DT 3 ) / M w
where T is the temperature in Kelvin, Mw is the molecular mass of water in kg/kmol, and A, B, C, and D are the regression coefficients as reported in Table 2.
(c) Density (ρ)
Water density as a function of the temperature was determined by using the following equation [25]:
ρ = A B ( 1 T / T c ) N
where Tc is the critical temperature, and the parameters A, B, and N are the regression coefficients, as reported in Table 3.
(d) Dynamic viscosity (μ)
Equation (34) was used to determine the dynamic viscosity of water [30].
μ = 10 ( A + B T + C T + D T 2 )
where T is the temperature in Kelvin, and the terms A, B, C, and D are the regression coefficients, as shown in Table 4.
The model validation was realized in a previous work [9], by comparing the numerical results obtained in the simulations with the data available in some works reported in the literature [23,24,25,26,31].

2.4. Studied Cases

Simulations were performed, by using different fluid inlet temperatures, and the other parameters were kept constant; these values were adopted based on the work developed by Magalhães et al. [9], Damak et al. [23], and Magalhães et al. [32]. Table 5 presents information about the parameters used in the simulations, and Table 6 shows all the cases studied.

3. Results and Discussion

3.1. Concentration Profile

Figure 5 shows the concentration profiles along the z coordinate (for y = 0.0225 m), for Cases 2 (two-dimensional) and 9 (three-dimensional), both with T0 = 55 °C. From the analysis of this figure, it is observed that the oil concentration increases along with the axial position, due to the particle accumulation near the membrane surface. Furthermore, it is observed that the geometry of the device considered in the analysis slightly affect the concentration profile, since the difference between the concentrations of the 3D and 2D analyses, along with the radial position, is only 3 × 10−3 kg/m3. When comparing the two-dimensional result with that obtained in a three-dimensional approach, where the permeation occurs on the outer tube surface of the shell and tube membrane separation module, it is verified that the 3D geometry presents higher concentrations, with variations close to the inlet (z = 0 m) and outlet (z = 3 m), presenting inlet and outlet solute concentration equal to 1.000 and 1.0089 kg/m3, respectively; this behavior is induced justly by feed inlet and outlet locations of the concentrate in the tridimensional geometry, which is not observed for 2D geometry.
Figure 6 shows the oil-concentration behavior as a function of the longitudinal position over four angles (0°, 90°, 180°, and 270°) on the membrane surface of the module, Cases 5 (T0 = 35 °C) and 8 (T0 = 95 °C). It illustrates that the oil-concentration distribution is affected by the variation temperature, especially at the outlet region (z = 3 m). In Figure 6a, it can be seen that the concentration profile, at 0° and 180° angular positions, have higher oil concentrations, as compared with the results observed in Figure 6b, for the angular positions 90° and 270°, due to the fluid flow behavior around the membrane. Higher concentrations were obtained, as expected, at the impact region between the fluid and membrane, provided by the feed inlet perpendicular position.

3.2. Permeation Velocity

Figure 7 shows the permeation velocity profile as a function of the axial position for Cases 1–4 (two-dimensional cases) and Cases 5–8 (three-dimensional cases). It is verified that, as the temperature is increased, a reduction in the permeation velocity is observed; this behavior is associated with the reduction of the pressure gradient within the equipment, due to the inlet velocity reduction, caused by the variation in the fluid properties. For example, by raising the temperature from 35 to 95 °C, the relationship between viscosity and density (μ/ρ) is reduced by 4.14 × 10−7 m2/s. Furthermore, it was verified that the geometry influences the permeation velocity; the characteristic curves for the three-dimensional cases present lower velocities than the two-dimensional curves. This behavior is associated with the enlargement of the runoff area under the same operating conditions imposed in the 2D analysis. Reinforcing this analysis, it is possible to visualize the temperature influence on permeate mass flow (Table 7).

3.3. Velocity Field

Figure 8 shows the velocity field on the YZ longitudinal plane of the separation module for the temperatures of 35, 55, 75, and 95 °C (Cases 5–8). By analyzing this figure, it is observed that the temperature influences the velocity fields: Higher velocity gradients are presented for 35 °C (Figure 8a). This behavior is expected, because, as the temperature increases from 35 to 95 °C the fluid viscosity and density are reduced; consequently, the velocity inside the equipment is also decreased. It is possible to analyze this behavior through the maximum velocities reached in each analysis (U35°C = 0.024 m/s, U55°C = 0.016 m/s, U75°C = 0.012 m/s, and U95°C= 0.010 m/s). There are larger variations in the fluid velocity at the inlet and outlet regions. Furthermore, we can see that velocity values are small, with laminar flow behavior.

3.4. Temperature Field

Figure 9, Figure 10, Figure 11 and Figure 12 show the temperature fields within the separation module for feed temperatures 35, 55, 75, and 95 °C (Cases 5–8). A temperature gradient is observed, at annular space between the shell module and the ceramic membrane, caused by the heat transfer between the heated membrane surface and atmospheric air, by natural convection, thereby cooling the fluid at membrane proximity. Furthermore, higher temperature gradients are verified for z = 2.25 m and 3 m (close to outlet region), due to the gradual heat losses along with the axial position.
In Figure 9, Figure 10, Figure 11 and Figure 12, significant temperature variations can be observed inside the module, due to the pre-established entry conditions, which are equal to ΔT ≈ 2 °C (Case 5), ΔT ≈ 7 °C (Case 6), ΔT ≈ 13 °C (Case 7), and ΔT ≈ 22 °C (Case 8). This behavior is associated with heat losses, by natural convection for the environment, since only the fluid mixture at the inlet has a higher temperature than the equipment wall.

3.5. Concentration Field

Figure 13, Figure 14, Figure 15 and Figure 16 show the concentration fields inside the separation module for feed temperatures of 35, 55, 75, and 95 °C (Cases 5–8). From the analysis of these figures, it is observed that the oil concentration increases along with the axial position, being reduced in 0.05 kg/m3 with the temperature variation from 35 to 95 °C. Furthermore, larger concentration gradients near the concentrate exit region were verified.
This behavior is associated with the oil transport to the proximity of the membrane surface by the shear forces, which is attenuated by increasing the temperature, due to the decrease in dynamic viscosity and density of the fluid, thereby reducing the accumulation of particles along the membrane surface.
This behavior can be seen more adequately in Figure 17, which represents the oil concentration profile along the z-axis, in y equal 0.017 m, for temperatures ranging from 35 to 95 °C.

3.6. Final Consideration and Limitations of the Model

In this work, the steady-state regime was used, being a relevant approach when it is desired to analyze the behavior of the separation process without the influence of time. As the main objective is to evaluate the influence of temperature on the flow-dynamics behavior, this approach is appropriate for the proposed study. However, we recommend a transient analysis for future works, to determine the time to change or cleaning the membrane.
The model developed was based on the work of Damak et al. [23,24,25] (used to validate this work). Therefore, we decided to work similarly. Then, the flow laminar regime was assumed. In general, this condition has been used in different experimental works. In addition, Equation (30), which was used to describe the concentration polarization, is only valid for R = 300–1000, i.e., under low Reynolds number condition. However, in industrial scale, outside of the membrane, turbulent flow is the dominant flow regime. Thus, studies under turbulent regime are recommended, since the turbulent flow could help to reduce the concentration polarization.
Furthermore, based on Damak et al. [23,24,25], it was decided to work with the series resistance model, considering the membrane resistance and the resistance due to the formation of the concentration polarization, as well as the Carman–Kozeny equation, to describe the pressure drop across the porous medium. Then, the effect of the membrane data and thickness were incorporated into the model, not being necessary to model the porous media explicitly.
Furthermore, when disregarding the pore obstruction and the fouling mechanisms, it is not possible to analyze their influence on the process. These effects are important, but the information is not found easily. However, within the proposed objective, these considerations do not compromise the results, since the model considers the rejection coefficient equal to one.
The temperature range for analysis was chosen based on the physics of the problem (temperatures below 100 °C), since high temperatures promote physical phenomena such as water evaporation and oil degradation. In industrial practice, all devices operate in isothermal conditions. Thus, new researches incorporating thermal effects play an important role.
With respect to the heat-transfer process, some aspects need to be explained. Based on the work of Damak and co-authors [23,24,25], the membrane was modeled entirely by boundary conditions. Then, it was necessary to use an appropriated thermal condition at the surface of the membrane. The fluid flows inside the system under action of one piece of pumping equipment; however, since the lower fluid velocity was found in the device (typically, less than 0.05) and Gr/Re2 >> 1 [28], the consideration of free convection to analyze heat transfer, it is very appropriate.
Furthermore, with respect to the choice of the numerical meshes, some important points must be explained. The number of elements used in 2D analysis was adequate. For this case, all predicted results were validated with the results reported by Damak and co-authors [23,24,25]. This information can be verified by the good agreement between the data, as shown in previous work [9]. It is expected that the number of elements in the 3D mesh would be higher than that used in this work. However, because of the simultaneous solution of the governing equations as applied to multiphase flow in porous and non-porous media and the limitation in hardware (computational equipment), larger computational efforts and time were needed, and numerical convergence problems occurred (overflow), when a larger element number was used. Thus, as an alternative numerical procedure, the number of elements was reduced, sequentially, until we obtained a convergent solution with the highest possible element number of the mesh (217,352 elements).

4. Conclusions

Based on the numerical results of the water/oil separation process through the tubular membrane, it is possible to conclude the following:
(a)
The mathematical model was able to predict the fluid behavior during the water/oil separation process with accuracy, providing a good understanding of the flow inside the system;
(b)
The increase in temperature provided modifications in the fluid properties, reducing the feed velocity and modifying the concentration, pressure, and velocity distributions;
(c)
The increase in temperature inside the equipment reduces the pressure gradient, due to a decrease in the inlet velocity, directly influencing the fluid viscosity and density ratio, causing a reduction in the permeation velocity.

Author Contributions

All the authors contributed to the development, analysis, writing, and revision of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CNPq, CAPES, and FINEP (Brazilian research agencies).

Acknowledgments

The authors thank the lab technical support (LPFI and LCTF) and the Brazilian Research Agencies for financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brasil, N.I.; Araújo, M.A.S.; Souza, E.C.M. Oil and Gas Processing Oil and its Derivatives; Primary Processing; Refining Processes; Petrochemicals; Environment; Publisher GEN National Editorial Group: Rio de Janeiro, Brazil, 2014. (In Portuguese) [Google Scholar]
  2. Aljuboury, D.; Palaniandy, P.; Abdul Aziz, H.B.; Feroz, S. Treatment of petroleum wastewater by conventional and new technologies—A review. Glob. Nest J. 2017, 19, 43–452. [Google Scholar]
  3. Santos, M.R.G.; Goulart, M.O.F.; Tonholo, J.; Zanta, C.L.P.S. The application of electrochemical technology to the remediation of oily wastewater. Chemosphere 2006, 64, 393–399. [Google Scholar] [CrossRef]
  4. Dezotti, M. Processes and Techniques for the Environmental Control of Liquid Effluents. Rio de Janeiro: E-Papers 2008, 5, 360. (In Portuguese) [Google Scholar]
  5. Etchepare, R.; Oliveira, H.; Azevedo, A.; Rubio, J. Separation of emulsified crude oil in saline water by dissolved air flotation with micro and nanobubbles. Sep. Purif. Technol. 2017, 186, 326–332. [Google Scholar] [CrossRef]
  6. Jamaly, S.; Giwa, A.; Hasan, S.W. Recent improvements in oily wastewater treatment: Progress, challenges, and future opportunities. J. Environ. Sci. 2015, 37, 15–30. [Google Scholar] [CrossRef]
  7. Padaki, M.; Murali, R.S.; Abdullah, M.S.; Misdan, N.; Moslehyani, A.; Kassim, M.A.; Hilal, N.; Ismail, A.F. Membrane technology enhancement in oil-water separation. A review. Desalination 2015, 357, 197–207. [Google Scholar] [CrossRef]
  8. Magalhães, H.L.F.; De Lima, A.G.B.; Farias Neto, S.R.; Alves, H.G.; De Souza, J.S. Produced water treatment by ceramic membrane: A numerical investigation by computational fluid dynamics. Adv. Mech. Eng. 2017, 9, 1–20. [Google Scholar] [CrossRef] [Green Version]
  9. Magalhães, H.L.F.; Lima, A.G.B.; Farias Neto, S.R.; De Almeida, A.F.; Andrade, T.H.F.; Brandão, V.A.A. Ceramic Membranes: Theory and Engineering Applications. Transp. Phenom. Multiph. Syst. 2018, 111–137. [Google Scholar] [CrossRef]
  10. Vasanth, D.; Pugazhenthi, G.; Uppaluri, R. Cross-flow microfiltration of oil-in-water emulsions using low cost ceramic membranes. Desalination 2013, 320, 86–95. [Google Scholar] [CrossRef]
  11. Zsirai, T.; Al-Jaml, A.K.; Qiblawey, H.; Al-Marri, M.; Ahmed, A.; Bach, S.; Watson, S.; Judd, S. Ceramic membrane filtration of produced water: Impact of membrane module. Sep. Purif. Technol. 2016, 165, 214–221. [Google Scholar] [CrossRef] [Green Version]
  12. Ding, D.; Mao, H.; Chen, X.; Qiu, M.; Fan, Y. Underwater superoleophobic-underoil superhydrophobic Janus ceramic membrane with its switchable separation in oil/water emulsions. J. Membr. Sci. 2018, 565, 303–310. [Google Scholar] [CrossRef]
  13. Frederic, E.; Guigui, C.; Jacob, M.; Machinal, C.; Krifi, A.; Line, A.; Schmitz, P. Modelling of fluid flow distribution in multichannel ceramic membrane: Application to the filtration of produced water. J. Membr. Sci. 2018, 567, 290–302. [Google Scholar] [CrossRef]
  14. Mota, T.V.; Alves, H.G.; Neto, S.R.F.; De Lima, A.G.B. Oily Water Treatment Using Ceramic Membrane in Presence of Swirling Flow Induced by a Tangential Inlet via CFD. Defect Diffus. Forum 2014, 348, 51–57. [Google Scholar] [CrossRef]
  15. Motin, A.; Tarabara, V.V.; Bénard, A. Numerical investigation of the performance and hydrodynamics of a rotating tubular membrane used for liquid-liquid separation. J. Membr. Sci. 2015, 473, 245–255. [Google Scholar] [CrossRef]
  16. Serra, C.A.; Wiesner, M.R. A comparison of rotating and stationary membrane disk filters using computational fluid dynamics. J. Membr. Sci. 2000, 165, 19–29. [Google Scholar] [CrossRef]
  17. Serra, C.A.; Wiesner, M.R.; Laîné, J.M. Rotating membrane disk filters: Design evaluation using computational fluid dynamics. Chem. Eng. J. 1999, 72, 1–17. [Google Scholar] [CrossRef]
  18. Geraldes, V.; Semião, V.; Pinho, M.N. Numerical modeling of mass transfer in slits with semi-permeable membrane walls. Eng. Comput. 2000, 17, 192–217. [Google Scholar] [CrossRef]
  19. Darcovich, K.; Dal-Cin, M.M.; Ballevre, S.; Wavelet, J.P. CFD assisted thin channel membrane characterization module. J. Membr. Sci. 1997, 124, 181–193. [Google Scholar] [CrossRef] [Green Version]
  20. Souza, J.S.; Santos Filho, S.J.; Farias Neto, S.R.; Lima, A.G.B.; Magalhães, H.L.F. Mass Transfer in Tubular Ceramic Membranes for Polluted Water Treatment - Numerical Simulation. Diffus. Found. 2019, 20, 16–33. [Google Scholar] [CrossRef]
  21. Cunha, A.L.; Farias Neto, S.R.; Lima, A.G.B.; Barbosa, E.S.; Santos, J.P.L.; Silva, G.F. Water-oil separation process using a concentric tubular ceramic membrane module: A numerical investigation. Braz. J. Pet. Gas 2016, 10, 205–219. [Google Scholar] [CrossRef]
  22. Alves, H.G.; Magalhães, H.L.F.; Santos, W.R.G.; Araújo, M.V.; Farias Neto, S.R.; Lima, A.G.B. Water/oil separation process via ceramic membranes: Modeling and Simulation. Int. J. Model. Simul. Pet. Ind. 2016, 9, 23–32. [Google Scholar]
  23. Damak, K.; Ayadi, A.; Zeghmati, B.; Schmitz, P. A new Navier-Stokes and Darcy’s law combined model for fluid flow in crossflow filtration tubular membranes. Desalination 2004, 161, 67–77. [Google Scholar] [CrossRef]
  24. Damak, K.; Ayadi, A.; Zeghmati, B.; Schmitz, P. Concentration polarization in tubular membranes—A numerical approach. Desalination 2005, 171, 139–153. [Google Scholar] [CrossRef]
  25. Damak, K.; Ayadi, A.; Schmitz, P.; Zeghmati, B. Modeling of cross-flow membrane separation processes under laminar flow conditions in tubular membrane. Desalination 2004, 168, 231–239. [Google Scholar] [CrossRef]
  26. Pak, A.; Mohammad, T.; Hosseinalipour, S.M.; Allahdinib, V. CFD modeling of porous membranes. Desalination 2008, 222, 482–488. [Google Scholar] [CrossRef]
  27. Minnikanti, V.S.; Dasgupta, S.; De, S. Prediction of mass transfer coefficient with suction for turbulent flow in cross flow ultrafiltration. J. Membr. Sci. 1999, 157, 227–239. [Google Scholar] [CrossRef]
  28. Bergman, T.L.; Lavine, A.S.; Incropera, F.P.; Dewitt, D.P. Fundamentals of Heat and Mass Transfer, 7th ed.; John Wiley & Sons, Inc.: South Tower, Singapore, 2007. [Google Scholar]
  29. Ramires, M.L.V.; Nieto De Castro, C.A.; Nagasaka, Y.; Nagashima, A.; Assael, M.J.; Wakeham, W.A. Standard reference data for the thermal conductivity of water. J. Phys. Chem. Ref. Data 1995, 24, 1377–1381. [Google Scholar] [CrossRef]
  30. Coker, A.K. Ludwig’s Applied Process Design for Chemical and Petrochemical Plants, 4th ed.; Elsevier and Gulf Professional Publishing: Houston, TX, USA, 2010. [Google Scholar]
  31. Cunha, A.L. Treatment of Effluents from the Petroleum Industry via Ceramic Membranes—Modeling and Simulation. Ph.D. Thesis, Federal University of Campina Grande, Campina Grande, Brazil, 2014; 201p. (In Portuguese). [Google Scholar]
  32. Magalhães, H.L.F.; Moreira, G.; Correia, B.R.B.; Gomez, R.S.; Lima, A.G.B.; Farias Neto, S.R. Thermo-Fluid Dynamics Analysis of the Oil-Water Separation Using Ceramic Membrane. Diffus. Found. 2019, 24, 37–60, ISSN 2296–3642. [Google Scholar] [CrossRef]
Figure 1. Geometric representation of the tubular membrane and the 2D transverse plane.
Figure 1. Geometric representation of the tubular membrane and the 2D transverse plane.
Energies 13 02092 g001
Figure 2. The two-dimensional mesh of the tubular-membrane longitudinal section.
Figure 2. The two-dimensional mesh of the tubular-membrane longitudinal section.
Energies 13 02092 g002
Figure 3. Geometric representation of the 3D separation module: (a) membrane in the external tube; (b) membrane in the internal tube.
Figure 3. Geometric representation of the 3D separation module: (a) membrane in the external tube; (b) membrane in the internal tube.
Energies 13 02092 g003
Figure 4. The three-dimensional mesh of the separation module (shell and tube membrane).
Figure 4. The three-dimensional mesh of the separation module (shell and tube membrane).
Energies 13 02092 g004
Figure 5. Oil concentration profile along the z-axis, in y = 0.0225 m (Cases 2 and 9).
Figure 5. Oil concentration profile along the z-axis, in y = 0.0225 m (Cases 2 and 9).
Energies 13 02092 g005
Figure 6. Concentration profile inside the separation module, at the angular positions (a) 0° and 180°, and (b) 90° and 270° (Cases 5 and 8).
Figure 6. Concentration profile inside the separation module, at the angular positions (a) 0° and 180°, and (b) 90° and 270° (Cases 5 and 8).
Energies 13 02092 g006
Figure 7. Permeation velocity as a function of axial position along the membrane (Cases 1–8).
Figure 7. Permeation velocity as a function of axial position along the membrane (Cases 1–8).
Energies 13 02092 g007
Figure 8. Velocity field for temperatures of (a) 35 °C, (b) 55 °C, (c) 75 °C and (d) 95 °C (Cases 5–8).
Figure 8. Velocity field for temperatures of (a) 35 °C, (b) 55 °C, (c) 75 °C and (d) 95 °C (Cases 5–8).
Energies 13 02092 g008
Figure 9. Temperature field inside the shell and tube module for T0 = 35 °C (Case 5).
Figure 9. Temperature field inside the shell and tube module for T0 = 35 °C (Case 5).
Energies 13 02092 g009
Figure 10. Temperature field inside the shell and tube module for T0 = 55 °C (Case 6).
Figure 10. Temperature field inside the shell and tube module for T0 = 55 °C (Case 6).
Energies 13 02092 g010
Figure 11. Temperature field inside the shell and tube module for T0 = 75 °C (Case 7).
Figure 11. Temperature field inside the shell and tube module for T0 = 75 °C (Case 7).
Energies 13 02092 g011
Figure 12. Temperature field inside the shell and tube module for T0 = 95 °C (Case 8).
Figure 12. Temperature field inside the shell and tube module for T0 = 95 °C (Case 8).
Energies 13 02092 g012
Figure 13. Oil concentration field inside the shell and tube module for T0 = 35 °C (Case 5).
Figure 13. Oil concentration field inside the shell and tube module for T0 = 35 °C (Case 5).
Energies 13 02092 g013
Figure 14. Oil concentration field inside the shell and tube module for T0 = 55 °C (Case 6).
Figure 14. Oil concentration field inside the shell and tube module for T0 = 55 °C (Case 6).
Energies 13 02092 g014
Figure 15. Oil concentration field inside the shell and tube module for T0 = 75 °C (Case 7).
Figure 15. Oil concentration field inside the shell and tube module for T0 = 75 °C (Case 7).
Energies 13 02092 g015
Figure 16. Oil concentration field inside the shell and tube module for T0 = 95 °C (Case 8).
Figure 16. Oil concentration field inside the shell and tube module for T0 = 95 °C (Case 8).
Energies 13 02092 g016
Figure 17. Oil concentration profile along the z-axis, at y = 0.017 m (Cases 5–8).
Figure 17. Oil concentration profile along the z-axis, at y = 0.017 m (Cases 5–8).
Energies 13 02092 g017
Table 1. Geometrical dimensions of the separation module.
Table 1. Geometrical dimensions of the separation module.
Geometrical DimensionsModule Dimensions (mm)
Length, L3000
Internal diameter, di30
External diameter, De60
a = b15
c20
Table 2. Regression coefficients for the specific heat of the water. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
Table 2. Regression coefficients for the specific heat of the water. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
A (J mol−1 K−1)B (J mol−1 K−2)C (J mol−1 K−3)D (J mol−1 K−4)
92.053−3.9953 × 10−2−2.1103 × 10−45.3469 × 10−7
Table 3. Regression coefficients for water density. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
Table 3. Regression coefficients for water density. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
A (kg)B (-)N (−)Tc (K)
0.34710 × 1030.274000.28571647.13
Table 4. Regression coefficients for the dynamic viscosity of water. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
Table 4. Regression coefficients for the dynamic viscosity of water. Reproduced from [30], Elsevier and Gulf Professional Publishing: 2010.
A (-)B (K)C (K−1)D (K−2)
−10.21581.7925 × 1031.7730 × 10−2−1.2631 × 10−5
Table 5. Parameters considered constant in simulations.
Table 5. Parameters considered constant in simulations.
Permeation Reynolds NumberRew0.1
Initial concentration (kg⋅m−3)C01
Schmidt numberSc1000
Axial Reynolds numberRe1000
Mean diameter of the oil droplet (m)dp51
Intrinsic retention of the oil by the membraneRr1
Membrane thickness (m)e0.01
External pressure (Pa)Pex101,325
Porosityεp0.3
Permeability (m2)K3.33 × 10−11
Table 6. Conditions used in the simulations.
Table 6. Conditions used in the simulations.
CasesT0 (°C)Membrane Location
Two-dimensional analysis
135-
255-
375-
495-
Three-dimensional analysis
535Internal tube (Figure 3b)
655Internal tube (Figure 3b)
775Internal tube (Figure 3b)
895Internal tube (Figure 3b)
955External tube (Figure 3a)
Table 7. Permeate mass flow for different temperatures.
Table 7. Permeate mass flow for different temperatures.
TemperaturePermeate Mass Flow (kg/s)
2D3D
35 °C4.28158 × 10−51.70724 × 10−4
55 °C2.94879 × 10−51.14345 × 10−4
75 °C2.18247 × 10−58.25664 × 10−5
95 °C1.70792 × 10−56.32973 × 10−5

Share and Cite

MDPI and ACS Style

Magalhães, H.L.F.; Moreira, G.; Gomez, R.S.; Porto, T.R.N.; Correia, B.R.B.; Silva, A.M.V.; Farias Neto, S.R.; Lima, A.G.B. Non-Isothermal Treatment of Oily Waters Using Ceramic Membrane: A Numerical Investigation. Energies 2020, 13, 2092. https://doi.org/10.3390/en13082092

AMA Style

Magalhães HLF, Moreira G, Gomez RS, Porto TRN, Correia BRB, Silva AMV, Farias Neto SR, Lima AGB. Non-Isothermal Treatment of Oily Waters Using Ceramic Membrane: A Numerical Investigation. Energies. 2020; 13(8):2092. https://doi.org/10.3390/en13082092

Chicago/Turabian Style

Magalhães, Hortência L. F., Gicelia Moreira, Ricardo S. Gomez, Túlio R. N. Porto, Balbina R. B. Correia, Anderson M. V. Silva, Severino R. Farias Neto, and Antonio G. B. Lima. 2020. "Non-Isothermal Treatment of Oily Waters Using Ceramic Membrane: A Numerical Investigation" Energies 13, no. 8: 2092. https://doi.org/10.3390/en13082092

APA Style

Magalhães, H. L. F., Moreira, G., Gomez, R. S., Porto, T. R. N., Correia, B. R. B., Silva, A. M. V., Farias Neto, S. R., & Lima, A. G. B. (2020). Non-Isothermal Treatment of Oily Waters Using Ceramic Membrane: A Numerical Investigation. Energies, 13(8), 2092. https://doi.org/10.3390/en13082092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop