===== Boltzmann equation ===== ==== Set ==== | @#55CCEE: context | @#55CCEE: $ {\bf K}:\mathbb R^3\times\mathbb R^3\times\mathbb t\to\mathbb R $ | | @#DDDDDD: range | @#DDDDDD: $ :: {\bf K}({\bf x},{\bf v},t)$ | | @#55CCEE: context | @#55CCEE: $ S\in\mathbb N $ | | $i,j\in\text{range}(S)$ | | @#55CCEE: context | @#55CCEE: $ m_i\in \mathbb R^* $ | | @#55CCEE: context | @#55CCEE: $ I_{ij}: \mathbb R^3\times[-\tfrac{\pi}{2},\tfrac{\pi}{2}]\times[0,2\pi]\to\mathbb R $ | | @#DDDDDD: range | @#DDDDDD: $ :: I_{ij}({\bf v},\vartheta,\varphi) $ | The integer $S$ denote the species, $m_i$ their respective masses and $I_{ij}$ the differential cross sections for two particle collisions. | @#55CCEE: context | @#55CCEE: $ {\bf v'},\ {\bf v'}_1 : \mathbb R^{2\times 3}\to \mathbb R^3 $ | | @#DDDDDD: range | @#DDDDDD: $ :: {\bf v'}({\bf v'},{\bf v}_1),\ {\bf V}_1({\bf v},{\bf v}_1) $ | | @#FFBB00: definiendum | @#FFBB00: $ {\bf f} \in \mathrm{it} $ | | @#55EE55: postulate | @#55EE55: $ f_i:\mathbb R^3\times\mathbb R^3\times\mathbb R\to\mathbb R_+ $ | | @#DDDDDD: range | @#DDDDDD: $ :: f_i({\bf x},{\bf v},t)$ | | @#DDDDDD: range | @#DDDDDD: $ J[f_i|f_j]({\bf x},{\bf v},t) \equiv \int\int g\ I_{ij}(g,\vartheta,\varphi)\ \left(f_i({\bf x},{\bf v'}({\bf v},{\bf v}_1),t)\cdot f_j({\bf x},{\bf v'}_1({\bf v},{\bf v}_1),t)-f_i({\bf x},{\bf v},t)\cdot f_j({\bf x},{\bf v}_1,t)\right)\ \mathrm d\Omega(\vartheta,\varphi)\ \mathrm d^3v_1 $ | | @#55EE55: postulate | @#55EE55: $ \left(\frac{\mathrm \partial}{\partial t}+{\bf v}\cdot\nabla_{\bf x}+\frac{1}{m_i}{\bf K}\cdot\nabla_{\bf v}\right)f_i = \sum_{j=1}^S J[f_i|f_j]$ | >I'm not sure about the summation "$\cdots = \sum_{j=1}^S J[f_i|f_j]$" here --- check that. ==== Discussion ==== There are different routes to obtain the Boltzmann equation as approximation to the one-particle distribution function equation in the [[BBGKY]] hierarchy. In that case one also has $ \int f_i({\bf x},{\bf v},t)\ \mathrm d^3x \mathrm d^3v \overset{!}{=} 1 $ Neglecting some arguments, we can read the "collision integral functional" $J$ as: $ J[f_i|f_j] \equiv \int\int g\ I_{ij}(g,\vartheta,\varphi)\ \left(f_i({\bf v'})\cdot f_j({\bf v'}_1)-f_i({\bf v})\cdot f_j({\bf v}_1)\right)\ \mathrm d\Omega(\vartheta,\varphi)\ \mathrm d^3v_1. $ The velocity-after-collision-${\bf V}$-independent term contains a factor $f$ which you can pull out, leaving an integral which is a velocity average over a collision rate $ -\ \left\langle \int g\ I_{ij}\ \mathrm d\Omega \right\rangle_{f_j({\bf v}_1)}\cdot f_i({\bf v}). $ The function $f_j$ arises as one-particle $({\bf x}_1,{\bf v}_1)$-density and so if $g\int I\ \mathrm d\Omega\equiv g\ \sigma_\text{tot}$ is volume transfer per time, then integrating out ${\bf v}$ makes the terms in bracket into a process characteristic time $\frac{1}{\tau}$. Hence the term is of the form of a decay rate $r:=-\frac{1}{\tau} f.$ at the fixed point velocity ${\bf v}$. The second term is more involved, since it doesn't represent the loss at ${\bf v}$, but the gain: It's the sum of processes of pairs particles with velocities ${\bf v'},{\bf v'}_1$, which end up with particles having the velocity ${\bf v}$ and any other velocity ${\bf v}_1$. === Note === $Assumptions = {kT > 0, n > -1, \[CapitalTheta] > 0}; f[EE_] = E^(-(EE/kT))* E^(-(EE^2/\[CapitalTheta]^2))/(E^(\[CapitalTheta]^2/(2 kT)^2) Sqrt[\[Pi]] /2 \[CapitalTheta] Erfc[\[CapitalTheta]/(2 kT)]); Integrate[f[EE] EE^n, {EE, 0, \[Infinity]}] === Reference === Wikipedia: [[https://en.wikipedia.org/wiki/Boltzmann_equation|Boltzmann equation]] ==== Parents ==== === Subset of === [[PDE system]] === Related === [[BBGKY hierarchy]]