El Coeficiente de Difusión: Cálculo y Medición

El coeficiente de difusión es una constante física fundamental que describe la velocidad a la que las partículas se dispersan de una región de mayor concentración a una de menor concentración. Este fenómeno, conocido como difusión, es omnipresente en la naturaleza y la ingeniería, afectando desde la forma en que el oxígeno llega a nuestras células hasta cómo se mezclan los contaminantes en la atmósfera. Comprender y poder calcular o medir este coeficiente es crucial en una vasta gama de disciplinas científicas y tecnológicas, incluyendo la química, la biología, la farmacia, la ingeniería de materiales y la ciencia ambiental. A lo largo de este artículo, exploraremos en profundidad qué es exactamente el coeficiente de difusión, cómo se calcula utilizando diversas fórmulas teóricas y empíricas, los métodos experimentales para su medición, los factores que influyen en su valor y sus múltiples aplicaciones prácticas.

¿Qué es el Coeficiente de Difusión (D)?

El coeficiente de difusión, denotado por la letra D, es el factor de proporcionalidad central en la Ley de Fick, una ley fundamental que rige el transporte de masa. Según la Ley de Fick, la masa de una sustancia (dM) que se difunde a través de una superficie (dF) normal a la dirección de difusión en un tiempo (dt) es directamente proporcional al gradiente de concentración (grad c) de dicha sustancia. Matemáticamente, esto se expresa como: dM = −D grad c dF dt. Físicamente, el coeficiente de difusión representa la masa de una sustancia que se difunde a través de una superficie unitaria en una unidad de tiempo, bajo un gradiente de concentración unitario. Su dimensión en el Sistema Internacional (SI) es metros cuadrados por segundo (m²/s).

Es importante destacar que el coeficiente de difusión no es una constante universal, sino una propiedad física que depende de múltiples factores. Entre ellos se incluyen el tamaño y otras propiedades de la molécula que se está difundiendo, así como las condiciones ambientales como la temperatura y la presión. Dada esta dependencia, los coeficientes de difusión de una sustancia en otra se determinan comúnmente de forma experimental y se recopilan en tablas de referencia, lo que facilita su uso en cálculos y diseños ingenieriles.

Cálculo del Coeficiente de Difusión: Fórmulas y Modelos

Aunque la determinación experimental es frecuente, en muchos casos, especialmente cuando los datos experimentales son escasos, el coeficiente de difusión puede ser calculado o estimado mediante diversas fórmulas teóricas y empíricas. Estas fórmulas varían significativamente dependiendo del estado de la materia del medio (gaseoso, líquido o sólido) y de la complejidad del sistema (binario o multicomponente).

Difusión en Gases

La difusión en gases es generalmente más rápida que en líquidos o sólidos debido a la mayor separación y movilidad de las moléculas. Para gases ideales, el coeficiente de difusión no depende de la concentración de la sustancia. La teoría cinética de los gases proporciona una base sólida para entender y modelar este fenómeno.

• Teoría Cinética de los Gases: Según esta teoría, la longitud de camino libre medio (l) de las moléculas es inversamente proporcional al área de la sección transversal media de la molécula (S) y a la densidad numérica de las moléculas (n) en una mezcla. Esta última es inversamente proporcional al espacio ocupado por la mezcla, es decir, T/p (temperatura/presión). La velocidad media de las moléculas (u) es proporcional a √(T/μ), donde μ es la masa molecular. Así, en el caso de la interdifusión de gases con la misma masa molecular (μ), o la autodifusión, la expresión es:

  D ∝ l * u ∝ (T/p) * √(T/μ)

• Difusión Binaria (Interacción Molecular): Si las moléculas de los tipos A y B interactúan, el coeficiente de interdifusión se expresa como:

  D_AB = Constante * (T^3/2 / (p * S_AB * √μ_AB))

  Donde S_AB es el valor medio de las secciones transversales de las moléculas de ambos tipos. La teoría cinética de los gases permite determinar la constante asumiendo moléculas esféricas y secciones transversales iguales a las de estas esferas.

