Кавитация

Решение

A1. Естественно предположить, что пузырек не схлопнется если \(p_\infty \le p_v\). Поэтому будем считать, что \(p_\infty > p_v\), а время схлопывания будет зависеть только от разности \(p_\infty - p_v\). Размерности используемых величин

\[[R_0] = \text{м}, \; [\rho] = \text{кг}/\text{м}^3, \; [p_\infty - p_v] = \text{кг}/(\text{м}\cdot \text{с}^2).\]Из них можно составить единственную комбинацию с размерностью времени

Ответ:

\[\tau \sim R_0 \sqrt{\frac{\rho}{p_\infty - p_v}} \approx 0.1 ~\text{мс}\]

A2. В равновесии давления должны быть связаны условием (давление внутри пузырька больше давления на бесконечности за счет поверхностного натяжения)

\[p_v + q = p_\infty + \frac{2 \sigma}{R},\]при \(R = R_0\) и \(q = q_0\) оно дает

\[q_0 = p_\infty - p_v + \frac{2 \sigma}{R_0} \approx 113 ~\text{кПа}.\] С другой стороны при постоянной температуре

\[q_0 R_0^3 = q R^3,\]поэтому

\[p_\infty = p_v + q_0 \frac{R_0^3}{R^3} - \frac{2 \sigma}{R}.\]Критическому значению давления отвечает минимум \(p_\infty\) как функции \(R\). Этот минимум получается при

\[R_c = R_0 \sqrt{\frac{3 q_0 R_0}{2 \sigma}} \approx 48~\text{мкм},\]соответствующее значение давления

\[p_c= p_v - \frac{4}{3} \sqrt{\frac{2 \sigma^3}{3 q_0 R_0^3}} = p_v - 2 q_0 \left( \frac{R_0}{R_c}\right)^3 \approx 332~\text{Па}.\]

Ответ:

\[q_0 = p_\infty - p_v + \frac{2 \sigma}{R_0} \approx 113 ~\text{кПа}, \quad p_c = p_v - 2 q_0 \left( \frac{R_0}{R_c}\right)^3 \approx 332~\text{Па}.\]

B1. Пусть \(v(r,t)\) – скорость элемента жидкости на расстоянии \(r\) от центра пузырька в момент времени \(t\), а \(u(t) = v(R, t)\) – скорость на границе пузырька

\[u(t) = \dot{R}(t).\]

Тогда из условия несжимаемости получим

\[v(r,t) = \frac{R^2}{r^2} u(t).\]

Первый способ. Использование законов Ньютона

Запишем второй закон Ньютона для элемента жидкости. Для этого нам нужно найти его ускорение. Сначала вычислим производную по времени:

\[\frac{\partial v}{\partial t} = \frac{\partial}{\partial t} \left( \frac{R^2}{r^2} u(t) \right) = \frac{2 R}{r^2} u^2 + \frac{R^2}{r^2} \dot{u}.\]

Скорость выбранного элемента жидкости меняется из-за двух причин. Во-первых распределение скоростей в пространстве меняется с течением времени. Во-вторых положение данного элемента жидкости в пространстве меняется с течением времени:

\[\Delta v = v(r + \Delta r, t + \Delta t) - v(r, t),\]где изменение радиальной координаты \(\Delta r = v \Delta t\). Тогда радиальное ускорение частицы жидкости

\[a = \frac{\partial v}{\partial t} + v \frac{\partial v}{\partial r} = \frac{R^2}{r^2} \dot{u} + \frac{2R}{r^2} u^2 - \frac{2 R^4}{r^5} u^2.\]

