Logo
Logo

Всплывающие пузыри

При наблюдении за воздушными пузырьками, всплывающими в жидкости, можно заметить, что их форма очень сильно зависит от объёма воздуха внутри пузырьков и может быть как близкой к сферической, так и напоминать тонкий диск. Данная задача посвящена изучению форм воздушных пузырьков, всплывающих в жидкости. При решении задачи примем следующую модель явления: 

  •  если не указано обратного, все движущиеся тела имеют форму, слабо отличающуюся от сферической; 
  •  сжимаемостью и вязкостью жидкости можно пренебречь; 
  •  гравитационным полем в частях A и B можно пренебречь.

Часть А. Хорошо забытое старое(5.0 балла)

Рассмотрим шар радиусом $R$, движущийся со скорость $\vec{v}_0$ в бесконечном объёме несжимаемой жидкости плотностью $\rho$ с нулевой вязкостью. Шар приводит жидкость в движение. Напомним, что гравитационным полем можно пренебречь. Движение жидкости характеризуется величиной $\mu$, называемой присоединённой массой. С присоединённой массой связаны следующие характеристики движущейся жидкости: 

  1.  Если шар движется в жидкости со скоростью $\vec{v}_0$, то кинетическая энергия $E_k$ жидкости, окружающей шар, определяется выражением: $$E_k=\cfrac{\mu v_0^2}{2}.  $$ 
  2.  Если шар движется в жидкости со скоростью $\vec{v}_0$, то импульс $\vec{p}$ жидкости, окружающей шар, определяется выражением $$ \vec{p}=\mu\vec{v_0}. $$ 

Пусть $\vec{r}$ — радиус-вектор, проведенный из центра шара в произвольную точку пространства, а $\vec{v}(\vec{r})$ — скорость соответствующего кусочка жидкости в лабораторной системе отсчёта. При решении задачи вам понадобится теорема Томсона, утверждающая, что если изначально неподвижную несжимаемую жидкость с нулевой вязкостью привести в движение, то циркуляция вектора скорости жидкости по любому замкнутому контуру, не пересекающему поверхности движущихся в жидкости тел, обращается в ноль, что можно записать следующим образом: $$\oint_l\vec{v}d\vec{l}=0{.} $$ В пунктах A1 - A9 будем решать задачу в в системе отсчёта, связанной с шаром. Скорость жидкости в этой системе отсчёта будем обозначать как $\vec{v}_{\text{отн}}$.

A1  0.30 Чему равен поток вектора относительной скорости жидкости через произвольную замкнутую поверхность, не пересекающую поверхность шара? Ответ обоснуйте.

 Примечание: потоком вектора относительной скорости жидкости через замкнутую поверхность называется величина $\Phi$, определяемая выражением: \begin{equation*} \Phi=\oint_S\vec{v}_\text{отн}\bigl(\vec{r}\bigr)d\vec{S}{.} \end{equation*}

A2  0.10 Чему равна нормальная компонента скорости жидкости ${v}_{\text{отн},n}$ в системе отсчёта, связанной с шаром?

A3  0.10 Чему равна относительная скорость жидкости $\vec{v}_{\infty}$ на бесконечном удалении от шара?

A4  0.40 Рассмотрим сверхпроводящее тело, помещенное в однородное магнитное поле с индукцией $\vec{B}_0$. Напишите выражение для циркуляции магнитного поля $\vec{B}(\vec{r})$ по контуру, не пересекающему поверхность сверхпроводника, и выражение для потока магнитного поля через замкнутую поверхность, не пересекающую поверхность сверхпроводника. Также укажите граничные условия магнитного поля на границе тела и на бесконечном удалении от него.

A5  0.80 Поместим сверхпроводящий шар радиусом $R$ в однородное магнитное поле с индукцией $\vec{B_0}$. Найдите индукцию магнитного поля $\vec{B}(\vec{r})$ во всем пространстве вне шара. Ответ выразите через $\vec{B_0}, \vec{r}$ и $R$.

A6  0.20 Пользуясь результатами предыдущих пунктов, найдите зависимость $\vec{v}_{\text{отн}}(\vec{r})$.

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

A7  0.40 Рассмотрим стационарное движение несжимаемой жидкости с нулевой вязкостью по тонкой трубке. Пусть $p$ — давление в жидкости, $\varphi$ — потенциальная энергия единицы массы жидкости, а $v$ — скорость течения в трубе.
Используя законы сохранения энергии и массы, покажите, что для всех точек трубки выполняется соотношение (уравнение Бернулли): $$\varphi + \frac{p}{\rho} + \frac{ v_\text{отн}^2}{2} = const$$

Пусть $\theta$ — угол между радиус-вектором $\vec{r}$ и скоростью $\vec{v_0}$.

A8  0.60 Считая, что давление на бесконечном удалении от шара равно $p_0$, найдите зависимость $p(\theta)$ на поверхности шара. Выразите ответ через $p_0$, $\rho$, $v_0$, $\theta$.

A9  0.40 Вычислите равнодействующую сил давления, действующих на поверхность шара. Выразите ответ через $\vec{v_0}$, $R$, $\rho$.

Вернёмся в лабораторную систему отсчёта. В лабораторной системе отсчёта скорость жидкости в точке с радиус-вектором $\vec{r}$ равна $\vec{v}(\vec{r})$

A10  0.10 Найдите распределение скорости жидкости $\vec{v}(\vec{r})$ в лабораторной системе отсчёта.

A11  0.80 Чему равна кинетическая энергия жидкости $E_k$ в лабораторной системе отсчёта? Выразите ответ через $\rho$, $R$, $v_0$.

A12  0.20 Найдите присоединённую массу $\mu$. Выразите ответ через $\rho$, $R$.

Если вы не смогли решить предыдущий пункт, в дальнейшем можете считать $\mu = \frac{4}{3}\pi R^3\rho$.

A13  0.30 Рассмотрим шар, движущийся со скоростью $\vec{v}_0$ и ускорением $\vec{a}_0$. Используя выражение для импульса $\vec{p}$ жидкости, получите выражение для силы $\vec{F}$, действующей на шар со стороны жидкости. Ответ выразите через $\rho$, $R$ и $\vec{a}_0$.

Поскольку пузырь не является твёрдым телом, в результате движения он будет деформироваться. В дальнейшем считайте его деформацию малой.

Рассмотрим движение пузыря с коэффициентом поверхностного натяжения $\sigma$ со скоростью $v_0$. Считайте применимыми все выражения, полученные в предыдущих пунктах. Давление внутри пузыря всюду одинаково и равно $p_{in}$.

A14  0.30 Получите зависимость кривизны $K$ поверхности пузыря от угла $\theta$. Используйте $p_{in}$, $p_0$, $\rho$, $\theta$ и $v_0$.

 Примечание: кривизной $K$ поверхности называется сумма обратных главных радиусов кривизны: $$ K=\frac{1}{R_1}+\frac{1}{R_2}$$

Она связана с Лапласовым давлением как:

$$\Delta p = \sigma K$$

 

Часть B. Влияние скорости на деформацию пузыря(2.5 балла)

Рассмотрим эллипсоид, получаемый из эллипса с большой полуосью $a$ и малой полуосью $b$ вращением относительно малой оси. Ось вращения обозначим за $y$.

Эту фигуру можно получить сжатием сферы радиусом $a$ вдоль оси $y$ с коэффициентом $k$, где $k = b/a$. Найдём кривизну поверхности полученной фигуры в произвольной точке $A'$, которой до сжатия соответствовала точка $A$ и угол $\theta$, отсчитываемый от оси $y$.

Можно показать, что в плоскости $(x, y)$ радиус кривизны эллипсоида равен:
$$R_1(\theta)= \frac{a}{k} \cdot (\cos^2{\theta} + k^2 \sin^2{\theta})^\frac{3}{2},$$
а в плоскости, содержащей нормаль к эллипсоиду в точке $A'$ и перпендикулярной плоскости $(x, y)$, радиус кривизны эллипсоида равен:
$$R_2(\theta) = \frac{a}{k} \sqrt{\cos^2{\theta} + k^2 \sin^2{\theta}}. $$

Так как $K(\theta)$ в части А было выведено в предположении, что капля является сферой, то отклонения от сферы должны быть маленькими. Это означает, что $k = 1-\varepsilon$, где $\varepsilon \ll1$. Также считайте, что $a \approx R$, где $R$ — радиус пузыря, который мы использовали в части А.

B1  0.60 С учётом малости $\varepsilon$ величины $R_1$ и $R_2$ можно представить в виде:
\begin{equation*}
R_1(\theta)\approx R(1+A_1\varepsilon-B_1\varepsilon\sin^2\theta)\qquad R_2\approx R(1+A_2\varepsilon-B_2\varepsilon\sin^2\theta){.}
\end{equation*}
Найдите $B_1$ и $B_2$. Далее слагаемыми $A_1\varepsilon$ и $A_2\varepsilon$ можно пренебречь.

Если вы не смогли решить данный пункт, то далее можно считать, что $B_1=B_2=1$. К потере последующих баллов это не приводит.

B2  0.20 Найдите зависимость кривизны эллипсоида $K(\theta)$ в первом порядке по $\varepsilon$. Выразите ответ через $R$, $\theta$.