• Corrección de Sutherland: Sutherland introdujo una corrección a la ecuación anterior, teniendo en cuenta las fuerzas de atracción intermoleculares que influyen en la longitud del camino libre de las moléculas:

  D_AB = D_AB,ideal * (1 + C_AB/T)^-1

  Donde C_AB es el coeficiente de Sutherland, que considera la atracción intermolecular.

• Potencial de Lennard-Jones: La teoría cinética contemporánea considera el carácter intrincado de la interacción molecular (repulsión a corta distancia y atracción a larga distancia). El potencial de Lennard-Jones es ampliamente reconocido para describir esta interacción:

  φ(r) = 4ε[(σ/r)¹² - (σ/r)⁶]

  Donde φ(r) es la energía potencial, r es la distancia entre los centros de las moléculas, y ε y σ son las constantes de interacción de Lennard-Jones, tabuladas para muchos gases.

• Ecuación de Bird, Hirshfelder y Curtiss: Para mezclas de gases diluidas, asumiendo colisiones moleculares binarias y elásticas, movimiento clásico y ausencia de efectos cuánticos, se obtiene la expresión:

  D_AB = (1.885 × 10^-2 * T^3/2 * √(1/M_A + 1/M_B)) / (P * σ_AB² * Ω)

  Aquí, Ω = f(kT/ε_AB) es el integral de colisión (un factor adimensional que depende de la energía cinética de las moléculas y la profundidad del pozo de potencial de Lennard-Jones), y k es la constante de Boltzmann. Los parámetros de interacción ε_AB y σ_AB se determinan para el sistema binario a partir de las constantes apropiadas para las sustancias puras:

  ε_AB = √(ε_Aε_B) y σ_AB = (σ_A + σ_B)/2

  Si los datos de ε y σ no están disponibles, pueden estimarse utilizando los parámetros críticos conocidos para una sustancia dada:

  ε/k ≈ 0.77 * T_c y σ ≈ 0.841 * V_c^1/3

  Donde T_c y V_c son la temperatura crítica (K) y el volumen molar crítico (cm³/mol), respectivamente.

• Corrección de Wilke y Lee: Wilke y Lee señalaron que el coeficiente 1.885 × 10^-2 no es constante en realidad y depende de las masas moleculares de los gases en difusión. Una definición más precisa resulta en un mejor acuerdo entre los datos experimentales y calculados.

  

• Fórmula de Correlación de Fuller, Schetter y Gittings: Esta fórmula empírica, obtenida mediante correlación asistida por ordenador de 340 puntos experimentales, es ampliamente conocida:

  D_AB = (1.00 × 10^-3 * T^1.75 * √(1/M_A + 1/M_B)) / (P * [(ΣV_A)^1/3 + (ΣV_B)^1/3]²)

  Donde (ΣV_A) y (ΣV_B) son los valores derivados de la suma de los volúmenes de difusión atómica para cada componente de la mezcla binaria (moléculas A y B). Los valores para algunos átomos y moléculas simples se presentan en la siguiente tabla:

  Tabla 1: Volúmenes de Difusión Atómica y Moléculas Simples

  Átomo/Molécula	Volumen de Difusión (cm³/mol)
  H	2.31
  He	2.88
  N	5.69
  O	6.78
  Cl	19.5
  S	17.0
  CO	18.0
  CO₂	26.9
  H₂O	13.1
  NH₃	14.9
  CH₄	25.8
  C₂H₆	39.7
  C₃H₈	54.0
  C₄H₁₀ (n-butano)	68.3
  C₄H₁₀ (i-butano)	67.0
  C₆H₆ (benceno)	73.0

• Sistemas Multicomponentes (Ecuación de Wilke): Para calcular el coeficiente de difusión en sistemas con múltiples componentes, Wilke utilizó la ecuación de Maxwell-Stefan para derivar la expresión:

  D'_A = (1 - Y_A) / (Σ_B≠A (Y_B / D_AB))

  Donde D'_A es el coeficiente de difusión del componente A en la mezcla, Y_A, Y_B, etc., son las fracciones molares de los componentes apropiados; y D_AB y D_AC son los coeficientes de difusión en los sistemas binarios AB y AC, respectivamente.

