Método de Chapman-Enskog
# 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 propongan simultáneamente e independientemente en 1916 y 1917 una solución a este problema. Más recientemente este método se ha extendido al caso de un gas en desequilibrio termodinámico, 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 perturbaciones que consiste en definir la solución bajo la forma de series de funciones de distribución en función 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 permite 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 $1,\mathbf {v} ,\mathbf {v} \otimes \mathbf {v} ,\mathbf {v} \otimes \mathbf {v} \otimes \mathbf {v} ...$ ($\mathbf {v}$ es la velocidad microscópica de la ecuación de Boltzmann y $\otimes$ el producto tensorial, e integrado en velocidad. En este tipo de método, el término n de la expansión queda en función del término (n+1), por lo que 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 entropía del sistema de fermiones que son las partículas del medio que él estudió. 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 macroscópica

### Nivel microscópico

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

La ecuación de Boltzmann se escribe como

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

dónde $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 ($j=i$) o heterogéneas ($j\neq i$).

**El choque elástico**

Las velocidades antes de la interacción en un Sistema de referencia inercial son $\mathbf {v} _{i}$ y $\mathbf {v} _{j}$. Los índices representan indistintamente una misma especie o dos especies diferentes. Las velocidades después de la interacción son $\mathbf {v} _{i}'$ y $\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 galileano, la velocidad inicial de la partícula $j$ es la velocidad relativa $\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 $\mathbf {g} _{ij}$. Elegimos un punto de referencia como: $\mathbf {\Omega } ={\frac {\mathbf {g} _{ij}}{||\mathbf {g} _{ij}||}}=[1,0,0]$ (ver figura). En esta referencia, el ángulo de desviación respecto al eje $\mathbf {\Omega } =[1,0,0]$ es $\theta _{ij}$. Definimos también el parámetro de impacto $b$ de la velocidad relativa $g_{ij}$ y el potencial de interacción que, se supone, depende sólo de la distancia de las dos partículas que interaccionan. Si esta hipótesis es rigurosa para la interacción entre dos átomos, se puede considerar utilizable también para dos moléculas. El potencial es entonces un potencial medio estadístico.

La dirección de salida de interacción está definida por: $\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

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

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

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

       $\mathbf {v} _{i}'(b,g)=\mathbf {v} _{i}-(\mathbf {g_{ij}} \cdot \mathbf {\Omega '} )\,\mathbf {\Omega '}$
       $\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 $b'=b$.

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

       $\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 $2\,\pi \,b\,\mathrm {d} b$ por unidad de tiempo es $f(\mathbf {v} _{j})\,g_{ij}\,2\,\pi \,b\,\mathrm {d} b\,\mathrm {d} \mathbf {v} _{j}$. Estas partículas interactúan con el número probable de partículas en el volumen elemental. $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 $\Theta _{ij}^{-}\,\mathrm {d} \mathbf {r} \,\mathrm {d} \mathbf {v} _{i}$ con

       $\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 por unidad de tiempo, tenemos

       $\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

       $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}$,

donde $\mathbf {v} _{i}'=\mathbf {v} _{i}'\left(\mathbf {v} _{i},\mathbf {v} _{j}\right)$ y $\mathbf {v} _{j}'=\mathbf {v} _{j}'\left(\mathbf {v} _{i},\mathbf {v} _{j}\right)$ dependen explícitamente de las velocidades de colisión iniciales $\mathbf {v} _{i}$ y $\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 $\sigma _{ij}$ y se define:

       $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

       $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 $i$ a nivel macroscópico esto se define:

| - Densidad de partículas | $n_{i}=\int _{\mathbf {v} }f_{i}\mathrm {d} \mathbf {v} _{i}$ | $n=\sum _{i}n_{i}$ |
| - | - | - |
| - Densidad | $\rho _{i}=n_{i}m_{i}$ | $\rho =\sum _{i}\rho _{i}$ |
| - Velocidad media | $\mathbf {V} _{i}={\frac {1}{n_{i}}}\int _{\mathbf {v} }\mathbf {v} _{i}f_{i}\mathrm {d} \mathbf {v} _{i}$ | |
| - Energía interna | $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

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

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

| - Fracción de volumen | $x_{i}={\frac {n_{i}}{n}}$ | tal que $\sum _{i}x_{i}=1$ |
| - | - | - |
| - Fracción de masa | $c_{i}={\frac {\rho _{i}}{\rho }}$ | tal que $\sum _{i}c_{i}=1$ |
| - Masa molar | ${\mathcal {M}}_{i}={\mathcal {N}}m_{i}$ | de valor medio ${\overline {\mathcal {M}}}=\sum _{i}x_{i}{\mathcal {M}}_{i}$ |

**Los flujos**

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

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

| - Flujo de masa | $\mathbf {J} _{i}=m_{i}\int _{\mathbf {v} }f_{i}\mathbf {U} _{i}\mathrm {d} \mathbf {U} _{i}$ | $=\rho _{i}{\mathbf {V} _{D}}_{i}$ |
| - | - | - |
| - Tensor de presiones | ${\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. |
| - Presión parcial | $p_{i}={\frac {1}{3}}\int _{\mathbf {v} }m_{i}U_{i}^{2}f_{i}d\mathbf {U} _{i}$ | se define a partir de la traza del tensor de presiones. |
| - Flujo de calor | $\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 $i$ así como la presión: $p=\sum _{i}p_{i}$.

Podemos entonces definir una temperatura a partir de la ecuación de estado

       $p=n\,k\,T$

**Ecuaciones de evolución**

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

       ${\frac {\partial \rho _{i}}{\partial t}}+\nabla \cdot (\rho _{i}\mathbf {V} +\mathbf {J} _{i})=0$
       $\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$
       ${\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

       $\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 $\rho$, $\mathbf {V}$ y $e$ en el que los flujos $\mathbf {J} _{i}$, ${\mathsf {P}}$ y $\mathbf {q}$ quedan por definir.

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

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

       ${\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:

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

## Formulación de la solución

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

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

Como $f_{i}$ verifica las leyes de conservación de los invariantes de colisión $m_{i},m_{i}\mathbf {v} ,{\frac {1}{2}}m_{i}v^{2}$, cada uno de los términos de desarrollo también deben de verificarlas. De donde se imponen las siguientes restricciones sobre la solución.

       $\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}]$.

Por último, insertamos esta aproximación en la ecuación de Boltzmann, y separamos los términos correspondientes a cada potencia de $\Lambda$ para estudiar los distintos órdenes.

### Al orden cero

Se obtiene simplemente

       $\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

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

Lo que implica

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

o

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

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

       $\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

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

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

### Al orden uno

A orden uno, se obtiene una ecuación integral de Fredholm para la incógnita $f_{i}^{(1)}$

       ${\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 denominado sistema de Stefan-Maxwell

       ${\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 ${\mathcal {D}}_{ij}(T,p)$ y el "coeficiente de difusión térmica multicomponente",$k_{i}^{T}(T,p,x_{i})$ (en realidad un número adimensional) ${\mathcal {D}}_{i}^{T}(T,p,x_{i})$ (que no son coeficientes de difusión y que pueden ser negativos) por:

       $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 $N$ especies el rango de este sistema es $N-1$ ya que $\sum _{i}\mathbf {J} _{i}=0$. Su solución

       $\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). $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 $N-1$ La solución no es única e implica ${\frac {N(N-1)}{2}}$ términos independientes. Generalmente elegimos $D_{ii}=0$ por simetría, pero 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 semejante a la ley de Fick, la cual no es exacta para una mezcla binaria.

**El tensor de presiones**

El tensor de presiones tiene una forma clásica

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

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

       ${\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 $-\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.

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

**El flujo de calor**

Está dado por

       $\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}$

$\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 $\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 la práctica, 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

       $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 $\mathbf {d} _{i}$ se da en el flujo de difusión

       $\mathbf {W} _{i}={\sqrt {\frac {m_{i}}{2kT}}}\mathbf {V} _{i}$
       $\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 $f_{i}^{(2)}$

       ${\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 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.

