Artículo de investigación/Research article

Influencia de la altura del lecho sobre el comportamiento dinámico de columna de lecho fijo en la biosorción de mercurio

Influence of the bed height on the dynamic behavior of a fixed-bed column during mercury biosorption

Candelaria N. Tejada-Tovar¹, Angel Villabona-Ortíz², Gemaly Alvarez-Bajaire³, Leidy Jattin-Torres⁴ y Clemente Granados-Conde⁵

¹ M.Sc. en Ingeniería Ambiental. M.Sc. Educación, Grupo de Investigación en Diseño de procesos y Aprovechamiento de Biomasas IDAB, Facultad de Ingeniería, Universidad de Cartagena, Cartagena-Colombia, ctejadat@unicartagena.edu.co

² M.Sc. en Ingeniería Ambiental, Grupo de Investigación en Diseño de procesos y Aprovechamiento de Biomasas IDAB, Facultad de Ingeniería, Universidad de Cartagena, Cartagena-Colombia, avillabonao@unicartagena.edu.co

³ Ingeniera Química, Grupo de Investigación en Diseño de procesos y Aprovechamiento de Biomasas IDAB, Facultad de Ingeniería, Universidad de Cartagena, Cartagena-Colombia, gemalita92@hotmail.com

⁴ Ingeniera Química, Grupo de Investigación en Diseño de procesos y Aprovechamiento de Biomasas IDAB, Facultad de Ingeniería, Universidad de Cartagena, Cartagena-Colombia, jattincita_3004@hotmail.com

⁵ M.Sc. en Ciencia y Tecnología de los Alimentos, Grupo de investigación Ingeniería, Innovación, Calidad Alimentaria y Salud -INCAS, Facultad de Ingeniería, Universidad de Cartagena, Cartagena-Colombia, cgranadosc@unicartagena.edu.co

Fecha de recepción: 12 de junio de 2017/ Fecha de aceptación: 30 de agosto de 2017

Como citar / How to cite

C. N. Tejada-Tovar, A. Villabona-Ortíz, G. Alvarez-Bajaire, L. Jattin-Torres y C. Granados-Conde, Influencia de la altura del lecho sobre el comportamiento dinámico de columna de lecho fijo en la biosorción de mercurio. TecnoLógicas, vol. 20, no. 40, pp. 71-81, 2017.

Resumen

Se determinó el impacto de la altura del lecho en la remoción de Hg (II), presente en solución sobre biomasa residual de cáscara de cacao en sistema continuo de lecho fijo, evaluando el rendimiento de la densidad de empaque. El trabajo experimental estuvo fundamentado en: preparación de la biomasa, diseño y montaje de la unidad de adsorción y modelado matemático, considerando como variables intervinientes en el proceso la concentración inicial del metal, pH, velocidad de flujo y tamaño de partícula y como variable independiente incidente la densidad de empaque mg biomasa/volumen del lecho, traducido en altura en cm de empaque. Los análisis FTIR de la cáscara muestran la presencia de grupos funcionales que favorecen la adsorción de iones metálicos. La concentración de Hg (II) se midió por espectroscopia UV/vis; se obtuvo una capacidad máxima de adsorción fue de 99.62 % para el lecho de 10g (7.5 cm). Además, el modelo de Thomas fue el que mejor ajustó los datos experimentales. A partir de estos resultados, se concluye que la cáscara de cacao tiene potencial para ser utilizado como bioadsorbente de Hg (II) en solución y que el aumento de la altura del lecho en el sistema continuo favoreció la remoción del contaminante.

Palabras clave: Bioadsorbente, cáscara de cacao, sistema continuo, adsorción molecular, modelado matemático.

Abstract

This study evaluated the packing density of residual biomass of cocoa bean husk as adsorbent of Hg (II) dissolved in an aqueous solution inside a continuous fixed-bed system. The effects of the height of the bed on the removal of the pollutant were evaluated. This experimental work was based on biomass preparation, adsorber design and assembly, and mathematical modelling. The variables considered in the process were initial concentration of the metal, pH, flow rate and particle size. The incident-independent variable was the packing density (mg biomass/bed volume), which translated into the height (cm) of the packing. The FTIR analyses of the husk revealed the presence of functional groups in the spectrum that favour the adsorption of the metal. The residual concentration of the solution was measured by UV/Vis spectroscopy; the maximum adsorption capacity was 99.62%, by the 10g (7.5 cm) bed. In addition, the Thomas model was the best-fitting for the experimental data. On the basis of these results, we concluded that cocoa bean husk has potential to be used as bioadsorbent of Hg (II) from aqueous solutions and that the increase in bed height in the continuous system favours the removal of the pollutant.

