Modelizacin y anlisis del transporte no isotrmico de

Modelización y análisis del transporte no isotérmico de compuestos orgánicos Josep Ma Gastó y Jordi Grifoll Departament d’Enginyeria Química Escola Tècnica Superior d’Enginyeria Química Universidad Rovira i Virgili (Tarragona)

Simulación del transporte de solutos • Descripción del suelo • Condiciones de contorno • Propiedades del CO • Condiciones de contorno Simulación del movimiento de agua Simulación del movimiento de solutos FLUJOS Y CONCENTRACIONES

Temperatura constante y homogénea como hipótesis habitual Evolución de los perfiles de temperatura • Descripción del suelo • Condiciones de contorno • Propiedades del CO • Condiciones de contorno Simulación del movimiento de agua Simulación del movimiento de solutos FLUJOS Y CONCENTRACIONES

Objetivos del trabajo • Elaboración de un modelo matemático del transporte no isotérmico de agua cerca de la superficie, implementar una solución numérica del mismo y analizar los efectos de los mecanismos considerados. • Elaboración de un modelo matemático del transporte no isotérmico de compuestos orgánicos incorporando el movimiento del agua, implementar una solución numérica y comparar los resultados respecto a la simulación bajo condiciones isotermas.

ECUACIONES BÁSICAS. TRANSPORTE AGUA • Balance de conservación del agua en fase líquida Flujo interfásico del L al G • Balance de conservación de la materia en fase gaseosa • Balance de conservación del vapor en fase gaseosa Flujo dispersivo y difusivo

Cálculo de flujos Ley de Darcy generalizada Dispersión hidráulica Para flujo sólo en la dirección vertical

Dispersividad longitudinal La dispersividad longitudinal en la zona no saturada varia típicamente entre 5 y 20 cm [Jury et al. , 1991]. A saturación, se ha tomado el valor experimental reportado por Biggar y Nielsen (1976) de a. LS = 7. 8 cm

TORTUOSIDAD g g A falta de medidas específicas, la tortuosidad puede estimarse mediante el denominado segundo modelo de Millington y Quirk, según prueban Jin y Jury (1996) en su recopilación de datos de laboratorio y de acuerdo con los datos de campo de Lahvis et al. (1999) que aquí se presentan.

ECUACIONES BÁSICAS. TRANSPORTE ENERGIA Hipótesis: equilibrio térmico local Mecanismos considerados • • Conducción Dispersión del vapor de agua Convección en la fase gaseosa Convección en la fase líquida

CONDICIONES DE CONTORNO TRANSPORTE AGUA • Rugosidad superficie • Velocidad del viento • Coeficiente de difusión Relación de lord Kelvin

CONDICIONES DE CONTORNO TRANSPORTE ENERGIA Intensidad de radiación en la superficie de la atmósfera • • • Dispersión Inclinación del sol latitud geográfica hora del día Absorción Emisión Nubes Vapor de agua Emisión de fondo Radiación reflejada Convección ENERGIA

DISCRETIZACIÓN • Principios de conservación aplicados a volúmenes de control Velocidad de acumulación = flujos de entrada - flujos de salida • Diferencias finitas para el cálculo de flujos 1 2

Resolución Numérica Para el paso de t a t + D t MATERIA LIQ Supuestos (t-Dt) Pl GAS ENERGIA CO VAP Tª CSM

Experimentos de campo descritos en la bibliografía • Procedimiento: ð Irrigar el suelo. ð Seguir la evolución del contenido en agua, por gravimetría, y de la temperatura a diferentes profundidades. ð Seguir la evolución de la evaporación mediante lisímetros. Jackson, R. D. (1973) Rose, C. W. (1968) Diurnal changes of soil water content during drying Water transport in soil with a daily temperature wave • • Suelo franco (Adelanto) 10 cm de irrigación inicial Suelo franco arenoso 30 cm de irrigación inicial

Contrastación con experimentos de campo (i)

Contrastación con experimentos de campo (ii) equivocat

Contrastación con experimentos de campo (iii)

Contrastación con experimentos de campo (iv)

Contrastación con experimentos de campo (v)

Procesos que rigen el transporte de agua en suelos secos (i) Profundidad Temperatura Contenido en agua Presión parcial del vapor de agua

Procesos que rigen el transporte de agua en suelos secos (ii) i-1 i • El punto de máxima presión parcial del vapor de agua delimita una zona donde el transporte es mayoritariamente en fase líquida de otra en que el transporte en fase vapor es significativo. • Los mecanismos de difusión y dispersión en fase gaseosa controlan este transporte. i+1

TRANSPORTE SOLUTOS • Mecanismos considerados ð dispersión hidrodinámica en las fases fluidas, ð convección en las fases fluidas, ð reparto entre fases de acuerdo con los coeficientes de distribución.

Simulaciones realizadas Suelo húmedo (inicio del experimento de Jackson) Concentración en la matriz porosa profundidad 0 2 escenarios 1 Suelo seco (final del experimento de Jackson) 10 cm Benceno (volátil) 2 compuestos Lindano (poco volátil) Sistema isotérmico 2 hipótesis Sistema no isotérmico

Volatilización del benceno

Volatilización del lindano

CONCLUSIONES • Se ha elaborado un modelo de transporte de agua, energía y compuestos orgánicos, válido para la zona no saturada del suelo, y se ha implementado un algoritmo para su resolución. • Se han contrastado favorablemente las simulaciones del proceso de secado cerca de la superficie. • Los flujos de volatilización son más sensibles a las variaciones de temperatura en compuestos con poca tendencia a la volatilización. • Para compuesto poco volátiles y en caso de suelos secos, la amplitud de las oscilaciones diarias de los flujos puede superar un orden de magnitud.

Condición contorno energía Is (intensidad onda corta corregida por el albedo) IL (intensidad onda larga emitida vapor agua) IIN (flujo hacia el interior del suelo) ILS (intensidad onda larga emitida por el suelo) ICONV (flujo convectivo desde el suelo)

Condición contorno energía Ángulo de incidencia Dispersión molecular Factor de turbidez Masa óptica del aire Efecto de las nubes Albedo

Conductividad térmica efectiva • Propuesta de Campbell (1994) W/(m K) es la densidad aparente. es la fracción de arcilla

Relaciones hidráulicas • Suelo utilizado por Rose (1968)

Relaciones hidráulicas Suelo utilizado por Jackson (1973)