Abstract
This study proposes the use of activated charcoal made from Umbaúba wood as an adsorbent for the removal of naphthenic acid in an aviation kerosene model mixture. The activated charcoal was characterised as mesoporous with a carbon graphite profile and presented pHpzc equal to 10.5. The best working conditions were obtained for activated charcoal levels of <0.09 mm and 300 r min−1. The system reached the equilibrium after 360 min, without significant statistical difference for the pseudo-first- and pseudo-second-order kinetic models. The Weber–Morris and Boyd models corroborated the conclusion that adsorption is not controlled only by the intraparticle diffusion step. For the equilibrium study, the adsorptive capacity obtained was of 1.1 g g−1, with the Brunauer–Emmett–Teller model better correlating with the experimental data. Given the results obtained, the activated charcoal demonstrated to have a remarkable potential for removing naphthenic acid in an aviation kerosene model mixture.
Introduction
The sources of crude oil in South America are the most acid in the world, consisting of heavy oils with a high relative density (<20° API), high viscosity and high levels of total acidity. This high acidity can lead to more severe corrosion due to the presence of certain compounds, such as naphthenic acids, which are also responsible for forming insoluble salts and deposits under heating (Farah, 2012; Gruber et al., 2012).
Naphthenic acids, when present in aviation kerosene (AVK), are the main source of corrosion in jet engines (Farah, 2012). For being chemically stable and non-volatile, these acids can act as emulsion surfactants. Furthermore, they are highly toxic for several organisms, with a harmful effect to the environment, mainly to aquatic environments (Clemente and Fedorak, 2005).
Several ways of dealing with this problem have already been suggested, such as liquid–liquid extraction (Shi et al., 2008), biodegradation (Quesnel et al., 2011) and catalytic esterification (Wang et al., 2014). However, these technologies are not entirely efficient or are not economically viable in the treatment of these compounds, due to their restricted use, associated with high costs, the destruction of the molecules of naphthenic acids and the treatment time.
Among the treatment processes suggested, adsorption has proved to be an efficient method for the removal of contaminants, as well as offering the possibility of regeneration, recovery and recycling of the adsorbent material (Eren et al., 2010).
The study of the removal of naphthenic acids, by adsorption, is relevant due to the non-destructive character of the process, thus enabling the purification and commercialisation of the acids removed. This allows the development of recycling methods of these removed acids, justifying the economic viability of the removal process (Gruber et al., 2012).
Among the adsorbents recommended, activated charcoal is one of the most used for the removal of industrial contaminants. The material is industrially produced by the oxidative pyrolysis of wood, mineral coal, bones and coconut shell, among others. Its adsorption capacity depends on the precursor material as well as on the activation and carbonisation methods used (Gonzalez and Cid, 2005).
The use of agro-industrial residues as adsorbents is considered to reduce the environmental impact in two different ways, with the polluting residual biomass accumulated being removed from its manufacturing plant or from the place it is deposited, and allowing the residual water to be treated with this biomass (Bhatnagar and Sillanpää, 2010).
In this context, researchers have investigated the development of new types of activated charcoal made from biomasses, such as nut shells and olive seeds (Martínez et al., 2006), mandioca shell (Sudaryanto et al., 2006), pine fruit (Royer et al., 2009), coconut shell (Cazetta et al., 2011), cotton fibres (Sun et al., 2012), bamboo (González et al., 2014), rice hulls (Muniandy et al., 2014), pine tree (Tonucci et al., 2015), potato residue (Zhang et al., 2015), papaya shell (Abbaszadeh et al., 2016) and guava seeds (Pezoti et al., 2016) and used in the removal of pollutants, such as heavy metals and organic compounds.
The use of Umbaúba wood as a precursor for the manufacturing of activated charcoal is advantageous due to the rapid growth of the tree, making it a renewable source of feedstock. The Umbaúba tree is characteristic of riparian forests, of large forest clearing, and is found from Mexico until Argentina, being dominant throughout all of the Brazilian territory (Hernández-Terrones et al., 2007). In the literature consulted, no reports were found of its use as a precursor in the preparation of activated carbon, nor of its application in the removal of naphthenic acids in AVK.
Taking this into account, this work was aimed at physically evaluating the use of activated carbon with CO2 made from the Umbaúba wood, as an adsorbent for the removal of naphthenic acids present in the model mixture of AVK and also determining the kinetic and isothermal parameters of the adsorption process.
Materials and methods
Preparation and quantification of the AVK model mixture
The n-dodecanoic carboxylic acid was chosen as the model compound to represent naphthenic acids, due to its physical properties and similarity with the acid compounds frequently found in AVK, according to Nascimento et al. (2014).
The stock solution of the 8% (v/v) AVK model mixture was prepared from 6 g of n-dodecanoic carboxylic acid (Merck, 99% purity) dissolved in 100 ml of n-dodecanoic solvent (Merck, 99% purity). The working solutions (1–6%) were obtained by diluting the stock solution, pH 6 (natural solution pH).
The concentration analyses of n-dodecanoic acid were carried out by gas chromatography, using a Varian Chromatograph CP-3380 (USA), with a Carbowax capillary column (30 m × 0.32 mm × 0.25 µm) and a flame ionisation detector, under the following conditions: oven temperature at 220°C, temperature of the vaporiser at 250°C, temperature of the detector set to 260°C, H2 used as the mobile phase, 2.0 ml min−1 outflow pressure of 5 psi and injected volume of 1 µl. The samples were analysed before and after each experiment, with blank experiments being carried out following the same procedures.
Charcoal preparation and activation
The charcoal made from Umbaúba wood, used in the activation process, was provided by Elephant Chemical Industry Ltd.
The physical activation process took place in two separate steps. In the first step, the Umbaúba wood was submitted to a 10°C min−1 heating ramp, starting at ambient temperature (25 ± 2°C) up until 600°C, in a N2 inert atmosphere at an outflow of 100 ml min−1, being kept at this temperature for 60 min. Subsequently, the N2 outflow was changed to CO2, at an outflow of 100 ml min−1 for an additional 60 min. After the material had reached room temperature, it was classified in a Tyler sieves series into three different groups of particle sizes <0.09, 0.09–0.15 and 0.15–0.21 mm.
Activated charcoal characterisation
The activated charcoal was characterised through the following techniques:
Point of zero charge. The value of pH at the point of zero charge (pHpzc) of the activated (Cc) charcoal was estimated using the pH measurement of water before and after the contact with the solids. A total of 0.1 g of charcoal for 25 ml of the solution with pH ranging from 2 to 11 was used in the trial, being continuously stirred for 24 h. The pH of the solutions was adjusted using a pH meter (Quimis, Q400AS) with HCl and NaOH (0.1 mol l−1), respectively.
Textural characterisation. The specific surface area of the material was determined through the adsorption of N2 at 350 ± 5°C in a BELSORP-MINI equipment, from Bel Japan Inc. Initially, 200 mg of the sample was pre-treated at 333 K under vacuum conditions (DEGASS) for 3 h.
X-ray diffraction. The charcoal was characterised in a BRUKER X-ray diffractometer (model D8 Advance), using a Cu-Kα radiation source, at a voltage of 30 kV and a current of 30 mA. The data were collected at a 2θ range from 5° to 80° with a spacing of 0.05° and a time of 2.0 s.
Fourier transform infrared spectroscopy. The absorption spectra were obtained from a BRUKER spectrometer (model VERTEX 70), using the attenuated total reflection technique. The data were collected in the infrared region between 4000 and 400 cm−1, with a spectral resolution of 2 cm−1.
Scanning electron microscopy with energy dispersive X-ray spectroscopy. The images were obtained using energy dispersive spectrophotometer connected to a SHIMADZU scanning electron microscope (Superscan SS-550 model) at 15 K. The charcoal was placed on a carbon tape and, subsequently, a thin layer of gold was added to further improve the sample handling.
Factorial design
The factorial design was used to evaluate the influences of the variables particle size (<0.09, 0.09–0.15 and 0.15–0.21 mm) and stirring speed (0, 150 and 300 r min−1) on the adsorption process. The trials were carried out in a random order, and the factorial design used was a 22 type with a central point run performed in triplicate, in order to guarantee the required repeatability of the experimental data. The response used to determine the efficiency of the process was the adsorptive capacity, q (g g−1).
The values of q were obtained from equation (1)
The 0.05 g of Cc that remained in contact with the 5 ml of the AVK model solution was used at 1% (v/v) during 360 min at a temperature of 28 ± 2°C, on a shaker (IKA, KS 130 control), according to Nascimento et al. (2014).
The calculations of the effects of the factors and interactions among them, as well as the respective standard errors, were made according to Barros Neto et al. (2007) using the Statistica 6.0 programme.
Kinetic adsorption study
Based on the conditions established by the factorial design, a kinetic study was carried out aiming at determining the system’s equilibrium time. The kinetic trials were carried out by adding the 0.05 g of Cc to the 1% (v/v) AVK model solution (5 ml), with the samples being filtered at time intervals from 5 to 540 min, at a temperature of 28 ± 2°C.
The kinetic pseudo-first-order (equation (2)) and pseudo-second-order (equation (3)) kinetic models, the Weber–Morris model (equation (4)) and the Boyd model (equation (5)) were adjusted to the experimental data obtained using a non-linear regression model with the Origin 8.0 computer programme. The models that best adjusted were compared using an F-test, according to Montgomery (2012)
Adsorption equilibrium study
Based on the conditions established by the factorial design, an equilibrium study was developed aimed at determining the maximum adsorptive capacities for the Cc charcoal.
It was utilised, in this study, 5 ml AVK model solution with concentration ranging from 0.5 to 5%, being in contact with 0.05 g of the adsorbent under constant stirring at 300 r min−1 at 28 ± 2° C. After 360 min, the samples were filtered and quantified.
The Langmuir (equation (6)), Freundlich (equation (7)), Fritz–Schlünder (equation (8)) and BET (equation (9)) adsorption models were adjusted to the experimental data by using the MATLAB version 7.8.0 programme, in order to evaluate the equilibrium data of the process
Results and discussion
Characterisation of the carbonaceous materials
The Cc charcoal exhibited a pHpzc equal to 10.5. Activated charcoals physically at high temperatures are generally characterised by their basic surface, according to Matos et al. (2011). A similar result was found by Largitte et al. (2016) who obtained a pHpzc of 9.6 and 10.3 for the activated charcoal made of guava seeds and almond shells, activated with CO2, respectively.
According to Ould-Idriss et al. (2011) below the pHpzc the material exhibits positive charges on its surface, favouring the adsorption of anions. Consequently, above the pHpzc, the material has negative charges on its surface, thus favouring the adsorption of cations. As the pHpzc > pH of the model mixture of AVK (4.1), there is a predominance of particles positively charged on the surface of the charcoal, favouring the adsorption of naphthenic acid.
The results obtained in the textural characterisation by adsorption/desorption of N2 for the Cc were surface area of 243.2 m2 g−1, pore volume of 0.13 cm3 g−1 and a pore diameter of 22 Å. The charcoal studied exhibited a pore diameter within the range of 20–500 Å, being classified as mesoporous, according to the International Union of Pure and Applied Chemistry (Thommes et al., 2015). The surface area found was within the range of values reported for activated charcoal, according to Habila et al. (2015) and Rodrigues et al. (2011).
In the X-ray diffraction for the Cc, peaks overlapping the two halos (approximately 24° and 43°) can be noticed in the diffractogram; such peaks demonstrate the presence of crystalline particles in the materials, which consist of residues produced in the carbonisation process. The halos in such positions are characteristic of the diffraction profile of carbon graphite and are attributed to the reflections at both planes (002) and (101). The latter represents the overlapping of reflections at the 100 and 101 planes. Similar results were obtained by Huang et al. (2011) for activated charcoal of stalks lotus.
The adsorption spectra in the infrared region, with Fourier transform, of the material studied are exhibited in Figure 1.