Keywords: Bioadsorbent, cocoa shell, continuous system, molecular adsorption, mathematical modelling.

1. Introducción

Uno de los productos agroalimentarios de mayor inserción en el mercado internacional es el cacao (Theobroma cacao L.), lo cual ha generado que su procesamiento produzca grandes volúmenes de desechos [1] . Esto debido a que, durante la explotación cacaotera sólo se usufructúa la semilla, que se traduce en el 10 % de la masa total del producto. Los subproductos generados, cáscara y la pulpa, al ser depositados en los suelos se consideran un foco para la propagación de Moniliophthora roreri y Phytophthora spp., los cuales ocasionan mermas mercantiles, relacionadas con la descomposición de la mazorca originada por estos hongos [2] . De igual forma, en Colombia las enfermedades moniliasis y pudrición negra del fruto, han sido reportadas en plantíos cacaoteros; los resultados de esta investigación ayudaron a la comprensión del ciclo espacial en la Orinoquía Colombiana de estas enfermedades en tres sembradíos, ligados a la inadecuada disposición final de estos restos [3] .

Los vertimientos de aguas contaminadas con metales pesados a los cuerpos de agua superficiales, provenientes de actividades antrópicas e industriales como metal, textiles, curtiembres, minería, papel, refinerías de petróleo, y manufacturas químicas; representa una amenaza para la salud pública y ambiente, debido a la toxicidad por bioacumulación en los organismos de sales hidrosolubles de iones metálicos como el mercurio (Hg), contaminando así las cadenas tróficas, suscitando padecimientos como: ceguera, amnesia, raquitismo, miastenia, cáncer, insuficiencia renal, incluso la muerte [4] [5] [6] .

Diferentes tecnologías se han usado para la eliminación de metales presentes en efluentes de origen industrial, entre las que se cuentan: la precipitación, ultrafiltración, ósmosis inversa, intercambio iónico, extracción por membranas y solventes, electrodiálisis, entre otras [7] [8] [9]. Estos métodos pueden remover hasta el 99 % de estos metales, sin embargo, su alto costo de instalación y mantenimiento, la generación de lodos y su poca eficiencia a bajas concentraciones dificulta su aplicación a las pequeñas y medianas empresas [10]. Por lo tanto, se han utilizado procesos sostenibles para este fin, empleando residuos agrícolas como bioadsorbentes, los cuales son eficientes, económicos, abundantes y reutilizables [11] .

Adicionalmente, en la remoción de iones en sistema batch y continuo, se han usado como adsorbente carbones activados a partir de agujas de pino, cáscara de Swietenia mahagoni , nanotubos de carbono de paredes múltiples con grupos amino y tiolados, residuos de cacao, bagazo de palma, cáscaras de ñame, entre otros, demostrando el potencial de las biomasas lignocelulósicas en la remediación de aguas [4] [12-15]. En contexto, el presente trabajo evaluó el desempeño de la densidad de empaque de biomasa residual de cáscara de cacao como adsorbente de Hg (II), disuelto en solución acuosa, en un sistema continuo de lecho fijo determinando el efecto de la altura del lecho en la remoción de este contaminante.

2. Materiales y métodos

2.1 Preparación del bioadsorbente

La cáscara de cacao fue recolectada de cultivos cacaoteros del departamento de Bolívar en Colombia. Este residuo se lavó con agua destilada para remover las impurezas hasta que el agua no presentó coloración alguna. Terminada esta etapa, se secó el material en el horno a 90°C durante 24 h con el fin de retirar todo el exceso de agua del residuo. El material seco se llevó a molienda y se tamizó para obtener el tamaño de partícula seleccionado (0.5 mm) [1].

Se realizó el análisis fisicoquímico a la biomasa por el método EAA (Espectroscopia de Absorción Atómica) para determinar los elementos químicos: calcio, sodio, potasio, hierro, cobre y magnesio. Y un análisis elemental para determinar el porcentaje de carbono (método AOAC 949.14), hidrogeno (método AOAC 949.14), porcentaje de nitrógeno (método AOAC 984.13), contenido de azufre (digestión), porcentaje de cenizas y pectina (termogravimetría y digestión acida), porcentaje de lignina, celulosa y hemicelulosa (Fotocolorimetria y digestión). Análisis FTIR a la biomasa antes y después de la remoción, para identificar los grupos funcionales intervinientes en la adsorción [9].

