La PizarraMatemática.

Máster y doctorado

Ecuaciones cinéticas y límites de muchos cuerpos

Saber cuántas partículas hay en un lugar no dice hacia dónde viajan. La teoría cinética conserva posición y velocidad en una misma distribución. Seguiremos primero partículas sin fuerza; después promediaremos sus velocidades y veremos qué información se pierde al pasar de muchas partículas a unas pocas variables.

1. Una distribución de fase contiene posición y velocidad

f(t,x,v)≥0 describe cantidad o probabilidad en regiones de posición y velocidad, según su normalización. La densidad espacial \rho(t,x)=\int f(t,x,v)\,dv elimina la información de velocidad. Dos poblaciones con igual ρ pueden evolucionar de formas distintas.

Sin fuerza ni colisiones, f_t+v\cdot\nabla_xf=0.

Cada partícula mantiene velocidad. En ejemplos discretos usamos masas puntuales, interpretadas como medidas, para visualizar la misma idea; no fingimos que esas masas sean funciones densidad ordinarias.

Mismo lugar inicial, futuro diferente

Dos poblaciones tienen masa total uno concentrada inicialmente en x=0.

  1. Primera

    Mitad viaja a v=−1 y mitad a v=1.
  2. Segunda

    Toda permanece a v=0.
  3. En t=2

    La primera queda en x=−2 y x=2; la segunda en x=0.
Resultado

La densidad espacial inicial sola no determina la evolución.

Lo comprobamos

La comparación se formula con medidas puntuales.

Un dominio puede perder partículas

Una partícula en x=1/2 viaja con v=1; contamos solo el intervalo [0,1].

  1. Trayectoria

    x(t)=1/2+t.
  2. Salida

    Tras t=1/2 atraviesa el extremo derecho.
  3. Balance

    La masa de toda la recta se conserva; la de este intervalo disminuye si no se reinserta ni refleja.

Resultado

Conservar masa total exige decir dónde la contamos.

Lo comprobamos

La palabra “libre” no especifica condiciones de frontera.

2. Las características transportan el dato sin cambiar su valor

La trayectoria libre cumple x(t)=x₀+tv y v(t)=v.
A lo largo de ella, \frac{d}{dt}f(t,x_0+tv,v)=f_t+v\cdot\nabla_xf=0.

Para hallar el valor actual retrocedemos hasta x₀=x−tv: f(t,x,v)=f_0(x-tv,v).

El transporte libre inclina una región de fase sin cambiar su área
(x₀,v) → (x₀+v,v)120123xvAzul: t=0; verde: t=1Misma área, otra forma.No hay fuerza ni colisiones.
Plano de fase de una posición y una velocidad. El rectángulo azul inicial contiene posiciones entre cero y uno y velocidades entre uno y dos. Tras tiempo uno, cada punto va a x más v manteniendo su velocidad, y forma el paralelogramo verde con vértices uno uno, dos uno, tres dos y dos dos. Ambos polígonos tienen la misma área y representan el soporte de una densidad transportada. Las alturas corresponden a velocidad, no al tiempo. Los ejes y las escalas se mantienen para comparar ambos conjuntos; la masa total se conserva en todo el espacio, no necesariamente en un intervalo espacial fijo.
El mapa (x,v)→(x−tv,v) tiene jacobiano uno.

Para una densidad integrable en todo el espacio conserva masa. La fórmula puede leerse clásicamente con regularidad suficiente o mediante transporte de medidas.

Retroceder antes de evaluar

Dato f₀(x,v)=1_[0,1](x)·g(v), t=1 y v=2.

  1. Origen

    x₀=x−2.
  2. Soporte

    0≤x−2≤1 equivale a 2≤x≤3.
  3. Evaluar

    En x=5/2, f=g(2); en x=1, f=0.
Resultado

La franja de velocidad dos se traslada dos unidades.

Lo comprobamos

No usamos x+tv para buscar el punto de partida.

Conservar área en fase

En una dimensión el mapa directo es (x₀,v)→(x₀+tv,v).

  1. Matriz

    \begin{pmatrix}1&t\\0&1\end{pmatrix}.

  2. Determinante

    Vale uno.

  3. Interpretar

    Un rectángulo se inclina en un paralelogramo sin cambiar su área.

Resultado

El perfil puede deformarse al mirar posición y velocidad juntas.

Lo comprobamos

La conservación de área no obliga a conservar la densidad espacial en cada punto.

3. Una velocidad media cero puede ocultar mucho movimiento

Los momentos básicos son \rho=\int f\,dv,\quad j=\int vf\,dv.

Donde ρ>0, la velocidad media es u=j/ρ.

La energía cinética con masa de partícula uno usa \tfrac12\int|v|^2f\,dv. Cada integral necesita el momento correspondiente finito.