Difusión en Líquidos

Como se ha señalado, la difusión en líquidos encuentra una resistencia mucho mayor que en gases, resultando en coeficientes de difusión de 10⁴ a 10⁵ veces menores. Esto se debe a la mayor densidad de moléculas y a su menor movilidad en la fase líquida.

• Ecuación de Stokes-Einstein: Una de las ecuaciones más antiguas para determinar el coeficiente de difusión en soluciones diluidas se basa en el modelo de movimiento de una partícula esférica de la sustancia difusora A en un continuo líquido viscoso B:

  D_AB = kT / (b * π * η_B * r₀)

  Donde k es la constante de Boltzmann, T es la temperatura, r₀ es el radio de la partícula (molécula) y η_B es la viscosidad del líquido. La constante b depende del tamaño de las moléculas en difusión: b = 6 para moléculas más grandes que las de la sustancia base; b = 4 para moléculas idénticas; y b puede ser menor que 4 para moléculas más pequeñas. Asumiendo que el diámetro de la molécula es d = (6V_A / (πN₀))^1/3, donde V_A es el volumen molar de la sustancia difusora y N₀ es el número de Avogadro, la ecuación se convierte en:

  D_AB = k T / (b * π * η_B * (3V_A / (4πN₀))^1/3)

  La comparación de esta fórmula con datos experimentales ha demostrado que, en la mayoría de los casos, la discrepancia es moderada, alcanzando solo el 40% en algunos casos.

• Fórmula de Wilke y Chang: En 1955, Wilke y Chang propusieron una fórmula más general basada en extensas investigaciones experimentales, que involucra valores empíricos:

  D_AB = (7.4 × 10^-8 * (φ * M_B)^0.5 * T) / (η_B * V_A^0.6)

  Donde D_AB es el coeficiente de interdifusión en una solución infinitamente diluida (cm²/s); φ es el parámetro de asociación del disolvente B (refleja cómo las moléculas del disolvente se agrupan, afectando la difusión); M_B es la masa molecular del disolvente B; V_A es el volumen molar del soluto A en el punto de ebullición bajo condiciones normales (cm³/mol); η_B es la viscosidad del disolvente (Ns/m²); y T es la temperatura (K). La introducción del parámetro de asociación se debe a que las moléculas asociadas se comportan como moléculas de gran tamaño y difunden a una velocidad menor. Wilke y Chang presentaron valores para los disolventes más comunes: para agua φ = 2.6; metanol, 1.9; etanol, 1.5; benceno, éster, heptano y disolventes no asociados, 1.0.

• Fórmula Semiempírica de Scheibel: Una fórmula semiempírica sugerida por Scheibel es digna de atención:

  D_AB = K * T / (η_B * V_A^1/3)

  Donde K es una constante. En algunos casos, parece ser más exacta que la anterior; pero para V_A/V_B ≤ (1 - 2), la desviación del experimento se vuelve importante y se recomiendan las siguientes relaciones:

  Si V_A/V_B < 1, K = 8.2 × 10^-8

  Si V_A/V_B > 1, K = 10.0 × 10^-8

• Ecuación de Reddy y Doraiswamy: Reddy y Doraiswamy sugirieron la ecuación:

  D_AB = K_RS * T / (η_B * V_A^1/3)

  Donde K_RS varía dependiendo de la relación de volúmenes molares: K_RS = 8.5 × 10^-8 para V_A/V_B < 1 y K_RS = 10^-7 para V_A/V_B > 1. La comparación de esta ecuación con 96 puntos experimentales ha mostrado un buen acuerdo, con una dispersión de puntos de aproximadamente ±15%.

  Es crucial señalar que todas las fórmulas mencionadas para el cálculo del coeficiente de difusión en líquidos son válidas para líquidos de baja viscosidad. Para disolventes de alta viscosidad, estas fórmulas presentan un gran error y, por lo tanto, son inaplicables.