2.3 Preparación del agua residual sintética

Para preparar el agua residual sintética se usó 0,135 g de cloruro de mercurio (II) (HgCl₂) disuelto en 1 L de agua desionizada obteniendo la concentración deseada de 100 ppm y se ajustó el pH a 6 usando NaOH [4].

2.4 Comportamiento dinámico de la columna de adsorción

Las pruebas de adsorción se realizaron en un equipo piloto que consta de cuatro columnas, con diámetro de 3.5 cm y altura de 18.5 cm, equipado con un tanque de alimentación con el agua sintética preparada, una bomba sumergible que alimenta el agua hacia el tanque de distribución donde se mantiene a un nivel constante por rebose, después por gravedad se distribuye el agua a las columnas, como se muestra en la Fig. 1 , con control de flujo manual. Se realizó un montaje experimental variando la altura del lecho durante 7.5 h, se utilizó para este estudio un lecho empacado de 4 y 7.5 cm, por el cual pasó la solución de agua sintética preparada a una velocidad de flujo de 6 mL/min [15]. Los ensayos se realizaron con réplica para la confiabilidad de los resultados obtenidos.

3. Resultados y Discusión

3.1 Caracterización: análisis proximal elemental y FTIR

En los resultados presentados en la Tabla 1, se observa un alto contenido de carbono en el adsorbente (50.35 %); y que los contenidos de celulosa, lignina, pectina y hemicelulosa tienen un porcentaje de 19.82 %, 12.66 %, 9.54 % y 9.45 % respectivamente.

Los espectros FTIR de la cáscara de cacao se localizaron entre 4000 y 400 cm^-1. La Fig. 2 representa el cambio del espectro de la biomasa tratada solo con agua destilada, con respecto al obtenido después de colocar la biomasa en contacto con la solución de mercurio.

En la Tabla 2 se muestra el estado de energía vibracional y el grupo funcional correspondiente a la longitud de onda del espectro, permitiendo relacionar estas tipologías con los diferentes componentes y establecer cualitativamente las especies que pueden propiciar la adsorción; y en la Tabla 3 los grupos funcionales que se mantuvieron luego de ser utilizada la biomasa como adsorbente.

3.2 Efecto de la altura del lecho empacado

Se valoró estableciendo el porcentaje de adsorción y el tiempo de ruptura de cada columna para la remoción de mercurio. De acuerdo con la Tabla 4, el porcentaje de remoción aumenta con la altura del lecho empacado de biomasa, siendo la altura de 7.5 cm la que favorece el proceso respecto a la de 4 cm.

Donde, t(5%) es el tiempo de servicio en minutos conseguido cuando la concentración de la solución en la salida es el 5 % de la concentración de entrada, C0 es la concentración en mg/L de metal en la entrada a la columna, Qv es el caudal (L/min) y m es la cantidad de biomasa de cascara de cacao. Los parámetros generados del análisis de los datos de las curvas de ruptura se muestran en la Tabla 5.

3.3 Modelado de la curvas de ruptura

La eficiencia de una columna de lecho empacado se establece mediante el concepto de curva de ruptura, la cual es una característica que permite determinar el tiempo de servicio y establecer la respuesta dinámica de la columna; el desempeño del ion retenido en lecho de cacao; se expresa en términos de la variación de la concentración inicial con respecto al tiempo, Ct/Co, la cual es una función del volumen del efluente a una altura dada, generando de esta manera la curva respectiva [15-19] . Con el fin de analizar el rendimiento de la columna se realizó el ajuste del modelo cinético de los datos experimentales a los siguientes modelos:

3.3.1 Modelo de Adams-Bohart (AB)

Este modelo es usado en la descripción de la parte inicial de la curva de ruptura, asumiendo que la velocidad de adsorción es proporcional a la capacidad de adsorción del adsorbente y a la concentración del soluto, y se enuncia en (2) [15, 17].

En esta expresión, kAB es la constante cinética de adsorción (L mg ^-1min^-1), Fes la velocidad lineal (cm min ^-1), Z es la altura del lecho empacado y N0 es la concentración de saturación del adsorbente (mg L ^-1); estos parámetros se establecen a partir del gráfico Ct/C0 vs t.