La misma velocidad media puede ocultar energías diferentes
Masa total 1 en ambasA1/21/2−202Media 0; energía 2B1−202Media 0; energía 0Media cero no significa reposo.
Dos distribuciones de velocidad de masa total uno. En la primera hay masa un medio en menos dos y un medio en dos, representada por dos barras de la misma altura. En la segunda toda la masa está en cero, representada por una barra del doble de altura con la misma escala de masa. Ambas tienen corriente y velocidad media cero. La primera tiene segundo momento cuatro y energía cinética dos para masa de partícula uno; la segunda tiene segundo momento y energía cero. Las barras representan masas discretas, no una densidad continua ni trayectorias espaciales. Mirar únicamente la media no permite reconstruir la distribución.

La dispersión alrededor de u se describe mediante \int(v-u)\otimes(v-u)f\,dv. En una dimensión es una cantidad no negativa relacionada con presión. Ni la energía ni esa dispersión se determinan solo por ρ y j.

Dos velocidades opuestas

Masa 1/2 en v=−2 y 1/2 en v=2.

  1. Masa y corriente

    ρ=1; j=−1+1=0; u=0.
  2. Energía

    (1/2)[(1/2)·4+(1/2)·4]=2.
  3. Comparar

    Masa uno en v=0 tiene la misma ρ y j, pero energía cero.

Resultado

Media cero no significa reposo de cada partícula.

Lo comprobamos

La varianza de velocidad de la primera población es cuatro.

Dividir corriente por masa

Masa 3/2 en v=1 y 1/2 en v=−1.

  1. Masa

    ρ=2; no es una distribución normalizada a uno.
  2. Corriente

    j=3/2−1/2=1; u=1/2.
  3. Dispersión

    (3/2)(1−1/2)²+(1/2)(−1−1/2)²=3/2.
Resultado

j y u tienen valores distintos porque la masa local no es uno.

Lo comprobamos

El segundo momento total es dos y la energía cinética uno.

4. Integrar la ecuación da continuidad, pero no siempre cierra el sistema

Con regularidad y momentos que permitan intercambiar operaciones, integrar en v da \rho_t+\operatorname{div}_xj=0. También puede formularse débilmente usando funciones de prueba. La continuidad conserva masa local mediante un flujo.

Multiplicar antes por v introduce el segundo momento \Pi=\int v\otimes vf\,dv. En general la ecuación siguiente introduce un momento mayor. El cierre necesita una hipótesis o un teorema, no solo seguir integrando formalmente.

Dos velocidades tienen un cierre especial

En una dimensión solo se permiten v=±1, con densidades f₊ y f₋.

  1. Transportes

    ∂ₜf₊+∂ₓf₊=0 y ∂ₜf₋−∂ₓf₋=0.
  2. Sumar

    ρ=f₊+f₋ y j=f₊−f₋ dan ρₜ+jₓ=0.
  3. Restar

    jₜ+ρₓ=0 porque v²=1 en ambas poblaciones.
Resultado

El cierre sale aquí de una restricción exacta de velocidades.

Lo comprobamos

No se aplica ese mismo cierre a velocidades reales arbitrarias.

Por qué falla el cierre general

Compara las poblaciones de media cero anteriores.

  1. Datos comunes

    Ambas tienen ρ=1 y j=0.
  2. Segundo momento

    La de velocidades ±2 tiene Π=4; la quieta Π=0.
  3. Conclusión

    Una fórmula universal Π=Π(ρ,j) no puede describir ambas sin más datos.

Resultado

La pérdida de información se comprueba con dos ejemplos explícitos.

Lo comprobamos

Una aproximación de equilibrio añade información que no estaba en esos dos momentos.

5. La fuerza exige revisar la forma conservativa en velocidades

Con aceleración a(t,x) independiente de v, f_t+v\cdot\nabla_xf+a\cdot\nabla_vf=Q. Si el borde en velocidad desaparece, la fuerza no crea masa y aporta ρa al balance de impulso.

Si a depende de v, la forma conservativa es f_t+\operatorname{div}_x(vf)+\operatorname{div}_v(af)=Q. Expandir el último término añade f·divᵥa. Una fricción puede comprimir volumen de fase y cambiar el valor de f a lo largo de trayectorias.

Aceleración constante

a=2, x₀=0, v₀=1.

  1. Velocidad

    v(t)=1+2t.
  2. Posición

    x(t)=t+t².
  3. En t=2

    v=5 y x=6; usar x=x₀+tv₀ habría dado dos y sería incorrecto.

Resultado

Las características también se modifican al añadir fuerza.

Lo comprobamos

La fórmula libre conserva velocidad; esta ya no lo hace.

Compresión por fricción

