Hear that, air?

An ideal gas of non-interacting helium atoms is in a thermal equilibrium at the temperature T = 3 0 C T = 30^\circ \mbox{C} . What's the probability that an atom in this gas has a speed near the speed of sound, i.e. has a speed in the range 340 ± 5 m/s 340 \pm 5~\mbox{m/s} ?

Details and assumptions

  • Mass of a helium atom is 6.646 × 1 0 27 kg 6.646 \times 10^{-27}~\mbox{kg}
  • The probability density as a function of energy for an atom in an ideal gas is p ( E ) = 2 k T π E k T exp ( E k T ) p(E) = \frac{2}{kT\sqrt{\pi}} \sqrt{\frac{E}{kT}} \exp({-\frac{E}{kT}}) .


The answer is 0.001683.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

10 solutions

Noah Fang
Jul 28, 2013

First, it is important to note two things:

  1. Ideal gas has no internal interactions between particles. This means that all of the energy in the particle is kinetic energy.

    1. exp(w) = e w e^{w}

\hspace{10mm}

Now, we will focus on energy range. Since the target speed is anywhere from 335 335 m/s to 345 345 m/s, we will have energy levels from

\hspace{5mm}

1 2 ( 6.646 × 1 0 27 ) ( 335 ) 2 3.729 × 1 0 22 J \frac{1}{2}(6.646 \times 10^{-27})(335)^{2} \approx 3.729 \times 10^{-22} J

to

1 2 ( 6.646 × 1 0 27 ) ( 345 ) 2 3.955 × 1 0 22 J \frac{1}{2}(6.646 \times 10^{-27})(345)^{2} \approx 3.955 \times 10^{-22} J

\hspace{10mm}

Next, we can integrate the probability density function, p ( E ) p(E) , from the lower energy level to the higher energy level.

\hspace{5mm}

3.729 × 1 0 22 3.955 × 1 0 22 ( 2 k T π × E k T × e E k T ) d E \int^{3.955 \times 10^{-22}}_{3.729 \times 10^{-22}} (\frac{2}{kT\sqrt{\pi}} \times \sqrt{\frac{E}{kT}} \times e^{-\frac{E}{kT}})\,dE

\hspace{10mm}

Solving this integral, we find the probability to be 0.001683 \fbox{0.001683}

\hspace{10mm}

Relevant Quantities

k k is the Boltzmann constant or 1.38 × 1 0 23 1.38 \times 10^{-23} m 2 m^{2} k g kg s 2 s^{-2} K 1 K^{-1}

T T is the temperature in Kelvins . In this case, it is 303.1 5 K 303.15 ^\circ K

Louie Tan Yi Jie
Jul 29, 2013

Energies in J corresponding to 335 m/s and 345 m/s are:

E 1 = 1 2 m ( 33 5 2 ) = 3.72924 × 1 0 22 E 2 = 1 2 m ( 34 5 2 ) = 3.95520075 × 1 0 22 E_1= \frac{1}{2} m (335^2) = 3.72924 \times 10^{-22}\\ E_2= \frac{1}{2} m (345^2) = 3.95520075 \times 10^{-22}

Using k = 1.38 × 1 0 23 k=1.38 \times 10^{-23} , T = 273.15 + 30 T=273.15+30 , m = 6.646 × 1 0 27 m=6.646 \times 10^{-27} , required probability is:

E 1 E 2 p ( E ) d E = 0.00168489 \int_{E_1}^{E_2} p(E) \, dE=0.00168489

Suhail Sherif
Jul 29, 2013

# from math import * p = lambda E: (2.0/((1.380649 pow(10,-23)) 303 sqrt(pi))) * sqrt(E/((1.380649 pow(10,-23)) 303)) * pow(e,(-E/((1.380649 pow(10,-23)) 303))) ll = 3.323 pow(10,-27) 335 335 ul = 3.323 pow(10,-27) 345 345 def integrate(f, l, u, dx): ans = 0.0 x=l while x<u: ans += f(x) x += dx return ans dx integrate(p,ll,ul,pow(10,-29))