• Efecto de la Temperatura: El efecto de la temperatura en el coeficiente de difusión ha sido poco estudiado hasta ahora. Dentro de un rango de temperatura estrecho (de 10 a 20°C), la dependencia de la temperatura del coeficiente de difusión puede asumirse como lineal:

  D = D₂₉₈ * (η₂₉₈ * ρ) / (η * ρ₂₉₈)

  Donde D₂₉₈ es el coeficiente de difusión a T = 298K, η₂₉₈ es la viscosidad del disolvente a T = 298K, y ρ es la densidad del disolvente. Teóricamente, esta debe ser una dependencia exponencial del tipo D = AT exp (-E/RT), donde A es una constante, E es la energía de activación para la difusión y R es la constante de los gases ideales. Datos experimentales de Wilke y Chang sugieren que la energía de activación varía de 12.6 a 28.1 kJ/mol.

• Efecto de la Concentración: La dependencia del coeficiente de difusión de la concentración de la sustancia difusora es, estrictamente hablando, una consecuencia de que el flujo de difusión depende de la diferencia (gradiente) del potencial termodinámico del sistema, más que de la concentración. Esto implica que la fórmula debe tener en cuenta la actividad de la sustancia difusora. Así, a V_B = constante:

  D = D₀ * (1 + d ln a / d ln c)

  Donde D₀ y D son los coeficientes de difusión en una solución infinitamente diluida y en una solución con concentración finita c, respectivamente; a y c son la actividad y la concentración de la sustancia difusora; y n es el coeficiente de actividad de esta sustancia. Las fórmulas semiempíricas presentadas anteriormente son más exactas que las teóricas porque estas últimas se derivaron haciendo suposiciones. Sin embargo, para evitar un error apreciable, es aconsejable realizar cálculos con varias fórmulas simultáneamente y comparar los resultados.

Difusión en Soluciones Electrolíticas

En soluciones electrolíticas, las sales se disocian y difunden como iones y moléculas, dependiendo del grado de disociación. La teoría de la difusión de sales se ha desarrollado principalmente para soluciones diluidas, donde el grado de disociación es cercano a uno.

• Ecuación de Nernst-Heckell: Para una sal simple infinitamente diluida, el coeficiente de difusión puede encontrarse utilizando la ecuación de Nernst-Heckell:

  D_AB = (RT / Fa²) * ((n+ + n-) / (n+ * n-)) * (λ+ * λ-) / (λ+ + λ-)

  Donde D_AB es el coeficiente de difusión (cm²/s); T es la temperatura (K); Fa es el número de Faraday; n+ y n- son las valencias del catión y el anión; y λ+ y λ- son las conductividades iónicas límite (bajo dilución infinita) del catión y el anión (cm²/Ω mol).

• Fórmula de Gordon: En soluciones electrolíticas, el coeficiente de difusión depende sustancialmente de la concentración de la sustancia difusora. Si su concentración no supera los 2N, Gordon sugiere la fórmula:

  D = D₀ * (V / V₁) * (x₁ / (ν * η)) * (η₁ / η)

  Donde D₀ y D son los coeficientes de difusión para una solución infinitamente diluida y una solución molar, respectivamente; V es el volumen de la solución; V₁ es el volumen parcial del disolvente; x₁ es el número de moles de disolvente en el volumen V; ν es la molaridad de la solución; η y η₁ son la viscosidad de la solución y del disolvente; y v es el coeficiente de actividad molar.

Medición del Coeficiente de Difusión: Métodos Experimentales

Aunque las fórmulas de cálculo son útiles, la determinación experimental del coeficiente de difusión sigue siendo el método más preciso, especialmente para sistemas complejos o cuando se requiere una alta exactitud. Un método común y eficaz implica el uso de una celda de difusión de tres cámaras.

