IUSUR - Instituto Universitario Sudamericano

Hidrogeno Verde

Numeros adimensionales y distribución de corriente en

ingeniería electroquimica

CARLOS F. ZINOLA
LABORATORIO DE ELECTROQUÍMICA
FUNDAMENTAL, GRUPO INTERDISCIPLINARIO DE
INGENIERÍA ELECTROQUÍMICA, UNIVERSIDAD DE LA
REPÚBLICA, MONTEVIDEO, URUGUAY.
 Ingeniería Electroquímica
La ingeniería electroquímica es la rama de la ingeniería química que se
ocupa del diseño y desarrollo de las reacciones electroquímicas con
aplicaciones tecnológicas, como la electrosíntesis de productos químicos
de interés, la electrodeposición y el refinado de metales (galvanoplastia),
baterías convencionales y de flujo, celdas de combustible y
ultracapacitores, estudio de la morfología de la superficie por
electrodeposición, separaciones electrolíticas y corrosión e inhibición.

La ingeniería electroquímica combina el estudio de la

transferencia de carga en la interfase de los electrodos: conductor
electronico/conductor iónico (electrolitos) con el desarrollo de
materiales y los fenómenos asociados de transferencia de calor y
masa. Las consideraciones fundamentales incluyen materiales de
electrodos y la cinética de las reacciones electroquímicas y
especies redox.

El desarrollo de la tecnología implica el estudio de los reactores

electroquímicos, su distribución potencial y corriente, las
condiciones de transporte de masa, la hidrodinámica, la
geometría y los componentes, así como la cuantificación de su
rendimiento general en términos de rendimiento de reacción,
eficiencia de conversión y faradaica o energética.
 Reactores Electroquímicos para celdas de combustible
Celdas Acido Fosfórico

La primera monocelda de Grove en 1939 incorporaba ácido sulfúrico como

electrolito pero la corrosión disminuyó su desarrollo asi que comenzó a
estudiarse la posibilidad de electrolitos sólidos menos corrosivas.
En 1961 G.V. Elmore y H.A. Tanner dan los primeros resultados alentadores para
esta tecnología. Durante los años 60 la Armada de los Estados Unidos se
muestra muy interesada en el desarrollo de esta tecnología. Las celdas fueron
construidas por Allis- Chalmers (que también trabajó en AFC), Engelhard
Industries (se encargaba del reformador) y Varo, Inc. (desarrolló el inversor de
corriente).

Utilizan ácido fosfórico concentrado como electrolito, retenido en una

matriz de carburo de silicio, con una temperatura de operación de 150-
250°C. Los electrodos son de platino soportado sobre carbón, con una
alta sensibilidad al envenenamiento por monóxido de carbono. Este tipo
de celdas son las que tienen un mayor nivel de desarrollo tecnológico,
con utilización en aplicaciones estacionarias en fase de demostración.

Esta tecnología presenta una tolerancia a contaminantes generados en

el proceso de reformado de hidrocarburos muy superior al resto de pilas
de combustible de baja temperatura, por lo que su implantación será
más fácil a corto plazo.
 Reactores Electroquímicos para celdas de combustible

Celdas Alcalinas

Bacon junto a su modelo AFC. Su trabajo fue un éxito

cuando lo que se buscaba era una alternativa a los
electrolitos ácidos y los electrodos enteros de Platino.

El gran avance se conseguiría con el Proyecto Apollo, posterior al

Gemini (en el que se usaron pilas de combustible polimérica con
poco éxito). Pratt & Whitney, compañía fabricante de aviones
fundada en 1925, adquiere las patentes de Bacon se hace con el
contrato de la NASA para este proyecto. Esta tecnología también
fue usada en los programas espaciales soviéticos (SOYUZ) y, más
tarde, por la Agencia Rusa del Espacio en el Programa BURAN.
La NASA selecciona las pilas de combustible alcalinas porque
alcanzan eficiencias del 70%. Científicos de la NASA trabajando en
uno de los tres stacks de la nave Apollo (1964).
Reacción electroquímica y geometría de electrodos

