Articulo de referencia

La ley de Darcy

La ley de Darcy es una ecuación que describe el flujo de un fluido a través de un medio poroso y a través de una celda de Hele-Shaw . La ley fue formulada por Henry Darcy basánd...

La ley de Darcy es una ecuación que describe el flujo de un fluido a través de un medio poroso y a través de una celda de Hele-Shaw . La ley fue formulada por Henry Darcy basándose en los resultados de experimentos [ 1 ] sobre el flujo de agua a través de lechos de arena , formando la base de la hidrogeología , una rama de las ciencias de la tierra . Es análoga a la ley de Ohm en electrostática, relacionando linealmente el caudal volumétrico del fluido con la diferencia de carga hidráulica (que a menudo es simplemente proporcional a la diferencia de presión) a través de la conductividad hidráulica . De hecho, la ley de Darcy es un caso especial de la ecuación de Stokes para el flujo de momento , que a su vez deriva de la ecuación de Navier-Stokes para el momento .

La ley de Darcy es análoga a la ley de Fourier en el campo de la conducción del calor , a la ley de Ohm en el campo de las redes eléctricas y a la ley de Fick en la teoría de la difusión .

Historia

Enrique Darcy

Desde la década de 1830, el ingeniero hidráulico francés Henry Darcy estudió el caudal de agua en filtros de arena con el fin de mejorar el sistema de abastecimiento de agua municipal de la ciudad de Dijon , Francia. Estableció su ley empíricamente y la publicó en un informe de construcción en 1856. [ 2 ]

Descripción

Diagrama que muestra las definiciones y direcciones de la ley de Darcy. A es el área de la sección transversal (m² ) del cilindro. Q es el caudal (m³ / s) del fluido que fluye a través del área A. El flujo de fluido a través de A es q = Q/A. L es la longitud del cilindro. Δp = p salida - p entrada = p b - p a .pag{\displaystyle \nabla p}= Δp/L = gradiente hidráulico aplicado entre los puntos a y b.

En su forma integral, la ley de Darcy, en ausencia de fuerzas gravitatorias y en un medio homogéneamente permeable, viene dada por una sencilla relación de proporcionalidad entre el caudal volumétrico.Q{\displaystyle Q}y la caída de presiónΔpag{\displaystyle \Delta p} a través de un medio poroso . La constante de proporcionalidad está ligada a la permeabilidad.k{\displaystyle k}del medio, la viscosidad dinámica del fluidoμ{\displaystyle \mu }, la distancia dadaL{\displaystyle L}sobre la cual se calcula la caída de presión y el área de la sección transversalA{\displaystyle A}, en la forma: Q=kAμLΔpag{\displaystyle Q={\frac {kA}{\mu L}}\Delta p}

La proporción:

R=μLkA{\displaystyle R={\frac {\mu L}{kA}}}

También se puede definir como la resistencia hidráulica de la ley de Darcy .

La ley de Darcy puede generalizarse a una forma local:

Ecuación constitutiva de Darcy ( medios porosos isotrópicos )
q=kμpag{\displaystyle \mathbf {q} =-{\frac {k}{\mu }}\nabla p}

dóndepag{\displaystyle \nabla p}es el gradiente hidráulico yq{\displaystyle \mathbf {q} }es el flujo volumétrico que aquí también se denomina velocidad superficial . Nótese que la relación:

σ=kμ{\displaystyle \sigma ={\frac {k}{\mu }}}

puede considerarse como la ley de conductividad hidráulica de Darcy .

En la forma integral (menos general), el flujo volumétrico y el gradiente de presión corresponden a las siguientes relaciones:

q=QA{\displaystyle q={\frac {Q}{A}}}

pag=ΔpagL{\displaystyle \nabla p={\frac {\Delta p}{L}}}.

En el caso de un medio poroso anisotrópico, la permeabilidad es un tensor de segundo orden , y en notación tensorial se puede escribir la ley más general:

Ecuación constitutiva de Darcy ( medios porosos anisotrópicos )
qi=kijμjpag{\displaystyle q_{i}=-{\frac {k_{ij}}{\mu }}\partial _{j}p}

Observe que la cantidadq{\displaystyle \mathbf {q} }El flujo volumétrico, a menudo denominado flujo de Darcy o velocidad de Darcy, no es la velocidad a la que el fluido se desplaza a través de los poros. Es el flujo volumétrico, o caudal por unidad de área (véase Flux Flummoxed, Ground Water , 44 (2): 125–128, 2006). La velocidad de flujo ( u ) se relaciona con el flujo ( q ) mediante la porosidad ( φ ) con la siguiente ecuación:

q=φ.{\displaystyle \mathbf {q} =\varphi \,\mathbf {u} .}

La ecuación constitutiva de Darcy, para el flujo monofásico (de fluidos), es la ecuación que define la permeabilidad absoluta (permeabilidad monofásica).

En referencia al diagrama de la derecha, la velocidad del flujo está en unidades del SI.(metro/s){\displaystyle \mathrm {(m/s)} }y dado que la porosidad φ es un número adimensional , el flujo de Darcyq{\displaystyle \mathbf {q} }, o descarga por unidad de área, también se define en unidades(metro/s){\displaystyle \mathrm {(m/s)} }; la permeabilidadk{\displaystyle k}en unidades(metro2){\displaystyle \mathrm {(m^{2})} }, la viscosidad dinámicaμ{\displaystyle \mu }en unidades(PAGas){\displaystyle \mathrm {(Pa\cdot s)} }y el gradiente hidráulico se expresa en unidades.(PAGa/metro){\displaystyle \mathrm {(Pa/m)} }.

En forma integral, la caída de presión totalΔpag=pagbpaga{\displaystyle \Delta p=p_{b}-p_{a}}está en unidades(PAGa){\displaystyle \mathrm {(Pa)} }, yL{\displaystyle L}es la longitud de la muestra en unidades(metro){\displaystyle \mathrm {(m)} }, el caudal volumétrico de DarcyQ{\displaystyle Q}, o descarga, también se define en unidades(metro3/s){\displaystyle \mathrm {(m^{3}/s)} }y el área de sección transversalA{\displaystyle A}en unidades(metro2){\displaystyle \mathrm {(m^{2})} }Varios de estos parámetros se utilizan en definiciones alternativas a continuación. Se utiliza un signo negativo en la definición del flujo siguiendo la convención física estándar de que los fluidos fluyen de regiones de alta presión a regiones de baja presión. Nótese que debe tenerse en cuenta la altura de elevación si la entrada y la salida están a diferentes alturas. Si el cambio de presión es negativo, entonces el flujo será en la dirección positiva del eje x . Se han propuesto varias ecuaciones constitutivas para la permeabilidad absoluta, y la más famosa es probablemente la ecuación de Kozeny (también llamada ecuación de Kozeny-Carman ).

Considerando la relación para la presión estática del fluido ( ley de Stevin ):

pag=ρgramoh{\displaystyle p=\rho gh} También se puede declinar la forma integral en la ecuación: Q=kAgramoνLΔh{\displaystyle Q={\frac {kAg}{\nu L}}\,{\Delta h}} donde ν es la viscosidad cinemática . Por lo tanto, la conductividad hidráulica correspondiente es:

K=kρgramoμ=kgramoν.{\displaystyle K={\frac {k\rho g}{\mu }}={\frac {kg}{\nu }}.}