Este método es adecuado para determinar los coeficientes de difusión de solutos pequeños en volúmenes limitados de fluido (aproximadamente 4-9 ml). La difusión se mide en una celda de difusión de tres cámaras, donde los compartimentos están separados por membranas de nitrocelulosa y un compartimento central no agitado. Los compartimentos externos están agitados para mantener una concentración uniforme.

El proceso implica derivar el gradiente de concentración instantáneo y el flujo instantáneo del soluto hacia el compartimento final diluido a partir de los cambios en la concentración del soluto en los dos compartimentos finales agitados a lo largo del tiempo. El coeficiente de difusión se calcula a partir de la pendiente de la línea de regresión de mínimos cuadrados que relaciona la magnitud del flujo de soluto instantáneo con la del gradiente de concentración instantáneo. El aparato se calibra con un soluto de difusividad conocida, como el KCl. Los coeficientes de difusión determinados con este método en agua a 25°C para diversas sustancias (por ejemplo, CaCl₂: 7.54 × 10⁻⁶ cm²/s; Na₂-ATP: 7.01 × 10⁻⁶ cm²/s; 2-desoxiglucosa: 5.31 × 10⁻⁶ cm²/s; D-Na-lactato: 5.62 × 10⁻⁶ cm²/s) difirieron en un promedio de solo 3.7% de los valores de la literatura, lo que demuestra la precisión y la relativa rapidez de este procedimiento.

Factores que Afectan la Difusividad Másica

El coeficiente de difusión no es una constante inmutable para una sustancia dada, sino que está influenciado por varios factores clave. Comprender estos factores es esencial para predecir y controlar los procesos de difusión en diversas aplicaciones:

• Temperatura: Generalmente, la difusividad de una sustancia aumenta con la temperatura. Esto se debe a que las moléculas adquieren mayor energía cinética a temperaturas más altas, moviéndose más rápidamente y colisionando con mayor frecuencia y energía, lo que facilita su dispersión.

• Viscosidad del Medio: Un medio más viscoso (como el jarabe) impone una mayor resistencia al movimiento de las moléculas, reduciendo significativamente la difusividad en comparación con un medio menos viscoso (como el agua). La fricción interna es mayor, dificultando el desplazamiento molecular.

• Tamaño Molecular: Las moléculas más grandes suelen difundir más lentamente que las más pequeñas. Su mayor masa y tamaño molecular les otorgan una menor movilidad y requieren más energía para desplazarse a través del medio circundante.

• Gradiente de Concentración: Aunque el gradiente de concentración es la fuerza impulsora de la difusión (como se ve en la Ley de Fick), su magnitud también influye en la tasa. Cuanto mayor sea la diferencia de concentración entre dos regiones, mayor será la tasa de difusión inicial. Sin embargo, a medida que el gradiente disminuye (la mezcla se acerca a la uniformidad), la velocidad de difusión se ralentiza.

  

• Naturaleza de la Sustancia Difusora: Las propiedades químicas intrínsecas de la sustancia, como su polaridad, solubilidad o capacidad para formar puentes de hidrógeno, juegan un papel crucial en cómo interactúa con el medio y, por lo tanto, en su difusividad.

• Propiedades del Medio: Además de la viscosidad, otras propiedades del medio, como su porosidad (en materiales sólidos o geles) o su densidad, también pueden afectar la facilidad con la que una sustancia se difunde a través de él. Por ejemplo, los gases tienen una difusividad intrínsecamente más alta que los líquidos debido a las menores fuerzas intermoleculares en la fase gaseosa.

Aplicaciones del Coeficiente de Difusión

El conocimiento y la capacidad para manipular el coeficiente de difusión son de vital importancia en una amplia gama de campos científicos e industriales. A continuación, se detallan algunas de sus aplicaciones más destacadas:

• Ingeniería Química: La difusión es fundamental en numerosos procesos unitarios como la mezcla de fluidos, la separación de componentes (por ejemplo, en destilación o extracción), y la cinética de reacciones. En los reactores químicos, la velocidad de difusión puede ser el paso limitante que afecta la eficiencia de las reacciones, especialmente en procesos catalíticos donde los reactivos deben difundirse hacia la superficie del catalizador.

