Abstract

In the present work, the equilibrium, thermodynamics, and kinetics of boron removal from aqueous solutions by the adsorption on commercial magnesium oxide powder were studied in a batch reactor. The adsorption efficiency of boron removal increases with temperature from 25°C to 50°C. The experimental results were fitted to the Langmuir, Freundlich, and Dubinin–Radushkevich (DR) adsorption isotherm models. The Freundlich model provided the best fitting, and the maximum monolayer adsorption capacity of MgO was 36.11 mg·g−1. In addition, experimental kinetic data interpretations were attempted for the pseudo-first-order kinetic model and pseudo-second-order kinetic model. The results show that the pseudo-second-order kinetic model provides the best fit. Such result suggests that the adsorption process seems to occur in two stages due to the two straight slopes obtained through the application of the pseudo-first-order kinetic model, which is confirmed by the adjustment of the results to the pseudo-second-order model. The calculated activation energy () was 45.5 kJ·mol−1, and the values calculated for G°, H°, and S° were −4.16 kJ·mol−1, 21.7 kJ·mol−1, and 87.3 kJ·mol−1, respectively. These values confirm the spontaneous and endothermic nature of the adsorption process and indicated that the disorder increased at the solid-liquid interface. The results indicate that the controlling step of boron adsorption process on MgO is of a physical nature.

1. Introduction

The exploration and production of petroleum from marine subsoil generate large amounts of liquid effluents, also known as produced water, extracted along with crude oil [1]. Produced water is the water found along with oil in the marine reservoir or the water injected into the reservoir in order to recover petroleum. The produced water flow is low at the early stages of production of the reservoir. Nevertheless, it can reach up to 80% of the crude oil extraction in the final years of exploitation of the well [2].

The produced water composition can vary according to the well and usually contains high salinity, organic and inorganic substances, and levels of dissolved solids over 40%, which makes it toxic to the environment [2, 3].

Boron is an element present in the produced water, and its concentration may vary from 4 mg·L−1 to 350 mg·L−1 [3, 4]. Although boron is considered a micronutrient essential to the development of microorganisms, plants, microalgae, and animals, this element can be toxic at concentration of 0.3 mg·L−1 to sensitive plants, at 2 mg·L−1 to semitolerant ones, and at 4 mg·L−1 to tolerant ones. Therefore, this is one of the reasons why boron is required to be removed from water and other effluents by environmental protection agencies [5, 6].

The maximum boron concentration recommended by the World Health Organization (WHO) for potable water was 0.5 mg·L−1 in 1998. However, this value was modified to 2.4 mg·L−1 in 2011. According to standards set by legislation, the limit allowed by the European Union, the UK, and Japan is 1.0 mg·L−1. In South American countries such as Brazil and Peru, the limits for fresh water and wastewater are 0.5 mg·L−1 and 5 mg·L−1, respectively. On the contrary, in the USA, the limit is not subject to federal regulations on this issue. The states of Minnesota, Florida, and California have allowed limits of 0.6, 0.63, and 1 mg·L−1, respectively [711].

Consequently, the produced water containing boron requires treatment before being released back into the sea or being employed as a hydrosource of potable water.

There are many possible treatments for the boron removal such as the electrocoagulation, process which achieved efficiencies over 98% from synthetic solutions and real produced waters with initial boron concentration between 10 and 30 mg·L−1 [12]. Other treatments are the adsorption by different types of adsorbents, including activated carbon, fly ash, clay, mesoporous silica, oxides, nanoparticles, biological material, layered double hydroxides, and natural minerals [1315]. Some additional methods studied at laboratory, for example, electrodialysis [16], phytoremediation [17, 18], and bioelectrochemical systems, presented removal efficiencies up to 90% [19, 20].

However, the most extensive processes employed for boron removal from produced water and seawater, on a larger scale, are reverse osmosis and ion exchange [21, 22].

