Método de Chapman-Enskog

Al finalizar el siglo XIX se conoce la ecuación de Boltzmann que rige la dinámica del medio gaseoso a la escala microscópica y las ecuaciones de Euler y de Navier-Stokes para el nivel macroscópico. Pasar de una escala a la otra constituye una parte del sexto problema de Hilbert. David Hilbert, autor de las declaraciones de los principales problemas considerados al finalizar el siglo XIX plantea las bases de un método bajo la forma de un desarrollo que lleva su nombre (1912). Hará falta esperar algunos años para que Sydney Chapman y David Enskog propogan simultáneamente e independientemente en 1916 y 1917 una solución a este problema.^[1]​^[2]​^[3]​ Más recientemente este método se ha extendido al caso de un gas en desequilibrio termodinámico,^[4]​ siendo este último aspecto un área de investigación muy activa en la actualidad.

El método de Chapman-Enskog es un método de perturbación que consiste en definir la solución bajo la forma de series de funciones de distribución en funciones de un « pequeño parámetro » asimilable al número de Knudsen. En orden cero se encuentra la distribución de Maxwell-Boltzmann y las ecuaciones de Euler. El orden uno permite conocer la expresión de los flujos de calor y de cantidad de movimiento y aquella de los coeficientes de transporte (los coeficientes de difusión por gradientes de concentración, de presión y de temperatura, las viscosidades dinámica, volumétricas, y la conductividad). De los potenciales de interacción molecular. Este enfoque permite encontrar las ecuaciones de Navier-Stokes y para justificar la difusión por gradientes térmico, desconocida en el tiempo en el que están publicados los trabajos de Chapman y de Enskog. Este método permitirá calcular todos estos coeficientes a partir del conocimiento de uno de ellos mediante la reconstitución a una medida (generalmente la viscosidad) de un potencial de interacción como el potencial de Lennard-Jones.