B3  0.60 Чему равна величина $\varepsilon$ для пузыря, рассмотренного в пункте A14? Ответ выразите через скорость пузыря $v_0$, радиус пузыря $R$, коэффициент поверхностного натяжения $\sigma$ и плотность жидкости $\rho$.

B4  0.30 Чему равен эксцентриситет $e$ пузыря, рассмотренного в пункте A14? Ответ выразите через скорость пузыря $v_0$, радиус пузыря $R$, коэффициент поверхностного натяжения $\sigma$ и плотность жидкости $\rho$. 

B5  0.80 Сферический пузырь с начальным диаметром $d = 2\text{мм}$ поднимается в воде плотностью $\rho = 10^3 \text{кг}/\text{с}$ и поверхностным натяжением $\sigma = 72 \cdot 10^{-3} ~\text{Н}/\text{м}$. Проведя необходимые измерения, вычислите скорость пузыря на последних трёх фотографиях. Также вычислите ускорение пузыря $a$.

Часть C. Динамика движения пузыря(1.5 балла)

В реальной жизни ускорение свободного падения $g$ создаст в сосуде вертикальный градиент давления, в следствие которого появится сила Архимеда $F_A$, разгоняющаяя пузырь вверх. Изначально пузырь является сферическим, но в процессе разгона, согласно теории из предыдущего пункта, он начинает сплющиваться.

Второй закон Ньютона для почти сферического пузыря будет выглядеть следующим образом:
$$m \vec{a} = m\vec{g} + \vec{F_A} + \vec{F_\mu} + \vec{F_{\text{тр}}},$$
где $m$ - масса пузыря, которой можно пренебречь, $\vec{F_A}$ — сила Архимеда, $\vec{F_\mu}$ — сила, действующая со стороны воды на пузырь, связанная с присоединённой массой воды, $\vec{F_\text{тр}}$ — сила трения, возникающая из-за эффектов турбулентности. В нашей модели $\vec{F_\text{тр}}$ хорошо описывается формулой $\vec{F_\text{тр}} = -\pi R^2 \rho v \vec{v}$, где $\vec{v}$ -скорость пузыря.

Скорость пузыря изначально равна нулю.

C1  0.60 Чему равно его ускорение $a$ в начальный момент времени? Ответ выразите через $g$.

C2  0.30 Найдите установившуюся скорость пузыря $v_{\max}$. Чему равно $\varepsilon$ при установившейся скорости? Ответы выразите через $R$, $g$, $\rho$ и $\sigma$.

С3  0.60 Рассчитайте численные значения $v_{\max}$ и $\varepsilon$ для $\rho = 10^3~ {\text{кг}}/{\text{м}^3}$, $\sigma = 72 \cdot 10^{-3} ~\text{Н}/\text{м}$, $g =10 ~\text{м}/ \text{с}^{2}$ и $R = 1~\text{мм}$. Применима ли выбранная модель для указанного радиуса?

Часть D. Крутой медузоподобный пузырь(2.0 балла)

В больших по объёму пузырях сила поверхностного натяжения слабее, чем в маленьких, поэтому их форма при достижении постоянной скорости далека от сферической.

В общем случае форма пузерей, приведенных на рисунке выше, моделируется при помощи компьютера. Однако в предыдущих частях мы смогли изучить форму, близкую к сферической: это пузыри $a$ и $b$. В этой части мы изучим модель, применимую к пузырям $m$ и $n$.

Опишем пузырь как сегмент сферы c радиусом основания $r$. Радиус кривизны пузыря, то есть радиус сферы, постоянен и равен $R$. На фотографиях $m$ и $n$ видно, что пузыри почти плоские, а значит можно считать, что $r \ll R$. Точку на поверхности пузыря будем характеризовать малым углом $\theta$, отсчитываемым от центра кривизны. Давление воздуха внутри пузыря равно $p_{in}$.

D1  0.20 Докажите, что давление воды в точках на поверхности пузыря постоянно.

Линии тока воды схематично обозначены на рисунке. Будем считать, что распределение скоростей жидкости $\vec{v}_{\text{отн}}(\vec{r})$ над пузырём точно такое же, как и в случае движения твёрдого шара радиусом $R$ со скоростью $v_0$. Ускорение свободного падения $g$.

Пусть пузырь всплывает в глубоком водоёме. Давление на поверхности водоёма равно атмосферному $p_0$. Расстояние от поверхности водоёма до точки $A$ равно $h \gg R$. 

D2  1.00 Найдите зависимость давления над поверхностью пузыря $p(\theta)$, используя $\rho$, $v_0$, $g$, $R$, $\theta$ и $h$.

D3  0.80 Получите выражение для установившейся скорости пузыря $v_0$. Рассчитайте значение для $R = 30~\text{см}$ и $g = 10~\text{м}/\text{с}^2$.