As far as reverse osmosis is concerned, boron removal can reach values of 98% at pH = 10.5 since, at this pH, the predominant boron specie is a borate ion B(OH)4− that has a negative charge and larger size when compared to the boric acid (H3BO3) [2, 22]. In the case of ion exchange, the removal of boron occurs via its adsorption in specific resins, which are usually synthesized by means of macroporous cross-linked polystyrene resins, and functionalized by the N-methyl-D-glucamine (1-amino-1-deoxy-D-glucitol; NMDG) group achieving removal efficiencies up to 99% [2326].

The most common commercial resins used in ion exchange process are Diaion CRB 02, Purolite S 108, and Amberlite IRA-743 [6, 2730].

In general, the initial concentration of boron in produced waters treated with resins can vary between 15 and 60 mg·L−1. However, after a certain period of time, this resin saturates, reaching the breakthrough, and needs to be regenerated to be reused [12, 31].

Regeneration step consists of passing an acid solution (usually H2SO4) through a saturated resin, resulting in an acid solution with high concentration of boron, usually between 350 and 700 mg·L−1, which must also be treated.

The volume of this effluent containing boron is relevant when it refers by example to oil waterway terminals, such as São Sebastião Waterway Terminal from Petrobras/Brazil, with a capacity of 1,585,345 m3 [32]. Due to their magnitude, they generate a large quantity of effluent, which is product of the elution of the resins used at wastewater treatment plants of these terminals.

The treatment of this effluent by conventional technologies is not adequate due to the high concentration of boron. Therefore, under this circumstance, technologies such as precipitation or adsorption seem to be viable alternatives.

Through these treatments, boron concentration could be reduced up to 15–30 mg·L−1 that may allow the mixing of treated effluent with production water and these could be subsequently treated by ion exchange resins or reverse osmosis. The resulting solid residue of precipitation/adsorption process can be disposed in landfills or be used as raw material in the manufacture of glass due to its high content of boron. In this context, magnesium oxide is provided as a good alternative as it is environmentally friendly, cost-effective, and nontoxic. Furthermore, it shows low solubility in water, and it is an effective sorbent for the removal of contaminants such as fluoride and toxic dyes [33, 34].

The aim of the present work was to evaluate the kinetics, thermodynamics, and equilibrium of boron removal from aqueous solutions with high content (350 mg·L−1) through adsorption in magnesium oxide.

2. Experimental

2.1. Materials and Methods

For each experiment, a synthetic solution of 350 mg·L−1 was prepared by dissolution of boric acid (H3BO3) PA in distilled water. The pH of the solution was adjusted with 1 M NaOH and/or 1 M HCl solutions provided by Sigma-Aldrich.

The magnesium oxide used in the adsorption process was supplied by Magnesita SA (MgO-500).

2.2. Characterization of MgO and Boron Concentration Analysis

The sizes of the MgO particles were analyzed by a dual laser liquid dry dispersion particle size analyzer, CILAS 1064L. Boron concentrations in aqueous solutions were determined by optical emission spectrometry analysis with inductively coupled plasma from PerkinElmer, Optima 4300DV. The morphology of the MgO was determined by scanning electron microscopy (MRV).

2.3. Experimental Procedure

The experiments were carried out in a batch reactor with 500 mL of solution. The use of borosilicate glass beaker was to prevent further contamination by dissolution of boron substances from the silicate material (beaker), mainly because the experiments were performed at alkaline pH. The solution was stirred at a speed of 150 rpm, 40 g·L−1 of MgO was added, and the pH was readily adjusted after starting the experiment. The pH of the solution was measured using a pH meter (Bel Engineering W3B) and adjusted with 1 M NaOH and/or 1 M HCl solutions (Sigma-Aldrich).

The study of reaction kinetic was performed by collecting the sample every 5 minutes. The samples were vacuum-filtered through a cellulose nitrate membrane of 0.4 μm pore diameter (Nalgene). After 240 minutes of reaction, the stirring was turned off.

The filtered samples were kept for boron concentration analysis by inductively coupled plasma spectrometer (Optima 4300DV, from PerkinElmer).

Simultaneously, a sludge sample was dried at 70°C and adequately stored for further analysis of the surface morphology by scanning electron microscopy (MRV).

The effect of temperature was evaluated in the process of boron removal at range of 5°C to 50°C, with initial concentration of boron of 350 mg·L−1. For constant temperature maintenance during the experiments, a cooling or heating system was used.