La ley de Darcy es una sencilla afirmación matemática que resume de forma concisa varias propiedades conocidas que presenta el agua subterránea que fluye en los acuíferos , entre ellas:

  • Si no hay gradiente de presión a lo largo de una distancia, no se produce flujo (estas son condiciones hidrostáticas ).
  • Si existe un gradiente de presión, el flujo se producirá desde la alta presión hacia la baja presión (en dirección opuesta al aumento del gradiente; de ​​ahí el signo negativo en la ley de Darcy).
  • cuanto mayor sea el gradiente de presión (a través del mismo material de formación), mayor será el caudal de descarga y
  • La tasa de descarga del fluido a menudo será diferente, ya sea a través de diferentes materiales de formación (o incluso a través del mismo material, en una dirección diferente), incluso si existe el mismo gradiente de presión en ambos casos.

Una ilustración gráfica del uso de la ecuación de flujo de agua subterránea en estado estacionario (basada en la ley de Darcy y la conservación de la masa) se encuentra en la construcción de redes de flujo , para cuantificar la cantidad de agua subterránea que fluye debajo de una presa .

La ley de Darcy solo es válida para flujos lentos y viscosos ; sin embargo, la mayoría de los casos de flujo de agua subterránea entran en esta categoría. Normalmente, cualquier flujo con un número de Reynolds menor que uno es claramente laminar, y sería válido aplicar la ley de Darcy. Las pruebas experimentales han demostrado que los regímenes de flujo con números de Reynolds de hasta 10 aún pueden ser darcianos, como en el caso del flujo de agua subterránea. El número de Reynolds (un parámetro adimensional) para el flujo en medios porosos se expresa típicamente como

Rmi=qdν,{\displaystyle \mathrm {Re} ={\frac {qd}{\nu }}\,,}

donde ν es la viscosidad cinemática del agua , q es el caudal específico (no la velocidad de los poros, con unidades de longitud por tiempo), d es un diámetro de grano representativo para el medio poroso (la opción estándar es math| d 30 , que es el tamaño que pasa al 30% de un análisis del tamaño de grano utilizando tamices, con unidades de longitud).

Derivación

Para un flujo estacionario, lento e incompresible , es decir , D ( ρu i ) / Dt ≈ 0 , la ecuación de Navier-Stokes se simplifica a la ecuación de Stokes , que al despreciar el término global es:

μ2iipag=0,{\displaystyle \mu \nabla ^{2}u_{i}-\partial _{i}p=0\,,}

donde μ es la viscosidad, u i es la velocidad en la dirección i y p es la presión. Suponiendo que la fuerza de resistencia viscosa es lineal con la velocidad, podemos escribir:

(k1)ijμφjipag=0,{\displaystyle -\left(k^{-1}\right)_{ij}\mu \varphi u_{j}-\partial _{i}p=0\,,}

donde φ es la porosidad y k ij es el tensor de permeabilidad de segundo orden. Esto da la velocidad en la dirección n ,

knortei(k1)ijj=δnortejj=norte=knorteiφμipag,{\displaystyle k_{ni}\left(k^{-1}\right)_{ij}u_{j}=\delta _{nj}u_{j}=u_{n}=-{\frac {k_{ni}}{\varphi \mu }}\partial _{i}p\,,}

lo que da la ley de Darcy para la densidad de flujo volumétrico en la dirección n ,

qnorte=knorteiμipag.{\displaystyle q_{n}=-{\frac {k_{ni}}{\mu }}\,\partial _{i}p\,.}

En medios porosos isotrópicos, los elementos fuera de la diagonal en el tensor de permeabilidad son cero, k ij = 0 para ij, y los elementos diagonales son idénticos, k ii = k , y la forma común se obtiene como se muestra a continuación, lo que permite determinar la velocidad del flujo del líquido resolviendo un conjunto de ecuaciones en una región dada. [ 3 ]

q=kμpag.{\displaystyle \mathbf {q} =-{\frac {k}{\mu }}\,{\boldsymbol {\nabla }}p\,.}

La ecuación anterior es una ecuación que rige el flujo de fluidos monofásicos en un medio poroso.

Aplicaciones

Ingeniería petrolera

Otra derivación de la ley de Darcy se utiliza ampliamente en la ingeniería petrolera para determinar el flujo a través de medios permeables; la más simple de las cuales es para una formación rocosa homogénea unidimensional con una sola fase fluida y viscosidad de fluido constante .