Disco rotante

Esfera rotante
 Tratamiento para superficies de revolución
ejesimétricas. Teorema de Lighthill-Acrivos.

x
Disco normal

Eje de simetría  y
R(x)

Elemento diferencial electródico

Placa normal

Electrodo
Este capilar permite medir el potencial de
un punto cercano al electrodo (cerca de 3
veces el diámetro del tip)
PRODUCCIÓN De PERÓXIDO De HIDRÓGENO PRODUCCIÓN De HIDRÓGENO y OXIGENO

Reactores electroquímicos con electrodos de cilindro Celdas Alcalinas

rotatorio: tratamiento de efluentes (PRELINE)

(a) Esquema del reactor electroquímico; (b) Corte transversal del El electrolito es hidróxido de potasio concentrado al
cátodo; (c) Vista inferior del cátodo; (d) Representación esquemática de 35-50% a temperaturas menores a 120°C, pudiendo
un compartimiento anódico para anolito ácido. 1-Estructura 3D, 2- operar a 250°C cuando la concentración es del 85%.
soporte de grafito, 3-compartimiento anódico, 4-ánodo, 5-ventana del Los electrodos contienen a) metales nobles, ánodos
separador, 6-alimentador de O2, 7-cámara de O2, 8-canales para O2, 9-eje de platino/paladio y cátodos de oro/platino, b)
del electrodo, 10-cilindro de Teflón, 11-recipiente de vidrio, 12-camisa electrocatalizadores, como níquel, plata, óxidos
metálicos y espinelas de ellos.
de calefacción.
ARMADO DE MULTICELDAS DE COMBUSTIBLE de Hidrógeno/oxígeno.
ARMADO FINAL DE 44 UNIDADES 500 W Proyecto PDT 2009
 Reactores electroquímicos para recuperar metales o producción de gases con reciclo de líquidos

CAS es la concentración a la salida

CAE es la concentración a la entrada

v es el caudal o velocidad de avance
de un volumen del líquido por unidad
de tiempo
VTK es el volumen del tanque

REQ es el reactor electroquimico

iLim  nF k d C A (t )

iLim es la corriente límite de la reacción

considerando control por transferencia de
materia. n el número de electrones, F la
constante de Faraday
Reactores electroquímicos con reciclo de líquidos
Una de las aplicaciones típicas es la recuperación de metales a partir de soluciones que los
contienen en forma ionizada (Cobre a partir de una solución de sulfato de cobre).

El REQ (reactor electroquímico) es un canal rectangular de espesor h y largo L, con placas planas
paralelas.

A t=0 la concentración de la especie electroactiva A, es CAi y a t=tF es CAf.

El factor de conversión para el compuesto en el

reactor, x

Ae o ae es el área específica de los electrodos

Centrada  Csalida  (k d ae L / v)
Kd es el coeficiente de transferencia de materia x  1 e
 L
Centrada
1 D
K d   K d ( x)dx   i
L0  Kd =10-5 - 10-2 cm s-1 donde L es la distancia que separa ambos
electrodos (x de 0 a L)
Asi existe la posibilidad de definir un Leff (es decir la separación mínima que produzca v
Leff 
una conversión x importante) a partir de la siguiente expression: k d ae
Asi mismo, dado un valor de L que pueda fabricarse (ya que esto es acotado) lograr un
valor óptimo de área de electrode, ae eff efectiva: v
ae eff 
kd L

Área de electrodo necesaria en función del

Variación de CAE(t) a la entrada del REQ caudal en el REQ
 Lo que sí es necesario recalcar es que ae = as (1 – ) con as la actividad específica por unidad de volumen de
electrodo (sólido solamente) y el  es la porosidad de la superficie del soporte donde se coloca el catalizador de
actividad as. Es decir que tanto as como ae se miden en cm2/cm3. Aca es donde Uds. entran a decir ¿de donde voy a
sacar esos valores de a?. Bueno esos valores se estiman por medidas de microscopía.