2.3.1. Adsorption Isotherms

The experiments for adsorption isotherm evaluation were carried out in a batch system with MgO concentration range of 8 g·L−1 to 64 g·L−1 at three different temperatures of 25°C, 40°C, and 50°C, during 240 min of reaction.

The initial concentration of boron was 350 mg·L−1, the stirring speed was 150 rpm, the pH was 10, the temperature was at 25°C, and the volume of the solution was 500 mL. After filtration, the samples were sent to boron concentration analysis.

3. Results and Discussions

3.1. Characterization of MgO

The chemical composition provided by the manufacturer is shown in Table 1.

Figures 1 and 2 show the surface morphology of MgO-500 before and after the adsorption process, respectively. The difference in structure presented in Figure 2 was due to the formation of a layer of H3BO3 on the surface of MgO and the hydration of MgO into Mg(OH)2 after the adsorption process.

3.2. Adsorption Isotherms

The adsorption isotherms describe the equilibrium between the concentration of a material in aqueous phase and its concentration on the surfaces of particle adsorbents. This study employed the Langmuir, Freundlich, and Dubinin–Radushkevich models to describe the equilibrium adsorption.

3.2.1. Dubinin–Radushkevich Isotherm

Dubinin–Radushkevich (DR) isotherm is a model that considers the adsorption in multilayers and in the heterogeneous surfaces.

DR model is expressed mathematically as follows:where is the amount of boron adsorbed per MgO (mg·g−1), is the maximum adsorption capacity, is the constant related to sorption energy (mol2·kJ−2), and ∈ is the Polanyi potential.and are expressed by (2) and (3), respectively [35]:where is the gas constant (J·mol−1 K−1) and is the temperature (K). (kJ·mol−1) is the sorption energy, and the magnitude of its value indicates if the adsorption is of chemical or physical nature.

3.2.2. Langmuir Isotherm

The Langmuir isotherm is a model that considers monolayer adsorption onto a uniform surface with a finite number of adsorption sites and uniform adsorption energy, and this model is given by the following equation [36]:

Equation (4) can be linearized as follows:where is the equilibrium concentration of boron in solution (mg·L−1), is the quantity of boron adsorbed onto the MgO (mg·g−1), is the maximum monolayer adsorption capacity of MgO (mg·g−1), and is the Langmuir adsorption constant related to the energy sorption (L·mg−1).

3.2.3. Freundlich Isotherm

This model is applied to the sorption processes on heterogeneous surfaces and reversible adsorption and admits multilayer adsorption [37].

The Freundlich equation and its linear form can be given as follows:where is the boron concentration adsorbed at equilibrium (mg·g−1), is the equilibrium concentration of boron in solution (mg·L−1), and (L·g−1) and n are the Freundlich sorption isotherm constants, related to the adsorption capacity and the adsorption intensity, respectively.

Figure 3 shows the adsorption isotherms of boron MgO-500 at 25°C, 40°C, and 50°C.

It can be observed that the adsorption isotherm at 25°C could be classified as an H4-type isotherm according to the classification of Giles et al. [38]. This type of isotherm suggests a high affinity between adsorbate-adsorbent, and the subgroup 4 suggests the formation of monolayers of adsorbate on the surface of adsorbent.

On the contrary, for the isotherms at 40°C and 50°C, it is observed that they are similar L-type isotherms according to the classification of Giles et al. [38].

This type of isotherm assumes the existence of an affinity between ion B(OH)4- and MgO, and if more sites of the adsorbent are filled, it will be more difficult to fill the empty sites with other solute molecules. This type of isotherm is commonly represented under the following mechanism: (1) molecules are adsorbed in layers and (2) there is competition for the active sites on the adsorbent surface between adsorbate molecules and solvent molecules [3941].

In the present work, the Langmuir, Freundlich, and Dubinin–Radushkevich models were evaluated, which are the most widely used models to describe an adsorption process.

According to the results shown in Table 2, the model that provides a better fit to the experimental data is the isotherm of Freundlich at 25°C, 40°C, and 50°C with 99% correlation. Figure 4 shows the data adjusting to the Freundlich isotherm, and Table 3 compares the Freundlich isotherm parameters between MgO used at present work and other adsorbents.