Casi todos los yacimientos de petróleo tienen una zona de agua debajo de la columna de petróleo, y algunos también tienen una capa de gas encima de la columna de petróleo. Cuando la presión del yacimiento cae debido a la producción de petróleo, el agua fluye hacia la zona de petróleo desde abajo, y el gas fluye hacia la zona de petróleo desde arriba (si existe la capa de gas), y se produce un flujo simultáneo y una mezcla inmiscible de todas las fases fluidas en la zona de petróleo. El operador del campo petrolífero también puede inyectar agua (o gas) para mejorar la producción de petróleo. Por lo tanto, la industria petrolera utiliza una ecuación de Darcy generalizada para el flujo multifásico desarrollada por Morris Muskat y otros . Debido a que el nombre de Darcy es tan conocido y está fuertemente asociado con el flujo en medios porosos, la ecuación multifásica se denomina ley de Darcy para el flujo multifásico , ecuación (o ley) de Darcy generalizada, o simplemente ecuación (o ley) de Darcy o ecuación de flujo si el contexto indica que el texto está tratando la ecuación multifásica de Muskat . El flujo multifásico en yacimientos de petróleo y gas es un tema extenso, y uno de los muchos artículos sobre este tema es la ley de Darcy para el flujo multifásico .

Preparación del café

Varios artículos han utilizado la ley de Darcy para modelar la física de la preparación del café en una cafetera moka , específicamente cómo el agua caliente se filtra a través del café molido bajo presión, comenzando con un artículo de 2001 de A. Varlamov y G. Balestrino, [ 4 ] y continuando con un artículo de 2007 de Concetto Gianino, [ 5 ] un artículo de 2008 de L. Navarini et al., [ 6 ] y un artículo de 2008 de W. King. [ 7 ] Los artículos tomarán la permeabilidad del café como constante como una simplificación o medirán el cambio a lo largo del proceso de preparación.

Formularios adicionales

Expresión diferencial

La ley de Darcy puede expresarse de forma muy general como:

q=Kh{\displaystyle \mathbf {q} =-K\nabla h}

donde q es el vector de flujo volumétrico del fluido en un punto específico del medio, h es la carga hidráulica total y K es el tensor de conductividad hidráulica en ese punto. La conductividad hidráulica a menudo se puede aproximar como un escalar . (Nótese la analogía con la ley de Ohm en electrostática. El vector de flujo es análogo a la densidad de corriente, la carga hidráulica es análoga al voltaje y la conductividad hidráulica es análoga a la conductividad eléctrica).

Ley cuadrática

Para flujos en medios porosos con números de Reynolds mayores que aproximadamente 1 a 10, los efectos inerciales también pueden volverse significativos. Generalmente, se agrega un término inercial a la ecuación de Darcy, conocido como término de Forchheimer . Este término permite tener en cuenta el comportamiento no lineal de la diferencia de presión en función del caudal. [ 8 ] [ 9 ]

pag=μkqρk1q2,{\displaystyle \nabla p=-{\frac {\mu }{k}}q-{\frac {\rho }{k_{1}}}q^{2}\,,}

donde el término adicional k 1 se conoce como permeabilidad inercial, en unidades de longitud(metro){\displaystyle \mathrm {(m)} }.

Cuando el número de Reynolds es mayor que 1, se observa el inicio del régimen de flujo inercial tras una transición desde el régimen de Darcy. La ley de Darcy deja de ser válida y se añade la corrección inercial no lineal. El origen de este comportamiento reside en el cambio de la forma de las líneas de corriente del flujo a escala de poro. Este cambio es inducido por las zonas de recirculación, que desempeñan un papel importante en la variación de la estructura del flujo. [ 9 ]

