Abstract

In this research, the adsorption potential of quaternized palm kernel shell (QPKS) to remove F from aqueous solution was investigated using fixed-bed adsorption column. Raw palm kernel shell waste was reacted with 3-chloro-2-hydroxypropyl trimethylammonium chloride (CHMAC) in order to modify the surface charge. The effects of inlet F concentrations (2–12 mg/l) and QPKS bed height (2–10 cm) with optimum pH (pH = 3) on the breakthrough characteristics of the adsorption system were determined. In the fixed-bed column, breakthrough time increases with increasing bed height due to increasing amount of active site on adsorbents to adsorb the fluoride ion. Decreasing trend of breakthrough values was obtained with increasing initial fluoride concentration due to greater driving force for the transfer process to overcome the mass transfer resistance in the column. The adsorptions were fitted to three well-established fixed-bed adsorption models, namely, Thomas, Yoon–Nelson, and Adams–Bohart models. The results fitted well to the Thomas and Yoon–Nelson models with correlation coefficient, R2 ≥ 0.96.

1. Introduction

Today, hazardous waste has become one of the environmental problems. Hazardous waste can be classified into different types of waste mainly nuclear waste, industrial waste, universal waste, medical waste, and construction waste. Hazardous waste can be in the form of solid, liquid, and gas which are not easily disposed due to its corrosiveness, toxicity, and ignitability. Fluoride is one of the hazardous wastes commonly produced by industries using hydrofluoric acid as cleaning agents or etchant. Industries manufacturing glass and ceramic, semiconductors, electroplating, coal-fired power stations, beryllium extraction plants, brick and iron works, and aluminium smelters discharge high concentration of fluoride [1]. High concentration of fluoride ions in water is mainly affected by natural minerals and industrial activities. Lv et al. [2] reported that concentration of fluoride in the wastewater from Zhejiang Juhua Fluorine Chemical Co. Ltd. exceeds 1000 mg/L (pH = 5–7). According to Ruan et al. [3], most natural rocks and minerals contain large quantities of fluoride which are released as fluoride ions into water contributing to levels of fluoride concentrations rise in water that exceed the acceptable limit set by the WHO which is 1.5 mg/L. Therefore, fluoride content in wastewater and drinking water should be reduced in order to improve the health quality.

The standards prescribed by various regulatory bodies for fluoride concentration in drinking water are different according to their climatic conditions. According to the WHO, the standard prescribed for fluoride ion concentration in drinking water is 1.5 mg/L. Fluoride in smaller doses (0.8–1.0 mg/L) helps to prevent dental caries particularly among children below eight years of age. Fluoride in higher concentration causes dental fluorosis (1.5–2.0 mg/L) and skeletal fluorosis (>3.0 mg/L) [4]. Chemically, fluorine is the most electronegative element and is always present in a combined state as fluoride due to its high chemical reactivity. Fluoride is a great calcium-seeking element and can affect the calcified structure of bones and teeth in the human body at higher concentration, resulting in dental fluorosis or skeletal fluorosis. Fluoride toxicity can also cause nonskeletal diseases such as aches and pain in the joints, nonulcer dyspepsia, polyurea (tendency to urinate more frequently) and polydipsia (excessive thrust), muscle weakness, fatigue, and anemia with very low haemoglobin levels [5].

Over the years, many techniques and processes have been developed and used in the treatment and removal of fluoride from the contaminated water effluent. These techniques include coagulation and flocculation, chemical precipitation, the use of membranes for membrane separation, aerobic and anaerobic degradation using various microorganisms, chemical oxidation, ion exchange, electrodialysis, reverse osmosis, foam flotation, electrolysis, and adsorption [6]. Some of these techniques have been proven to be effective, although they displayed some limitations such as the excess amount of chemical usage or accumulation of highly concentrated sludge, causing disposal problems, lack of effective fluoride reduction, and sensitivity towards various wastewater inputs [7].

Amongst the processes mentioned above, the use of adsorption has been the most prominent and widely used because of its cost-effectiveness, efficiency, and technology readiness due to the fact that it produces effluents containing very low levels of dissolved organic compounds [8]. In addition, many adsorption techniques have been employed for the treatment of drinking water, as reviewed in the work reported by Tikki [9].