• Industria Farmacéutica: El coeficiente de difusión es crucial en el desarrollo y diseño de sistemas de liberación de fármacos. Las formulaciones de liberación controlada dependen de una comprensión precisa de cómo los medicamentos difunden a través de membranas biológicas o barreras poliméricas en el cuerpo, asegurando que el fármaco se libere a la velocidad y en el lugar adecuados.

• Sistemas Biológicos: En biología, la difusividad másica es esencial para comprender procesos vitales como el transporte de oxígeno y nutrientes a través de las membranas celulares y los tejidos, la eliminación de productos de desecho metabólico, y la difusión de moléculas señalizadoras que regulan las funciones celulares y orgánicas.

• Ciencias Ambientales: La difusión juega un papel crítico en la dispersión de contaminantes en el aire, el agua y el suelo. La capacidad de modelar cómo las sustancias se difunden en sistemas naturales ayuda a predecir los impactos ambientales de los derrames o emisiones, y a diseñar estrategias de remediación efectivas para mitigar la contaminación.

• Ciencia de Materiales: La difusividad es importante en procesos como la sinterización (la compactación de polvos mediante calor), la aplicación de recubrimientos protectores o funcionales, y la fabricación de materiales donde las sustancias se difunden en matrices sólidas para modificar sus propiedades (por ejemplo, dopado en semiconductores o endurecimiento superficial de metales).

Valores Típicos del Coeficiente de Difusión

Los valores del coeficiente de difusión (D) varían ampliamente, dependiendo de la sustancia que se difunde y del medio en el que lo hace. Generalmente, la difusión es mucho más rápida en gases que en líquidos, y extremadamente lenta en sólidos.

• En Agua: La difusividad de sustancias comunes en agua suele estar en el rango de 10^-9 a 10^-6 m²/s.
• En Aire: La difusividad de gases como el oxígeno o el dióxido de carbono en el aire es típicamente más alta, en el rango de 10^-5 a 10^-4 m²/s.
• En Sólidos: La difusividad en sólidos es generalmente mucho menor, en el rango de 10^-15 a 10^-10 m²/s, debido a la estructura rígida y la baja movilidad atómica/molecular.

A continuación, se presentan algunas tablas con valores típicos del coeficiente de difusión para diversas sustancias y medios:

Preguntas Frecuentes (FAQs) sobre el Coeficiente de Difusión

1. ¿Cuál es la diferencia entre difusividad másica y difusividad térmica?

La difusividad másica (o coeficiente de difusión) se refiere a la velocidad a la que las partículas o moléculas de una sustancia se dispersan dentro de otro medio, es decir, cómo se transfiere la masa. Por otro lado, la difusividad térmica describe la rapidez con la que el calor se propaga a través de una sustancia. Ambas son fenómenos de transporte, pero la difusividad másica se centra en la transferencia de masa, mientras que la térmica se enfoca en la transferencia de energía térmica.

2. ¿Cómo afecta el peso molecular de una sustancia a su difusividad?

Generalmente, las moléculas más pesadas difunden más lentamente que las más ligeras. Esto se debe a que su mayor tamaño y masa reducen su movilidad y aumentan la inercia, lo que dificulta su movimiento a través del medio.

3. ¿Puede el coeficiente de difusión ser constante en un sistema?

En muchos casos, especialmente bajo condiciones de estado estacionario o en soluciones diluidas y homogéneas, el coeficiente de difusión puede tratarse como una constante. Sin embargo, en sistemas no homogéneos, con grandes gradientes de concentración o temperatura, o en medios muy viscosos o complejos (como polímeros o soluciones electrolíticas concentradas), el coeficiente de difusión puede variar significativamente y no ser constante.

Si quieres conocer otros artículos parecidos a El Coeficiente de Difusión: Cálculo y Medición puedes visitar la categoría Física.