Harold Grad ha propuesto un enfoque alternativo que consiste en buscar la solución por los métodos de momentos de la función de distribución (1949). La ecuación de Boltzmann está multiplicada por ${\displaystyle 1,\mathbf {v} ,\mathbf {v} \otimes \mathbf {v} ,\mathbf {v} \otimes \mathbf {v} \otimes \mathbf {v} ...}$ (${\displaystyle \mathbf {v} }$ es la velocidad microscópica de la ecuación de Boltzmann y la ${\displaystyle \otimes }$el productor tensorial e integrado en velocidad. En este tipo de método de la ecuación en el n° de tiempo muestra la (n+1)°. Por lo tanto, debemos hacer una hipótesis para "cerrar" el sistema. Grad asume la solución expresada por una serie truncada de polinomios de Hermite. David Levermore ha propuesto más recientemente (1996) un cierre que usa una propiedad general: la solución maximiza la entropie del sistema de cerrábamos que son las partículas del medio.^[5]​ De los códigos de cálculo basado en estos métodos están quedado en la propiedad del laboratorio porque no aportando una ganancia notable en términos de propiedad de validez (en términos de número de Knudsen) por informe a los códigos estándares que resuelven las ecuaciones de Navier-Stokes, los cuales han hecho el objeto de desarrollos considerables.

Ecuaciones de evolución microscópica y macroscopia

Nivel microscópico

Apuntamos la función de distribución estadística de la velocidad. ${\displaystyle f_{i}(\mathbf {x} ,\mathbf {v} ,t)}$ . ${\displaystyle \mathbf {v} }$  en el momento ${\displaystyle t}$  en el punto ${\displaystyle \mathbf {x} }$  para la partícula (átomo o molécula) perteneciente a la especie ${\displaystyle i}$ . El número probable de partículas en el volumen ${\displaystyle [\mathbf {x} ,\mathbf {x} +\mathrm {d} \mathbf {x} ]}$  velocidades ${\displaystyle [\mathbf {v} ,\mathbf {v} +\mathrm {d} \mathbf {v} ]}$  en este momento es ${\displaystyle f_{i}\,\mathrm {d} \mathbf {x} \,\mathrm {d} \mathbf {v} }$ . La distribución estadística ${\displaystyle f_{i}}$ , por lo tanto se mide en s³⋅m^-6.

La ecuación de Boltzmann se escribe

${\displaystyle {\frac {\partial f_{i}}{\partial t}}+\mathbf {v} \cdot \nabla _{x}f_{i}=\sum _{j}Q(f_{i},f_{j})}$ 

dónde ${\displaystyle Q}$ , es el operador (o núcleo) de colisión, es un operador integral cuadrático que se describe a continuación, dando el efecto de las colisiones que se supondrán elásticas para simplificar el problema: no hay intercambio entre los grados de libertad interna y traslación, no hay reacción química. Por lo tanto se excluye el volumen de viscosidad que resulta de este tipo de intercambio.

Hay tantas distribuciones como de especies presentes en el medio. A cada una corresponde una ecuación de Boltzmann acoplada a las demás por los segundos miembros que representan colisiones homogéneas (${\displaystyle j=i}$ ) o heterogéneas (${\displaystyle j\neq i}$ ).

El choque elástico

 

Esquema de una interacción molecular elástica en el sistema relacionado al baricentro.

Las velocidades antes de la interacción en un Sistema de referencia inercial son ${\displaystyle \mathbf {v} _{i}}$  y ${\displaystyle \mathbf {v} _{j}}$ . Los índices representan indistintamente una misma especie o dos especies diferentes. Estas velocidades valen ${\displaystyle \mathbf {v} _{i}'}$  et ${\displaystyle \mathbf {v} _{j}'}$  después nos colocamos en un sistema centrado en el centro de gravedad que tiene una velocidad constante a causa de la conservación de la cantidad de movimiento. En este sistema que es por lo tanto galileo, la velocidad inicial de la partícula ${\displaystyle j}$  es la velocidad relativa ${\displaystyle \mathbf {g} _{ij}=\mathbf {v} _{i}-\mathbf {v} _{j}}$ . Por simetría se puede afirmar que la trayectoria estará contenida en el plano que contiene el origen y ${\displaystyle \mathbf {g} _{ij}}$ . Elegimos un punto de referencia como: ${\displaystyle \mathbf {\Omega } ={\frac {\mathbf {g} _{ij}}{||\mathbf {g} _{ij}||}}=[1,0,0]}$  (ver figura). En esta referencia la desviación es ${\displaystyle \theta _{ij}}$ . Según el parámetro de impacto ${\displaystyle b}$  de la velocidad relativa ${\displaystyle g_{ij}}$  y el potencial de interacción que se supone depende solo de la distancia en las dos partículas que interactúan. Si esta hipótesis es rigurosa para la interacción entre dos átomos, se puede considerarla utilizable para dos moléculas: el potencial es entonces un potencial medio estadístico.

La dirección de salida de interacción está definida por: ${\displaystyle \mathbf {\Omega '} ={\frac {\mathbf {g'} _{ij}}{||\mathbf {g'} _{ij}||}}}$ . Las velocidades finales se pueden calcular a partir de las siguientes consideraciones:

• la conservación de la cantidad de movimiento en la interacción implica

${\displaystyle \mathbf {v} _{i}'+\mathbf {v} _{j}'=\mathbf {v} _{i}+\mathbf {v} _{j}}$ 

• la velocidad relativa ${\displaystyle \mathbf {g} _{ij}}$  tiene un módulo constante debido a la conservación de la energía, por lo que ${\displaystyle g_{ij}'=g_{ij}}$  o

${\displaystyle |\mathbf {v} _{i}'-\mathbf {v} _{j}'|=|\mathbf {v} _{i}-\mathbf {v} _{j}|}$ 

Las velocidades después de la interacción son por lo tanto:

${\displaystyle \mathbf {v} _{i}'(b,g)=\mathbf {v} _{i}-(\mathbf {g_{ij}} \cdot \mathbf {\Omega '} )\,\mathbf {\Omega '} }$ 

${\displaystyle \mathbf {v} _{j}'(b,g)=\mathbf {v} _{j}+(\mathbf {g_{ij}} \cdot \mathbf {\Omega '} )\,\mathbf {\Omega '} }$ 

Además la conservación de momento cinético durante la interacción conduce a ${\displaystyle b'=b}$ .

El sistema que describe la colisión es reversible. El Teorema de Liouville (mecánica hamiltoniana) permite escribir

${\displaystyle \mathrm {d} \mathbf {v} _{i}'\mathrm {d} \mathbf {v} _{j}'=\mathrm {d} \mathbf {v} _{i}\,\mathrm {d} \mathbf {v} _{j}}$ 

La colisión del núcleo

El número probable de partículas que atraviesan el área ${\displaystyle 2\,\pi \,b\,\mathrm {d} b}$  por unidad de tiempo es ${\displaystyle f(\mathbf {v} _{j})\,g_{ij}\,2\,\pi \,b\,\mathrm {d} b\,\mathrm {d} \mathbf {v} _{j}}$ . Éllas interactúan con el número probable de partículas en el volumen elemental. ${\displaystyle f(\mathbf {v} _{i})\,\mathrm {d} \mathbf {r} \,\mathrm {d} \mathbf {v} _{i}}$  El número de partículas que desaparecen de la estadística por unidad de tiempo es ${\displaystyle \Theta _{ij}^{-}\,\mathrm {d} \mathbf {r} \,\mathrm {d} \mathbf {v} _{i}}$  con;

${\displaystyle \Theta _{ij}^{-}=2\pi \int _{\mathbf {v} }\int _{0}^{\infty }f(\mathbf {v} _{i})f(\mathbf {v} _{j})\,g_{ij}\,b\,\mathrm {d} b\,\mathrm {d} \mathbf {v} _{j}}$ 

De la misma manera, contamos la cantidad de partículas que aparecen

${\displaystyle \Theta _{ij}^{+}=2\pi \int _{\mathbf {v} }\int _{0}^{\infty }f(\mathbf {v} _{i}')f(\mathbf {v} _{j}')\,g_{ij}'\,b'\,\mathrm {d} b'\,\mathrm {d} \mathbf {v} _{j}'}$ 

Dadas las relaciones anteriores para la colisión. el operador de colisión se escribe

${\displaystyle Q(f_{i},f_{j})=\Theta _{ij}^{+}-\Theta _{ij}^{-}=2\pi \int _{\mathbf {v} }\int _{0}^{\infty }[f(\mathbf {v} _{i}')f(\mathbf {v} _{j}')-f(\mathbf {v} _{i})f(\mathbf {v} _{j})]\,g_{ij}\,b\,\mathrm {d} b\,\mathrm {d} \mathbf {v} _{j}}$ 

Esta ecuación está nombrada como la ecuación de Wang Chang y Uhlenbeck.

Se puede dar una formulación equivalente que introduce la sección transversal diferencial ${\displaystyle \sigma _{ij}}$  y se define:

${\displaystyle 2\,\pi \,b\,\mathrm {d} b=2\,\pi \,\sigma _{ij}\sin \theta _{ij}\,\mathrm {d} \theta _{ij}=\sigma _{ij}\,\mathrm {d} \Omega }$ 

por lo tanto

${\displaystyle Q(f_{i},f_{j})=\int _{\mathbf {v} }\int _{4\pi }[f(\mathbf {v} _{i}')f(\mathbf {v} _{j}')-f(\mathbf {v} _{i})f(\mathbf {v} _{j})]\,g_{ij}\,\sigma _{ij}\,\mathrm {d} \mathbf {\Omega } }$ 

Nivel macroscópico

Las variables

La ecuación de Boltzmann describe la evolución de partículas a nivel macroscópico. Para describir cada una de las especies ${\displaystyle i}$  a nivel macroscópico esto se define:

- la densidad de partículas	${\displaystyle n_{i}=\int _{\mathbf {v} }f_{i}\mathrm {d} \mathbf {v} _{i}}$	${\displaystyle n=\sum _{i}n_{i}}$
- la densidad	${\displaystyle \rho _{i}=n_{i}m_{i}}$	${\displaystyle \rho =\sum _{i}\rho _{i}}$
- la velocidad media	${\displaystyle \mathbf {V} _{i}={\frac {1}{n_{i}}}\int _{\mathbf {v} }\mathbf {v} _{i}f_{i}\mathrm {d} \mathbf {v} _{i}}$	- la energía interna	${\displaystyle e_{i}=\int _{\mathbf {v} }{\frac {1}{2}}m_{i}v_{i}^{2}f_{i}\mathrm {d} \mathbf {v} _{i}}$

Podemos entonces definir valores para el conjunto de las especies

- la velocidad total	${\displaystyle \mathbf {V} ={\frac {1}{\rho }}\sum _{i}\rho _{i}\mathbf {V} _{i}}$	(velocidad baricéntrica de todas las partículas)
- la velocidad de difusión	${\displaystyle {\mathbf {V} _{D}}_{i}=\mathbf {V} _{i}-\mathbf {V} }$	(cantidad cinématique)

Algunas variables auxiliares (${\displaystyle {\mathcal {N}}}$  es el número de Avogadro)

- la fracción de volumen	${\displaystyle x_{i}={\frac {n_{i}}{n}}}$	telle que ${\displaystyle \sum _{i}x_{i}=1}$
- la fracción masa	${\displaystyle c_{i}={\frac {\rho _{i}}{\rho }}}$	telle que ${\displaystyle \sum _{i}c_{i}=1}$
- la masa molar	${\displaystyle {\mathcal {M}}_{i}={\mathcal {N}}m_{i}}$	de valeur moyenne ${\displaystyle {\overline {\mathcal {M}}}=\sum _{i}x_{i}{\mathcal {M}}_{i}}$

Los flujos

El flujo de cantidad ${\displaystyle \psi _{i}}$  es por definición la cantidad ${\displaystyle \int _{\mathbf {v} }\psi _{i}f_{i}\mathbf {U} _{i}\mathrm {d} \mathbf {U} _{i}}$  donde ${\displaystyle \mathbf {U} _{i}=\mathbf {v} _{i}-\mathbf {V} }$ 

esto se define observando el producto producto diádico ${\displaystyle \otimes }$ 

- el flujo de masa	${\displaystyle \mathbf {J} _{i}=m_{i}\int _{\mathbf {v} }f_{i}\mathbf {U} _{i}\mathrm {d} \mathbf {U} _{i}}$	${\displaystyle =\rho _{i}{\mathbf {V} _{D}}_{i}}$
- el tensor de presión	${\displaystyle {\mathsf {P}}_{i}=\int _{\mathbf {v} }m_{i}\mathbf {U} _{i}\otimes \mathbf {U} _{i}f_{i}d\mathbf {U} _{i}}$	que representa el flujo de cantidad de movimiento. Es simétrico por construcción.
- la presión parcial	${\displaystyle p_{i}={\frac {1}{3}}\int _{\mathbf {v} }m_{i}U_{i}^{2}f_{i}d\mathbf {U} _{i}}$	se define a partir del rastro del tensor de presión.
- el flujo de calor	${\displaystyle \mathbf {q} _{i}=\int _{\mathbf {v} }{\frac {1}{2}}m_{i}U_{i}^{2}\mathbf {U} _{i}f_{i}\mathrm {d} \mathbf {v} }$	representa el flujo de energía interna.

Los flujos globales se obtienen simplemente sumando ${\displaystyle i}$  así como la presión:${\displaystyle p=\sum _{i}p_{i}}$ . Podemos entonces definir una temperatura a partir de la ecuación de estado

${\displaystyle p=n\,k\,T}$ 

Ecuaciones de evolución

Multiplicando cada una de las ecuaciones de Boltzann sucesivamente por cada uno de los invariantes de colisión ${\displaystyle m_{i},m_{i}\mathbf {v} ,{\frac {1}{2}}m_{i}v^{2}}$  y al integrarse en las velocidades y, si es necesario, en la especie, se obtienen las ecuaciones de evolución macroscópicas llamadas ecuaciones del producto contraído.

${\displaystyle {\frac {\partial \rho _{i}}{\partial t}}+\nabla \cdot (\rho _{i}\mathbf {V} +\mathbf {J} _{i})=0}$ 

${\displaystyle \rho \,{\frac {\partial \mathbf {V} }{\partial t}}+\rho \,(\mathbf {V} \cdot \nabla )\mathbf {V} +\nabla \cdot {\mathsf {P}}={\frac {\partial (\rho \,\mathbf {V} )}{\partial t}}+\nabla \cdot (\rho \mathbf {V} \otimes \mathbf {V} )+\nabla \cdot {\mathsf {P}}=0}$ 

${\displaystyle {\frac {\partial (\rho \,e)}{\partial t}}+\nabla \cdot (\rho \,e\,\mathbf {V} )+\nabla \cdot \mathbf {q} +\nabla \cdot ({\mathsf {P}}{\textbf {:}}\mathbf {V} )=0}$ 

Todos los segundos miembros son nulos a causa de las leyes de conservación

${\displaystyle \int _{\mathbf {v} }\psi _{i}\,Q\,\mathrm {d} \mathbf {v} =0,~~\forall i,~~\forall \psi \in [m_{i},m_{i}\mathbf {v} ,{\frac {1}{2}}m_{i}v^{2}]}$ 

Así obtenemos un sistema de evolución ${\displaystyle \rho }$ , ${\displaystyle \mathbf {V} }$  y ${\displaystyle e}$  en el que los flujos ${\displaystyle \mathbf {J} _{i}}$ , ${\displaystyle {\mathsf {P}}}$  y ${\displaystyle \mathbf {q} }$  quedan por aclarar.

Ecuación de Boltzmann adimensional

Se supone que es un medio homogéneo (una sola especie presente).

Con el fin de estimar la aportación de cada término en la ecuación de Boltzmann es necesario modificar este. Para eso se define las cantidades de referencia siguiente:

• la temperatura ${\displaystyle T^{*}}$  de la cual deducimos una velocidad de referencia que es la velocidad del sonido para una gas perfecto ${\displaystyle a^{*}={\sqrt {\gamma RT^{*}}}}$ 
• la longitud de referencia macroscópica ${\displaystyle L^{*}}$  Relacionado con el problema que se define un tiempo de referencia ${\displaystyle t^{*}={\frac {L^{*}}{a^{*}}}}$  Relacionados con la propagación de ondas acústicas
• presión de referencia ${\displaystyle p^{*}}$ , con la que se define una densidad de referencia ${\displaystyle \rho ^{*}={\frac {p^{*}}{RT^{*}}}}$  y la cantidad ${\displaystyle f^{*}={\frac {\rho ^{*}}{{a^{*}}^{3}}}}$ 
• una velocidad de referencia microscópica que corresponde a la velocidad mediana de una partícula ${\displaystyle c^{*}={\sqrt {\frac {8kT^{*}}{m\pi }}}}$ 
• otra velocidad de referencia microscópica que corresponde a la velocidad relativa promedia entre dos partículas ${\displaystyle g^{*}=c^{*}{\sqrt {2}}}$ 
• una longitud de referencia microscópica correspondiente a la trayectoria libre promedia ${\displaystyle l^{*}={\frac {1}{{\sqrt {2}}\sigma n}}}$  donde ${\displaystyle \sigma }$  es la sección transversal diferencial.

Si ahora definimos las variables reducidas ${\displaystyle {\tilde {f}}={\frac {f}{f^{*}}}}$ , ${\displaystyle {\tilde {t}}={\frac {t}{t^{*}}}}$ , ${\displaystyle {\tilde {v}}={\frac {v}{c^{*}}}}$ , ${\displaystyle {\tilde {g}}={\frac {g}{g^{*}}}}$ , ${\displaystyle {\tilde {b}}={\frac {b}{\sqrt {\frac {\sigma }{\pi }}}}}$  y ${\displaystyle {\tilde {x}}={\frac {x}{L^{*}}}}$ , la ecuación de Boltzmann esta escrita

${\displaystyle {\mathcal {S}}r{\frac {\partial {\tilde {f}}}{\partial {\tilde {t}}}}+{\tilde {\mathbf {v} }}\cdot \nabla _{\tilde {x}}{\tilde {f}}={\frac {1}{Kn}}{\tilde {Q}}({\tilde {f}},{\tilde {f}})}$ 

donde:

• ${\displaystyle {\tilde {Q}}({\tilde {f}},{\tilde {f}})={\frac {\tau ^{*}}{f^{*}}}Q(f,f)}$ 

• número de Strouhal ${\displaystyle {\mathcal {S}}r={\frac {L^{*}}{c^{*}t^{*}}}}$  que pesa la variación temporal del sistema,
• número de Knudsen ${\displaystyle Kn={\frac {l^{*}}{L^{*}}}}$  cuyo peso inverso es el término de colisión, por lo que la tendencia a volver al equilibrio. En la práctica, para que las ecuaciones del continuo sean válidas, es necesario que ${\displaystyle Kn<0.01}$ , típicamente.

Formulación de la solución

Escribimos la solución como una serie usando un parámetro ${\displaystyle \Lambda }$  del mismo orden de magnitud que el número de Knidsen

${\displaystyle f_{i}=f_{i}^{(0)}+\Lambda \,f_{i}^{(1)}+\Lambda ^{2}\,f_{i}^{(2)}+...}$ 

${\displaystyle f_{i}}$  Respetando las leyes de conservación, cada uno de los términos de desarrollo también deben de respetarlas. De donde las restricciones sobre la solución.

${\displaystyle \int _{\mathbf {v} }\psi \,f_{i}^{(k)}\,\mathrm {d} \mathbf {v} =0,~~\forall i,~~\forall k,~~\forall \psi \in [m,m\mathbf {v} ,{\frac {1}{2}}mv^{2}]}$ 

usamos esta aproximación en la ecuación de Boltzmann y se separa los términos correspondientes a cada potencia de ${\displaystyle \Lambda }$ .

Al orden cero

Se obtiene simplemente

${\displaystyle \sum _{j}Q(f_{i}^{(0)},f_{j}^{(0)})=0}$ 

Esta ecuación está verificada si todos los términos que la componen son nulos, así que en particular

${\displaystyle Q(f_{i}^{(0)},f_{i}^{(0)})=0}$ 

Lo que implica

${\displaystyle f_{i}^{(0)}(\mathbf {v'} )f_{i}^{(0)}(\mathbf {w'} )=f_{i}^{(0)}(\mathbf {v} )f_{i}^{(0)}(\mathbf {w} )}$ 

o

${\displaystyle f_{i}^{(0)}(\mathbf {v'} )+f_{j}^{(0)}(\mathbf {w'} )=\log f_{i}^{(0)}(\mathbf {v} )+\log f_{j}^{(0)}(\mathbf {w} )}$ 

${\displaystyle \log f}$  Es un variante colisional. Por lo tanto se escribe como una combinación lineal de las invariantes colisionales canónicos

${\displaystyle \log \,f_{i}^{(0)}=A_{i}\,m_{i}+\mathbf {B} _{i}\cdot (m_{i}\mathbf {v} )+C_{i}\,{\frac {1}{2}}m_{i}v^{2}}$ 

Introduciendo esta expresión en las ecuaciones que definen las variables microscópicas se identifica los parámetros de este desarrollo y se encuentra la ley de distribución de las velocidades de Maxwell

${\displaystyle f_{i}^{(0)}=n_{i}\left({\frac {\beta _{i}}{\pi }}\right)^{\frac {3}{2}}e^{-\beta _{i}||\mathbf {V} _{i}-\mathbf {v} ||^{2}}}$ 

con. ${\displaystyle \beta _{i}={\frac {m_{i}}{2kT}}}$ . Los flujos de difusión ${\displaystyle \mathbf {J} _{i}}$  se anulan, así como el flujo de calor ${\displaystyle \mathbf {q} }$ . El tensor de presión está reducido a su rastro ${\displaystyle {\mathsf {P}}=p\,{\mathsf {I}}}$  donde ${\displaystyle {\mathsf {I}}}$ . Las ecuaciones macróscopicas correspondientes son las ecuaciones de Euler.

Al orden uno

El orden uno mustra una ecuación integral de Fredholm para la incógnita ${\displaystyle f_{i}^{(1)}}$ 

${\displaystyle {\frac {\partial f_{i}^{(0)}}{\partial t}}+\mathbf {v} \cdot \nabla f_{i}^{(0)}=\sum _{j}\left[Q(f_{i}^{(0)},f_{j}^{(1)})+Q(f_{i}^{(1)},f_{j}^{(0)})\right]}$ 

La difícil resolución de esta ecuación permite dar las diversas cantidades desconocidas de las ecuaciones de Enskog que se puede entonces asimilar a las ecuaciones de Navier-Stokes.

El flujo de difusión

Se obtiene bajo forma de un sistema lineal llamado sistema de Stefan-Maxwell

${\displaystyle {\begin{array}{rcl}\sum _{j\neq i}{\frac {x_{i}x_{j}}{\rho {\mathcal {D}}_{ij}}}\left({\frac {\mathbf {J} _{j}}{c_{j}}}-{\frac {\mathbf {J} _{i}}{c_{i}}}\right)&=&\nabla x_{i}+(x_{i}-c_{i})\nabla \log p+k_{i}^{T}\nabla \log T\\[0.6em]&=&\mathbf {d} _{i}+k_{i}^{T}\nabla \log T\end{array}}}$ 

Donde vemos la aparición coeficiente de difusión binaria ${\displaystyle {\mathcal {D}}_{ij}(T,p)}$  y el "coeficiente de difusión térmica multicomponente",${\displaystyle k_{i}^{T}(T,p,x_{i})}$  (en realidad un número adimensional) ${\displaystyle {\mathcal {D}}_{i}^{T}(T,p,x_{i})}$  (que no son coeficientes de difusión y que pueden ser negativos) por:

${\displaystyle k_{i}^{T}=\sum _{j\neq i}{\frac {x_{i}x_{j}}{\rho {\mathcal {D}}_{ij}}}\left({\frac {{\mathcal {D}}_{j}^{T}}{c_{j}}}-{\frac {{\mathcal {D}}_{i}^{T}}{c_{i}}}\right)}$ 

Para un medio con ${\displaystyle N}$  especies el rango de este sistema es ${\displaystyle N-1}$  ya que ${\displaystyle \sum _{i}\mathbf {J} _{i}=0}$ . Su solución

${\displaystyle \mathbf {J} _{i}=-{\frac {\rho }{{\overline {\mathcal {M}}}^{2}}}\sum _{j\neq i}{\mathcal {M}}_{i}{\mathcal {M}}_{j}D_{ij}[\nabla x_{i}+(x_{i}-c_{i})\nabla \log p]-{\mathcal {D}}_{i}^{T}\nabla \log T}$ 

Encontramos los términos de difusión por gradiente de concentración, de presión y de temperatura (efecto Soret).${\displaystyle D_{ij}(T,p,x_{i})}$  Es el coeficiente de difusión multicomponente, solución de un sistema lineal que involucra los coeficientes binarios. Este sistema también es de rango términos ${\displaystyle N-1}$  La solución no es única e implica ${\displaystyle {\frac {N(N-1)}{2}}}$  términos independientes. Generalmente elegimos ${\displaystyle D_{ii}=0}$  por el bien de la simetría. Esta elección es arbitraria.

Hay diversas soluciones aproximadas del sistema de Stefan-Maxwell que permite obtener una expresión explícita del flujo de difusión bajo una forma vecina de la ley de Fick, la cual no es exacta para una mezcla binari.

El tensor de presión

El tensor de presión tiene una forma clásica

${\displaystyle {\mathsf {P}}=p{\mathsf {I}}-2\mu (T,p,x_{i}){\mathsf {S}}}$ 

donde ${\displaystyle {\mathsf {I}}}$  es el tensor unitario ${\displaystyle {\mathsf {S}}}$  El tensor de estrés viscoso

${\displaystyle {\mathsf {S}}={\frac {1}{2}}[\nabla \otimes \mathbf {V} +(\nabla \otimes \mathbf {V} )^{t}]-{\frac {1}{3}}(\nabla \cdot \mathbf {V} )\,{\mathsf {I}}}$ 

Un término adicional ${\displaystyle -\eta (\nabla \cdot \mathbf {V} )\,{\mathsf {I}}}$  aparece en la presión cuando se tienen en cuenta las interacciones inelásticas. Su influencia es débil incluso totalmente despreciable para los gases poco densos.^[4]​

Se observa que la hipótesis de Stokes se justifica naturalmente por este enfoque.

El flujo de calor

Está dado por

${\displaystyle \mathbf {q} =-\lambda \nabla T+\sum _{i}\rho _{i}e_{i}{\mathbf {V} _{D}}_{i}+p\sum _{i}k_{i}^{T}{\mathbf {V} _{D}}_{i}}$ 

${\displaystyle \lambda (T,p,x_{i})}$  Esta es la conductividad térmica. El último término de la ecuación es el corolario del efecto Soret y está nombrado efecto Dufour.

Los coeficientes de transporte

Los coeficientes de transporte se expresan en forma de sistemas lineales que involucran cantidades del tipo ${\displaystyle \int _{\mathbf {V} }F(V_{i},f_{i}^{(0)})d\mathbf {V} }$ , las cuales se desarrollan en polinomios de Sonine-Laguerre. Los coeficientes del desarrollo se expresan en funciones integrales de colisión. En lo práctico estamos satisfechos con el primer orden para el desarrollo y las integrales de colisión son funciones de la temperatura tabuladas por varios autores. Además hay soluciones aproximadas de los sistemas lineales que dan los diversos coeficientes bajo forma explícita.

La función de distribución

La función de distribución es

${\displaystyle f_{i}^{(1)}=f_{i}^{(0)}\left[1+{\frac {1}{x_{i}}}\mathbf {V} _{i}\cdot \mathbf {d} _{i}+\mathbf {b} _{i}:(\mathbf {\nabla } \otimes \mathbf {V} )+\left(W_{i}^{2}-{\frac {5}{2}}\right)\mathbf {V} _{i}\cdot \nabla \log T\right]}$ 

donde ${\displaystyle \mathbf {d} _{i}}$  se da en el flujo de difusión

${\displaystyle \mathbf {W} _{i}={\sqrt {\frac {m_{i}}{2kT}}}\mathbf {V} _{i}}$ 

${\displaystyle \mathbf {b} _{i}=2\left(\mathbf {W} _{i}\otimes \mathbf {W} _{i}-W_{i}^{2}\,{\mathsf {I}}\right)}$ 

Esta función de distribución es necesaria para el cálculo de la capa de Knudsen que da las condiciones a la pared para la ecuación de Navier-Stokes.

Al orden dos

Aquí nuevamente obtenemos una integral de Fredholm para la incógnita ${\displaystyle f_{i}^{(2)}}$ 

${\displaystyle {\frac {\partial f_{i}^{(1)}}{\partial t}}+\mathbf {v} \cdot \nabla f_{i}^{(1)}=\sum _{j}\left[Q(f_{i}^{(0)},f_{j}^{(2)})+Q(f_{i}^{(1)},f_{j}^{(1)})+Q(f_{i}^{(2)},f_{j}^{(0)})\right]}$ 

David Burnett propuso en 1935 una solución de esta ecuación. Esta tiene el inconveniente de no respetar el théorème^[6]​ H. Parece que la subida en orden constituye un callejón sin salida, todas las variantes propuestas hasta este día no resuelven este problema.^[7]​

Notas y referencias

Notas

Referencias

1. (en) Sydney Chapman et Thomas George Cowling, The Mathematical Theory of Non-uniform Gases, Cambridge University Press, 1991 (ISBN 0-521-40844-X)
2. (en) Joseph Oakland Hirschfelder, Charles Francis Curtiss et Robert Byron Bird, Molecular Theory of Gases and Liquids, John Wiley and Sons, 1966 (ISBN 978-0-471-40065-3)
3. (en) Gilberto Medeiros Kremer, « The Methods of Chapman-Enskog and Grad and Applications », RTO-EN-AVT 194, 2011 [1]
4. Raymond Brun, Introduction à la dynamique des gaz réactifs, Cépaduès, 2015 (ISBN 978-2-364-93190-9)
5. (en) Charles David Levermore, « Moment Closure Hierarchies for Kinetic Theories », Journal of Statistical Physics, vol. 83, 1996 (lire en ligne)
6. Byung Chan Eu, Nonequilibrium Statistical Mechanics. Ensemble Method., Springer, 1998 (ISBN 978-90-481-5007-6)
7. Cédric Villani, Limites hydrodynamiques de l'équation de Boltzmann, coll. « Séminaire Bourbaki » (n^o 893 [2]), 2001

Véase también

Véase también

• Choque élastique
• Colisión inélastique
• Desarrollo asymptotique
• Teoría de las perturbaciones

•   Datos: Q4133108