Selected final answers: Partial answer sheet, PDF page 181. This sheet does not cover every question.
Assume a hollow square electromagnetic waveguide with perfectly conducting walls, filled with vacuum or air. For side , the waveguide cutoff frequency is
A transverse electric mode can have one zero index, but not both. A transverse magnetic mode requires both indices positive. Thus the lowest cutoff belongs to the two degenerate transverse electric modes and , with
The corresponding cutoff vacuum wavelength is . Strictly, a travelling mode carrying power along the guide needs : at the cutoff itself, the axial wavenumber and group velocity are zero. The number above is therefore the lower frequency threshold, not a propagating mode with nonzero longitudinal energy flux exactly at that frequency. A dielectric filling would replace by the wave speed in that medium.
Treat the lake as a large thermal reservoir at , and take the copper's heat capacity constant over the cooling interval from . Although the actual cooling is irreversible, entropy is a state function, so its change can be calculated along a hypothetical reversible path:
The lake receives at constant temperature, giving
Therefore, neglecting other changes,
The positive total is the entropy generated by transferring heat across a finite temperature difference. The copper's entropy decreases, but the colder lake gains a greater amount.
Let be the length of the free straight segment from the current contact point to the mass, and let measure the contact point's rotation from first contact. For an inextensible taut rope, winding through angle uses length :
The point mass follows an involute of the pole's circular cross-section. Its instantaneous velocity is perpendicular to the free rope. To verify the geometry explicitly, use planar polar unit vectors at the contact point:
The rope's tension therefore does no work, and the speed stays at its initial value . The geometry gives . Consequently
The distance from the pole's axis is , so a point mass first reaches the pole when . The winding time around a circular pole is thus
This has units of time. The model assumes a light inextensible rope, no energy loss at contact and a point mass. The required tension diverges at the ideal final contact, so finite strength and size matter very close to impact.
The symmetry is electromagnetic reciprocity applied to mutual inductance. For thin closed circuits in vacuum, the magnetic vector potential generated by circuit B is
Stokes theorem converts the magnetic flux through circuit A into its magnetic vector potential circulation. Dividing by gives Neumann's mutual-inductance formula,
The dot product and separation are unchanged by interchanging A and B, so
Each has units . The winding orientations must be used consistently to determine the common sign. The same reciprocity holds in ordinary linear reciprocal magnetic media; it is not an unrestricted statement about biased nonreciprocal materials. Equivalently, a reversible magnetic energy function has equal mixed derivatives with respect to the two currents.
Assume all rings have the same circulation and orientation, and that there is no imposed background flow. Vortex lines describe the direction of vorticity, not the direction of velocity. For an incompressible fluid, the velocity induced by a thin vortex filament has the same Biot-Savart integral form as the magnetic field of a current filament:
A closely spaced coaxial array of vortex rings is therefore analogous to a solenoid. The velocity field passes along the common axis through the rings and returns outside a finite array. Its orientation reverses if every circulation is reversed. In a long array, far from its ends and vortex cores, the velocity field of a dense coaxial vortex-ring array is approximately
The exterior field is exactly zero only in the ideal continuous infinitely long limit, with the background velocity set to zero. In a finite array the streamlines spread out at one end, run back outside, and converge into the other. A large number of closely spaced rings alone does not guarantee that the total array length greatly exceeds ; if it does not, the end and return flows remain substantial. Different unspecified ring circulations would require superposing their actual signed contributions.
Figure 1.
Induced velocity through and around a finite array of co-oriented vortex rings
. Original Biot-Savart calculation for 81 equal rings with spacing 0.1 radii and a small numerical core smoothing. Streamlines lie in a plane containing the common axis; the grey markers sample intersections with the vortex rings, rather than fluid trajectories. Flow returns outside the finite array. Details.
Maxwell's displacement current is the electric-field contribution needed to make Ampère's law consistent with charge conservation. In macroscopic form,
The displacement current density has units . Since , it contains the vacuum field term and, in matter, the polarization-current contribution. It need not represent free charge moving through a material.
Taking the divergence gives . Gauss's law, , then yields the continuity equation . Without the displacement term, Ampère's original law would incorrectly require every current distribution to be divergence-free even while charge accumulates.
A charging parallel-plate capacitor makes this concrete. Ignoring fringing, , so the displacement current across its dielectric gap is , equal to the conduction electric current in its wire. The integral law
therefore gives the same answer whether the spanning surface cuts the wire or passes through the gap. Together with Faraday's law, this field term also permits self-consistent electromagnetic waves in vacuum, with speed .
For a cold, unmagnetized, collisionless plasma, take an equilibrium Electron density , Electrons of mass and charge , and approximately stationary ions. Use fields proportional to . The Electron displacement satisfies
Its electric polarization is , hence the cold collisionless plasma dielectric response
The negative susceptibility arises from the inertial Electron response, rather than a restoring force binding an Electron to a particular atom.
A transverse electromagnetic wave then obeys
For , the plasma supports propagation, with refractive index . The phase velocity is greater than , while the group velocity is ; their product is , and no superluminal signal is implied. Below the plasma frequency, the bulk wavenumber is imaginary and the field is evanescent. A sufficiently thick plasma largely reflects such radiation in the lossless model.
Including an Electron momentum-relaxation rate replaces the denominator by :
Its imaginary part describes absorption with the stated time convention. Finite Electron temperature also produces plasma screening of slowly varying electrostatic fields, with . That spatial, longitudinal screening problem is distinct from simply taking in the cold transverse-wave formula. Magnetic bias, ion motion and kinetic resonances require more general dielectric responses.
For weak fields in a linear isotropic material, magnetization is , with magnetic susceptibility . Paramagnetism has , while diamagnetism has . In either case, . These weak responses generally disappear when the applied field is removed, unlike hysteretic ferromagnetism.
A paramagnet can contain permanent magnetic moments from unpaired Electrons or atomic angular momentum. Alignment lowers the energy , while thermal disorder opposes it. For independent classical moments of magnitude , their weak-field average projection is , so
This is Curie's law; it assumes weak magnetization and no important interactions between moments. With quantum angular momentum , replace by . Strong-field magnetization saturates instead of growing linearly without limit. Conduction Electrons also exhibit Pauli paramagnetism: only occupations near the Fermi surface readily redistribute between spin states, giving a nearly temperature-independent susceptibility when .
Diamagnetism is an induced orbital response opposing the applied magnetic field. For bound Electrons, a weak field changes orbital motion and induces an opposing magnetic moment. For an isotropic collection of atoms, the elementary orbital estimate is
where the sum runs over the Electrons of one atom. It is usually small and weakly temperature dependent. Delocalized Electrons also have orbital Landau diamagnetism; for an ideal three-dimensional free Electron gas with a quadratic isotropic dispersion and electron spin with , . Both spin and orbital contributions must be included to determine the total response. Diamagnetism is not a consequence of permanent dipoles choosing to align against the field; that orientation would increase their energy.
The Poynting vector is the electromagnetic energy flux density,
In vacuum it equals . For a linear nondispersive medium with time-independent constitutive coefficients, Maxwell's equations and the identity give Poynting's theorem,
The source term is work done by the fields on free charges; for an Ohmic conductor it is the positive local Joule heating rate . Integrating over a volume states that loss of stored field energy plus outward energy flow equals energy transferred to matter. In dispersive media, such as a plasma, material motion contributes to stored energy and the simple nondispersive expression cannot be used indiscriminately.
For a vacuum plane wave, and , and the propagation direction form a right-handed triad. With real peak field , the time-averaged wave intensity is
More generally, complex peak amplitudes give . Along a transmission line, the cross-section integral of the longitudinal Poynting vector is the transmitted electrical power . Energy flows through the electromagnetic fields in the surrounding dielectric, even though the conductors guide the currents.
Distributed circuit model. A transmission line must be treated as a distributed system when its propagation time is appreciable compared with the signal timescale. Divide a uniform two-conductor line into elements of length . Each has series resistance and inductance , and shunt electrical conductance and capacitance . The primes denote quantities per unit length, with units , , and , respectively.
Kirchhoff's voltage and current laws give the Telegrapher's equations:
For a lossless line, , differentiation gives and the same wave equation for . The propagation speed and characteristic impedance are
A general lossless signal has and . Thus characteristic impedance is the voltage/current ratio of one travelling wave, not simply the resistance of the wires. It can be real even when the line dissipates no energy, because energy is carried away by propagation.
For sinusoidal quantities with convention , the general lossy-line formulas are
Choose the square-root branches for forward attenuation and positive forward power. The real part of is attenuation per unit length; its imaginary part is the phase wavenumber. For small losses, attenuation is approximately , using the lossless real in this estimate.
Fields and geometry. Two separated conductors support a transverse electromagnetic mode, with electric and magnetic fields transverse to the propagation axis. For a homogeneous dielectric and negligible conductor penetration, both examples below have and hence . The distributed circuit model follows from the field energy and charges in each slice.
For a coaxial cable with inner radius and inner surface of the outer conductor at , Gauss's law and Ampère's law give
The electric field is radial and the magnetic field circles the axis. Integration of voltage and field energy gives
The equal and opposite return current on the outer conductor confines the ideal external field to the annulus.
For wide parallel plates of width separated by , neglect fringing. The electric field crosses the gap, , while the magnetic field lies in the gap parallel to the plates and perpendicular to propagation, . Therefore
For a forward wave, points along the line. Integrating this Poynting vector over either cross-section gives instantaneous power . For complex peak amplitudes, average power is , or for a single forward wave on a lossless line. The conductors guide an energy flux predominantly in the dielectric between them.
Figure 1.
Electric and magnetic fields in coaxial and parallel-plate transmission lines
. Original transverse cross-sections for a forward TEM wave. Blue arrows show the electric field and orange arrows the magnetic field strength. Dots denote current out of the page and crosses current into it. Parallel-plate edge fields and coaxial conductor penetration are neglected. Details.
Terminations and reflection. At a load , require . If and are the incident and reflected voltage amplitudes there, the total current is . Thus
A matched load, , absorbs the incident power with no reflection. An open circuit has , a current node and a doubled voltage at the end. A short circuit has , a voltage node and a doubled current there. A purely reactive load on a lossless line has and accepts no average power. For real , the reflected power fraction is and the accepted fraction is .
Incident and reflected waves produce a standing-wave pattern with voltage standing wave ratio . The input impedance of a lossless length is
This gives for a short and for an open. A quarter-wave section gives , enabling impedance transformation. After a source launches a pulse, a distant termination affects it only after the reflected pulse returns; the return time is . Real conductor resistance, dielectric loss and departures from a TEM field profile limit the ideal model at high frequency.
An Euler strut is a slender elastic column under axial compression. Let its bending stiffness about the weaker section axis be and its length be . With small lateral deflection and a dead compressive load , the quadratic total potential energy is
The first term is bending energy; the negative term is work released by the load as bending shortens the end-to-end distance. Stationarity gives . For pinned ends, at , and requires . The first instability is Euler's critical load,
Below this load, the straight state has positive quadratic energy against these perturbations. At the threshold, the lowest bending mode has zero quadratic restoring energy; above it the straight state is unstable. For a rectangular section, about the appropriate bending axis, so the smaller second moment of area controls the first buckling direction.
Endpoint constraints matter. Write : for pinned-pinned, for fixed-free and for fixed-fixed conditions. This elastic instability differs from compressive material failure. Imperfections cause bending before the ideal threshold, and a short column may yield before Euler buckling occurs.
A gravitational assist changes a spacecraft's velocity relative to the Sun by scattering it in the gravitational field of a moving planet. In the planet's rest frame, ignore other forces during the encounter and take the spacecraft mass negligible. The incoming and outgoing asymptotic velocities have equal magnitude,
because the encounter is conservative, but their directions differ along a hyperbolic orbit. If the planet has approximately constant heliocentric velocity , the heliocentric spacecraft velocities are . Its specific energy change is
Thus a flyby can increase or decrease heliocentric energy depending on its geometry. The planet exchanges momentum and orbital energy with the spacecraft; its speed change is extremely small because its mass is much greater. There is no energy creation and no increase of asymptotic speed in the planet's own frame.
For a point-mass planet, conservation of energy and angular momentum gives the hyperbolic eccentricity and deflection angle
where is periapsis distance and the impact parameter. The velocity-vector change has magnitude , so the maximum energy-change magnitude for this fixed deflection is . A small periapsis increases the turning angle, but collision, atmosphere and tidal forces limit how close a real spacecraft may pass. Multiple encounters allow large changes without expending the propellant that equivalent direct velocity changes would require.
Take axes fixed in a rigid body, chosen along its principal axes, with moments of inertia . The angular momentum components are . Newton's rotational equation is , while the rotating-frame transport rule gives
Thus Euler's equations are
The nonlinear cross terms arise because the principal axes rotate. Even with no external torque, the body-frame components of angular momentum can vary while the angular momentum vector stays fixed in inertial space.
With zero torque, the equations conserve
Uniform rotation about any principal axis is a torque-free solution. If , linearization about the intermediate-axis solution gives
The positive coefficient gives exponential growth: intermediate-axis rotation is unstable, whereas rotation about the largest or smallest moment is linearly stable. For a symmetric body, Euler's equations also describe the precession of the transverse angular velocity around its symmetry axis. Torques from gravity or an applied couple enter through the right-hand sides rather than being absorbed into the inertial cross terms.
Archimedes' principle states that the resultant hydrostatic pressure force on a body equals the weight of the fluid displaced, directed opposite to gravity. Let be the displaced fluid region, with outward normal from the body into the fluid. Since the fluid pressure force is inward,
Hydrostatic balance is , with downward gravitational acceleration . Therefore
for uniform fluid density and gravity. Uniform atmospheric pressure gives no net force on a closed surface. For a floating body, the same displaced-volume construction includes the free surface; any significant air buoyancy can be included as a second displaced fluid.
A freely floating body of mass satisfies , so its displaced volume is . A homogeneous body of density floats with submerged fraction ; if , hydrostatic buoyancy alone cannot support it. The buoyant force acts through the centre of buoyancy, the centroid of the displaced volume for constant density, while weight acts through the centre of mass. Their relative positions determine rotational equilibrium and stability. For a fully submerged body, a centre of mass below the centre of buoyancy gives a restoring torque; a floating body's changing displaced shape additionally requires the metacentric height criterion. Surface tension, externally driven flow and contact forces can alter the simple force balance.
Following the moving fluid. The Eulerian velocity field assigns a velocity to each position. A fluid particle satisfies , so the chain rule gives the material derivative, also called the convective derivative,
Even in steady flow a particle can accelerate because it moves between regions with different velocity. For example, flow around a curved obstacle changes direction despite . Mass conservation is
For a constant-density incompressible fluid, .
Euler equations for an inviscid fluid and Bernoulli's equation. Neglect viscosity and take the body force per unit mass to be . Newton's second law for a fluid element gives
For steady flow, dot with displacement along a streamline. The tangential component of is the derivative of along that line. For constant density, integration gives
In uniform gravity, . Thus pressure, kinetic energy per unit mass and gravitational potential energy per unit mass exchange along a streamline. For a barotropic compressible fluid, replace by . Different streamlines need not have the same Bernoulli constant.
Vorticity and potential flow. Vorticity is . A small fluid element rotates locally with angular velocity after its deformational part is separated. The vector identity
gives for steady constant-density inviscid flow. Thus an irrotational region has one common Bernoulli constant on each connected region. In a simply connected irrotational domain, introduce a velocity potential . Incompressibility then gives Laplace's equation, . For unsteady irrotational flow,
An impermeable stationary wall fixes , but ideal potential flow does not enforce no slip. In a multiply connected domain, irrotational flow may still have nonzero circulation around an excluded vortex core, so a globally single-valued potential requires care.
For example, uniform two-dimensional flow past a circular cylinder has outside . Its normal velocity vanishes at the cylinder, and its surface tangential velocity is . Bernoulli's equation predicts a symmetric pressure distribution and zero drag. This is D'Alembert's paradox: the assumptions omit the viscous boundary layers, separation and wake responsible for real drag.
For an incompressible fluid of constant kinematic viscosity , the curl of the motion equation gives
The first term stretches and tilts vortex lines; the second diffuses vorticity. In two-dimensional incompressible inviscid flow the stretching term is absent and vorticity is carried with the fluid. In a three-dimensional inviscid flow it can be amplified by stretching. Density gradients not parallel to pressure gradients provide an additional source in a general compressible fluid.
Viscous flow. A Newtonian fluid's shear stress is proportional to rate of strain. For constant viscosity and incompressibility the Navier-Stokes equation is
with dynamic viscosity and kinematic viscosity . At an ordinary solid wall, the no-slip boundary condition sets fluid velocity equal to the wall velocity. Viscosity transfers momentum between neighbouring layers and converts mechanical energy into heat; unlike the pressure force, its shear contribution depends on velocity gradients.
For plane Couette flow between plates , with the upper plate moving at , the steady solution is , with shear stress . For pressure-driven flow through a pipe of radius , a constant pressure gradient gives the parabolic Poiseuille profile
The Reynolds number compares convective inertia with viscous force. Small Reynolds number often permits the Stokes-flow approximation; large Reynolds number can produce thin boundary layers and, under suitable conditions, turbulence. It does not justify neglecting the wall conditions. Bernoulli's inviscid constant is generally not conserved along a viscous streamline, and pipe flow instead shows a pressure loss corresponding to dissipation.
A Carnot cycle consists of four reversible stages: isothermal expansion at the hot-reservoir temperature , adiabatic expansion to , isothermal compression at , and adiabatic compression back to the starting state. The adiabatic stages exchange no heat. For one mole of an ideal gas, the two isothermal stages have heats and in magnitude. The reversible adiabatic relations ensure .
More generally, the total entropy change over any reversible cycle is zero. The two heat exchanges therefore satisfy
With work output , its thermal efficiency is
The temperatures must be absolute. The second law makes this the maximum efficiency for an engine exchanging heat only with those two reservoirs; irreversible entropy production lowers the efficiency.
Label the initial state 1, the end of compression 2, the end of constant-pressure heating 3, and the end of adiabatic expansion 4. The volume sequence is , with . The Diesel cycle has a compression adiabat , a horizontal constant-pressure heat-input segment , an expansion adiabat , and a vertical constant-volume cooling segment .
Figure 1. . Original dimensionless example with V2/V1=0.2, V3/V1=0.4 and gamma=1.4. The curved segments are reversible adiabats, the upper horizontal segment adds heat at fixed pressure and the right vertical segment rejects heat at fixed volume. The enclosed area is the positive net work output. Details.
For an ideal gas, a reversible adiabat obeys , so it falls more steeply than an isotherm. The expansion adiabat lies above the compression adiabat at equal volume between and , giving a clockwise cycle with positive work output.
For one mole, , and constant heat capacities satisfy . The specified diatomic-gas model has
A reversible adiabatic process satisfies . Constant-pressure expansion satisfies . Applying these in order gives
Every volume ratio is dimensionless, so the expressions have units of temperature. For the intended volume ordering, the minimum temperature is and the maximum is .
Use heat into the gas as positive. The two adiabatic stages have . During constant-pressure heating and constant-volume cooling,
The first law over a complete cycle gives , since the internal energy returns to its initial value. Therefore
Substitute the temperatures:
Hence
Equivalently, define compression ratio and cutoff ratio :
These formulas use the ideal constant-heat-capacity model and reversible adiabatic strokes. The heat-addition and heat-rejection strokes occur over changing gas temperatures, unlike the isothermal heat exchanges of a Carnot engine.
For , Taylor expansion gives . The vanishing-cutoff limit of the Diesel cycle is
The maximum temperature tends to , while the minimum stays at . A Carnot engine between those limiting extrema has
so the limiting efficiencies coincide. At exactly , however, both heat input and net work vanish, and the cycle retraces its compression adiabat. The efficiency there is an undefined ; the boxed value is the continuous limit for neighbouring nonzero cycles, not finite work at zero heat input.
For a finite cutoff ratio , the Carnot comparison uses , giving . The Diesel cycle has smaller efficiency. For example,
which makes its rejected/input heat ratio larger than the Carnot ratio.
An equation of state relates equilibrium state variables such as pressure, volume, temperature and material amount. At fixed composition, it reduces the number of independent thermodynamic variables. It does not by itself specify all caloric properties: information such as heat capacity or internal energy is also needed.
At a liquid-gas first-order transition, an equilibrium isotherm can contain a constant-pressure interval in which changing total volume changes the phase fractions. A metastable single-phase extrapolation can instead show a region with ; that region has negative isothermal compressibility and is mechanically unstable. Thermodynamic phase coexistence requires equal temperature, pressure and chemical potential, and a consistent free-energy description replaces the unstable region by coexistence. For an appropriate analytic fluid equation of state, a liquid-gas critical point occurs where .
These tests diagnose the phase behaviour represented by that equation of state. They do not guarantee detection of every possible transition involving an additional order parameter which the mechanical equation omits.
For the given one-mole model,
At fixed positive temperature its isotherm is an ideal-gas hyperbola shifted right by : pressure diverges as and tends to zero as . Higher temperatures give higher pressure at any fixed allowed volume.
Figure 1. . Original dimensionless isotherms with a common vertical asymptote at V=b. They are monotonically decreasing and have no coexistence plateau or critical inflection. Details.
For positive , the volume is the effective volume available to molecular motion. Thus describes excluded-volume repulsion or an effective incompressible molecular core. It is a molar volume in this one-mole equation, not necessarily the literal sum of molecular geometric volume. Setting recovers the ideal gas law. No attractive-energy term is included.
Start with the fundamental relation at fixed material amount, , and expand entropy in independent variables :
The last equality uses the supplied Maxwell relation. Consequently
For the excluded-volume ideal gas, , so the bracket vanishes:
Once depends only on temperature, cannot depend on volume either. Its stipulated temperature independence therefore makes it a constant, giving with arbitrary energy reference . Excluded volume changes the pressure without adding the volume-dependent attractive internal energy that a more general interacting-gas model might have.
Use the enthalpy . The equation of state gives , so
At fixed pressure, is constant. Therefore
Since and are constants, this ratio is independent of both temperature and pressure. Although the gas is not mechanically ideal when , its constant-pressure/constant-volume heat-capacity difference remains in this model.
For a reversible adiabatic change, , equivalently . Substitute :
Integration, with fixed positive reference values , gives
Since and ,
The relation applies to the quasistatic reversible adiabat. Adiabatic heat isolation alone does not enforce it for an irreversible process which generates entropy. The shifted volume , rather than itself, is the appropriate coordinate in this model.
At every , ,
There is no constant-pressure coexistence interval, unstable loop or finite-temperature critical inflection. The isothermal compressibility is positive:
Thus this equation of state predicts no liquid-gas phase transition in its physical domain. Its divergence at the excluded-volume boundary is a repulsive constraint, not a liquid-gas coexistence mechanism. This conclusion concerns the stated model and cannot rule out an ordering transition requiring additional physics or an additional thermodynamic variable. Extrapolation to absolute zero is also outside the classical gas model's reliable regime.
For a classical dilute ideal gas in thermal equilibrium, the probability of a single-particle velocity is proportional to the Boltzmann factor for kinetic energy, . Velocities with speed between and fill a spherical shell of velocity-space volume . The normalized Maxwell speed distribution is therefore
The factor counts states, rather than representing an extra Boltzmann energy factor.
Set and use . Normalization cancels in ratios, so the given recurrence gives
Consequently
The root mean square speed is larger than the mean speed, consistent with the positive variance of the speed distribution.
Choose a flat wall and let measure the angle from the direction towards it. Particles of speed in solid angle occupy fraction of an isotropic gas. In time , those within normal distance reach the wall. Their number per area per time is therefore .
Only the incoming hemisphere contributes. Integrating gives the one-way collision rate per unit area,
Its units are . This is not : selecting the incoming half of the particles must also account for their mean normal velocity. For effusion through a thin hole whose diameter is much less than the mean free path, the same one-way rate applies because collisions within the opening can be neglected.
Wall arrivals are not an unbiased sample of the gas: faster atoms arrive more frequently. The normalized speed distribution among incident atoms is
The angular factor in the arrival rate cancels when calculating the mean kinetic energy. Thus
Using gives
By contrast, an atom chosen randomly from the bulk has mean kinetic energy . The arrival weighting raises the normal-direction contribution to while each tangential direction still contributes . This is also the mean energy carried by atoms effusing through a thin aperture from an equilibrated reservoir.
Denote the heated and cold vessels by subscripts . For argon, the opposing effusive particle rates are equal in the steady state:
Because and ,
This is thermal transpiration in the free-molecular limit. It fixes a ratio; the initial conditions determine the absolute pressures.
For the stated two identical vessels, take their fixed volume to be and the combined system sealed. Initially it contains atoms. Hence
Let . Combining this inventory with gives the closed-vessel thermal transpiration result
Heating makes the hot side less dense but at higher pressure. Both expressions reduce to when . For a small imposed temperature rise , the pressure difference starts as .
If instead the cold vessel is interpreted as an externally maintained pressure reservoir with , the answer becomes . That alternative requires the fixed-pressure reservoir assumption; it does not conserve the initial particle inventory of two sealed equal-volume vessels. Merely describing the vessels as large does not keep their pressures fixed through the redistribution. The pressures in the boxed answer assume the natural closed-vessel interpretation.
Let L be the vessel initially containing argon and R the vessel initially containing helium. Both have volume and fixed temperature , and the aperture has area . The incoming and outgoing rates of each species must be considered separately. For species ,
Thus its density difference relaxes exponentially with time constant
The two-vessel effusive mixing solution for the four densities is
Using atomic masses of approximately and atomic mass units gives . Helium therefore transfers faster. Each species eventually reaches in each vessel, but its transient is different. In particular, the total density temporarily rises in the initially argon vessel and falls in the initially helium vessel, before both return to .
Figure 1.
Species densities during effusive mixing of initially separated argon and helium
. Original solutions of the equal-volume, equal-temperature rate equations, with the helium rate faster by sqrt(10). The dashed total densities show the temporary pressure imbalance even though both temperatures are equal. Each species ultimately has density n0/2 in each vessel. Details.
For unequal maintained temperatures , replace the rate equation by
With each species' total density sum still , its final densities and relaxation time are
Every species has the same final density ratio and partial-pressure ratio . Thus there is no permanent composition separation in this ideal steady state, although the rates remain mass dependent. A maintained hot side has lower density and higher pressure. Equal particle rates do not imply thermal equilibrium: atoms arriving from the hotter side carry more energy, so external heating and cooling are required to hold different temperatures. The model assumes thin-hole effusion, local equilibrium within each vessel and negligible aperture volume.
Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path small compared with the gradient length, and relaxation times short compared with the observation time.
For a dilute hard-sphere gas,
where is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.
A simple kinetic estimate makes the common structure explicit. For a quantity carried per particle, a molecule moving at angle to brings approximately . Expanding to first order and averaging an isotropic velocity distribution gives
since and . This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.
Momentum transport and viscosity. For a shear flow with mean velocity , the transported directed momentum per molecule is . The net flux of momentum in the direction is
Here is dynamic viscosity in . The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is , with units . Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.
Heat transport. Let be the thermal heat capacity per molecule, in . Substituting gives Fourier's law,
Thermal conductivity has units . For a monatomic gas, translational equipartition gives . Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.
Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,
The diffusivity has units ; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.
Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,
Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as ; at fixed pressure, and
Kinematic viscosity and thermal diffusivity also grow roughly as at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change . Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.
Heat diffusion equation. For a stationary gas with constant volumetric heat capacity , constant thermal conductivity and negligible compression work, the energy balance in any fixed volume is
Here is a possible heat generation rate per volume. The divergence theorem and Fourier's law yield
Thermal diffusivity has units and is of order . In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.
Representative solutions. With no source, a steady one-dimensional slab between temperatures at has and uniform heat flux . Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.
For an infinite medium and an initial sinusoidal temperature perturbation,
Short spatial wavelengths decay faster, and the characteristic smoothing length after time is . A finite slab with both ends held at has a sine expansion with factors ; its slowest relaxation time is .
An instantaneous heat input per unit transverse area at in an otherwise uniform infinite medium gives the Gaussian solution
Its integral conserves the deposited heat. The peak falls as while its width grows as ; diffusion redistributes energy instead of transporting a sharp front at a single velocity.
Finally, a semi-infinite body initially at and held at surface temperature from has
It has and approaches far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

Articles by others on the same topic (0)

There are currently no matching articles.