IPhO 2025, theory — Задача 1. Champagne!

Автор: Olympiads XYZ · транскрипция на официалните материали

Проверена срещу оригинала на 13.9.2026 от същия модел, който я е транскрибирал (без независима проверка)

Съдържание

УсловиеРешение

Theory — International Physics Olympiad FRANCE 2025 — Q3 — Champagne! (10 points) — English (Official) · 17 юли 2025 г. · 10 т.

Внимание: Бележка към темата

Фрагмент от решение; първата страница на прозореца (стр. 6) е покрита от предишния прозорец.

Условие

Warning: Excessive alcohol consumption is harmful to health and drinking alcohol below legal age is prohibited.

Champagne is a French sparkling wine. Fermentation of sugars produces carbon dioxide (CO2_2) in the bottle. The molar concentration of CO2_2 in the liquid phase cc_\ell and the partial pressure PCO2P_{\mathrm{CO_2}} in the gas phase are related by c=kHPCO2c_\ell = k_{\mathrm{H}} P_{\mathrm{CO_2}}, known as Henry's law and where kHk_{\mathrm{H}} is called Henry's constant.

Data

  • Surface tension of champagne σ=47×103 Jm2\sigma = 47 \times 10^{-3}\ \mathrm{J \cdot m^{-2}}
  • Density of the liquid ρ=1.0×103 kgm3\rho_\ell = 1.0 \times 10^{3}\ \mathrm{kg \cdot m^{-3}}
  • Henry's constant at T0=20CT_0 = 20\,{}^\circ\mathrm{C}, kH(20C)=3.3×104 molm3Pa1k_{\mathrm{H}}(20\,{}^\circ\mathrm{C}) = 3.3 \times 10^{-4}\ \mathrm{mol \cdot m^{-3} \cdot Pa^{-1}}
  • Henry's constant at T0=6CT_0 = 6\,{}^\circ\mathrm{C}, kH(6C)=5.4×104 molm3Pa1k_{\mathrm{H}}(6\,{}^\circ\mathrm{C}) = 5.4 \times 10^{-4}\ \mathrm{mol \cdot m^{-3} \cdot Pa^{-1}}
  • Atmospheric pressure P0=1 bar=1.0×105 PaP_0 = 1\ \mathrm{bar} = 1.0 \times 10^{5}\ \mathrm{Pa}
  • Gases are ideal with an adiabatic coefficient γ=1.3\gamma = 1.3

Part A. Nucleation, growth and rise of bubbles

Immediately after opening a bottle of champagne at temperature T0=20CT_0 = 20\,{}^\circ\mathrm{C}, we fill a glass. The pressure in the liquid is P0P_0 and its temperature stays constant at T0T_0. The concentration cc_\ell of dissolved CO2_2 exceeds the equilibrium concentration and we study the nucleation of a CO2_2 bubble. We note aa its radius and PbP_\mathrm{b} its inner pressure.

In the liquid, the concentration of dissolved CO2_2 depends on the distance to the bubble. At long distance we recover the value cc_\ell and we note cbc_\mathrm{b} the concentration close to the bubble surface. According to Henry's law, cb=kHPbc_\mathrm{b} = k_{\mathrm{H}} P_\mathrm{b}. We furthermore assume in all the problem that bubbles contain only CO2_2.

Since ccbc_\ell \neq c_\mathrm{b}, CO2_2 molecules diffuse from areas of high to low concentration. We assume also that any molecule from the liquid phase reaching the bubble surface is transferred to the vapour.

In practice, bubbles mainly grow from pre-existing gas cavities. Consider then a bubble with initial radius a040 μma_0 \approx 40\ \mathrm{\mu m}. The number of moles of CO2_2 transferred at the bubble's surface per unit area and time is noted jj. Two models are possible for jj.

  • model (1) j=Da(ccb)j = \dfrac{D}{a}(c_\ell - c_\mathrm{b}) where DD is the diffusion coefficient of CO2_2 in the liquid.
  • model (2) j=K(ccb)j = K(c_\ell - c_\mathrm{b}) where KK is a constant here.