En una dimensión, a=−γv con γ>0, sin colisiones.

  1. Velocidad

    v=e⁻ᵞᵗv₀.

    El jacobiano completo de fase es e⁻ᵞᵗ.

  2. Densidad

    La forma conservativa da df/dt=γf sobre la característica.
  3. Balance

    f crece por eᵞᵗ mientras el volumen se contrae por e⁻ᵞᵗ: la masa se conserva.

Resultado

Confundir a·∇ᵥf con divᵥ(af) perdería este factor.

Lo comprobamos

La fuerza y el dominio deben permitir las operaciones y el balance declarados.

6. Una colisión elástica cambia velocidades conservando ciertos momentos

El operador de Boltzmann Q compara ganancia y pérdida de partículas en una velocidad por colisiones binarias. En el marco elástico y con integrabilidad y simetrías apropiadas, \int Q\,dv=\int vQ\,dv=\int |v|^2Q\,dv=0. Son conservación de masa, impulso y energía de colisiones.

Para soluciones positivas suficientemente regulares y un sistema cerrado, el teorema H hace no creciente la integral de f·log f. La entropía física tiene el signo opuesto. Núcleo de colisión, fronteras y régimen importan; no toda redistribución arbitraria es una colisión de Boltzmann.

Dos partículas de igual masa

Velocidades iniciales v=(2,0), w=(0,0), con normal n=(1,1)/√2.

  1. Regla

    v′=v−[(v−w)·n]n y w′=w+[(v−w)·n]n.
  2. Calcular

    El producto es √2; v′=(1,−1), w′=(1,1).
  3. Comprobar

    Impulso antes y después: (2,0). Suma de módulos cuadrados: cuatro en ambos casos.

Resultado

Cambian velocidades individuales sin cambiar esos totales.

Lo comprobamos

La energía cinética añade el mismo factor de masa por un medio.

Ver el signo de una producción de H

Modelo ilustrativo de dos casillas positivas: f₁′=f₂−f₁ y f₂′=f₁−f₂.

  1. Masa

    La suma de derivadas es cero.

  2. Derivar H

    H′=(f₂−f₁)(log f₁−log f₂).
  3. Signo

    El logaritmo es creciente: los factores tienen signos opuestos y H′≤0.

Resultado

La cuenta muestra un mecanismo de relajación de entropía.

Lo comprobamos

Es un modelo de dos casillas, no una afirmación de conservar todos los momentos de Boltzmann.

7. Una Maxwelliana cierra momentos solo bajo un modelo apropiado

En d dimensiones, M(v)=\rho(2\pi\theta)^{-d/2}e^{-|v-u|^2/(2\theta)}, con ρ>0 y θ>0, tiene masa ρ, velocidad media u y dispersión térmica. En el operador elástico estándar está en equilibrio de colisiones.

Equilibrio de colisiones local no implica que todo el transporte esté en equilibrio. Usar Maxwellianas locales para cerrar momentos es una aproximación que necesita justificación. Límites hacia fluidos requieren escalas y estimaciones; escribir una gaussiana no demuestra por sí solo un límite de Navier–Stokes.

Leer sus momentos en una dimensión

ρ=2, u=1 y θ=4.

  1. Masa y corriente

    ∫M=2 y j=ρu=2.
  2. Segundo momento

    ∫v²M=ρ(u²+θ)=2·5=10.
  3. Energía

    La energía cinética es cinco; la presión térmica ρθ=8.
Resultado

La energía incluye movimiento medio y dispersión.

Lo comprobamos

θ es la varianza normalizada de velocidad, no su raíz.

Colisiones equilibradas, transporte pendiente

f(x,v)=ρ(x)·(2π)⁻¹ᐟ²e⁻ᵛ²ᐟ², independiente de t, sin fuerza.

  1. En cada x

    Es una Maxwelliana: Q=0.
  2. Transportar

    v·∂ₓf=vρ′(x)·(2π)⁻¹ᐟ²e⁻ᵛ²ᐟ².
  3. Si ρ′≠0

    Ese término no se anula para toda v: no es una solución estacionaria libre.

Resultado

Equilibrio de Q y equilibrio de toda la ecuación son condiciones distintas.

Lo comprobamos

Con una fuerza adecuada podría haber otro balance, que habría que especificar.

8. Campo medio y Boltzmann–Grad usan escalas diferentes

En campo medio, una fuerza agregada depende de la distribución, como a=K*ρ. En sistemas de N partículas aparece a menudo un factor 1/N en la suma de interacciones. Bajo hipótesis sobre fuerza y datos se estudia su aproximación por Vlasov y la propagación del caos.

Boltzmann–Grad hace crecer el número de partículas mientras disminuye su tamaño con una sección eficaz total controlada. Es otro régimen, con correlaciones y recolisiones que deben analizarse. Los teoremas especifican condiciones e intervalo temporal; no se presume validez universal de cualquier cierre para todos los tiempos.

