Una implementaci n num rica sencilla para simular procesos de transferencia de masa multitasa
Download
1 / 23

Una implementación numérica sencilla para simular procesos de transferencia de masa multitasa . - PowerPoint PPT Presentation


  • 147 Views
  • Uploaded on

Una implementación numérica sencilla para simular procesos de transferencia de masa multitasa. Orlando Silva, Jesús Carrera, Sireesh Kumar Abril 2008. Haggerty and Gorelick WRR (1995). Adsorción en superficie de minerales. Problema. Inclusiones de líquido. Zonas estancadas.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Una implementación numérica sencilla para simular procesos de transferencia de masa multitasa .' - elgin


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Una implementaci n num rica sencilla para simular procesos de transferencia de masa multitasa

Una implementación numérica sencilla para simular procesos de transferencia de masa multitasa.

Orlando Silva, Jesús Carrera, Sireesh Kumar

Abril 2008


Problema

Haggerty and Gorelick WRR (1995) de transferencia de masa

Adsorción en superficie de minerales

Problema

Inclusiones de líquido

Zonasestancadas

Poros/fracturasintrapartícula

Difusión a/de lajas, vainas de arcilla

Agregados de granos

INTERCAMBIO ZONA MÓVIL-ZONA INMÓVIL


Ecuaciones de transferencia de masa de conservación

1. Ecuación de flujo con MRMT

(1.1)

Ecuaciónquegobierna la transferencia con la j-ésimazonainmóvil

(1.2)

2. Ecuación de transporte de soluto con MRMT

(2.1)

Ecuaciónquegobierna la transferencia con la j-ésimazonainmóvil

(2.2)

3

3


Esto es lo que se supone que hace el presente método de transferencia de masa

Representaciónnumérica

  • 2 maneras de representar numéricamente el problema:

  • Con una malla apropiada de nodos im.

  • Eliminando im como variable de estado explícita, i.e. expresando SVim como función de SVm)


Implementación de transferencia de masa numérica

  • - El nodo m (zona móvil) está conectado a sus nodos adyacentes en la malla 1, 2 o 3D.

  • - El “nodo” im (zona inmóvil) sólo está conectado con el nodo m.

  • Geométricamente, el “nodo” im se solapa con el nodo m.

  • Numéricamente, la variable de estado (presión, concentración) en el “nodo” im puede ser resuelta de manera explícita como función de la variable de estado en la zona móvil.

  • En otras palabras, el “nodo” im puede considerarse como uno nodo de dimensión cero. (0D)

5


Ecuación de transferencia de masa de transportegeneralizada

1. Balance de la fase móvil

(3.1)

2. Transferencia de masa dede/hacia la zona inmóvil

(3.2)

3. Relaciones adicionales

(3.3)


Solución de transferencia de masa numérica

1 Se asume un esquema de Euler para la integración en el tiempo de la variable um:

2. Reemplazando en la ecuación (3.2) se obtienen N ecuaciones diferenciales de primer orden

3. Hay una solución analítica (variación de parámetros) para la condición inicial

4. Combinando lo anterior con las ecuaciones (3.2), (3.3) se obtiene el flujo a tiempo tk+q.


Contribución de transferencia de masa del proceso MRMT

Sistema numérico


Equivalencia de transferencia de masa con otrasrepresentacionesnuméricas

1. Modelosdifusionales (e.g., Harmon et al., 1989).

2. Formulaciónintegro-diferencial y función de memoria (e.g., Carrera et al., 1998).

3. Comportamientoasintótico (late-time) (e.g., Willmann et al., 2007)

g(t) sigue una ley de potencia (mg) distribución de aj, bj.

4. CTRW (e.g., Dentz & Berkowitz, 2003)

y(t): distribución de tiempos de transición

9


Módulo de transferencia de masa Fortran 90: mod_MRMT.f90

1. Tipot_immobilecontieneatributosquecaracterizan a unazonainmóvil (porosidad, coeficientes de transferencia de masa, concentraciones en la faseinmóvil, etc.).

2. Subrutinasquecalculan el aporte del proceso a las matrices numéricasD (ContriToMatrices_) y b (ContriToSink_).

  • 3. Entrada de datosmediantearchivos XML

  • Parameters fim,j, Lim,j, Dim,j, Rim,j

  • Coefficients aim,j y bim,j

  • ExpansionTermsgeometría, aim,j y bim,j, y número de términos de expansión

  • LateTime mg, t1 y t2

4. Acoplamiento con códigosestándaresprogramados en Fortran.


Acoplamiento de transferencia de masa de mod_MRMT.f90 con TRACONF


Ejemplos de transferencia de masa de verificación. Flujo radial convergente


Verificación de transferencia de masa 1. Ensayo de trazadoresconvergente (TRANSIN).

Pulso de masa de trazador radioactivo sin adsorción. La matriz o fase inmóvil se conceptualiza en forma de lajas.

10 términos de expansión


Parámetros de transferencia de masa de simulación

Tiempo total de simulación: 4 días


Curvas de transferencia de masa de llegada


Perfiles de transferencia de masa de concentración de trazador


Sensibilidad de transferencia de masa al número de términos de expansión


Verificación de transferencia de masa 2. Flujo radial hacia un pozo de bombeo (Haggerty & Gorelick, 1995).

Q* = 1 L/s, b = 20 m, rw = 0.1 m, aL = 0.1 m, fm = 0.25, Rm = ?


(a) Borden Sand. de transferencia de masa Restauración de un acuíferohomogéneo

Rm = 1.6 (PCE en acuífero Borden, Rivet & Allen-King, 2003)

50 términos de expansión por cada zona inmóvil


Evolución de transferencia de masa de la fracción de masaremanente de PCE. 500 días


(b) de transferencia de masa Casohipotético con unamezcla de procesos de transferencia de masa. Restauración de un acuíferoheterogéno

50 términos de expansión por cada zona inmóvil


Evolución de transferencia de masa de la fracción de masaremanente de PCE. 20000 días


CONCLUSIONES de transferencia de masa

1. Modelo MRMT permiterepresentarunamezclaheterogénea de procesos de transferencia de masa (partículas de acuífero, lajas de arcilla, reacciones de superficie de primer orden).

2. Presenteaproximaciónnuméricaessencilla y equivalente a otrasformulaciones de MRMT.

3. Implementación en módulo Fortran fácil de acoplar a códigosestándares de flujo y transporte.

4. Predicciónbuena de resultadosobtenidos con otraformulación y con soluciones semi-analíticas.


ad