Por ejemplo, para electrodos de lecho fluidizado (que son los más pobres) los as
andan en 15.000 mientras que ae en 10.000. Para el caso de granos (catalizador) fijos,
los as andan en 25.000 mientras que ae en 15.000.
Pero esto además cambia con el tiempo de operación continua.
Potencial AC Potencial AC
(1.50 a 0 V) (1.35 a 0.70 V)

Potencial DC
(-3.0 - 0.50 V)
Estudios electroquímicos en metanol
Una microscopía de
transmisión de la superficie del
electrodo muestra la
1.5 composición entramada del
platino depositado sobre
carbono.
1.0
-2
j / mA cm

0.5

0.0

0.0 0.4 0.8 0.0 0.4 0.8 0.0 0.4 0.8 1.2

E / V vs. rhe
Reactores “Swiss-roll cell” pilas con o
sin reciclo de liquidos
 Diseño de electrodos concéntricos estáticos y rotatorios
Flujo

[Link]í
rotaciones mica
Perfil de Couette, es el flujo tipo laminar de un fluido en el
espacio entre dos planos o esferas o.. en el que uno esta en
movimiento relativo con respecto al otro.
f .centrifugas 4 2 r 2
Ta  
f .vis cos as 2
Número de Taylor (Ta) es un número Papel
adimensional que caracteriza la
influencia de las fuerzas centrífugas Aluminio
(fuerzas inerciales debidas a la rotación
de un electrolito alrededor de un eje Aplicaciones en
vertical) respecto a las fuerzas viscosas. impresiones,
polímeros y primers
Componentes de las celdas PEM
 Placas bipolares de grafito

Una microscopía de
transmisión de la
superficie del electrodo
muestra la composición
entramada del platino
depositado sobre carbono.

Placas colectoras de corriente y difusoras de gases

Membrana polimérica y catalizadores

soportados en carbono (MEA)
Algunos trabajos técnicos desarrollados por nuestro Grupo para el
armado de Celdas de Combustible PEM
Placas bipolares de grafito
Tallado de las placas bipolares de grafito
 Esquema de estudio y diseño superficies homogéneas y columnares
Canal de Flujo
Flujo de
reactivos Capa de difusión de gases
homogénea 3Dai Sc1 / 6 w
Da  1/ 2
Capa catalítica de platino liso WaGz xy
de microesferas adyacentes
Soporte de carbono

Modelo de Superficie Columnar de Platino

0.8

0.6

coordenada vertical
0.4

0.2

d h

2 4 6 8 10 12
L
coordenada tangencial

3Dai Sc1 / 6 w  (T   sin T )

Da(T ) 
Canal de flujo WaGz1t ,/ohm
2
Flujo de
reactivos Capa de difusión de gases

Capa catalitica de platino columnar

Soporte de carbono
Curva de desempeño de una monocelda de hidrógeno/oxígeno

1.0 Celda Original hidrogeno

Patente No.43651 ‫ل‬nodo
Patente No.43651 c‫ل‬todo
0.8

E/V 0.6

0.4

0.2

0.0
0.0 0.2 0.4 0.6 0.8 1.0 1.2
I/A
ARMADO DE MULTICELDAS DE COMBUSTIBLE de Hidrógeno/aire.
ARMADO FINAL de 2500 W con AJUSA 2011.
Baterias o celdas multielectrodos. En un reactor electroquímico en general se
requieren varios electrodos para conseguir la
producción o consumo de sustancia que se
Electrodos en conexión (a) monopolar y (b) bipolar. desea alcanzar.

En la conexión en serie el potencial obtenido

es la suma de todos ellos y en algunos casos
no sería tan conveniente puesto que lo que se
desea es lograr altas corrientes (velocidades de
producción).

En la conexión en paralelo el potencial es

único y queda determinado por el del electrodo
que opere en mejores condiciones, y la
corriente igual sería la contribución de todos.

