Distribución de Maxwell-Boltzmann para ecuaciones de transporte

Tengo que calcular los coeficientes de transporte para la distribución de Maxwell-Boltzmann. Pero no estoy seguro de qué distribución tengo que usar. Hasta donde yo sé, no debería ser la distribución de MB para v -espacio (Velocidad) o mi -eje (Energía), ya que eso me dará las dimensiones incorrectas al final. Tengo que usar la distribución por estado.

Pero no estoy seguro de cómo se ve esto. La integral que tengo que resolver, para obtener la conductividad eléctrica (primer coeficiente de transporte) que necesito, está dada por:

L ( 0 ) = ( 2 metro 2 ) 3 / 2 mi 2 τ π 2 metro ( F METRO B ε ) ε 3 / 2 d ε ,

al menos, de nuevo, al intentar calcular la conductividad eléctrica, que al final debería acabar siendo la fórmula de Drudes σ = norte mi 2 τ metro .

Básicamente, no es difícil. Pero tengo que hacer bien la función de distribución.

Hasta donde yo sé, la distribución de MB viene dada por:

F METRO B ( ε ) = C mi ε / k B T ,

dónde C es lo que necesito averiguar, ya que eso determinará las dimensiones de mis coeficientes.

Según mi libro, la función de distribución de MB normalizada es:

norte ¯ = norte ¯ Z 1 ( T , V ) mi ε / k B T ,

dónde:

Z 1 ( T , V ) norte ¯ = V norte ¯ ( 2 π metro k B T h 2 ) Z En t ( T ) ,

y Z En t ( T ) = 1 en mi caso.

Pero no estoy muy seguro de cómo hacer esto. Por lo que puedo ver, no se trata simplemente de insertar el término inverso de esto en C - al menos no por lo que puedo ver. Tal vez sea el V / norte No estoy seguro.

Bueno, ¿alguien que me pueda dar una pista?

Respuestas (1)

Factor de normalización

Definamos una función de densidad de probabilidad gaussiana generalizada (PDF) como:

(0) F s ( X ) = A o   mi ( X X o ) 2 2 σ 2
dónde A o es la constante de normalización, X es el argumento y s denota el conjunto de distribuciones (por ejemplo, especies de partículas), X o es el desplazamiento del pico desde X = 0 , y σ es la varianza de la distribución.

El factor de normalización A o se encuentra usando la siguiente restricción:

(1) +   d X   F s ( X ) = 1

La solución analítica a esta integral se puede encontrar en cualquier tabla integral estándar o usando algo como Mathematica , donde uno encuentra:

(2) A o = 1 2   π   σ 2

Distribución de velocidad maxwelliana

Para convertir de la PDF gaussiana 1D en la Ecuación 0 anterior a una distribución de Maxwell-Boltzmann , o Maxwellian, simplemente convertimos las variables de la siguiente manera:

  • X v , dónde v es el argumento de la velocidad de F s que van desde a +
  • X o v o , dónde v o es la velocidad de deriva o la velocidad de flujo a granel
  • 2   σ 2 V T s 2 , dónde V T s es la velocidad térmica (aquí la velocidad más probable )

Entonces podemos ver que el Maxwelliano 1D está dado por:

(3) F s ( v ) = 1 π   V T s   mi ( v v o V T s ) 2

La conversión a una distribución 3D completa es lo suficientemente simple siempre que cada componente de velocidad, v j , no está correlacionado con ningún otro componente. Luego (dejando caer el subíndice s por brevedad) podemos definir:

(4a) F ( v X , v y , v z ) = F ( v X )   F ( v y )   F ( v z ) (4b) = k = X , y , z   A k   mi ( v k v o k V T k ) 2
donde el factor de normalización total viene dado por:
(5) A X   A y   A z   = 1 π 3 / 2   V T X   V T y   V T z

Conversión a Energía

Para hacer esto correctamente, se debe usar un análogo de cantidad de movimiento de la versión espacio-velocidad descrita anteriormente.

En el límite no relativista, la conversión de cantidad de movimiento a energía es mi = pag 2 2   metro o pag = 2   metro   mi . Por lo tanto, podemos ver que d pag / d mi mi 1 / 2 o pag 1 . Dado que la energía es un escalar, los límites de integración (p. ej., al encontrar el factor de normalización) cambian de v j + a 0 mi + .

Luego, la versión 3D (p. ej., la Ecuación 4b anterior) se convierte en:

F ( mi ) = 1 Z   mi ( mi k B   T )
dónde Z es la función de partición , k B es la constante de Boltzmann y T es la temperatura

En el límite relativista , la situación se vuelve increíblemente complicada como se explica en ¿Cuál es la función de distribución relativista correcta? .

Coeficientes de transporte

Se pueden calcular estos mediante el uso de los momentos de la distribución. Escribí una respuesta detallada para una distribución de velocidad no relativista en https://physics.stackexchange.com/a/218643/59023 .

Para hacer que las distribuciones de probabilidad que se muestran arriba sean consistentes con los momentos de velocidad discutidos en mi otra respuesta, puede redefinir F s tal que cuando se integra sobre todo el espacio de velocidad se obtiene la densidad numérica del conjunto s .

Creo que su conversión de energía ha omitido la velocidad de deriva.
@LuKas: es un desplazamiento constante, así que agregue un mi o después de la mi en el exponente y si lo desea. O podría transformarse en un marco de referencia donde la velocidad de deriva sea cero.
No creo que puedas simplemente agregar un mi 0 término, porque entonces te faltan los términos cruzados del ( v v 0 ) 2 .
@LuKas: es una compensación. Así que transfórmate en el v o el marco luego se convierte en energía, luego se transforma de nuevo en el marco de laboratorio y tendrá un mi o término allí.
Puede que esté haciendo algo mal, pero cuando calculo rápidamente la distribución de energía 1D a partir de la distribución de velocidad usando sustitución, obtengo F ( mi ) = 1 4 π k B T mi Exp ( metro ( 2 mi metro 2 v 0 ) 2 2 k B T ) . En el exponencial obtendrá un término cruzado que implica la raíz cuadrada de la variable aleatoria, que por lo tanto no se puede escribir como un término de compensación constante. Obtendrás lo mismo (o peor) en 3D.
Necesitas transformarte en el v o marco primero para que desaparezca en el exponente.
Pero si luego se transforma de nuevo al marco de referencia en movimiento, no solo obtiene el término de compensación constante, ¿verdad? Esta no es una transformación lineal en la energía. Entonces, sí, puede simplemente ir al marco de descanso y el término desaparece, pero si quiere trabajar en cualquier otro marco, el término de deriva lo complica. (También noté el pequeño error tipográfico en mi F ( mi ) arriba: el 2 delante del v 0 no debería estar allí)
¿Por qué no? Ha eliminado la información direccional al convertirla en energía, por lo que tendrá que ser una compensación constante en energía (las unidades deben coincidir).
Por supuesto, las unidades coincidirán, el término cruzado obviamente tiene las mismas unidades que los cuadrados de los términos que se están sumando (y tienen unidades coincidentes).