Promediar una fuerza sencilla

En una dimensión, K(z)=−z y a_{i}=(1/N)Σ_jK(x_{i}−x_{j}).

  1. Sumar

    a_{i}=−x_{i}+(1/N)Σ_jx_j= x̄−x_{i}.
  2. Dos posiciones

    Para −1 y 1, x̄=0 y las aceleraciones son 1 y −1.
  3. Centro

    La suma de aceleraciones es cero: el centro no se acelera.

Resultado

El factor 1/N mantiene una escala de fuerza promediada.

Lo comprobamos

Este cálculo finito no prueba el teorema de límite para fuerzas arbitrarias.

Una escala geométrica en dimensión tres

Como ilustración, fija Nε²=1 para partículas de tamaño ε.

  1. Primer caso

    N=100, ε=1/10: Nε²=1 y Nε³=1/10.
  2. Segundo

    N=10000, ε=1/100: Nε²=1 y Nε³=1/100.
  3. Leer

    La escala de sección permanece mientras la de volumen ocupado disminuye.

Resultado

Este escalado no es el mismo que debilitar cada fuerza por 1/N.

Lo comprobamos

Factores geométricos, datos y evolución siguen formando parte de la derivación física.

9. Una marginal no determina las correlaciones de dos partículas

La distribución conjunta de N partículas contiene dependencias. La marginal de una sola no determina la de dos. La jerarquía BBGKY refleja que la evolución de una marginal puede depender de la siguiente; sustituir una distribución conjunta por un producto es una afirmación adicional.

Propagación del caos significa, en modelos y límites apropiados, que marginales de un número fijo de partículas se aproximan a productos a partir de datos caóticos. No afirma independencia exacta de todas las partículas interactuantes en todo tiempo finito.

Tres parejas con las mismas marginales

Cada velocidad individual vale ±1 con probabilidades 1/2.

  1. Independientes

    Las cuatro parejas tienen probabilidad 1/4; E[v₁v₂]=0.
  2. Correlacionadas

    Solo ++ y −−, cada una con 1/2; E[v₁v₂]=1.
  3. Opuestas

    Solo +− y −+, cada una con 1/2; E[v₁v₂]=−1.
Resultado

La misma marginal permite correlaciones completamente diferentes.

Lo comprobamos

No se puede reconstruir la tabla conjunta multiplicando sin hipótesis de independencia.

La suma revela la dependencia

S=v₁+v₂ en esas tres distribuciones.

  1. Media

    E[S]=0 en las tres.
  2. Segundo momento

    E[S²]=E[v₁²]+E[v₂²]+2E[v₁v₂].
  3. Resultados

    Son dos, cuatro y cero, respectivamente.

Resultado

Las correlaciones afectan a fluctuaciones aunque no cambien la media individual.

Lo comprobamos

Un límite que suprima correlaciones debe justificar esa supresión.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Seguir una partícula

Con x₀=2,v=3, calcula x en t=4.

  1. Datos

    x₀=2, v=3 y t=4.
  2. Modelo

    Transporte libre sin aceleración.

  3. Característica

    x(t)=x₀+tv.
  4. Producto

    tv=4·3=12.
  5. Sumar

    x=2+12=14.
  6. Comprobar

    La velocidad sigue siendo tres: no se ha añadido fuerza.

Resultado

La posición es 14.

Lo comprobamos

La velocidad sigue siendo 3.

EJEMPLO RESUELTO 2

Dos velocidades discretas

Una distribución tiene masas 1/4 en v=−1 y 3/4 en v=1. Calcula velocidad media.

  1. Distribución

    Masas 1/4 en menos uno y 3/4 en uno.

  2. Normalizar

    La masa total es uno, así que la corriente coincide con la media.

  3. Aporte negativo

    (1/4)(−1)=−1/4.
  4. Aporte positivo

    (3/4)(1)=3/4.
  5. Sumar

    u=−1/4+3/4=1/2.
  6. Comprobar

    La media está entre menos uno y uno; no dice que cada partícula viaje a un medio.

Resultado

La velocidad media es 1/2.

Lo comprobamos

La masa total es uno.

Los símbolos que hemos usado

f(t,x,v)≥0
Densidad en espacio de fases, normalizada como probabilidad o cantidad según el modelo.
ρ=∫f dv
Densidad espacial.
j=∫vf dv
Flujo o primer momento de velocidad.
u=j/ρ
Velocidad media local, donde ρ>0.
Q(f,f)
Operador de colisiones del modelo de Boltzmann.
Marginal
Distribución de parte de las partículas obtenida integrando las restantes.
Propagación del caos
Convergencia de marginales finitas hacia productos en un régimen especificado.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar