Logo
Logo

Торнадо

A1  0.80 Считайте, что ускорение свободного падения и температура постоянны, а воздух можно считать идеальным газом. Покажите, что давление зависит от высоты $z$ следующим образом:
\begin{equation*}
P(z) = P_0 e^{-\alpha z},
\end{equation*}где $\alpha$ – некоторая константа. Найдите $\alpha$. Выразите ответ через постоянные, приведенные в части "Физические константы".

Рассмотрим небольшой объем воздуха площади $dS$ и высоты $dz$. Вдоль вертикальной оси на него действуют силы давления $P(z) \cdot dS$ вверх и $P(z+dz) \cdot dS$ вниз, а также сила тяжести $\rho g \, dV$ вниз. Т.к. ветер отсутствует и воздух покоится, эти силы уравновешены:
\begin{equation*}
P(z) \cdot dS - P(z+dz) \cdot dS - \rho g \, dS \cdot dz = 0.
\end{equation*}Сокращая на $dS \cdot dz$, получим дифференциальное уравнение на давление:
\begin{equation*}
\frac{dP}{dz} = - \rho g.
\end{equation*}Выражая плотность через давление и температуру, получим
\begin{equation*}
\frac{dP}{P} = - \frac{\mu_{\text{air}} g}{R T_0} dz.
\end{equation*}После интегрирования находим $P$:
\begin{equation*}
P = P_0 \exp \left( - \frac{\mu_{\text{air}} g}{R T_0} z \right).
\end{equation*}Сравнивая с выражением $P = P_0 e^{-\alpha z}$, находим $\alpha$:

Ответ: \begin{equation*}
\alpha = \frac{\mu_{\text{air}} g}{R T_0} = 1.17\cdot10^{-4} \text{м}^{-1}.
\end{equation*}
A2  1.00 Пусть теперь плотность воздуха постоянна и равна $\rho_{\text{air}}$. Найдите зависимость давления от высоты $P(z)$.
В такой ситуации температура зависит от высоты линейно. Найдите градиент температуры $b = \frac{dT}{dz}$.

С учетом $\rho = \rho_{\text{air}}$, аналогично предыдущему пункту, находим
\begin{equation*}
\frac{dP}{dz} = - \rho_{\text{air}} g.
\end{equation*}После интегрирования получаем
\begin{equation*}
P = P_0 - \rho_{\text{air}} g z.
\end{equation*}Выразим температуру через давление:
\begin{equation*}
T = \frac{\mu_{\text{air}}}{R \rho_{\text{air}}} \left( P_0 - \rho_{\text{air}} g z \right).
\end{equation*}Дифференцируя, находим $b$:

Ответ: \begin{equation*}
b = -\frac{\mu_{\text{air}} g}{R} = -0.0342 \text{К}/\text{м}.
\end{equation*}
A3  0.20 Используя результат пункта A2, найдите давление в точке $B$ на основании кучево-дождевого облака. Высота $h = 1~\text{км}$.

Ответ: \begin{equation*}
P_B = P_0 - \rho_{\text{air}} g h = 88.2 \text{кПа}.
\end{equation*}
B1  0.40 Покажите, что в обеих областях I и II выполняется
\begin{equation*}
\frac{dP}{dr} = \rho_{\text{air}} \frac{v^2}{r}.
\end{equation*}

Силы, действующие на объем воздуха (показано сечение плоскостью z = сonst)

