stringtranslate.com

ecuación de Boltzmann

El lugar de la ecuación cinética de Boltzmann en la escalera de la reducción de modelos de la dinámica microscópica a la dinámica macroscópica del continuo (ilustración del contenido del libro [1] )

La ecuación de Boltzmann o ecuación de transporte de Boltzmann ( BTE ) describe el comportamiento estadístico de un sistema termodinámico que no se encuentra en estado de equilibrio ; fue ideado por Ludwig Boltzmann en 1872. [2] El ejemplo clásico de tal sistema es un fluido con gradientes de temperatura en el espacio que causan que el calor fluya de regiones más calientes a otras más frías, mediante el transporte aleatorio pero sesgado de las partículas que lo componen. líquido. En la literatura moderna, el término ecuación de Boltzmann se utiliza a menudo en un sentido más general, refiriéndose a cualquier ecuación cinética que describe el cambio de una cantidad macroscópica en un sistema termodinámico, como la energía, la carga o el número de partículas.

La ecuación surge no analizando las posiciones individuales y los momentos de cada partícula en el fluido, sino más bien considerando una distribución de probabilidad para la posición y el momento de una partícula típica, es decir, la probabilidad de que la partícula ocupe una región muy pequeña del espacio. (matemáticamente el elemento de volumen ) centrado en la posición , y tiene un impulso casi igual a un vector de impulso dado (ocupando así una región muy pequeña del espacio de impulso ), en un instante de tiempo.

La ecuación de Boltzmann se puede utilizar para determinar cómo cambian las cantidades físicas, como la energía térmica y el momento , cuando un fluido está en transporte. También se pueden derivar otras propiedades características de los fluidos, como la viscosidad , la conductividad térmica y la conductividad eléctrica (tratando los portadores de carga de un material como un gas). [2] Véase también la ecuación de convección-difusión .

La ecuación es una ecuación integrodiferencial no lineal , y la función desconocida en la ecuación es una función de densidad de probabilidad en un espacio de seis dimensiones de la posición y el momento de una partícula. El problema de la existencia y unicidad de las soluciones aún no está completamente resuelto, pero algunos resultados recientes son bastante prometedores. [3] [4]

Descripción general

El espacio de fase y la función de densidad.

El conjunto de todas las posiciones posibles r y momentos p se denomina espacio de fase del sistema; en otras palabras, un conjunto de tres coordenadas para cada coordenada de posición x, y, z , y tres más para cada componente de impulso p x , p y , p z . Todo el espacio es de 6 dimensiones : un punto en este espacio es ( r , p ) = ( x, y, z, p x , p y , p z ) , y cada coordenada está parametrizada por el tiempo t . El volumen pequeño (" elemento de volumen diferencial ") se escribe

Dado que la probabilidad de N moléculas, todas las cuales tienen r y p dentro de , está en duda, en el centro de la ecuación hay una cantidad f que da esta probabilidad por unidad de volumen del espacio de fase, o probabilidad por unidad de longitud al cubo por unidad de momento al cubo. , en un instante de tiempo t . Esta es una función de densidad de probabilidad : f ( r , p , t ) , definida de modo que,

tienenrde espacio de impulsopt[5] La integración

que es una integral séxtuple . Mientras que f está asociado con varias partículas, el espacio de fases es para una partícula (no para todas, lo que suele ser el caso de los sistemas deterministas de muchos cuerpos ), ya que solo se trata de una r y una p . No es parte del análisis usar r 1 , p 1 para la partícula 1, r 2 , p 2 para la partícula 2, etc. hasta r N , p N para la partícula N .

Se supone que las partículas del sistema son idénticas (por lo que cada una tiene una masa idéntica m ). Para una mezcla de más de una especie química , se necesita una distribución para cada una, ver más abajo.

Declaración principal

La ecuación general se puede escribir entonces como [6]

donde el término "fuerza" corresponde a las fuerzas ejercidas sobre las partículas por una influencia externa (no por las partículas mismas), el término "diff" representa la difusión de partículas y "coll" es el término de colisión , que representa las fuerzas. actuando entre partículas en colisiones. A continuación se proporcionan expresiones para cada término del lado derecho. [6]

Tenga en cuenta que algunos autores utilizan la velocidad de la partícula v en lugar del momento p ; están relacionados en la definición de impulso por p = m v .

Los términos de fuerza y ​​difusión.

Considere las partículas descritas por f , cada una de las cuales experimenta una fuerza externa F que no se debe a otras partículas (consulte el término de colisión para este último tratamiento).

Supongamos que en el momento t un número de partículas tienen todas la posición r dentro del elemento y un momento p dentro . Si una fuerza F actúa instantáneamente sobre cada partícula, entonces en el tiempo t + Δ t su posición será y el momento p + Δ p = p + F Δ t . Entonces, en ausencia de colisiones, f debe satisfacer

Tenga en cuenta que hemos utilizado el hecho de que el elemento volumen del espacio de fases es constante, lo que se puede demostrar utilizando las ecuaciones de Hamilton (consulte la discusión bajo el teorema de Liouville ). Sin embargo, dado que ocurren colisiones, la densidad de partículas en el volumen del espacio de fases cambia, por lo que

donde Δ f es el cambio total en f . Dividiendo ( 1 ) por y tomando los límites Δ t → 0 y Δ f → 0 , tenemos

El diferencial total de f es:

donde es el operador de gradiente , · es el producto escalar ,

ê xê yê zvectores unitarios cartesianos

Declaración final

Dividiendo ( 3 ) por dt y sustituyendo en ( 2 ) se obtiene:

En este contexto, F ( r , t ) es el campo de fuerza que actúa sobre las partículas en el fluido y m es la masa de las partículas. El término del lado derecho se añade para describir el efecto de las colisiones entre partículas; si es cero entonces las partículas no chocan. La ecuación de Boltzmann sin colisiones, donde las colisiones individuales se reemplazan con interacciones agregadas de largo alcance, por ejemplo, interacciones de Coulomb , a menudo se denomina ecuación de Vlasov .

Esta ecuación es más útil que la principal anterior, pero aún está incompleta, ya que f no se puede resolver a menos que se conozca el término de colisión en f . Este término no se puede encontrar tan fácil o generalmente como los demás: es un término estadístico que representa las colisiones de partículas y requiere conocimiento de las estadísticas que obedecen las partículas, como las distribuciones de Maxwell-Boltzmann , Fermi-Dirac o Bose-Einstein .

El término de colisión (Stosszahlansatz) y el caos molecular

Término de colisión de dos cuerpos

Una idea clave aplicada por Boltzmann fue determinar el término de colisión resultante únicamente de colisiones de dos cuerpos entre partículas que se supone no estaban correlacionadas antes de la colisión. Boltzmann se refirió a esta suposición como " Stosszahlansatz " y también se la conoce como "suposición del caos molecular ". Bajo este supuesto, el término de colisión se puede escribir como una integral de momento-espacio sobre el producto de funciones de distribución de una partícula: [2]

p Ap BABp′ Ap′ B
velocidad relativaI ( g , Ω)sección transversal diferencialθángulo sólido d Ω

Simplificaciones del término de colisión

Dado que gran parte del desafío para resolver la ecuación de Boltzmann se origina en el término de colisión complejo, se han realizado intentos de "modelar" y simplificar el término de colisión. El modelo de ecuación más conocido se debe a Bhatnagar, Gross y Krook. [7] El supuesto en la aproximación BGK es que el efecto de las colisiones moleculares es forzar una función de distribución de no equilibrio en un punto del espacio físico a volver a una función de distribución de equilibrio de Maxwell y que la velocidad a la que esto ocurre es proporcional a la frecuencia de colisión molecular. Por tanto, la ecuación de Boltzmann se modifica a la forma BGK:

donde es la frecuencia de colisión molecular y es la función de distribución local de Maxwell dada la temperatura del gas en este punto del espacio. Esto también se denomina "aproximación del tiempo de relajación".

Ecuación general (para una mezcla)

Para una mezcla de especies químicas etiquetadas por índices i = 1, 2, 3, ..., n la ecuación para la especie i es [2]

donde f i = f i ( r , p i , t ) y el término de colisión es

donde f′ = f′ ( p′ i , t ) , la magnitud de los momentos relativos es

y I ij es la sección transversal diferencial, como antes, entre las partículas i y j . La integración se realiza sobre los componentes del momento en el integrando (que están etiquetados como i y j ). La suma de integrales describe la entrada y salida de partículas de la especie i dentro o fuera del elemento del espacio de fases.

Aplicaciones y extensiones

Ecuaciones de conservación

La ecuación de Boltzmann se puede utilizar para derivar las leyes de conservación de la dinámica de fluidos para masa, carga, momento y energía. [8] : 163  Para un fluido que consta de un solo tipo de partícula, la densidad numérica n está dada por

El valor promedio de cualquier función A es