The Freundlich isotherms suggest a heterogeneous adsorbent surface and a reversible adsorption process, which considers the formation of multilayer. However, the adsorption of ion B(OH)4- may involve different mechanisms, such as ion exchange, microprecipitation, complexation/chelation, and electrostatic attraction [4042].

The Freundlich isotherms were obtained for values of 0.33, 0.51, and 0.57 at temperatures of 25°C, 40°C, and 50°C, respectively. The values of n, which are related to the distribution of ion B(OH)4− linked to the active sites on the adsorbent, were 1.392, 1.476, and 1.481, respectively, for 25°C, 40°C, and 50°C (Table 2).

The value of n is a constant that indicates the strength of adsorption and is also known as a measure of linearity since n is equal to one. In this case, the adsorption process is linear; therefore, the adsorption sites are homogeneous concerning energy, and there would not be any interactions between the adsorbate and the adsorbent. For values of n smaller than the unit, the connection between the adsorbate and the adsorbent is very weak, and consequently, the process is not favorable for the adsorption and its adsorptive capacity decreases.

However, for values of n greater than the unit, the adsorption process is favorable; therefore, the adsorption capacity increases, indicating that the adsorption of ions is favorable under experimental conditions studied [4143]. As seen in Table 3, all values of n are greater than the unit for the Freundlich isotherm, and the temperature increases from 25°C to 50°C, indicating that the adsorption of boron MgO-500 is favored by the temperature.

It can also be noticed from Table 3 that the adsorption constants () for boron removal process at different temperatures are higher than those obtained for boron removal processes with activated charcoal, activated alumina, Pural (76% Al2O3), Siral (28% SiO2, 72% Al2O3), tannin gel, and mineral wastes. It is clear that the capacity of boron adsorption by MgO is higher than other adsorbents with exception to N-methyl-D-glucamine onto SPC and neutralized red mud.

3.3. Kinetic Analysis

The kinetics of boron removal by adsorption onto MgO was studied in the present work. The rate in which boron is removed by adsorption is the most important parameter for the design of a continuous or batch system in an effluent treatment plant. Therefore, it is essential to establish the dependency of time in the adsorption process for different operational conditions.

The kinetics of adsorption depends on the adsorbate-adsorbent interaction. The rate of adsorbate removal determines the residency time required for complete adsorption and could be calculated by kinetic analysis.

3.3.1. Equilibrium Kinetics

During an adsorption process, the sorbed solute tends to desorb and return to the solution, and vice versa. The process occurs continuously until, at a given time, the adsorption and desorption rates reach an equilibrium state. Hence, in this stage, there will not be any additional adsorption of the pollutant in the solution [4547].

Tests were conducted to determine the equilibrium of the reaction under the following conditions: 500 mL solution with a concentration of 350 mg·L−1 of boron was prepared; then, MgO was added to obtain a concentration of 40 g·L−1, maintaining the pH at 10 and stirring speed of 150 rpm. Samples were collected at different time periods (5, 10, 15, 20, 25, 30, 40, 50, 60, 90, 140, 190, and 240 min).

Figure 5 shows the effect of temperature on boron removal.

The maximum adsorption capacities were 7.2, 7.7, and 8.0 mg·g−1 at temperatures of 25°C, 40°C, and 50°C, respectively.

It can also be observed, according to Figure 5, that the rate of the process is greatly affected by an increase in temperature, and at 25°C, 40°C, and 50°C, the equilibrium of reaction was reached around 40 min, 30 min, and 15 min, respectively. From this stage, the variations in adsorption capacities were insignificant. Figure 5 also shows that the period of time for reaching equilibrium at temperatures of 40°C and 50°C is relatively short, which could be an advantage for future industrial applications.

3.3.2. Kinetic Models

In 1898, Lagergren presented the kinetic model known as “pseudo-first-order equation” for a liquid-solid system based on the adsorption capacity of the solid.

Later, Ho and Mckay [46] presented a modification of the Lagergren model based on the concentration of the pollutant in the solution, known as “pseudo-second-order equation.”