Рассмотрим небольшой объем воздуха, занимающий цилиндрические координаты от $r$ до $r+dr$, от $\varphi$ до $\varphi + d\varphi$ и от $z$ до $z+dz$ (см. рис. 1). Площади участков равны $S_{AC} = r \, d\varphi \, dz$, $S_{AB} = S_{CD} = dr \, dz$, $S_{BD} = d\varphi \, dz \cdot \left(r + dr\right)$. Найдем проекцию силы, действующей на рассматриваемый объем, на ось $x$, удерживая слагаемые до третьего порядка включительно:
\begin{multline*}
dF_x = dF_{AC} - dF_{BD} + \frac{d\varphi}{2} \left( dF_{CD} + dF_{AB} \right) = P(r) \cdot S_{AC} - P(r+dr) \cdot S_{BD} + \frac{d\varphi}{2} \cdot P(r+dr/2) \cdot \left( S_{AB} + S_{CD} \right) = - r \, \frac{dP}{dr} dr \, d\varphi \, dz.
\end{multline*}Удобно перейти к элементу объема $dV = dr \cdot r\,d\varphi \cdot dz$, тогда сила перепишется как
\begin{equation*}
d F_x = - \frac{dP}{dr} dV.
\end{equation*}(Альтернативно, можно было рассмотреть малый прямоугольник и показать, что действующая на него со стороны давления сила равна $d\vec{F} = - \nabla P \, dV$)
Перейдем к ускорению. Т.к. радиальной скоростью мы пренебрегаем, то остается лишь центростремительное ускорение:
\begin{equation*}
a_x = -\frac{v^2}{r}.
\end{equation*}Записывая второй закон Ньютона в проекции на ось $x$ для рассматриваемого объема (с учетом $dm = \rho \, dV$), получаем искомое соотношение.

B2  0.50 Найдите зависимость тангенциальной скорости от радиуса $v(r)$ в области I. Выразите ответ через $v_G$ и $r_G$ (скорость и радиус на границе воронки).

Силы, действующие на объемы воздуха, радиальны, поэтому сохраняется момент импульса, откуда $r v(r) = \text{const}$. Используя условие $v(r_G) = v_G$, находим $\text{const} = r_G v(r_G)$. Тогда для скорости получаем следующее выражение:

Ответ: \begin{equation*}
v(r) = \frac{v_G r_G}{r}.
\end{equation*}
B3  0.50 Найдите скорость $v_G$.

Подставим найденную скорость в уравнение, полученное в B1:
\begin{equation*}
\frac{dP}{dr} = \rho_{\text{air}} \frac{v_G^2 r_G^2}{r^3}.
\end{equation*}После интегрирования получим
\begin{equation*}
P(r) = \text{const} - \frac{\rho_{\text{air}} v_G^2}{2} \frac{r_G^2}{r^2}.
\end{equation*}Используя условие $P_0 = P_A \approx P(\infty) = \text{const}$, получаем зависимость давления от радиуса (на нулевой высоте):
\begin{equation*}
P(r) = P_0 - \frac{\rho_{\text{air}} v_G^2}{2} \frac{r_G^2}{r^2}.
\end{equation*}Подставляя $r=r_G$, получим выражение для давления в точке $G$:
\begin{equation*}
P_G = P_0 - \frac{\rho_{\text{air}} v_G^2}{2}.
\end{equation*}С другой стороны, давление в точках $G$ и $B$ совпадает (т.к. они лежат на изобарном граничном слое), т.е. $P_G = P_B = P_0 - \rho_{\text{air}} g h$. Сравнивая это с предыдущим выражением, находим $v_G$:

Ответ: \begin{equation*}
v_G = \sqrt{2 g h}.
\end{equation*}
B4  2.00 Найдите зависимость $r_C(z)$. Выразите ответ через $r_G$, $v_G$, $z$ и $g$. Качественно нарисуйте форму торнадо в безразмерных осях $z/h$ и $r/r_G$ (где $h$ изображено на рис. 2).

Подставим в зависимость $P(r)$ из предыдущего пункта выражение для $v_G$:
\begin{equation*}
P(r) = P_0 - \rho_{\text{air}} g h \frac{r_G^2}{r^2}.
\end{equation*}С высотой давление падает как $\frac{dP}{dz} = -\rho_{\text{air}} g$ (т.к. ускорение воздуха по вертикали отсутствует), поэтому для давления в произвольной точке области I получим следующее выражение:
\begin{equation*}
P(r, z) = P_0 - \rho_{\text{air}} g \left( h \frac{r_G^2}{r^2} + z \right).
\end{equation*}На пограничном слое давление совпадает с $P_B = P_G = P_0 - \rho_{\text{air}} g h$. Сравнивая с полученным выражением для $P(r,z)$, получаем уравнение на $r_C$:
\begin{equation*}
h \frac{r_G^2}{r_C^2} + z = h.
\end{equation*}Отсюда находим $r_C$:

Ответ: \begin{equation*}
r_C = \frac{r_G}{\sqrt{1 - \frac{z}{h}}}.
\end{equation*}

Форма торнадо, полученная в B4.

B5  0.50 Большинство торнадо выглядят как воронка (радиус увеличивается при увеличении высоты), но некоторые имеют более постоянный диаметр, как труба. Если другие параметры одинаковы, какой из торнадо имеет бóльшую скорость $v_G$?

Ясно, что чем больше $h$, тем меньше $\frac{d r_C}{dz}$. Тогда более однородными по диаметру кажутся торнадо с бóльшим $h$ и, соответственно, бóльшей скоростью.

C1  1.10 Считая, что в области II ($r < r_C$) торнадо движется как твердое тело, найдите (тангенциальную) скорость $v(r)$ в этой области. Постройте график профиля скорости от $r=0$ до $r=\infty$.

Т.к. мы считаем, что в области II торнадо движется как твердое тело, то там скорость будет линейно зависеть от расстояния: $v_{II} = \omega r = \frac{v_G}{r_G} r$. В области I, согласно пункту B2, $v \propto \frac{1}{r}$. Профиль скорости показан на рис. 3.

Профиль скорости. $v_C$ и $r_C$ -- скорости и радиус точки на границе.

C2  1.20 Найдите давление в центре торнадо (точка $C$, на той же высоте, что и точка $G$).

Подставим полученное выражение для скорости в уравнение из пункта B1:
\begin{equation*}
\frac{dP}{dr} = \rho_{\text{air}} \frac{v_G^2}{r_G^2} r.
\end{equation*}Проинтегрируем:
\begin{equation*}
P = \text{const} + \frac{\rho_{\text{air}} v_G^2}{2} \frac{r^2}{r_G^2}.
\end{equation*}С учетом условия $P(r_G) = P_G = P_0 - \rho_{\text{air}} g h$ и $v_G^2 = 2 g h$ получаем
\begin{equation*}
P = P_0 - 2 \rho_{\text{air}} g h + \rho_{\text{air}} g h \cdot \frac{r^2}{r_G^2}.
\end{equation*}Соответственно, давление в центре равно

Ответ: \begin{equation*}
P_C = P_0 - 2 \rho_{\text{air}} g h = 76.5 \text{кПа}.
\end{equation*}
C3  0.50 Найдите температуру в центре торнадо (точка $C$).

Ответ: \begin{equation*}
T_C = \frac{\mu_{\text{air}}}{R \rho_{\text{air}}} \cdot \left( P_0 - 2 \rho_{\text{air}} g h \right) = -17 ^{\circ}\text{С}.
\end{equation*}
C4  0.30 Опираясь на результаты пункта C3, предположите возможный источник огромной энергии торнадо.

Низкая температура торнадо приводит к конденсации водяного пара, в результате чего выделяется большое количество энергии в виде теплоты парообразования.

D1  0.80 Найдите отношение подъемной силы, действующей на крышу, к ее весу.

Используя выражение для $P(r)$ из пункта B4 находим давление:
\begin{equation*}
P(r) = P_0 - \frac{\rho_{\text{air}} g h}{4}.
\end{equation*}Тогда разница давлений на крышу составит $\Delta P = \frac{\rho_{\text{air}} g h}{4}$, что приведет к подъемной силе\begin{equation*}
F_{\text{под}} = \Delta P \cdot S,
\end{equation*}где $S$ – площадь крыши.
Обозначив толщину крыши через $d$, найдем ее вес:
\begin{equation*}
F_{\text{вес}} = \rho g S d,
\end{equation*}откуда искомое отношение равно

Ответ: \begin{equation*}
\frac{F_{\text{под}}}{F_{\text{вес}}} = \frac{\rho_{\text{air}} h}{4 \rho d} = 3.75.
\end{equation*}
D2  0.20 Лучше открыть или закрыть окна? Для определенности считайте, что крыша может выдержать силы, в пять раз большие ее веса.

Окна лучше закрыть.