I wasn't sure whether this would work. I tried finding the integral mathematically but couldn't get anywhere. So I thought I'd try it some other day and just wrote the above out of curiosity to see if it'd work. It's not very accurate, but it works for upto 5 decimal places.

Indefinite integral (in exact terms) is e k T ( k T erf ( e k T ) e 2 e e k T π ) \sqrt{\frac{e}{k T}} \left(\frac{\sqrt{k} \sqrt{T} \text{erf}\left(\frac{\sqrt{e}}{\sqrt{k} \sqrt{T}}\right)}{\sqrt{e}}-\frac{2 e^{-\frac{e}{k T}}}{\sqrt{\pi }}\right)

But I guess we only need a 3sf answer?

Louie Tan Yi Jie - 7 years, 10 months ago
Aaron Aquino
Jul 29, 2013

The probability can be determined by evaluating the integral of the probability density function p ( E ) p(E) with lower and upper limits E 1 E_1 and E 2 E_2 , respectively.

E 1 = 1 2 m v 1 2 = 1 2 ( 6.646 × 1 0 27 k g ) ( 335 m / s ) 2 = 3.7292 × 1 0 22 J . E_1 = \frac12mv_1^2 = \frac12(6.646 \times 10^{-27} kg)(335 m/s)^2 = 3.7292 \times 10^{-22} J. E 2 = 1 2 m v 2 2 = 1 2 ( 6.646 × 1 0 27 k g ) ( 345 m / s ) 2 = 3.9552 × 1 0 22 J . E_2 = \frac12mv_2^2 = \frac12(6.646 \times 10^{-27} kg)(345 m/s)^2 = 3.9552 \times 10^{-22} J. Evaluating the integral, we get E 1 E 2 p ( E ) d E = 0.001683. \int\limits_{E_1}^{E_2} p(E) dE = 0.001683.

Enric Boix
Aug 3, 2013

E = 1 2 m v 2 E = \frac{1}{2} m v^2

Minimum energy is when v = 335 m / s v = 335 m/s : E = 3.72923675 1 0 22 E = 3.72923675*10^-22

Maximum energy is when v = 345 m / s v = 345 m/s : E = 3.95520075 1 0 22 E = 3.95520075*10^-22

Integrate probability distribution function from minimum energy to maximum energy. I did it numerically.

Francisco Rivera
Jul 31, 2013

The lowest possible energy that the atom can have to be in our acceptable range happens with the lowest possible speed and equals E min = 1 2 m v min 2 = 3.729 1 0 22 J E_{\text{min}} = \frac{1}{2} mv_{\text{min}}^2 = 3.729 \cdot 10^{-22} \text{ J}

Similarly the maximum happens with the maximum speed E max = 1 2 m v max 2 = 3.955 1 0 22 J E_{\text{max}} = \frac{1}{2} m v_{\text{max}}^2 = 3.955 \cdot 10^{-22} \text{ J}

Now that we have the energy bounds, we must simply integrate the PDF function in the range to find the probability that the energy (and by extension the speed) is in the desired range: p = E min E max 2 k T π E k T exp ( E k t ) d E p = \int_{E_{\text{min}}}^{E_{\text{max}}} \frac{2}{kT\sqrt{\pi}} \sqrt{\frac{E}{kT}} \exp\left( -\frac{E}{kt} \right) \, dE where k = 1.38 1 0 23 k=1.38\cdot 10^{-23} J \cdot K 1 ^{-1} and T = 273 + 30 T = 273+30 K.

Evaluating gives us that p 0.0001683 p \approx \boxed{0.0001683} .

Vasileios Tsiolis
Jul 31, 2013