Dado que las ecuaciones de conservación involucran tensores, se utilizará la convención de suma de Einstein cuando índices repetidos en un producto indiquen suma sobre esos índices. Por tanto y , ¿dónde está el vector de velocidad de las partículas? Defina como alguna función únicamente del momento , que se conserva en una colisión. Supongamos también que la fuerza es función únicamente de la posición y que f es cero para . Multiplicar la ecuación de Boltzmann por A e integrar sobre el momento produce cuatro términos que, usando la integración por partes, se pueden expresar como

donde el último término es cero, ya que A se conserva en una colisión. Los valores de A corresponden a momentos de velocidad (y momento , ya que son linealmente dependientes).

momento cero

Dejando , la masa de la partícula, la ecuación integrada de Boltzmann se convierte en la ecuación de conservación de la masa: [8] : 12, 168 

Primer momento

Dejando , el momento de la partícula, la ecuación integrada de Boltzmann se convierte en la ecuación de conservación del momento: [8] : 15, 169 

donde es el tensor de presión (el tensor de tensión viscosa más la presión hidrostática ).

Segundo momento

Dejando , la energía cinética de la partícula, la ecuación integrada de Boltzmann se convierte en la ecuación de conservación de la energía: [8] : 19, 169 

donde es la densidad de energía térmica cinética y es el vector de flujo de calor.

Mecánica hamiltoniana

En mecánica hamiltoniana , la ecuación de Boltzmann a menudo se escribe de manera más general como

Loperador de LiouvilleCL

Teoría cuántica y violación de la conservación del número de partículas.

Es posible escribir ecuaciones cuánticas relativistas de Boltzmann para sistemas cuánticos relativistas en los que el número de partículas no se conserva en las colisiones. Esto tiene varias aplicaciones en cosmología física , [9] incluida la formación de elementos ligeros en la nucleosíntesis del Big Bang , la producción de materia oscura y la bariogénesis . No está claro a priori que el estado de un sistema cuántico pueda caracterizarse por una densidad espacial de fase clásica f . Sin embargo, para una amplia clase de aplicaciones existe una generalización bien definida de f que es la solución de una ecuación de Boltzmann efectiva que puede derivarse de los primeros principios de la teoría cuántica de campos . [10]

Relatividad general y astronomía.

La ecuación de Boltzmann es útil en la dinámica galáctica. Una galaxia, bajo ciertas suposiciones, puede aproximarse a un fluido continuo; su distribución de masa está entonces representada por f ; En las galaxias, las colisiones físicas entre las estrellas son muy raras y el efecto de las colisiones gravitacionales puede despreciarse durante tiempos mucho más largos que la edad del universo .

Su generalización en la relatividad general es [11]

Γ α βγsímbolo de Christoffel( xi , p i ) espacio de( xi , p i )[12] [13]

En cosmología física, se ha utilizado el enfoque totalmente covariante para estudiar la radiación cósmica de fondo de microondas. [14] De manera más genérica, el estudio de los procesos en el universo temprano a menudo intenta tener en cuenta los efectos de la mecánica cuántica y la relatividad general . [9] En el medio muy denso formado por el plasma primordial después del Big Bang , las partículas se crean y aniquilan continuamente. En tal entorno, la coherencia cuántica y la extensión espacial de la función de onda pueden afectar la dinámica, lo que hace cuestionable si la distribución clásica del espacio de fase f que aparece en la ecuación de Boltzmann es adecuada para describir el sistema. Sin embargo, en muchos casos es posible derivar una ecuación de Boltzmann eficaz para una función de distribución generalizada a partir de los primeros principios de la teoría cuántica de campos . [10] Esto incluye la formación de los elementos ligeros en la nucleosíntesis del Big Bang , la producción de materia oscura y la bariogénesis .

Resolviendo la ecuación

En algunos casos se ha demostrado que existen soluciones exactas a las ecuaciones de Boltzmann; [15] este enfoque analítico proporciona información, pero generalmente no se puede utilizar en problemas prácticos.

En cambio, los métodos numéricos (incluidos los métodos de elementos finitos y de celosía de Boltzmann ) se utilizan generalmente para encontrar soluciones aproximadas a las diversas formas de la ecuación de Boltzmann. Los ejemplos de aplicaciones van desde la aerodinámica hipersónica en flujos de gas enrarecido [16] [17] hasta flujos de plasma. [18] Una aplicación de la ecuación de Boltzmann en electrodinámica es el cálculo de la conductividad eléctrica; el resultado es, en orden principal, idéntico al resultado semiclásico. [19]