Adsorption is the process through which a substance, originally present in one phase, is removed from that phase of accumulation at the interface between that phase and a separate (solid) phase. It can also be referred to as the separation of components in a fluid mixture by the transfer of one or more components (known as the adsorbate) to the internal surface of a porous solid (known as the adsorbent) where they are held by intermolecular forces. However, adsorption is different from absorption in the sense that, in absorption, a substance diffuses into a liquid or solid to form a solution. Thus, the connecting term “sorption” encompasses both processes, whereas desorption is the reverse process.

Generally, the solid material that adsorbs the adsorbate is referred to the adsorbent. There are various types and classifications of adsorbents; but in general, adsorbents are classified to be either organic or inorganic adsorbents. The physical properties of adsorbents include their shape, porous structure, and specific surface area. In order to determine the type of adsorbent to be selected for a particular adsorption process, the following criteria must be taken into consideration:(a)High capacity to reduce the amount of adsorbent needed(b)Favourable kinetic and transport properties to allow for rapid sorption(c)Should be thermally and chemically stable in order to allow for the preservation of both quantity and properties(d)Good hardness and material strength(e)Resistance to fouling(f)Regenerative capability to allow for reuse [10]

The adsorption process is based on three steps as explained in the following by Viegas et al. [11]:(1)Mass transfer of solute from the bulk solution through the stagnant film surroundings to the particle external surface (external or film mass transport).(2)Mass transfer of solute within the particle (internal or intraparticle diffusion). Intraparticle diffusion may be due to pore volume diffusion (diffusion in the fluid-filled pores), surface diffusion (migration along the pore surface in which an adsorbate hops from one to another available adsorption site in a series of adsorption-desorption reactions), or a combination of both.(3)Solute attachment onto the surface of the adsorbent surface site, adsorption.

Step (3), adsorption, is usually very rapid compared to the first two steps, and the overall rate of adsorption is therefore often controlled by the first or the second steps, whichever is slower, or a combination of both.

Adsorption process is widely used due to its effectiveness giving satisfactory results and also in terms of cost, simplicity of design, and operations [12, 13]. Many research studies have focused on the low-cost adsorbent in the removal of wastewater pollutant such as heavy metal by carbonized green mussel shell [14], modified agriculture wastes [15], aerogel adsorbent from combination of waste office paper and chitosan [16], and many more. Meanwhile, low-cost adsorbent application in fluoride removal has been studied by Singh et al. [17] who used sugarcane bagasse as an adsorbent, Ayu Haslija et al. [18] used modified palm kernel shell in batch process adsorption, Mukherjee et al. [19] synthesized zeolite adsorbent from rice husk, and Hu et al. [20] used chitosan/aluminum hydroxide bead adsorbent, and many more.

Malaysia is one of the largest producers and exporters of palm oil in the world producing a large quantity of by-products which is the palm kernel shell. According to the data obtained from the Malaysian Palm Oil Board (MPOB) which is represented in Table 1, 4.5 million tonnes of palm kernel shell (PKS) has been produced as of November 2017. This big amount could lead to a serious environmental problem if no proper waste management system is carried out.

The shell is practically removed by steam treatment and mechanically crushed from the oil palm fruit. These by-products are usually burned off emitting dark smoke as well as partially carbonized fibrous particulates due to incomplete combustion leading to pollutions. Thus, to prevent the pollutions, the PKS can be further processed as the adsorbent for adsorption process.

In this research, palm kernel shell is chemically modified to act as an adsorbent for F removal from aqueous solution by using the fixed-bed column. Characterization analysis was done to analyze the properties of the adsorbent produced. The crucial design parameters such as inlet concentration of F solution and column bed height were investigated. The breakthrough curves for the adsorption of F were analyzed using Adams–Bohart, Thomas, and Yoon–Nelson models.

2. Materials and Methods

2.1. Materials

Palm kernel shell was obtained from Seri Ulu Langat Palm Oil Mill Sdn. Bhd., Selangor, Malaysia. 60 wt. % of sodium hydroxide (NaOH) was used for the mercerization process of palm kernel shell while 3-chloro-2-hydroxypropyl trimethylammonium chloride (CHMAC) solution used for the quaternization process. Stock solutions of F were prepared by dissolving appropriate amount of NaF in the deionized distilled water.

2.2. Quaternization and Experimental Setup

The quaternization process was carried out by impregnating mercerized PKS with NaOH solution and CHMAC at a ratio (w/w) of 5 H2O : 1 sample : 1 NaOH : 4 CHMAC [22]. All chemicals were purchased from Sigma-Aldrich. The sample was then shaken at 100 rpm on a digital orbital shaker, Wise Shake SHO-2D (Wise Laboratory Instruments, Korea), for 24 hours to accelerate the reaction. After quaternization, the quaternized PKS (QPKS) was washed with 2% acetic acid solution, followed by distilled water until the solution was neutral (pH 7) to stop the reaction. The washed QPKS was dried in the oven at 60°C overnight and cooled in a desiccator for further experimental usage.