Aquí el efecto de caída óhmica es importante

asi como la posibilidad de “bypass” de algunos
electrodos que no funcionen adecuadamente en
Corriente las celdas bipolares.
bypass
 PROYECTO NACIONES UNIDAS GEF
Projecto EMTU Brasil Săo Paulo, Marcio Schetino.
Colaborador y Asesor: Dr. Fernando Zinola
PNUD (Programa das Nações Unidas para o Desenvolvimento) y
del GEF (Global Environmental Facility)

Garage de EMTU/SP
durante 4 años, para
obtener 1.000.000 km
rodados de experiencia
(julio de 2004)
Galvanizado

Distribución de corriente y
potencial

Cromado
 Definiciones y conceptos
La distribución de corriente y potencial es la formulación del perfil correspondiente en función del espacio
al cual se desarrolla y eventualmente en función del tiempo de progreso.
Los factores principales de los cuales depende esta distribución en un electrodo son;
Geometría del sistema en estudio, conductividad superficial de los electrodos y de la solución electrolítica,
sobrepotencial de activación (propiedades cinéticas de la reacción en cuestión) y de transporte de masa
(definida por la configuración para cada caso).
Dependiendo de las condiciones de operación del reactor y de su
estado de desarrollo (régimen estacionario o fuera del mismo), la
influencia de cada uno de los factores anteriores en la distribución de
corriente y potencial puede tener mas peso uno que otro.

Líneas de potencial

Líneas de
corriente
 Definiciones y conceptos
De acuerdo con ello se distinguen 3 tipos de distribución de corriente;

• distribución de corriente primaria, cuando solo depende de la geometría de

los electrodos y del reactor.

• distribución de corriente secundaria, cuando además de los factores

geométricos es necesario considerar la cinética de la reacción por medio
del sobrepotencial de activación. La homogeneidad de la distribución
establece un aumento del número adimensional llamado Número de Líneas isospotenciales
Wagner, Wa;
(d / dj ) 
Wa 
L
donde  es el sobrepotencial de activación definido como el apartamiento del
potencial de operación, E respecto del de circuito abierto, Ej=0 ocasionado por el
avance de la propia reacción;  es la conductividad de la solución cercana al
electrodo y L la longitud característica de la geometría del reactor de la celda
de combustible.

• distribución de corriente terciaria, cuando el proceso se encuentra

controlado por fenómenos de transporte del índole que sean por medio del Líneas de equiflujos -
sobrepotencial de concentración, o sea cuando estamos operando cerca corriente
de la corriente límite.
Introducción

En la revisión de la literatura se pueden ver los números importantes y la variedad de modelos de celdas
de combustible PEM. Una explicación que obstaculiza la comparación entre los modelos surge de la
ausencia de números adimensionales en la mayoría de las celdas PEM.
Gyenge publicó una revisión exhaustiva de todos los adimensionales que pueden encentrarse en las
celdas PEM, y enumera no menos de 25 números.
E.L. Gyenge, J. Power Sources 152 (2005) 105-121.
Xuan et al. centraron su estudio de celdas hidrógeno/aire usando el número de Damkhoeler, Da, en el
electrodo y el número de Graetz, Gz, en el canal. El impacto de estos números se traduce en la curva de
operación y de potencia.
J. Xuan, H. Wang, D.Y.C. Leung, M.K.H. Leung, H. Xu, L. Zhang, Y. Shen, J. Power Sources 231 (2013) 1-5.

(dm / dt )
Gz 
rDL donde el numerador se refiere al flujo másico, dm/dt, r es la densidad, D la difusividad, L la
longitud característica.
d d
Gz  Re Sc  Pe
L L donde d es el diámetro hidráulico, Pe el número de Peclet, Pe = ReSc, Re el número
de Reynolds y Sc el número de Schmidt.
En este caso veremos que la ecuación adimensional queda independiente de Gz o sea que actuaría
indistinta de una relación difusiva/convectiva usando en su lugar Pe.
 p