Infrared spectrum for the Cc.
For the Cc, the 3725–2950 cm−1 band was obtained indicating the presence of functional groups associated with the phenolic, alcohol and carboxylic acids O–H. The band at 1640–1510 cm−1 can be attributed to the oxygen functional groups, such as the highly conjugated C = O bonds, the stretching of the C–O bonds, in carboxylic, aldehyde, ketone, ester and carboxylate salt groups (Kumar and Jena, 2015). The band ranging between 1510 and 1335 cm−1 can refer to the C–C stretching vibration in aromatic rings. Peaks at 877 and 750 cm−1 are indicative of the presence of C–H in aromatic ring (Huang et al., 2015).
The Cc presented a wrinkled, heterogeneous and irregular surface, as well as a wide range of particle sizes. Similar results were found by Habila et al. (2015), who prepared activated charcoal from agro-industrial solid residues, obtaining a wrinkled surface with similar pores and cavities.
Factorial design
Based on the experimental data, and using the Statistica 6.0 computer programme, the main effects and interactions, with their respective standard deviations, were calculated, and the Pareto chart and response surfaces of the model in study were drawn.
The effects that were considered as being statistically significant at a confidence level of 95% on the levels studied are those placed after the line (p = 0.05) of the Pareto chart in Figure 2(a). The results of the effects of the factors and the interactions between them demonstrated that the effect Stirring Speed and the interaction between the effects Particle Size and Stirring Speed were statistically significant. The effect of the interaction between both effects can be better visualised through the surface presented in Figure 2(b), in which a higher adsorption capacity can be observed for the level < 0.09 mm and 300 r min−1.