A transparent cylindrical glass column of 2 cm internal diameter and 35 cm height, as shown by Figure 1, was used to conduct the continuous column adsorption experiments. A glass wool of 0.5 mm was attached at the bottom of the column. The column performance of F adsorption onto QPKS was studied at different F concentration (2–12 mg/l) and bed height (2–10 cm) at pH 3. The pH was adjusted by adding 0.1 M HCl into the solution. Effluents were collected at the bottom of the column at regular intervals and analyzed for fluoride using UV/VIS double beam spectrophotometer (GENESYS 10 UV, USA) at a wavelength of 570 nm.

2.3. Adsorbent Analysis
2.3.1. Functional Group Analysis

Functional groups of samples and infrared (IR) spectra were determined and recorded on the Fourier-transform infrared (FTIR) spectrometer (PerkinElmer 1750X, PerkinElmer, USA). Attenuated total reflectance (ATR) technique was used for all sample measurements whereby spectra were collected based on measurements of the temporal coherence of a radiative source. The percentage of transmittance against wave number (cm−1) within the range of 400–4000 cm−1 at room temperature was plotted for all spectra.

2.3.2. Surface Morphology and Elemental Analysis

Surface morphology and elements in of raw palm kernel shell (PKS) were analyzed and compared with quaternized PKS and PKSAC by using scanning electron microscopy with energy-dispersive X-ray (Hitachi Model S-3400N, Hitachi, Japan). SEM-EDX is a microscope that uses electrons to form images.

2.3.3. Porosity and Surface Area Analysis

Brunauer–Emmett–Teller (BET) method was used to analyze the surface area of samples. The BET surface area analyser used was AutoSorb-1 (Quantachrome Instruments, USA). Pore size distribution, average pore size, and pore volume of the samples were determined using the Barrett–Johner–Halenda (BJH) equation.

2.3.4. Crystallinity Analysis

XRD analysis was performed by PANalytical XRD instrument with Cu radiation 1.542 Å in order to identify the adsorbents’ crystalline materials (e.g., minerals and inorganic compounds). The X-ray tube was operated at 45 kV and 40 mA where the 2θ angle of the diffractometer was stepped from 10° to 80° scan rate with a step size of 0.02°.

2.4. Breakthrough Profile

Breakthrough profiles were plotted for each parameter to analyze the best condition to efficiently remove fluoride by plotting Ct/C0 against time, t (min). The maximum column capacity, q (mg/g), for a given parameter is equal to the area under the breakthrough curve.

2.5. Kinetic Models

In this research, Thomas, Yoon–Nelson, and Adams–Bohart models were used to analyze the kinetic studies. The Thomas model is based on the assumption that the process follows Langmuir kinetics of adsorption-desorption with no axial dispersion. It describes that the rate driving force obeys the second-order reversible reaction kinetics [23].

The linearized form of the model is given as follows:where is the Thomas rate constant (mL/min·mg), is the equilibrium adsorbate uptake per g of adsorbent, is the inlet adsorbate concentration (mg/L), is the outlet concentration at time t (mg/L), is the mass of adsorbent (g), is the filtration velocity (mL/min), and is the filtration time (min).

The Adams–Bohart model is used for the description of the initial part of a breakthrough curve but its validity is limited in the concentration of C < 0.5C0 range [24]. The mathematical equation of the model can be written as follows:where (cm/min) is the linear flow rate, (mg/L) is the saturation concentration, (L/mg·min) is the kinetic constant, and is the bed depth of column (cm). The linear flow rate was calculated from the equation where is the volumetric flow rate (ml/min) and is the transversal area of the column (cm2).

A prediction model bed depth service time (BDST) which focuses on the theory of surface reaction rate and prediction of column efficiency between bed height (Z) and service time (t) is based on the Adams–Bohart model. This BDST model can be used to achieve the required breakthrough at the given column condition [25]. The Adams–Bohart equation was modified by Hutchins [26] who proposed a linear relationship between bed height and service time aswhere is the initial F concentration (mg/L), is the breakthrough F ion concentration (mg/L), is the linear velocity (cm/min), is the adsorption capacity of bed (mg/L), is the rate constant of the BDST model (L/mg/min), is the time (min), and is the bed height (cm) of the column. Equation (3) can be linearized into equation (4):where the critical bed depth, , occurs at time which is the theoretical value of bed depth of the adsorbent which the effluent concentration does not exceed the breakthrough concentration (Cb). Therefore, when , equation (3) changes to equation (5):