Da 
kC o L2 donde k es la cte de velocidad de la reacción, p es el orden de la reacción.
CoD
donde jo es la densidad de corriente de intercambio, f=F/RT siendo F la cte de
f ( E  E ) Faraday, R la cte de Regnault, E y Ej=0 son los potenciales operativos y de circuito
abierto, n el número de electrones intercambiados en la reacción,  el espesor
j0
jo e
Da 
nFDC o /  hidráulico o de difusión según corresponda y  el coeficiente de transferencia de
carga.

El uso de estos números presenta la ventaja de permitir la comparación de sus resultados con otros
estudios, independientemente de las condiciones de funcionamiento o de las propiedades de los
materiales de la pila de combustible.

La otra ventaja del análisis dimensional es definir regímenes operativos. Basado en el valor de algunos
parámetros adimensionales, es posible definir los procesos limitantes para un rango de condiciones de
funcionamiento (densidad de corriente, caudal de aire, difusión de oxígeno en la capa difusora, etc)
donde se pueden encontrar modelos simplificados. Los resultados presentados por Iranzo et al. ilustran
esta metodología.

A. Iranzo, M. Muñoz, F.J. Pino, F. Rosa, J. Power Sources 196 (2011) 4264-4269.
Desarrollo

Siguiendo el enfoque anterior, este estudio muestra un análisis adimensional con centro en un pseudo
modelo de celda de combustible 2D, el cual toma en cuenta el transporte de oxígeno en el canal y en
la capa difusora, el transporte de carga en la membrana y la reacción electroquímica en el cátodo.
Basado en este análisis dimensional, se define el régimen de funcionamiento de la celda de
combustible en los que se simplifica el modelo y se resuelven las ecuaciones analíticamente.

Distribuciones de densidad de corriente a lo largo de los canales se presentan para 4 rangos de

condiciones de funcionamiento definidas por los valores de tres números adimensionales solamente.

Para validar estos resultados teóricos se comparan los mismos utilizando datos experimentales para
distribuciones de densidad de corriente informadas anteriormente por nuestro Grupo de Trabajo y
mediciones de la presión parcial de oxígeno local a lo largo del canal, ya que es el reactivo limitante
por el alto sobrepotencial de activación.

Vamos a considerar el aspecto particular de la relación propia geométrica de una

celda de combustible, o sea de pequeño espesor y longitudes muy largas.
Asimismo, consideraremos el caso propio del transporte del oxígeno en el canal que
será solo a lo largo de la coordenada y, o sea 1 D pero que será a lo largo de la
coordenada x en la capa difusora de gases. Este modelo se conoce como pseudo
2D o 1D + 1D.
 Ecuaciones fundamentales para un canal estrecho y lineal.

Flujo de aire
canal
o
U aire

Capa difusora de gases

Capa catalítica
Soporte de celda

Transporte de Oxígeno en estado estacionario

Transporte de oxígeno en el canal

Donde CO2,canal es la concentracion molar de oxígeno en el
 CO2 ,canal ( y )  DO2
capa , dif
 CO2 , gas ( x, y )  canal, CO2,gas es la concentracion molar de oxígeno en la capa
U o
     difusora, Uoaire es la velocidad promedio del oxígeno en el aire
y x
aire
  hc   x 0 supuesta constante para un canal, hc es la altura del canal y
DO2capadif es la difusividad del oxígeno efectiva en la capa
Transporte de oxígeno en la capa difusora que es cte a lo largo de y.
difusora de gases

  2CO2 , gas ( x, y ) 
 0 DOcapa , dif
DOcapa , dif

 x 2  2
 2

1 1

1
  o
U aire hc  o
U aire hc Sc Re Pe
 Ley de Tafel modificada por transporte de oxígeno
CO2 , gas ( x  hd , y )  F ( E  E j  0 )  Donde E y Ej=0 son los potenciales en la capa catalítica
j( y)  j catodo
exp  para la operación de la celda y en circuito abierto,
CO2 ,canal ( y  0)
o
 RT  respectivamente, jocatodo es la densidad de corriente de