Es utilizado para graficar y analizar los datos experimentales derivados de ensayos en columnas, con el fin de determinar la máxima capacidad de adsorción del material y su constante cinética; y se describe mediante (3) [15, 18]:

Donde, kTh es la constante cinética (mL min ^-1 mg ^-1), q0 es la cantidad de soluto adsorbida por gramo de adsorbente (mg g ^-1), x es la masa de adsorbente en la columna (g), C0 es la concentración inicial del efluente (mg L ^-1), Ct es la variación de la concentración del efluente en el tiempo t(mg L-1) y v es el caudal de entrada (mL min ^-1). El coeficiente cinético kThy la capacidad de adsorción de la columna, q0, pueden ser determinados por regresión no-lineal del gráfico Ct/C0 vs t,a una velocidad de entrada dada.

3.3.3 Modelo de Yoon-Nelson (YN)

El modelo linealizado para un sistema monocomponente se expresa en (4); y asume que la tasa de disminución en la probabilidad de adsorción para cada molécula de adsorbato es proporcional a la probabilidad de adsorción de adsorbato y la probabilidad de ruptura de adsorbato en el adsorbente. No requiere de datos detallados sobre las características del adsorbato, el tipo de adsorbente, y las propiedades físicas del lecho de adsorción [15, 17]:

Donde kYN es la constante de velocidad (min^-1) y τ (min) el tiempo necesario para una penetración del adsorbato del 50 %.

En la Fig. 3 (a) se muestra el comportamiento de la curva de ruptura de la adsorción de Hg (II) sobre una columna de 4 cm de cacao y el modelado de los datos experimentales para un tiempo de 450 min. Por su parte, en la Fig. 3 (b), se observa el desempeño de la curva de ruptura de la columna de 7.5 cm para el mismo tiempo de operación y adsorbato.

4. Discusión

El análisis químico proximal de la Tabla 1 presenta los principales componentes propiciadores de la remoción de iones metálicos, puesto que en su estructura cuentan con gran presencia grupos OH (hidroxilos y fenólicos); por lo tanto es de esperar que la cáscara de cacao sea un adsorbente con gran capacidad de remoción, ya que su naturaleza lignocelulósica confirma la presencia de compuestos con largas cadenas poliméricas (celulosa, hemicelulosa, pectina y lignina), que son reconocidos por la gran cantidad de grupos hidroxilos y carboxilos que presentan en su disposición [4-6] , [15] , [16] . En el análisis del espectro FTIR de la Fig. 2, se identificó la presencia de: hidrocarburos alifáticos, aromáticos e insaturados con una energía vibracional variable, compuestos orgánicos, halógenos, metilo, sílica y nitrógeno y grupos hidroxilo, de acuerdo con la Tabla 2. La vibración en la región entre 3900 a 3200 cm^-1 confirma la presencia de grupos OH, los cuales benefician la remoción de metales pesados, debido a que los grupos de carga aniónica en la superficie del sorbente atraen electrostáticamente a los cationes. El espectro presenta una tendencia decreciente en este rango, porque el enlace de hidrógeno tiende a ampliar los picos y desplazarlos hacia números de onda más bajos. Por tanto, la forma ancha y de gran intensidad del pico sobre 3627.29- 3167.25 cm^-1 es típica [9]. Con base en lo reportado en la Tabla 3, se puede decir que el espectro de la biomasa luego de la adsorción de mercurio y el espectro de la biomasa natural presenta la misma tendencia, debido a que después se conservaron varias de las vibraciones características de los grupos funcionales principales presentes en la cadena.

Teniendo en cuanta la Tabla 4, la relación entre el tiempo y el porcentaje de remoción es inversamente proporcional, esto puede ser ocasionado porque el fluido con un contenido constante de contaminante, que es alimentado a la columna, es inicialmente removido en la zona inferior del lecho, por lo que la solución a la salida no contiene adsorbato en este punto. Sin embargo, a medida que aumenta el tiempo de residencia el ion se va ligando al adsorbente en zonas cada vez más altas, lo cual se traduce en el movimiento de la zona de transferencia de masa con el tiempo; por consiguiente, una vez la zona de transferencia de masa se ubica en la parte más alta de la columna, quiere decir que se ha saturado el lecho [10, 15, 19].