According to Mohan et al. [25], a prediction of the concentration-time profile or breakthrough curve and adsorption capacity for the effluent under given specific operating conditions is required in order to design a column for the adsorption process.

Yoon and Nelson [27] developed a model known as Yoon–Nelson model, which assumed that the rate of decrease in the probability of sorption for each adsorbate molecule was proportional to the probability of adsorbate sorption and the probability of adsorbate breakthrough on the adsorbent [28]. The model is described by the following equation:where (min−1) is the Yoon–Nelson rate constant and (min) is the time required for 50% sorbate breakthrough.

Tofan et al. [23] stated that, according to the Yoon–Nelson model, the amount of adsorbate ion sorbed in a fixed bed is half of the total adsorbate ion entering the adsorption bed within 2τ period. Results obtained for adsorption of Cr (III) by modified hemp fibers showed that the values of the Yoon–Nelson rate constant, , increase as the initial concentration increases. This is probably due to the increases of competition between the sorbate species on the sorption sites by increasing initial concentration of the metal ion, which results in a higher rate of retention. It is known that and are inversely related; therefore, the time required for 50% breakthrough, , decreases with increasing the influent concentration of Cr(III).

3. Results and Discussion

3.1. Adsorbent Characteristic
3.1.1. Surface Morphology Analysis

Surface morphology of raw PKS and QPKS was analyzed using SEM. Figure 2(a) presents the surface of raw PKS which showed small size, uneven, and low amount of pores while the surface seems to be covered by a layer of impurities. On the other hand, Figure 2(b) shows that QPKS had more pores compared to raw PKS (Figure 2(a)) which resulted from the mercerization and quaternization process which remove impurities under alkaline condition in NaOH [29] and therefore created more pores. The quantitative data on PKS and QPKS pore size were determined using BET analysis in the subsequent part.

3.1.2. Surface Elemental Analysis

Energy-dispersive X-ray analysis is a surface elemental analysis used in conjunction with SEM. The results for PKS and QPKS are listed in Table 2.

A high percentage of nitrogen was detected in QPKS as the result of quaternization process, but was not detected in PKS. Hence, this analysis again verifies that the quaternary ammonia group () reacted on the surface of QPKS, resulting in a successful synthesis.

3.1.3. Porosity and Surface Area Analysis

The adsorbent’s pore characteristics are given in Table 3, which show that the average pore diameters of PKS and QPKS were 2.2 nm and 35.2 nm, respectively. Due to the size of these pores, they are classified as mesopores. Enlargement of the pore size in QPKS was due to the dissolution of impurities in NaOH solution during the pretreatment mercerization and quaternization processes.

3.1.4. Functional Group Analysis

The fundamental of FTIR analysis is to identify the functional groups exist in an organic compound. Figure 3 shows FTIR spectra of raw PKS and QPKS, where the existence of –CN, –COC (ether), and –Cl was shown in the peaks at wavenumbers of 1331 cm−1, 1105 cm−1, and 682 cm−1, respectively, which were not present in raw PKS. In order to prove the designated peaks such as –CN, –COC, and –Cl peaks, software from Bio-Rad, KnowItAll-Academic Edition, has been applied. Figures 3(b)–3(d) show the peaks extracted from software. The existence of the –CN peak showed that the surface of PKS was successfully quaternized () with CHMAC, and this is one of the main factors to enhance the adsorption of fluoride ions.

These peaks can attribute quaternization mechanism process as illustrated in Figure 4 adapted from de Lima et al. There are three reaction steps in a quaternization reaction: the interaction between cellulose and sodium hydroxide, epoxy formation from CHMAC, and the reaction between epoxide and the hydroxyl group of cellulose to produce ether with the presence of acid to stop the reaction, as shown in Figure 4.

3.1.5. Crystallinity Analysis