intercambio en el cátodo,  es el coeficiente de
O2 + 4H+ + 4e- ->2H2O transferencia de carga, F y R ctes universales.

La vinculación entre E y Ej=0 es la definición de sobrepotencial de activación para el propio

transporte de carga en la capa catalítica , (y), de donde surge el número Wa.

E  E j 0   ( y )    ( y )  (y), es el sobrepotencial óhmico que en general se expresa como

una relación lineal de la densidad de corriente por la ley de Ohm,
donde R es la resistencia óhmica siendo la suma de la resistencias
iónica y electrónica del soporte PEM y la capa difusora de gases,
  ( y )  j ( y ) R respectivamente.
 Ecuaciones fundamentales en números adimensionales
Aplicando el Teorema P de Buckingham obtenemos 3 números adimensionales:
dJ ( z ) J ( z)
 Pe 
dz 1  J ( z)   J ( z) 
1 1   exp 
Da  Wa   Wa 

Pe compara la razón entre el oxígeno transportado a lo largo

del canal y su difusión en la capa difusora de gases.

Da compara la cte cinética de la reacción catódica con la

difusión en la capa difusora de gases.

Wa compara la polarización por la transferencia de carga con

la resistencia interna de la celda PEM.
 Ecuación general en la superficie de los electrodos adimensional
La ecuación del reactor PEM es de carácter general y muestra la distribución de corriente a lo largo de
los canales simples en el Modelo 1D + 1D. Para el caso de al superficie de un electrodo, z=0 obtenemos:

Distribución de Corriente

Wa
Da
 Da 
1

Pe cte a   Da  e Wa
 Pe 

1 Potencial de Celda 0

Esta Figura representa la distribución de corriente vs. potencial de celda tomando la solución para
valores máximos iguales, J(0)=1, a 1 en z=0. Pe es cte e independiente del potencial.
Ecuación general de distribución de J(z)
corriente adimensional
J(z)
dJ ( z ) J ( z)
 Pe 
dz 1  J ( z)   J ( z)  Wa
1 1   exp 
Da  Wa   Wa 
J(z)
Wa
z
J(z)
z
Wa
J(z)

z
J(z) Wa

z
Wa
z Distribución de corriente adimensional
en función de la coordenada
z adimensional para diferentes Wa.
Estabilidad de funcionamiento a largo
plazo de las membranas Nafion
Esquema de 5 capas en un dispositivo típico de celda
polimérica de hidrógeno/oxígeno PEMFC.
Capa catódica
electrocatalítica

Entrada y salida de
oxígeno húmedo Capa difusora de
Placas gases catódica
bipolares

Placas
bipolares
Capa difusora de gases
anódica Entrada y salida de
hidrógeno húmedo

Capa anódica
electrocatalítica Membrana de Nafion de
electrolito polimérico
Diseño experimental para la evaluación de los perfiles de velocidad a lo largo de la capa porosa de Nafion. (a) flujímetros de
control; (b) tanques de agua acida; (c) Reservorio central; (d) Reservorio de recolección; (e) bomba; (f) filtro de agua; (g)
conector/inversor; (h) juntas planas de tope ; (i) soporte de plato; (j) n-microelectrodos de pH; (k) soportes de presión; (l)
medio poroso de Nafion; (m) electrode auxiliar; (n) válvulas; (o) potenciostato-galvanostato; (p) electrodos de referencia.
Micrografías SEM para la membrana de Nafion 117 virgen: (a) luego del método de purificación, (b) corte en seccion del
ensamble de membrana/catalizador, (c) envejecido por 2 años a 0.1 A cm-2, (d) envejecido por 2 años a 1 A cm-2.

10 micrones
(a) 0.2 micrones (b) 150 micrones

(c) 50 micrones (d) 5 micrones

Tabla Parámetros y propiedades geométricas, físicas y operacionales para la membrana
de Nafion 117 usada en la celda “home made” de hidrógeno/oxígeno PEMFC.