El flujo en el centro de un yacimiento de arenisca es tan lento que la ecuación de Forchheimer generalmente no es necesaria, pero el flujo de gas hacia un pozo de producción puede ser lo suficientemente alto como para justificar su uso. En este caso, los cálculos de rendimiento de entrada para el pozo, y no para la celda de la malla del modelo 3D, se basan en la ecuación de Forchheimer. Esto implica que aparece un factor de daño adicional dependiente del caudal en la fórmula de rendimiento de entrada.

Algunos yacimientos carbonatados presentan numerosas fracturas, por lo que la ecuación de Darcy para el flujo multifásico se generaliza para describir tanto el flujo en las fracturas como el flujo en la matriz (es decir, la roca porosa tradicional). La superficie irregular de las paredes de las fracturas y el elevado caudal en ellas podrían justificar el uso de la ecuación de Forchheimer.

Corrección para gases en medios finos (difusión de Knudsen o efecto Klinkenberg)

Para el flujo de gas en dimensiones características pequeñas (por ejemplo, arena muy fina, estructuras nanoporosas, etc.), las interacciones partícula-pared se vuelven más frecuentes, dando lugar a una fricción de pared adicional (fricción de Knudsen). Para un flujo en esta región, donde están presentes tanto la fricción viscosa como la de Knudsen , es necesario utilizar una nueva formulación. Martin Knudsen presentó un modelo semiempírico para el flujo en régimen de transición basado en sus experimentos en capilares pequeños. [ 10 ] [ 11 ] Para un medio poroso, la ecuación de Knudsen se puede expresar como [ 11 ]

norte=(kμpaga+pagb2+DKmiFF)1RTpagbpagaL,{\displaystyle N=-\left({\frac {k}{\mu }}{\frac {p_{a}+p_{b}}{2}}+D_{\mathrm {K} }^{\mathrm {eff} }\right){\frac {1}{RT}}{\frac {p_{\mathrm {b} }-p_{\mathrm {a} }}{L}}\,,}

donde N es el flujo molar, R es la constante de los gases , T es la temperatura, D eff K es la difusividad de Knudsen efectiva del medio poroso. El modelo también puede derivarse del modelo de fricción binaria basado en primeros principios (BFM). [ 12 ] [ 13 ] La ecuación diferencial del flujo de transición en medios porosos basada en BFM se da como [ 12 ]

pagincógnita=RT(kpagμ+DK)1norte.{\displaystyle {\frac {\partial p}{\partial x}}=-RT\left({\frac {kp}{\mu }}+D_{\mathrm {K} }\right)^{-1}N\,.}

Esta ecuación es válida tanto para capilares como para medios porosos. La terminología del efecto Knudsen y la difusividad de Knudsen es más común en ingeniería mecánica y química . En ingeniería geológica y petroquímica, este efecto se conoce como efecto Klinkenberg . Utilizando la definición de flujo molar, la ecuación anterior se puede reescribir como

pagincógnita=RT(kpagμ+DK)1pagRgramoTq.{\displaystyle {\frac {\partial p}{\partial x}}=-RT\left({\frac {kp}{\mu }}+D_{\mathrm {K} }\right)^{-1}{\dfrac {p}{R_{\mathrm {g} }T}}q\,.}

Esta ecuación se puede reorganizar en la siguiente ecuación.

q=kμ(1+DKμk1pag)pagincógnita.{\displaystyle q=-{\frac {k}{\mu }}\left(1+{\frac {D_{\mathrm {K} }\mu }{k}}{\frac {1}{p}}\right){\frac {\partial p}{\partial x}}\,.}

Comparando esta ecuación con la ley de Darcy convencional, se puede dar una nueva formulación como:

q=kmiFFμpagincógnita,{\displaystyle q=-{\frac {k^{\mathrm {eff} }}{\mu }}{\frac {\partial p}{\partial x}}\,,}

dónde

kmiFF=k(1+DKμk1pag).{\displaystyle k^{\mathrm {eff} }=k\left(1+{\frac {D_{\mathrm {K} }\mu }{k}}{\frac {1}{p}}\right)\,.}