Pattern of XRD peaks for raw PKS and QPKS is shown in Figure 5. Figure 5(a) shows raw PKS’ peak before modification in which two main peaks can be observed at 2θ = 22° and 35° indicating the existence of carbon and oxygen. The peak of QPKS is represented by Figure 5(b) which shows different pattern due to the quaternization process. The highest peak of QPKS is at 2θ = 12° which is not present in raw PKS including peak at 2θ = 24°, 26°, and 33°. These peaks indicate the presence of nitrogen and chloride element in the QPKS base on reference database in Match software. Figure 5(c) is resulted after removing raw PKS peak from QPKS peak and clearly proves the existence of nitrogen and chloride on QPKS due to the quaternization process. This characteristic is the main factor on the adsorption of fluoride from solution.

3.1.6. Point of Zero Charge

The point of zero charge was performed to determine the surface chemistry characterisation of adsorbent. By conducting this test, the pH of the point of zero charge, pHpzc, of QPKS was determined. The pHpzc of the adsorbent is a very important characteristic that determines the pH when the electrical charge density on a surface is zero. From Figure 6, the pHpzc of QPKS, which was the point when ∆pH = 0 (pHinitial = pHfinal) was close to pH = 4.1. At ∆pH < 0, there was a net positively charge on the surface of adsorbents, which was likely to attract anions. The most negative value of ∆pH for QPKS was found at pH 3.2 (−1.01), which indicated that the surface of the adsorbent was most positively charged at this pH and thus, exerted the greatest electrostatic forces of attraction for anions.

According to Railsback [31], pH below pHpzc had positively charged surface due to H+ ions being adsorbed on the surface rather than entering the low pH (H+ rich) solution, resulting in a positively charged (H+ rich) surface, which is a very important characteristic in adsorbing fluoride ion (F). In this study, the electrolyte solution used was 0.1 M NaNO3, which consisted of Na+ as cations and as anions, which may be adsorbed. The following equations derived from Cardenas-Peña et al. [32] showed Q+, which represents the surface of QPKS. Hence, by combining with Railsback’s [31] statement, the following equations can best describe the mechanism of the adsorption process at optimum pH:

Meanwhile, in the condition of high pH, which contributed of more ions in the solution, hydrofluoric acid might be produced. On the other hand, displacement of ion from the quaternized adsorbent is replaced by F producing HCl in the solution. The following equation shows the possibility of the reactions to occur:

3.2. Effect of Inlet Initial F Concentration

Effect of initial fluoride concentration on breakthrough profile at 2, 8, and 12 mg/L with optimum pH, bed height of 6 cm, and flow rate at 4 ml/min was studied, and the results are shown in Figure 7 and Table 4. The breakthrough curve clearly showed the transfer mechanism at different initial fluoride concentrations.

Using 2 mg/l as the lowest initial concentration, breakthrough and exhausted time for each adsorbent were analyzed. The curve from Figure 7 clearly shows that QPKS achieved breakthrough time at 17.5 minute and exhausted at 23.75 minute. At initial fluoride concentration of 8 mg/L, the breakthrough time and exhausted time of QPKS were at 9.3 and 15 minutes. Meanwhile, for 12 mg/l of fluoride concentration, QPKS breakthrough at a faster rate was 5.3 minutes due to the vacant site filled up faster by fluoride ion and exhausted at 8.75 minutes.

From the analysis, it can be concluded that breakthrough time decreased with increased initial fluoride concentration. As the concentration increased, amount of adsorbate also increased; therefore, active sites were occupied at a faster rate. Adsorbate ions will saturate on adsorbent faster (exhaust) and in a shorter time to achieve the breakthrough time with increasing influent fluoride concentration due to greater mass transfer driving force for the transfer process to overcome the mass transfer resistance in the column [3335]. Similar trend was reported by Chen et al. [33] for fluoride removal using Kanuma mud and hydrous ferric oxide by Nur et al. [34]. Table 1 summarizes all the breakthrough and exhaust time for QPKS.

3.3. Effects of Bed Height

Figure 8 shows the breakthrough profile of bed heights at different bed heights. The breakthrough curve for 2 cm bed height had steeper curve as compared to curve at 6 and 10 cm bed height. The curve trend of bed height effect was observed to decrease with higher bed height but the higher the bed height, the longer the breakthrough and exhaust time. Breakthrough occurred when Ct/C0 was equal to 0.75. Similar trend was observed in fluoride adsorption on manganese oxide coated alumina [36] and on hydrous ferric oxide (HFO) [34]. The area above the breakthrough curve also represented the adsorption capacity of adsorbent [37]. Exhaustion occurred when all active sites on adsorbents were occupied or saturated by adsorbate ions which in this study were the fluoride ions. At this point, final concentration was the same as initial concentration.