Next, a kinetic analysis based on those models will be presented, which are the most usual and well-established ones for this type of adsorption process.

(1) Pseudo-First-Order Kinetic Model. The Lagergren model was developed for solid-liquid adsorption systems and is based on the adsorption capacity of the adsorbent. Equation (7) shows how this model is usually expressed, being one of the most used models to study the kinetics of adsorption processes:

Figures 6 and 7 present the application of the experimental data to the pseudo-first-order model at 25°C, 40°C, and 50°C, respectively.

It is observed in Figure 6 that there is an apparent multilinearity of the kinetic curves, which may lead us to think that the process happens in two stages. Figure 7 shows the good fit of the experimental kinetic data at 25°C in the pseudo-first-order model with R2 = 0.96.

(2) Pseudo-Second-Order Kinetic Model. The pseudo-second-order model is based on the concentration of the pollutant in the solution. The adsorption rate equation can be written as follows:where is the adsorption rate constant (g·mg−1·min−1), is the concentration of boron adsorbed at equilibrium (mg·g−1), and is the concentration of boron on the surface of the adsorbent at time (mg·g−1). The constants can be determined experimentally by plotting versus .

Figure 8 illustrates the adjustment of the experimental data to the pseudo-second-order kinetic model. Table 4 presents the values of the kinetic parameters obtained from this setting.

It can be observed from the values presented in Table 4 that the pseudo-second-order model is the one with the best fit of the experimental data of the adsorption process, since R2 values of 0.99 were obtained for this model at the three temperatures studied (25°C, 40°C, and 50°C) compared to the R2 values obtained for the pseudo-first-order model of 0.96, 0.83, and 0.43 at the same temperatures, respectively.

The values of obtained at 25°C, 40°C, and 50°C by the pseudo-second-order model were 7.81 mg·g−1, 8.12 mg·g−1, and 8.17 mg·g−1, respectively, while the experimental values obtained at the same temperatures were 7.23 mg·g−1, 7.75 mg·g−1, and 7.98 mg·g−1, respectively. It can be observed from the comparison of the values of obtained experimentally and by the pseudo-second-order model that the values of are very similar, which confirms that the process of removal of boron by adsorption with MgO follows a pseudo-second-order kinetics.

Studies available in the literature show that many models provide a clarification for the overall adsorption process. However, in many cases, that is not possible when graphs present multilinear feature [48].

In order to understand these cases, it is usual to divide the graph into two or more straight lines and to suggest that the adsorption mechanism is controlled by each straight line. This practice may help to understand the adsorption mechanism to some extent [48].

From Figure 6, which represents the application of the experimental data to a pseudo-first-order kinetic model, it is observed that the adsorption process at the three temperatures seems to occur in two stages due to the two straight slopes obtained from the application of the model, which is confirmed by the adjustment of the results to the pseudo-second-order model that describes a two-step reaction occurring consecutively.

3.3.3. Apparent Activation Energy

The apparent activation energy of the adsorption process was calculated. The variation in the rate constant according to an increase in temperature can be described by the Arrhenius equation (9): where is the sorption rate constant (g·mg−1·min−1), is the independent factor of temperature (g·mg−1·min−1), is the adsorption activation energy (kJ·mol−1), is the gas constant (8.314 J·mol−1·K−1), and is the temperature of the solution (K).

Equation (10) can be expressed in linear form as follows:

As shown in Figure 8, there is a linear correlation between the rate constant of the pseudo-second-order model and its corresponding absolute temperature, with a correlation coefficient of 0.98. Figure 9 shows the correlated experimental data in the linearized Arrhenius equation. From Figure 9, the relation between and can be represented as follows:

From (11), we can observe that the frequency factor is 3.02 × 106 g·mg−1·min−1, and the activation energy is 45.54 kJ·mol−1, which is very slightly outside of the suggested normal range (8–40 kJ·mol−1) for a typical physical adsorption processes [4952].

3.4. Thermodynamics of the Process