Experimentally, the bubble radius a(t)a(t) is found to depend on time as shown in Fig. 2. Here c4c0c_\ell \approx 4c_0, and since bubbles are large enough to be visible, the excess pressure due to surface tension can be neglected and PbP0P_\mathrm{b} \approx P_0.

Eventually bubbles detach from the bottom of the glass and continue to grow while rising. Fig. 3. shows a train of bubbles. The bubbles of the train have the same initial radius and are emitted at a constant frequency fb=20 Hzf_\mathrm{b} = 20\ \mathrm{Hz}.

For the range of velocities studied here, the drag force FF on a bubble of radius aa moving at velocity vv in a liquid of dynamic viscosity η\eta is given by Stokes' law F=6πηavF = 6\pi\eta a v. Measurements show that at any moment in time, the bubble can be assumed to be travelling at its terminal velocity.

The quasi-stationary growth of bubbles with rate qa=dadtq_a = \frac{\mathrm{d}a}{\mathrm{d}t} still applies during bubble rise.

There are NbN_\mathrm{b} nucleation sites of bubbles. Assume that the bubbles are nucleated at a constant frequency fbf_\mathrm{b} at the bottom of a glass of champagne (height HH_\ell for a volume VV_\ell), with a0a_0 still negligible. Neglect diffusion of CO2_2 at the free surface.

Part B. Acoustic emission of a bursting bubble

Small bubbles are nearly spherical as they reach the free surface. Once the liquid film separating the bubble from the air thins out sufficiently, a circular hole of radius rr forms in the film and, driven by surface tension, opens very quickly (Fig. 4. left). The hole opens at constant speed vfv_\mathrm{f} (Fig. 4. right). The film outside the rim remains still, with constant thickness hh.

Due to dissipative processes, only half of the difference of the surface energy between tt and t+dtt + \mathrm{d}t of the rim and the accumulated liquid is transformed into kinetic energy. We further assume that the variation of the surface of the rim is negligible compared to that of the film.

When the film bursts, it releases internal pressure and emits a sound. We model this acoustic emission by a Helmholtz resonator: a cavity open to the atmosphere at P0P_0 through a bottleneck aperture of area SS (Fig. 5. left). In the neck, a mass mpm_\mathrm{p} makes small amplitude position oscillations due to the pressure forces it experiences as the gas in the cavity expands or compresses adiabatically. The gravity force on mpm_\mathrm{p} is negligible compared to pressure forces. Let V0V_0 be the volume of gas under the mass mpm_\mathrm{p} for P=P0P = P_0 as z=0z = 0.

The Helmholtz model may be used for a bubble of radius aa. V0V_0 is the volume of the closed bubble. From litterature, the mass of the equivalent of the piston is mp=8ρgr3/3m_p = 8\rho_\mathrm{g} r^3/3 where rr is the radius of the circular aperture and ρg=1.8 kgm3\rho_\mathrm{g} = 1.8\ \mathrm{kg \cdot m^{-3}} is the density of the gas (Fig. 5. right). During the bursting process, rr goes from 0 to rcr_\mathrm{c}, given by rc=23a2ρg0σr_\mathrm{c} = \dfrac{2}{\sqrt{3}} a^2 \sqrt{\dfrac{\rho_\ell g_0}{\sigma}}. At the same time, the frequency of emitted sound increases until a maximum value of 40 kHz40\ \mathrm{kHz} and the bursting time is tb=3×102 mst_b = 3 \times 10^{-2}\ \mathrm{ms}.

Part C. Popping champagne

In a bottle, the total quantity of CO2_2 is nT=0.2 moln_\mathrm{T} = 0.2\ \mathrm{mol}, either dissolved in the volume VL=750 mLV_\mathrm{L} = 750\ \mathrm{mL} of liquid champagne, or as a gas in the volume VG=25 mLV_\mathrm{G} = 25\ \mathrm{mL} under the cork (Fig. 6. left). VGV_\mathrm{G} contains only CO2_2. The equilibrium between both CO2_2 phases follows Henry's Law. We suppose that the fast gaseous CO2_2 expansion when the bottle is opened, is adiabatic and reversible. Ambient temperature T0T_0 and pressure P0=1 barP_0 = 1\ \mathrm{bar} are constant.

Another step of champagne production (not described here) leads to the following values of PiP_\mathrm{i} that we will use for the next questions: Pi=4.69 barP_\mathrm{i} = 4.69\ \mathrm{bar} at T0=6CT_0 = 6\,{}^\circ\mathrm{C} and Pi=7.45 barP_\mathrm{i} = 7.45\ \mathrm{bar} at T0=20CT_0 = 20\,{}^\circ\mathrm{C}.

During bottle opening, two different phenomena can be observed, depending on T0T_0 (Fig. 6. right).

  • either a blue fog appears, due to the formation of solid CO2_2 crystals (but water condensation is inhibited);
  • or a grey-white fog appears, due to water vapor condensation in the air surrounding the bottleneck. In this latter case, there is no formation of CO2_2 solid crystals.

The saturated vapor pressure PsatCO2P_\mathrm{sat}^{\mathrm{CO_2}} for the CO2_2 solid/gas transition follows : log10(PsatCO2P0)=ABT+C\log_{10}\left(\dfrac{P_\mathrm{sat}^{\mathrm{CO_2}}}{P_0}\right) = A - \dfrac{B}{T + C} with TT in K, A=6.81A = 6.81, B=1.30×103 KB = 1.30 \times 10^{3}\ \mathrm{K} and C=3.49 KC = -3.49\ \mathrm{K}.

During bottle opening, the cork stopper pops out. We now determine the maximum height HcH_\mathrm{c} it reaches. Assume that the friction force FF due to the bottleneck on the cork stopper is F=αAF = \alpha A where AA is the area of contact and α\alpha is a constant to determine. Initially, the pressure force slightly overcomes the friction force. The cork's mass is m=10 gm = 10\ \mathrm{g}, its diameter d=1.8 cmd = 1.8\ \mathrm{cm} and the length of the cylindrical part initially stuck in the bottleneck is 0=2.5 cm\ell_0 = 2.5\ \mathrm{cm}. Once the cork has left the bottleneck, you can neglect the net pressure force.

Photo of a champagne flute filled with sparkling wine with rising bubbles.
Fig. 1. A glass filled with champagne.
Graph of bubble radius a in micrometres from about 40 to 270 versus time t in seconds from 0 to about 0.9, shown as data points with error bars growing roughly linearly.
Fig. 2. Time evolution of CO2_2 bubble radius in a glass of champagne (adapted from [1]).
Horizontal photo of a train of progressively larger bubbles; below the photo an axis z points right with arrows for velocity u_z, gravity g_0 pointing left, and a scale bar of 1 mm.
Fig. 3. A train of bubbles. The photo is rotated horizontally for the page layout (adapted from [1]).
Left column: four sketches (α to δ) of a bubble at the liquid surface with a film hole of radius r(t) opening, film thickness h, rim moving at speed v_f. Right column: top view of the pierced film showing growing hole radius r(t) and rim annulus of width v_f dt; below, a cross-section with rim radius R_m, film thickness h, hole radius r(t) and velocity v_f.
Fig. 4. (Left) (α) Bubble at the surface: (1) liquid, (2) air at pressure P0P_0 and (3), CO2_2 at pressure PbP_\mathrm{b}, (β) and (γ) retraction of the liquid film, where the rim is in dark blue, (δ) bubble collapse. (Right) Retraction of the liquid film at time tt. Top: sketch of the pierced film seen from above. Bottom: cross-section of the rim and the retracting film. During dt\mathrm{d}t the rim accumulates nearby liquid (dotted).
Left: a cavity with pressure P(t) under a neck where a piston of mass m_p of area S is displaced by z above equilibrium position z=0, open to pressure P_0; the z axis points up. Right: a spherical bubble of radius a below a free surface with a circular hole of diameter 2r.
Fig. 5. (Left) a Helmholtz resonator. (Right) a bubble as an oscillator.
Left: schematic of a bottle neck with cork of diameter d, embedded length ℓ_0, below it a gas region labelled P_i, V_G and liquid labelled c_ℓ = k_H P_i, V_L, with pressure P_0 above. Right: two photos of a champagne bottle being opened with a fog plume at the neck.
Fig. 6. Left: traditional bottleneck: (1) surrounding air, (2) cork stopper, (3) headspace, (4) liquid champagne. Right: Two phenomena observed while opening the bottle at two different temperatures (adapted from [2]).

