UAMI13771
UAMI13771
Presenta:
M. en C. Francisco José Valdés Parada
Asesor:
Dr. J. Alberto Ochoa Tapia
If it is true, there is a proof
La realización de este trabajo no habría sido posible sin el financiamiento otorgado por el
Consejo Nacional de Ciencia y Tecnología. Al mismo tiempo, le estoy sumamente
agradecido a mi asesor, el Dr. J. Alberto Ochoa-Tapia, por haberme aceptado para llevar
a cabo este trabajo y por siempre compartir parte de su tiempo para discutir mis inquietudes y
apoyarme en todo momento.
Aprecio sinceramente las revisiones de las versiones preliminares del trabajo que muy
amablemente hicieron los doctores Benoît Goyeau y Michel Quintard, quienes además me
permitieron visitarlos y me propusieron importantes mejoras en aspectos específicos del trabajo. En
particular, deseo expresar mi sincero agradecimiento al profesor Stephen Whitaker, no sólo por
sus valiosas correcciones y recomendaciones, sino también por su amable disponibilidad para
platicar conmigo sobre una gran variedad de temas en los días previos a la defensa pública de este
trabajo. Asimismo, deseo agradecer a los profesores José Álvarez-Ramírez, Antonio del
Río-Portilla, Felipe López Isunza y Rosa María Velasco-Belmont, quienes aportaron
importantes comentarios y sugerencias al trabajo durante las diversas etapas de evaluación. En
particular, le estoy plenamente agradecido al profesor Alberto Soria-López por las valiosas
horas de discusión que tuvimos, lo cual me ayudó y motivó considerablemente a mejorar el trabajo.
Como siempre, estoy eternamente agradecido a todas las personas que, en un momento u otro, han
creído en mí.
Lo último que me resta es manifestar que ha sido un honor y un placer llevar a cabo este trabajo.
Resumen
Las condiciones de frontera adecuadas, ya sea de transferencia de cantidad de
Lo anterior ha despertado el interés de varios autores que, en las tres últimas décadas,
permite obtener condiciones de salto cerradas, es decir aquellas en las que el coeficiente
cerradura correspondiente.
En este trabajo, se utilizan las ideas de estos trabajos (Ochoa-Tapia y Whitaker 1995a y
Wood y col., 2000) para desarrollar una metodología que lleve a condiciones de salto
correspondientes. Los resultados se comparan con los obtenidos con otras metodologías.
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Índice
Introducción 1
Capítulo 1: Revisión bibliográfica 4
1.1.- Condición de salto empírica 5
1.2.- Condición de salto metodológica 9
1.3.- Predicción del coeficiente de salto 14
1.4.- Condiciones de salto cerradas 16
1.5.- Estimación de la permeabilidad del medio poroso homogéneo 17
1.6.-Modelos capilares 18
1.7.-Modelo basado en el radio hidráulico 19
1.8.-Modelos de arrastre para estructuras periódicas 20
Objetivos 24
Objetivo general 24
Objetivos particulares 24
Parte I: Transferencia de masa difusiva
Capítulo 2: Deducción de la condición de salto de transferencia de masa en la 26
inter-región microporo-fluido
2.1.- Planteamiento del problema 26
2.2.- Promedio volumétrico 27
2.3.- Ecuación de transferencia de masa generalizada 30
2.4.- Ecuaciones de transferencia de masa en cada región 34
2.5.- Condición de salto 36
2.6.- Relación entre 〈RA〉s y RAs 43
Capítulo 3: Planteamiento y solución del problema de cerradura
3.1.- Problemas puntual y promedio 47
3.2.- Planteamiento del problema de cerradura 48
3.3.- Forma cerrada de la condición de salto 59
3.4.- Solución del problema de cerradura 60
3.5.- Resultados 64
3.6.- Efecto de la microestructura 66
3.7.- Resumen de la metodología 71
Capítulo 4: Ecuación de transferencia de masa generalizada cerrada para 74
difusión pasiva en la frontera fluido-medio poroso
I
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
II
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Bibliografía 175
Apéndice A-1: Prueba de simetría y antisimetría de los problemas de A-1.1
cerradura
A-1.1. Problema bidimensional para la inter-región A-1.1
A-1.2. Problema de cerradura adimensional A-1.5
Apéndice A-2: Solución numérica del problema de cerradura A-2.1
A-2.1.- Planteamiento del problema A-2.1
A-2.2.- Discretización de las ecuaciones diferenciales A-2.2
A-2.3.- Esquema numérico A-2.3
A-2.4.- Independencia de parámetros numéricos y del número de A-2.10
cuadrados
Apéndice A-3: Prueba de simetría y antisimetría para el modelo de medio A-3.1
poroso centrado en el cuerpo
A-3.1.- Medio poroso centrado en el cuerpo horizontal A-3.1
A-3.2.- Medio poroso centrado en el cuerpo vertical A-3.6
Apéndice A-4: Gráficos de la variable de cerradura sγ A-4.1
III
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Índice de figuras
región ω − η .
Figura 2-4: Distancias características en la inter-región. 45
Figura 2-5: Esquema de distribución de velocidad de reacción en la inter- 45
región para el transporte dominado por a) reacción y b) difusión.
Figura 3-1: Zona representativa de la inter-región. 52
Figura 3-2: Modelo de medio poroso centrado en las caras y celdas unitarias 63
para la inter-región y la región homogénea − ω .
Figura 3-3: Coeficiente efectivo vs. φ como función de la razón de aspecto A; 65
vista bidimiensional.
Figura 3-4: Coeficiente efectivo vs. φ como función de la razón de aspecto A; 66
vista tridimiensional.
Figura 3-5: Modelo de medio poroso centrado en el cuerpo horizontal. 68
Figura 3-6: Comparación del coeficiente efectivo obtenido para un medio 68
poroso centrado en el cuerpo horizontal ( − − − ) y centrado en las
caras ( ).
Figura 3-7: Modelo de medio poroso centrado en el cuerpo vertical. 69
Figura 3-8: Comparación del coeficiente efectivo obtenido para un medio 69
poroso centrado en el cuerpo vertical ( − − − ) y centrado en las
caras ( ).
Figura 3-9: Modelo de medio microporoso bidisperso, centrado en el cuerpo 70
horizontal.
Figura 3-10: Coeficiente efectivo vs. φ para un medio poroso centrado en el 71
cuerpo horizontal con distribución bidispersa de diámetros de
IV
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( ).
Figura 4-1: Región de promediado, longitudes y regiones (homogéneas y 76
heterogéneas) características del sistema.
Figura 4-2: Zona representativa de la inter-región. 80
Figura 4-3: Campos de la componente X de la variable de cerradura para a) 83
ξ / 2h = 0.4 , b) ξ / 2h = 0.0 y c) ξ / 2h = −0.5 .
Figura 4-4: Dependencia espacial de las componentes tangencial ( ) y 84
ε γω .
Figura 4-5: Región de promediado para la inter-región y vectores unitarios 88
asociados.
Figura 4-6: Influencia de la porosidad en los perfiles de concentración como 94
función del tiempo. a) ε γω = 0.2 , b) ε γω = 0.4 , c) ε γω = 0.6 , d)
(− − − −) Modelo-II.
Figura 5-1: Flujo de un fluido sobre un medio poroso en un sistema tipo BJ. 100
Figura 5-2: Región de promediado en las regiones homogéneas y en la inter- 101
región.
Figura 5-3: Región de promediado V∞ y vectores unitarios para la inter-región. 107
V
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Figura 6-4: Dependencia espacial de los esfuerzos de Brinkman para los 135
materiales usados por BJ. La línea punteada horizontal indica el
valor de B ( 0 ) / Bmax para cada material. Bmax = max ( B ) A 2 .
Figura 6-5: Estimación del coeficiente de permeabilidad para los datos de BJ 137
usando los modelos de a) ( ) Rumpf y Gupte ( − − − ) tubos
(− − −) cilindros.
Figura 6-6: Estimación del coeficiente de permeabilidad para los datos de BJ 138
usando los modelos de a) Happel y Brenner para flujo
perpendicular ( ) y paralelo ( − − − ) , b) Larson y Higdon
Figura 6-8: Dependencia de los esfuerzos de globales con la posición para los 141
materiales estudiados por BJ γ 1 ( ) y γ 2 (− − −) .
Figura 6-9: Dependencia del coeficiente de salto con la porosidad usando 143
γ = γ 1 , γ 2 y la metodología de Goyeau y col. (2003).
Figura 6-10: Dependencia del coeficiente K −1ε βω / avs con ε βω y el tamaño de 144
poro. Usando a) γ = γ 1 y b) γ = γ 2 .
Figura 7-1: Región de promediado V, regiones y longitudes características del 151
sistema.
Figura 7-2: Región de promediado V∞ y vectores unitarios asociados. 156
VI
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Figura A-2.6: Coeficiente efectivo vs. inverso del número de cuadrados para A-2.12
A = 0.5 y φ = 3 .
Figura A-3.1: Celda unitaria para el problema en la inter-región. A-3.3
Figura A-3.2: Celda unitaria simplificada para la inter-región. A-3.6
Figura A-3.3: Celda unitaria para el problema en la inter-región. A-3.8
Figura A-3.4: Celda unitaria simplificada para la inter-región. A-3.11
Figura A-4.1: Campo de la variable de cerradura sγ para A = 0.1 , φ = 0.5 A-4.1
VII
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
VIII
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Índice de tablas
0.8.
con ε γω .
Tabla 6-1: Datos de los materiales estudiados por Beavers y Joseph (BJ) 133
Tabla 6-2: Estimación del coeficiente de salto para los materiales usados por 142
BJ
Tabla A-2.1: Dependencia del coeficiente efectivo con la tolerancia escogida A-2.10
origen.
IX
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Nomenclatura
β
Aβ tensor de segundo orden que mapea μ β ∇ v β en pβ .
B( z) esfuerzos de Brinkman, m −2 .
β
Bβ tensor de segundo orden que mapea v β en v β .
β
Bω tensor de segundo orden que mapea v β
ω
en vω .
β
bβ vector que mapea μ β v β en pβ , m −1 .
β
bω vector que mapea μ β v β ω
en pβ ,ω , m −1 .
β
Cβ tensor de tercer orden que mapea ∇ v β en v β , m.
γ β
c Aγ j
, c Aβ j
promedio intrínseco de la especie A en la región homogénea j , j = η , ω ,
mol/m3 .
γ β
c Aγ , cAβ promedio intrínseco de la especie A, mol/m3 .
j = β , σ , mol/m3 .
X
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Dγ , Dβ difusión molecular, m 2 /s .
I tensor identidad.
K desviaciones espaciales de K, m / s.
k desviaciones espaciales de k, m / s.
k ηω
promedio en área de k, m / s.
XI
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
CK y LH, m.
j =κ o σ .
β
pβ promedio intrínseco de la presión, N/m 2 .
β
pβ
j
promedio intrínseco de la presión en la región − j ; j = ω ,η , N/m 2 .
pβ
j
promedio superficial de la presión en la región − j ; j = ω ,η , N/m 2 .
r vector de posición, m.
γ
sγ variable de cerradura que mapea cAγ en c Aγ .
γ
sω variable de cerradura que mapea c Aγ
ω
en c Aω .
T s
esfuerzo superficial, N/m 2 .
Tβ
s
esfuerzo global, N/m 2 .
t tiempo, s.
XII
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
V región de promediado.
j = β , γ ,σ ,κ .
β
vβ promedio intrínseco del vector de velocidad, m/s .
β
vβ
j
promedio intrínseco del vector de velocidad en la región − j ; j = ω ,η ,
m/s .
m/s .
vβ j
componente tangencial del promedio superficial del vector de velocidad
en la región − j ; j = ω ,η , m/s .
m.
j = β , σ , m.
XIII
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Símbolos griegos
Subíndices
A relativo a la especie A.
C relativo a la concentración.
v relativo a la velocidad.
XIV
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Introducción
El estudio de los fenómenos de transporte en sistemas multifásicos, y en particular en
XIX (una descripción detallada de la evolución de este tema hasta nuestros días se
Ruthven, 1985; Dalvie y col., 1990; Goto y McCoy, 2000). Los resultados de este
sistemas de tres fases (como son las emulsiones dobles) así como en biopelículas, entre
Karniadakis y Beskok, 2005) las cuales pueden derivarse a partir de la metodología aquí
propuesta.
1
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
respecto, Wood y col. (2000) propusieron una metodología para desarrollar condiciones
Wood y col. (2000) para así obtener condiciones de frontera cerradas para el transporte
sistema, es decir tanto las ecuaciones de medio efectivo como sus correspondientes
los trabajos recientes que se han dedicado a obtener estimaciones aproximadas del
En este capítulo se incluye, además, una revisión de algunos de los modelos reportados
plantean los casos concretos a tratar. Estos, se orientan a tres situaciones: transferencia
masa.
2
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Con base en estos problemas, el resto del trabajo se presenta en las siguientes tres
3
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 1:
Revisión Bibliográfica
Este capítulo está dedicado a la revisión de los antecedentes del presente trabajo. Para
ello, primeramente se expone una breve revisión de los trabajos más relevantes
con la revisión del trabajo de Wood y col. (2000) donde se desarrollan condiciones de
salto cerradas, es decir, que no dependen de parámetros ajustables, las cuales involucran
clasifican en dos clases, aquellos en los que se considera el flujo a través de conductos
tortuosos y por otro lado, los que consideran el flujo alrededor de objetos sólidos
inmersos en un fluido.
4
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
tema de interés en el campo de la Ingeniería Química desde hace más de medio siglo.
Sin embargo, como se mencionó en la Introducción, una buena parte de los trabajos se
estudiar los efectos interfaciales que ocurren al hacer pasar un fluido sobre un medio
poroso, saturado con el mismo fluido, mediante una caída de presión determinada.
Consideraron dos tipos de materiales, uno granular (aloxita) y otro con una distribución
externos de los poros superficiales) no necesariamente debía ser continua; por lo que
dvβ
dy y =O+
=
α
K βω
(v
β y =O
+
)
− vω , (1.1)
donde α es un parámetro adimensional que se fijó de tal forma que se logren reproducir
del material en la vecindad con el fluido más no de la viscosidad del fluido; los valores
raíz cuadrada de la permeabilidad del medio poroso, K βω . Relacionado con esto último,
5
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
determinaron experimentalmente que el espesor es más bien del orden del diámetro de
grano.
rectángulos paralelos, con una caída de presión despreciable. De esta forma, fueron
constante.
Más tarde, Saffman (1971) utilizó un análisis de capa límite para encontrar una
medio poroso, además no reporta datos que comparen su metodología con los resultados
estar sujeta a restricciones de escala. Esta idea será retomada más adelante.
Por otro lado Jones (1973) sugirió agregar un término extra a la condición de BJ con el
dvβ
dy y =O+
+
dvβ
dx y
=
α
K βω
(v
β y =O
+
)
− vω . (1.2)
6
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Por otro lado, Neale y Nader (1974) utilizando la ecuación de Darcy-Brinkman en lugar
análogos a los de BJ. La razón para este cambio recae en el argumento que la ley de
Darcy es incompatible con la existencia de una capa límite en el medio poroso dado que
debe incluirse al extrapolar la ley de Darcy a la frontera. Bajo esta línea de pensamiento,
hecho, de acuerdo a Whitaker (1999) la viscosidad efectiva está dada explícitamente por
1988; Givler y Altobelli, 1994; Martys y col., 1994; Starov y Zhdanov, 2001) han
Altobelli (1994), la razón μeff / μ β puede variar entre 5.1 y 10.9 dependiendo de ε β . El
uso de este tipo de modelo se recomienda para flujo, de bajo número de Reynolds, en
Sobre lo anterior, Haber y Mauri (1983) usando ideas muy similares a las de Saffman
(1971) demuestran la diferencia que existe entre usar la ecuación de Brinkman y la ley
vuelve despreciable, recuperándose así la ley de Darcy. Entre sus conclusiones recalcan
7
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
“algo más que la intuición física”, que permitan que los resultados obtenidos con la ley
A partir del trabajo de Whitaker (1969) sobre la deducción de la ley de Darcy; Ross
Marle, 1967; Slattery, 1967; Whitaker, 1967, Bear, 1972) a la frontera medio poroso-
fluido, obteniendo una ecuación general que es válida tanto en la inter-región como en
las regiones homogéneas. Para el caso de un medio isótropo, obtiene una expresión del
una ecuación que sea válida en todo el dominio fue usada posteriormente por otros
Por otro lado, Prat (1989, 1990, 1992) propuso una metodología para la determinación
representativa con la cual hace el promediado y de esta manera obtiene los valores en la
frontera. Para probar la pertinencia del método, se calculan las desviaciones entre las
8
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
El modelo de celda unitaria y otras de las ideas expresadas en sus trabajos serán
cambiar la estructura del medio poroso (por medio de distintos arreglos de cilindros) los
Motivado por el trabajo previo de Stone (1990), Whitaker (1992) desarrolló de manera
nivel de escala puntual usando una metodología que puede extenderse para desarrollar
condiciones de salto a otros niveles de escala. Entre las aplicaciones que tiene la
9
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
desarrollar condiciones de salto cuando las ecuaciones promedio sean idénticas a las
ecuaciones puntuales en las dos fases. Una discusión detallada de los fenómenos de
dominio.
10
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Se estimó que dicho parámetro tiene orden de magnitud de la unidad. De esta forma,
de la velocidad y del flux*. Notaron que, en este caso, la comparación con datos
experimentales deja mucho que desear a pesar de ser un modelo más complicado que el
que utiliza la condición de salto. Para estimar los cambios de la porosidad con la
son relativamente suaves por lo que pudieron representarlos por medio de polinomios
Una versión del trabajo de Ochoa-Tapia y Whitaker (1995a) que incluye los efectos
inerciales en la condición de salto fue presentada más tarde por los mismos autores
esfuerzos inerciales.
*
A lo largo del texto se usará la palabra flux para denotar densidad de flujo de alguna propiedad.
11
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Esta metodología fue aplicada más tarde para obtener la condición de salto de
1995) y con reacción química (Soria y col., 1996). Los resultados obtenidos fueron
medio poroso y la ecuación de Stokes para el fluido. Notó que los perfiles de velocidad
Srivastava (1999), para el flujo en medios porosos inducido por un disco que gira cerca
fluido del medio poroso a lo largo de una capa límite, cuyo espesor depende de las
propiedades del medio poroso y no del parámetro β . Por otro lado, recientemente
esfuerzo cortante en el flujo entre esferas no sólo se ven afectados por cambios en la
permeabilidad sino que también el coeficiente de salto ejerce una influencia importante.
12
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Más adelante, Hager y Whitaker (2000), desarrollaron las condiciones de salto de masa
y energía entre un medio poroso y dos fluidos. Debido a que la suposición de equilibrio
local térmico falla en las fronteras entre estas tres fases, un frente de evaporación tuvo
momento entre un medio poroso y un fluido fue presentada por Alazmi y Vafai (2001).
Encontraron que los modelos 1 y 2 dan resultados similares a los del modelo 5 para
valores altos de la porosidad, mientras que para valores bajos, el modelo 2 se aproxima
a los modelos 3 y 4. Por lo que propusieron que los modelos 1 y 5 son muy parecidos, al
igual que los modelos 3 y 4, mientras que el modelo 2 cae entre estos dos pares de
y col. (2003) considerando los efectos de adsorción entre un medio poroso y un fluido.
13
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
sea por técnicas analíticas o numéricas. Entre ellos, cabe mencionar a Goyeau y col.
(2003) quienes desarrollaron una metodología para obtener una expresión semi-analítica
dominios, en el primero se utiliza una sola ecuación para describir el transporte en los
dos medios y en el segundo se usan las ecuaciones de Brinkman y de Stokes, junto con
la información suficiente para obtener una expresión analítica del parámetro ajustable.
Sin embargo, a pesar de que sus resultados son muy similares a los que se habían
propuesto ad hoc para ajustar los datos experimentales, admiten que su trabajo es un
paso intermedio hacia la estimación de dicho parámetro; ya que para ello es necesario
metodología se utilizan más adelante en el trabajo (Capítulo 6). Esta metodología fue
de curvas de las soluciones del problema en uno y dos dominios, sólo que para esto
Por otro lado, Min y Kim (2005) estudiaron los problemas de transferencia de momento
metodología no considera los efectos de la geometría del medio poroso, por ejemplo, se
metodología predice valores negativos y cercanos a cero del coeficiente de salto para
menor de 0.2. Estos resultados no concuerdan con los reportados por Ochoa-Tapia y
Whitaker (1995b) donde los valores de dicho parámetro están entre −1 y 1.47 para los
trabajos de investigación en diversos campos, son una muestra de que el tema de las
de lo anterior, son claros tanto el interés como la necesidad de contar con una
comentan algunos detalles del trabajo de Wood y col. (2000), el cual permite desarrollar
15
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Una de las contribuciones más importantes para el presente trabajo es la de Wood y col.
cambios espaciales.
coeficiente en la misma.
16
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
por ejemplo, considerando una reacción química con cinética no lineal (Bletzacker y
col., 2003), o bien entre un medio poroso y un fluido como se hará en este trabajo.
Por otro lado, Lichter y col. (2003) usando análisis estocástico determinaron la
dependencia del coeficiente efectivo de velocidad de reacción, definido por Wood y col.
(2000), con el tiempo en un reactor batch. Encontraron que los valores del coeficiente
Hasta este punto se llevará la revisión de trabajos relacionados con condiciones de salto,
β
d pβ
μ β vβ ω
= − K βω ω
, (1.4)
dx
donde la permeabilidad ( K βω ) es la medida de la conductividad de un medio poroso al
flujo de fluidos. Sin embargo, dada la complejidad para modelar los fenómenos de
clasifican en dos tipos, los modelos capilares donde el flujo a través de la matriz porosa
es visto como el flujo a través de conductos tortuosos y los modelos de arrastre donde
17
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
porosidades bajas e intermedias, los modelos capilares son más apropiados, mientras
que para altas porosidades la mejor opción se considera que son los modelos de arrastre
modelos mencionados.
1.6.-Modelos capilares.
desarrollado. De esta forma, suponiendo que el medio poroso está constituido por n
β
l 2 d pβ
vβ =− σ ω
. (1.5)
ω, p 32 μ β dx
β
nπ lσ4 d pβ
vβ = ε βω vβ =− ω
. (1.6)
ω ω, p 128μ β dx
nπ lσ4 ε βω lσ
2
K βω = = . (1.7)
128 32
18
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
(Carman, 1937) se define como la razón de cuatro veces el volumen vacío entre el área
superficial, o bien
4ε βω
lβ = . (1.8)
av ,ω
En este caso,
av ,ω =
4
dc
(1 − ε βω ) , (1.9)
factor de forma k0 ,
β
l β2 d pβ
vβ =− ω
, (1.10)
ω, p 16τ k0 μ β dx
como sigue,
ε βω
vβ = vβ . (1.11)
ω τ ω, p
ε βω l β2 ε βω
3
K βω = = , (1.12)
16k K k K av2,ω
19
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
6 (1 − ε βω )
lσ = . (1.13)
av ,ω
3
De manera que para el caso de cilindros lσ = d c y para esferas lσ = d s , donde d s es el
2
modificada a,
ε βω
3 2
lσ
K βω = . (1.14)
180 (1 − ε βω )
2
Para un lecho empacado con partículas esféricas con una baja distribución de tamaños,
Rumpf y Gupte (1971) muestran que la permeabilidad se puede calcular por medio de,
ε βω
5.5 2
lσ
K βω = , (1.15)
5.6
la cual ofrece un buen ajuste con resultados experimentales para valores de porosidad
Las ideas detrás de estos modelos son muy parecidas a las empleadas en la solución del
problema de cerradura (Eidsath y col., 1983; Barrère y col. 1992), es decir, se resuelven
en una cela unitaria periódica las ecuaciones de Navier-Stokes para obtener de esta
perpendicular a los cilindros. Note que δ puede considerarse como el diámetro de poro.
el perfil de velocidad, con el que se puede obtener la velocidad de Darcy, de manera que
ε βω
3 2
lσ
K βω = , (1.16)
16kk (1 − ε βω )
2
2ε βω
3
kk = , (1.17)
⎡ ⎛ ⎞ 2⎤
(1 − ε βω ) ⎢2 ln ⎜⎜ 1 −1ε ⎟⎟ − 3 + 4 (1 − ε βω ) − (1 − ε βω ) ⎥
⎣⎢ ⎝ βω ⎠ ⎦⎥
δ 2 + 2 lσ δ
ε βω = . (1.18)
( lσ + δ )
2
(l +δ ) ⎡ ⎛ l + δ ⎞ 1 ( lσ + δ )4 − lσ4 ⎤
2
K βω = σ ⎢ln ⎜ σ ⎟− ⎥. (1.19)
⎢⎣ ⎝ lσ ⎠ 2 ( lσ + δ ) + lσ ⎥⎦
4 4
16
Es importante mencionar que al usar otros arreglos de cilindros como los que
21
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Figura 1-1: Modelo periódico de medio poroso y celda unitaria usando cilindros
Sin embargo, debido a la baja capacidad de estos modelos simplificados para reproducir
valores bajos de porosidad. Por ello, se han resuelto numéricamente las ecuaciones de
Por otro lado, Larson y Higdon (1989) consideraron un modelo de medio poroso
constituido por esferas, el cual se muestra junto con la correspondiente celda unitaria en
22
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
K βω ε βω
3
=A , (1.20)
(1 − ε βω )
2
l2
l 1
= . (1.21)
(1 − ε )
1/ 3
lσ βω
Para finalizar este capítulo se presenta la siguiente expresión propuesta por Kaviany
(1999) para ajustar los resultados numéricos de Sahraoui y Kaviany (1992) para arreglos
de cilindros,
K βω π ε βω
5.1
23
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Objetivos
Objetivo general:
Objetivos particulares:
primer orden.
comparar los resultados con los que se obtienen al usar una condición de salto
fluido.
24
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Parte I:
25
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 2:
inter-región microporo-fluido.
En este capítulo se desarrollarán condiciones de salto cerradas (es decir, sin involucrar
medio poroso compuesto de microporos. Para ubicar el nivel de escala al que se dirige
el análisis, en la Figura 2-1 se muestran cinco niveles que van desde la escala reactor
por Wood y col. (2000). De esta manera, el problema que se desea estudiar se ubica en
el nivel de Escala III. Este capítulo está dedicado a la deducción de las condiciones de
salto cerradas para este problema, mientras que en el Capítulo 3 se presenta y resuelve
resultante.
Problema puntual.
Las condiciones de salto puntual y promedio para la Escala V fueron reportadas por
Wood y col. (2000) suponiendo que la reacción química sólo tiene lugar en la interfase
sólido-fluido. Para este análisis, se referirá a la fase sólida como la fase −κ mientras que
Escala V son
26
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂c Aγ
= ∇ i(Dγ ∇c Aγ ) , en la fase − γ , (2.1)
∂t
Elemental Representativo o V.E.R.; Bear, 1972) en cada punto del sistema de la Escala
*
A lo largo de este documento, se usarán las fuentes Times New Roman cursiva y Shelley Allegro
cursiva para referirse a espacios (o superficies) geométricos y sus volúmenes (o áreas), respectivamente.
27
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
V =V γ ( x ) +V κ ( x ) , (2.3)
Por otro lado, en este trabajo se maneja el concepto de región homogénea de la misma
forma que Ochoa-Tapia y Whitaker (1995a); es decir, es una porción del sistema
está influenciada por los rápidos cambios en la estructura que ocurren en la frontera. En
otras palabras, una región homogénea es aquella porción del sistema donde las
inter-región o frontera macroscópica, como la porción del sistema que separa dos
por las fases −κ y −γ , siendo esta última la única fase presente en la región −η ; en
partícula catalítica.
La región de promediado (V ) debe ser lo suficientemente grande para que las no-
sigue,
lγ r0 L. (2.4)
28
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1
ψγ
x
=
V ∫
Vγ ( x )
ψγ x+y γ
dV , (2.5)
1
ψγ =
V ∫
Vγ ( x )
ψ γ dV . (2.6)
γ 1
V γ ( x ) Vγ∫( x ) γ
ψγ = ψ dV , (2.7)
29
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
dado que la integral en el lado derecho de las ecuaciones (2.6) y (2.7) es la misma,
1 Vγ (x)
εγ ( x) = ∫ dV = . (2.9)
V Vγ ( x )
V
Dado que en este trabajo se estudian los cambios espaciales al pasar de una región
consiste en obtener una ecuación de medio efectivo que sea válida en cualquier punto
∂c Aγ
= ∇ i(Dγ ∇cAγ ) . (2.10)
∂t
(2.8), como:
γ
∂c Aγ ∂ cAγ ∂ c Aγ
= = εγ ( x) . (2.11)
∂t ∂t ∂t
30
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Por otro lado, mediante el teorema del promediado espacial (Howes y Whitaker, 1985),
∇ i(Dγ ∇c Aγ ) = ∇ i Dγ ∇c Aγ +
1
V ∫
Aγκ ( x )
nγκ iDγ ∇c Aγ dA . (2.12)
expresar el primer término del lado derecho de la ecuación anterior, con una segunda
⎡ ⎛ 1 ⎞⎤
∇ i Dγ ∇cAγ = ∇ i ⎢Dγ
⎢
⎜ ∇ c Aγ +
⎜ V ∫ γκ Aγ ⎟⎟⎥⎥ .
n c dA (2.13)
⎣ ⎝ Aγκ ( x ) ⎠⎦
⎡ ⎛ ⎞⎤ 1
∇ i(Dγ ∇c Aγ ) = ∇ i ⎢Dγ
1
⎢
⎜ ∇ c Aγ +
⎜ V ∫ γκ Aγ ⎟⎟⎥⎥ −V
n c dA ∫ kcAγ dA . (2.14)
⎣ ⎝ Aγκ ( x ) ⎠⎦ Aγκ ( x )
γ
c Aγ x + yγ
= c Aγ + c Aγ x + yγ
. (2.15)
x + yγ
⎡ ⎛ε x ∇ c
∇ i(Dγ ∇c Aγ ) = ∇ i ⎢Dγ
γ γ
⎜ γ( ) Aγ + c Aγ ∇ε γ ( x )
⎣ ⎝ x x
⎞⎤ 1 (2.16)
1 γ 1
+
V ∫ nγκ c Aγ
x + yγ
dA +
V Aγκ∫( x )
n γκ c Aγ d A ⎟⎥ −
⎟⎥ V ∫ kc Aγ dA.
Aγκ ( x ) ⎠⎦ Aγκ ( x )
31
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ 1 γ
cAγ
x
∇ε γ = −
V ∫
Aγκ ( x )
nγκ c Aγ
x
dA . (2.17)
1
dado que ∇ε γ = −
V ∫
Aγκ ( x )
nγκ dA . Por lo que la ecuación (2.16) se escribe ahora como,
(
∇ i(Dγ ∇c Aγ ) = ∇ i ε γ ( x ) D ( x )i∇ c Aγ
γ
) −V1 ∫
Aγκ ( x )
kc Aγ dA , (2.18)
γ ⎡ γ 1 ⎛ γ γ ⎞
D ( x )i∇ c Aγ
V γ ( x ) ∫( ) γκ ⎝
= Dγ ⎢∇ cAγ + n ⎜ c Aγ − cAγ ⎟ dA
⎢⎣ Aγκ x
x + yγ x⎠
(2.19)
1 ⎤
V γ ( x ) Aγκ∫( x ) γκ Aγ ⎥
+ n c dA ⎥.
⎦
y como consecuencia
γ 1 ⎛ γ γ ⎞
∇ c Aγ ∫ nγκ ⎜ c Aγ
V γ ( x h ) Aγκ ( xh ) ⎝ xh + yγ
− c Aγ ⎟ dA,
xh ⎠ (2.21)
en las regiones homogéneas,
Por otro lado, debe remarcarse que el uso de las restricciones de escala dadas en la
ecuación (2.20), está limitado a las regiones homogéneas dado que existe cierta
cuenta la transición de las propiedades del tensor desde el fluido homogéneo (Dγ I) a
32
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
para estimar esta transición, se requiere conocer, entre otras propiedades, los cambios
detalladamente en el Capítulo 4.
Por otro lado, al despreciar los cambios espaciales del coeficiente de reacción k en
(
∇ i(Dγ ∇c Aγ ) = ∇ i ε γ ( x ) D ( x )i∇ c Aγ
γ
) − ka ( x) c
v Aγ γκ
, (2.22)
1
A γκ ( x ) Aγκ∫( x )
cAγ γκ
= cAγ dA , (2.23)
1 A γκ ( x )
av ( x ) = ∫ dA = . (2.24)
V Aγκ ( x )
V
sustituir ahora las ecuaciones (2.11) y (2.22) en la ecuación (2.10) resulta finalmente,
γ
( ) − ka ( x) c
∂ c Aγ γ
εγ ( x) = ∇ i ε γ ( x ) D ( x )i∇ cAγ Aγ . (2.25)
∂t
v γκ
estudiado en la Escala III, la cual carece de restricciones de escala y es, por lo tanto,
válida en todo el dominio. A este tipo de modelos se les suele llamar modelos de un solo
dominio (one domain approach) (Goyeau y col., 2003). Debe notarse que los
donde ε γω y avω son los valores de la porosidad y el área superficial por unidad de
33
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
No debe perderse de vista que la ecuación (2.25) no está cerrada y por lo tanto no es
γ
posible utilizarla para obtener el campo de c Aγ . Esta dificultad puede superarse, ya
sea estimando las propiedades necesarias para calcular los cambios espaciales de D ( x )
Región homogénea- ω .
lγ LC , (2.27)
lγ r0 , (2.28)
se cumple que,
γ
cAγω c Aγ . (2.29)
ω x +y
ω γ
γ
promedio. De acuerdo a la desigualdad (2.29), es posible remplazar a c Aγ γκ
por c Aγ
ω
34
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
ε γ ( xω ) D ( xω )i∇ cAγ = ε γωDω i∇ cAγ ω
, (2.30)
xω
donde,
γ ⎛ γ 1 ⎞
V γ ( xω ) Aγκ ∫( xω ) γκ Aγω ⎟
Dω i∇ c Aγ = Dγ ⎜ ∇ c Aγ + n c d A ⎟. (2.31)
ω ⎜ ω
⎝ ⎠
fue llevada a cabo por Ryan y col. (1981), obteniendo un buen ajuste con datos
experimentales.
( )−a
∂ c Aγ γ γ
ε γω ω
= ∇ i ε γωDω i∇ c Aγ vω k cAγ , en la región homogénea − ω . (2.32)
∂t ω ω
Región homogénea-η .
En este caso, ( x → xη ) las variables promedio son iguales a sus valores puntuales
ε γ ( xη ) D ( xη )i∇ cAγ
γ γ
= Dγ ∇ c Aγ . (2.34)
η η
35
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ
∂ cAγ
∂t
η
(
= ∇ i Dγ ∇ c Aγ
γ
η ) , en la región homogénea −η . (2.35)
Cabe mencionar que las ecuaciones (2.32) y (2.35) son análogas a las ecuaciones (2-21)
y (2-29) de Whitaker (1983); quién dedujo el modelo de una sola ecuación para un
A π RT
Dω = Dω I = I, (2.36)
avω 8M A
orden que depende sólo de la geometría del medio poroso. Los detalles que llevan a la
Como se mencionó anteriormente, las restricciones de escala, bajo las cuales las
ecuaciones (2.32) y (2.35) son válidas, pueden no ser satisfechas en la inter-región. Por
Con estas condiciones de frontera, se contará con un modelo completo para predecir el
campo de concentración que consiste en dos ecuaciones de medio efectivo válidas en las
modelos, algunos autores les llaman modelos de dos dominios (two domain approach)
V ∞ =V ω +V η , (2.37)
descomponerse en,
A ∞ = A ω +Aη , (2.38)
movimiento; sin embargo, existen otros casos donde el flujo del fluido debe
Capítulo 7).
γ
∂ c Aγ
Región−η: ∫
Vη
∂t
η
dV = ∫ nη i(Dγ ∇
Aη
c Aγ
γ
η ) dA + ∫ n i(D ∇ c ) dA ,
*
Aωη
ηω γ Aγ η
γ
(2.40)
Todo el dominio:
γ γ
∫ nω i(ε γ ( x ) D ( x )i∇ ) dA
∂ cAγ ∂ cAγ γ
∫ εγ ( x)
Vω
∂t
dV + ∫ ε γ ( x )
Vη
∂t
dV =
Aω
c Aγ
(2.41)
Aη
(
+ ∫ nη i ε γ ( x ) D ( x )i∇ c Aγ
γ
) dA − ∫ ka ( x ) c
Vω
v Aγ γκ
dV ,
37
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
∂ c Aγ ∂ c Aγ
∫ ε γω dV + ∫
ω η
dV
Vω
∂t Vη
∂t
= ∫ nω i(ε γωDω i∇
Aω
cAγ
γ
ω ) dA + ∫ n i(D ∇ c ) dA
Aη
η γ Aγ η
γ
(2.42)
+ ∫*
Aωη
(
nωη i ε γωDω i∇ c Aγ
γ
ω
−D ∇ c γ ) dA − ∫ a k c
Aγ η
γ
Vω
vω Aγ
γ
ω
dV .
expresión
γ
⎡ ∂ c Aγ
γ
∂ c Aγ
γ
⎤ ⎡ ∂ c Aγ
γ
∂ c Aγ ⎤
⎥ dV + ⎢ε γ ( x ) ⎥ dV
∫V ⎢⎢ε γ ( x ) ∂t ∫
η
− ε γω ω
−
∂t ⎥ ⎢ ∂t ∂t ⎥
ω
⎣ ⎦ Vη
⎢⎣ ⎦⎥
= ∫ nω i(ε γ ( x ) D ( x )i∇
Aω
c Aγ
γ
− ε γωDω i∇ cAγ
γ
ω ) dA
(
+ ∫ nη i ε γ ( x ) D ( x )i∇ cAγ
Aη
γ
− Dγ ∇ cAγ ) dA
γ
η
(2.43)
38
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
∂ cAγ ⎡ ∂ cAγ
γ
∂ c Aγ ω ⎤
∫ εγ s
dA = ∫ ⎢ε γ ( x ) − ε γω ⎥ dV
s
∂t ⎢ ∂ t ∂ t ⎥
⎣ ⎦
* Vω
Aωη
γ
(2.44)
⎡ ∂ c Aγ
γ
∂ c Aγ η ⎤
⎢
+ ∫ εγ ( x) − ⎥ dV ,
⎢ ∂t ∂ t ⎥
⎣⎢ ⎦⎥
Vη
∫ n iε γ D i∇
C
s s s s c Aγ
γ
s
dσ = ∫ nω i(ε γ ( x ) D ( x )i∇
Aω
c Aγ
γ
− ε γωDω i∇ cAγ
γ
ω ) dA
( ) dA,
γ γ
(2.45)
+ ∫ nη i ε γ ( x ) D ( x )i∇ cAγ − Dγ ∇ c Aγ η
Aη
1978; Slattery, 1980; Ochoa-Tapia y col., 1993; Slattery, 2007) en la ecuación (2.45) se
obtiene que
∫ ∇ iε γ D i∇
*
Aωη
s s s s c Aγ
γ
s
dA = ∫ nω i(ε γ ( x ) D ( x )i∇
Aω
c Aγ
γ
− ε γωDω i∇ c Aγ
γ
ω ) dA
(2.46)
(
+ ∫ nη i ε γ ( x ) D ( x )i∇ c Aγ
Aη
γ
− Dγ ∇ c Aγ
γ
η ) dA.
En este momento es conveniente introducir un tensor de proyección P de acuerdo a
(Whitaker, 1992),
términos de P y ∇ , como
39
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
entonces,
1
RA s
=
Aηω ∫R As dA = RAs ηω
= kc Aγ ηω
, (2.50)
Aηω
(2.43), se obtiene, tras agrupar todos los términos en una sola integral y tomando en
*
cuenta que la superficie Aωη es arbitraria (lo que permite igualar a cero el integrando),
la siguiente expresión
)
∂ c Aγ
εγ s
∂t
s
(
− ∇ s i ε γ sDs i∇ s c Aγ
γ
s
=
exceso de acumulación superficial exceso de transporte por difusión interna
(2.51)
(
−nωη i ε γωDω i∇ cAγ
γ
ω
− Dγ ∇ cAγ
γ
η )− kcAγ ηω
.
transporte difusivo neto hacia la inter-región exceso en la
tasa de reacción
En este momento conviene notar que, si los cambios espaciales de los integrandos en los
dos términos del lado derecho de la ecuación (2.46) se suponen despreciables en las
∫ ∇ iε γ D i∇
*
Aωη
s s s s c Aγ
γ
s
dA = ∫
Aω ( xω )
(
nω i ε γ ( x ) D ( x )i∇ cAγ
γ
− ε γωDω i∇ c Aγ
γ
ω ) dA
( ) dA,
γ γ
(2.52)
+ ∫ nη i ε γ ( x ) D ( x )i∇ c Aγ − Dγ ∇ c Aγ η
( )
Aη xη
las restricciones de escala bajo las cuales la ETMG y las propiedades promedio
40
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
lado derecho de la ecuación (2.52) se puede suponer que es mucho menor que la unidad
( )−
∂ c Aγ γ γ
εγ s s
= −nωη i ε γωDω i∇ c Aγ − Dγ ∇ c Aγ kc Aγ . (2.53)
∂t ω η ηω
Por otro lado, al llevar a cabo un análisis de orden de magnitud en la ecuación (2.44) se
obtiene que
) t VA
∂ c Aγ
(
⎡ γ γ ⎤
εγ s s
= O ⎢ ε γ Δ t c Aγ − ε γω Δ t cAγ ω
⎥
∂t ⎣⎢
ω * *
⎥
ωη ⎦
(2.54)
⎡
⎢⎣
(
+O ⎢ ε γ Δ t c Aγ
γ
− Δ t c Aγ
γ
η ) Vη ⎤
* ⎥
t * Aωη ⎥⎦
,
γ γ
donde Δ t ψ γ representan los cambios en ψ γ que ocurren durante un tiempo
característico t * , el cual, como se mostrará más adelante (Sección 3.2), está restringido
a t * = O ( LC LC1 / Dγ ) . Por otro lado, de acuerdo a la Figura 2-3, se pueden proponer los
siguientes estimados
⎛ V ⎞ ⎛ BAωη
*
⎞
O ⎜ *j
⎜A ⎟⎟ = O ⎜⎜ * ⎟⎟ = O ( B ) , j = ω ,η , (2.55)
⎝ ωη ⎠ ⎝ Aωη ⎠
donde, B es la mitad del espesor de la zona de transición (Figura 2-4). De esta forma, la
γ
∂ c Aγ γ γ γ γ
εγ s s
= O ⎡⎢ε γ c Aγ − ε γω c Aγ ⎤ B + O ⎡ε c − cAγ ⎤B. (2.56)
∂t ⎣ ω ⎥⎦ t * ⎢⎣ γ Aγ η ⎥⎦ t *
hacia la inter-región es
(
nωη i ε γωDω i∇ c Aγ
γ
ω
− Dγ ∇ c Aγ
γ
η ) = O ⎛⎜⎜⎝ ε DD
γω
γ
ω
cAγ
γ
ω
− c Aγ
γ⎞ Dγ
⎟
η ⎟ L
⎠ C
. (2.57)
41
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
⎛ ⎛ε D γ γ ⎞ ⎞
⎜ O ⎜ γω ω c Aγ ω − c Aγ η ⎟ ⎟
⎜ ⎟
⎜ ⎝ Dγ ⎟ Dγ t
*
⎠
( ) ( )
⎜ γ γ γ γ ⎟ BL 1, (2.58)
⎜ O ε γ cAγ − ε γω c Aγ ω + O ε γ c Aγ − cAγ η ⎟ C
⎜ ⎟
⎝ ⎠
⎛ε D γ γ ⎞
O ⎜ γω ω c Aγ − c Aγ ⎟⎟
⎜ D ω η
⎝ γ ⎠ = O (1) ,
( ) + O (ε )
γ γ γ γ
(2.59)
O ε γ c Aγ − ε γω cAγ ω γ c Aγ − cAγ η
Dγ t *
la restricción de (2.58) puede expresarse, tomando en cuenta que 1 , como
BLC
LC1
1 , proceso cuasiestacionario. (2.60)
B
γ
asociada con los cambios de ∇ cAγ (la cual puede suponerse del orden de la longitud
puede considerarse despreciable respecto a los términos del lado derecho de la ecuación
(
−nηω i Dγ ∇ c Aγ η
γ
) + n i(ε D i∇ c ) = kc
ηω γ ω Aγ
γ
ω Aγ ηω , en la superficie divisoria , (2.61)
Conviene remarcar que el uso de este resultado sólo puede considerarse aceptable
42
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
sin embargo, este análisis aún requiere más discusión y se presentará en un trabajo
Hasta este punto, la condición de salto dada por la ecuación (2.61) no está cerrada, pues
medibles en las regiones homogéneas. Para ello, es necesario el desarrollo del problema
superficial como,
γ
∫
*
RA s dA = ∫
V∞
RAγ γκ
dV − ∫ RAγ
Vω
ω
dV , (2.63)
Aωη
definiciones,
RAγ γκ
= av ( x ) k c Aγ γκ
, (2.64)
γ γ
y RAγ = avω k c Aγ . (2.65)
ω ω
+B 0 γ
RA =∫ RAγ γκ
dz − ∫ RAγ ω
dz . (2.66)
s −B −B
43
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
medio poroso hasta la posición del centroide del volumen de promediado donde se
encuentra alguna región homogénea. Por simplicidad, esta distancia se supuso igual
para ambas regiones, pero el análisis puede extenderse sin mayores problemas a
situaciones cuando no lo son. A diferencia del caso estudiado por Wood y col. (2000),
en este caso la reacción química tiene lugar tanto en la superficie divisoria como dentro
del medio poroso. A la escala puntual, la reacción química tiene lugar (dondequiera que
exista superficie catalítica) a lo largo de una distancia igual a 2Δ* y, para el caso de la
frontera, se incluye en la Figura 2-4 una comparación con las otras longitudes
características del sistema. Cabe mencionar que la distancia 2Δ* será igual a 2Δ (Wood
y col., 2000) si los efectos de curvatura locales pueden suponerse despreciables, como
De acuerdo a la Figura 2-4, es posible descomponer el primer sumando del lado derecho
+B γ − r0 + r0 +B
∫−B
RAγ dz = ∫
−B
RAγ γκ
dz + ∫
− r0
RAγ γκ
dz + ∫
+ r0
RAγ γκ
dz . (2.67)
44
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ
RAγ
ω
, de manera que la ecuación (2.67) se reduce a,
+B γ − r0 γ + r0
∫−B
RAγ dz = ∫
−B
RAγ ω
dz + ∫
− r0
RAγ γκ
dz . (2.68)
Por otro lado, el segundo sumando en la ecuación (2.66) puede descomponerse en,
0 γ − r0 γ 0 γ
∫−B
RAγ ω
dz = ∫
−B
RAγ ω
dz + ∫
− r0
RAγ ω
dz (2.69)
a) b)
45
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( )
0 γ r0 γ
. (2.70)
=∫ RA γκ
− RAγ ω
dz + ∫ RAγ ω
dz
− r0 0
1
RA s
=
Aηω ∫R As dA = RAs ηω
= kc Aγ ηω
, (2.71)
Aηω
46
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 3:
Planteamiento y solución del problema de cerradura
∂c Aγ
= ∇ i(Dγ ∇c Aγ ) , en la fase − γ , (3.1)
∂t
C.F. −nγκ iDγ ∇c Aγ = kcAγ , en la interfase γ − κ . (3.2)
Por otro lado, al usar el método del promedio volumétrico (Whitaker, 1999) fue posible
tiene restricciones de escala y es por tanto válida tanto en las regiones homogéneas
como en la inter-región,
( ) − ka ( x ) c
∂ c Aγ γ
εγ ( x) = ∇ i ε γ ( x ) D ( x )i∇ c Aγ Aγ , (3.3)
∂t
v γκ
A partir de la ecuación (3.3) fue posible obtener, al imponer las restricciones de escala
)−a
∂ c Aγ
ε γω
∂t
ω
(
= ∇ i ε γωDω i∇ c Aγ
γ
ω vω k c Aγ
γ
ω
, en la región homogénea − ω , (3.5)
γ
∂ cAγ
∂t
η
(
= ∇ i Dγ ∇ c Aγ
γ
η ), en la región homogénea − η . (3.6)
Por otro lado, siguiendo la metodología de Ochoa-Tapia y Whitaker (1995a) fue posible
47
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
(
−nηω i Dγ ∇ c Aγ
γ
η ) + n i(ε D i∇ c ) = kc
ηω γ ω Aγ
γ
ω Aγ ηω , en la inter-región ω − η . (3.7)
concentración,
γ γ
c Aγ = c Aγ , en la inter-región ω − η . (3.8)
η ω
Debe mencionarse que una alternativa a este planteamiento consiste en imponer una
condición de continuidad del flux de masa y desarrollar una condición de salto para la
Por otro lado, la condición de (3.7) no está cerrada, por lo que en este capítulo se
col., 2000),
γ
cAγ = c Aγ + c Aγ , (3.9)
k= k ηω
+k . (3.10)
obtiene,
γ γ
kcAγ ηω
= k ηω
c Aγ + c Aγ k ηω
+ k c Aγ + kcAγ . (3.11)
ηω
48
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Como se puede notar, lo que se necesita para cerrar esta última ecuación, es una
la inter-región.
( )
∂ cAγ γ γ
= Dγ ∇ 2 c Aγ + ε γ−1 ( x ) ∇ε γ ( x )i Dγ ∇ cAγ
∂t
⎡D ⎛ γ γ ⎞ ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − c Aγ ⎟ dA⎥ (3.13)
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
⎡D ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ n c
γκ Aγ dA ⎥ − ε γ−1 ( x ) kav ( x ) cAγ .
⎢⎣V Aγκ ( x ) ⎥⎦ γκ
⎡D ⎛ γ γ ⎞ ⎤
−ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ cAγ − c Aγ ⎟ dA⎥ , en la fase γ . (3.14)
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
⎡D ⎤
−ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ c Aγ dA⎥ + ε γ−1 ( x ) kav ( x ) cAγ
⎢⎣V Aγκ ( x ) ⎥⎦ γκ
Dγ ∇ 2 c Aγ = ε γ−1 ( x ) ∇ε γ ( x )i⎛⎜ Dγ ∇ c Aγ ⎞
γ
⎟
⎝ x ⎠
⎡D ⎛ γ γ ⎞ ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − c Aγ ⎟ dA⎥ (3.16)
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
⎡D ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ n γκ c γ dA ⎥ − ε γ−1 ( x ) kav ( x ) cAγ .
⎢⎣V
A γκ
Aγκ ( x ) ⎥⎦
49
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Dγ t *
= O (1) . (3.17)
LC LC1
γ ⎛ γ ⎞
−nγκ iDγ ∇ cAγ − nγκ iDγ ∇c Aγ = k ⎜ c Aγ + c Aγ ⎟ . (3.18)
x + yγ ⎝ x + yγ ⎠
Al hacer un análisis de órdenes de magnitud a esta ecuación, resulta (Whitaker, 1999),
⎡ ⎤
⎢ ⎥
⎢ klγ / Dγ lγ / LC ⎥ γ
cAγ =O⎢ , c Aγ . (3.19)
⎛ kl ⎞ ⎛ klγ ⎞ ⎥ x + yγ
⎢1 + O ⎜ γ ⎟⎟ 1 + O ⎜⎜ ⎟⎟ ⎥
⎢⎣ ⎜D
⎝ γ ⎠ ⎝ Dγ ⎠ ⎥⎦
γ
Aquí se supuso que la longitud característica asociada con los cambios de cAγ es LC .
klγ
1, (3.21)
Dγ
lγ
1. (3.22)
LC
La restricción sobre el término de reacción se satisface, en general, debido a la
kavω rp2
≤ 10 . (3.23)
ε γω Dω
50
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
−nγκ iDγ ∇ c Aγ − nγκ iDγ ∇c Aγ = k c Aγ , en la interfase γ -κ . (3.24)
x + yγ x + yγ
1 γ
A γκ ( x ) ∫( )
cAγ γκ
= cAγ dA . (3.25)
x + yγ
Aγκ x
Además, la restricción de escala dada por la ecuación (3.22) permite afirmar que
⎡ 1 ⎤
∇ 2c Aγ ε γ−1 ( x ) ∇ i ⎢
V Aγκ∫( x ) γκ Aγ ⎥
n c dA ⎥. (3.26)
difusión
⎢
⎣ ⎦
difusión no-local
⎞ − ε γ ( x ) kav ( x )
−1
= ε γ ( x ) ∇ε γ ( x )i⎛⎜ Dγ ∇ cAγ
γ γ
A γκ ( x ) Aγκ∫( x ) Aγ
2 −1
Dγ ∇ cAγ ⎟ c dA
⎝ x⎠ x + yγ
(3.27)
⎡D ⎛ γ γ ⎞ ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − cAγ ⎟ dA ⎥ .
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
región como la que se muestra en la Figura 3-1. Dicha celda unitaria debe ser
⎞ − ε γ ( x ) kav ( x )
−1
= ε γ ( x ) ∇ε γ ( x )i⎛⎜ Dγ ∇ cAγ
γ γ
A γκ ( x ) Aγκ∫( x ) Aγ
2 −1
Dγ ∇ cAγ ⎟ c dA
⎝ x⎠ x + yγ
fuente difusiva volumétrica
fuente reactiva non-local
(3.28)
⎡D ⎛ γ γ ⎞ ⎤
+ ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − c Aγ ⎟ dA⎥ , en la fase − γ ,
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
fuente difusiva volumétrica no-local
γ γ
−nγκ iDγ ∇ c Aγ − nγκ iDγ ∇c Aγ = k c Aγ , en la interfase γ -κ . (3.29)
x + yγ x + yγ
51
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
La condición (3.30) implica que la celda unitaria debe ser espacialmente periódica a lo
difusivas, de las cuales cabe distinguir entre las fuentes difusivas superficiales y
volumétricas, cuyos órdenes de magnitud pueden no ser los mismos. Si este es el caso,
⎧orden de magnitud ⎫
⎪ ⎪ 1 γ ⎛D γ ⎞
⎨ de la fuente ⎬= ∫ nγκ iDγ ∇ c Aγ dA = O ⎜ γ ∇ c Aγ
⎜l ⎟⎟ , (3.31)
⎪ difusiva superficial ⎪ V Aγκ ( x ) ⎝ γ ⎠
⎩ ⎭
⎧ orden de magnitud ⎫
⎪
⎨ de la fuente
⎪ 1
⎬=
⎪difusiva volumétrica ⎪ V
∫ (ε γ ( x ) )
Vγ ( x )
−1
(
∇ε γ ( x )i Dγ ∇ cAγ
γ
) dV
⎩ ⎭ . (3.32)
⎛ Δε D γ ⎞
= O ⎜ γ γ ∇ c Aγ ⎟⎟ .
⎜ ε r
⎝ γ 0 ⎠
52
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Las ecuaciones anteriores permiten hacer el estimado de la razón entre las fuentes
fuente volumétrica ⎛ Δε l ⎞
= O⎜ γ γ ⎟⎟ . (3.33)
fuente superficial ⎜ ε r
⎝ γ 0 ⎠
ε γ−1 ( x ) kav ( x ) γ
A γκ ( x ) Aγκ∫( x ) Aγ
Dγ ∇ 2 c Aγ = − c dA
x + yγ
γ γ
C.F.1 −nγκ iDγ ∇ c Aγ − nγκ iDγ ∇c Aγ = k c Aγ , en Aγκ ( x ) , (3.36)
x + yγ x + yγ
Periodicidad: c Aγ ( r + li ) = c Aγ ( r ) , i = 1, 2. (3.37)
Por otro lado, como se indica en la Figura 3-1, el alto de la celda unitaria debe ser tal
frontera,
ahora,
53
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
ε γ−1 ( x ) kav ( x ) γ
A γκ ( x ) Aγκ∫( x ) Aγ
Dγ ∇ 2 c Aγ = − c dA
x + yγ
, en la fase- γ , (3.40)
⎡D ⎛ γ γ ⎞ ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − c Aγ ⎟ dA ⎥
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
γ γ
C.F.1 −nγκ iDγ ∇ c Aγ − nγκ iDγ ∇cAγ = k cAγ , en Aγκ ( x ) , (3.41)
x + yγ x + yγ
Periodicidad: c Aγ ( r + li ) = c Aγ ( r ) , i = 1, 2. (3.44)
En este punto es conveniente notar que, mientras que la geometría en la que está
γ
r → r + li , el problema de valor a la frontera no lo será a menos que cAγ y su
x + yγ
gradiente puedan tratarse como constantes. Para ello se consideran las siguientes
γ γ γ 1 γ
cAγ = c Aγ + y γ i∇ c Aγ + y γ y γ : ∇∇ c Aγ + ... , (3.45)
x + yγ x x 2 x
γ γ γ
∇ c Aγ = ∇ c Aγ + y γ i∇∇ c Aγ + ... (3.46)
x + yγ x x
γ
Dγ ∇ 2 cAγ = −ε γ−1 ( x ) kav ( x ) c Aγ
x
ε γ−1 ( x ) kav ( x ) ⎡ ⎤ γ
− ⎢ ∫ y γ dA⎥ i∇ cAγ
A γκ ( x ) ⎢ Aγκ ( x ) ⎥⎦
⎣
x
ε γ−1 ( x ) kav ( x ) ⎡ 1 ⎤ γ
− ⎢ ∫ y γ y γ dA⎥ : ∇∇ cAγ , en la fase- γ , (3.47)
A γκ ( x ) ⎢ Aγκ ( x ) 2 ⎥
⎣ ⎦
x
γ
−ε γ−1 ( x ) ∇ i ⎡⎣Dγ ∇ y γ ⎤⎦ i∇ cAγ
x
⎡D ⎤ γ
−ε γ−1 ( x ) ∇ i ⎢ γ ∇ y γ y γ ⎥ : ∇∇ cAγ + ...
⎣ 2 ⎦ x
54
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ γ
−nγκ iDγ ∇cAγ = k cAγ + nγκ iDγ ∇ cAγ + nγκ iDγ y γ i∇∇ cAγ
x x x
(3.48)
γ k γ
+ k y γ i∇ c Aγ + y γ y γ : ∇∇ c Aγ + ..., en Aγκ ( x ) .
x 2 x
1
V ∫
Aγκ ( x )
nγκ y γ y γ dA = −∇ y γ y γ . (3.50)
γ ε γ−1 ( x ) kav ( x ) ⎡ ⎤ γ
−ε γ
−1
( x ) kav ( x ) cAγ − ⎢ ∫ y γ dA⎥ i∇ c Aγ
A γκ ( x ) ⎢ Aγκ ( x ) ⎥⎦
⎣
x x
ε γ−1 ( x ) kav ( x ) ⎡ 1 ⎤ γ γ
− ⎢ ∫ y γ y γ dA⎥ : ∇∇ c Aγ − ε γ−1 ( x ) ∇ i ⎡⎣Dγ ∇ y γ ⎤⎦ i∇ c Aγ
A γκ ( x ) ⎢ Aγκ ( x ) 2 ⎥
⎣ ⎦
x x
γ γ (3.51)
⎡D ⎤ γ
⎛ ka c ⎞ ⎛ ka r Δ c ⎞
−ε γ ( x ) ∇ i ⎢ γ ∇ y γ y γ ⎥ : ∇∇ c Aγ = O⎜ ⎟ + O⎜ v 0 ⎟
−1 v Aγ Aγ
⎣ 2 ⎦ x ⎜ εγ ⎟ ⎜ LC ε γ ⎟
⎝ ⎠ ⎝ ⎠
( ) ( )
γ γ
⎛ ka r 2 Δ c γ
⎞ ⎛D Δ ∇ y Δ c ⎞ ⎛D Δ ∇ y y Δ c ⎞
+O ⎜
v 0 Aγ
⎟ + O⎜ γ γ Aγ
⎟ + O⎜ γ γ γ Aγ
⎟.
⎜ ε γ LC LC1 ⎟ ⎜ ε γ r0 LC ⎟ ⎜ ε γ r0 LC LC1 ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
r0 r02
1, 1. (3.52)
LC LC LC1
γ ε γ−1 ( x ) kav ( x ) ⎡ ⎤ γ
ε γ−1 ( x ) kav ( x ) c Aγ ⎢ ∫ y γ dA⎥ i∇ c Aγ , (3.53)
A γκ ( x ) ⎢ Aγκ ( x ) ⎥⎦
⎣
x x
γ ε γ−1 ( x ) kav ( x ) ⎡ 1 ⎤ γ
ε γ−1 ( x ) kav ( x ) c Aγ ⎢ ∫ y γ y γ dA⎥ : ∇∇ c Aγ . (3.54)
A γκ ( x ) ⎢ Aγκ ( x ) 2 ⎥
⎣ ⎦
x x
55
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
Dγ ∇ 2 c Aγ = −ε γ−1 ( x ) kav ( x ) cAγ − ε γ−1 ( x ) ∇ i ⎡⎣Dγ ∇ y γ ⎤⎦ i∇ cAγ
x x
(3.55)
⎡D ⎤ γ
−ε γ ( x ) ∇ i ⎢ γ ∇ y γ y γ ⎥ : ∇∇ c Aγ
−1
, en la fase − γ .
⎣ 2 ⎦ x
Nótese que los términos de orden superior no fueron incluidos en la ecuación anterior.
Por otro lado, el llevar a cabo un análisis de orden de magnitud en cada término en el
γ γ γ γ
k c Aγ + nγκ iDγ ∇ cAγ + nγκ iDγ y γ i∇∇ c Aγ + k y γ i∇ c Aγ
x x x x
γ
⎛D Δ c ⎞
k
+ y γ y γ : ∇∇ c Aγ
2
γ
x
(
= O k cAγ
γ
+ O⎜)⎜
γ
LC
Aγ
⎟
⎟
(3.56)
⎝ ⎠
γ
⎛D r Δ c ⎞ ⎛ kr Δ c γ ⎞ ⎛ kr 2 Δ c γ
⎞
+O ⎜ ⎟ + O⎜ 0 ⎟ + O⎜ 0 ⎟.
γ 0 Aγ Aγ Aγ
⎜ LC LC1 ⎟ ⎜ LC ⎟ ⎜ LC LC1 ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
γ γ γ k γ
k c Aγ k y γ i∇ cAγ ; k c Aγ y γ y γ : ∇∇ cAγ ;
x x x 2 x
γ γ
nγκ iDγ ∇ cAγ nγκ iDγ y γ i∇∇ cAγ . (3.57)
x x
γ γ
−nγκ iDγ ∇c Aγ = k c Aγ + nγκ iDγ ∇ c Aγ , at Aγκ ( x ) . (3.58)
x x
difusiva en la ecuación (3.58) con los dos últimos términos en la ecuación (3.55). Estos
γ
γ
⎛D Δ c ⎞
1
dA = O ⎜ ⎟,
γ Aγ
V ∫ nγκ iDγ ∇ c Aγ
x ⎜ LC lγ ⎟
(3.59)
Aγκ
⎝ ⎠
56
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( )
γ
⎛D Δ ∇ y Δ c ⎞
1 γ
dV = O ⎜ ⎟,
γ γ Aγ
∫ ε γ ( x ) ∇ i ⎡⎣Dγ ∇ y γ ⎤⎦ i∇ cAγ
−1
(3.60)
V Vγ ( x )
x ⎜ ε γ r0 LC ⎟
⎝ ⎠
( )
γ
⎛D Δ ∇ y y Δ c ⎞
1 ⎡D ⎤ γ
⎜ γ γ γ Aγ
⎟ . (3.61)
∫ ε γ ( x ) ∇ i ⎢ γ ∇ y γ y γ ⎥ : ∇∇ c Aγ
−1
dV = O
V Vγ ( x ) ⎣ 2 ⎦ x ⎜ ε γ r0 LC LC1 ⎟
⎝ ⎠
Por otro lado, al usar los estimados desarrollados por Pérez-Córdova y Ochoa-Tapia
(1995),
(
O yγ ) = O ⎛⎜⎝ 323 r (1 − ε ) ⎞⎟⎠ ∼ O ⎛⎜⎝ 101 r (1 − ε ) ⎞⎟⎠ ,
0 γω 0 γω (3.62)
⎛ r02 ⎞
O yγ yγ ( ) = O ⎜ (1 + ε γω ) ⎟ , (3.63)
⎝ 10 ⎠
γ
⎛ (1 − ε ) D r Δ c γ
⎞
1
ε γ ( x ) ∇ i ⎡⎣Dγ ∇ y γ ⎤⎦ i∇ c Aγ dV = O ⎜ ⎟ , (3.64)
γω γ 0 Aγ
∫
−1
V Vγ ( x )
x ⎜ 10ε γ 2
Lε LC ⎟
⎝ ⎠
⎡D ⎤ γ
⎛ (1 + ε ) D r 2 Δ c γ
⎞
1
ε γ ( x ) ∇ i ⎢ γ ∇ y γ y γ ⎥ : ∇∇ c Aγ dV = O ⎜ ⎟.
γω γ 0 Aγ
∫
−1
(3.65)
V Vγ ( x ) ⎣ 2 ⎦ x ⎜ 10ε γ 2
Lε LC LC1 ⎟
⎝ ⎠
γ
Dγ ∇ 2 c Aγ = −ε γ−1 ( x ) kav ( x ) c Aγ , en la fase- γ , (3.67)
x
fuente reactiva volumétrica
γ γ
C.F.1 −nγκ iDγ ∇ c Aγ − nγκ iDγ ∇c Aγ = k cAγ , en Aγκ ( x ) , (3.68)
x x
fuente difusiva superficial fuente reactiva superficial
57
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ γ
cAγ = sγ c Aγ + bγ i∇ cAγ , (3.72)
x x
frontera,
Problema I
kav ( x )
Dγ ∇ 2 sγ = − , en la fase- γ , (3.73)
εγ ( x)
k
−nγκ i∇sγ = , en Aγκ ( x ) , (3.74)
Dγ
sγ = 0, en y = 2h , (3.75)
Periodicidad sγ ( r + li ) = sγ ( r ) , i = 1, 2. (3.76)
Problema II
bγ = 0, en y = 2h , (3.79)
Periodicidad bγ ( r + li ) = bγ ( r ) , i = 1, 2. (3.80)
Junto con,
γ γ
sγ cAγ + bγ i∇ c Aγ = c Aγω , en y = 0 . (3.81)
determinar si es necesario resolver los dos problemas arriba presentados o sólo uno de
ellos.
58
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ ⎛l γ ⎞
kbγ i∇ c Aγ = O ⎜ γ k c Aγ s ⎟
. (3.83)
ηω
⎝L ⎠
s
bγ = O ( lγ ) . Por otro lado, el estimado de orden de magnitud del primer término del
k ηω
c Aγ
γ
s
= O k c Aγ ( γ
s ). (3.84)
salto se reduce a,
(
−nηω i Dγ ∇ c Aγ
γ
η ) + n i(ε D i∇ c ) = K
ηω γ ω Aγ
γ
ω eff cAγ
γ
ω
, (3.85)
K eff = k ηω
+ ksγ . (3.86)
ηω
59
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
por Aris (1975) a partir de la ponderación del flux en las bocas de los poros y la
Por otro lado, el salto en el flux en la ecuación (3.85) puede despreciarse cuando se
(
nηω i ε γ Dω i∇ c Aγ
γ
ω ) K eff c Aγ
γ
ω
(3.88)
La expresión anterior puede ser resultado de suponer que la velocidad de reacción total
que se lleva a cabo dentro de los microporos es mucho mayor que la velocidad total de
1 γ 1 γ
V ∫ a ωk
Vγ
v c Aγ ω
dV
V ∫K
Aηω
eff c Aγ ω
dA . (3.89)
Este es el caso en algunas aplicaciones prácticas en medios porosos (Aris, 1975). Bajo
(
nηω i Dγ ∇ c Aγ
γ
η ) = n i(ε D i∇ c ) .
ηω γ ω Aγ
γ
ω
(3.90)
Por último, como puede notarse de la ecuación (3.86), para obtener el coeficiente de
Para llevar a cabo la solución del Problema I, se debe considerar la condición de (3.81),
región homogénea −ω . Para ello, conviene recordar que en la Sección 3.2 se obtuvo la
ecuación (3.67) a partir de las ecuaciones válidas en la fase −γ y en todo el dominio, por
necesarias para obtener esta ecuación son más fáciles de satisfacer en esta región. De
kavω γ
∇ 2 c Aγω = − cAγ , en la fase − γ. (3.91)
ε γωDγ ω
Por otro lado, la condición de frontera (3.68) es válida dondequiera que exista una
interfase sólido-fluido, por lo que se puede utilizar con seguridad en esta región,
γ k γ
−nγκ i∇ c Aγ − nγκ i∇c Aγω = c Aγ , en Aγκ ( xω ) . (3.92)
ω Dγ ω
periodicidad,
Este problema es idéntico al expresado en las ecuaciones (1.4-51) a (1.4-53) del texto de
Whitaker (Cap. 1, 1999). Para resolver este problema se propone la siguiente solución,
γ γ
cAγω = sγω cAγ + bγω i∇ c Aγ . (3.94)
ω ω
γ γ
debe recordarse que en esta frontera es razonable suponer que cAγ = c Aγ , de
ω
sγ = sγω , en y = 0 . (3.96)
Por otro lado, dado que esta frontera es artificial, es aceptable imponer la condición,
61
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂c Aγ ∂cAγω
= , en y = 0 . (3.97)
∂y ∂y
∂sγ ∂sγω
= , en y = 0 . (3.98)
∂y ∂y
Problema III,
kavω
∇ 2 sγω = − , en la fase – γ, (3.99)
ε γωDγ
k
−nγκ i∇sγω = , en Aγκ ,ω , (3.100)
Dγ
Los problemas I y III no están restringidos, hasta el momento, a algún tipo de geometría
mostrado en la Figura 3-1 sería una labor por demás laboriosa y probablemente
mostrado en la Figura 3-2, donde además se identifican las celdas unitarias a utilizar
para resolver los problemas de cerradura. En adelante se denotará a este modelo como
62
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂sγω
= 0, en y = 0 . (3.102)
∂y
∂sγ
= 0, en y = 0 . (3.103)
∂y
Figura 3-2: Modelo de medio poroso centrado en las caras y celdas unitarias para la
planteamiento del Problema I para la mitad de la celda unitaria de la Figura 3-2 (las
63
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂ 2 sγ ∂ 2 sγ f1 (Y0 ) K1 2
+ =− φ , en la fase fluida, (3.104)
∂X 2
∂Y 2
f 2 (Y0 )
∂sγ φ2 A
C.F.1 = K1 , en X = , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , (3.105)
∂X avω l 2
∂sγ
= 0 , en X = 0, Y2 ( k − 1) < Y ≤ Y1 ( k ) , k = 1,..., nc, Y2 ( nc ) ≤ Y ≤ 2 H ;
∂X
1
X= , 0 ≤ Y ≤ 2H , (3.106)
2
∂sγ φ2 A
= K1 , en Y = H , 0 ≤ X ≤ , (3.107)
∂Y avω l 2
∂sγ φ2 A
− = K1 , en Y = Y1 ( k ) , 0 ≤ X ≤ , (3.108)
∂Y avω l 2
1
C.F.2 sγ = 0 , en Y = 2 H , 0 ≤ X ≤ , (3.109)
2
∂sγ 1
C.F.3 = 0 , en Y = 0 , 0 ≤ X ≤ , (3.110)
∂Y 2
donde nc es el número de cuadrados a usar en la celda unitaria; A es la relación del largo
del lado de cada cuadrado de la celda unitaria entre el largo de la base de la misma y el
k ηω
avω l 2
φ2 = . (3.111)
Dγ
En la ecuación (3.104) el subíndice “0” indica que las propiedades están evaluadas en la
como el lugar donde la porosidad adquiere el valor del promedio de las porosidades de
3.5.- Resultados.
El problema de cerradura descrito por las ecuaciones (3.104) a (3.110) fue resuelto
utilizando el método de diferencias finitas con inversión línea por línea (ver detalles en
64
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
una tolerancia de 1x10−4 y once cuadrados en la celda unitaria, es decir, once celdas
De esta forma se obtuvieron los resultados que se muestran en las Figuras 3-3 y 3-4. El
comportamiento de los resultados es similar al de los obtenidos por Wood y col. (2000),
efectivo para la inter-región disminuye hasta llegar a cero; este efecto se ve retardado
Figura 3-3: Coeficiente efectivo de reacción vs. φ como función de la razón de aspecto
A; vista bidimiensional.
65
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
inter-región entre un sistema formado por microporos y fluido. Con estos resultados, es
posible evaluar los efectos que tiene la condición de salto sobre el flux interfacial
convección, el sistema estudiado hasta este punto, puede tomarse como análogo al de
mezclado. Esta analogía será válida únicamente cuando las restricciones de escala
asociadas con los microporos pueden aplicarse a los poros de las partículas catalíticas.
Figura 3-4: Coeficiente efectivo de reacción vs. φ como función de la razón de aspecto
A; vista tridimiensional.
proponen dos modelos de celda unitaria, mostrados en las figuras 3-5 y 3-7,
l / 2 , una de cada dos filas de cuadrados del modelo de la Figura 3-2; mientras que el
66
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
anterior, una de cada dos columnas de cuadrados del modelo de la Figura 3-2. De esta
vertical.
al caso mostrado en la Figura 3-2 para resolver estos dos problemas; los detalles de las
Las rutinas numéricas para resolverlos difieren de las del caso anterior, primordialmente
tomó de 1x10−5 sobre el coeficiente efectivo. Una comparación de los resultados que se
obtienen con los tres tipos de modelos se incluye en el Apéndice A-4 para cada valor de
variable de cerradura para los tres modelos de celda unitaria discutidos aquí. Las
continuación se presentan algunos comentarios sobre los resultados de las figuras 3-6 y
3-8
en el cuerpo vertical, la frontera se vuelve más irregular. Esto puede ser la razón
67
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Los resultados presentados en esta parte del trabajo confirman las observaciones
efectuadas para el caso del coeficiente efectivo en las regiones homogéneas (Kim y col.,
1987), es decir, que el efecto del tipo de microestructura no tiene un efecto significativo
variable de cerradura para los tres modelos de medio poroso utilizados en esta parte del
5000 puntos.
Figura 3-6: Comparación del coeficiente efectivo obtenido para un medio poroso
Figura 3-8: Comparación del coeficiente efectivo obtenido para un medio poroso
el coeficiente efectivo obtenido no difiere en esa misma proporción, esto se debe al paso
69
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Para cerrar este capítulo, se presentan en la Figura 3-10, los resultados que se obtienen
al tratar al medio microporoso como un medio bidisperso, es decir con dos distintos
diámetros de microporo (ver Figura 3-9). Para elaborar este modelo de medio poroso, se
horizontal entre los cuadrados de medio sólido. Este modelo es similar al propuesto
vertical entre los cuadrados sea distinto que el horizontal fue tomada del trabajo de
esta forma una diferencia significativa con los resultados que se obtienen de usar el
modelo de medio poroso consistente en un medio de frontera regular y, por ende, del
70
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Esta es una muestra más de la versatilidad que ofrece el concepto de celda unitaria para
1.0
0.8
0.6
K eff εγω=0.99 0.91 0.75 0.51 0.19
k ηω 0.4
0.2
0.0
0.1 1 10 40
Figura 3-10: Coeficiente efectivo vs. φ para un medio poroso centrado en el cuerpo
Para concluir estos dos capítulos, se presentan a continuación los pasos clave de la
metodología utilizada:
escala apropiadas.
71
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
región.
9. Se propone (si la linealidad del problema lo permite) una solución del problema
cerradura.
problema de cerradura.
el coeficiente de salto.
Como puede notarse, los pasos 1-5 fueron originalmente propuestos por Ochoa-Tapia y
72
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
col. (2000) para completar el análisis a partir del planteamiento y solución del problema
de cerradura correspondiente.
Los resultados de estos capítulos pueden extenderse a casos más complicados, como es
fluido (Escala II, Figura 2-1) esto se lleva a cabo en el Capítulo 7 incluyendo efectos de
comprender el efecto que tienen las suposiciones que permiten calcular los coeficientes
73
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 4:
generalizada (la cual no está cerrada) que constituye una pieza clave en el desarrollo de
los valores del flux de masa en la inter-región predichos con este modelo en la frontera
entre un medio poroso y un fluido. Para ello, se plantean y resuelven los problemas de
Este tipo de modelo ha sido llamado por algunos autores (Arquis y Caltagirone, 1984;
Beckermann y col., 1987, 1988; Gobin y col., 1998; Goyeau y col., 2003) como modelo
de un solo dominio (One domain approach). Mientras que el modelo que involucra el
A pesar que el uso de modelos que implican una sola ecuación para describir el
interfacial). Por ello, en este capítulo se predicen dichos cambios a partir de la solución
74
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
el transporte difusivo de masa de la especie A desde el seno del fluido (fase −γ ) hacia
un medio poroso
∂c Aγ
= ∇ i(Dγ ∇c Aγ ) , en la fase − γ , (4.1)
∂t
acuerdo a la ecuación (4.2). Al igual que en el caso anterior, se supondrá que el mismo
fluido se encuentra dentro (región −ω ) y fuera (región −η ) del medio poroso (Figura 4-
1). Bajo esta suposición, la ecuación (4.1) es válida dondequiera que esté presente la
fase −γ . Note que las ecuaciones (4.1) y (4.2) son idénticas a las ecuaciones (2.1) y
problema local. Sin embargo, dado que este análisis es análogo al desarrollado en el
γ
∂ c Aγ γ
εγ ( x) = ∇ i ⎡⎢ε γ ( x ) D ( x )i∇ c Aγ ⎤.
⎥⎦ (4.3)
∂t ⎣
75
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
(2.19) como
γ
D ( x )i∇ cAγ γ 1 ⎛ γ γ ⎞
V γ ( x ) ∫( ) γκ ⎜⎝
= ∇ cAγ + n c Aγ − c Aγ ⎟ dA
Dγ Aγκ x
x + yγ x⎠
(4.4)
1
V γ ( x ) Aγκ∫( x )
+ nγκ cAγ dA.
γ
cAγ x + yγ
= c Aγ + c Aγ x + yγ
. (4.5)
x + yγ
76
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
En este momento, conviene retomar la ecuación (2.26), donde se presentan los valores
Caltagirone, 1984; Beckermann y col., 1987, 1988; Gobin y col., 1998; Goyeau y col.,
2003). De acuerdo a la ecuación (4.4), este modelo no está cerrado, pues se encuentra en
col. (2003) propusieron usar relaciones tipo Heaviside, lineales, sinusoidales y basadas
en funciones error para describir estos cambios espaciales. Sin embargo, no hay garantía
de que dichos modelos heurísticos puedan ofrecer una descripción adecuada del
transporte cerca de la frontera entre un medio poroso y un fluido. Por ello, en este
para describir los cambios especiales de las propiedades geométricas (porosidad, área
77
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( )+ε ( )
∂ cAγ γ γ
= ∇ i Dγ ∇ c Aγ γ
−1
( x ) ∇ε γ ( x )i Dγ ∇ cAγ
∂t
⎡D ⎛ γ γ ⎞ ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ nγκ ⎜ c Aγ − c Aγ ⎟ dA ⎥ (4.7)
⎢⎣V Aγκ ( x ) ⎝ x + yγ x⎠ ⎥⎦
⎡D ⎤
+ε γ−1 ( x ) ∇ i ⎢ γ ∫ γκ Aγ ⎥⎥ .
n c d A
⎢⎣V Aγκ ( x ) ⎦
∂c Aγ
= ∇ i(Dγ ∇c Aγ )
∂t
− ε γ−1 ( x ) ∇ε γ ( x )i Dγ ∇ cAγ ( γ
) , en la fase − γ , (4.8)
fuente difusiva volumétrica
⎡D ⎤
− ε γ−1 ( x ) ∇ i ⎢ γ
⎢⎣V
∫
Aγκ ( x )
(
nγκ Δ cAγ
γ
+ cAγ dA⎥)
⎥⎦
fuente difusiva volumétrica no-local
γ γ γ
donde se usó Δ cAγ = c Aγ − c Aγ .
x + yγ x
γ
−nγκ iDγ ∇c Aγ = nγκ iDγ ∇ c Aγ , en la interfase γ − κ . (4.9)
x + yγ
Sin embargo, como se mostró en el Capítulo 3 (ver Sección 3.2), si se cumplen las
siguientes restricciones de escala en la inter-región
Δε γ lγ
1;
(1 − ε ) l r γω γ 0
1;
(1 + ε ) γω r02lγ
1, (4.10)
r0ε γ 10ε γ L2ε 10ε γ L2ε LC1
γ
Δ c Aγ se puede considerar despreciable y como consecuencia
γ γ
c Aγ − c Aγ cAγ , de manera que el último término en el lado derecho de la
x + yγ x
78
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
⎡D ⎤ ⎡D ⎤
ε γ−1 ( x ) ∇ i ⎢ γ
⎣ ⎢V
∫
Aγκ ( x )
(
nγκ Δ c Aγ
γ
)
+ c Aγ dA⎥ = ε γ−1 ( x ) ∇ i ⎢ γ
⎥⎦ ⎢⎣V
∫ γκ Aγ ⎥⎥ . (4.11)
Aγκ ( x )
n c d A
⎦
γ
D ( x )i∇ c Aγ γ 1
V γ ( x ) Aγκ∫( x )
= ∇ c Aγ + nγκ c Aγ dA . (4.12)
Dγ
Δε γ lγ
1, (4.13)
ε γ r0
momento conviene hacer notar que, de las restricciones de escala arriba impuestas, las
más difíciles de satisfacer son las de primer orden. Esto puede restringir la validez de
t *Dγ
1, (4.14)
lγ2
γ
C.F.1 −nγκ iDγ ∇c Aγ = nγκ iDγ ∇ c Aγ , en la interfase γ − κ , (4.16)
x
fuente difusiva superficial
∂c Aγ ∂c Aγω
C.F.3 cAγ = c Aγω ; = , en Y = 0 , (4.18)
∂Y ∂Y
Periodicidad: c Aγ ( r + li ) = c Aγ ( r ) , i = 1, 2. (4.19)
79
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
El cual está referido a una celda unitaria como la mostrada en la Figura 4-2. En las
tal que las desviaciones espaciales corresponden a las de las regiones homogéneas. Cabe
mencionar que en la parte inferior de la celda unitaria se supuso que tanto el campo de
Dada la linealidad del problema de valor a la frontera, definido por las ecuaciones (4.15)
como sigue
γ
cAγ = bγ i∇ cAγ . (4.20)
80
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
C.F.2 bγ = 0, en Y = 2h , (4.23)
γ ∂bγ γ ∂c Aγω
C.F.3 bγ i∇ c Aγ = c Aγω ; i∇ c Aγ = , en Y = 0 , (4.24)
Y =0 ∂Y Y =0 ∂Y
Periodicidad: bγ ( r + li ) = bγ ( r ) , i = 1, 2. (4.25)
Por otro lado, a partir del trabajo de Ryan y col. (1981), en la región homogénea−ω se
tiene que
γ
c Aγω = bγω i∇ c Aγ . (4.26)
ω
γ γ
De manera que, suponiendo que ∇ c Aγ = ∇ c Aγ ω
, la condición de frontera (4.24)
Y =0
∂bγ ∂bγ ω
C.F.3 bγ = bγω ; = , en Y = 0 . (4.27)
∂Y ∂Y
∂bγω X
Componente X = 0, en Y = 0 , (4.28)
∂Y
∂bγ X
C.F.3a = 0, en Y = 0 , (4.30)
∂Y
C.F.3b bγ ,Y = 0, en Y = 0 . (4.31)
Note que, debido a las ecuaciones (4.30) y (4.31), no es necesario resolver los
81
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
De esta forma, al sustituir la ecuación (4.20) en el segundo sumando del lado derecho de
como sigue
εγ ( x)D ( x) ⎡ 1 ⎤
= ε γ ( x ) ⎢I + ∫ nγκ bγ dA⎥
Dγ ⎢⎣ V γ ( x ) Aγκ ( x ) ⎥⎦
(4.32)
1
= εγ ( x) I + ∫ nγκ bγ dA.
V Aγκ ( x )
X y Y de la variable de cerradura bγ .
usando diferencias finitas con el método de inversión línea por línea. Al igual que en el
caso anterior, los cálculos se hicieron usando diferentes tamaños de malla y valores de
3-2. Sin embargo, dado que se desean calcular los cambios espaciales del coeficiente de
diez cuadrados se obtienen resultados aceptables (es decir, con una desviación menor al
en el lado derecho de la ecuación (4.32) adquiere valores cada vez menores respecto al
primer sumando.
disminuye a la vez que aumenta la diferencia entre los valores de las componentes X y Y
1.0
εγω = 0.9
0.8
0.8
0.6
0.7
εγ Deff / Dγ
0.6
0.4 0.5 0.4
0.3
0.2 0.2
0.1
0.0
-0.50 -0.25 0.00 0.25 0.50
ξ /2h
Para tener una idea más clara de la magnitud de las diferencias de las componentes del
tensor de difusión efectiva, en la Tabla 4-1 se presenta una comparación de los cambios
espaciales de ambas componentes para ε γω = 0.2 y 0.8. Como puede notarse, para
ε γω = 0.2 , la diferencia en los resultados es inferior a 0.2%, mientras que para ε γω = 0.8
el valor máximo de la diferencia es 1.63%. Dado que estas diferencias son menores al
5%, puede concluirse que el uso de las componentes normal o tangencial del tensor de
difusión da lugar a predicciones similares de, por ejemplo, los perfiles de concentración
(Sección 4.5).
84
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
ξ ε γω = 0.2 ε γω = 0.8
2h Componente Componente % Componente Componente %
tangencial normal Diferencia* tangencial normal Diferencia*
-0.5 0.10875 0.10875 0.00000 0.63324 0.63312 0.01200
-0.4 0.19789 0.19759 0.03000 0.67004 0.66591 0.41300
-0.3 0.28701 0.28643 0.05800 0.70683 0.69927 0.75600
-0.2 0.37614 0.37529 0.08500 0.74348 0.73332 1.01600
-0.1 0.46532 0.46416 0.11600 0.78017 0.76816 1.20100
0.0 0.55439 0.55306 0.13300 0.81684 0.80056 1.62800
0.1 0.64351 0.64200 0.15100 0.85351 0.84045 1.30600
0.2 0.73266 0.73100 0.16600 0.89019 0.87822 1.19700
0.3 0.82176 0.82013 0.16300 0.92686 0.91726 0.9600
0.4 0.91089 0.90956 0.13300 0.96353 0.95777 0.57600
0.5 1.00000 1.00000 0.00000 1.00000 1.00000 0.00000
* % Diferencia = Componente normal − Componente tangencial x100%
Por otro lado, los resultados de la Figura 4-4 se pueden aproximar mediante las
siguientes expresiones
ε γ Deff , X ⎛ξ ⎞ ξ
Componente tangencial, = AX + BX ⎜ ⎟ , ∀ ∈ [ −0.5, 0.5] , (4.34)
Dγ ⎝ 2h ⎠ 2h
ε γ Deff ,Y ⎛ξ ⎞ ξ
Componente normal, = AY + BY ⎜ ⎟ , ∀ ∈ [ −0.5, 0.5] . (4.35)
Dγ ⎝ 2h ⎠ 2h
que en todos los ajustes el coeficiente de correlación fue siempre superior a 0.99.
85
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Tabla 4-2: Dependencia de los coeficientes de las ecuaciones (4.34) y (4.35) con ε γω .
ε γω AX BX AY BY
0.1 0.52604 0.94791 0.52565 0.94762
0.2 0.55439 0.89125 0.55294 0.89038
0.3 0.58551 0.82900 0.58199 0.82800
0.4 0.61996 0.76014 0.61450 0.75998
0.5 0.65874 0.68200 0.65130 0.68041
0.6 0.70299 0.59422 0.69313 0.59184
0.7 0.75471 0.49087 0.74272 0.48865
0.8 0.81679 0.36679 0.80307 0.36530
0.9 0.89336 0.21350 0.88258 0.21294
desigualdad
γ 1 ⎛ γ γ ⎞
∇ c Aγ ω ∫ nγκ ⎜ cAγ
V γ ( xω ) Aγκ ( xω ) ⎝ xω + y γ
− c Aγ
xω
⎟ dA .
⎠
(4.37)
86
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ
Dω i∇ cAγ γ 1
V γ ( xω ) Aγκ ∫( xω ) γκ Aγ
ω
= ∇ cAγ + n c dA , (4.38)
Dγ ω
γ
∂ c Aγ γ
ε γω ω
= ∇ i ⎡⎢ε γωDω i∇ c Aγ ⎤ , en la región homogénea − ω .
⎥⎦ (4.39)
∂t ⎣ ω
Por otro lado, en la región ocupada por el fluido homogéneo (región homogénea −η ), la
única fase presente es la fase −γ , por lo que si se satisfacen las restricciones de escala
γ
∂ cAγ
∂t
η
(
= ∇ i Dγ ∇ c Aγ
γ
η ) , en la región homogénea −η . (4.41)
de manera que
V ∞ =V ω +V η . (4.42)
87
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
A ∞ = A ω +Aη . (4.43)
Al repetir esta operación para las ecuaciones de medio efectivo en las regiones
a la siguiente expresión
γ γ
∂ c Aγ ∂ c Aγ
∫ ε γω dV + ∫
ω η
dV
Vω
∂t Vη
∂t
= ∫ nω i⎡⎢⎣ε γωDω i∇
Aω
c Aγ
γ
ω ⎥⎦
A η
(
⎤ dA + n i D ∇ c
∫ η γ Aγ
γ
η ) dA (4.45)
γ γ
+ ∫ nωη i⎡⎣⎢ε γωDω i∇
*
cAγ ω
− Dγ ∇ cAγ η
⎤ dA.
⎦⎥
Aωη
*
En la ecuación anterior, Aωη representa la porción de la superficie divisoria entre las dos
regiones asociada con V∞ . Debe recalcarse que el uso de las ecuaciones (4.39) y (4.41)
88
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ
⎡ ∂ c Aγ
γ
∂ c Aγ
γ
⎤ ⎡ ∂ c Aγ
γ
∂ c Aγ ⎤
⎢ε γ ( x ) ⎥ dV + ⎢ε γ ( x ) ⎥ dV
∫⎢ ∫⎢
η
− ε γω ω
−
∂t ∂t ⎥ ∂t ∂t ⎥
Vω
⎣ ⎦ Vη
⎢⎣ ⎦⎥
γ γ
= ∫ nω i⎡⎣⎢ε γ ( x ) D ( x )i∇
Aω
c Aγ − ε γωDω i∇ c Aγ ω
⎤ dA
⎦⎥ (4.46)
γ γ
+ ∫ nη i ⎡⎢ε γ ( x ) D ( x )i∇ c Aγ − Dγ ∇ c Aγ ⎤ dA
A
⎣ η ⎦⎥
η
γ γ
− ∫ nωη i⎡⎣⎢ε γωDω i∇
*
c Aγ ω
− Dγ ∇ c Aγ η
⎤ dA.
⎦⎥
Aωη
γ
(4.47)
⎡ ∂ cAγ
γ
∂ cAγ ⎤
⎢
+ ∫ εγ ( x) − η ⎥
dV .
⎢ ∂t ∂t ⎥
Vη
⎢⎣ ⎥⎦
Exceso de transporte superficial
γ γ γ
∫ n iε γ D i∇
C
s s s s c Aγ s
dσ = ∫ nω i⎡⎣⎢ε γ ( x ) D ( x )i∇
Aω
c Aγ − ε γωDω i∇ c Aγ ω
⎤ dA
⎦⎥
γ γ
(4.48)
+ ∫ nη i ⎡⎢ε γ ( x ) D ( x )i∇ cAγ − Dγ ∇ cAγ ⎤ dA,
η⎥
A
⎣ ⎦
η
donde C representa la curva cerrada que rodea la superficie divisoria (Figura 4-5) y n s
89
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
γ
∂ c Aγ γ
εγ s s
− ∇ s iε γ sD s i∇ s c Aγ
∂t s
exceso de transporte por difusión interna
exceso de acumulación superficial (4.49)
γ γ
= − nωη i ⎡⎢ε γωDω i∇ c Aγ − Dγ ∇ c Aγ ⎤ , en la superficie divisoria.
⎥⎦
⎣ ω η
γ γ
nωη iε γωDω i∇ cAγ ω
= nωη iDγ ∇ cAγ η
, en la superficie divisoria . (4.50)
Esto puede resultar en una simplificación excesiva del problema dado que, como se
errores causados por usar las ecuaciones (4.39) y (4.41) fuera de las regiones
homogéneas. Por ello, resulta más conveniente usar la ecuación (4.49) en lugar de la
difusión interna pueden calcularse usando las definiciones (4.47) y (4.48) suponiendo
γ γ
∂ c Aγ ⎡
ξ =0 ∂ cAγ
γ
∂ c Aγ ω ⎤
εγ s s
=∫ ⎢ε (ξ ) − ε γω ⎥ dξ
∂t ξ =− r0 ⎢ γ ∂t ∂t ⎥
⎣ ⎦
γ
(4.51)
⎡ ∂ c Aγ
γ
∂ c Aγ η ⎤
ξ = r0
+∫ ⎢ ε (ξ ) − ⎥ dξ ,
ξ =0 ⎢ γ ∂t ∂t ⎥
⎣⎢ ⎦⎥
γ γ γ
∫ nω i⎡⎣⎢ε γ ( x ) D ( x )i∇ ⎤ dA
1
∇ s iε γ sDs i∇ s c Aγ = c Aγ − ε γωDω i∇ cAγ
s A ωη *
Aω
ω ⎦⎥
(4.52)
1 γ γ
+
A ωη* ∫ nη i⎡⎢⎣ε γ ( x ) D ( x )i∇ cAγ
A
− Dγ ∇ c Aγ ⎤ dA.
⎥
η⎦
η
Debe notarse que para calcular estos términos se requiere de las soluciones tanto de la
90
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
medio efectivo de las regiones homogéneas. Dado que este es un proceso transitorio,
las ecuaciones (4.51) y (4.52) en un instante previo al usado para discretizar las
Por otro lado, como se expuso anteriormente, para propósitos de este trabajo se
γ γ
cAγ η
= cAγ ω
, en la superficie divisoria . (4.53)
cual se discutirá más a fondo en un trabajo futuro. De esta forma, este modelo está
∂U ∂ ⎛ ∂U ⎞
ε γ (ζ ) = ⎜ D (ζ ) ⎟ , ∀ζ ∈ ( −1,1) , (4.54)
∂τ ∂ζ ⎝ ∂ζ ⎠
∂U
= 0, en ζ = − H , (4.55)
∂ζ
91
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
U = 1, en ζ = 1 , (4.56)
γ γ
c Aγ cAγ ξ tDγ ε γ (ζ ) Deff (ζ )
U= , U0 = 0
,ζ = ,τ= , D (ζ ) = , (4.58)
( L / 2)
γ γ
cAγ c Aγ L/2 2
Dγ
∞ ∞
γ
donde cAγ ∞
es el valor de la concentración lejos de la frontera en la fase fluida,
γ
la cual se supone constante y c Aγ es el valor inicial de la concentración. En la
0
ecuación (4.55), H es una distancia relacionada con el largo del medio poroso, la
cual se fijó como tres veces la distancia adimensional desde la frontera hasta el
valores de H.
Debe remarcarse que, para usar este modelo, se debe definir una región de
2h 20lκ 20 lκ r0 1
δ= = = = . (4.59)
L / 2 L / 2 1 − ε γω 1 − ε γω r0 L / 2 10 1 − ε γω
difusión efectiva.
⎧ε , ζ ≤0
ε γ (ζ ) = ⎨ γω , (4.60)
⎩1, ζ >0
92
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
⎧ ε γω Deff ω
⎪ ,ζ ≤ 0
D (ζ ) = ⎨ Dγ . (4.61)
⎪ 1, ζ >0
⎩
• Modelo II. Por último se considera el siguiente problema de valor a la frontera
∂U ω ∂ 2U ω
ε γω = Dω , ∀ζ ∈ ( −1, 0 ) , (4.62)
∂τ ∂ζ 2
∂Uη ∂ 2Uη
= , ∀ζ ∈ ( 0,1) , (4.63)
∂τ ∂ζ 2
∂U ω
= 0, en ζ = −1, (4.64)
∂ζ
ε γ s ∂U s 1 ∂ ⎛ ε γ s Ds ∂U s ⎞ ∂U ω ∂Uη
− ⎜⎜ ⎟⎟ = − Dω + , en ζ = 0 , (4.65)
L / 2 ∂τ L / 2 ∂ζ ⎝ D γ ∂ζ ⎠ ∂ζ ∂ζ
Uη = U ω , en ζ = 0 , (4.66)
Uη = 1, en ζ = 1 , (4.67)
U ω = Uη = U 0 , cuando τ = 0 . (4.68)
se muestra en la Figura 4-6. Como se puede notar, para valores bajos de la porosidad de
la región homogénea −ω ( ε γω < 0.2 ) se presentan las mayores diferencias en los perfiles
pueden atribuirse a una mayor resistencia que presenta el medio poroso, dado que para
un retardo para el cálculo de los términos del lado izquierdo de la ecuación (4.65). Por
otro lado, para τ > 0.5 las predicciones obtenidas con el Modelo-I* son
considerablemente distintas a las resultantes con los otros modelos. Cabe mencionar
que, para valores extremos de τ ( τ < 0.1 y τ 50 ) todos los modelos ofrecen los
93
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
el seno del fluido. Note que, en todos los resultados, se puede observar una inflexión en
los perfiles cerca de ζ = 0 , además, todos los modelos predicen perfiles semejantes
1.0 1.0
τ = 50 τ = 50
0.8
10
0.8 10
0.6 5 0.6
5
U U
0.4 0.4 1
1 0.5
0.2 0.5 0.2
a) 0.1 b) 0.1
0.01 0.01
0.0 0.0
-3 -2 -1 0 1 -3 -2 -1 0 1
ζ ζ
1.0
1.0 τ = 50
τ = 50
0.8 10
0.8 10
0.6
0.6 5
5 U
U 0.5
0.4
0.4
1 1
0.5 0.1
0.2 0.1 0.2
d) 0.01
c) 0.01
0.0
0.0 -3 -2 -1 0 1
-3 -2 -1 0 1
ζ
ζ
( i i i i i ) Modelo-I* y ( − − − − ) Modelo-II.
Por otro lado, en la Figura 4-7 se muestra la respuesta transitoria del flux de masa en
este caso, las predicciones del Modelo I* son inferiores a las obtenidas con los otros
modelos. Al igual que en la Figura 4-6, al aumentar la porosidad las diferencias entre
los modelos disminuyen. Además, para valores extremos del tiempo ( τ < 0.1 y τ 10 )
94
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
las predicciones de los modelos se asemejan entre ellas por las razones explicadas
arriba.
0.5 0.5
0.3 0.3
0.2 0.2
0.1 a) 0.1 b)
0.0 0.0
0.01 0.1 1 10 50 0.01 0.1 1 10 50
τ τ
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
c) 0.1 d)
0.1
0.0
0.0 0.01 0.1 1 10 50
0.01 0.1 1 10 50
τ τ
Por último, debido a la aproximación adoptada en la solución del Modelo II, se observa
embargo, las áreas bajo las curvas resultantes de los modelos I y II son similares, por lo
Para concluir este capítulo, es pertinente presentar los siguientes comentarios sobre los
resultados obtenidos:
95
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
ser también una fuente de discrepancia en los resultados. Esto puede deberse a
la solución de este tipo de modelos constituye una alternativa para el cálculo del
movimiento.
96
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Parte II:
97
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 5:
fluido ha sido un tema de intensa investigación desde hace al menos tres décadas. Al
igual que en el caso de transferencias de masa, actualmente existen dos corrientes para
un solo dominio (one domain approach) en el cual el dominio (medio poroso, fluido y
transición continua de las propiedades efectivas. Para ello se requieren usar ecuaciones
que sean por lo tanto válidas en todo el dominio (Arquis y Caltagirone, 1984;
Beckermann y col., 1987, 1988; Gobin y col., 1998). La otra línea de análisis es el
modelo de dos dominios (two domain approach) en donde se usan ecuaciones de medio
efectivo en las regiones homogéneas las cuales se acoplan por medio de sus
Sobre esta última, Beavers y Joseph (BJ) (1967) estudiaron teórica y experimentalmente
promedio volumétrico
∂ vβ ∂ vβ η β
ε βω
−1 ω
− = vβ , (5.1)
∂y ∂y K βω ω
98
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
coeficiente (Goyeau y col., 2003; Deng y Martínez, 2005; Min y Kim, 2005; Chandesris
cual no sólo debe tener un significado físico claro, sino que además, sus valores se
puedan obtener a partir del problema de cerradura correspondiente. Dado que una buena
Tapia y Whitaker (1995a), en las siguientes secciones se presentan algunos de los pasos
más importantes.
análisis más sencillo y contar con un problema idéntico al estudiado por BJ (1967) y
99
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Por otro lado, en este contexto se denotará a la fase fluida, dentro y fuera del medio
poroso, como la fase − β , mientras que la fase sólida será referida como la fase −σ . Las
∇ i v β = 0, en la fase − β (5.2)
0 = −∇ p β + ρ β g + μ β ∇ 2 v β , en la fase − β (5.3)
Figura 5-1: Flujo de un fluido sobre un medio poroso en un sistema tipo BJ.
Al igual que en los casos anteriores, el análisis se inicia definiendo una región de
respectivamente como,
1
Promedio superficial ψβ
x
=
V ∫
Vβ ( x )
ψβ x+y β
dV , (5.5)
β 1
V β ( x ) Vβ∫( x ) β
Promedio intrínseco ψβ = ψ x+yβ
dV . (5.6)
x
100
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∇ivβ = 0 . (5.9)
101
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∇i vβ = 0 . (5.10)
siguiente expresión
0 = − ∇ pβ + ε β ( x ) ρ β g + μ β ∇ 2 v β , (5.11)
dentro del volumen de promediado. El uso del teorema del promediado espacial en el
1
∇ pβ = ∇ pβ +
V ∫
Aβσ ( x )
n βσ pβ dA . (5.12)
Por otro lado, al usar la descomposición espacial (Gray, 1975) para la presión,
β
pβ = pβ + pβ . (5.13)
x+y β
β 1 ⎛ β β ⎞
∇ pβ = ε β ( x ) ∇ pβ + ∫ n βσ ⎜ pβ − pβ ⎟ dA
V Aβσ ( x ) ⎝ x+y β x⎠
(5.14)
1
+
V ∫
Aβσ ( x )
n βσ pβ dA,
β 1 β
pβ ∇ε β ( x ) = − ∫ n βσ pβ dA . (5.15)
x V Aβσ ( x )
x
Por otro lado, utilizando el teorema del promediado espacial en el último término de la
1
∇ 2 v β = ∇ i∇v β = ∇ i ∇v β +
V ∫
Aβσ ( x )
n βσ i∇v β dA . (5.16)
102
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1
∇2 v β = ∇2 v β +
V ∫
Aβσ ( x )
n βσ i∇v β dA , (5.17)
β
vβ = vβ + vβ , (5.18)
x+y β x+y β x+yβ
1 ⎛ β β ⎞
∇2 v β = ∇2 v β +
V Aβσ ( x )∫n βσ i ⎜ ∇ v β
⎝ x+yβ
− ∇ vβ ⎟ dA
x⎠
(5.19)
β 1
−∇ε β ( x )i∇ v β + ∫ n βσ i∇v β dA.
V Aβσ ( x )
β β
0 = −ε β ( x ) ∇ pβ + ε β ( x ) ρ β g + μ β ∇ 2 v β − μ β ∇ε β ( x )i∇ v β
⎡ ⎞ ⎤ dA
+ μ β ∇ ⎛⎜ Δ v β
1 β β
+ ∫ n i −IΔ pβ
V Aβσ ( x ) βσ ⎢⎣ x ⎝ x
⎟⎥
⎠⎦
(5.20)
n βσ i( −Ipβ + μ β ∇v β ) dA,
1
+
V ∫
Aβσ ( x )
β β β β β β
donde se usaron Δ pβ = pβ − pβ y Δ vβ = vβ − vβ
x x+y β x x x+yβ x
β
0 = −∇ p β + ρ β g + μ β ε β−1 ( x ) ∇ 2 v β
(5.21)
(
− μ β ε β−1 ( x ) ∇ε β ( x )i∇ ε β−1 ( x ) v β )−μ K β
−1
β ( x )i vβ ,
como
n βσ i( Ipβ − μ β ∇v β ) dA
1
μ β K −β1 ( x )i v β =
V β ( x ) Aβσ∫( x )
(5.22)
⎡ ⎞ ⎤ dA.
− μ β ∇ ⎛⎜ Δ v β
1 β β
+ ∫
V β ( x ) Aβσ ( x )
n βσ i ⎢IΔ pβ
⎣ x ⎝
⎟⎥
x ⎠⎦
103
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
promediado de la velocidad.
constituyen el modelo de un solo dominio (Goyeau y col., 2003). Para utilizar este
modelo Goyeau y col. (2003) supusieron que los cambios espaciales de las propiedades
Por otro lado, conviene hacer notar que el término μ β K −β1 ( x )i v β corresponde al
referencia), el cual es el componente clave del término de exceso global (ecuación 68,
reportada del modelo de un solo dominio, en este trabajo se prefiere denotar al vector
⎧0 en la región homogénea − η
K −β1 ( x ) = ⎨ −1 , (5.23)
⎩K βω en la región homogénea − ω
Región homogénea−η
En esta región (x → x ) ,
η como ya lo han mostrado previamente Ochoa-Tapia y
104
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
r0 Lv , Lp , (5.24)
los valores de las variables promediadas corresponden a los de las variables puntuales,
de manera que
β β
pβ = pβ = pβ ,
xη xη + y β
(5.26)
β β
vβ = vβ = vβ ,
xη xη + y β
β
∇i vβ η
= 0, en la región homogénea − η , (5.27)
β
0 = −∇ p β + ρ β g + μβ ∇ 2 v β , en la región homogénea − η . (5.28)
η η
Debe notarse que esta última es idéntica a la ecuación de Stokes (5.3), siempre que se
β
cumpla que v β = vβ = v β , de la misma forma que la ecuación (5.27) es idéntica a
η η
Región homogénea−ω
⎡ ⎛ ⎞⎤
n βσ i( −Ipβ + μ β ∇v β ) dA
β β
∫
Aβσ ,ω
∫
Aβσ ,ω
n βσ i ⎢ −I Δ pβ
⎣ xω
+ μβ ∇ ⎜ Δ v β
⎝ xω
⎟ ⎥ dA .
⎠⎦
(5.29)
n βσ i( −Ipβ + μ β ∇v β ) dA .
1
μ β K −βω1 i v β ω
=−
V β ,ω ∫
Aβσ ,ω
(5.30)
Por otro lado, en esta región la porosidad se supone constante de manera que se
(
O ∇2 v β ω ) (
O ∇ ( ε βω )i∇ ε βω
−1
(vβ ω )) . (5.31)
(
∇ i ε βω v β
β
ω ) = 0, en la región homogénea − ω , (5.32)
β
0 = −∇ p β ω
+ ρ β g + ε βω
−1
μβ ∇ 2 v β ω
(5.33)
− μ β K βω i v β
−1
ω
, en la región homogénea − ω ,
esta última puede reacomodarse para obtener la ley de Darcy con la primera corrección
de Brinkman,
vβ ω
=−
K βω
μβ (
i ∇ pβ
β
ω
− ρ β g − ε βω
−1
μβ ∇ 2 v β ω ) , en la región homogénea − ω . (5.34)
ecuaciones (5.27), (5.28), (5.32) y (5.33) no son válidas en la inter-región. Por lo que,
siguiente sección.
V ∞ =V ω +V η . (5.35)
106
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
A ∞ = A ω +Aη . (5.36)
Ecuación de continuidad
la divergencia, se obtienen
∫ ∇i vβ
V∞
dV = ∫ ni v β
A∞
dA = ∫ nω i v β
Aω
dA + ∫ nη i v β dA = 0 ,
Aη
(5.37)
∫ nη i v β
Aη
η
dA + ∫ nηω i v β
*
η
dA = 0 , (5.38)
Aηω
∫ nω i v β
Aω
ω
dA + ∫ nωη i v β
*
ω
dA = 0 , (5.39)
Aωη
*
donde Aωη representa la porción de la superficie divisoria entre las dos regiones
contenida en V∞ (ver Figura 5-3). Al sumar las ecuaciones (5.38) y (5.39) y restando el
resultado a la ecuación (5.37) se obtiene
∫ nω i(
Aω
vβ − vβ ω ) dA + ∫ n i( v
Aη
η β − vβ η ) dA = ∫ n i( v
*
Aωη
ωη β ω − vβ η ) dA . (5.40)
∫ n i(δ
C
s vβ s
) dσ = ∫ n i( v
Aω
ω β − vβ ω ) dA + ∫ n i( v
Aη
η β − vβ η ) dA , (5.41)
107
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
donde C representa una curva cerrada sobre la superficie divisoria (Figura 5-3);
∫ ∇ i(δ
*
Aωη
s vβ s
) dA = ∫ n i( v
Aω
ω β − vβ ω ) dA + ∫ n i( v
Aη
η β − vβ η ) dA , (5.42)
nωη i v β ( ω
− vβ η ) = ∇ i(δs vβ s
) en la superficie divisoria . (5.43)
Por otro lado (como se mostró en el Capítulo 2), cabe mencionar que si los cambios
∫ ∇ i(δ
*
Aωη
s vβ s
) dA = ∫ n i( v
( )Aω xω
ω β − vβ ω ) dA + ∫
( )
Aη xη
(
nη i v β − v β η ) dA , (5.44)
satisfacen las restricciones de escala, bajo las cuales las ecuaciones de transporte
(
nωη i v β ω
− vβ η ) (
∇s i δ vβ s
). (5.45)
108
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
(
nωη i v β ω
− vβ η ) = 0, en la superficie divisoria , (5.46)
vβ η
= vβ ω
= vβ ηω
, en la superficie divisoria. (5.47)
Ecuación de movimiento
( )
β
0 = −∇ p β + ρ β g + μ β ∇ i ε β−1 ( x ) ∇ v β
(5.49)
+ μ β ε β−3 ( x ) ( ∇ε β ( x )i∇ε β ( x ) ) v β − μ β K −β1 ( x )i v β .
obtienen, respectivamente
β β
0 = − ∫ n ω pβ dA − ∫ nη pβ dA + ∫ ρ β g dV + ∫ ρ β g dV
Aω Aη Vω Vη
(
+ ∫ μ β nω i ε β−1 ( x ) ∇ v β ) dA + ∫ μ n i(ε β η
−1
β (x) ∇ vβ ) dA (5.50)
Aω Aη
+ ∫ μ β ε β ( x ) ( ∇ε β ( x )i∇ε β ( x ) ) v β dV − ∫ μ β K −β1 ( x )i v β dV ,
−3
V∞ V∞
β β
0 = − ∫ nη pβ η
dA − ∫ nηω pβ η
dA + ∫ ρ β gdV
Aη * Vη
Aηω
(5.51)
+ ∫ μ β nη i∇ v β η
dA + ∫ μ β nηω i∇ v β η
dA,
Aη *
Aηω
109
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β β
0 = − ∫ n ω pβ ω
dA − ∫ nωη pβ ω
dA + ∫ ρ β gdV + ∫ ε βω μβ nω i∇
−1
vβ ω
dA
Aω * Vω Aω
Aωη
(5.52)
∫ ε βω μ β nωη i∇ v β dA − ∫ μ β K −βω1 i v β
−1
+ ω ω
dV .
* Vω
Aωη
∫ nωη i⎡⎢⎣ −I (
*
Aωη
pβ
β
ω
− pβ η) + μ (ε ∇ v − ∇ v )⎤⎥⎦ dA
β
β
−1
βω β ω β η
(p − p ) dA − ∫ n ( p − p ) dA
β β β β
= − ∫ nω β β ω η β β η
Aω Aη
Aω
(
+ ∫ μ β nω i ε β−1 ( x ) ∇ v β − ε βω
−1
∇ vβ ω ) dA
(5.53)
Aη
(
+ ∫ μ β nη i ε β ( x ) ∇ v β − ∇ v β
−1
η ) dA
− ∫ μ β K −β1 ( x )i v β dV + ∫ μ β K −βω1 i v β ω
dV
V∞ Vω
+ ∫ μ β ε β−3 ( x ) ( ∇ε β ( x )i∇ε β ( x ) ) v β dV .
V∞
La estructura del lado derecho de esta última ecuación sugiere la introducción de los
∫ n i(δ
C
s T s )dσ = ∫ nω i⎡⎢⎣ −I (
Aω
pβ
β
− pβ
β
ω ) + μ (ε ( x ) ∇ v − ε
β
−1
β β
−1
βω ∇ vβ ω )⎤⎥⎦ dA
(5.54)
+ ∫ nη i ⎡ −I pβ ( p ) + μ (ε ( x ) ∇ v ) ⎤ dA,
β β −1
− − ∇ vβ
Aη
⎣⎢ β η β β β η ⎦⎥
∫ nωη i Tβ ∫ μ β K β ( x )i v β dV − ∫ μ β K −βω1 i v β
−1
s
dA = ω
dV , (5.55)
* V∞ Vω
Aωη
∫ nωη i Tβ ∫ μ β K β ( x )i v β dV − ∫ μ β K −βω1 i v β
−1
s
dA = ω
dV . (5.56)
* V∞ Vω
Aωη
condición de salto,
110
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
−nωη (p β ω
β
− pβ
β
η ) + μ n i(ε
β ωη
−1
βω ∇ vβ ω
− ∇ vβ η ) = ∇ iδ
s T s
μβ μβ (5.57)
+ K −β1
ηω
ivβ + ε β−3 ( ∇ε β i∇ε β ) vβ , en la superficie divisoria.
avs avs ηω
ηω ηω
μβ
nωη i Bβ = ε β−3 ( x ) ( ∇ε β ( x )i∇ε β ( x ) ) vβ . (5.59)
s s avs ηω ηω
Por otro lado, al igual que en la definición de la velocidad superficial, la ecuación (5.54)
∫ ∇ i(δ
*
Aωη
s T s ) dA = ∫
(Aω xω )
nω i ⎡ −I pβ
⎣⎢ ( β
− pβ
β
ω) + μ (ε ( x ) ∇ v − ε ∇ v )⎤⎦⎥ dA
β
−1
β β
−1
βω β ω
(5.60)
nη i ⎡ −I pβ ( p ) + μ ( ε ( x ) ∇ v − ∇ v ) ⎤ dA.
β β
∫
−1
+ −
( )
Aη xη
⎣⎢ β η β β
⎦⎥ β β η
consecuencia
μ β nωη i ε βω
−1
(∇ vβ ω
− ∇ vβ η ) ∇ s iδ T s . (5.61)
111
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
−nωη (p β
β ω − pβ
β
η ) + μ n i(ε
β ωη
−1
βω ∇ vβ ω
− ∇ vβ η ) = μavs
β
K −1 i v β ω (5.62)
en la superficie divisoria,
para obtener la ecuación anterior se despreciaron los cambios espaciales de la
microestructura del medio poroso (de aquí en adelante se les denotará como esfuerzos
*
como constante en el pequeño dominio de integración Aηω (ver Sección 5.5). La cual es
divisoria. No debe olvidarse que se está analizando el sistema estudiado por Beavers y
vβ λ
= vβ λ
ii , λ = ω ,η . (5.64)
112
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂ vβ ∂ vβ K −1
η
ε βω
−1 ω
− = vβ , en la superficie divisoria , (5.65)
∂y ∂y avs ω
K β ,ω K −1
β= . (5.67)
avs
puede obtenerse directamente a partir de una expresión que describa los cambios
llevarse a cabo a partir de la ecuación (5.22). Esto implica el conocer los campos de las
objetivo del siguiente capítulo. Por el momento se presentan los detalles que permiten
Relación entre t β s
y tβs .
113
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∫
*
tβ s
dA = ∫
V∞
t β dV − ∫ t β
Vω
ω
dV . (5.69)
Aωη
μβ y =+ r0
tβs = K −β1 i v β = ∫ t β dy . (5.70)
avs s y =− r0
donde avs se refiere al valor del área interfacial por unidad de volumen ( av = A βσ /V )
evaluada en la superficie.
Para cumplir el objetivo es conveniente, para el caso en que los efectos de curvatura de
y =+ B y =0
tβ =∫ t β dy − ∫ tβ ω
dy . (5.71)
s y =− B y =− B
(las longitudes características asociadas con la inter-región son las mismas que en los
casos anteriores).
y =0 y =+ r0
∫ y =− r0
tβ ω
dy y representa ∫ y =− r0
t β dy .
114
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
tβ
ω
, de manera que la ecuación (5.71) se puede simplificar a
y =+ r0 y =0
tβ =∫ t β dy − ∫ tβ ω
dy . (5.72)
s y =− r0 y =− r0
1 μβ
− tβ s
=
Aηω * ∫ tβ
*
s dA = t β s ηω
=
avs
K −β1
ηω
ivβ
ηω
. (5.73)
Aηω
(5.72) es negativa (ver Figura 5-4). Cuando este no sea el caso, se debe cambiar el signo
Relación entre b β s
y bβ s .
∫
Aωη
bβ s
dA = ∫
V∞
b β dV , (5.74)
bβ s
= nωη i Bβ s
; b β = μ β ε β−3 ( x ) ( ∇ε β ( x )i∇ε β ( x ) ) v β . (5.75)
Es importante mencionar que estos esfuerzos, a diferencia de los globales, sólo toman
115
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Debe recordarse que el término b β viene de una combinación de las dos correcciones
interpretar como esfuerzos generados por los cambios espaciales de la porosidad cerca
Al igual que en el caso anterior, el primer paso consiste en obtener la forma simplificada
despreciables,
y =+ B
bβ =∫ bβ dy . (5.76)
s y =− B
se puede simplificar a,
y =+ r0
bβ =∫ b β dy . (5.77)
s y =− r0
μβ
ε β−3 ( x ) ( ∇ε β ( x )i∇ε β ( x ) )
1
bβ s
=
Aηω * ∫*
b β s dA = b β s ηω
=
avs ηω
vβ
ηω
. (5.78)
Aηω
116
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 6:
región. Estos últimos pueden estimarse a partir de las expresiones desarrolladas por
una expresión que permita predecir los esfuerzos globales en la frontera. Para resolver
este problema, conviene recordar que en el capítulo anterior se obtuvo una definición de
los esfuerzos globales dependiente de la posición K −β1 ( x ) que puede usarse para
compararán los resultados con los reportados por Goyeau y col. (2003).
117
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
fluido ( fase − β ) a través de un medio poroso (compuesto por las fases fluida y sólida
(fase −σ )),
∇ i v β = 0, en la fase − β , (6.1)
0 = −∇ p β + ρ β g + μ β ∇ 2 v β , en la fase − β , (6.2)
v β = 0, en la interfase β − σ . (6.3)
(
∇i ε β ( x) vβ
β
)=0, (6.4)
β
0 = −∇ p β + ρ β g + μ β ε β−1 ( x ) ∇ 2 v β
. (6.5)
(
− μ β ε β−1 ( x ) ∇ε β ( x )i∇ ε β−1 ( x ) v β )−μ K β
−1
β ( x )i vβ .
n βσ i( Ipβ − μ β ∇v β ) dA
1
μ β K −β1 ( x )i v β =
V β ( x ) Aβσ∫( x )
(6.6)
1 ⎡⎛ β β ⎞ ⎛ β β ⎞⎤
+ ∫
V β ( x ) Aβσ ( x )
n βσ i ⎢I ⎜ pβ
⎣⎝ x+yβ
− pβ ⎟ + μβ ⎜ ∇ v β
x⎠ ⎝ x
− ∇ vβ
x+ yβ
⎟ ⎥ dA,
⎠⎦
⎧ 0, en la región homogénea − η
. K −β1 ( x ) = ⎨ −1 . (6.7)
⎩ K βω , en la región homogénea − ω
118
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
∇i vβ = 0, en la región homogénea − η , (6.8)
η
β
0 = −∇ p β + ρ β g + μβ ∇ 2 v β , en la región homogénea − η , (6.9)
η η
(
∇ i ε βω v β
β
ω ) = 0, en la región homogénea − ω , (6.10)
β
0 = −∇ p β ω
+ ρ β g + ε βω
−1
μβ ∇ 2 v β ω
(6.11)
− μ β K βω i v β
−1
ω
, en la región homogénea − ω.
Por otro lado, al usar la metodología de Ochoa-Tapia y Whitaker (1995a) junto con la
(
nωη i v β ω
− vβ η )=0 en la superficie divisoria , (6.12)
−nωη (p β
β ω − pβ
β
η ) + μ n i(εβ ωη
−1
βω ∇ vβ ω
− ∇ vβ η ) = aμ β
vs
K −1 i v β ω (6.13)
en la superficie divisoria.
119
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
pβ = pβ − pβ , (6.15)
x+y β x+ y β x+yβ
β
vβ x+y β
= vβ x+ y β
− vβ . (6.16)
x+y β
El uso de las ecuaciones anteriores en la resta de las ecuaciones (6.1) y (6.4) así como
β
0 = −∇ p β + μ β ∇ 2 v β − μ β v β ε β−1 ( x ) ∇ 2ε β ( x )
fuente
β β
(6.18)
+ μ β ε β−1 ( x ) ∇ε β ( x )i∇ v β + μ β ε β ( x ) K −βσ1 ( x )i v β , en la fase − β .
fuente fuente
β
C.F.1 vβ = − vβ , en la interfase β − σ . (6.19)
fuente
Debe ser claro que no se desea resolver el problema de valor a la frontera para las
representativa como la que se muestra en la Figura 6-1. Donde se puede notar que a una
determinada distancia 2h, las desviaciones espaciales corresponden a las de las regiones
C.F.2 En ζ = − h v β = vω ; pβ = pβ ,ω , (6.20)
C.F.3 En ζ = h v β = vη = 0; pβ = pβ ,η = 0 . (6.21)
En esta última ecuación se aprovechó el hecho que las ecuaciones de medio efectivo en
Periodicidad: pβ ( r + li ) = pβ ( r ) , v β ( r + li ) = v β ( r ) , i = 1, 2 . (6.22)
En este momento es conveniente notar que el problema de valor a la frontera para las
β β β 1 β
vβ = vβ + y β i∇ v β + y β y β : ∇∇ v β + ... , (6.23)
x+yβ x x 2 x
121
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β β β 1 β
∇ vβ = ∇ vβ + y β i∇∇ v β + y β y β : ∇∇∇ v β + ... (6.24)
x+yβ x x 2 x
β β
∇ i v β = ε β−1 ( x ) ∇ε β ( x )i v β + ε β−1 ( x ) ∇ε β ( x )iy β i∇ v β (6.25)
x x
β
0 = −∇pβ + μ β ∇ 2 v β − μ β ε β−1 ( x ) ∇ 2ε β ( x ) v β
x
β
− μ β ε β−1 ( x ) ∇ 2ε β ( x ) y β i∇ v β
x
β
+ μ β ε β−1 ( x ) ∇ε β ( x )i∇ v β
x
(6.26)
β
+ μ β ε β−1 ( x ) ∇ε β ( x )iy β i∇∇ v β
x
β
+ μ β ε β ( x ) K −β1 ( x )i v β
x
β
+ μ β ε β ( x ) K β ( x )iy β i∇ v β
−1
,
x
β β
C.F.1 vβ = − vβ − y β i∇ v β . (6.27)
x x
de escala,
r0 r0
1, 1, (6.28)
Lv Lv1
β
donde Lv es la longitud característica asociada con v β , mientras que Lv1 es la
x
β
longitud característica correspondiente a ∇ v β . Bajo estas condiciones, el problema
x
β
∇ i v β = ε β−1 ( x ) ∇ε β ( x )i v β , en la fase − β , (6.29)
x
fuente
β
0 = −∇ p β + μ β ∇ 2 v β − μ β v β ε β−1 ( x ) ∇ 2ε β ( x )
x
fuente
(6.30)
β β
+ μ β ε β ( x ) ∇ε β ( x )i∇ v β
−1
+ μ β ε β ( x ) K β ( x )i v β−1
, en la fase − β ,
x x
fuente fuente
122
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
C.F.1 vβ = − vβ , en la interfase β − σ , (6.31)
x
fuente
C.F.2 En ζ = − h v β = vω , (6.32)
C.F.3 En ζ = h v β = vη = 0 , (6.33)
Periodicidad: pβ ( r + li ) = pβ ( r ) , v β ( r + li ) = v β ( r ) , i = 1, 2 . (6.34)
β β
pβ = μ β A β : ∇ v β + μβ b β i v β , (6.36)
donde los coeficientes que multiplican a las fuentes son las variables de cerradura.
Problema I:
C.F.2 En ζ = h Bβ = 0 , (6.40)
C.F.3 En ζ = −h B β = Bω , (6.41)
Periodicidad: b β ( r + li ) = b β ( r ) , Bβ ( r + li ) = Bβ ( r ) , i = 1, 2 . (6.42)
Problema II:
∇ iC β = 0, en la fase − β , (6.43)
123
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
C.F.2 En ζ = h Cβ = 0 , (6.46)
C.F.3 En ζ = −h Cβ = 0 , (6.47)
Periodicidad: A β ( r + li ) = A β ( r ) , Cβ ( r + li ) = Cβ ( r ) , i = 1, 2 . (6.48)
β β
En las ecuaciones (6.41) y (6.47) se supuso que en ζ = − h , vβ − vβ 1 . Además,
ω
β
vω = Bω i v β , (6.49)
ω
β
pβ ,ω = μ β bω i v β ω
, (6.50)
Por otro lado, en su trabajo Eidsath y col. (1983) propusieron imponer condiciones de
simetría en las posiciones de la celda unitaria (representativa del seno del medio poroso)
∂Bω
En ζ = − h = 0. (6.51)
∂ζ
De esta forma, al suponer continuidad de las derivadas del campo de las desviaciones
por,
∂Bβ
En ζ = − h = 0. (6.52)
∂ζ
homogénea −ω .
124
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Bβ = O (1) . (6.53)
⎡ lβ2 Δε β ⎤
O ( v β ) = O ⎢1,
β
⎥ vβ , (6.56)
⎢⎣ r0 Lε β ⎥⎦
β
⎡ l 2 Δε ⎤ μβ v β
O ( pβ ) = O ⎢1, β β ⎥ . (6.57)
⎢⎣ r0 Lε β ⎥⎦ lβ
sólo es necesario resolver el Problema I para tener una estimación de los campos de
Cuando este no sea el caso, la definición del esfuerzo global y, por tanto, la estructura
125
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Hasta este punto se ha conservado cierto grado de generalidad al no fijar la posición del
resolver este problema en cada posición de la inter-región y, de esta forma, calcular los
de un solo dominio (Goyeau y col., 2003). Sin embargo, para obtener una forma cerrada
superficie divisoria entre el medio poroso y el fluido (ver Sección 5.5). En este
lugar donde ε β =
1
2
(1 + ε βω ) y será referida, de aquí en adelante, como la posición
donde el centroide está ubicado en x = x 0 .
No debe olvidarse que el propósito de estimar los campos de las variables de cerradura
se obtiene
μβ ⎡ ⎤
⎢ ∫ n βσ i( Ib β − ∇Bβ ) dA⎥ i v β
β
μ β K −β1 ( x0 )i v β =
x0 V β ( x0 ) ⎢ Aβσ ( x0 ) ⎥⎦
⎣
x0
(6.61)
1 ⎡ β ⎛ β ⎞⎤
+ ∫
V β ( x0 ) Aβσ ( x0 )
n βσ i ⎢IΔ pβ
⎣ x0
− μβ ∇ ⎜ Δ v β
⎝
⎟ ⎥ dA.
x0 ⎠ ⎦
Con el fin de simplificar la ecuación anterior, se sustituyen las series de Taylor para
β β
pβ y ∇ vβ alrededor de x0 en el último término de la ecuación (6.61) para
obtener,
126
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
μβ ⎡ ⎤
⎢ ∫ n βσ i( −Ib β + ∇Bβ ) dA⎥ i v β
β
μ β K −β1 ( x0 )i v β =−
x0 V β ( x0 ) ⎢ Aβσ ( x0 ) ⎥⎦
⎣
x0
ε β−1 ( x0 ) ⎡⎛ β ⎞ ⎛ β ⎞⎤
+
V ∫
Aβσ ( x )
n βσ i ⎢I ⎜ y β i∇ pβ
⎣⎝
⎟ − μ β ⎜ y β i∇∇ v β
x0 ⎠ ⎝
⎟ ⎥ dA
x0 ⎠ ⎦
(6.62)
0
ε β ( x0 )
−1
⎡⎛ β ⎞ ⎛ β ⎞⎤
+
2V ∫
Aβσ ( x )
n βσ i ⎢I ⎜ y β y β : ∇∇ pβ
⎣⎝
⎟
x0 ⎠
− μ β ⎜ y β y β : ∇∇∇ v β
⎝
⎟ ⎥ dA + ...
x0 ⎠ ⎦
0
1
V ∫
Aβσ ( x )
n βσ y β dA = −∇ y β , (6.63)
1
V ∫
Aβσ ( x )
n βσ y β y β dA = −∇ y β y β . (6.64)
resulta que
μβ ⎡ ⎤
⎢ ∫ n βσ i( −Ib β + ∇B β ) dA⎥ i v β
β
μ β K −β1 ( x 0 )i v β =−
x0 V β ( x 0 ) ⎢ Aβσ ( x0 ) ⎥⎦
⎣
x0
⎡ β β ⎤
−ε β−1 ( x 0 ) ∇ y β : ⎢I∇ pβ − μ β ∇∇ v β ⎥ (6.65)
⎣ x0 x⎦
ε −1 ( x ) ⎡ β β ⎤
− β 0 ∇ y β y β ⎢I∇∇ pβ − μ β ∇∇∇ v β ⎥ + ...
2 ⎣ x0 x⎦
⎛ 3 r0 ⎞ ⎛ 1 r0 ⎞
ji∇ y β i j
ηω
= O⎜ ( ε βη − ε βω ) ⎟ ∼ O ⎜ (1 − ε βω ) ⎟ , (6.66)
⎝ 32 Lε ⎠ ⎝ 10 Lε ⎠
⎛ ⎞
(
ji∇ ji y β y β i j ) r2
= O ⎜ (1 + ε βω ) 0 ⎟. (6.67)
ηω
⎝ 10 Lε ⎠
De manera que los estimados de orden de magnitud de los tres términos del lado
μβ ⎡ ⎤ ⎡1 ⎤
⎢ ∫ n βσ i( −Ib β + ∇Bβ ) dA⎥ i v β
β β
= O ⎢ 2 ε β−1 ( x0 ) μ β v β ⎥ , (6.68)
V β ( x0 ) ⎢ Aβσ ( x0 ) ⎥⎦ ⎣⎢ lβ ⎦⎥
⎣
x0
127
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
⎡ β β ⎤
ε β−1 ( x0 ) ∇ y β : ⎢I∇ pβ − μ β ∇∇ v β ⎥
⎣ x0 x⎦
⎡ (1 − ε βω ) r −1
(6.69)
β ⎤
=O⎢ 0
ε β ( x0 ) μβ v β ⎥,
⎢⎣ 10 Lv Lv1 Lε ⎥⎦
ε β−1 ( x 0 ) ⎡ β β ⎤
∇ yβ yβ ⎢I∇∇ pβ − μ β ∇∇∇ v β ⎥
2 ⎣ x0 x⎦
(6.70)
⎡⎛ r 2 ⎞ (1 + ε βω ) −1 β ⎤
= O ⎢⎜ 0 ⎟ ε β ( x0 ) μβ v β ⎥.
⎢⎣⎝ Lε L ⎠ 20 Lv Lv1 ⎥⎦
(1 − ε ) l
βω β
2
r0 ⎛ r 2 ⎞ (1 + ε βω ) lβ
1; ⎜ 0 ⎟
2
1, (6.71)
10 Lv Lv1 Lε ⎝ Lε L ⎠ 20 Lv Lv1
las cuales son más fáciles de cumplir que la ecuación (6.58). Entonces, la ecuación
−1
ε β ( x0 ) K −β1 ( x0 ) = n i( −Ib β + ∇Bβ ) dA, en la superficie divisoria .(6.72)
V β ( x0 ) Aβσ∫( x0 ) βσ ηω
Para completar esta parte del análisis, en la siguiente sección se presenta el problema
bidimensional de valor ala frontera a resolver para calcular la forma escalar de este
coeficiente.
velocidad que es tangente a la superficie del medio poroso; se necesita la forma escalar
del tensor de esfuerzos globales, la cual se obtiene a partir de la ecuación (6.72) como
⎛ ∂Bβ ,ξξ ⎞
⎜⎜ − ( i in βσ ) bβ ,ξ +
1
ε β ( x0 ) K β−1 ( x0 ) = − ∫ ⎟⎟ dA , (6.73)
V β ( x0 ) Aβσ ( x0 ) ⎝ ∂n βσ ⎠ ηω
donde
128
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
bβ ,ξ = b β ii , (6.76)
condición de salto es
2
⎛ dε ⎞
K −1
= Kβ −1
+ εβ ⎜ β
−3
⎟⎟ . (6.77)
ηω ⎜ dy
⎝ ⎠ ηω
K β ,ω K −1
β= . (6.78)
avs
de resolver la forma escalar del Problema I para una celda unitaria bidimensional
dy 2 y =0
129
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂ vβ ∂ vβ ⎡ ⎛ ∂Bβ ,ξξ ⎞
∂y
ω
− ε βω
∂y
η
= ⎢−
⎢ V
1
∫ ⎜ ⎜ − ( i in βσ ) bβ ,ξ +
∂n βσ
⎟⎟ dA
Aβσ ( x0 ) ⎝ ⎠ x0
⎣
(6.87)
2 ⎤
⎛ dε ⎞ ⎥ ε βω
+ ε β−2 ⎜ β ⎟⎟ v , en la superficie divisoria.
⎜ dy ⎥ ε β ( x 0 ) avs β ω
⎝ ⎠ ηω ⎥
⎦
Para resolver el problema de valor a la frontera arriba planteado se pueden usar modelos
simples de medios porosos como los mostrados en el Capítulo 3 (ver Figuras 3-2, 3-5,
3-7 y 3-9). El problema de las variables de cerradura tiene una estructura similar al
problema de Stokes. Sin embargo, los últimos términos en las ecuaciones (6.80) y (6.81)
se puede superar, siguiendo las ideas de Barrère y col. (1992), mediante una
redefinición de variables. Una metodología alternativa para tratar con este problema
numérica del problema de cerradura está más allá de los objetivos de este trabajo y será
centroide se localiza en dicha región). Bajo esta idea es que, en la siguientes secciones
130
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
globales.
⎛ y ⎞ ⎡⎛ y ⎞ ⎤
2
εβ ( y) =
1
( ε βω + 1) + ( ε βω − 1) ⎜ ⎟ ⎢⎜ ⎟ − 3⎥ , para − r0 ≤ y ≤ + r0 ,
1
(6.88)
2 4 ⎝ r0 ⎠ ⎢⎣⎝ r0 ⎠ ⎥⎦
av ( y ) 1 1 ⎛ y ⎞ ⎡⎛ y ⎞ ⎤
2
= + ⎜ ⎟ ⎢⎜ ⎟ − 3⎥ , para − r0 ≤ y ≤ + r0 , (6.89)
av ,ω 2 4 ⎝ r0 ⎠ ⎢⎝ r0 ⎠ ⎥⎦
⎣
los cuales se obtuvieron suponiendo que el medio poroso está compuesto de esferas o
r0 = 20 lσ . (6.90)
dε β ( y ) ⎡⎛ y ⎞ 2 ⎤
=
3
( βω ) ⎢⎜ r ⎟ − 1⎥ , para − r0 ≤ y ≤ + r0 .
ε − 1 (6.91)
dy 4r0 ⎢⎣⎝ 0 ⎠ ⎥⎦
131
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
muestran los resultados obtenidos a partir de las ecuaciones (6.88) y (6.91) para
a los materiales estudiados por BJ (1967) cuyas propiedades se resumen en la Tabla 6-1.
mientras que lo inverso se muestra en la Figura 6-2b para dε β dy . Debe notarse que
1.0
1.0
εβω=0.9
0.8 0.8
0.8 εβω=0.1
0.7
0.6 EA,EB,EC 0.2 A2
dε β
0.6
0.6 0.3
r0 b) A1
εβ 0.5 A1
a) dy 0.4
0.4 0.4 0.5 EA,EB,EC
0.4 A2
0.3 0.6
0.2 0.7
0.2 0.2 0.8
0.1
0.9
0.0 0.0
-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0
y / r0 y / r0
ambos casos se identifican las curvas correspondientes a los materiales estudiados por
BJ (1967).
⎡ 3 ( ε βω − 1) ⎛ y 2 ⎞ ⎤
2
⎛ dε β ⎞
2
B ( y) = εβ ( y)⎜
−3
⎟ = εβ ( y) ⎢
−3
⎜ 2 − 1⎟ ⎥ , (6.92)
⎝ dy ⎠ ⎢⎣ 4 r0 ⎝ r0 ⎠ ⎥⎦
la cual está restringida a las suposiciones hechas para obtener la ecuación (6.88). La
132
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
9 ( ε βω − 1)
2
B ( 0) = . (6.93)
2r02 ( ε βω + 1)
3
Tabla 6-1: Datos de los materiales estudiados por Beavers y Joseph (BJ)
Figura 6-3. Puede notarse que, al aumentar la fracción porosa de la región −ω , los
esta forma, se aprecia que sólo para ε βω ≥ 0.9 el valor del esfuerzo es menor o igual al
esfuerzo en la superficie.
En la Figura 6-4 se muestran los resultados para el caso particular de los materiales
empleados por BJ. En este caso se puede apreciar la gran sensibilidad de los resultados
se necesita conocer el valor del área interfacial en la superficie ( avs ) , el cual se puede
av ,ω
avs = av ( 0 ) = , (6.94)
2
133
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
av ,ω =
6
lσ
(1 − ε βω ) . (6.95)
0.8
B/B(0)
6 0.6
0.2
4 0.4
0.3
2 0.2
0.4
0 0.0
-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0
y/r0 y/r0
avs =
3
lσ
(1 − ε βω ) . (6.96)
esfuerzos globales en la frontera. Para ello, primero se comparan varios modelos para el
acuerdo a su capacidad para reproducir los datos reportados por BJ. Y posteriormente se
A2
1.0
0.8
A1
0.6
B/Bmax
0.4
0.2
EA, EB, EC
0.0
-1.0 -0.5 0.0 0.5 1.0
y/r0
Figura 6-4: Dependencia espacial de los esfuerzos de Brinkman para los materiales
usados por BJ. La línea punteada horizontal indica el valor de B ( 0 ) / Bmax para cada
siguientes,
• Tubos capilares
ε βω lσ2
K βω = . (6.97)
32
ε βω
3 2
lσ
K βω = . (6.98)
180 (1 − ε βω )
2
135
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
ε βω
5.5 2
lσ
K βω = . (6.99)
5.6
ε βω
3 2
lσ
o Flujo paralelo K βω = . (6.100)
16kk (1 − ε βω )
2
o Flujo perpendicular
(l +δ ) ⎡ ⎛ l +δ ⎞ 1 ( lσ + δ ) − lσ4 ⎤
2 4
K βω = σ ⎢ln ⎜ σ ⎟− ⎥, (6.101)
⎠ 2 ( lσ + δ ) + lσ ⎥⎦
4
⎢⎣ ⎝ lσ
4
16
donde,
2ε βω
3
kk = . (6.102)
⎡ ⎛ ⎞ 2⎤
(1 − ε βω ) ⎢2 ln ⎜⎜ 1 −1ε ⎟⎟ − 3 + 4 (1 − ε βω ) − (1 − ε βω ) ⎥
⎣⎢ ⎝ βω ⎠ ⎦⎥
K βω ε βω
3
l 1
=A , = . (6.103)
(1 − ε ) 3 1− ε
2 2
l βω
lσ βω
K βω
= 0.0606 ε βω
5.1
; 0.4 ≤ ε βω ≤ 0.8 . (6.104)
lσ2
Como criterio de selección entre estos modelos, se decidió que sólo aquellos modelos
que reproduzcan satisfactoriamente los datos de permeabilidad de la Tabla 6-1, serán los
muestran en las figuras 6-5 y 6-6, las observaciones para los modelos capilares son las
siguientes:
datos experimentales.
136
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.00E-007 2.00E-007
a) b)
8.00E-008 1.60E-007
4.00E-008 8.00E-008
2.00E-008 4.00E-008
0.00E+000 0.00E+000
1 2 3 4 5 1 2 3 4 5
Material Material
Figura 6-5: Estimación del coeficiente de permeabilidad para los datos de BJ usando
satisfactoriamente todos los datos, mientras que su versión para cilindros sólo es
reproducen mejor los datos que para flujo paralelo, sin embargo los resultados
todos los datos de BJ, sin embargo en el caso de las dos especies de aloxita, el
Larson y Higdon no siendo éste el caso para las tres especies de espuma de
metal.
137
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.20E-007
1.80E-007 b)
1.60E-007
a) 1.05E-007
1.40E-007 9.00E-008
1.20E-007 7.50E-008
Kβω Kβω
1.00E-007 6.00E-008
8.00E-008
4.50E-008
6.00E-008
3.00E-008
4.00E-008
1.50E-008
2.00E-008
0.00E+000 0.00E+000
1 2 3 4 5 1 2 3 4 5
Material Material
Figura 6-6: Estimación del coeficiente de permeabilidad para los datos de BJ usando
para las tres especies de espuma de metal es recomendable usar el modelo de Happel y
Brenner para flujo perpendicular y para las dos especies de aloxita el modelo de Rumpf
medio poroso en su región homogénea. Sin embargo, no se encontró una relación que
que permite reducir los esfuerzos a cero en el límite cuando ε β → 1 . Esta propiedad la
138
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
mencionó, los únicos modelos (de los que se discutieron aquí) que reproducen
Higdon. Ambos permiten una transición suave en los esfuerzos al pasar de una región
K βω γε βω
3
= , (6.105)
(1 − ε βω )
2
lσ2
γ = γ 1 = (180 )
−1
donde lleva a la ecuación (6.98), mientras que si
Kβ ( y ) γε β3 ( y )
= , (6.106)
(1 − ε ( y ) )
2
lσ2 β
γ lσ2 (1 + ε βω )
3
Kβ = Kβ ( 0) = . (6.107)
2 (1 − ε βω )
ηω 2
valor en la superficie del medio poroso, se obtuvieron los resultados de la Figura 6-7. Al
máximo de los esfuerzos globales crece y llega a ser hasta 660 veces su valor en la
superficie, (para ε βω = 0.1 ). Por otro lado, su valor se vuelve despreciable para
y ≥ 0.5r0 .
139
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
700
14 εβω= 0.5
600 εβω= 0.1
12
500
10 0.6
K(0) / K
400
K(0) / K
8
0.7
300
6 0.8
200
0.4 4 0.9
100
0.2
2
0.3
0 0
-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0
y / r0 y / r0
Para los materiales empleados por BJ, los resultados se muestran en la Figura 6-8. La
ecuación (6.106) muestra que los esfuerzos globales son función tanto de la fracción
porosa como del diámetro de la esfera (o cilindro) que componen el medio poroso. El
normalizados de acuerdo a este valor. Al igual que en la Figura 6-4, los resultados son
bastante sensibles a cambios en las propiedades de los materiales. Cabe mencionar que
140
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0
A1
0.8
0.6
K / Kmax
0.4
EA EC
0.2
A2
EB
0.0
-1.0 -0.5 0.0 0.5 1.0
y / r0
Figura 6-8: Dependencia de los esfuerzos de globales con la posición para los
Con los resultados obtenidos en las dos secciones anteriores, se puede obtener la
2 (1 − ε βω ) 9 ( ε βω − 1)
2 2
−1
K = + , (6.108)
γ lσ2 (1 + ε βω ) 800 lσ2 ( ε + 1)3
3
βω
donde se usó la relación de (6.90). Debe ser claro, de la ecuación anterior, que los
2 (1 − ε βω )
2
−1
K = . (6.109)
γ lσ2 (1 + ε βω )
3
2ε βω
3/ 2
β= , (6.110)
3 γ (1 + ε βω )
3
141
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
La estimación del coeficiente de salto, para los materiales estudiados por BJ, se presenta
en la Tabla 6-2. Donde se muestran los resultados que se obtienen al usar los modelos
que los tres modelos ofrecen estimaciones del orden de la unidad para el coeficiente de
(1995a,b). Sin embargo, los resultados del modelo γ 2 se acercan más a los reportados
por Goyeau y col. (2003) para A1 y A2. Para los materiales EA, EB y EC los modelos
Tabla 6-2: Estimación del coeficiente de salto para los materiales usados por BJ
del coeficiente de salto con la porosidad y además la intersección de las curvas de los
límite cuando ε βω → 0 ; no siendo este el caso del modelo de Goyeau y col. (2003), en
ε βω → 1 , tanto los resultados que se obtienen con el método de Goyeau y col. (2003),
142
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
velocidad.
∂ vβ ∂ vβ K −1
η
ω
− ε βω = ε βω vβ , en la inter-región η − ω .
∂y ∂y avs ω
se aprecia que el término responsable del salto es K −1ε βω / avs , el cual puede expresarse,
K −1ε βω 2 (1 − ε βω ) ε βω
≈ . (6.111)
3γ lσ (1 + ε βω )
3
avs
100
10
Goyeau y col. (2003)
β 1
β (γ1)
0.1
β (γ2)
0.01
0.0 0.2 0.4 0.6 0.8 1.0
εβω
diámetro de las esferas (o cilindros) que componen el medio poroso. Esta última,
143
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
K −1ε βω / avs con la porosidad y l β se muestra en la Figura 6-10 para los resultados
frontera aumenta con l β . Además, la comparación de las Figuras 6-10a y 6-10b muestra
que la forma de las curvas difiere para ε βω > 0.3 . Como puede observase, para un
tamaño de poro dado, los modelos γ 1 y γ 2 predicen valores máximos de K −1ε βω / avs
ecuación (6.110). No está de más remarcar que estas observaciones están basadas en
7 −4 16 −4
lβ = 2 x 10 m lβ = 2 x 10 m b)
6 a) 14
10 x K εβω / avs ( m )
10 x K εβω / avs ( m )
12
−1
−1
5
−4
3 x 10 A1 10 3 x 10
−4
4 −4
4 x 10 A2 −4
8 4 x 10 A1 A2
3
−4
6 x 10 FA 6
6 x 10
−4
−1
−1
2 −4 FB −4 FA
8 x 10 4 8 x 10
FC FB
4
1
4
−3
2 −3
1.2 x 10 1.2 x 10 FC
0 0
0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0
εβω εβω
Figura 6-10: Dependencia del coeficiente K −1ε βω / avs con ε βω y el tamaño de poro.
Usando a) γ = γ 1 y b) γ = γ 2 .
A partir de lo desarrollado en estas secciones se puede concluir que no existe una forma
literatura ya que no se cuenta con una expresión general que relacione la permeabilidad
con la porosidad. Sin embargo, es posible confirmar que el coeficiente de salto tiene
y Whitaker (1995a,b).
término.
Con los resultados obtenidos, el siguiente paso en este trabajo consiste en estudiar el
en el siguiente capítulo.
145
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Parte III:
146
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 7:
como en la frontera del medio poroso. Un estudio similar a este ya ha sido reportado por
embargo, varios de los pasos del desarrollo de la condición de salto se harán siguiendo
las propuestas de Valencia-López y col. (2003). Por otro lado, el estudio de este
problema en el seno del medio poroso, empleando el método del promedio volumétrico,
fue llevado a cabo previamente por Whitaker (1986) y varias de las ideas de ese trabajo
serán retomadas en este capítulo. Asimismo, una buena parte de los desarrollos
partir del uso de estas ecuaciones, y las correspondientes a las regiones homogéneas, en
plantean los problemas de las variables de cerradura que permiten obtener la forma
cerrada de la condición de salto. El capítulo concluye con una breve discusión de los
resultados.
147
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
dentro y fuera de un medio poroso. Al igual que en los dos capítulos anteriores, se
denota a la fase fluida como la fase − β , mientras que la fase sólida es la fase −σ . La
diluida ( X A 1) es
∂c Aβ
+ ∇ i( c Aβ v β ) = ∇ i(Dβ ∇c Aβ ) , en la fase − β . (7.1)
∂t
∂c As
= −n βσ i(Dβ ∇c Aβ ) − ks cAs , en la interfase β − σ . (7.2)
∂t
donde c As es la concentración molar superficial de la especie A y k s es el coeficiente de
velocidad de reacción superficial. Para obtener una forma útil de la ecuación anterior, se
requiere contar con más información sobre las concentraciones molares interfacial y la
(k + k )
2
ks *
ε βω −1 s 1 ; ( k−1 + ks ) t 1, (7.4)
avω k−1k1 k −1
bajo las cuales es válida la suposición de equilibrio de adsorción local. De esta forma, la
concentración superficial puede expresarse en términos de la concentración en el seno
de la fase − β por medio de,
148
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
cAs = Kc Aβ , (7.5)
K eq k1
donde, K = , K eq = . Al sustituir la ecuación anterior en la ecuación (7.2)
1 + ks / k−1 k−1
resulta en
∂cAβ
K = −n βσ i(Dβ ∇c Aβ ) − kr c Aβ , en la interfase β − σ , (7.6)
∂t
en esta última, kr = ks K , es un coeficiente de pseudo-reacción (Whitaker, 1986). En
flujo (dentro y fuera del medio poroso) se supone estacionario, incompresible y con
las siguientes
∇ i v β = 0, en la fase − β , (7.7)
v β = 0, en la interfase β − σ . (7.9)
lleva a cabo el proceso de promediado, sin imponer restricciones de escala, para obtener
superficial,
1
ψβ
x
=
V ∫ ψβ
Vβ ( x )
x+y β
dV , (7.10)
149
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
las regiones homogéneas del sistema como en la inter-región (Figura 7-1); de esta
∂ cAβ
+ ∇ i( c Aβ v β ) = ∇ i(Dβ ∇c Aβ ) . (7.11)
∂t
Para obtener esta ecuación se supuso que Vβ ( x ) es independiente del tiempo. Como se
promediado.
Una aplicación del teorema del promediado espacial (Howes y Whitaker, 1985), en los
∂ cAβ
n βσ i(Dβ ∇c Aβ − v β c Aβ ) dA . (7.12)
1
∂t
+ ∇ i c Aβ v β = ∇ i Dβ ∇cAβ +
V ∫
Aβσ ( x )
Sin embargo, con base en las condiciones en la interfase sólido-fluido (7.6) y (7.9), la
∂ c Aβ ∂ c Aβ
+ ∇ i c Aβ v β = ∇ iDβ ∇c Aβ − kr av ( x ) c Aβ − Kav ( x ) βσ
, (7.13)
∂t βσ ∂t
área
1
A βσ ( x ) Aβσ∫( x )
c Aβ βσ
= c Aβ dA , (7.14)
A βσ ( x )
av ( x ) = . (7.15)
V
150
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Una segunda aplicación del teorema del promediado espacial, en el primer término del
∂ cAβ ⎡ ⎛ 1 ⎞⎤
∂t
+ ∇ i c Aβ v β = ∇ i ⎢Dβ ⎜ ∇ c Aβ +
⎢ ⎝ ⎜ V ∫ βσ Aβ ⎟⎟ ⎥⎥
n c d A
⎣ Aβσ ( x ) ⎠⎦ (7.16)
∂ c Aβ
−kr av ( x ) c Aβ − Kav ( x ) βσ
.
βσ ∂t
β 1
V β ( x ) Vβ∫( x ) β
ψβ = ψ x+yβ
dV , (7.17)
x
151
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( ) = ∇i⎡⎣⎢ε
∂ c Aβ β β
ε β (x) + ∇ i ε β ( x ) c Aβ v β ( x ) D ( x )i∇ cAβ ⎤
∂t
β
⎦⎥
(7.20)
∂ c Aβ
−kr av ( x ) cAβ − Kav ( x ) βσ
,
βσ ∂t
donde se introdujo el siguiente tensor de difusión dependiente de la posición,
β
D ( x )i∇ cAβ β 1 ⎛ β ⎞ dA ,
V β ( x ) ∫( ) βσ ⎝ β
= ∇ c Aβ + n ⎜c − c Aβ ⎟ (7.21)
Dβ ⎠
A x + yγ x
Aβσ x
(Gray, 1975)
β
c Aβ = c Aβ + c Aβ , (7.22)
x+y β x+y β x+y β
β
vβ x+y β
= vβ + vβ x+yβ
, (7.23)
x+y β
β
para poder expresar al término cAβ v β como sigue,
β β β β β
cAβ v β = vβ c Aβ + c Aβ v β + c Aβ v β ex
, (7.24)
β β
(7.25)
β β
+ c Aβ v β + c Aβ vβ .
x+y β x+yβ
( ) = ∇i⎡⎣⎢ε
∂ c Aβ
( x ) (D ( x ) + D*β ( x ) )i∇
β β β
ε β (x) + ∇i ε β ( x) vβ c Aβ cAβ ⎤
∂t
β
⎦⎥
(7.26)
∂ cAβ
−kr av ( x ) c Aβ − Kav ( x ) βσ
,
βσ ∂t
152
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Dβ ( x ) = D ( x ) + D*β ( x ) . (7.28)
( ) = ∇i⎡⎣⎢ε
∂ c Aβ β β β
ε β (x) + ∇i ε β ( x) vβ c Aβ ( x ) Dβ ( x )i∇ c Aβ ⎤
∂t
β
⎦⎥
(7.29)
∂ cAβ
−kr av ( x ) c Aβ − Kav ( x ) βσ
.
βσ ∂t
(
∇i ε β vβ
β
) = 0, (7.30)
β
0 = −∇ pβ + ρ β g + μ β ε β−1 ( x ) ∇ 2 v β
(7.31)
(
− μ β ε β−1 ( x ) ∇ε β ( x )i∇ ε β−1 ( x ) v β )−μ K β
−1
β ( x )i vβ ,
ecuación (5.22)). Conviene hacer notar que en la ecuación (7.31) se prefirió expresar al
153
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
y como consecuencia,
β 1 ⎛ β β ⎞
∇ c Aβ ∫
V β ( x h ) Aβσ ( xh )
n βσ ⎜ c Aβ
⎝ xh + y β
− c Aβ ⎟ dA,
xh ⎠ (7.33)
en las regiones homogéneas,
β β β β β
v β cAβ vβ cAβ , vβ c Aβ , en las regiones homogéneas . (7.34)
Región homogénea−ω
lβ LC , lβ r0 , (7.35)
β
La cual permite reemplazar cAβ
βσ
→ c Aβ ω
en la ecuación (7.29). Por otro lado, con
β ⎡ β 1 ⎤ β
Dβ ( xω )i∇ c Aβ = Dβ ⎢∇ c Aβ + ∫ n βσ cAβ dA⎥ − v β c Aβ . (7.37)
ω
⎢⎣ ω V β ( xω ) Aβσ ( xω ) ⎥⎦ ω
β
⎛ Ka
ε β ,ω ⎜1 + v ,ω
⎜
⎝ ε β ,ω
⎞ ∂ c Aβ
⎟⎟
⎠ ∂t
ω
(
+ ∇ i ε β ,ω v β
β
ω
c Aβ
β
ω ) (7.38)
β β
= ∇ i ⎡⎢ε β ,ωDβ ( xω )i∇ c Aβ ⎤−k a c
⎥⎦ r v ,ω Aβ , en la región homogéna − ω.
⎣ ω ω
154
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Región homogénea−η
Dado que en esta región ( x → xη ) la única fase presente es la fase fluida (Figura 7-1),
D β ( xη )i∇ c Aβ
β β
= Dβ ∇ c Aβ , (7.39)
η η
β
∂ cAβ
∂t
η
(
+ ∇i vβ
β
η
c Aβ
β
η ) = ∇i(D ∇ c ) , en la región homogénea −η .
β Aβ η
β
(7.40)
La cual tiene la misma estructura que la ecuación (7.1). El procedimiento para obtener,
puede proceder al desarrollo de las condiciones de salto que acoplan las ecuaciones
(7.38) y (7.40).
V ∞ =V η +V ω , (7.41)
A ∞ = Aη +A ω . (7.42)
β β
∂ cAβ ∂ cAβ
∫ ε β ( x)
Vω
∂t
dV + ∫ ε β ( x )
Vη
∂t
dV
Aω
(
+ ∫ nω i ε β ( x ) v β
β
cAβ
β
) dA + ∫ n i(εAη
η β (x) vβ
β
c Aβ
β
) dA
(7.43)
β β
= ∫ nω i ⎡⎢ε β ( x ) Dβ ( x )i∇ c Aβ ⎤ dA + n i ⎡ε ( x ) D ( x )i∇ c
∫ η ⎢⎣ β ⎤ dA
Aω
⎣ ⎦⎥ A
β Aβ
⎦⎥
η
∂ cAβ
− ∫ kr av ( x ) c Aβ dV − ∫ Kav ( x ) βσ
dV .
V∞
βσ
V∞
∂t
β β
∂ c Aβ
) dA
∂ cAβ
∫ nω i(ε β ω
β β
∫ εβ ω dV + ∫ η
ω
dV + vβ c Aβ
∂t ∂t
, , ω ω
Vω Vη Aω
+ ∫ nη i v β
Aη
( β
η
c Aβ
β
η ) dA + ∫ n i(ε
*
Aωη
ωη β ,ω vβ
β
ω
c Aβ
β
ω
− vβ
β
η
cAβ
β
η ) dA
= ∫ nω i⎡⎣⎢ε β ωDβ ( xω )i∇
Aω
, cAβ
β
⎤ dA + n i D ∇ c
⎥
ω⎦ ∫ η β Aβ Aη
( β
η ) dA (7.44)
β β
+ ∫ nωη i⎡⎢⎣ε β ωDβ ( xω )i∇
*
, c Aβ ω
− Dβ ∇ c Aβ η
⎤ dA
⎥⎦
Aωη
β
β ∂ c Aβ
− ∫ kr av ,ω cAβ dV − ∫ Kav ,ω ω
dV .
Vω
ω
Vω
∂t
156
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
⎛ ∂ cAβ
β
∂ cAβ
β
⎞ ⎛ ∂ cAβ
β
∂ c Aβ ⎞
⎟ dV + ⎜ ε ( x ) ⎟ dV
∫V ⎜⎜ ε β ( x ) ∂t ∫
η
− ε β ,ω ω
−
∂t ⎟ ⎜ β
∂t ∂t ⎟⎟
Vη ⎜
ω
⎝ ⎠ ⎝ ⎠
(
+ ∫ nω i ε β ( x ) v β
Aω
β
cAβ
β
− ε β ,ω v β
β
ω ) dAcAβ
β
ω
(
+ ∫ nη i ε β ( x ) v β
Aη
β
c Aβ
β
− v c
β η ) dA
β
Aβ η
β
= ∫ nωη i(ε β ω
*
Aωη
, vβ
β
ω
cAβ ω
β
− v cβ η ) dA
β
Aβ η
β
+ ∫ nωη i(Dβ ∇
*
Aωη
c Aβ
β
η
− ε β ,ω D ( x )i∇ c
β ω ) dA Aβ ω
β
β β
+ ∫ nω i ⎡⎢ε β ( x ) Dβ ( x )i∇ cAβ − ε β ,ωD β ( xω )i∇ c Aβ ⎤ dA
⎥⎦
A
⎣ ω
ω
β β
+ ∫ nη i ⎡⎢ε β ( x ) Dβ ( x )i∇ c Aβ − Dβ ∇ c Aβ ⎤ dA
A
⎣ η ⎦⎥
η
β
− ∫ kr av ( x ) cAβ βσ
dV + ∫ kr av ,ω c Aβ ω
dV
V∞ Vω
(7.45)
β
∂ cAβ ∂ c Aβ
− ∫ Kav ( x ) βσ
dV + ∫ Ka ω ω
dV .
∂t ∂t
v,
V∞ Vω
y col., 2003),
β β
∂ c Aβ ⎡ ∂ cAβ
β
∂ cAβ ω ⎤
∫ εβ s
dA = ∫ ε β ( x )
⎢ − ε β ,ω ⎥ dV
s
∂t ⎢ ∂t ∂t ⎥
⎣ ⎦
* Vω
Aωη
β
(7.46)
⎡ ∂ c Aβ
β
∂ cAβ η ⎤
⎢
+ ∫ ε β (x) − ⎥ dV ,
⎢ ∂t ∂t ⎥
Vη
⎢⎣ ⎥⎦
∫ n i(ε β
C
s s vβ
β
s
c Aβ
β
s ) dσ = ∫ n i(ε ( x ) v
Aω
ω β β
β
cAβ
β
− ε β ,ω v β ω
β
cAβ
β
ω ) dA
+ ∫ n i( ε ( x ) v ) dA,
β β β β
(7.47)
η β β c Aβ − vβ η
cAβ η
Aη
157
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
∫ n i⎡⎢⎣ε β D i∇
C
s s s s c Aβ s
⎤ dσ
⎥⎦
β β
= ∫ nω i⎡⎣⎢ε β ( x ) Dβ ( x )i∇
Aω
c Aβ − ε β ,ω Dβ ( xω )i∇ cAβ ω
⎤ dA
⎦⎥
(7.48)
β β
+ ∫ nη i ⎡⎢ε β ( x ) Dβ ( x )i∇ c Aβ − Dβ ∇ c Aβ ⎤ dA,
⎥⎦
A
⎣ η
η
β
∂ c Aβ ∂ c Aβ
∫ ∫ Ka ( x ) dV − ∫ Kav ,ω
βσ
RAds dA = ω
dV . (7.50)
∂t ∂t
v
* V∞ Vω
Aωη
y Ds es el tensor de exceso de dispersión total superficial. Por otro lado, a partir de las
la Sección 2.6)
1
RRx s
=
A ωη ∫ kcβ
Aωη
r A s
dA = kr cAβ ηω
, (7.51)
1 ∂c Aβ ∂c Aβ
RAds s
=
A ωη ∫ K
∂t
dA = K
∂t
. (7.52)
Aωη s ηω
(Deemer y Slattery, 1978; Slattery, 1980; Ochoa-Tapia y col., 1993; Slattery, 2007), la
158
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
( )
∂ c Aβ β β β
εβs s
+ nωη i ε β s v β cAβ − ε β sDs i∇ s c Aβ
∂t s s s
(
−nωη i ε β ,ω v β
β
ω
c Aβ
β
ω
− vβ
β
η
cAβ
β
η )
transporte convectivo
β β
(7.53)
= −nωη i ⎡⎢ε β ,ωDβ ( xω )i∇ cAβ − Dβ ∇ c Aβ ⎤
η ⎥
⎣ ω ⎦
transporte difusivo
∂c Aβ
− kr cAβ − K , en la superficie divisoria.
ηω ∂t ηω
reacción superficial
adsorción superficial
Por otro lado, si se aceptan las suposiciones y restricciones de escala introducidas en las
β β
−nωη i ⎡⎢ε β ,ωD β ( xω )i∇ c Aβ − Dβ ∇ cAβ ⎤
⎣ ω η ⎦⎥
(
+nωη i ε β ,ω v β
β
ω
c Aβ
β
ω
− vβ
β
η
cAβ
β
η ) (7.54)
∂cAβ
= kr c Aβ + K , en la superficie divisoria.
ηω ∂t ηω
No está de más remarcar que, a diferencia del trabajo de Valencia-López y col. (2003),
de concentraciones promedio,
β β β
cAβ = cAβ = c Aβ , en la superficie divisoria . (7.55)
ω η s
159
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
(
nωη i v β ω
− vβ η ) = 0, en la superficie divisoria , (7.56)
−nωη (p β
β ω − pβ η
β
) + μ n i(ε
β ωη
−1
βω ∇ vβ ω
− ∇ vβ η ) = aμvs
β
K −1 i v β ω (7.57)
en la superficie divisoria.
ecuaciones (7.55) y (7.56) en la condición de salto (7.54), permite que se cancelen los
reduce a,
∂cAβ β β
K = −nωη i ⎡⎢ε β ,ωDβ ( xω )i∇ c Aβ − Dβ ∇ cAβ ⎤
∂t ηω
⎣ ω η ⎦⎥ (7.58)
− kr c Aβ ηω
, en la superficie divisoria.
obtenida en el Capítulo 2 para el caso de difusión y reacción (ver ecuación (2.47)). Sin
kr = kr ηω
+ kr , (7.59)
K= K ηω
+K , (7.60)
160
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
que permiten descomponer los términos de reacción y adsorción superficial como sigue,
β β
kcAβ ηω
= k ηω
c Aβ + k ηω
c Aβ + k c Aβ + kcAβ , (7.61)
ηω
β β
∂cAβ ∂ cAβ ∂c Aβ ∂ c Aβ ∂c Aβ
K = K ηω
+ K ηω
+K +K . (7.62)
∂t ηω
∂t ∂t ∂t ∂t
ηω
β
∂c Aβ ∂ c Aβ ∂cAβ
K = K ηω
+ K . (7.64)
∂t ηω
∂t ∂t ηω
concentración. Por lo que se debe obtener el problema de valor a la frontera, a partir del
(
= ∇ i Dβ ∇ cAβ
β
)+ε −1
β ( x ) ∇ε β ( x )i Dβ ∇ ( c Aβ
β
)
∂ cAβ
−ε β−1 ( x ) kr av ( x ) c Aβ − ε β−1 ( x ) Kav ( x ) βσ
(7.65)
βσ ∂t
⎡D ⎤
+ε β−1 ( x ) ∇ i ⎢ β ∫ n βσ ⎛⎜ c Aβ − c Aβ ⎞ dA ⎥
β
⎟
⎢⎣V Aβσ ( x ) ⎝ x⎠ ⎥⎦
−ε β−1 ( x ) ∇ i ⎡⎢ cAβ v β + c Aβ v β ex ⎤⎥ .
⎣ ⎦
∂c Aβ
= ∇ i(Dβ ∇c Aβ )
β β
+ v β i∇c Aβ − v β ∇ i c Aβ
∂t
(
= ∇ i Dβ ∇ cAβ
β
)−ε −1
β ( x ) ∇ε β ( x )i Dβ ∇ ( c Aβ
β
)
161
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂ cAβ
+ε β−1 ( x ) kr av ( x ) c Aβ + ε β−1 ( x ) Kav ( x ) βσ
βσ ∂t
⎡D ⎤
n βσ ⎛⎜ c Aβ − c Aβ ⎞ dA ⎥
β
−ε β−1 ( x ) ∇ i ⎢ β ∫ ⎟
⎢⎣V Aβσ ( x )
⎝ x⎠ ⎥⎦ (7.66)
+ε β−1 ( x ) ∇ i ⎡⎢ c Aβ v β + c Aβ v β ⎤ , en la fase − β .
⎥⎦
⎣ ex
Dβ t *
1 (7.67)
δ β2
v β i∇c Aβ + v β ∇ i cAβ
β
= Dβ ∇ 2 cAβ − ε β−1 ( x ) ∇ε β ( x )i Dβ ∇ c Aβ ( β
)
∂ c Aβ
+ε β−1 ( x ) kr av ( x ) c Aβ + ε β−1 ( x ) Kav ( x ) βσ
βσ ∂t
⎡D ⎤
−ε β−1 ( x ) ∇ i ⎢ β ∫ n βσ cAβ dA⎥ (7.68)
⎢⎣V Aβσ ( x ) ⎥⎦
⎡D ⎛ β β ⎞ ⎤
−ε β−1 ( x ) ∇ i ⎢ β ∫ n βσ ⎜ c Aβ − c Aβ ⎟ dA⎥
⎢⎣V Aβσ ( x ) ⎝ x+y β x⎠ ⎥⎦
+ε β−1 ( x ) ∇ i ⎡⎢ c Aβ v β + c Aβ v β ex ⎤⎥ , en la fase − β .
⎣ ⎦
( )−n
∂ c Aβ
i(Dβ ∇c Aβ )
β
K = −n βσ i Dβ ∇ c Aβ βσ
∂t (7.69)
β
− kr cAβ − kr cAβ , en la interfase β − σ ,
162
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
⎡ ⎤
⎢ ⎥
⎢ k δ /D δ β / LC K δ β / Dβ t * ⎥ β
cAβ =O⎢ r β β , , c Aβ . (7.70)
⎛kδ ⎞ ⎛ krδ β ⎞ ⎛ krδ β ⎞ ⎥
⎢1 + O ⎜ r β ⎟⎟ 1 + O ⎜⎜ ⎟⎟ 1 + O ⎜⎜ ⎟⎟ ⎥
⎢⎣ ⎜
⎝ Dβ ⎠ D
⎝ β ⎠ ⎝ β ⎠ ⎥⎦
D
β
cAβ x+y β
cAβ , (7.71)
x+y β
krδ β δβ Kδ β
1; 1; 1. (7.72)
Dβ LC Dβ t *
( )−n
∂ c Aβ
i(Dβ ∇c Aβ )
β
K = −n βσ i Dβ ∇ cAβ βσ
∂t (7.73)
β
− kr c Aβ , en la interfase β − σ .
siguientes desigualdades
⎡ 1 ⎤
∇ 2 cAβ ε β−1 ( x ) ∇ i ⎢ ∫ βσ Aβ ⎥⎥ ,
n c d A (7.74)
⎢⎣V Aβσ ( x ) ⎦
β β
∫ n βσ i∇ cAβ dA ∫ ε β−1 ( x ) ∇ε β ( x )i∇ cAβ dV , (7.75)
Aβσ ( x ) Vβ ( x )
permite que v β = O v β ( β
) . Por lo que, la ecuación (7.68) se reduce, al tomar en
cuenta la desigualdad en (7.71), a
163
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β
v β i∇c Aβ + v β ∇ i cAβ = Dβ ∇ 2 cAβ
x
⎡ β ∂ β ⎤
+ε β−1 ( x ) av ( x ) ⎢ kr cAβ +K c Aβ ⎥
⎣⎢ x+yβ
βσ
∂t x+ yβ ⎥
βσ ⎦
(7.77)
⎡D ⎛ β β ⎞ ⎤
−ε β−1 ( x ) ∇ i ⎢ β ∫ n βσ ⎜ c Aβ − cAβ ⎟ dA⎥
⎢⎣V Aβσ ( x ) ⎝ x+yβ x⎠ ⎥⎦
+ε β−1 ( x ) ∇ i c Aβ v β ex
, en la fase − β .
En este punto, conviene encontrar qué tipo de restricciones de escala deben satisfacerse
β
para que la propiedades ψ β puedan suponerse constantes en el volumen de
x+y β
de x,
β β β 1 β
ψβ = ψβ + y β i∇ ψ β + y β y β : ∇∇ ψ β + ... (7.78)
x+yβ x x 2 x
El siguiente paso consiste en sustituir esta expansión en las ecuaciones (7.73) y (7.77);
analizar término por término, de esta forma, los términos de reacción y adsorción se
β ∂ β
kr c Aβ +K c Aβ
x+yβ
βσ
∂t x+y β
βσ
β
β β β
∂ cAβ
k
= k cAβ + ky β i∇ c Aβ + y β y β : ∇∇ c Aβ +K x
(7.79)
x x 2 x ∂t
∂⎛ β ⎞ + K ∂ ⎛ 1 y y : ∇∇ c β ⎞
+K ⎜ y β i∇ c Aβ ⎟ ⎜ β β Aβ ⎟ + ...
∂t ⎝ x⎠ ∂t ⎝ 2 x⎠
⎡D ⎛ β β ⎞ ⎤
ε β−1 ( x ) ∇ i ⎢ β ∫ n βσ ⎜ c Aβ − c Aβ ⎟ d A ⎥=
⎢V
⎣ Aβσ ( x ) ⎝ x+y β x⎠ ⎥⎦
( )
β
−ε β−1 ( x ) ∇ i Dβ ∇ y β i∇ cAβ (7.80)
x
⎡D ⎤ β
−ε β−1 ( x ) ∇ i ⎢ β ∇ y β y β ⎥ : ∇∇ c Aβ + ...
⎣ 2 ⎦ x
164
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
y Whitaker (1994)
1
V ∫
Aβσ ( x )
n βσ y β dA = −∇ y β , (7.81)
1
V ∫
Aβσ ( x )
n βσ y β y β dA = −∇ y β y β . (7.82)
Por último, los términos que componen el término de transporte dispersivo no-local se
β β ⎛1 ⎞ β
cAβ vβ = vβ c Aβ + v β ⎜ y β y β ⎟ : ∇∇ c Aβ + ... , (7.83)
x+yβ x ⎝ 2 ⎠ x
β β ⎛1 ⎞ β
vβ cAβ = cAβ vβ + c Aβ ⎜ y β y β ⎟ : ∇∇ v β + ... , (7.84)
x+yβ x ⎝2 ⎠ x
β β ⎛ β β 1 β ⎞ β
cAβ vβ = ⎜ c Aβ + y β i∇ cAβ + y β y β : ∇∇ cAβ ⎟ vβ
x+y β x+yβ ⎝ x x 2 x⎠ x
⎛ β β 1 β ⎞⎛ β ⎞
+ ⎜ c Aβ + y β i∇ cAβ + y β y β : ∇∇ cAβ ⎟ ⎜ y β i∇ v β ⎟ (7.85)
⎝ x x 2 x ⎠⎝ x⎠
⎛ β β 1 β ⎞⎛ 1 β ⎞
+ ⎜ c Aβ + y β i∇ cAβ + y β y β : ∇∇ cAβ ⎟⎜ y β y β : ∇∇ v β ⎟ + ...
⎝ x x 2 x ⎠⎝ 2 x⎠
restricciones de escala,
2 2
10ε β L2ε 10ε β L2ε LC1
r0 L, r L, 1; 1, (7.86)
0
(1 − ε ) δβω r
β 0 (1 + ε ) r δ
βω
2
0 β
β β
v β i∇c Aβ + v β ∇ i c Aβ = Dβ ∇ 2 c Aβ + ε β−1 ( x ) kr av ( x ) cAβ
x x
fuente convectiva fuente de reacción volumétrica
β (7.87)
∂ c Aβ
+ ε β ( x ) Kav ( x )
−1 x
, en la fase − β .
∂t
fuente de adsorción volumétrica
165
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Debe ser claro que no se pretende resolver el problema de las desviaciones espaciales de
la concentración en toda la superficie divisora entre las regiones, sino en una porción
representativa como la mostrada en la Figura 7-3. Como en los casos anteriores, se debe
seleccionar la altura h de la celda unitaria de manera tal que, en sus extremos, las
fuentes como,
β
β β ∂ c Aβ
cAβ = sβ cAβ + b β i∇ cAβ + dβ . (7.93)
∂t
La cual, da lugar a los siguientes tres problemas de valor a la frontera para las variables
de cerradura,
Problema I
kr
C.F.1 −n βσ i∇sβ = , en la interfase β − σ , (7.95)
Dβ
C.F.2 En ζ = h, sβ = 0 , (7.96)
∂sβ ∂sω
C.F.3 En ζ = − h, sβ = sω , = , (7.97)
∂ζ ∂ζ
Periodicidad: sβ ( r + li ) = sβ ( r ) , i = 1, 2 . (7.98)
Problema II
C.F.2 En ζ = h, bβ = 0 , (7.101)
∂b β ∂bω
C.F.3 En ζ = − h, b β = bω , = , (7.102)
∂ζ ∂ζ
Periodicidad: b β ( r + li ) = b β ( r ) , i = 1, 2 . (7.103)
Problema III
167
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
C.F.2 En ζ = h, dβ = 0 , (7.106)
∂d β ∂dω
C.F.3 En ζ = − h, d β = dω , = , (7.107)
∂ζ ∂ζ
Periodicidad: d β ( r + li ) = d β ( r ) , i = 1, 2 . (7.108)
Nótese que los problemas I y III tienen la misma estructura y que los tres problemas de
valor a la frontera dependen del campo de velocidades. Para obtener las condiciones de
β
β β ∂ c Aβ
cAω = b βω i∇ c Aβ + sβω c Aβ + dω ω
. (7.109)
ω ω ∂t
β β
β β ∂ c Aβ ∂ c Aβ β β
En ζ = −h, c Aβ = c Aβ ; = ω
; ∇ c Aβ = ∇ c Aβ . (7.110)
ω ∂t ∂t ω
sección.
En la sección anterior se obtuvo una ecuación que permite expresar el campo de las
168
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
β β
kr c Aβ ηω
= kr ηω
cAβ + kr b β i∇ cAβ
s ηω s
β (7.111)
β ∂ c Aβ
+ k r sβ c Aβ + kr d β s
,
ηω s ηω ∂t
K
∂C Aβ
= K
∂ cAβ
β
s
+ Kb β i
(
∂ ∇ cAβ
β
s )
ηω
∂t ηω
∂t ηω ∂t
(7.112)
β β
∂ c Aβ
2
∂ cAβ
+ Kd β s
+ Ksβ s
.
ηω ∂t 2 ηω ∂t
Sin embargo, a partir de los problemas II y III se pueden obtener los siguientes
b β = O (δ β ) , (7.113)
⎛ Kδ β ⎞
dβ = O ⎜ ⎟⎟ . (7.114)
⎜ D
⎝ β ⎠
Kδ β
De manera que, si las restricciones de escala δ β L y 1 se satisfacen,
Dβ t *
β
kr c Aβ = K effrx c Aβ , (7.115)
ηω ω
β
∂C Aβ ∂ cAβ
K =K ads ω
, (7.116)
∂t ∂t
eff
ηω
K effrx = kr ηω
+ k r sβ , (7.117)
ηω
K effads = K ηω
+ Ksβ . (7.118)
ηω
169
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
A partir de las ecuaciones (7.117) y (7.118) se deduce que sólo es necesario resolver el
Problema I para calcular los coeficientes efectivos. Nótese que en el caso en que no se
K effads = K ηω
, caso no reactivo. (7.119)
Región homogénea−ω
β
⎛ Ka
ε β ,ω ⎜⎜1 + v ,ω
⎝ ε β ,ω
⎞ ∂ c Aβ ω
⎟⎟
⎠ ∂t
+ ∇ i v β ω cAβ ω (
β
) (7.120)
β β
= ∇ i ⎡⎢ε β ,ωDβ ( xω )i∇ cAβ ω ⎤⎥ − kr av ,ω cAβ ,
⎣ ⎦ ω
∇i vβ = 0, (7.121)
ω
vβ ω
=−
K βω
μβ (
i ∇ pβ
β
ω
− ρ β g − ε βω
−1
μβ ∇ 2 v β ω ). (7.122)
Región homogénea−η
β
∂ cAβ
∂t
η
+ ∇i vβ ( η
cAβ
β
η ) =D ∇ β
2
cAβ
β
η
, (7.123)
∇i vβ η
= 0, (7.124)
β
0 = −∇ p β η
+ ρβ g + μβ ∇2 v β η
. (7.125)
Inter-región ω−η
β β
cAβ ω
= cAβ η
, (7.126)
β
∂ c Aβ β β β
K effads ω
= −nωη i ⎡⎢ε β ,ωDβ ( xω )i∇ cAβ − Dβ ∇ c Aβ ⎤ − K rx c
⎥⎦ Aβ , (7.127)
∂t ⎣ ω η eff ω
(
nωη i v β ω
− vβ η ) = 0, (7.128)
−nωη (p β ω
β
− pβ
β
η ) + μ n i(ε
β ωη
−1
βω ∇ vβ ω
− ∇ vβ η ) = μa β
vs
K −1 i v β ω
. (7.129)
170
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Con esto se termina el análisis. No sin antes enfatizar que, a pesar que la estructura, en
depende del campo de velocidad local. Estas diferencias dejarán de ser significativas
sólo en el caso en que los efectos convectivos sean despreciables respecto a los
difusivos (es decir, cuando el número de Péclet sea mucho menor que la unidad).
Para finalizar este capítulo, es oportuno comentar que este tipo de análisis puede
poroso-fluido. El análisis de estos y otros sistemas, así como la solución del Problema I,
171
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Conclusiones
En este trabajo, se desarrolló una metodología que permite obtener condiciones de salto
Joseph (1967).
Whitaker (1995a) y Wood y col. (2000). Los pasos más importantes que conforman esta
iii) Determinar la relación entre las propiedades de exceso y las locales, en este
El contar con esta metodología permite plantear, usando el método del promedio
172
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Por otro lado, en cada capítulo se incluye una sección para conclusiones y discusión de
resultados, sin embargo, resulta conveniente retomar aquí las conclusiones a las que se
la frontera. Las diferencias entre estos resultados y los obtenidos con un modelos
Esto permitió comparar los resultados con los de Goyeau y col. (2003),
173
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Capítulo 3.
columnas cromatográficas, sistemas de tres fases, entre otras. Esto se llevará a cabo en
trabajos futuros.
174
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Bibliografía
Alazmi, B., Vafai, K., 2001. “Analysis of fluid flow and heat transfer interfacial
Anderson, T.B.; Jackson, R., 1967. “A fluid mechanical description of fluidized beds”
Arquis, E., Caltagirone, J.P., 1984. “Sur les conditions hydrodynamiques au voisinage
Barrère, J., Gipouloux, O., Whitaker, S. 1992. “On the closure problem for Darcy's law”
Beavers, G.S., Joseph, D.D., 1967. “Boundary conditions at a naturally permeable wall”
Beckermann, C., Ramadhyani, S., Viskanta, R., 1987. “Natural convection flow and
heat transfer between a fluid and a porous layer inside a rectangular enclosure”
Beckermann, C., Viskanta, R., Ramadhyani, S., 1988. “Natural convection in vertical
Bhattacharyya, A., Sekhar, G.P.R., 2005. “Stokes flow inside a porous spherical shell:
175
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Bletzacker, L.; Prat, M.; Quintard, M. 2003., “Transport réactif en milieux poreux:
Brusseau, M.L., Larsen, T., Christensen, T.H., 1991. “Rate-limited sorption and non-
Carman, P.C., 1937. “The determination of the specific surface area of powder I”
Chang, H.C., 1983, “Effective diffusion and conduction in two-phase media: a unified
Coppens, M.O. (1999). “The effect of fractal surface roughness on diffusion and
Crapiste, G.H., Rotstein, E., Whitaker, S., 1986. “A general closure scheme for the
method of volume averaging” Chemical Engineering Science, Vol. 41, No. 2, 227-
235.
Dalvie, S.K., Gajiwala, K.S., Baltus, R.E., 1990. “Mathematical model of a rotating
Deemer, A.R., Slattery, J.C., 1978. “Balance equations and structural models for phase
176
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Deng, C., Martinez, D.M., 2005. “Viscous flow in a channel partially filled with a
porous medium and with wall suction” Chemical Engineering Science, Vol. 60,
329-336.
DeVault, D., 1943. “The theory of chromatography” Journal of the American Chemical
Dullien, F. A. L., 1979, Porous media: Fluid transport and pore structure. Academic
press.
Eidsath, A.B., Carbonell, R.G., Whitaker, S., Herrmann, L.R., 1983, “Dispersion in
pulsed systems III: Comparison between theory and experiments for packed beds”
Fletcher, C.A.J., 1991. Computational Techniques for Fluid Dynamics, second edition,
Gibbs, J.W., 1928. The Collected works of J. Willard Gibbs, vol. 1, Yale University
Givler, R.C., Altobelli, S.A., 1994. “A determination of the effective viscosity for the
370.
Gobin, D., Goyeau, B., Songbe, J.-P., 1998. “Double diffusive natural convection in a
composite fluid-porous layer” ASME Journal of Heat Transfer, Vol. 120, 234-242.
Goto, M., McCoy, B.J. 2000. “Inverse size-exclusion chromatography for distributed
Goyeau, B., Lhuillier, D., Gobin, D., Velarde, M.G., 2003. “Momentum transport at a
fluid-porous interface” International Journal of Heat and Mass Transfer, Vol. 46,
4071-4081.
177
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
masa entre un fluido y un medio poroso” Memorias del XXV Encuentro nacional
de la AMIDIQ.
Goharzadeh, A., Saidi, A., Khalili, A., 2005. “Small scale measurement and theoretical
Gray, W.G., 1975. “A derivation of the equations for multiphase transport” Chemical
Haber, S., Mauri, R., 1983. “Boundary conditions for Darcy´s flow through porous
Hager, J., Whitaker, S., 2000. “Vapor-liquid jump condition within a porous medium:
results for mass and energy” Transport in Porous Media, Vol. 40, 73-111.
Happel, J., Brenner, H., 1991. Low Reynolds number hydrodynamics. Kluwer academic
publishers.
Howes, F.A, Whitaker, S., 1985. “The spatial averaging theorem revisited” Chemical
Jones, I.P., 1973. “Low Reynolds number flow past a porous spherical shell”
Kakac, S., Vasiliev, L. L., Yener, Y., 2005. Microscale heat transfer. Fundamentals
Karniadakis, G.E., Beskok, A., 2005. Microflows and Nanoflows: fundamentals and
Kaviany, M., 1999. Principles of heat transfer in porous media, 2nd edition, Springer-
New York.
Kim, J.H., Ochoa-Tapia, J.A., Whitaker, S., 1987. “Diffusion in anisotropic porous
178
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Kolodziej, J. A., 1988. “Influence of the porosity of a porous medium on the effective
Kuznetsov, A.V., 1996. “Analytical investigation of the fluid flow in the interface
region between a porous medium and a clear fluid in channels partially filled with
Kuznetsov, A.V., 1997. “Influence of the stress jump condition at the porous-
Lichtner, P.C., Tartakovsky, D.M., 2003. “Stochastic analysis of effective rate constant
Martys, N., Bentz, D.P., Garboczi, E.J., (1994). “Computer simulation study of the
Min, Y.J., Kim, S.J., 2005. “A novel methodology for thermal analysis of a composite
system consisting of a porous medium and an adjacent fluid layer” ASME Journal
Neale, G., Nader, W., 1974. “Practical significance of Brinkman´s extension of Darcy´s
law: Coupled parallel flows within a channel and a bounding porous medium”
Ochoa-Tapia, J.A., Del Río, J.A., Whitaker, S., 1993. “Bulk and surface diffusion in
179
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Ochoa-Tapia, J.A., 1994. “Un estudio comparativo de dos modelos para el transporte de
101-106.
Ochoa-Tapia J.A., Whitaker, S., 1995a. “Momentum transfer at the boundary between a
International Journal of Heat and Mass Transfer, Vol. 38, No. 14, 2635-2646.
Ochoa-Tapia, J.A., Whitaker, S., 1995b. “Momentum transfer at the boundary between
International Journal of Heat and Mass Transfer, Vol. 38, No. 14, 2647-2655.
Ochoa-Tapia, J.A., Soria, A., 1995. “Emulsiones dobles: condiciones de salto para el
Ochoa-Tapia, J.A., Whitaker, S. 1997. “Heat transfer at the boundary between a porous
Ochoa-Tapia, J.A., Whitaker, S. 1998. “Heat transfer at the boundary between a porous
Partha, M.K., Murthy, P.S.V.N., Sekhar, R. 2005. “Viscous flow past a porous spherical
180
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Tecnológico de Celaya.
Prat, M., 1989. “On the boundary conditions at the macroscopic level” Transport in
Prat, M., 1990. “Modeling of heat transfer by conduction in a transition region between
a porous medium and an external fluid” Transport in Porous Media, Vol. 5, 71-95
Prat, M., 1992. “Some refinements concerning the boundary conditions at the
Quintard, M., Whitaker, S., 1994. “Transport in ordered and disordered porous media II:
Reis, J.F.G., Lightfoot, E.N., Noble, P.T., Chiang, A.S., 1979. “Chromatography in a
Richardson, S., 1971. “A model for the boundary condition of a porous material. Part 2”
181
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Ryan, D., Carbonell, R.G., Whitaker, S., 1981. “A theory of diffusion and reaction in
Saffman, P.G., 1971. “On the boundary condition at the surface of a porous medium”
Sahraoui, M., Kaviany, M., 1992. “Slip and no-slip velocity boundary conditions at the
Sahraoui, M., Kaviany, M., 1993. “Slip and no-slip temperature boundary conditions at
Sahraoui, M., Kaviany, M., 1994. “Slip and no-slip temperature boundary conditions at
Scheidegger, A.E., 1974. The physics of flow through porous media. 3rd edition. Univ.
Slattery, J.C., 1967. “Flow of viscoelastic fluids through porous media” AIChE Journal,
Slattery, J.C., Sagis, L., Oh, E.-S., 2007. Interfacial transport phenomena. 2nd edition.
Springer-Verlag.
Soria-López, A., Whitaker, S., Ochoa-Tapia, J.A., 1997. “The mass transfer jump
condition at the boundary between the liquid membrane and the external regions
Sparrow, E.M., Loeffler, A.L., 1959. “Longitudinal laminar flow between cylinders
182
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
disk near its surface” Zeitschrift für angewandte Mathematik und Physik, Vol. 50,
529-545.
Starov V.M.; Zhdanov V.G. (2001). “Effective viscosity and permeability of porous
192, 363-375.
A2, 111-112.
Taylor, G.I., 1971. “A model for the boundary condition of a porous material. Part 1”.
Wakao, N., Smith, J.M., 1964. “Diffusion and reaction in porous catalysts”. Industrial
Whitaker, S., 1967. “Diffusion and dispersion in porous media”. AIChE Journal, Vol.
13, 420-427.
Whitaker, S., 1969. “Advances in Theory of Fluid Motion in Porous Media” Industrial
Whitaker, S., 1986. “Transient diffusion, adsorption and reaction in porous catalysts:
183
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Whitaker, S., 1991. “Improved constraints for the principle of local thermal
Whitaker, S. 1992. “The species mass jump condition at a singular surface” Chemical
Whitaker S., 1999. The method of volume averaging, Kluwer academic publishers.
boundaries: the catalytic surface” Chemical Engineering Science, Vol. 55, 5231-
5245.
184
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Apéndice A-1:
dominio bidimensional, es
Problema I
∂ 2 sγ ∂ 2 sγ av ( y ) k
+ =− , en la fase fluida, (A-1.1)
∂x 2
∂y 2
Dγ ε γ ( y )
∂sγ k 1
C.F. 1 = , en x = a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.2)
∂x Dγ 2
∂sγ k 1
− = , en x = − a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.3)
∂x Dγ 2
∂sγ k
, en y = y2 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
= (A-1.4)
∂y Dγ 2 2
∂sγ k
, en y = y1 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
− = (A-1.5)
∂y Dγ 2 2
C.F.2 sγ = 0 , en y = 2h , (A-1.6)
∂sγ
C.F.3 = 0 , en y = 0 , (A-1.7)
∂y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2. (A-1.8)
se puede notar, se tienen suficientes condiciones de frontera en las dos direcciones, sin
A-1.1
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
numérico. Por ello, para analizar la posibilidad de modificar estas últimas condiciones
de frontera por otras de tipo Dirichlet o Neumann, se lleva a cabo un análisis de simetría
y = 2h
l
.
.
.
y=h
.
.
.
y2 ( k)
y1 (k)
.
.
.
y
x
x=a/2
Problema I
∂ 2 sγ ∂ 2 sγ av ( y ) k
+ =− , en la fase fluida, (A-1.9)
∂x 2
∂y 2
Dγ ε γ ( y )
∂sγ k 1
C.F. 1 = , en x = a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.10)
∂x Dγ 2
∂sγ k 1
− = , en x = − a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.11)
∂x Dγ 2
A-1.2
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂sγ k
, en y = y2 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
= (A-1.12)
∂y Dγ 2 2
∂sγ k
, en y = y1 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
− = (A-1.13)
∂y Dγ 2 2
C.F.2 sγ = 0 , en y = 2h , (A-1.14)
∂sγ
C.F.3 = 0 , en y = 0 , (A-1.15)
∂y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2. (A-1.16)
x → −x , (A-1.17)
y → y, (A-1.18)
se obtiene,
∂ 2 sγ ∂ 2 sγ av ( y ) k
+ =− , en la fase fluida, (A-1.19)
∂x 2
∂y 2
Dγ ε γ ( y )
∂sγ k 1
C.F. 1 − = , en x = − a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.20)
∂x Dγ 2
∂sγ k 1
= , en x = a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.21)
∂x Dγ 2
∂sγ k
, en y = y2 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
= (A-1.22)
∂y Dγ 2 2
∂sγ k
, en y = y1 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
− = (A-1.23)
∂y Dγ 2 2
C.F.2 sγ = 0 , en y = 2h , (A-1.24)
∂sγ
C.F.3 = 0 , en y = 0 , (A-1.25)
∂y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2. (A-1.26)
A-1.3
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Como se puede notar, para que el problema anterior sea idéntico al Problema I basta
sγ ( − x, y ) = sγ ( x, y ) . (A-1.27)
∂sγ
= 0, en x = 0 . (A-1.28)
∂x
l
Esta condición se puede extender para x = ± gracias a la condición de periodicidad.
2
∂ 2 sγ ∂ 2 sγ av ( y ) k
+ =− , en la fase fluida, (A-1.29)
∂x 2
∂y 2
Dγ ε γ ( y )
∂sγ k 1
C.F. 1 = , en x = a ; y1 ( k ) ≤ y ≤ y2 ( k ) , k = 1,..., nc , (A-1.30)
∂x Dγ 2
∂sγ k
, en y = y2 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
= (A-1.31)
∂y Dγ 2 2
∂sγ k
, en y = y1 ( k ) , − ≤ x ≤ , k = 1,..., nc ,
a a
− = (A-1.32)
∂y Dγ 2 2
C.F.2 sγ = 0 , en y = 2h , (A-1.33)
∂sγ
C.F.3 = 0 , en y = 0 , (A-1.34)
∂y
∂sγ
C.F.4 = 0, en x = 0 , (A-1.35)
∂x
∂sγ l
C.F.5 = 0, en x = . (A-1.36)
∂x 2
A-1.4
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
x y h a A 4a
X= , Y = , H = , A = , avω = γκ = 2 ,
l l l l V l
k k k ηω avω l 2
K1 = , K1 = ,φ =
2
. (A-1.37)
k ηω k ηω Dγ
Por otro lado, los cambios espaciales de la porosidad y el área interfacial serán descritos
por medio de las siguientes expresiones propuestas por Pérez Córdova y Ochoa-Tapia
(1995)
av ( y ) ⎡⎛ Y ⎞ ⎤
2
1⎛ Y ⎞
= 1+ ⎜ ⎟
4 ⎝ 20 A ⎠ ⎢⎜ 20 A ⎟ − 3⎥ = f1 (Y ) , 0 ≤ Y ≤ 40 A , (A-1.38)
avω ⎣⎝ ⎠ ⎦
⎛ Y ⎞ ⎡⎛ Y ⎞ ⎤
2
ε γ ( y ) = ε γω +
1
4
( ε γω − 1) ⎜ ⎟ ⎢⎜ ⎟ − 3⎥ = f 2 (Y ) , 0 ≤ Y ≤ 40 A . (A-1.39)
⎝ 20 A ⎠ ⎣⎝ 20 A ⎠ ⎦
Por lo que, el problema de cerradura se puede escribir acorde con la Figura A-1.2, como
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-1.40)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2 A
C.F.1 = K1 , en X = , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , (A-1.41)
∂X avω l 2
∂sγ φ2 A
= K1 , en Y = Y2 ( k ) , 0 ≤ X ≤ , (A-1.42)
∂Y avω l 2
∂sγ φ2 A
− = K1 , en Y = Y1 ( k ) , 0 ≤ X ≤ , (A-1.43)
∂Y avω l 2
A-1.5
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1
C.F.2 sγ = 0 , en Y = 2 H , 0 ≤ X ≤ , (A-1.44)
2
∂sγ 1
C.F.3 = 0 , en Y = 0 , 0 ≤ X ≤ , (A-1.45)
∂Y 2
∂sγ
C.F.4 = 0 , en X = 0, Y2 ( k − 1) < Y ≤ Y1 ( k ) , Y2 ( nc ) ≤ Y ≤ 2 H ; (A-1.46)
∂X
∂sγ 1
C.F.5 = 0 , en X = , 0 ≤ Y ≤ 2H . (A-1.47)
∂X 2
.
.
.
.
.
.
.
.
.
y
región.
A-1.6
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Apéndice A-2:
∂ 2 sγ ∂ 2 sγ f1 (YNX / 2 ) K1 2
+ =− φ , en la fase fluida, (A-2.1)
∂X 2
i, j
∂Y 2
i, j
f 2 (YNX / 2 )
∂sγ φ2
C.F.1 = K1 , en i = MX , OY ( k ) ≤ j ≤ MY ( k ) , k = 1,..., nc , (A-2.2)
∂X MX , j
avω l
∂sγ φ2
= K1 , en j = MY ( nc ) , 1 ≤ i ≤ MX , (A-2.3)
∂Y i, j
avω l
∂sγ φ2
− = K1 , en j = OY ( k ) , 1 ≤ i ≤ MX , (A-2.4)
∂Y i, j
avω l
C.F.2 sγ i , NY
= 0 , en j = NY , 1 ≤ i ≤ NX , (A-2.5)
∂sγ
C.F.3 = 0 , en j = 1 , 1 ≤ i ≤ NX , (A-2.6)
∂Y i ,1
∂sγ
C.F.4 = 0 , en i = 1, MY ( k − 1) ≤ j ≤ OY ( k ) , MY ( nc ) ≤ j ≤ NY ; (A-2.7)
∂X i, j
∂sγ
C.F.5 = 0 , en i = NX , 1 ≤ j ≤ NY , (A-2.8)
∂X i, j
A 1
X = 0 ≡ i =1; X = ≡ i = MX ; X = ≡ i = NX , (A-2.9)
2 2
1+ A
Y = 0 ≡ j = MY ( 0 ) = 1 ; Y2 ( k ) = ( k − 1) + ≡ j = MY ( k ) ; k = 1,.., nc ;
2
A-2.1
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1− A
Y1 ( k ) = ( k − 1) + ≡ j = OY ( k ) ; k = 1,..., nc ;
2
Y = H ≡ j = MY ( nc ) , Y = 2 H ≡ j = NY . (A-2.10)
Para resolver el problema se usará el método de diferencias finitas con inversión línea
f1 (YNX / 2 ) K1
sγ i −1, j
− 2 (1 + r ) sγ i, j
+ sγ i +1, j
=−
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i , j −1
+ sγ i , j +1 ), (A-2.11)
donde,
⎛ ΔX ⎞
2
r =⎜ ⎟ . (A-2.12)
⎝ ΔY ⎠
Ai sγ i −1, j
+ Bi sγ i, j
+ Ci sγ i +1, j
= Di , (A-2.13)
donde,
f1 (YNX / 2 ) K1
Ai = 1, Bi = −2 (1 + r ) , Ci = 1, Di = −
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i , j −1
+ sγ i , j +1 ) . (A-2.14)
Para discretizar las condiciones de frontera se usarán diferencias hacia delante o hacia
A-2.2
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
La dirección del barrido se es de arriba hacia abajo, como se indica en la Figura A-2.1;
sγ i , NY
= 0 , para 1 ≤ i ≤ NX . (A-2.15)
Posteriormente se barre toda la sección ocupada por el fluido que se encuentra por
A1 sγ 0, j
+ B1 sγ 1, j
+ C1 sγ 2, j
= D1 , (A-2.16)
donde,
A1 = 0, B1 = −1, C1 = 1, D1 = 0 . (A-2.17)
Ai sγ i −1, j
+ Bi sγ i, j
+ Ci sγ i +1, j
= Di , (A-2.18)
donde,
f1 (YNX / 2 ) K1
Ai = 1, Bi = −2 (1 + r ) , Ci = 1, Di = −
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i , j −1
+ sγ i , j +1 ) . (A-2.19)
Por último, para i = NX la ecuación es
ANX sγ NX −1, j
+ BNX sγ NX , j
+ CNX sγ NX +1, j
= DNX , (A-2.20)
donde,
Para resolver las matrices tridiagonales que se generan con este esquema, se utiliza el
algoritmo de Thomas (el cual está disponible en monografías como la de Fletcher, Cap.
6, 1991).
A-2.3
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Como se puede notar, en la Figura A-2.1 se eligió una celda con la que se pueden
calcular tantos cuadrados como se desee. La primera parte de la celda barre desde
Para 1 ≤ i ≤ MX , la ecuación es
Ai sγ i −1,OY ( k +1)
+ Bi sγ i ,OY ( k +1)
+ Ci sγ i +1,OY ( k +1)
= Di , (A-2.22)
donde,
φ2
Ai = 0, Bi = 1, Ci = 0, Di = sγ i ,OY (
− K1ΔY . (A-2.23)
)
k +1 −1 avω l
A-2.4
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Ai sγ i −1,OY ( k +1)
+ Bi sγ i ,OY ( k +1)
+ Ci sγ i +1,OY ( k +1)
= Di , (A-2.24)
donde,
Ai = 1, Bi = −2 (1 + r ) , Ci = 1,
Di = −
f1 (YNX / 2 ) K1
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i ,OY ( )
k +1 −1
+ sγ i ,OY ( k +1) +1 ). (A-2.25)
donde,
por j = MY ( k ) , de 1 ≤ i ≤ MX ,
Ai sγ i −1, MY ( k )
+ Bi sγ i , MY ( k )
+ Ci sγ i +1, MY ( k )
= Di , (A-2.28)
donde,
φ2
Ai = 0, Bi = 1, Ci = 0, Di = sγ i , MY (
− ΔYK1 . (A-2.29)
)
k +1 avω l
Ai sγ i −1, MY ( k )
+ Bi sγ i , MY ( k )
+ Ci sγ i +1, MY ( k )
= Di , (A-2.30)
donde,
Ai = 1, Bi = −2 (1 + r ) , Ci = 1,
Di = −
f1 (YNX / 2 ) K1
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i , MY ( k ) −1
+ sγ i , MY ( k ) +1 ). (A-2.31)
ANX sγ NX −1, MY ( k )
+ BNX sγ NX , MY ( k )
+ C NX sγ NX +1, MY ( k )
= DNX , (A-2.32)
A-2.5
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
donde,
Para i = MX la ecuación es
AMX sγ MX −1, j
+ BMX sγ MX , j
+ CMX sγ MX +1, j
= DMX , (A-2.34)
donde,
φ2
AMX = 0, BMX = −1, CMX = 1, DMX = K1ΔX . (A-2.35)
avω l
Ai sγ i −1, j
+ Bi sγ i, j
+ Ci sγ i +1, j
= Di , (A-2.36)
donde,
f1 (YNX / 2 ) K1
Ai = 1, Bi = −2 (1 + r ) , Ci = 1, Di = −
f 2 (YNX / 2 )
( ΔX ) φ 2 − r sγ
2
( i , j −1
+ sγ i , j +1 ) . (A-2.37)
Finalmente, para i = NX , la ecuación es
ANX sγ NX −1, j
+ BNX sγ NX , j
+ CNX sγ NX +1, j
= DNX , (A-2.38)
donde,
Es de notarse que para el cuadrado número nc, no se emplearán las ecuaciones (A-2.22)
sγ i ,1
= sγ i ,2
. (A-2.40)
Con esto se completa el barrido. El programa repetirá las operaciones mostradas en esta
sección hasta que se satisfaga una determinada tolerancia sobre el coeficiente efectivo,
definida como,
A-2.6
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
K eff − K effsupuesta
Tolerancia = , (A-2.41)
K effsupuesta
donde K effsupuesta se refiere al valor obtenido con el campo que de la iteración anterior y
para la primera iteración es el obtenido con el campo supuesto para iniciar el programa.
3/8 (ver por ejemplo, Nakamura 1992). A continuación (Figuras A-2.2-4) se muestran
diagrama de la subrutina con la que se resuelven las matrices tridiagonales ya que está
Inicio
Lee
datos
Fin
A-2.7
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Inicio
Fin
A-2.8
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
A-2.9
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
cuadrados.
Una vez alcanzada la convergencia del coeficiente efectivo con el esquema numérico
Respecto al número de puntos, lo que se hizo fue correr el programa con diferentes
números de puntos de manera tal que se usen celdas autocontenidas. Por ejemplo, al
usar 100, 200, 400, 800 y 1600 puntos para A = 0.2 , φ = 1 y once cuadrados, a partir de
la obtención de la función que ajuste a una gráfica del inverso de número de puntos
contra K eff / k ηω
se calcula la ordenada al origen (en este caso de 0.56), como se
muestra en la Figura A-2.5. Cabe mencionar que si sólo se hubieran tomado los últimos
A-2.10
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
0.55
0.50
K eff
0.45
k ηω
0.40
0.35
0.000 0.002 0.004 0.006 0.008 0.010
1/No. de puntos
Figura A-2.5: Coeficiente efectivo contra el inverso del número de puntos para A = 0.2 ,
φ = 1 y once cuadrados.
celda unitaria, para ello se resolvió el problema para A = 0.5 y φ = 3 usando desde uno
que ajuste adecuadamente (con un coeficiente de correlación superior a 0.95) los datos
como,
K eff 0.33298
= 1.0079 − 3.84764
, (A-2.42)
k ηω ⎛ ( No. de cuadrados )−1 ⎞
1+ ⎜ ⎟
⎜ 0.23883 ⎟
⎝ ⎠
calculó el porcentaje de error absoluto de cada uno de los resultados con respecto a este
valor, obteniendo los resultados presentados en la Tabla A-2.2, donde se puede observar
que con once cuadrados es suficiente para obtener un error inferior al 0.5%.
A-2.11
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.00
0.95
0.90
K eff 0.85
k ηω 0.80
0.75
0.70
0.65
0.04 0.1 1
1/No. de cuadrados
Figura A-2.6: Coeficiente efectivo vs. inverso del número de cuadrados para A = 0.5 y
φ = 3.
No. de K eff
− 0.67492
(100% )
cuadrados
k ηω
0.67492
1 48.16571
2 43.72074
3 37.79411
4 27.2536
5 16.33823
6 9.67522
7 6.04812
8 4.01233
9 2.53067
10 2.27434
11 0.20891
12 2.06099
13 0.42375
14 1.79577
15 0.36597
16 0.65637
17 0.60007
18 0.10223
19 0.19113
20 0.15261
siguiente.
A-2.12
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
A-2.13
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Apéndice A-3:
Al igual que en el caso del modelo de medio poroso centrado en las caras, es
medio poroso, con el fin de reemplazar las condiciones de periodicidad por condiciones
Problema Ih
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.1)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = − X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.2)
∂sγ φ2
= K1 , en X = − X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.3)
A-3.1
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂sγ φ2 1
= K1 , en Y = Y4 ( k ) , − X 2 ≤ X ≤ X 2 ; Y = Y2 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.4)
∂Y avω l 2
∂sγ φ2 1
− = K1 , en Y = Y3 ( k ) ; − X 2 ≤ X ≤ X 2 ; Y = Y5 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.5)
∂Y avω l 2
C.F.2 sγ = 0, en Y = 2 H , (A-3.6)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.7)
∂Y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.8)
Para llegar a la condición de frontera (A-3.7) se hicieron las mismas suposiciones que
en el caso del medio porosos centrado en las caras y, además, se tomó en cuenta que aún
∂sγ ,ω
= 0, en Y = 0 . (A-3.9)
∂Y
Problema Ih
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1
+ =− φ 2 , en la fase fluida, (A-3.10)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = − X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.11)
∂sγ φ2
= K1 , en X = − X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
A-3.2
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
X = X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.12)
∂sγ φ2 1
= K1 , en Y = Y4 ( k ) , − X 2 ≤ X ≤ X 2 ; Y = Y2 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.13)
∂Y avω l 2
∂sγ φ2 1
− = K1 , en Y = Y3 ( k ) ; − X 2 ≤ X ≤ X 2 ; Y = Y5 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.14)
∂Y avω l 2
Y=2H
.
.
.
Y=H
.
.
.
Y6 ( k)
Y5 ( k)
Y4 ( k)
Y3 ( k)
Y2 ( k)
. Y1 ( k)
.
.
X
X1
X2
C.F.2 sγ = 0, en Y = 2 H , (A-3.15)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.16)
∂Y
A-3.3
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.17)
X → −X , (A-3.18)
Y →Y . (A-3.19)
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.20)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 = K1 , en − X = X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.21)
∂sγ φ2
− = K1 , en X = X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
− X = X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.22)
∂sγ φ2 1
= K1 , en Y = Y4 ( k ) , − X 2 ≤ X ≤ X 2 ; Y = Y2 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.23)
∂Y avω l 2
∂sγ φ2 1
− = K1 , en Y = Y3 ( k ) ; − X 2 ≤ X ≤ X 2 ; Y = Y5 ( k ) , ± X 1 ≤ X ≤ ± , (A-3.24)
∂Y avω l 2
C.F.2 sγ = 0, en Y = 2 H , (A-3.25)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.26)
∂Y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.27)
Como se puede notar, para que el problema anterior sea idéntico al Problema Ih se debe
cumplir que,
sγ ( − X , Y ) = sγ ( X , Y ) . (A-3.28)
A-3.4
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1
las condiciones de periodicidad, este resultado se puede extender a X = ± . De esta
2
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.29)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 1 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ; (A-3.30)
∂X avω l
∂sγ φ2
= K1 , en X = X 2 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) ; (A-3.31)
∂X avω l
∂sγ φ2 1
= K1 , en Y = Y4 ( k ) , 0 ≤ X ≤ X 2 ; Y = Y2 ( k ) , X 1 ≤ X ≤ , (A-3.32)
∂Y avω l 2
∂sγ φ2 1
− = K1 , en Y = Y3 ( k ) ; 0 ≤ X ≤ X 2 ; Y = Y5 ( k ) , X 1 ≤ X ≤ , (A-3.33)
∂Y avω l 2
C.F.2 sγ = 0, en Y = 2 H , (A-3.34)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.35)
∂Y
∂sγ
C.F.4 = 0, en X = 0, 0 ≤ Y ≤ 2 H , (A-3.36)
∂X
∂sγ 1
C.F.5 = 0, en X = , 0 ≤ Y ≤ 2 H . (A-3.37)
∂X 2
Por otro lado, el esquema computacional a seguir para resolver este problema es similar
al caso anterior. Es decir, se inicia el barrido de arriba hacia abajo; una vez que se llega
A-3.5
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
computacional usado para este problema comparado con el caso del medio poroso
.
.
.
. Celda unitaria
.
. computacional
.
.
.
una tolerancia de 1x10−5 . Es decir, se usaron prácticamente la mitad de las celdas usadas
A-3.6
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Problema Iv
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.38)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = − X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.39)
∂sγ φ2
= K1 , en X = − X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.40)
∂sγ φ2
= K1 , en Y = Y4 ( k ) , − X 1 ≤ X ≤ X 1 ; Y = Y2 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.41)
∂Y avω l
∂sγ φ2
− = K1 , en Y = Y3 ( k ) ; − X 1 ≤ X ≤ X 1 ; Y = Y5 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.42)
∂Y avω l
C.F.2 sγ = 0, en Y = 2 H , (A-3.43)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.44)
∂Y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.45)
Problema Iv
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1
+ =− φ 2 , en la fase fluida, (A-3.46)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = − X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.47)
A-3.7
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂sγ φ2
= K1 , en X = − X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.48)
∂sγ φ2
= K1 , en Y = Y4 ( k ) , − X 1 ≤ X ≤ X 1 ; Y = Y2 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.49)
∂Y avω l
Y=2H
.
.
.
Y=H
.
.
.
Y6 (k)
Y4 ( k) Y5 (k)
Y2 (k)
Y3( k) Y1 (k)
.
.
.
Y
X
X1
X2
∂sγ φ2
− = K1 , en Y = Y3 ( k ) ; − X 1 ≤ X ≤ X 1 ; Y = Y5 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.50)
∂Y avω l
C.F.2 sγ = 0, en Y = 2 H , (A-3.51)
A-3.8
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.52)
∂Y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.53)
X → −X , (A-3.54)
Y →Y , (A-3.55)
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.56)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 = K1 , en − X = X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
X = X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.57)
∂sγ φ2
− = K1 , en X = X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ;
∂X avω l
− X = X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.58)
∂sγ φ2
= K1 , en Y = Y4 ( k ) , − X 1 ≤ X ≤ X 1 ; Y = Y2 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.59)
∂Y avω l
∂sγ φ2
− = K1 , en Y = Y3 ( k ) ; − X 1 ≤ X ≤ X 1 ; Y = Y5 ( k ) , ± X 2 ≤ X ≤ ±1 , (A-3.60)
∂Y avω l
C.F.2 sγ = 0, en Y = 2 H , (A-3.61)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.62)
∂Y
Periodicidad: sγ ( r + li ) = sγ ( r ) , i = 1, 2 . (A-3.63)
Al comparar el problema anterior con el Problema Iv, se puede notar que basta con
sγ ( − X , Y ) = sγ ( X , Y ) , (A-3.64)
A-3.9
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
para que ambos problemas sean idénticos. De lo anterior se concluir que sγ es simétrico
∂ 2 sγ ∂ 2 sγ f1 (Y ) K1 2
+ =− φ , en la fase fluida, (A-3.65)
∂X 2
∂Y 2
f 2 (Y )
∂sγ φ2
C.F.1 − = K1 , en X = X 2 , Y1 ( k ) ≤ Y ≤ Y2 ( k ) , Y5 ( k ) ≤ Y ≤ Y6 ( k ) ; (A-3.66)
∂X avω l
∂sγ φ2
= K1 , en X = X 1 , Y3 ( k ) ≤ Y ≤ Y4 ( k ) , (A-3.67)
∂X avω l
∂sγ φ2
= K1 , en Y = Y4 ( k ) , 0 ≤ X ≤ X 1 ; Y = Y2 ( k ) , X 2 ≤ X ≤ 1 , (A-3.68)
∂Y avω l
∂sγ φ2
− = K1 , en Y = Y3 ( k ) ; 0 ≤ X ≤ X 1 ; Y = Y5 ( k ) , X 2 ≤ X ≤ 1 , (A-3.69)
∂Y avω l
C.F.2 sγ = 0, en Y = 2 H , (A-3.70)
∂sγ
C.F.3 = 0, en Y = 0 , (A-3.71)
∂Y
∂sγ
C.F.4 = 0, en X = 0 , (A-3.72)
∂X
∂sγ
C.F.5 = 0, en X = 1 . (A-3.73)
∂X
radica en el tipo de celda unitaria computacional a usar, la cual es más complicada que
las anteriores. Nótese que, a diferencia de los casos anteriores, el ancho de la celda
unitaria computaciones es de la unidad. Por último, al igual que en el caso del medio
poroso centrado en las caras, al considerar once celdas y una tolerancia de 1x10−5 se
A-3.10
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
.
.
.
.
Celda unitaria
. computacional
.
.
.
.
j
A-3.11
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
Apéndice A-4:
1.0 M1
M2
M3
0.8
0.6
K eff
k eff ηω 0.4
d)
0.2
0.0
0.1 1
con los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el
A-4.1
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff ηω 0.4
d)
0.2
0.0
0.1 1 2
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.2
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff ηω 0.4
d)
0.2
0.0
0.1 1 3
φ
con los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el
A-4.3
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff
ηω 0.4
d)
0.2
0.0
0.1 1 4
φ
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.4
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff ηω 0.4
d)
0.2
0.0
0.1 1 6
φ
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.5
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff ηω 0.4
d)
0.2
0.0
0.1 1 8
φ
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.6
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff
ηω 0.4
d)
0.2
0.0
0.1 1 10 20
φ
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.7
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff
ηω 0.4
d)
0.2
0.0
0.1 1 10 20
φ
los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el cuerpo
A-4.8
Condiciones de frontera para el transporte Francisco J. Valdés-Parada
entre medios multifásicos
1.0 M1
M2
M3
0.8
0.6
K eff
k eff
ηω 0.4
d)
0.2
0.0
0.1 1 10 40
φ
con los modelos del medio poroso, a) centrado en las caras (M1), b) centrado en el
A-4.9