Если \(p = p(r)\) – давление на расстоянии \(r\) (здесь мы явно не указываем возможную зависимость давления от времени), второй закон Ньютона имеет вид (штрих здесь – радиальная производная) \[\rho a = - p',\]

или

\[\frac{R^2}{r^2} \dot{u} +\left( \frac{2R}{r^2} - \frac{2 R^4}{r^5} \right)u^2 = - \frac{p'}{\rho}.\]

Проинтегрируем это уравнение от \(r = R\) до \(r = \infty\) и получим

\[R \dot{u} + \frac{3}{2} u^2 = \frac{p(R) - p_\infty}{\rho}.\]

Поскольку давление внутри пузырька связано с давлением снаружи соотношением \[p = p(R) + \frac{2\sigma}{R},\]

мы получим \[R \ddot{R} + \frac{3}{2} \dot{R}^2 + \frac{2\sigma}{\rho R} = \frac{p- p_\infty}{\rho}.\] Наконец учтем, что давление в пузырьке складывается из давления пара \(p_v\) и давления воздуха, а значит \[p = p_v + q = p_v + q_0 \left( \frac{R_0}{R}\right)^{3 \gamma},\]откуда уравнение движения

\[\rho R\ddot{R} + \frac{3}{2} \rho \dot{R}^2 + \frac{2\sigma}{ R} - q_0\frac{ R_0^{3 \gamma}}{R^{3 \gamma}} = p_v - p_\infty.\]

Второй способ. Энергетический подход

Найдем кинетическую энергию движущейся жидкости. Масса жидкости на расстоянии от \(r\) до \(r + dr\) от центра пузырька равна

\[dm = 4\pi r^2 \rho dr,\]откуда ее кинетическая энергия

\[dE_k = \frac{dm }{2} v(r, t)^2 = 2 \pi \rho r^2 u^2 \frac{R^4}{r^4} dr .\]Тогда полная кинетическая энергия

\[E_k = \int_R^{\infty} 2 \pi \rho u^2 R^4 \frac{dr}{r^2} = 2 \pi \rho R^3 u^2\]Также нужно учесть поверхностную энергию

\[E_s = 4 \pi R^2 \sigma.\]Энергия жидкости и поверхности меняется за счет того, что давление воздуха внутри пузырька и жидкости на бесконечности совершают работу. За время \(dt\) радиус пузырька меняется на величину \(dR = u dt = \dot{R} dt\), поэтому газ внутри совершает работу

\[dW_1 = 4\pi R^2 p dR = 4\pi R^2 p u dt.\]

Для вычисления работы жидкости на бесконечности выберем сферу очень большого радиуса \(R_1 \gg R\). Скорость элементов жидкости на этой поверхности

\[u_1 = u \frac{R^2}{R_1^2},\]поэтому работа над жидкостью внутри этой сферы равна

\[d W_2 = - 4 \pi R_1^2 p(R_1) u_1 dt = - 4\pi R^2 p(R_1) u dt.\]Когда радиус сферы стремится к бесконечности, мы можем заменить \(p(R_1) = p_\infty\), и тогда радиус сферы исчезнет из выражения для работы. Тогда закон сохранения энергии имеет вид

\[d(E_k + E_s) = dW_1 + dW_2.\]Вычисляя дифференциалы и сокращая \(dt\), получим

\[\frac{d}{dt} \left( 2 \pi \rho R^3 \dot{R}^2 + 4 \pi \sigma R^2\right) = 4\pi R^2 \dot{R}(p - p_\infty),\]откуда (после деления на \(4\pi\) и сокращения общего множителя \(R^2 \dot{R}\))

\[\frac{3}{2} \rho \dot{R}^2 + \rho R \ddot{R} + \frac{2 \sigma}{R} = p- p_\infty.\]Подставляя выражения для \(p\), получим то же уравнение

Ответ:

\[\rho R\ddot{R} + \frac{3}{2} \rho \dot{R}^2 + \frac{2\sigma}{ R} - q_0\frac{ R_0^{3 \gamma}}{R^{3 \gamma}} = p_v - p_\infty\]

B2. Используя то, что \(\ddot{R} \to 0\) и \(R \to \infty\), получим

\[\dot{R}(\infty) = \sqrt{\frac{2 (p_v -p_\infty)}{3 \rho}} = \sqrt{\frac{2 p_v }{3 \rho}} \approx 1.24~\text{м}/\text{с}.\]Из уравнения движения найдем начальное ускорение

\[\ddot{R}(0) = \frac{q_0 + p_v - 2\sigma/R_0}{\rho R_0} = \frac{p^-_\infty}{\rho R_0} \approx 10^7~\text{м}/\text{с}^2,\]и поэтому время установления скорости можно оценить как

\[t \approx \frac{\dot{R}(\infty)}{\dot{R}(0)} \approx 0.1~\text{мкс}.\]

Ответ:

\[\dot{R}(\infty) =\sqrt{\frac{2 p_v }{3 \rho}} \approx 1.24~\text{м}/\text{с}, \quad t \approx \frac{\dot{R}(\infty)}{\dot{R}(0)} \approx 0.1~\text{мкс}\]

B3. Умножим уравнение из B1 на \(2 R^2 \dot{R}\), получим

\[\rho \frac{d}{dt}\left[ R^3 \dot{R}^2\right]+ 4 \sigma R \dot{R} - 2 q_0 R_0^{3 \gamma} R^{2 - 3 \gamma}\dot{R} = 2 (p_v - p_\infty) R^2 \dot{R}.\]Проинтегрируем это уравнение от \(R(0) = R_0\) и \(\dot{R}(0) = 0\) до \(R(t) = R\) и \(\dot{R}(t)\), получим

\[\rho R^3 \dot{R}^2 + 2 \sigma (R^2 - R_0^2) + \frac{2}{3 (\gamma -1)} q_0 R_0^{3 \gamma} \left( R^{3 - 3\gamma} - R_0^{3 - 3\gamma}\right) = \frac{2}{3} (p_v - p_\infty)(R^3 - R_0^3),\]

или

\[\rho \dot{R} ^2 = \frac{2 (p_v -p_\infty)}{\rho} - \frac{2 \sigma}{R} + \frac{2 R_0^3}{R^3} \left( \frac{\sigma}{R} + \frac{p_\infty - p_v}{3} + \frac{q_0}{3 (\gamma - 1)} - \frac{q_0 R_0^{3 \gamma -3}}{3 (\gamma -1) R^{3 \gamma -3}}\right).\]Это же соотношение является законом сохранения энергии, записанным в интегральной форме. Чтобы получить его непосредственно из закона сохранения энергии, подставим в выражение для работы из B1 явную зависимость давления от радиуса:

\[dW = dW_1 + dW_2 = 4 \pi R^2 dR \left( p_v + q_0 \frac{R_0^{3 \gamma}}{R^{3 \gamma}} - p_\infty\right).\]Интегрируя от \(R_0\) до \(R\), найдем полную работу

\[W = \frac{4\pi}{3} (R^3 - R_0^3)(p- p_\infty) + \frac{4\pi q_0 R_0^{3 \gamma}}{3 - 3 \gamma} \left( R^{3 - 3 \gamma} - R_0^{3- 3 \gamma} \right).\]Она равна изменению энергии жидкости и поверхностной энергии

\[\Delta (E_k + E_s) = 2\pi \rho R^3 \dot{R}^2 + 4 \pi \sigma (R^2 - R_0^2) = W.\]Подставляя выражения для работы и деля обе части на \(2\pi R^3\), получаем прежнее выражение для \(\rho \dot{R}^2\).

Когда радиус пузырька достигнет минимального значения, \(\dot{R} = 0\). Если считать, что при этом \(R \ll R_0\), то можно учитывать только наиболее быстро растущие слагаемые. Поскольку слагаемое с \(1/R^3\) растет быстрее остальных, для оценки можно считать, что выражение в круглых скобках в формуле для \(\rho \dot{R}^2\) должно быть равно нулю:

\[\frac{\sigma}{R} + \frac{p_\infty - p_v}{3} + \frac{q_0}{3 (\gamma - 1)} - \frac{q_0 R_0^{3 \gamma -3}}{3 (\gamma -1) R^{3 \gamma -3}} =0,\]откуда

\[\left( \frac{R_0}{R}\right)^{3 \gamma -3} = \frac{\gamma -1}{q_0} \left( p_\infty - p_v + \frac{3 \sigma}{R_0}\right) +1.\]Аналогично A2 найдем начальное значение давления воздуха внутри пузырька:

\[q_0 = p_\infty^- - p_v +\frac{2 \sigma}{R_0} = 13.8~ кПа.\]Тогда численное значение минимального радиуса

\[R \approx 2.88~ мкм.\]Точное решение уравнения \(\dot{R}^2 = 0\) дает \(R = 2.95 ~ мкм\).

Ответ:

\[R = 2.95 ~ мкм\]

B4. При \(R \ll R_0\) из уравнения движения следует (наиболее быстро растущее слагаемое в выражении для \(\dot{R}^2\) ведет себя как \(1/R^3\))

\[R' \sim - R^{-3/2}.\]Подставим в это уравнение решение вида \(R \sim (T- t)^\alpha\), получим

\[(T- t)^{\alpha - 1} \sim (T- t)^{-3 \alpha/2},\]откуда

\[\alpha - 1 = - \frac{3 \alpha}{2},\]

Ответ:

\[\alpha = \frac{2}{5}.\]

Может показаться, что если использовать вместо закона сохранения непосредственно уравнение движения, получится другой показатель степени. Действительно, подставив в уравнение движения решение вида \(R = A (T-t)^\alpha\), получим

\[A (T-t)^\alpha A \alpha (\alpha -1) (T- t)^{\alpha - 2} + \frac{3}{2} A^2 \alpha^2 (T- t)^{2 \alpha} + \frac{2 \sigma}{\rho A} (T-t)^{-\alpha} = \frac{p_v - p_\infty}{\rho}.\]Тогда первые два слагаемых пропорциональны \((T- t)^{2\alpha-2}\), третье слагаемое порядка \((T- t)^{-\alpha}\), правая часть не зависит от времени, а слагаемое с \(q_0\) по условию не учитываем. Тогда чтобы первые три слагаемых могли сократиться, нужно потребовать \(2 \alpha - 2 = -\alpha\), \(\alpha = 2/3\).

Однако это рассуждение работает только при условии, что коэффициент перед первыми двумя слагаемыми не равен нулю

\[A^2\left( \alpha(\alpha -1) + \frac{3}{2} \alpha^2 \right) = A^2 \left( \frac{5}{2} \alpha^2 - \alpha\right).\]Однако он обращается в 0 как раз при правильном значении показателя степени \(\alpha = 2/5\), так что никакого противоречия нет. Последующие слагаемые сократятся за счет поправок к лидирующему решению:

\[R(t) = A(T-t)^{2/5} (1 + f (T- t) + \dots),\]при правильном выборе постоянной \(f\).

B5. Для описания поведения пузырька введем безразмерную переменную \(x\) согласно \(R = R_0(1 + x)\), тогда уравнение движения

\[\rho R_0^2 (1 + x)\ddot{x} + \frac{3}{2} R_0^2 \dot{x}^2 + \frac{2 \sigma}{R_0(1 + x)}- \frac{q_0}{(1 + x)^{3 \gamma}} = p_v - p_\infty.\]Считая колебания малыми, оставим только линейный по \(x\) слагаемые. Тогда в первом слагаемом можно пренебречь \(x\), вторым слагаемым можно пренебречь полностью, а оставшиеся разложить до первого порядка по \(x\). Не зависящие от \(x\) слагаемые сократятся в силу условия равновесия пузырька. Окончательно получим

\[\rho R_0^2 \ddot{x} + \frac{2 \sigma}{R_0} (1 - x) - q_0 (1 - 3 \gamma x) =p_v - p_\infty, \quad \frac{2\sigma}{R_0} - q_0 = p_v - p_\infty,\]\[\rho R_0^2 \ddot{x} + \left(3 \gamma q_0 - \frac{2 \sigma}{R_0} \right) x= 0.\]Тогда пузырек будет устойчив при условии \(3 \gamma q_0 > 2 \sigma /R_0\), а частота колебаний (при вычислении считаем, что пузырек находится в равновесии при радиусе из условия и давлении \(p_\infty = 101~ кПа\), откуда \(q_0 \approx 100~ кПа\))

Ответ:

\[f = \frac{1}{2\pi} \sqrt{\frac{3 \gamma q_0 - 2 \sigma/R_0}{\rho R_0^2}} \approx 32.6 ~ кГц\]

B6. Будем использовать то же линеаризованное уравнение колебаний пузырька, но теперь учтем, что давление на бесконечности меняется по закону из условия

\[\rho R_0^2 \ddot{x} + \left(3 \gamma q_0 - \frac{2 \sigma}{R_0} \right) x= A \sin\left(\frac{2\pi f}{c}a\right) \sin(2\pi ft).\]Установившееся решение имеет вид \(x = B \sin (2\pi f t)\), амплитуду решения можно найти, подставив его в уравнение. Для упрощения записи ответов выразим коэффициент перед \(x\) через \(f_0^2\):

\[\rho R_0^2 4\pi^2 (f_0^2 - f^2) B \sin (2\pi f t) = A \sin\left(\frac{2\pi f}{c}a\right) \sin(2\pi ft),\]откуда решение имеет вид

\[x(t) = \frac{A \sin\left(\frac{2\pi f} {c}a\right)}{4\pi^2 \rho R_0^2 (f_0^2 -f^2)} \sin (2\pi f t).\]Давление меняется с координатой, поэтому на пузырек действует сила, аналогичная силе Архимеда

\[F_x = - \rho g_x V,\]где \(g_x\) – эффективное ускорение свободного падения, которое можно найти, зная градиент давления:

\[\rho g_x = \frac{\partial p}{\partial x}.\]Тогда среднее по времени значение силы

\[F = -\left\langle V \frac{\partial p}{\partial x}\right\rangle.\]Объем пузырька

\[V = \frac{4 \pi }{3 } R^3 = \frac{4 \pi }{3 } R^3_0 (1 + x)^3 \approx \frac{4 \pi }{3 } R^3_0(1 + 3x),\]постоянное слагаемое сократится при усреднении, тогда получим

\[F = - 4\pi R_0^3 \left\langle x \frac{\partial p}{\partial x}\right\rangle = - 4\pi R_0^3 \frac{A \sin\left(\frac{2\pi f} {c}a\right)}{4\pi^2 \rho R_0^2 (f_0^2 -f^2)} \left\langle\sin (2 \pi ft) \frac{2 \pi f}{c} A \cos\left( \frac{2\pi f}{c} f\right) \sin(2\pi ft)\right\rangle.\]Используя среднее от квадрата синуса \(\langle \sin^2 (2\pi ft) = 1/2\) и формулу для синуса двойного угла, получим окончательно

Ответ:

\[F = - \frac{R_0 f A^2 }{2 \rho c (f_0^2 - f^2)} \sin \left( \frac{4\pi f}{c} a\right)\]

C1. Концентрация растворенного воздуха на большом расстоянии до пузырька

\[u_i = H p_\infty,\]а концентрация вблизи пузырька

\[u = H q > u_i,\]поэтому возникнет диффузионный поток растворенного воздуха, направленный от пузырька. Из-за этого масса воздуха внутри пузырька будет убывать, и со временем пузырек полностью растворится.

Оценим поток воздуха из пузырька. Для оценки примем, что концентрация растворенного воздуха существенно изменяется на расстоянии порядка радиуса пузырька, тогда полный поток

\[J = 4 \pi R^2 \kappa \frac{u - u_i}{R} = 4\pi R \kappa H(q - p_\infty).\]Поскольку давлением водяного пара можно пренебречь, давление воздуха в пузырьке

\[q = p_\infty + \frac{2 \sigma}{R}.\]Тогда окончательное выражение для потока массы из пузырька

\[J = 4 \pi R \kappa H \frac{2\sigma}{R} = 8 \pi \kappa H \sigma.\]

Видим, что это выражение не зависит от радиуса, поэтому для оценки времени можно просто разделить полную массу воздуха в пузырьке на поток. Эта масса равна

\[m = \frac{4\pi}{3} R^3_0 \frac{\delta}{p_\infty} \left( p_\infty + \frac{2 \sigma}{R_0}\right).\]Здесь \(\delta\) – плотность воздуха при атмосферном давлении \(p_\infty\).

Тогда время жизни пузырька

\[t \approx \frac{m}{J} = \frac{\frac{4\pi}{3} R_0^3 \delta \left( 1 + \frac{2 \sigma}{p_\infty R_0}\right)}{8 \pi \kappa H \sigma} = \frac{R_0^3 \delta}{6 \kappa \sigma H} \left( 1 + \frac{2 \sigma}{p_\infty R_0} \right)\]

Ответ:

\[t \approx \frac{R_0^3 \delta}{6 \kappa \sigma H} \left( 1 + \frac{2 \sigma}{p_\infty R_0} \right) \approx 7~ с\]

В основной части решения мы использовали оценку для потока воздуха из пузырька

\[J = 4 \pi R^2 \kappa \frac{u - u_i}{R}.\]

Из качественных соображений невозможно получить коэффициент в формуле. Рассмотрев решение уравнения диффузии покажем, что коэффициент действительно выбран правильно (на олимпиаде такое вычисление в действительности требовалось от участников).

Рассмотрим сферу радиуса \(r\) и толщины \(dr\). Скорость изменения массы воздуха в этой сфере можно выразить через производную концентрации:

\[\frac{dm}{dt} = 4\pi r^2 dr \frac{\partial u}{\partial t}.\]С другой стороны, ее можно выразить как разность втекающих и вытекающих потоков

\[\frac{dm}{dt} = 4\pi (r + dr)^2 J(r + dr) - 4\pi r^2 J(r) = 4\pi (r + dr)^2 \kappa\frac{\partial u}{\partial r} (r + dr) - 4\pi r^2 \kappa \frac{\partial u}{\partial r} (r).\]Используя соотношение

\[\frac{\partial u}{\partial r} (d + dr) = \frac{\partial u}{\partial r} (r) + \frac{\partial ^2 u}{\partial r^2}(r) dr,\]получим

\[\frac{dm}{dt} = 4\pi r^2 \kappa dr \frac{\partial^2 u}{\partial r^2} + 8 \pi r \kappa dr \frac{\partial u} {\partial r}.\]Объединив эти два уравнения, получим

\[\frac{\partial u}{\partial t} = \kappa \frac{\partial^2 u}{\partial r^2} + \frac{2 \kappa}{r} \frac{\partial u} {\partial r}.\]В начальный момент времени во всем объеме воды концентрация растворенного воздуха одинакова и равна \(u = u_i\). Во все моменты времени концентрация воздуха в непосредственной окрестности пузырька равна \(H q\):

\[\begin{cases} u(r, 0) = u_i \quad \text{ при } r >R,\\ u(R, t) = Hq, \quad \text{ при } t >0. \end{cases}\]

Введем новую переменную \(v = r(u - q H)\), и новый параметр времени \(\tau = \kappa t\). Тогда уравнение диффузии перейдет в

\[\frac{\partial v}{\partial \tau} = \frac{\partial ^2 v}{\partial r^2},\]с граничными условиями

\[\begin{cases} v(r, 0) = r(u_i -q H) \quad \text{ при } r >R,\\ v(R, t) = 0 \quad \text{ при } \tau >0. \end{cases}\]Также введем переменную \(\xi = r - R\) и функцию \(w(\xi, \tau) = v(R + \xi, \tau)\), тогда

\[\frac{\partial w}{\partial \tau} = \frac{\partial ^2 w}{\partial \xi^2},\]с граничными условиями

\[\begin{cases} w(\xi, 0) = (R + \xi)(u_i -q H) \quad \text{ при } \xi > 0,\\ w(0, t) = 0 \quad \text{ при } \tau >0. \end{cases}\]Решение этого уравнения имеет вид

\[w(\xi, \tau) = \frac{u_i - qH}{\sqrt{4\pi \tau}} \int_0^{\infty} \left(e^{- (\xi - \eta)^2/(4 \tau)} - e^{- (\xi + \eta)^2/(4 \tau)}\right) (\eta + R) d\eta.\]Поскольку

\[u(r,t) = q H + \frac{v(r, \kappa t)}{r} = qH + \frac{w(r - R, \kappa t)}{r},\]производная имеет вид

\[\frac{\partial u}{\partial r} = \frac{1}{r} \frac{\partial w}{\partial \xi} - \frac{w}{r^2},\]откуда

\[\left. \frac{\partial u}{\partial r}\right|_{r = R} = \frac{1}{R} \left. \frac{\partial w}{\partial \xi} \right|_{\xi = 0} = (u_i - qH) \left( \frac{1}{R} + \frac{1}{\sqrt{\pi \kappa t}}\right).\]Здесь знаменатель второго слагаемого порядка характерного размера диффузионной области, в которой существенно меняется концентрация. Мы считаем, что размер этой области много больше размера пузырька, поэтому вторым слагаемым можно пренебречь. С учетом этого поток воздуха через поверхность пузырька

\[J = - 4\pi R^2 \kappa \left. \frac{\partial u}{\partial r}\right|_{r = R} = 4\pi R^2 \kappa \frac{u - u_i}{R}.\]Это и есть выражение, которое мы использовали в основной части решения.

C2. Условия равновесия имеют вид

\[p_v + q = p_\infty \pm \frac{2\sigma}{R}, \quad u = Hq.\]Выбор знака зависит от направления выпуклости границы. Пусть начальное давление воздуха в пузырьке равно \(q\). При больших значениях давления воздуха в пузырьке \(q\), это давление будет убывать, а пузырек будет уменьшаться. В некоторый момент система придет в равновесие. Объем воздуха не исчезнет полностью, если поверхность воды выпуклая (и давление внутри может быть меньше давления на бесконечности, так что диффузия остановится):

Ответ:

\[\pi + \alpha < 2 \theta.\]

Перевод: pho.rs.