Esto es equivalente a la formulación de permeabilidad efectiva propuesta por Klinkenberg: [ 14 ]

kmiFF=k(1+bpag).{\displaystyle k^{\mathrm {eff} }=k\left(1+{\frac {b}{p}}\right)\,.}

donde b se conoce como el parámetro de Klinkenberg, que depende de la estructura del gas y del medio poroso. Esto resulta evidente al comparar las formulaciones anteriores. El parámetro de Klinkenberg b depende de la permeabilidad, la difusividad de Knudsen y la viscosidad (es decir, de las propiedades tanto del gas como del medio poroso).

Ley de Darcy para escalas de tiempo cortas

Para escalas de tiempo muy cortas, se puede agregar una derivada temporal del flujo a la ley de Darcy, lo que da como resultado soluciones válidas en tiempos muy pequeños (en transferencia de calor, esto se llama la forma modificada de la ley de Fourier ).

τqt+q=kh,{\displaystyle \tau {\frac {\partial q}{\partial t}}+q=-k\nabla h\,,}

donde τ es una constante de tiempo muy pequeña que hace que esta ecuación se reduzca a la forma normal de la ley de Darcy en tiempos "normales" (> nanosegundos ). La razón principal para hacer esto es que la ecuación regular de flujo de agua subterránea ( ecuación de difusión ) produce singularidades en los límites de carga constante en tiempos muy pequeños. Esta forma es matemáticamente más rigurosa, pero conduce a una ecuación hiperbólica de flujo de agua subterránea, que es más difícil de resolver y solo es útil en tiempos muy pequeños, generalmente fuera del ámbito de la aplicación práctica.

La versión de Brinkman de la ley de Darcy

Otra extensión a la forma tradicional de la ley de Darcy es el término de Brinkman, que se utiliza para explicar el flujo de transición entre límites (introducido por Henri Brinkman en 1949 [ 15 ] ),

β2q+q=kμpag,{\displaystyle -\beta \nabla ^{2}q+q=-{\frac {k}{\mu }}\nabla p\,,}

donde β es un término de viscosidad efectiva . Este término de corrección tiene en cuenta el flujo a través de un medio cuyos granos son porosos, pero es difícil de usar y, por lo general, se omite.

flujo multifásico

En 1936, basándose en los resultados experimentales de sus colegas Ralph Dewey Wyckoff y Holbrook Gorham Boset, Morris Muskat y Milan W. Meres generalizaron la ley de Darcy para abarcar el flujo multifásico de agua, petróleo y gas en el medio poroso de un yacimiento petrolífero . [ 16 ] La ley de Darcy para el flujo multifásico, propuesta por Muskat y otros, constituye la base analítica de la ingeniería de yacimientos que perdura hasta nuestros días.

Validez

La ley de Darcy es válida para el flujo laminar a través de sedimentos . En sedimentos de grano fino, las dimensiones de los intersticios son pequeñas; por lo tanto, el flujo es laminar. Los sedimentos de grano grueso también se comportan de manera similar, pero en sedimentos de grano muy grueso, el flujo puede ser turbulento . [ 17 ] Por lo tanto, la ley de Darcy no siempre es válida en dichos sedimentos. Para el flujo a través de tuberías circulares comerciales, el flujo es laminar cuando el número de Reynolds es menor que 2000 y turbulento cuando es mayor que 4000, pero en algunos sedimentos, se ha encontrado que el flujo es laminar cuando el valor del número de Reynolds es menor que 1. [ 18 ]

Véase también