Estudiando el comportamiento dinámico de la columna de lecho fijo, se obtuvieron trup de 0 y 120 min para la altura del lecho de 4 cm y 7.5 cm un tiempo de ruptura de 0 y 120 min, respectivamente, como se muestra en la Tabla 5. Por lo tanto, el aumento de la altura del lecho favorece el tiempo de ruptura, ya que aumenta la cantidad de sitios de transferencia de masa. Además, se estableció que el cambio del trup es proporcional a la altura del lecho, lo cual significa que se adsorbe mayor cantidad de soluto por gramo de adsorbente, lo que es coherente con los porcentajes de remoción presentados. Con base en este desempeño, es posible afirmar que el proceso de adsorción de iones de Hg (II), sobre un lecho de cacao, es dependiente de la altura, ya que el incremento en la altura del lecho también significa un aumento el área superficial disponible y el tiempo de contacto entre los iones de Hg y el sorbente [13, 15].

Al modelar los datos experimentales de las curvas de ruptura se obtuvo que el modelo de Thomas es el que mejor los ajusta; estos resultados son coherentes con resultados obtenidos por otros investigadores en estudios similares; este ajuste es apropiado para flujo constante en columna sin dispersión axial, y su comportamiento coincide con el modelo de isoterma de Langmuir (adsorción en monocapa) y modelo cinético reversible de pseudo-segundo orden; además, se observa una relación proporcional entre la altura del lecho y la constante de velocidad de Thomas (kTH), lo cual describe un aumento en velocidad de transporte del ion de la solución acuosa al adsorbente y un aumento de kTH y qe, como se muestra en la Tabla 6 [15-18].

5. Conclusiones

Teniendo en cuenta el análisis FTIR se estableció que, el grupo funcional característico en la adsorción de iones metálicos presentes en solución acuosa es el grupo hidroxilo (OH-), el cual se encuentra en la superficie de la cáscara de cacao. Los ensayos de adsorción evidenciaron que la biomasa residual de cacao es eficiente como adsorbente de Hg (II), removiéndose hasta el 99.62 % del contaminante usándose una columna con un lecho de 7.5 cm. Se determinó que la altura del lecho es una variable que influye de manera proporcional con el tiempo de servicio, siendo el modelo de Thomas el que mejor ajustó los datos experimentales de adsorción.

6. Agradecimientos

Referencias

[1] M. F. Castillo, M. Ramírez, R. S. García, M. Bernal, B. Espinosa, J. A. Solís, C. Durán,“Reaprovechamiento integral de residuos agroindustriales: cáscara y pulpa de cacao para la producción de pectinas,” Rev. Latinoam. el Ambient. y las Ciencias, vol. 1, no. 2, pp. 45-66, 2010.         [ Links ]

[2] J. G. Ramírez Gil, “Pérdidas económicas asociadas a la pudrición de la mazorca del cacao causada por Phytophthora spp., y Moniliophthora roreri (Cif y Par) Evans et al., en la hacienda Theobroma, Colombia,” Rev. Protección Veg., vol. 31, no. 1, pp. 42-49, 2016.         [ Links ]

[3] N. J. Cárdenas Pardo, A. E. Darghan Contreras, M. D. Sosa Rico, and A. Rodríguez, “Análisis espacial de la incidencia de enfermedades en diferentes genotipos de cacao (Theobroma cacao L.) en El Yopal (Casanare), Colombia,” Acta Biológica Colomb., vol. 22, no. 2, pp. 209-220, 2017.         [ Links ]

[4] C. Tejada, A. Herrera, and E. Ruiz, “Kinetic and isotherms of biosorption of Hg ( II ) using citric acid treated residual materials,” Ing. y Competividad, vol. 18, no. 1, pp. 117-127, 2016.         [ Links ]

[5] C. Tejada, Á. Villabona, and J. Núñez, “Uso de biomasas para la adsorción de plomo , níquel , mercurio y cromo,” Univ. Cart. Colomb., vol. 9, pp. 41-51, 2015.         [ Links ]

[6] C. Tejada-Tovar, Á. Villabona-Ortiz, and L. Garcés-Jaraba, “Adsorción de metales pesados en aguas residuales usando materiales de origen biológico,” Tecno Lógicas, vol. 18, no. 34, pp. 109-123, 2015.         [ Links ]

