Source: http://docplayer.es/133448-Projecte-o-tesina-d-especialitat.html
Timestamp: 2018-09-25 04:08:36+00:00

Document:
PROJECTE O TESINA D ESPECIALITAT - PDF
Download "PROJECTE O TESINA D ESPECIALITAT"
Pablo Álvarez Vidal
1 PROJECTE O TESINA D ESPECIALITAT Títol TRANSPORT DE NANOPARTÍCULES EN MEDI PORÓS Autor/a HUGO OLMO SÁEZ Tutor/a XAVIER SÁNCHEZ VILA Departament Intensificació HIDROLOGÍA SUBTERRÀNEA Data GENER 2014
3 AGRADECIMIENTOS En primer lugar un especial agradecimiento a Xavier Sánchez Vila por su dedicación y ayuda durante la realización del trabajo. También muy importante han sido durante toda la carrera el apoyo y ayuda de mis padres y hermano Xavi. También merece un personal reconocimiento Ángel Guasch, amigo y compañero de clase, el cual me ha ayudado en la realización del código del programa MATLAB. 3
4 RESUMEN En el presente trabajo se describe de forma general el transporte de nanopartículas en medio poroso. El objetivo principal del trabajo es la resolución de la ecuación que rige dicho transporte y conocer el comportamiento de las nanoparticulas a lo largo de su recorrido a través del medio poroso. Cabe destacar el papel importante que tienen las nanoparticulas en la actualidad, con importantes y variables aplicaciones. En concreto, en el campo de la ingenieria sanitaria tiene una aplicación muy importante que permite reducir la concentración de diversos contaminantes en el agua. En esta tesina se presentan diversos aspectos relacionados con las propiedades de las nanoparticulas en relación al medio poroso y en base a los conceptos desarrollados, se plantea una ecuación matemática que permita analizar el transporte y comportamiento de las mismas. Para la resolución de la ecuación diferencial en derivadas parciales propuesta, se utiliza el método de las diferencias finitas y se desarrolla una herramienta computacional para dicho modelo mediante el programa MATLAB. 4
5 ABSTRACT The present work aims at describing transport of nanoparticles in porous media. The main contribution of this work is to present the governing equation for the fate of nanoparticles in porous media and to provide a numerical solution to this equation. It is necessary to emphasize the importance that nanoparticles have at present, with important and variable applications. In particular, in the field of sanitary engineering they can be used as a way to reduce the concentration of certain pollutants in water. This thesis deals with different aspects related to the properties of nanoparticles as related to porous medium, and based on these properties we postulate an equation that governs the fate of nanoparticles in aquifers. The resulting partial differential equation is solved using finite differences, in a code written in MATLAB. 5
6 ÍNDICE AGRADECIMIENTOS... 3 RESUMEN... 4 ABSTRACT... 5 ÍNDICE DE FIGURAS... 8 ÍNDICE DE TABLAS ) INTRODUCCIÓN Y OBJETIVOS Introducción y objetivos Metodología ) NANOPARTÍCULAS Definiciones previas Origen Características Clasificación Nanopartículas de origen natural Nanopartículas generadas de forma involuntaria por la actividad humana Nanopartículas artificiales Aplicaciones Eliminación de contaminantes del agua mediante el uso de nanopartículas Carácter tóxico de las nanoparticulas ) MEDIO POROSO Definición Fases del medio poroso Parámetros importantes del medio poroso Porosidad Área específica Tortuosidad Ecuación del flujo en medio poroso Ley de Darcy Principio de conservación de la masa Ecuación del flujo
7 ÍNDICE DE FIGURAS Figura 1: Imagen fullenero Figura 2: Imagen nanotubo de carbono Figura 3: Imagen nanohilos Figura 4: Imagen grafeno Figura 5: Imagen nanoespuma de carbono Figura 6: Geometría del medio poroso Figura 7: Fases del medio poroso Figura 8: Representación de la tortuosidad Figura 9: Experimento de Darcy Figura 10: Representación de los elementos de la ecuación de Bernoulli Figura 11: Representación del gradiente hidráulico entre dos líneas equipotenciales Figura 12: Fenómeno de advección en transporte de solutos Figura 13: Fenómeno de difusión en partículas Figura 14: Fenómeno de dispersión mecánica en suelos porosos Figura 15: Malla de puntos en el método de las diferencias finitas Figura 16: Representación del método implícito Figura 17: Gráfica 3D de la evolución de la concentración de nanopartículas en el proceso Figura 18: Gráfica 3D de la evolución de la concentración de partículas retenidas en el proceso Figura 19: Corte de las gráficas 3D al inicio y al final del medio poroso en el proceso Figura 20: Grafica 3D de la evolución de la concentración de nanopartículas en el proceso Figura 21: Gráfica 3D de la evolución de la concentración de nanopartículas retenidas en el proceso Figura 22: Corte de las gráficas 3D al inicio y al final del medio porosoen el proceso ÍNDICE DE TABLAS Tabla 1: Aplicaciones nanopartículas Tabla 2: Valores de la permeabilidad en un suelo (m/día) Tabla 3: Valores del coeficiente de difusión molecular
8 CAPÍTULO 1: INTRODUCCIÓN Y OBJETIVOS 9
9 1) INTRODUCCIÓN Y OBJETIVOS 1.1 Introducción y objetivos Las nanopartículas (NP) se utilizan en un conjunto de aplicaciones en ciencias de materiales, medicina y diversas aplicaciones industriales. Además pueden interaccionar entre ellas y con el medio que la envuelve. En estos momentos, se empieza a pensar en utilizar las nanopartículas para analizar el movimiento de solutos en medio poroso, con implicaciones potenciales en descontaminación de recursos hídricos y en evaluación de los servicios ecosistémicos. La hidrología subterránea nos proporciona modelos teóricos que describen, por un lado el movimiento del agua en el acuífero, y por otro el movimiento de un determinado soluto en el agua subterránea. En nuestro caso nos permitirá conocer el movimiento de las partículas en aguas subterráneas. El objetivo de este trabajo es realizar una recopilación de los procesos que gobiernan el transporte de nanopartículas en medios porosos saturados (acuíferos), y formular estos procesos en forma de ecuaciones diferenciales en derivadas parciales de gobierno, así como proponer un método para la resolución numérica (diferencias finitas) de la ecuación de transporte resultante. 1.2 Metodología Las distintas etapas que ha seguido el trabajo se pueden resumir de la siguiente manera: 1) Estudio de las características de las nanopartículas y de sus potenciales aplicaciones, con énfasis en aquellas relacionadas con el medio poroso. 2) Explicación de las características del medio poroso y ecuaciones que rigen el transporte de solutos (sustancias disueltas) en medios poroso. 3) Formulación de la ecuación de transporte de nanopartículas en medio poroso. 10
10 4) Resolución de la ecuación de transporte de nanoparticulas mediante el método de las diferencias finitas. 5) Adquirir unos conocimientos básicos de programación en MATLAB y escritura de un código numérico que permita resolver la ecuación de transporte de nanopartículas mediante el método de las diferencias finitas. 6) Valoración de los resultados obtenidos. 11
11 CAPÍTULO 2: NANOPARTÍCULAS 12
12 2) NANOPARTÍCULAS 2.1 Definiciones previas Nanopartícula: Partícula con una o más dimensiones del orden de 100 nm. Nanómetro: Unidad de longitud equivalente a una millonésima de milímetro (10-9 m). Símbolo: nm. Nanotecnología: Diseño, caracterización, producción y aplicación de estructuras, dispositivos y sistemas mediante el control del tamaño y la forma a una escala nanométrica. Nanociencia: El estudio de los fenómenos y la manipulación de materiales a escala atómica, molecular y macromolecular. Nanoescala: Tener una o más dimensiones representativas del orden de 100 nm. Nanomateriales: Material con una o más dimensiones externas o una estructura interna a nanoescala. Nanoestructura: Estructura con una o más dimensiones en la nanoescala. Nanocompuestos: Compuestos en los que al menos una de las fases tiene una dimensión de la nanoescala. 2.2 Origen Aunque esta ciencia puede sonar relativamente nueva, las nanopartículas existen en el planeta desde hace siglos. Ejemplo de ello son las partículas de humo y los virus
13 Antiguamente algunas civilizaciones utilizaron nanopartículas aprovechando sus propiedades ópticas y medicinales, desconociendo obviamente el concepto moderno de las mismas. Un ejemplo es la civilización egipcia, en la que se usaron nanopartículas de oro como coloides medicinales. Estamos rodeados de miles de millones de nanopartículas, incluyendo el polvo de arena, sustancias químicas generadas por el plancton oceánico, los productos de combustión (principalmente los gases de combustión de los motores diésel), procesos naturales como volcanes e incendios. Todos ellos han generado nanopartículas y otras nanoestructuras que permanecen en el ambiente a lo largo del tiempo. El ambiente habitual puede contener nanopartículas por centímetro cúbico. En un bosque, esta cifra puede elevarse a 50000, y en una calle de la ciudad a nanopartículas. La proporción de nanopartículas artificiales procedentes de la producción industrial es todavía muy pequeña, pero en el futuro se espera que las nanopartículas se fabriquen en mayores cantidades, para obtener nuevos productos, de tal manera que a medida que sus usos se extiendan, las posibilidades de exposición tanto del ser humano como del medioambiente se irán incrementando. 2.3 Características Una de las características de las nanopartículas es que la relación entre el número de átomos superficiales y el tamaño de la partícula es de carácter exponencial. Por ello, las propiedades relacionadas con la superficie, como las eléctricas, mecánicas, magnéticas, ópticas o químicas de los nanomateriales son diferentes a las de los mismos materiales a escala no nanométrica. Un ejemplo de ello es el color. Cuanto más pequeña sea una partícula de cualquier material, mayor es su superficie en relación con su masa. Con frecuencia, su reactividad y por tanto su toxicidad se ven 14
14 también aumentadas (ya que el número de átomos en contacto con la superficie es mayor). Si, además, la superficie de una partícula se ha modificado para lograr un determinado comportamiento, esto puede dar lugar a interacciones inesperadas con moléculas biológicas. En los sistemas nanométricos, las superficies pasan a jugar un papel más importante que el que tenían en los sistemas macroscópicos. Esto es así porque, al reducir la escala, el número de átomos que tenemos en la superficie con respecto al número de átomos que tenemos en volumen aumenta drásticamente. Los principales parámetros característicos de las nanopartículas son su forma, tamaño, características de la superficie y la estructura interna, aunque en la escala nanométrica, las propiedades de las partículas pueden cambiar de manera impredecible. Las nanopartículas tienen diferentes maneras de interactuar entre sí. Éstas pueden permanecer libres o agruparse en función de las fuerzas de interacción atractiva o repulsiva entre ellas. Estas interacciones siguen siendo difíciles de caracterizar. 2.4 Clasificación Las nanopartículas pueden clasificarse en tres grandes categorías según su origen Nanopartículas de origen natural Algunas son de origen biológico, como por ejemplo muchos virus y bacterias y otras son de origen mineral o medioambiental como las que contiene el polvo de arena o las nieblas y humos derivados de la actividad volcánica o de los fuegos forestales Nanopartículas generadas de forma involuntaria por la actividad humana Las nanopartículas como consecuencia de la actividad humana pueden ser generadas de forma involuntaria o deliberada. 15
15 Las nanopartículas producidas de forma involuntaria son las que se producen en ciertos procesos industriales bien conocidos, tales como la producción de materiales a gran escala por procedimientos a altas temperaturas (como el humo de sílice, partículas ultrafinas de óxido de titanio y metales ultrafinos), procesos de combustión (diésel, carbón), obtención de pigmentos, procesos mecánicos (lijado, mecanizado, pulido, etc.) o en procesos domésticos (barbacoas, humos de aceite, etc.) Nanopartículas artificiales Las nanopartículas generadas deliberadamente se producen mediante las llamadas nanotecnologías. Los métodos para la obtención de nanopartículas son, a grandes rasgos, de dos tipos: los llamados top-down, en los que se llega a nanomateriales sometiendo los materiales convencionales a diversos procesos y los bottom-up en los que se construyen nanopartículas a partir de átomos o moléculas. Se puede hacer otra clasificación en función del número de dimensiones que en la estructura considerada tengan carácter nanométrico, los nanomateriales los podemos encontrar en: Tres dimensiones a escala nanométrica: nanocristales, fullerenos y puntos cuánticos. Dos dimensiones a escala nanométrica: nanotubos y los nanohilos. Una dimensión a escala nanométrica: estructuras que se utilizan en los recubrimientos de superficies o películas finas en los que sólo su grosor es de orden nanométrico. A continuación se explican los tipos de nanopartículas artificiales Fullerenos Son estructuras cerradas formadas por átomos de carbono dispuestos en forma de pentágonos y hexágonos a modo de nanobalones de fútbol. Entre sus propiedades 16
16 físicas destaca la de que son capaces de resistir presiones extremas y recuperar su forma original cuando cesa la presión. Figura 1: Imagen fullenero. Extraída de imágenes de Google Puntos cuánticos Los puntos cuánticos son un conjunto de materiales entre los 2 a 10 nm, que pueden estar compuestos de metales, óxidos metálicos o materiales semiconductores y ofrecen propiedades electrónicas, ópticas, magnéticas y catalíticas Nanotubos de carbono Pueden describirse como un tubo cuya pared es una malla de agujeros hexagonales. También es posible obtener nanotubos de varias capas, a modo de varios tubos concéntricos. Son muy destacables sus propiedades eléctricas y mecánicas. Son grandes superconductores capaces de resistir el paso de corrientes elevadísimas, de densidades de corriente de hasta mil millones de amperios por metro cuadrado, y, por otra parte su resistencia mecánica es sesenta veces superior a la de los mejores aceros; a su vez son ligeros y flexibles ya que su peso específico es más de seis veces inferior a la de los mencionados aceros. 17
17 Existen principalmente dos tipos de nanotubos de carbono: con una única pared y con múltiples paredes. Los nanotubos de carbono de pared simple son cilindros de grafito de una sola capa que tiene diámetros del orden de unos pocos nanómetros, mientras que los múltiples paredes forman parte de 2 a 30 cilindros concéntricos que tienen diámetros exteriores entre 30 y 50 nm. Figura 2: Imagen nanotubo de carbono. Extraída de imágenes de Google Nanohilos Son también estructuras alargadas que destacan por sus propiedades como conductores o semiconductores. Se han obtenido nanohilos de diversos materiales como silicio, cobalto, oro y cobre. Sus aplicaciones más importantes pertenecen al campo de la nanoelectrónica. 18
18 Figura 3: Imagen nanohilos. Extraída de imágenes de Google Grafeno En el grafeno los átomos de carbono se unen en láminas planas de un átomo de espesor, formando un panel de abejas hexagonal (con un átomo en cada vértice). El grafeno es un semiconductor que puede operar a escala nanométrica y a temperatura ambiente, con propiedades que ningún otro semiconductor ofrece. Este material ha llegado para revolucionar la electrónica permitiendo fabricar dispositivos electrónicos mucho más pequeños que cualquiera de los que existen actualmente. Figura 4: Imagen grafeno. Extraída de imágenes de Google. 19
19 Sector Aplicaciones Automoción, aeronáutica Materiales reforzados y más ligeros; pinturas más brillantes, con efectos de color, resistentes a arañazos, anticorrosivas y antiincrustantes; aditivos para diésel que permiten una mejor combustión; neumáticos más duraderos y reciclables. Electrónica y comunicación Memorias de alta densidad y procesadores miniaturizados, células solares, tecnología inalámbrica, pantallas planas. Química y materiales Pinturas anticorrosión e ignífugas, catalizadores, textiles con recubrimientos antibacterianos y ultrarresistentes. Materiales más ligeros y resistentes. Farmacia, biomedicina y biotecnología Cosmética Medicamentos a medida liberados en órganos específicos, kits de autodiagnóstico, biosensores, prótesis, implantes. Cremas solares transparentes, nanoesferas con antioxidantes. Salud Dispositivos de diagnóstico y detección miniaturizados, destrucción de tumores por calor, terapia génica, microcirugía y medicina reparadora: nanoimplantes y prótesis, membranas para diálisis. Energía Generación de energía fotovoltaica, nuevos tipos de baterías, ventanas inteligentes, materiales aislante más eficaces. Medio ambiente y ecología Disminución de las emisiones de CO 2, producción de agua ultrapura a partir de agua de mar, pesticidas y fertilizantes más eficaces y menos dañinos, sensores para la limpieza del medio ambiente. 21
20 son esperanzadores, con eliminaciones próximas al 100% a concentraciones inferiores a 50 mg/l de fosfato y concentración de NP-CeO 2 de 250 mg/l. Por otro lado, también se ha procedido a anclar cadenas orgánicas sobre las nanopartículas con el objetivo de que sean estas las que retengan el contaminante. En este sentido, actualmente se están efectuando experimentos con nanopartículas de oro funcionalizadas con ciclodextrinas para eliminar pesticidas de varias formulaciones. En todos los casos se ha prestado especial atención a determinar el efecto que las nanopartículas pueden causar en el medio donde son liberadas. En este sentido se han realizado diferentes ensayos de toxicidad estandarizados con las nanopartículas. Los resultados obtenidos muestran, por lo tanto, que las nanopartículas presentan un potencial real para ser utilizadas en el tratamiento de aguas residuales. En algunos casos nos encontramos ante nanopartículas con cierta toxicidad ambiental pero con un elevado potencial descontaminante, por lo que hay que encontrar otras formas de utilizarlas, como por ejemplo la inmovilización de las NP sobre soportes porosos a través de los cuales se haría circular el agua contaminada. Esto facilitaría su aplicación. 2.7 Carácter tóxico de las nanoparticulas Un aspecto importante de las importantes de las nanopartículas y que muchas veces se debe tener en cuenta es su carácter tóxico. Una de las preocupaciones principales derivadas del uso y la manipulación de los nanomateriales es que su toxicidad parece deberse a ciertas propiedades fisicoquímicas que actualmente no están incluidas en ensayos estándar de toxicidad. También causa gran preocupación que el uso comercial de estas partículas puede producir una amplia dispersión de las mismas mediante su liberación a través de efluentes industriales y domésticos, y productos de higiene personal u otros productos ampliamente utilizados.las 23
21 implicaciones que dicha dispersión puede tener sobre la salud ambiental y humana son desconocidas en la actualidad. La información toxicológica del material en escala normal puede servir como base para las estimaciones preliminares sobre la toxicidad de una nanopartícula, junto con los estudios toxicológicos in vivo e in vitro y los estudios epidemiológicos que hayan podido realizarse. Los estudios experimentales en animales han mostrado que la respuesta biológica a ciertas nanopartículas puede ser mayor que la encontrada para la misma masa de partículas más grandes de composición química similar, debido probablemente al aumento del área superficial. Además del número de partículas y del área superficial, otras características podrían influir en la respuesta biológica, incluyendo la solubilidad, forma, carga, superficie química, propiedades catalíticas, contaminantes adsorbidos, así como el grado de aglomeración. Los principales factores que pueden determinar los efectos toxicológicos de los nanomateriales en el organismo son: 1) Factores que dependen de la exposición (vía de penetración, duración de la exposición, concentración). 2) Factores que dependen del organismo expuesto (susceptibilidad individual, realización de una actividad física en el lugar de trabajo, lugar de depósito, ruta que siguen las nanopartículas una vez que se han introducido en el organismo). 3) Factores relacionados con los nanomateriales (toxicidad intrínseca de la sustancia). 24
22 CAPÍTULO 3: MEDIO POROSO 25
23 Donde: = Densidad aparente. Definida como la relación entre la masa de suelo (excluyendo el agua y aire) dividida por el volumen total (M/L 3 ). = Densidad real o de partículas. Mide la relación entre la masa de las partículas que conforman el suelo y el volumen que ocupan, es decir, el volumen total, excluyendo el volumen ocupado por los poros (M/L 3 ). La porosidad depende de la composición, de la textura y de la estructura del suelo. Por lo general, varía de 0,4 a 0,6 y puede ser mayor en suelos con un contenido de materia orgánica (su forma irregular produce una escasa compactación). Un aumento en el contenido de arcilla favorece la formación de agregados del suelo e incrementa la porosidad Área específica El área específica de un sólido,, se define como el área intersticial de la superficie de poros por unidad de masa y es uno de los principales parámetros para evaluar la capacidad de absorción de los suelos. Varía desde unos pocos cm 2 /gr para rocas, llega a valores de entre 600 y 1000 m 2 /gr para aerogeles y alcanza valores de hasta m 2 /gr para carbones activados Tortuosidad El flujo en los poros del suelo no sigue un camino recto. El agua se mueve alrededor de partículas individuales y de poros de diferentes tamaños dentro del medio poroso, resultando un patrón mucho más largo que una línea recta. Se define y calcula la tortuosidad como la relación entre la longitud real que debe recorrer una partícula de fluido para unir dos puntos en el seno del medio poroso y la distancia en línea recta entre dichos puntos. (3.3) 28
24 entre el volumen de agua que atraviesa la columna por unidad de tiempo,, y las alturas piezométricas medidas por un manómetro en los extremos de la columna, y. La ley resultante se puede escribir como: Figura 9: Experimento de Darcy. Extraída de imágenes de Google. Donde: (3.4) = Caudal (L 3 /T). = Longitud en metros de la muestra (L). = Conductividad hidráulica o coeficiente de permeabilidad (L/T). = Área de la sección transversal de la muestra (L 2 ). = Altura, sobre el plano de referencia que alcanza el agua en un tubo colocado a la entrada de la capa filtrante. = Altura que alcanza el agua en un tubo colocado a la salida de la capa filtrante. 30
25 punto de la vertical de A varían y pero la suma ( ) puede permanecer constante (hipótesis de Dupuit). En este caso (aunque no necesariamente el más habitual en hidrogeología) la línea vertical se corresponde con una línea equipotencial. Matemáticamente: Donde: (3.7) = Altura piezométrica (L). = Altura de presión (L). = Altura de elevación (L) Figura 10: Representación de los elementos de la ecuación de Bernoulli. Extraída de imágenes de Google. Si se consideran dos puntos se obtiene el siguiente esquema. 33
26 Figura 11: Representación del gradiente hidráulico entre dos líneas equipotenciales. Extraída de imágenes de Google. En la figura 11 se puede observar que el flujo del agua se dirige de la línea equipotencial 1 a la línea equipotencial 2 y perpendicularmente a las mismas. El gradiente hidráulico está definido por: (3.8) Donde: =Diferencia de alturas piezométricas (L). Se puede observar que el gradiente hidráulico se puede definir también como la pendiente de la superficie piezométrica entre los puntos 1 y 2. 34
27 3.4.2 Principio de conservación de la masa La conservación de masa sobre cualquier volumen del acuífero implica: (3.9) Donde: = Densidad del agua (M/L 3 ). = Términos fuente/sumidero. Tiempo (T). En esta ecuación, el primer miembro representa la variación de masa de agua almacenada por unidad de volumen de acuífero. El primer sumando del segundo miembro representa las entradas menos las salidas de agua sobre cada volumen desde/hacia los volúmenes de acuífero adyacentes (expresadas como masa de agua por unidad de volumen de acuífero y por unidad de tiempo). Por último, el segundo sumando representa las entradas menos las salidas desde/hacia el exterior. Por otro lado, es habitual aproximar el primer sumando linealmente como: (3.10) Donde: = Coeficiente de almacenamiento específico (volumen de agua extraído por unidad de volumen de acuífero y por unidad de descenso de nivel). 35
28 3.5.1 Advección La advección es el arrastre de soluto por el agua. Si sólo existiera este proceso, el soluto viajaría a la misma velocidad del agua y la extensión ocupada por el soluto sería constante. La advección simplemente transporta los solutos. Figura 12: Fenómeno de advección en transporte de solutos. Extraída del documento de Javier Sánchez, Transporte de contaminantes, Universidad de Salamanca. En un medio poroso, el flujo de masa a través de una unidad perpendicular al flujo es igual a: (3.13) Donde: = Flujo de masa advectiva por unidad de sección y tiempo (M/L 2 T). = Concentración de partículas suspendidas (M/L 3 ) Difusión Movimiento molecular causado por el movimiento aparentemente aleatorio de las partículas (movimiento Browniano). Produce un flujo neto de partículas desde las zonas donde hay una mayor concentración hacia las zonas con menor concentración. Ello provoca un efecto de expansión y disolución. 37
29 Figura 13: Fenómeno de difusión en partículas. Extraída del documento de Javier Sánchez, Transporte de contaminantes, Universidad de Salamanca. En la difusión, las sustancias disueltas se mueven por un gradiente de concentraciones. En un medio poroso, el flujo de masa por difusión está regido por la primera ley de Fick y presenta el siguiente aspecto. Donde: (3.14) = Flujo de masa por difusión por unidad de área y masa (M/L 2 T). = Coeficiente de difusión molecular (L 2 /T). = Tensor de difusión molecular = A continuación se muestra algunos valores que tiene el coeficiente de difusión molecular. 38
30 Tabla 3: Valores del coeficiente de difusión molecular. Extraída del documento de Jesús Carrera, Transporte de contaminantes, UPC Dispersión mecánica Variación de velocidades en el espacio debido a la geometría de los poros y sus interconexiones. Esta dispersión se produce en dirección del flujo (longitudinalmente) y lateramente (transversalmente). Figura 14: Fenómeno de dispersión mecánica en suelos porosos. Extraída del documento de Javier Sánchez, Transporte de contaminantes, Universidad de Salamanca. 39
31 La dispersión longitudinal es debida a: Las moléculas que se encuentran en caminos más tortuosos se retrasan. Las moléculas que encuentran caminos más anchos avanzan más rápido. La dispersión transversal es debida a la constante bifurcación de los caminos que entran en el fluido. La capacidad de un medio poroso para dispersar mecánicamente de un fluido que circula por él se refleja en un coeficiente denominado dispersividad dinámica (L), en el influirá la porosidad, tortuosidad, forma de los granos, etc. Se distingue la dispersividad dinámica longitudinal (en el sentido del flujo) y dispersividad dinámica transversal. La dispersión mecánica presenta la siguiente expresión. (3.15) Donde: = Flujo de masa por dispersión mecánica por unidad de tiempo y área (M/L 2 T). = Tensor de dispersión mecánica. Además en la dispersión mecánica encontramos la dispersión hidrodinámica, que es la acción conjunta de la difusión y la dispersión mecánica; ambos fenómenos no pueden considerarse aisladamente. Para tomarlos en consideración de un modo conjunto, se establece el coeficiente de dispersión hidrodinámica. Desglosando este concepto en dirección del flujo (longitudinal) y perpendicular al flujo (transversal), queda expresado en términos del módulo de la velocidad según: 40
32 CAPÍTULO 4: TRANSPORTE DE NANOPARTÍCULAS EN MEDIO POROSO 43
33 CAPÍTULO 5: RESOLUCIÓN DE LA ECUACIÓN DE TRANSPORTE DE NANOPARTICULAS MEDIANTE EL MÉTODO DE LAS DIFERENCIAS FINITAS 47
34 5) RESOLUCIÓN DE LA ECUACIÓN DE TRANSPORTE DE NANOPARTÍCULAS MEDIANTE EL MÉTODO DE LAS DIFERENCIAS FINITAS Una vez se ha obtenido la ecuación unidimensional de transporte de nanopartículas en medio poroso se procederá a su resolución con condiciones de contorno adecuadas. Una vez planteado el método se resolverá numéricamente con un código ad-hoc escrito en MATLAB. A continuación se presentan los posibles métodos numéricos de resolución de ecuaciones diferenciales en derivadas parciales acopladas. 5.1 Tipos de métodos numéricos Los métodos para la simulación de la ecuación de transporte de nanopartículas pueden ser clasificados de acuerdo con la forma de la ecuación que pretenden resolver, como Eulerianos, Lagrangianos, y mixtos. Seguidamente se hace una pequeña revisión de los tres tipos de métodos Métodos Eulerianos Los métodos Eulerianos se basan en la discretización de la ecuación de transporte en malla fija. Son los más antiguos y comunes para la solución de esta ecuación. Entre éstos, una de las formulaciones más ampliamente usadas es la basada en diferencias finitas. La otra gran formulación es la basada en elementos finitos. A pesar de que ambas formulaciones dan lugar a sistemas de ecuaciones similares, sus filosofías son muy diferentes. 48
35 Agrupando términos: ( ) ( ) ( ) (5.16) Donde: ( ) (5.17) ( ) (5.18) ( ) (5.19) (5.20) Finalmente la ecuación presenta el siguiente aspecto: (5.21) Por otra parte aplicando el método implícito a la segunda ecuación se obtiene: (5.22) 54
36 Donde: (5.23) (5.24) Por lo tanto se obtiene: (5.25) Finalmente se expresa el conjunto de ecuaciones resultante en forma matricial. (5.26) ( ) ( ) ( ) 55
37 CAPÍTULO 6: MATLAB PARA LA RESOLUCIÓN DE LA ECUACIÓN DE TRANSPORTE 56
38 Neumann. (6.4) Operando sobre la condición de Neumann se obtiene: ( ) ( ) ( ) (6.5) Donde: (6.6) (6.7) Por lo tanto: (6.8) La primera condición de contorno de Dirichlet muestra la concentración de entrada en un tiempo t al inicio de la columna. 58
39 Es importante que para no obtener oscilaciones y dispersión en los resultados, los valores de los elementos de la ecuación de transporte deben cumplir: Número de Grid-Courant. (6.9) Número de Grid-Peclet. (6.10) Finalmente el código desarrollado presenta el siguiente aspecto: function tesinalog() %Variables condición contorno close all D= ; L=0.05; v=0.003; n=50; Ax=L/n; At=0.9*Ax/v; porosidad=0.40; ro=1.3; ka=0.43/3600; kb=6.43/3600; TiempoFinal=40*At; q=20/3600; %Condiciones contorno iniciales condicioninicial=zeros(1,n*2); Ci=1*10^-5; C0=3/1000; valorcondicioninicialc=ci; condicioninicial(1,1:n)=valorcondicioninicialc; 59
40 %Matrices sistema AA=(1+(2*D*At)/(Ax.^2)); AB=(-D*At/(Ax.^2)+v*At/(2*Ax)); AC=(-D*At/(Ax.^2)-v*At/(2*Ax)); MatrizA=zeros(n,n); BD=ro/porosidad; MatrizB=zeros(n,n); CF=(porosidad*ka/ro)*At; MatrizC=zeros(n,n); DE=(1+At*kb); MatrizD=zeros(n,n); %Condición de contorno de Neumann M1=(-0.5*q/Ax)-(D/((Ax).^2)); M2=(0.5*q/Ax)+(D/((Ax).^2)); %Interior de matrices for i=1:1:n for j=1:1:n if i==j MatrizA(i,j)=AA;%AA MatrizB(i,j)=BD; MatrizC(i,j)=CF; MatrizD(i,j)=DE; elseif i==j-1 && j-1>0 MatrizA(i,j)=AB;%AB elseif i==j+1 MatrizA(i,j)=AC;%AC end end end %Interior matriz sistema MatrizSistema=zeros(2*n,2*n); MatrizSistema(1:n,1:n)=MatrizA(1:n,1:n); MatrizSistema(n+1:2*n,1:n)=MatrizC(1:n,1:n); MatrizSistema(1:n,n+1:2*n)=MatrizB(1:n,1:n); MatrizSistema(n+1:2*n,n+1:2*n)=MatrizD(1:n,1:n); 60
41 %Resolución sistema k=1; MatrizSolucion(1,2*n)=condicionInicial(1,2*n); for t=0:at:tiempofinal %Iteramos para cada uno de los tiempos k=k+1; if k==2 %Iteración inicial, correspondiente al vector de condiciones iniciales vectorcondiciones(1,1:n)=condicioninicial(1,1:n); vectorcondiciones(1,1)=c0; vectorcondiciones(1,n+1:2*n)=condicioninicial(1,1:n)*ka*porosidad/(kb*ro) ; SolucionC(1,1:n)=condicionInicial(1,1:n); SolucionS(1,1:n)=vectorCondiciones(1,n+1:2*n); else %Las demás iteraciones vectorcondiciones(1,1:n)=matrizsolucion(k- 1,1:n)+MatrizSolucion(k-1,n+1:2*n)*ro/porosidad; %Procesos posibles examinados %%Entra en un primer momento y a la mitad del tiempo deja de entrar concentración if t<1/2*tiempofinal vectorcondiciones(1,1)=c0; end vectorcondiciones(1,n+1:2*n)=matrizsolucion(k-1,n+1:2*n); %%Entra al principio, para, y después vuelve a entrar concentración % if t<=tiempofinal/2 && t>(1/4)*tiempofinal % vectorcondiciones(1,1)=c0; % end end MatrizSistema(n,n-1)=M1*At; MatrizSistema(n,n)=1+M2*At; vectorcondiciones=vectorcondiciones'; vectorresultado=linsolve(matrizsistema,vectorcondiciones); vectorresultado=vectorresultado'; MatrizSolucion(k,1:n*2)=vectorResultado; %Guardado de los resultados separados de C y S, pues en la matriz de %solución estos se modifican según Ck+(ro/porosidad)*Sk vectorcondiciones=0; SolucionC(k,1:n)=MatrizSolucion(k,1:n); SolucionS(k,1:n)=MatrizSolucion(k,n+1:2*n); 61
42 end %Generación de gráficas figure; yl=0:ax:l-ax; xl=0:at:tiempofinal+at; surf(yl(:),xl(:),solucionc); title('concentración C'); xlabel('longitud'); ylabel('tiempo'); zlabel('concentración C'); figure; yl=0:ax:l-ax; xl=0:at:tiempofinal+at; surf(yl(:),xl(:),solucions); title('concentracion S'); ylabel('tiempo'); xlabel('longitud'); zlabel('concentración S'); tiempo=0:at:tiempofinal+at; logt=log(tiempo(:)); logci=log(solucionc(:,1)); logcf=log(solucionc(:,n)); logsi=log(solucions(:,1)); logsf=log(solucions(:,n)); figure; subplot(1,2,1); plot(tiempo(:),solucionc(:,1)); title('concentración C') xlabel('tiempo'); ylabel('concentración C'); hold on plot(tiempo(:),solucionc(:,n),'r'); hold off; subplot(1,2,2); plot(tiempo(:),solucions(:,1)); title('concentración S') xlabel('tiempo'); ylabel('concentración S'); hold on plot(tiempo(:),solucions(:,n),'r'); hold off; 62
43 figure subplot(1,2,1); plot(logt(:),logci(:)); hold on; plot(logt(:),logcf(:),'r'); hold off; title('concentración C log') xlabel('log(tiempo)'); ylabel('log(concentración C)'); subplot(1,2,2); plot(logt(:),logsi(:)); title('concentración S log'); xlabel('log(tiempo)'); ylabel('log(concentración S)'); hold on; plot(logt(:),logsf(:),'r'); hold off; title('concentración S log') xlabel('log(tiempo)'); ylabel('log(concentración S)'); end 63
44 CAPÍTULO 7: RESULTADOS 64
45 En la figura 17 se puede observar como al principio hay entrada de flujo de nanopartículas y a la mitad del tiempo deja de entrar flujo de concentración de nanopartículas y en consecuencia baja la concentración de nanopartículas, obteniéndose un comportamiento esperado. También se observa un comportamiento diferente al inicio respecto al final del periodo, donde se perciben los cambios de entrada de concentración de manera más suave. Figura 18: Gráfica 3D de la evolución de la concentración de partículas retenidas en el proceso 1. Extraída del MATLAB. En esta figura 18 se observa la evolución de las partículas retenidas. Como es lógico en el instante inicial hay más partículas retenidas que en el instante final y el número va disminuyendo progresivamente. 66
46 Figura 19: Corte de las gráficas 3D al inicio y al final del medio poroso en el proceso 1. Extraída del MATLAB. En esta graficas se ha realizado un corte de las gráficas 3D de C y S al inicio i al final del medio poroso. En color azul se muestra el corte al inicio y en color azul al final. El la gráfica de C se puede observar una gráfica con dos curvas diferentes ya que al final padece los cambios de concentración menos brusca y de una forma más lenta y progresiva. En la gráfica de las partículas retenidas S se pude observar como tiene un comportamiento parecido al inicio y al final. 67
47 Figura 21: Gráfica 3D de la evolución de la concentración de nanopartículas retenidas en el proceso 2. Extraída del MATLAB. De la misma manera que en el proceso 1 la concentración de partículas retenidas disminuye progresivamente, ya que cada vez van saliendo más nanoparticulas y hay menos retenidas. 69
48 CAPÍTULO 8: CONCLUSIONES 71
49 8) CONCLUSIONES Es importante recalcar los siguientes aspectos, a manera de conclusiones generales respecto al trabajo realizado y el programa desarrollado para el cálculo de la ecuación de transporte de nanopartículas en medio poroso. Las nanoparticulas tienen muchas aplicaciones en la actualidad y tienen mucha posibilidad de progreso en el campo de la ingeniería sanitaria como medio para reducir la cantidad de contaminantes en el agua y por lo tanto mejorar la calidad del agua, Se ha empezado actualmente a conocer el carácter sanitario de las nanoparticulas y se debe investigar más para obtener un mejor rendimiento en su uso. Se debe tener en cuenta el carácter tóxico de las nanopartículas, ya que en muchas ocasiones puede llegar a ser perjudicial para la salud humana y el medio ambiente. El modelo de las diferencias finitas es un modelo adecuado y sencillo para el cálculo de ecuaciones diferenciales, concretamente la ecuación de transporte de nanopartículas en medio poroso. Existen más modelos para el cálculo de ecuaciones diferenciales, aunque se ha adoptado desarrollar el método de las diferencias finitas por su fácil comprensión y manejo. Se ha desarrollado un programa que permite resolver la ecuación transporte de nanopartículas en medio poroso, mediante técnicas de programación. Se comprueba que los resultados y graficas obtenidas son satisfactorios y cumplen con lo esperado. 72
50 CAPÍTULO 9: BIBLIOGRAFÍA 73
La seguridad y salud en la exposición a NANOPARTÍCULAS Instituto Riojano de Salud Laboral Logroño 2011 Realizado por el Área de Higiene Industrial Depósito Legal: LR-101-2011 Impreso en España - Printed

References: resolución 
 resolución 
 resolución 
 Resolución 
 RESOLUCIÓN 
 RESOLUCIÓN 
 resolución 
 resolución 
 RESOLUCIÓN