Logo
Logo

Прилипание струи жидкости

A1  1.30 Экспериментально изучите зависимость $z(x)$. Сделайте не менее 15 измерений. Подробно опишите методику.

Для получения зависимости $z(x)$ необходимо снять координаты струи. Это можно сделать, поместив миллиметровку параллельно плоскости движения воды и ручкой проставить точки.

A2  0.80 Определите $U_{0}$.

Получим линеаризацию для зависимости $z(x)$:
$$x(t) = U_0\cdot t \\
z(t) = \frac{gt^2}{2} \\
z = \frac{gx^2}{2U_0^2}$$
С помощью данной линаеризации мы можем определить $U_0$ с помощью углового коэффициента зависимости $z(x^2)$.

A3  0.50 Определите площадь внутреннего сечения $A_0$ насадки и диаметр струи на выходе из нее $D_j$.

С помощью мерного цилиндра и секундомера можно измерить объёмный расход $Q$ при выбранной скорости.
$$A = \frac{Q}{U_0} = 1.96 \; мм^2$$
$$D = 2\sqrt{\frac{A}{\pi}} = 1.58 \;мм$$

B1  1.50 Экспериментально получите зависимость $z(X)$, сделайте не менее 30 измерений. Прочитайте предупреждение ниже.

Сделаем установку, которая будет создавать необходимую спираль.
Поместим внутрь пробирки миллиметровку и снимем зависимость по этой миллиметровке.

При снятии зависимости будем учитывать, что внешний и внутренний диаметры пробирки различаются, поэтому скорректируем значения $x$:

$$x_{ист} = \frac{D_{внешн}}{D_{внутр}}x$$

B2  0.80 Постройте график зависимости $z(X)$.

B3  1.30 Рассчитайте $\psi (s)$ и постройте график этой зависимо-
сти.

$\psi$ можно определить несколькими способами:

1.) С помощью построения сглаживающей кривой на графике $z(x)$ и определение $\psi$ с помощью касательных.
2.) Определение $\psi$ с помощью двух соседних точек: $ctg(\psi) = \frac{z_i-z_{i-1}}{x_i-x_{i-1}}$

$s$ определяется итеративно:
$s_{i+1} = s_i + \sqrt{(z_{i+1}-z_i)^2+(x_{i+1}-x_i)^2}$

С1  0.80 Для всех пробирок измерьте $\lambda (R)$ при одинаковых начальных условиях.

При измерении $\lambda$ на различных пробирках важно сохранять одинаковые условия измерения, в них входят $\psi_0$, скорость струи при попадании на пробирку, степень прижатости струи к пробирки, а так же метод измерения $\lambda$.

C2  0.50 Поставьте точки, соответствующие экспериментальным данным, полученным в предыдущем пункте, на графике $z(X)$. Основываясь на этом сделайте вывод относительно справедливости утверждения: «Радиус пробирки $R$ не влияет на зависимость $z(X)$ при одинаковых начальных условиях $(U_{0}, \psi_0)$».

При нанесении точек, соответствующих различным пробиркам, на изначальный график $z(x)$ можно убедится, что с точностью до погрешности они попадают на зависимость. Это говорит о том, что радиус пробирки не влияет на зависимость $z(x)$ при одинаковых начальных условиях.

D1  2.00 Получите коэффициенты $a$, $b$ и $c$ для параболического приближения зависимости $\psi (s) = as^2 + bs + c$.

Продифференциируем зависимость $\sigma(a,b,c,{s_i,\psi_i})$ по $a,b,c$:

$$\left\{\begin{array}{l}\frac{\partial \sigma}{\partial a}=\sum 2\left(a s_{i}^{2}+b s_{i}+c-\psi_{i}\right) x_{i}^{2}=0 \\ \frac{\partial \sigma}{\partial b}=\sum 2\left(a s_{i}^{2}+b s_{i}+c-\psi_{i}\right) x_{i}=0 \\ \frac{\partial \sigma}{\partial c}=\sum 2\left(a s_{i}^{2}+b s_{i}+c-\psi_{i}\right)=0\end{array}\right.$$

Упростим данную систему:

$$\left\{\begin{array}{l}a \sum s_{i}^{4}+b \sum s_{i}^{3}+c \sum s_{i}^{2}-\sum s_{i}^{2} \psi _{i}=0 \\ a \sum s_{i}^{3}+b \sum s_{i}^{2}+c \sum s_{i}-\sum s_{i} \psi _{i}=0 \\ a \sum s_{i}^{2}+b \sum s_{i}+c\cdot n-\sum \psi _{i}=0\end{array}\right.$$

Посчитав все необходимые суммы решим её.

D2  0.70 Допустим, что по плоской поверхности твердого тела стационарно и однородно течет слой воды высотой $H$. Можно считать, что разность давлений между точками, лежащими на разных торцах и одной линии тока не зависят от этой линии тока, т.к. иначе они бы искривились. Найдите распределение скорости $u(h)$ внутри этого слоя, где $h$ — высота над поверхностью. Ответ выразите через скорость воды у поверхности $U_{\text{пов}}$, $h$ и $H$. Выразите силу вязкого трения $F_{\text{вяз, 1}}$, действующую на участок слоя ширины $W$ и длины $L$ со стороны поверхности, через $U_{\text{пов}}$, $A$, $W$, $L$, $\eta$, где $A = WH$ — площадь сечения.

Рассмотрим слой, заключённый между высотами $h$ и $H$. В силу того, что течение носит установившийся характер, силы действующие на слой уравновешены.

Очевидно, что силы вязкого трения действуют только на нижнюю часть слоя, так как верхняя поверхность является свободной. Тогда уравнение баланса сил запишется следующим образом:
$$
\Delta P \cdot (H-h) W + \eta W L \frac{\partial u }{\partial h} = 0,
$$
где $\Delta P$ — постоянная разность давлений между торцами слоя. Отсюда несложно заключить, что зависимость $u(h)$ — квадратичная.

Подберём коэффициенты с помощью подстановки известных начальных условий:
$$
u(0) = 0 \quad u(H) = U_{\text{пов}} \\ F_{\text{тр}} (H) \propto \frac{\partial u }{\partial h} (H) = 0
$$
Решая систему из трёх линейных уравнений на параметры зависимости
$$
u(h) = C_1 h^2 + C_2 h + C_3,
$$
получим в итоге

Ответ: $$
u(h) = 2U_{пов} \left( \frac{h}{H}-\frac{h^2}{2H^2} \right)
$$

Выражение для силы вязкого трения можно получить, подставив зависимость $u(h)$ в явное выражение для $F_{\text{тр}}$:

Ответ: $$F_{вяз} = \frac{2\eta LW^2}{A}U_{пов}$$

D3  0.60 Выразите $U_{\text{пов}}$ через $U$.

$\textit{Напомним,}$ что $U$ — некоторая характерная величина, определяемая через выражение $Q = AU$.

Фактически, задача сводится к выражению суммарного потока жидкости через $U_{\text{пов}}$ с помощью элементарного интегрирования.
$$
Q = \int \limits_0^H u(h) W dh = 2U_{пов}WH \int \limits_0^1 \left( x-\frac{x^2}{2} \right) dx = \frac{2}{3} U_{пов}WH
$$
Из чего в итоге очевидно следует

Ответ: $$
U_{пов} = \frac{3}{2} U
$$

D4  0.20 Сформулируйте закон сохранения массы для данного кусочка. Выразите ответ через $U$, $A$, $ds$, $\rho$.

Ответ: $$
dm = \rho A ds = const
$$

D5  1.50 Запишите векторный II закон Ньютона в импульсной форме. Возьмите в качестве базиса ортонормированную тройку векторов $\left(\vec{t} (s), \vec{n} (s), \vec{b} (s) \right)$: $\vec{t}$ направлен вдоль струи; $\vec{n}$ — направлен перпендикулярно поверхности пробирки от нее; $\vec{b}$ — их ортогональное дополнение (в касательной плоскости перпендикулярно $\vec{t}$). При записи изменения импульса струи пренебрегайте распределением скоростей внутри нее и считайте ее равномерным по сечению потоком, двигающимся со скоростью $U(s)$, это приближение можно использовать в силу того, что учет неравномерности скорости в данном члене дает поправку, очень близкую к 1. Воспользуйтесь пунктом $D4$ для упрощения данного выражения.

Довольно очевидно можно получить, что импульс исследуемого кусочка
$$
dm \cdot \vec{U} = \rho A ds \cdot U \vec{t}
$$
Используя предыдущий пункт, $dm$ можно вынести из под знака дифференцирования по времени:
$$
\frac{d}{dt} \left[ \rho A ds \cdot U \vec{t} \right] = \rho A ds \cdot \frac{d}{dt} \left[ U \vec{t} \right]
$$
Осталось корректно записать силы, действующие на кусочек. Дополнительных преобразований может потребовать только сила вязкого трения: объединяя результаты $D2$ и $D3$ получим
$$
\vec F_{вяз} = - \frac{3C\eta W^2}{A} ds \cdot U \vec{t}.
$$
Комбинируя все полученные результаты, запишем закон в виде:

Ответ: $$
\rho Ads\cdot\frac{d}{dt}\left[U\vec{t}\right]=\rho A g \vec{z}ds-\Delta PW\vec{n}ds-\frac{3C\eta W^2}{A}U\vec{t}ds
$$

D6  0.20 Исключите время из уравнения, воспользовавшись тождественностью операторов
$$
\frac{d}{dt} \equiv U \cdot \frac{d}{ds}
$$

Ответ: $$
\rho AU\cdot\frac{d}{ds}\left[U\vec{t}\right]=\rho A g \vec{z}-\Delta PW\vec{n}-\frac{3C\eta W^2}{A}U\vec{t}
$$

D7  0.60 Найдите проекции $\cfrac{d}{ds} ( U \vec{t} ) $ на $\vec t$ и $\vec{b}$. Выразите ответы через функции от $U$, $\psi$ и их производные по $s$.

Дифференцируя произведение по правилу Лейбница, получим:
$$
\frac{d}{ds}\left[U\vec{t}\right] = \frac{dU}{ds} \vec t + U \frac{d \vec t}{ds}
$$
В силу того, что $\vec t$ — единичный, то
$$
\vec t \cdot \frac{d \vec t}{ds} = 0
$$
Проекцию $\vec b \cdot \frac{d \vec t}{ds}$ можно найти исходя из того, что угол $\psi$ определяет ориентацию векторов $\vec t, \vec b$ в касательной к трубке плоскости, что аналогично простому вращению пары векторов в $\mathbb{R}^2$ c использованием полярных координат:
$$
\vec b \cdot \frac{d \vec t}{ds} = - \frac{d \psi}{ds}
$$
Аналогично, последнее выражение можно получить дифференцированием тождества $\vec t \cdot \vec z = \cos \psi$, где $\vec z$ — постоянный единичный вектор, направленный вниз.
Итого:

Ответ: $$\left[\frac{d}{ds}\left(U\vec{t}\right)\right]\cdot\vec{t}=\frac{dU}{ds} \\
\qquad \: \left[\frac{d}{ds}\left(U\vec{t}\right)\right]\cdot\vec{b}=-\frac{d\psi}{ds}\cdot U$$

D8  0.60 Воспользовавшись результатами пунктов D4-D7, получите систему из двух диф.уравнений на $U$ и $\psi$ такого вида:
$$
\begin{cases}
\cfrac{dU}{ds} = - \alpha_1 U^2 + \cfrac{\alpha_2 \cos \psi}{U} \\
\cfrac{d \psi}{ds} = - \cfrac{\beta \sin \psi}{U^2},
\end{cases}
$$
где $\alpha_1$, $\alpha_2$ и $\beta$ — некоторые положительные константы. Выразите их через следующие постоянные: $\eta$, $U_0$, $D_j$, $\rho$, $g$, $C$.

Спроецировав уравнение $D6$ на $\vec t$ и $\vec b$, а затем воспользовавшись $D7$, получим следующие уравнения:
$$
\rho AU \cdot \frac{dU}{ds}=\rho A g \cos \psi -\frac{3C\eta W^2}{A}U \\
-\rho AU^2 \cdot \frac{d\psi}{ds}=\rho A g \sin \psi
$$
Разделив уравнения на множители перед производными для получения искомых коэффициентов остаётся вспомнить, что, строго говоря, $A$ является переменной величиной, зависящей от $U$. Считая жидкость несжимаемой, можно привести закон сохранения массы к закону сохранения объёмного потока, что будет иметь вид:
$$
Q = AU = \frac{\pi D_j^2}{4} U_0 = const
$$
Заменив с помощью этого выражения $A^2$ в последнем члене первого уравнения, получим:
$$
\frac{dU}{ds}=\frac{g \cos \psi}{U} -\frac{48C\eta}{\rho \pi^2D_j^2U_0^2} U^2 \\
\frac{d\psi}{ds}=-\frac{g \sin \psi}{U^2},
$$
из чего искомые коэффициенты

Ответ: $$
\alpha_1 = \frac{48C\eta}{\rho \pi^2D_j^2U_0^2} \\
\alpha_2= \beta = g
$$

D9  1.20 Рассчитайте значение $U$ в экспериментальных точках опираясь на параболическое приближение $\psi (s)$ (пункт D1). Погрешность $U$ можно не оценивать. Постройте график $U(s)$.

Из второго уравнения части $D8$ выразим $U$:
$$
U = \sqrt{-\frac{g \sin \psi}{d \psi /ds}}
$$
Производную $\frac{d \psi}{ds} $ рассчитаем, используя параболическое приближение:
$$
\frac{d \psi}{ds} = 2 a s + b
$$

D10  0.80 Оцените значение $C$.

Теперь воспользуемся первым уравнением пункта $D8$. Для оценки достаточно построить касательную к графику $U(s)$ в нуле, и, используя найденный коэффициент наклона, получить $C$:
$$
C = \frac{\rho \pi^2D_j^2}{48C\eta} \left[ \frac{g \cos \psi}{U_0} - \frac{dU}{ds}(0) \right]
$$

D11  1.20 Приближенно численно решите дифференциальные уравнения. Найдите вид «теоретической» зависимости $U_t(s)$, $\psi_t(s)$, $X_t(z)$ в не менее чем 15-ти точках.

D12  0.90 На графиках $z(X)$, $\psi (s)$, $U(s)$ отметьте точки теоретического моделирования $z_t(X)$, $\psi_t(s)$, $U_t(s)$.

E1  1.50 Экспериментально исследуйте зависимость $\alpha (U_0)$, постройте ее график. Сделайте не менее 13 измерений. Подробно опишите методику.

Методика измерения скорости воды:

Зафиксируем положение крана и иглы, что бы начальная скорость вытекания сохранялась постоянной. Определим эту скорость так же, как и в пункте A (однако теперь мы знаем $A$ и нам достаточно измерить объёмный расход).

Теперь будем изменять наклон струи, что бы изменить высоту, на которой струя становится горизонтальной. С помощью измерения этой высоты мы можем определить скорость на горизонтальном участке.

Методика определения угла поворота:

Поместим самую маленькую пробирку в отверстие экрана и прозрачной пластины с параллельными линиями. Повернём установку так, что бы изначальное направление струи совпадало с линиями на экране. С помощью транспортира определим необходимый угол, как показано на фотографии.

E2  0.50 При каком максимальном $U_{крит}$ наблюдается спираль?

$$U_{крит} = 0.60 \pm 0.05 \frac{м}{с}$$