|
Expansion of an Ideal Gas Uri Lachish, guma science Contents:
1. Introduction 2. Elastic collisions
Figure-1: Elastic collision between two particles of mass m and M, moving toward each other along a straight line. The velocities of the particles before the collisions are v and
vM, and after the collision they are v' and
vM'. m v + M vM = m v' + M vM' (1) and the law of conservation of energy states that: m v2+ M vM2 = m v'2 + M vM'2 (2) Division of equation (2) by equation (1) yields: v - vM = -(v' - vM') (3) That is, the relative velocity between the particles is conserved during the collision, and the direction is
reversed. (1 + m/M) vM' - (1 - m/M) vM = (m/M) 2v (4)(1 + m/M) v' + (1 - m/M) v = 2vM (5) Assume that m is a microscopic particle and M is a macroscopic body. In the limit (m/M) --> 0: vM' = vM (6)v' = -(v - 2vM) (7)The direction of the microscopic particle is reversed and its velocity changes by 2vM, two times the velocity of the macroscopic body. The velocity of the macroscopic body is almost not affected by the collision. If the velocity of the macroscopic body is extremely slower than that of the microscopic particle, vM << v, then the kinetic energy, ΔEk, transferred from body to body will be, by equations (6) - (7): Δ Ek = (m/2) (v'2 - v2) = 2m vM (v + vM) (8)2m vM v, the dominant term on the righthand side of equation (8), is a reversible energy transfer between the microscopic and macroscopic bodies. This energy is proportional to the velocity of the macroscopic body and is independent of its mass. In a front collision, when the two bodies move toward each other, the energy will transfer from the macroscopic to the microscopic body, and when they move at the same direction, the microscopic body will transfer energy to the macroscopic body. The other term, 2m vM2, is an irreversible transformation of kinetic energy into heat, for both cases. This term becomes important in processes of drag within an ideal gas (section 4). 3. The gas pressure
Figure-2: Elastic collision of a perticle with a wall. It will collide with the wall during a time Δt, if its distance from the wall is less than |v| Cosφ Δt. The number ΔNcol of particle collisions with a unit area of the wall during Δt is: Δ Ncol = [n ∫f(v) |v| Cosφ d3v] Δt (9)n is the particle density, f(v) is the probability that the particle velocity is v, normalized by ∫4π v2 f(v) dv = 1, and d3v = 4π v2 dv Sinφ dφ . Substituting in equation (9): ΔNcol = n (∫4π v2 f(v) v dv) (∫Cosφ Sinφ dq) Δt (10) The three-dimensional average velocity <v> of particles that collide with the wall is ΔNcol = n (<v>/2) Δt (11) Upon a collision with the wall the particle will transfer a momentum 2mv Cosφ to the wall. The overall momentum ΔP, delivered to a unit area of the wall during Δt, is the sum of momenta delivered by all the colliding particles that come from all directions of the half sphere. Similarly to equations (9) - (10): ΔP = n 2m (∫4π v2 f(v) v2 dv) (∫Cos2q Sinφ dφ) Δt (12) The three-dimensional average is
<v2> = ∫4π
v2 f(v) v2 dv, p = n m <v2>/3 (13) The thermal average is calculated by applying Maxwell-Boltzmann distribution: f(v) = (m / 2πkT)3 / 2 exp(-mv2 / 2kT) (14) to the calculation of <v> and <v2> (see appendix): <v> = 4π ∫f(v) v3dv = 2 (2kT / πm)1/2 (15) <v2> = 4π ∫f(v) v4 dv = 3kT / m (16) Substituting <v2> in equation (13) yields the pressure: p = n kT (17) Equation (17) is the state equation of an ideal gas. 4. Drag within an ideal gas A particle moves with a velocity vin
toward a moving disc at a velocity vM (figure-3 ).
The particle's velocity is first transformed from a static reference to coordinates that move together with the disc (a).
Then the particle is mirror reflected by an angle
φ in the moving coordinates (b).
Finally, the particle's velocity is transformed back to the static reference (c).
If the disc moves upward, an amount of 2vM will be added to
the vertical component of the reflected velocity, according to equation (7).
The vectorial change of a particle's velocity, will be
Figure-3: Collision of a particle with a moving disc: Particles that collide with the front surface gain a kinetic energy: Δ+ =
(m/2) ((vin Cos(β) + 2vM)2 +
vin2 Sin2(β)
- vin2) and particles that collide with the back surface lose a kinetic energy: Δ- =
(m/2) ((2vM - vin Cos(β))2 +
vin2 Sin2(β) - vin2) So that the average energy given to two particles by the disc is: Δ+ + Δ- = (m/2) 8vM2 (20) The overall energy given by the disc to the particles during a time interval Δt is
ΔE =
(Δ+
+ Δ-) A
ΔNcol,
where, by (11) and (15), ΔE = (m/2) 8vM2 A n (2kT/πm)1/2 Δt (21) The work of pushing the disc a distance Δx through the gas is pex A Δx = pex A vM Δt, where pex is the pressure (force per unit area) that pushes the disc. This work is equal to the energy given to the particles, equation (21), therefore: pex A vM Δt = 4vM2 A n (2mkT/π)1/2 Δt (22) Defining the drag coefficient γ by
pex = γ
vM, its value will be:
γ
= 4n (2mkT/π)1/2 The faster particles, after colliding with the disc's front surface, and the slower particles, after colliding with its back surface, will thermalize with the other particles in the gas, and the drag energy will be transformed into heat. Thus, the isothermal movement of any surface within a gas, including pistons, is an irreversible thermodynamic process. It becomes reversible only in the limit of slow motion. b. Drag of a tilted disc A particle with a velocity vin moves at an angle
β toward a tilted moving disc.
The disc's velocity is vM
Figure-4: Collision of a particle with a tilted moving disc: Similarly to equation (18)-(20), particles that collide with the front surface gain a kinetic energy: Δ+ =
(m/2) ((vin Cos(β) + 2vM
Cos(α))2 +
vin2 Sin2(β) -
vin2) Similarly, particles that collide with the back surface lose a kinetic energy: Δ- =
(m/2) ((2vM
Cos(α) - vin Cos(β))2 +
vin2 Sin2(β) -
vin2) So that the average energy given to two particles by the disc is: Δ+ + Δ- = (m/2) 8vM2 Cos2(α) (26) so that, similarly to equations (21)-(23), ΔE = Cos2(α) 4vM2 A n (2mkT/π)1/2 Δt (27) pex A vM Δt = Cos2(α) 4vM2 A n (2mkT/π)1/2 Δt (28) and by pex = γ vM, the drag coefficient is: γ = 4n (2mkT/π)1/2 Cos2(α) (29) c. Drag of a sphere
Figure-5: The sphere parameters and coordinates. The upper figure is a front look and the lower figure is a look from above. The sphere moves upward. R is the sphere radius, φ is the elevation angle observed
in the front look, and θ is the rotation angle observed from
above, and r = R sin(φ). dF = pex r dπ
R Cos2(φ) dφ The integration over the azimuth angle θ is 2π,
and the integration over the elevation angle φ, between 0
and π/2 is 1/3. Therefore, the drag force is: F = pex A = γ vM A = (8n/3) (2mkT/π)1/2 vM (πR2) (31) where A = πR2 is the sphere cross section.
The drag coefficient is: γ = (8n/3) (2mkT/π)1/2 (32) 5. Expansion ΔE = vM n 2m (∫4π v2 f(v) v2 dv) (∫Cos2q Sinq dq) Δt (33) By using equation (12): ΔE = p vM Δt = p ΔV (34) V Δt is the distance traveled by the partition and it is also the change of the volume ΔV of the box per unit area of the moving partition.The calculation demonstrates that elastic particle collisions with a moving partition is the transformation mechanism of microscopic thermal energy of the particles into macroscopic mechanical work of the partition. 6. Isothermal expansion Δ E = N kT (ΔV/V) (35)n = N/V , V is the box volume and N is the number of particles in it. The energy is integrated in an isothermal system:∫dE = N kT ∫dV/V (36) 7. Adiabatic expansion of a mono-atomic gas A single particle collides ΔNcol / N = (v/2) Δt / V times, equation (11), with the partition during Δt. During each collision the particle loses (or gains) a velocity 2vM, equation (7), twice the partition velocity in the direction of its movement. A partition that moves in the x-direction will change the vx component of the particle velocity. However, the velocity change will eventually spread to the vy and vz components. Therefore, the change of a particle velocity during Δt will be: dv = -(1/3) 2vM (v/2) Δt / V = -(v/3) dV/V (37)since vMΔt = dV. The equation of adiabatic expansion of a mono-atomic gas, without internal degrees of freedom of gas particles, is calculated by integration of equation (37): v/v0 = (V0/V)1/3 (38)and by equation (13): p/p0 = (V0/V)5/3 (39)or, by the state equation (16), T/T0 = (V0/V)2/3, or, p/p0 = (T/T0)5/2. 8. Expansion of a quantum well ψ (x, t) = (i/2L) (e-i(kx + wt) - e+i(kx - wt)) (40)where k = π/L and w = hπ2 / 2mL2 . If the boundary at L is moving slowly with a velocity vM , then the waves reflected from it will be Doppler shifted: k' = k (1 - 2vM/c) (41)where c = w/k is the velocity of the waves. During a time Δt each wave is reflected c Δt / 2L times by the moving boundary, so the overall shift of the wave number will be: Δ k = - k vM Δt / L = - π ΔL/L2 (42)When the well expands from L1 to L2 the wave number will shift, by integration of equation (42), from π/L1 to π/L2. Therefore, the ground state wave function related to L1 will be shifted to the ground state function related to L2. See: Appendix: integrals 0∫ ∞xn exp(-a2x2) dx = Γ((n+1)/2) / 2an+1 Γ is the gamma function: Γ(n) = (n – 1)! <Γ(x + 1) = x Γ(x) Γ(x) Γ(x + 1/2) = <Γ(2 x) π1/2 / 22 x - 1 Γ(1/2) = π1/2 References On the net: December 2002, section on drag added July 2006, drag of a sphere added June 2008. By the author:
Links:
|