Для получения зависимости $z(x)$ необходимо снять координаты струи. Это можно сделать, поместив миллиметровку параллельно плоскости движения воды и ручкой проставить точки.
Получим линеаризацию для зависимости $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)$.
С помощью мерного цилиндра и секундомера можно измерить объёмный расход $Q$ при выбранной скорости.
$$A = \frac{Q}{U_0} = 1.96 \; мм^2$$
$$D = 2\sqrt{\frac{A}{\pi}} = 1.58 \;мм$$
Сделаем установку, которая будет создавать необходимую спираль.
Поместим внутрь пробирки миллиметровку и снимем зависимость по этой миллиметровке.
При снятии зависимости будем учитывать, что внешний и внутренний диаметры пробирки различаются, поэтому скорректируем значения $x$:
$$x_{ист} = \frac{D_{внешн}}{D_{внутр}}x$$
$\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}$
При измерении $\lambda$ на различных пробирках важно сохранять одинаковые условия измерения, в них входят $\psi_0$, скорость струи при попадании на пробирку, степень прижатости струи к пробирки, а так же метод измерения $\lambda$.
При нанесении точек, соответствующих различным пробиркам, на изначальный график $z(x)$ можно убедится, что с точностью до погрешности они попадают на зависимость. Это говорит о том, что радиус пробирки не влияет на зависимость $z(x)$ при одинаковых начальных условиях.
Продифференциируем зависимость $\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.$$
Посчитав все необходимые суммы решим её.
Рассмотрим слой, заключённый между высотами $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)$ в явное выражение для $F_{\text{тр}}$:
Фактически, задача сводится к выражению суммарного потока жидкости через $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
$$
Из чего в итоге очевидно следует
Довольно очевидно можно получить, что импульс исследуемого кусочка
$$
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}.
$$
Комбинируя все полученные результаты, запишем закон в виде:
Дифференцируя произведение по правилу Лейбница, получим:
$$
\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$ — постоянный единичный вектор, направленный вниз.
Итого:
Спроецировав уравнение $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},
$$
из чего искомые коэффициенты
Из второго уравнения части $D8$ выразим $U$:
$$
U = \sqrt{-\frac{g \sin \psi}{d \psi /ds}}
$$
Производную $\frac{d \psi}{ds} $ рассчитаем, используя параболическое приближение:
$$
\frac{d \psi}{ds} = 2 a s + b
$$
Теперь воспользуемся первым уравнением пункта $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]
$$
Методика измерения скорости воды:
Зафиксируем положение крана и иглы, что бы начальная скорость вытекания сохранялась постоянной. Определим эту скорость так же, как и в пункте A (однако теперь мы знаем $A$ и нам достаточно измерить объёмный расход).
Теперь будем изменять наклон струи, что бы изменить высоту, на которой струя становится горизонтальной. С помощью измерения этой высоты мы можем определить скорость на горизонтальном участке.
Методика определения угла поворота:
Поместим самую маленькую пробирку в отверстие экрана и прозрачной пластины с параллельными линиями. Повернём установку так, что бы изначальное направление струи совпадало с линиями на экране. С помощью транспортира определим необходимый угол, как показано на фотографии.
$$U_{крит} = 0.60 \pm 0.05 \frac{м}{с}$$