(a) Pareto chart for the effects calculated and (b) response surface for the adsorption capacity (pure error: 0.00005).
Kinetic study of adsorption
The adsorption kinetics of n-dodecanoic carboxylic acid for the solid phase of the Cc as an adsorbent and the adjustment to the non-linear pseudo-first- and pseudo-second-order models are shown in Figure 3.

Adsorption kinetics in solid phase using a finite bath system and the adjustment to the pseudo-first- and pseudo-second-order models.
By observing Figure 3, the adsorption capacity, q, increased with the increase of contact time with the Cc, and it can be observed that the removal mostly took place at around 180 min, reaching the equilibrium after 360 min, with the experimental adsorption capacity (qexp) being equal to 0.57 g g−1.
This trend is caused due to the fact that, at the start, there is a great number of empty sites available for adsorption, which, with time, decreases; other than that, the repulsion forces of the adsorbed molecules tend to inhibit the adsorption at other sites.
Comparing these models using the F-test, it was found that Fcal (1) <Ftab (2.42); thus, both models did not exhibit any significant statistical difference at a confidence level of 95%.
The parameters of the kinetic models calculated are presented in Table 1.
Kinetic parameters calculated for the pseudo-first- and pseudo-second-order models.
Although by using the F-test, according to Montgomery (2012), no significant difference is found between the models evaluated, the qexp was similar to the qe for the pseudo-first-order model, but qe obtained through the pseudo-second-order model is not in agreement with the qexp. Therefore, it may be suggested that the adsorption of naphthenic acid on activated carbon from Umbaúba wood has a behaviour of first-order reaction according to Ho and Mckay (1999).
The adjustment of the experimental data to the Weber–Morris intraparticle diffusion model is presented in Figure 4.

