DINÁMICA DE FLUIDOS / MECÁNICA DE MEDIOS POROSOS

Ley de Darcy en Espresso: Resistencia Hidráulica, Permeabilidad ($k$) y Fracción $D_{10}$

Última revisión: 01/10/2026 · changelog: corrección de longitud de metadatos y estandarización de contenedores responsive
Esquema vectorial de permeabilidad de Darcy y distribución de finos D10 en la pastilla de espresso

La extracción de un espresso a 9 bar no es un proceso de disolución estática, sino un fenómeno de percolación forzada a través de un medio poroso compresible y soluble. El caudal de salida ($Q$) queda determinado de manera unívoca por la Ley de Darcy y la ecuación de Kozeny-Carman, donde la fracción granulométrica más fina ($D_{10}$) gobierna la resistencia hidráulica total del sistema.

1. Formulación Matemática del Caudal en Lecho Fijo

La ley de Darcy clásica para flujo unidireccional en medios porosos se expresa como:

$$Q = \frac{-k \cdot A}{\mu} \cdot \frac{\Delta P}{L}$$

Donde:

  • $Q$: Caudal volumétrico ($\text{m}^3/\text{s}$ o $\text{ml}/\text{s}$).
  • $k$: Permeabilidad intrínseca del lecho granular ($\text{m}^2$).
  • $A$: Área de la sección transversal de la cesta ($58.4\text{ mm} \rightarrow A \approx 2.68 \times 10^{-3}\text{ m}^2$).
  • $\mu$: Viscosidad dinámica del agua a $93\text{ }^\circ\text{C}$ ($\mu \approx 0.304\text{ mPa}\cdot\text{s} = 3.04 \times 10^{-4}\text{ Pa}\cdot\text{s}$).
  • $\Delta P$: Caída de presión a través de la pastilla ($\text{Pa}$, ej. $9\text{ bar} = 9 \times 10^5\text{ Pa}$).
  • $L$: Altura compactada del lecho de café ($\approx 10\text{ a }14\text{ mm}$).

2. Permeabilidad Intrínseca ($k$) y Ecuación de Kozeny-Carman

Para relacionar la permeabilidad intrínseca $k$ con la distribución de tamaño de partícula (PSD) y la porosidad intergranular ($\epsilon$), aplicamos la formulación de Kozeny-Carman:

$$k = \frac{\epsilon^3}{C_{KC} \cdot S_v^2 \cdot (1 - \epsilon)^2} = \frac{\epsilon^3 \cdot d_{\text{eff}}^2}{180 \cdot (1 - \epsilon)^2}$$

Donde $d_{\text{eff}}$ es el diámetro de partícula hidráulicamente efectivo. En distribuciones polidispersas de café molido, $d_{\text{eff}}$ no corresponde a la mediana ($D_{50} \approx 350\text{ }\mu\text{m}$), sino que está fuertemente sesgado hacia el percentil inferior $D_{10}$ ($< 100\text{ }\mu\text{m}$) debido a que el área superficial específica ($S_v$) es inversamente proporcional al diámetro de partícula ($S_v \propto 1/d$).

3. Matriz de Parámetros Hidráulicos

Rango Granulométrico $D_{10}$ Estimado ($\mu\text{m}$) Porosidad $\epsilon$ Permeabilidad $k$ ($\text{m}^2$) Caudal Resultante a 9 Bar ($\text{ml}/\text{s}$) Régimen Hidrodinámico
Molienda Gruesa (Turbo) $85 - 110$ $0.44$ $4.8 \times 10^{-13}$ $3.8 - 4.5$ Flujo laminar rápido
Molienda Tradicional $45 - 65$ $0.38$ $1.6 \times 10^{-13}$ $1.2 - 1.6$ Darcy estándar estable
Exceso de Finos (Muelas desalineadas) $< 30$ $0.31$ (migración) $0.4 \times 10^{-13}$ $< 0.4$ Colapso y canalización

4. Consecuencias Físicas en el Flujo de Extracción

Una reducción del $20\%$ en la fracción $D_{10}$ altera $d_{\text{eff}}^2$ en un factor de $0.64$, lo que provoca una caída del $36\%$ en la permeabilidad $k$. Si la bomba suministra un flujo volumétrico constante, la contrapresión hidrodinámica se disparará para compensar la resistencia, acelerando el colapso poroelástico del lecho.