[7] V. R. García, N. Borja, E. Guzmán, A. G. Yipmantin, and H. Maldonado, “Equilibrio de biosorción de plomo (II) y caracterización mediante FT-IR y SEM-EDAX en pectina reticulada proveniente de cáscaras de naranja,” Rev. Soc. Química Peru, vol. 79, no. 3, pp. 256-264, 2013.         [ Links ]

[8] C. Tejada Tovar, Á. Villabona Ortiz, and V. Ruiz Rangel, “Biomasa residual para remoción de mercurio y cadmio: una revisión,” Ingenium, vol. 6, no. 14, pp. 11-21, 2012.         [ Links ]

[9] C. T. Tovar, Á. V. Ortiz, and M. J. Villadiego, “Remoción de cromo hexavalente sobre residuos de cacao pretratados químicamente,” Rev. U.D.C.A Act. Div. Cient., vol. 20, no. 1, pp. 139-147, 2017.         [ Links ]

[10] S. Mädler, F. Sun, C. Tat, N. Sudakova, P. Drouin, R. J. Tooley, E. J. Reiner, T. A. Switzer, R. Dyer, H. M. Skip Kingston, M. Pamuku and V. I. Furdui “Trace-Level Analysis of Hexavalent Chromium in Lake Sediment Samples Using Ion Chromatography Tandem Mass Spectrometry,” J. Environ. Prot. (Irvine,. Calif)., vol. 7, no. February, pp. 422-434, 2016.         [ Links ]

[11] X. Li, H. Zhou, W. Wu, S. Wei, Y. Xu, and Y. Kuang, “Studies of heavy metal ion adsorption on Chitosan/Sulfydryl-functionalized graphene oxide composites,” J. Colloid Interface Sci., vol. 448, pp. 389-397, 2015.         [ Links ]

[12] R. Kumar, D. K. Arya, N. Singh, and H. K. Vats, “Removal of Cr (VI) Using Low Cost Activated Carbon Developed By Agricultural Waste,” IOSR J. Appl. Chem., vol. 10, no. 1, pp. 76-79, 2017.         [ Links ]

[13] S. Rangabhashiyam and N. Selvaraju, “Adsorptive remediation of hexavalent chromium from synthetic wastewater by a natural and ZnCl2 activated Sterculia guttata shell,” J. Mol. Liq., vol. 207, pp. 39-49, 2015.         [ Links ]

[14] M. Hadavifar, N. Bahramifar, H. Younesi, and Q. Li, “Adsorption of mercury ions from synthetic and real wastewater aqueous solution by functionalized multi-walled carbon nanotube with both amino and thiolated groups,” Chem. Eng. J., vol. 237, pp. 217-228, 2014.         [ Links ]

[15] J. Lara, C. Tejada, A. Villabona, and A. Arrieta, “Adsorción de plomo y cadmio en sistema continuo de lecho fijo sobre residuos de cacao,” Rev. Ion., vol. 29, no. 2, pp. 111-122, 2016.         [ Links ]

[16] C. Tejada Tovar, A. Villabona Ortiz, and E. Ruiz Paternina, “Cinética de adsorción de Cr (VI) usando biomasas residuales modificadas químicamente en sistemas por lotes y continuo,” Rev. Ion, vol. 28, no. 1, pp. 29-41, 2015.         [ Links ]

[17] S. M. El-Bahy and Z. M. El-Bahy, “Immobilization of 2-amino pyridine onto poly(acrylonitrile-co-N,N´-methylenebisacrylamide) nanoparticles for the removal of Hg(II), Cd(II) and Cr(III): Batch and column techniques,” J. Environ. Chem. Eng., vol. 5, no. 4, pp. 3560-3571, 2017.         [ Links ]

[18] J. L. Gong, L. Zhang, Y. Jiang, G. M. Zeng, Z. H. Cui, K. Liu, C. H. Deng, Q. Y. Niu, J. H. Deng and S. Y. Huan, “Continuous adsorption of Pb(II) and methylene blue by engineered graphite oxide coated sand in fixed-bed column,” Appl. Surf. Sci., vol. 330, pp. 148-157, 2015.         [ Links ]

[19] O. Allahdin, J. Mabingui, M. Wartel, and A. Boughriet, “Removal of Pb 2+ ions from aqueous solutions by fixed-BED column using a modified brick: (Micro)structural, electrokinetic and mechanistic aspects,” Appl. Clay Sci., vol. 148, no. July, pp. 56-67, 2017.         [ Links ]