A.1 Express the pressure PbP_\mathrm{b} in terms of P0P_0, aa and σ\sigma. [0,2 т.]

A.2 Express the critical radius aca_\mathrm{c} above which a bubble is expected to grow in terms of P0,σ,cP_0, \sigma, c_\ell and c0c_0 where c0=kHP0c_0 = k_\mathrm{H} P_0. Calculate numerically aca_\mathrm{c} for c=4c0c_\ell = 4c_0. [0,5 т.]

A.3 Express the number of CO2_2 moles in the bubble ncn_\mathrm{c} in terms of aa, P0P_0, T0T_0 and ideal gas constant RR. Find a(t)a(t) for both models. Indicate which model explains the experimental results in Fig. 2. Depending on your answer, calculate numerically KK or DD. [1,2 т.]

A.4 Give the expression of the main forces exerted on a vertically rising bubble. Obtain the expression of v(a)v(a). Give a numerical estimate of η\eta using ρ\rho_\ell, g0g_0 and quantities measured on Fig. 3. [0,8 т.]

A.5 Express the radius aHa_{H_\ell} of a bubble reaching the free surface in terms of height travelled HH_\ell, growth rate qa=dadtq_a = \frac{\mathrm{d}a}{\mathrm{d}t}, and any constants you may need. Assume aHa0a_{H_\ell} \gg a_0 and qaq_a constant, and give the numerical value of aHa_{H_\ell} with H=10 cmH_\ell = 10\ \mathrm{cm} and qaq_a corresponding to Fig. 2. [0,5 т.]

A.6 Write the differential equation for c(t)c_\ell(t). Obtain from this equation the characteristic time τ\tau for the decay of the concentration of dissolved CO2_2 in the liquid. [1,1 т.]

B.1 Express vfv_\mathrm{f} in terms of ρ,σ\rho_\ell, \sigma and hh. [1,1 т.]

B.2 Express the frequency of oscillation f0f_0 of mpm_\mathrm{p}. Hint: for ϵ1\epsilon \ll 1, (1+ϵ)α1+αϵ(1+\epsilon)^{\alpha} \approx 1 + \alpha\epsilon. [1,1 т.]

B.3 Find the radius aa and the thickness hh of the champagne film separating the bubble from the atmosphere. [1,1 т.]

C.1 Give the numerical value of the pressure PiP_\mathrm{i} of gaseous CO2_2 in the bottle for T0=6CT_0 = 6\,{}^\circ\mathrm{C} and T0=20CT_0 = 20\,{}^\circ\mathrm{C}. [0,4 т.]

C.2 Give the numerical value TfT_\mathrm{f} of the CO2_2 gas at the end of the expansion, after opening a bottle, if T0=6CT_0 = 6\,{}^\circ\mathrm{C} and if T0=20CT_0 = 20\,{}^\circ\mathrm{C}, if no phase transition occured. Choose which statements are true (several statements possible):

  1. At T0=6CT_0 = 6\,{}^\circ\mathrm{C} a grey-white fog appears while opening the bottle.
  2. At T0=6CT_0 = 6\,{}^\circ\mathrm{C} a blue fog appears while opening the bottle.
  3. At T0=20CT_0 = 20\,{}^\circ\mathrm{C} a grey-white fog appears while opening the bottle.
  4. At T0=20CT_0 = 20\,{}^\circ\mathrm{C} a blue fog appears while opening the bottle. [0,7 т.]

C.3 Give the numerical value of HcH_\mathrm{c} if the external temperature is T0=6CT_0 = 6\,{}^\circ\mathrm{C}. [1,3 т.]

Решение

Покажи официалното решение


Оригинал в Архива: IPhO_2025_Q3.pdf · официални решения: IPhO_2025_S3.pdf