In order to understand the thermodynamics of the process of boron removal by adsorption with MgO, some thermodynamic parameters were determined. Adsorption tests at different temperatures were performed using a thermostatic bath, under the following conditions: initial boron concentration = 350 mg·L−1, pH = 10, MgO concentration = 40 g·L−1, and stirring speed = 150 rpm.

The temperatures studied were 25°C, 40°C, and 50°C, and the parameters calculated were standard Gibbs free energy variation (∆G°), standard enthalpy variation (∆H°), and standard entropy variation (∆S°).

The standard free energy variation in the adsorption process is related to the equilibrium constant () and can be calculated according to the following equation:where is the gas constant (8.314 J·mol−1·K−1), is the absolute temperature, and is the equilibrium constant, which can be estimated from the following equation:where  is the concentration of boron at equilibrium in the solution (mg·L−1) and  is the concentration of boron at equilibrium in the adsorbent (mg·L−1).

The variations in standard enthalpy (∆H°) and standard entropy (∆S°) can be calculated according to the following Van’t Hoff equation:

Figure 9 represents the plot of versus 1/T from the Van’t Hoff equation, where the slope of the line is and the intercept with y-axis is , from which the corresponding thermodynamic parameters were obtained and are presented in Table 5.

From Table 5, it can be observed that when the temperature of the adsorption process increases from 25°C to 40°C and from 40°C to 50°C, the ∆G° values become increasingly negative varying from −4.16 kJ·mol−1 to −5.88 kJ·mol−1 and from −5.88 kJ·mol−1 to −6.26 kJ·mol−1, respectively. That means the process is spontaneous at all temperatures due to the negative value of ∆G°. Likewise, as the temperature increases, the B(OH)4- ion has higher affinity to be adsorbed by MgO.

The ∆G° for physical sorption is between −20 and 0 kJ·mol−1 and is ranged from −80 to −400 kJ·mol−1 for chemical sorption [52, 53]. The values of ∆G° presented in Table 5 indicate that the sorption process is controlled by physical adsorption.

The adsorption enthalpy is also a parameter used to indicate the intensity of the interaction between the adsorbate and the adsorbent. In the phenomenon of physisorption, this parameter has low values (up to about 40 kJ·mol−1) since it is characterized by a low degree of interaction, being the forces involved of the order of magnitude of van der Waals forces. The chemisorption phenomenon is characterized by a high degree of interaction between the adsorbate and the surface of the adsorbent; the enthalpy values, in this case, are around 800 kJ·mol−1 [51]. The ∆H° calculated in this work was 21.75 kJ mol−1·K−1, indicating that the sorption process is controlled by physical adsorption.

Positive value of ∆S° (87.33 kJ·mol−1) indicates that the degrees of freedom increase at the solid-liquid interface during the adsorption of boron onto magnesium oxide particles.

4. Conclusions

Equilibrium, kinetics, and thermodynamic studies for the adsorption of boron onto magnesium oxide powder were carried out.

The experimental data of adsorption were conducted by the Langmuir, Freundlich, and DR models, wherein the Freundlich isotherm was found to have a better fit for the equilibrium data for adsorption of boron under ambient temperature (25°C) and at higher temperatures (40°C and 50°C), respectively.

Additionally, the investigation of the kinetics of the overall adsorption process was conducted for the pseudo-first-order kinetic and pseudo-second-order kinetic models. Results show that the pseudo-second-order kinetic model generates the best fit to all experimental data. Such result suggests that the adsorption process at these temperatures seems to occur in two stages due to the two straight slopes obtained through the application of the pseudo-first-order kinetic model, which is confirmed by the adjustment of the results to the pseudo-second-order model.

The calculated activation energy () was 45.54 kJ·mol−1, which is very slightly outside the normal range (8–40 kJ·mol−1) for a typical physical adsorption process.

The values calculated for ∆G° and ∆H° were −4161.43 kJ·mol−1 and 21.75 kJ·mol−1, respectively. These values confirm the spontaneous and endothermic nature of the process and also suggest a physical adsorption process.

Data Availability

The boron analysis data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The funding for this research was provided by the Research and Development Center of Petrobras and the National Council for Scientific and Technological Research of Brazil (CNPq). The authors thank Mariana Lima for assistance during the experimental work.