Cantidad y magnitud Dominio de valores/ unidades

Longitud del canal / altura / ancho L/h/d 15.0-30.0 cm/0.5-1.0 mm/0.5-3.0 mm

Espesor del Difusor de gases / capa catalítica / membrana 0.2-0.3/0.05/0.183 mm

de Nafion117 , 

Presiones anódica / catódica, P 1.0-2.0/1.0-2.0 atm

Porosidades de las muestras virgen y envejecida,  0.32/0.46

Permeabilidades para las muestras virgen y envejecida, k 3x10−12/3x10−9 m2

Relación entre el número de Re normalizado para la muestra porosa (usando caidas de presión con velocidad
tangencial) vs. el Re local viscoso real (velocidad tangencial).

Micrografías SEM para la

evolución de los electro-
catalizadores con el tiempo
de operción.

5000

4000
-  dP/dx ) k /   v x )

3000

2000

1000

0 1000 2000 3000 4000 5000 6000

(Lvx)/
CASO 1.- GRAN DISTORSIÓN EN LA MEMBRANA EN LOS CAMINOS DE
POROSIDAD (ENVEJECIMIENTO) O UNA GRAN TORTUOSIDAD INTRÍNSECA.

Perfil de velocidad adimensional laminar, Ux(x) para Re = 100 en distintas relaciones /Da con  = 0.32. Los valores
experimentales fueron obtenidos de una muestra preparada de Nafion 117 utilizada en una celda monopolar de
hidrogeno/oxigeno PEMFC durante 2 años de operación a 1 A cm-2.

Re Coordenada adimensional X.
0.8
Velocidad adimensional Ux = vx/Uo relacionada con el
valor límite, Uo;
0.6 P= Dp/rUo2 es la presión adimensional;
Re= UoL/ es el Número de Reynolds y
U X (X)

0.4 Da= kd2 es el Número de Darcy (relaciona la

/Da permeabilidad hidraulica y su camino cuadrático)
313
0.2
8300
125000

0.0

0.0 0.2 0.4 0.6 0.8 1.0

X-Re(dPdY)
 Caso 2.- Membrana de Nafion virgen de un medio poroso uniforme de porosidad media .

Perfil de velocidad adimensional para flujo laminar, Vx*(X) con Re* = 100 empleando diferentes valores de /Da con =0.32. Los
puntos marcados son valores experimentales obtenidos para el diseño mostrado anteriormente con la membrana virgen de
Nafion.

(2Re )(dPdY)
f

1.0

0.8

0.6
V X (X)
*

0.4
/Da

313
8300
0.2 125000

0.0

0.0 0.2 0.4 0.6 0.8 1.0

X
 Diagrama de Moody Es comun definir el factor de fricción f como la relación entre las diferencias de
presión y la pérdida de energía cinética para el fluido:

d DP 2dP
f  
L / 2 rv i 2 / 2 U x 2
Factor de fricción f para la membrana de Nafion añejada por 2
años usando /Da = 313 (puntos y líneas verdes) y para la
muestra virgen (puntos y lineas rojas) en función del Numero
6
de Gz calculado con Dh, L, Sc y Re.
5 El gráfico dentro muestra la curva doble logarítmica similar al
1
Diagrama de Moody. Comportamientos micro-turbulentos
4
Aged Nafion /Da=313
Prisitine Nafion
dentro de esos caminos en la membrana
3
f

2
0.1

1 100 1000

0
0 100 200 300 400 500 600 700

Gz
 Fluidodinámica de reactivos gaseosos en
placas bipolares rugosas y
electrocatalizadores de morfología diferente
Vortices en
Electrocatalizadores
Hemisféricos

Vortices de Foeppl
Vortices en
Electrocatalizadores
Hemisféricos
Vortices de Burger y Rankine
En Electrocatalizadores Columnares
Vortices de Burger y Rankine
En Electrocatalizadores Columnares
Preguntas!!???