Referencias

  1. ^ Darcy, H. (1856). Les fontaines publiques de la ville de Dijon . París: Dalmont.
  2. Brown, GO (2002). "Henry Darcy y la creación de una ley" . Water Resources Research . 38 (7): 1106. Bibcode : 2002WRR....38.1106B . doi : 10.1029/2001WR000727 . ISSN 0043-1397 . 
  3. Adaptación de medios porosos para el flujo capilar controlable. Revista de Ciencia de Coloides e Interfaces 539 (2019) 379–387
  4. ^ A. Varlamov y G. Balestrino, "La fisica di un buon caffè", Il Nuovo Saggiatore 17°3-4°, 59–66°2001°.
  5. Gianino, Concetto. Análisis experimental de la cafetera italiana "moka". American Journal of Physics (2007)
  6. "Investigación experimental de la extracción de café a presión de vapor en una cafetera de estufa" L. Navarini, E. Nobile, F. Pinto, A. Scheri, F. Suggi-Liverani
  7. King, Warren. "La física de una máquina de café espresso de estufa". American Journal of Physics (2008)
  8. Bejan, A. (1984). Transferencia de calor por convección . John Wiley & Sons.
  9. 1 2 Agnaou, Mehrez; Lasseux, Didier; Ahmadi, Azita (2017). "Origen de la desviación inercial de la ley de Darcy: una investigación a partir de un análisis de flujo microscópico en estructuras modelo bidimensionales" . Physical Review E. 96 ( 4) 043105. American Physical Society. Bibcode : 2017PhRvE..96d3105A . doi : 10.1103/PhysRevE.96.043105 . hdl : 10985/12154 . PMID 29347623 . 
  10. Cunningham, RE; Williams, RJJ (1980). Difusión en gases y medios porosos . Nueva York: Plenum Press.
  11. 1 2 Carrigy, N.; Pant, LM; Mitra, SK; Secanell, M. (2013). "Difusividad de Knudsen y permeabilidad de capas de difusión de gas recubiertas microporosas de PEMFC para diferentes cargas de politetrafluoroetileno" . Journal of the Electrochemical Society . 160 (2): F81–89. doi : 10.1149/2.036302jes .
  12. 1 2 Pant, LM; Mitra, SK; Secanell, M. (2012). "Mediciones de permeabilidad absoluta y difusividad de Knudsen en capas de difusión de gas y capas microporosas de PEMFC". Journal of Power Sources . 206 : 153–160 . Bibcode : 2012JPS...206..153P . doi : 10.1016/j.jpowsour.2012.01.099 .
  13. Kerkhof, P. (1996). "Un modelo de Maxwell-Stefan modificado para el transporte a través de membranas inertes: El modelo de fricción binaria" . Chemical Engineering Journal y Biochemical Engineering Journal . 64 (3): 319– 343. doi : 10.1016/S0923-0467(96)03134-X .
  14. Klinkenberg, LJ (1941). "La permeabilidad de los medios porosos a líquidos y gases". Drilling and Production Practice . American Petroleum Institute. pp. 200–213 . 
  15. Brinkman, HC (1949). "Un cálculo de la fuerza viscosa ejercida por un fluido en movimiento sobre un enjambre denso de partículas". Applied Scientific Research . 1 (1) 27: 27– 34. Bibcode : 1949FTC.....1...27B . CiteSeerX 10.1.1.454.3769 . doi : 10.1007/BF02120313 . 
  16. III, Myron B. Allen (22 de junio de 2021). Las matemáticas del flujo de fluidos a través de medios porosos . John Wiley & Sons. ISBN 978-1-119-66384-3.
  17. Jin, Y.; Uth, M.-F.; Kuznetsov, AV; Herwig, H. (2 de febrero de 2015). "Investigación numérica de la posibilidad de turbulencia macroscópica en medios porosos: un estudio de simulación numérica directa". Journal of Fluid Mechanics . 766 : 76–103 . Bibcode : 2015JFM...766...76J . doi : 10.1017/jfm.2015.9 . S2CID 119946306 . 
  18. Arora, KR (1989). Mecánica de suelos e ingeniería de cimentaciones . Standard Publishers.
Obtenido de " https://en.wikipedia.org/w/index.php?title=Darcy%27s_law&oldid=1355639613 "