Cerca del equilibrio local , la solución de la ecuación de Boltzmann puede representarse mediante una expansión asintótica en potencias del número de Knudsen (la expansión de Chapman-Enskog [20] ). Los dos primeros términos de esta expansión dan las ecuaciones de Euler y las ecuaciones de Navier-Stokes . Los términos superiores tienen singularidades. El problema de desarrollar matemáticamente los procesos limitantes que conducen desde la visión atomista (representada por la ecuación de Boltzmann) a las leyes del movimiento de los continuos, es una parte importante del sexto problema de Hilbert . [21]

Limitaciones y usos adicionales de la ecuación de Boltzmann

La ecuación de Boltzmann es válida sólo bajo varios supuestos. Por ejemplo, se supone que las partículas son puntuales, es decir, sin tamaño finito. Existe una generalización de la ecuación de Boltzmann que se llama ecuación de Enskog . [22] El término de colisión se modifica en las ecuaciones de Enskog de modo que las partículas tienen un tamaño finito; por ejemplo, pueden modelarse como esferas que tienen un radio fijo.

Para las partículas no se suponen más grados de libertad que el movimiento de traslación. Si hay grados de libertad internos, la ecuación de Boltzmann debe generalizarse y podría poseer colisiones inelásticas . [22]

Muchos fluidos reales, como líquidos o gases densos, además de las características mencionadas anteriormente, tienen formas de colisión más complejas; no solo habrá colisiones binarias, sino también ternarias y de orden superior. [23] Estos deben derivarse utilizando la jerarquía BBGKY .

Para el movimiento de las células también se utilizan ecuaciones tipo Boltzmann . [24] [25] Dado que las células son partículas compuestas que tienen grados de libertad internos, las ecuaciones de Boltzmann generalizadas correspondientes deben tener integrales de colisión inelásticas. Estas ecuaciones pueden describir invasiones de células cancerosas en tejidos, morfogénesis y efectos relacionados con la quimiotaxis .

Ver también

