Editing Maxwell speed distribution
Jump to navigation
Jump to search
The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.
Latest revision | Your text | ||
Line 1: | Line 1: | ||
The '''Maxwellian | The '''Maxwellian velocity distribution''' <ref> [http://books.google.com/books?id=hYIBOMxIuvEC&source=gbs_similarbooks_r&cad=2 James C. Maxwell, "The scientific papers of James Clerk Maxwell", Edited by W.D. Niven, paper number XX, Dover Publications, Vol. I,II, New York, USA (2003)]</ref> provides probability that the speed of a molecule of mass ''m'' lies in the range ''v'' to ''v+dv'' is given by | ||
:<math>P(v)dv = 4 \pi v^2 dv \left( \frac{m}{2 \pi k_B T} \right)^{3/2} \exp (-mv^2/2k_B T) </math> | :<math>P(v)dv = 4 \pi v^2 dv \left( \frac{m}{2 \pi k_B T} \right)^{3/2} \exp (-mv^2/2k_B T) </math> | ||
Line 14: | Line 14: | ||
and the root-mean-square speed by | and the root-mean-square speed by | ||
:<math>\sqrt{\overline{v^2}} = \sqrt | :<math>\sqrt{\overline{v^2}} = \sqrt \frac{3}{2} v_{\rm max}</math> | ||
==Derivation== | ==Derivation== | ||
According to the '''Shivanian and Lopez-Ruiz model''' <ref>[http:// | According to the '''Shivanian and Lopez-Ruiz model''' <ref>[http://arxiv.org/abs/1105.4813 Elyas Shivanian and Ricardo Lopez-Ruiz "A New Model for Ideal Gases. Decay to the Maxwellian Distribution", arXiv:1105.4813v1 24 May (2011)]</ref>, consider an [[ideal gas]] composed particles having a mass of unity in the three-dimensional (<math>3D</math>) space. As long as there no privileged direction when in equilibrium, we can take any direction in space and study the discrete time evolution of the velocity distribution in that direction. Let us call this axis <math>U</math>. We can complete a Cartesian system with two additional orthogonal axis <math>V,W</math>. If <math>p_n(u){\mathrm d}u</math> represents the probability of finding a particle of the gas with velocity component in the direction <math>U</math> comprised between <math>u</math> and <math>u + {\mathrm d}u</math> at time <math>n</math>, then the probability to have at this time <math>n</math> a particle with a <math>3D</math> velocity <math>(u,v,w)</math> will be <math>p_n(u)p_n(v)p_n(w)</math>. The particles of the gas collide between them, and after a number of interactions of the order of system size, a new velocity distribution is attained at time <math>n+1</math>. Concerning the interaction of particles with the bulk of the gas, we make two simplistic and realistic assumptions in order to obtain the probability of having a velocity <math>x</math> in the direction <math>U</math> at time <math>n+1</math>: (1) Only those particles with an energy greater than <math>x^2</math> at time <math>n</math> can contribute to this velocity <math>x</math> in the direction <math>U</math>, that is, all those particles whose velocities <math>(u,v,w)</math> verify <math> u^2+v^2+w^2\ge x^2</math>; (2) The new velocities after collisions are equally | ||
distributed in their permitted ranges, that is, particles with velocity <math>(u,v,w)</math> can generate maximal velocities <math>\pm U_{max}=\pm\sqrt{u^2+v^2+w^2}</math>, then the allowed range of velocities <math>[-U_{max},U_{max}]</math> measures <math>2|U_{max}|</math>, and the contributing probability of these particles to the velocity <math>x</math> will be <math>p_n(u)p_n(v)p_n(w)/(2|U_{max}|)</math>. Taking all together we finally get the expression | |||
for the evolution operator <math>T</math>. This is: | |||
:<math> | :<math> | ||
p_{n+1}(x)= | p_{n+1}(x)=Tp_n(x) = \iiint_{u^2+v^2+w^2\ge x^2}\,{p_n(u)p_n(v)p_n(w)\over 2\sqrt{u^2+v^2+w^2}} \; {\mathrm d}u~{\mathrm d}v~{\mathrm d}w\,. | ||
</math> | </math> | ||
Let us remark that we have not made any supposition about the type of interactions or collisions | Let us remark that we have not made any supposition about the type of interactions or collisions | ||
between the particles and, in some way, the equivalent of the Boltzmann hypothesis of ''molecular chaos'' would be the two simplistic assumptions we have stated on the interaction of particles with the bulk of the gas. In fact, the operator <math> | between the particles and, in some way, the equivalent of the Boltzmann hypothesis of ''molecular chaos'' would be the two simplistic assumptions we have stated on the interaction of particles with | ||
the bulk of the gas. In fact, the operator <math>T</math> conserves the energy and the null momentum of the gas over time. Moreover, for any initial velocity distribution, the system tends towards its equilibrium, i.e. towards the Maxwellian Velocity Distribution (MVD). This means that | |||
:<math> | :<math> | ||
\lim_{n\rightarrow\infty} | \lim_{n\rightarrow\infty} T^n \left(p_0(x)\right) \rightarrow p_f(x)= \mathrm{MVD}\;(1D\;case)\,. | ||
</math> | </math> | ||
Line 45: | Line 47: | ||
===Theorem 1=== | ===Theorem 1=== | ||
For any <math>p</math> with <math>||p||=1</math>, we have <math>|| | For any <math>p</math> with <math>||p||=1</math>, we have <math>||Tp||=||p||</math>. | ||
This can be interpreted as the conservation of the number of particles, or in an equivalent way, the total mass of the gas. | This can be interpreted as the conservation of the number of particles, or in an equivalent way, the total mass of the gas. | ||
Line 51: | Line 53: | ||
===Theorem 2=== | ===Theorem 2=== | ||
The mean value of the velocity in the recursion <math>p_n= | The mean value of the velocity in the recursion <math>p_n=T^np_0</math> is conserved in time. | ||
In fact, it is null for all <math>n</math>: | In fact, it is null for all <math>n</math>: | ||
:<math> | :<math> | ||
\langle x, | \langle x,Tp \rangle = \langle x,T^2p \rangle = \langle x,T^3p \rangle=\cdots= \langle x,T^np \rangle =\cdots=0\,, | ||
</math> | </math> | ||
Line 64: | Line 66: | ||
</math> | </math> | ||
It means that the zero total momentum of the gas is conserved in its time evolution under the action of <math> | It means that the zero total momentum of the gas is conserved in its time evolution under the action of <math>T</math>. | ||
===Theorem 3=== | ===Theorem 3=== | ||
Line 71: | Line 73: | ||
:<math> | :<math> | ||
\langle x^2,p \rangle= \langle x^2, | \langle x^2,p \rangle= \langle x^2,Tp \rangle= \langle x^2,T^2p \rangle= \langle x^2,T^3p \rangle =\cdots= \langle x^2,T^np \rangle=\cdots \,. | ||
</math> | </math> | ||
Line 78: | Line 80: | ||
===Theorem 4=== | ===Theorem 4=== | ||
The one-parametric family of normalized Gaussian functions <math>p_{\alpha}(x)=\sqrt{\alpha\over\pi}e^{-\alpha x^2}</math>, <math>\alpha\ge 0</math>, <math>||p_{\alpha}||=1</math>, are fixed points of the operator <math> | The one-parametric family of normalized Gaussian functions <math>p_{\alpha}(x)=\sqrt{\alpha\over\pi}e^{-\alpha x^2}</math>, <math>\alpha\ge 0</math>, <math>||p_{\alpha}||=1</math>, are fixed points of the operator <math>T</math>. In other words, <math>Tp_{\alpha}=p_{\alpha}</math>. | ||
===Conjecture=== | ===Conjecture=== | ||
Line 88: | Line 86: | ||
As a consequence of the former theorems, and by simulation of many examples, the following conjecture can be stated: | As a consequence of the former theorems, and by simulation of many examples, the following conjecture can be stated: | ||
For any <math>p</math> with <math>||p||=1</math>, with finite <math> \langle x^2,p \rangle </math> and verifying <math>\lim_{n\rightarrow\infty} || | For any <math>p</math> with <math>||p||=1</math>, with finite <math> \langle x^2,p \rangle </math> and verifying <math>\lim_{n\rightarrow\infty} ||T^np(x)-\mu(x)||=0</math>, the limit <math>\mu(x)</math> is the fixed point <math>p_{\alpha}(x)=\sqrt{\alpha\over\pi}\,e^{-\alpha x^2}</math>, with <math>\alpha=(2\, \langle x^2,p \rangle)^{-1}</math>. | ||
===Conclusion=== | ===Conclusion=== | ||
In physical terms, it means that for any initial velocity distribution of the gas, it decays to the Maxwellian distribution, which is just the fixed point of the dynamics. Recalling that <math> \langle x^2,p \rangle= | In physical terms, it means that for any initial velocity distribution of the gas, it decays to the Maxwellian distribution, which is just the fixed point of the dynamics. Recalling that <math> \langle x^2,p \rangle=k\tau</math>, with <math>k</math> the Boltzmann constant and <math>\tau</math> the temperature of the gas, and introducing the mass <math>m</math> of the particles, let us observe that the MVD (above presented) is recovered in its <math>3D</math> format: | ||
:<math> | :<math> | ||
\mathrm{MVD} = p_{\alpha}(u)p_{\alpha}(v)p_{\alpha}(w)=\left({m\alpha\over\pi}\right)^{3\over 2}\,\exp^{-m\alpha (u^2+v^2+w^2)} \;\;\; with \;\;\; \alpha=( | \mathrm{MVD} = p_{\alpha}(u)p_{\alpha}(v)p_{\alpha}(w)=\left({m\alpha\over\pi}\right)^{3\over 2}\,\exp^{-m\alpha (u^2+v^2+w^2)} \;\;\; with \;\;\; \alpha=(2k\tau)^{-1}. | ||
</math> | </math> | ||