According to kinetic gas theory and since we are dealing with an ideal gas, the total energy for an atom equals its kinetic energy. Therefore E = ( m v 2 ) / 2 E=(mv^{2})/2 , where m m is the mass of a helium atom. The probability that an atom has a speed near the speed of sound, which means a speed between v 1 = 335 m / s v_1=335m/s and v 2 = 345 m / s v_2=345m/s is the probability that the atom's energy is between E 1 = ( m v 1 2 ) / 2 = 3.729 × 1 0 22 J E_1=(mv_1^{2})/2=3.729\times10^{-22}J and E 2 = ( m v 2 2 ) / 2 = 3.955 × 1 0 22 J E_2=(mv_2^{2})/2=3.955\times10^{-22}J . Given that the probability as a function of energy for an atom in an ideal gas is p ( E ) p(E) (where p ( E ) p(E) is the function given in details and assumptions) the desired probability is the integral of p ( E ) p(E) from E 1 E_1 to E 2 E_2 . By doing so we find that the probability is 0.00168 0.00168 , keeping 5 decimal digits. And now just a few notes:
The temperature we substitute in the p ( E ) p(E) formula needs to be in S.I. units so T = 30 + 273.15 = 303.15 K T=30+273.15=303.15K . That k k in p ( E ) p(E) formula is Boltzmann's constant and it is equal to 1.3806488 × 1 0 23 m 2 k g ( s ) 2 K 1 1.3806488\times10^{-23}m^{2}kg(s)^{-2}K^{-1} .

Daniel Chiu
Jul 30, 2013

The energy of an atom with speed 335 is a = 1 2 6.646 1 0 27 33 5 2 a=\dfrac{1}{2}6.646\cdot10^{-27}\cdot 335^2 The energy of an atom with speed 345 is b = 1 2 6.646 1 0 27 34 5 2 b=\dfrac{1}{2}6.646\cdot10^{-27}\cdot 345^2 If we integrate p ( E ) p(E) from a a to b b we get our answer, 0.001683 \boxed{0.001683} .

Ricky Escobar
Jul 30, 2013

To find the probability that that the particle has an energy in a certain range, we simply integrate the probability density function over that range.The probability density function is a function of energy, so we must find the kinetic energies associated with the given velocities v 1 = 335 m / s , v 2 = 345 m / s \ v_1=335 \ \mathrm{m/s}, \ v_2=345 \ \mathrm{m/s} .

$$K 1 = \frac{1}{2}m v 1^2=3.746 \times 10^{-22} \ \mathrm{J}$$

$$K 2 = \frac{1}{2}m v 2^2=3.973 \times 10^{-22} \ \mathrm{J}$$

So then the probability P P is

$$P=\int {K 1}^{K_2} {\frac{2}{kT\sqrt{\pi}}\sqrt{\frac{E}{kT}}\exp{\left(-\frac{E}{kT}\right)}} \ dE$$

which (taking k = 1.38 × 1 0 23 J / K k=1.38 \times 10^{-23} \ \mathrm{J/K} and T = 303 K T=303 \ \mathrm{K} ) can be evaluated numerically to be

$$ \fbox{0.001697}. $$

The probability for an atom to have a speed in a given range is P ( v 1 , v 2 ) = v 1 v 2 p ( v ) d v P(v_1,v_2)=\int_{v_1}^{v_2}p(v)\,\mathrm{d}v . As the speed can be bijectively mapped to the kinetic energy, we can also calculate the sought probability as P ( v 1 , v 2 ) = P ( E 1 , E 2 ) = E 1 E 2 p ( E ) d E P(v_1,v_2)=P(E_1,E_2)=\int_{E_1}^{E_2}p(E)\,\mathrm{d}E .

On the other hand, energies are expressed through speeds as E = 1 2 m v 2 E=\frac{1}{2}mv^2 . Thus P ( E 1 , E 2 ) = 2 k T π E k T E 1 E 2 exp ( E / k T ) d E P(E_1,E_2)=\frac{2}{kT\sqrt{\pi}}\sqrt{\frac{E}{kT}}\int_{E_1}^{E_2}{\exp(-E/kT)}\,\mathrm{d}E , which, substituting the numbers, gives us P ( 335 , 345 ) = 0.168 % P(335,345)=0.168\% .

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...