Notas

  1. ^ Gorban, Alexander N.; Karlin, Ilya V. (2005). Colectores invariantes para cinética física y química. Apuntes de conferencias de física (LNP, vol. 660). Berlín, Heidelberg: Springer. doi :10.1007/b98103. ISBN 978-3-540-22684-0.URL alternativa
  2. ^ abcd Encyclopaedia of Physics (segunda edición), RG Lerner , GL Trigg, editores VHC, 1991, ISBN (Verlagsgesellschaft) 3-527-26954-1, ISBN (VHC Inc.) 0-89573-752-3.
  3. ^ DiPerna, RJ; Leones, P.-L. (1989). "Sobre el problema de Cauchy para las ecuaciones de Boltzmann: existencia global y estabilidad débil". Ana. de Matemáticas . 2. 130 (2): 321–366. doi :10.2307/1971423. JSTOR  1971423.
  4. ^ Philip T. Gressman y Robert M. Strain (2010). "Soluciones clásicas globales de la ecuación de Boltzmann con interacciones de largo alcance". Procedimientos de la Academia Nacional de Ciencias . 107 (13): 5744–5749. arXiv : 1002.3639 . Código Bib : 2010PNAS..107.5744G. doi : 10.1073/pnas.1001185107 . PMC 2851887 . PMID  20231489. 
  5. ^ Huang, Kerson (1987). Mecánica estadística (Segunda ed.). Nueva York: Wiley. pag. 53.ISBN _ 978-0-471-81518-1.
  6. ^ ab Enciclopedia de física de McGraw Hill (segunda edición), SP Parker, 1993, ISBN 0-07-051400-3
  7. ^ Bhatnagar, PL; Bruto, EP; Krook, M. (1 de mayo de 1954). "Un modelo para procesos de colisión en gases. I. Procesos de pequeña amplitud en sistemas monocomponente cargados y neutros". Revisión física . 94 (3): 511–525. Código bibliográfico : 1954PhRv...94..511B. doi : 10.1103/PhysRev.94.511.
  8. ^ abcd de Groot, SR; Mazur, P. (1984). Termodinámica de no equilibrio . Nueva York: Dover Publications Inc. ISBN 978-0-486-64741-8.
  9. ^ ab Edward Kolb y Michael Turner (1990). El Universo Temprano . Prensa de Westview. ISBN 978-0-201-62674-2.
  10. ^ ab M. Drewes; C. Weniger; S. Mendizábal (8 de enero de 2013). "La ecuación de Boltzmann de la teoría cuántica de campos". Física. Letón. B . 718 (3): 1119-1124. arXiv : 1202.1301 . Código Bib : 2013PhLB..718.1119D. doi :10.1016/j.physletb.2012.11.046. S2CID  119253828.
  11. ^ Ehlers J (1971) Relatividad general y cosmología (Varenna), RK Sachs (Academic Press NY); Thorne KS (1980) Rev. Mod. Phys., 52, 299; Ellis GFR, Treciokas R, Matravers DR, (1983) Ann. Física, 150, 487}
  12. ^ Debbasch, Fabrice; Willem van Leeuwen (2009). "Ecuación relativista general de Boltzmann I: tratamiento covariante". Física A. 388 (7): 1079-1104. Código Bib : 2009PhyA..388.1079D. doi :10.1016/j.physa.2008.12.023.
  13. ^ Debbasch, Fabrice; Willem van Leeuwen (2009). "Ecuación relativista general de Boltzmann II: tratamiento manifiestamente covariante". Física A. 388 (9): 1818–34. Código Bib : 2009PhyA..388.1818D. doi :10.1016/j.physa.2009.01.009.
  14. ^ Maartens R, Gebbie T, Ellis GFR (1999). "Anisotropías de fondo cósmico de microondas: dinámica no lineal". Física. Rev. D. 59 (8): 083506
  15. ^ Philip T. Gressman; Robert M. cepa (2011). "Soluciones clásicas globales de la ecuación de Boltzmann sin corte angular". Revista de la Sociedad Matemática Estadounidense . 24 (3): 771. arXiv : 1011.5441 . doi :10.1090/S0894-0347-2011-00697-8. S2CID  115167686.
  16. ^ Evans, Ben; Morgan, Ken; Hassan, Oubay (1 de marzo de 2011). "Una solución discontinua de elementos finitos de la ecuación cinética de Boltzmann en formas sin colisiones y BGK para flujos de gas macroscópicos". Modelado Matemático Aplicado . 35 (3): 996–1015. doi : 10.1016/j.apm.2010.07.027 .
  17. ^ Evans, B.; Walton, SP (diciembre de 2017). "Optimización aerodinámica de un vehículo de reentrada hipersónico basada en la solución de la ecuación de Boltzmann-BGK y optimización evolutiva". Modelado Matemático Aplicado . 52 : 215–240. doi : 10.1016/j.apm.2017.07.024 . ISSN  0307-904X.
  18. ^ Pareschi, L.; Ruso, G. (1 de enero de 2000). "Solución numérica de la ecuación de Boltzmann I: aproximación espectralmente precisa del operador de colisión". Revista SIAM de Análisis Numérico . 37 (4): 1217-1245. CiteSeerX 10.1.1.46.2853 . doi :10.1137/S0036142998343300. ISSN  0036-1429. 
  19. ^ HJW Müller-Kirsten, Conceptos básicos de la mecánica estadística, Capítulo 13, 2ª ed., World Scientific (2013), ISBN 978-981-4449-53-3
  20. ^ Sídney Chapman; Thomas George Cowling La teoría matemática de los gases no uniformes: una descripción de la teoría cinética de la viscosidad, la conducción térmica y la difusión de los gases, Cambridge University Press, 1970. ISBN 0-521-40844-X 
  21. ^ "Número temático 'El sexto problema de Hilbert'". Transacciones filosóficas de la Royal Society A. 376 (2118). 2018. doi : 10.1098/rsta/376/2118 .
  22. ^ ab "Ecuación de Enskog: descripción general | Temas de ScienceDirect". www.sciencedirect.com . Consultado el 10 de mayo de 2022 .
  23. ^ van Noije, TPC; Ernst, MH (3 de junio de 1997). "Teoría cinética de anillos para un gas granular idealizado". arXiv : cond-mat/9706020 .
  24. ^ Chauvière, A.; Hillen, T.; Preziosi, L. (2007). "Modelado del movimiento celular en tejidos de redes anisotrópicas y heterogéneas". Instituto Americano de Ciencias Matemáticas . 2 (2): 333–357. doi :10.3934/nhm.2007.2.333.
  25. ^ Conte, Martina; Loy, Nadia (12 de febrero de 2022). "Modelo cinético de múltiples señales con detección no local para la migración celular en una red de fibra con quimiotaxis". Boletín de Biología Matemática . 84 (3): 42. doi :10.1007/s11538-021-00978-1. ISSN  1522-9602. PMC 8840942 . PMID  35150333. 

Referencias

enlaces externos