At a bed height of 2 cm, the breakthrough time and exhausted time of QPKS were 8.75 minutes and 17.5 minutes. This is due to the porosity effect of each adsorbent in which QPKS has larger pore size as compared to PKSAC which resulted in lower amount of active site on PKSAC to adsorb the fluoride ion. Efficiency and capacity of column adsorption was largely dependent on the amount of adsorbent in the column and porosity. Table 5 shows the breakthrough and exhaustion time for all bed heights.

3.4. Kinetic Models
3.4.1. Yoon–Nelson Model

Figure 9 shows the Yoon–Nelson plot on fluoride adsorption at different bed heights and initial concentrations. Values for KYN were determined from the slope, and τ was calculated from the intercept. All the values obtained are compiled in Table 6. KYN values were found to be decreased with the increase of the bed height, and τ values, which are the time required for breakthrough, showed an increasing trend with increasing bed height. It indicated the adsorbate breakthrough was taking longer duration with increasing bed height. These findings were in agreement with the observations reported by Chowdhury et al. [38] and Yadav et al. [39]. Meanwhile, values of KYN were observed to be increased, and τ values decreased with increasing initial fluoride concentration. Chowdhury et al. [38], Nwabanne and Igbokwe [40], and Mohan et al. [25] also reported a similar trend in their work.

3.4.2. Thomas Model

The Thomas model is one of the most general and widely used models to describe the performance theory of the adsorption process in fixed-bed columns. The Thomas model assumed the adsorption process follows Langmuir kinetic. Using equation (1), KTh and qe were calculated from slope and intercepts of linear plots of versus t which are showed in Figure 10 for different bed heights and initial concentrations. Table 7 shows the results of effect of different bed heights and initial fluoride concentrations. KTh values decreased with increased parameter values but qe values showed increasing trend with increased parameter values. Higher mass transfer driving force was produced by increasing concentration of fluoride ions, therefore increases the concentration gradient and the adsorption capacity [35].

3.4.3. Adams–Bohart and Bed Depth Service Time (BDST) Model

The Adams–Bohart model constants, KAB and N0, were determined from the plot in Figure 11(a) at different bed heights and at different initial concentrations (Figure 11(b)). All the values are presented in Table 8. The values of R2 were between 0.59 and 0.98. Table 8 clearly shows that the values of KAB increased with increasing bed height; meanwhile, N0 decreased as the bed height increased. In contrast, initial fluoride concentration gave the different values for KAB and N0. As the initial concentration increased, KAB showed the decreasing values but an increasing trend for N0. The trend is in agreement with the findings reported by the study carried out by Ataei-Germi and Nematollahzadeh [41] who removed methylene blue using nanoparticles and Yadav et al. [39] on fixed-bed fluoride removal by using a commercial Bio-F adsorbent.

According to equation (4), t vs. Z was plotted and is shown in Figure 12 for each adsorbent, where t values were considered for breakthrough times and the exhaust times. The slope of the BDST line represents the time required for the adsorption zone to travel a unit length through the adsorbent under the selected experimental conditions at a given concentration. This can be used to predict the performance of the bed if there is a change in the initial solute concentration. The values of the constants are presented in Table 9, and zo which is the minimum column height to produce 1.5 mg/L fluoride was found to be 3.45 cm for QPKS.

Different adsorbents produced different adsorption capacities on removing fluoride due to their different properties. Table 10 shows the values of adsorption capacity by using different adsorbent from previous research findings.

4. Conclusions

From this study, QPKS prepared was found to be suitable for F removal from aqueous solution using continuous adsorption column. Breakthrough time increases with increasing bed height due to increasing amount of the active site on adsorbents to adsorb the fluoride ion. Decreasing trend of breakthrough values was obtained with increasing initial fluoride concentration due to greater driving force for the transfer process to overcome the mass transfer resistance in the column. Kinetic models plotted for fixed-bed column were Yoon–Nelson, Thomas, Adams–Bohart, and BDST models. Experimental data fit well with Yoon–Nelson with R2 values between 0.96 and 0.99, a range from 0.96 to 0.99 for the Thomas model, and 0.85 to 0.96 for the Adams–Bohart model. Thus, according to the performance, adsorption capacity, and sustainability on waste management, column adsorption using waste material as adsorbent is best to be used on removing fluoride in both drinking water and wastewater treatment.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors express their sincere gratitude to Chemical and Environmental Engineering Department of University Putra Malaysia and Faculty of Engineering, Technology and Built Environment, UCSI University, for the support given to this research project.