Weber–Morris intraparticle diffusion model for the adsorption of naphthenic acid by Cc.
According to this model, when the graph of qt in relation to t1/2 results in a straight line passing through the origin, the rate of the sorption process can be considered as being limited by the intraparticle diffusion process. However, as the results obtained were multilinear (Figure 4), the adsorption process was considered as not being controlled only by the intraparticle diffusion step, with two or more steps controlling the process.
The data are represented by three linear phases, with the initial phase (kdif1) representing the effect of the boundary layer, with external mass transfer, in which the n-dodecanoic carboxylic acid was rapidly adsorbed by the Cc. After 60 min, the adsorption rate was reduced, resulting in the second phase, which extended until 300 min. According to Kumar and Porkodi (2007), this phase refers to the diffusion of molecules to the most internal adsorption sites of the adsorbent (kdif2).
The equilibrium was reached after 360 min, leading to the decline of intraparticle diffusion, due to less availability of adsorption sites (kdif3). Comparing the three diffusion constants of the Weber–Morris model (kdif), it is concluded that kdif1>kdif2>kdif3, as observed in Table 2. The molecules diffusion in the interior of the adsorbent is the determining step of the adsorption process, being confirmed by the decrease of kdif3 value.
Intraparticle diffusion constants of the Weber–Morris model.
According to Bhattacharyya and Sharma (2005), this is another indication that the intraparticle diffusion mechanism is not the only limiting step of the adsorption process, and that other interaction mechanisms (adsorption on the external surface and diffusion into the interior) might be acting simultaneously to it.
The Boyd model assumes, however, that the diffusion is the limiting step of the adsorption process, due to the boundary layer that involves the adsorbent particle. Aiming at ratifying the results obtained by the Weber–Morris model, the experimental data were adjusted to the Boyd kinetic expression (Figure 5).

Boyd kinetic model for the adsorption of naphthenic acid by Cc.
It was therefore noticed that the experimental data did not follow a linear trend, not passing through the origin. This demonstrates that the adsorption process is controlled both by the effects of intraparticle diffusion and external diffusion, corroborating with the results obtained by the Weber and Morris model, presented previously.
Equilibrium study of adsorption
The adsorption isotherm of the Cc and the adjustment to the experimental data to the BET equation for types IV and V isotherms are presented in Figure 6.

Equilibrium adsorption curve of naphthenic acid in the AVK model using Cc.
By analysing Figure 6, it can be observed that the Cc exhibited a maximum value of qexp equal to 1.1 g g−1. A similar result was obtained by Silva (2007), who obtained q equal to 1.48 g g−1 for the adsorption of naphthenic acid using an MgMCM-41 molecular sieve.
The experimental data applied to the Langmuir, Freundlich and Fritz–Schlünder models are presented in Table 3.
Parameters of the Langmuir, Freundlich and Fritz–Schlünder models.
It is therefore concluded that the Langmuir, Freundlich and Fritz–Schlünder models did not adjust to the experimental data, with r < 0.9; with only the BET model having adjusted to the equilibrium data in the adsorption of n-dodecanoic acid in AVK using Cc. Since the obtained isotherm (Figure 6) was type IV, the adsorption of n-dodecanoic acid using CC occurred in multilayers. According to Do (1998), the BET isotherm model assumes the same hypotheses presented by the Langmuir mechanism. However, while the Langmuir model only describes the adsorption in a single layer, the BET one can describe the occurrence of a multilayer adsorption. A comparison between the experimental data and the data obtained from the adjustment curve for the adsorption isotherm of Cc is shown in Figure 7.

Comparison between the experimental and calculated equilibrium data of the adsorption of naphthenic acid in AVK using Cc.
The value of the correlation coefficient (r) of 0.9598 demonstrates an agreement between the experimental and calculated data.
The parameters of this model were obtained from the adjustment of the equilibrium curve of the adsorbent material and are presented in Table 4.
Parameters of the BET model for the types IV and V isotherms for the adsorption of n-dodecanoic acid by Cc.
BET: Brunauer–Emmett–Teller.
The comparison between the parameters obtained in this work for the BET model with the ones found by using other adsorbents found in the literature for the adsorption of n-dodecanoic acid in an AVK model mixture is presented in Table 4.
According to Table 4, the material studied exhibited a higher Cs in equilibrium with qs, with that being an indication of higher adsorption on the first layer if compared to the remaining adsorbents, despite the lower number of layers.
Regarding the affinity between the adsorbate and the adsorbent surface (b), the material studied had a similar value to the MCM-41 molecular sieve – a remarkable result, given it is an activated charcoal. The activated charcoal exhibited more than three adsorption layers (nc) due to the heterogeneity of the process.
The extra energy required to dissolve an adsorbate (Q) for the activated charcoals MgMCM-41, MCM-41 and Sr-MCM-41 was calculated from the values of g presented in Table 4, being equal to 20.1, 31.5, 11.5 and 19.6 kJ mol−1, respectively. The higher the value of Q, the higher was the affinity between the adsorbent and adsorbate. The Cc had a higher energy for a higher desorption, if compared to both MCM-41 and Sr-MCM-41 indicating that the adsorbent studied had an affinity for naphthenic acid.
Conclusion
The activated charcoal produced from the Umbaúba wood is a mesoporous adsorbent, with a heterogeneous surface, and presented a wide range of positive charges on its surface, which favours the adsorption of naphthenic acid.
The best condition for the adsorption of naphthenic acid in the Cc was obtained for a particle size < 0.09 mm and for 300 r min−1, at the levels studied. The equilibrium time of the adsorption process was of 180 min, and there was no significant difference between the pseudo-first- and pseudo-second-order models. The diffusion was considered as not being the limiting step of the adsorption process studied, according to the Weber–Morris and Boyd models. An adsorption capacity of 1.1 g g−1 was obtained for the adsorption equilibrium, with only the BET being adjusted to the experimental data.
The adsorption study of naphthenic acids in an AVK model mixture using activated charcoal made from Umbaúba was considered as being significantly efficient. It is noteworthy that even though the charcoal studied was prepared from biomass, it presented an adsorption capacity similar to the ones found in the literature.
Footnotes
Acknowledgements
The authors are also thankful to Elephant Chemical Industry Ltda. for donating the charcoal. The authors are also thankful to CETENE, LACOM-UFPB and Analytical Center/DQF- UFPE for analysing the results of characterisation.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors would like to thank PETROBRAS for providing the financial resources needed for this work.
