Скорость течения жидкости постоянна, поэтому сумма силы давления и силы вязкого трения, действующей на слой жидкости толщины $dx$ должна быть равна нулю, откуда получаем уравнение
$$
\eta \frac{d^2 v}{dx^2}= -\frac{\Delta P}{L},
$$
его решение, удовлетворяющее правильным граничным условиям
$$
v(x) = \frac{\Delta P}{2 \eta L} \left( \frac{d^2}{4} -x^2\right).
$$
Интегрируя, находим среднюю скорость
$$
\langle v \rangle =\frac{1}{d} \int_{-d/2}^{d/2}dx\, v(x) = \frac{\Delta P d^2}{12 \eta L}
$$
и поток
$$
q = dw\langle v \rangle = \frac{d^3 w}{12 \eta L} \Delta P.
$$
Приравнивая нулю силу, действующую на цилиндр жидкости радиуса $\rho$, находим уравнение
$$
2 \pi \rho \frac{dv}{d\rho}\eta L = -\pi \rho ^2 \Delta P, \quad
\frac{dv}{d\rho} =- \rho \frac{\Delta P}{2 \eta L}.
$$
Его решение
$$
v(\rho) = \frac{\Delta P}{4 \eta L} (r^2 - \rho^2),
$$
откуда поток равен
$$
q = \int_0^r 2\pi \rho d\rho\, v(\rho) = \frac{\pi r^4 \Delta P}{8 \eta L}.
$$
Используя результаты предыдущих пунктов, получим
Работа сил вязкого трения равна по модулю и противоположна по знаку работе сил давления (жидкость движется равномерно, поэтому полная работа равна нулю).
Сила давления, действующая на узкий цилиндр жидкости площади $dS$ (ось цилиндра параллельна направлению движения жидкости) равна $F =\Delta P dS $. Соответствующий вклад в мощность $dW = \Delta P dS v$, а полная мощность сил трения
$$
W = -\int \Delta P v dS = - \Delta P q.
$$
$$W =-q\Delta P =-\Delta P^2/R= - \frac{\pi r^4}{8 \eta L}\Delta P^2$$
Альтернативно можно использовать факт, что в бесконечно малом объеме мощность сил вязкого трения $dW = - dV \eta (dv/d\rho)^2$.
Интегрируя по объему, найдем
$$
W = \int dV \eta (dv/d\rho)^2= - \eta L \int_0^r d\rho\, 2\pi \rho (dv/d\rho)^2 = - \frac{\pi r^4}{8 \eta L}\Delta P^2.
$$
Напишем выражения для потоков жидкостей через их средние скорости и ширины потоков:
$$
q_1 = \langle v_1 \rangle d w_1, \quad q_2 = \langle v_2 \rangle d w_2= \langle v_2 \rangle d (w-w_1).
$$
Решая полученные уравнения, находим
$$
\frac{q_2}{q_1} = \frac{\langle v_2 \rangle}{\langle v_1 \rangle}\left( \frac{w}{w_1}-1\right), \quad \frac{w_1}{w} = \frac{1}{1+ q_2 \langle v_1 \rangle/q_1\langle v_2 \rangle}.
$$
Из формул для средних скоростей
$$
\frac{\langle v_1 \rangle}{\langle v_2 \rangle} = \frac{\eta_2}{\eta_1}
$$
$$
\frac{q_1}{q_2} = \frac{\rho_1 g h_1/R_1}{\rho g h_2/R_2} = \frac{\rho_1 h_1}{\rho_2 h_2} \frac{R_2}{R_1} = \frac{\rho_1 h_1 \eta_2}{\rho_2 h_2 \eta_1 }
$$
Без учета гравитации
$$
\langle v \rangle =\frac{ d^2}{12 \eta L} \Delta P.
$$
За счет гравитации возникает дополнительный вклад в разность давлений
$$
\Delta P \to \Delta P + \rho g L.
$$
Отсюда скорость
Запишем выражения для средних скоростей двух жидкостей
$$
\langle v_1 \rangle =\frac{ d^2}{12 \eta_1 } \left(\frac{\Delta P}{L}+ \rho_1 g\right),
\quad
\langle v_2 \rangle =\frac{ d^2}{12 \eta_2 } \left(\frac{\Delta P}{L}+ \rho_2 g\right).
$$
Умножим эти уравнения на соответствующие вязкости и вычтем одно из другого
$$
\eta_1 \langle v_1 \rangle - \eta_2 \langle v_2 \rangle = \frac{gd^2}{12 } (\rho_1 - \rho_2).
$$
Разделив на $\eta_1 \langle v_2 \rangle$ получим
$$
\frac{\langle v_1 \rangle}{\langle v_2 \rangle} - \frac{\eta_2}{\eta_1} = \frac{g d^2}{12 \eta_1 \langle v_2 \rangle} (\rho_1-\rho_2)
$$
откуда
Используем соотношение из B1
$$
\frac{w}{w_1}=1+ q_2 \langle v_1 \rangle/q_1\langle v_2 \rangle.$$
Подставим в него найденное ранее отношение скоростей
$$
\frac{w}{w_1}=1+ \frac{q_2 \eta_2}{q_1 \eta_1} \left( 1+\frac{g d^2}{12 \eta_2 \langle v_2 \rangle} (\rho_1-\rho_2) \right).
$$
Выразим среднюю скорость $\langle v_2 \rangle$ через соответствующий поток:
$$
\langle v_2 \rangle = \frac{q_2}{d w_2} = \frac{q_2}{d(w-w_1)}.
$$
Тогда
$$
\frac{w}{w_1}=1+ \frac{q_2 \eta_2}{q_1 \eta_1} \left( 1+\frac{g d^3}{12 \eta_2 q_2} (\rho_1-\rho_2) (w-w_1)\right).
$$
Преобразуем это выражение:
$$
\frac{w}{w_1}=1+ \frac{q_2 \eta_2}{q_1 \eta_1} \left( 1+\frac{g d^3 w}{12 \eta_2 q_2} (\rho_1-\rho_2) \frac{w-w_1}{w}\right) =
$$
$$
1 + \alpha \left( 1+ \frac{q_G (\rho_1 - \rho_2)}{q_2 \rho_2}\left(1- \frac{w_1}{w}\right)\right) = 1+ \alpha \left( 1+ \beta \left(1- \frac{w_1}{w} \right)\right).
$$
Получили квадратное уравнение
$$
\alpha \beta\frac{w_1^2}{w^2} - (1 + \alpha + \alpha \beta)\frac{w_1}{w} +1 = 0.
$$
Его решение (нужный нам корень должен быть меньше 1, поэтому берем знак -)
Подставляем приведенные выше числа в формулы для сопротивления и $q_G$
$$
R = \frac{8 \eta L}{\pi r^4}, \quad q_G = \frac{\rho_2 g d^3 w}{12 \eta_2}
$$
Подставим числа в формулы, полученные в C4. Подставим числа в формулы из B4. Получим значения для $h= 0.25~м, \,0.5~м, \,1.0~м$ соответственно:
Подставим числа в формулы из B4. Получим значения для $h= 0.25~м, \,0.5~м, \,1.0~м$ соответственно:
На рисунке приведены графики для случаев a), b), c) соответственно. Обратите внимание, что все три зависимости монотонно убывают, а также на порядок расположения кривых (когда плотность центральной жидкости больше, поток становится уже после учета гравитации, а когда меньше — поток наоборот расширяется).
Скорость уменьшения уровня воды можно оценить из
$$
\pi r_r^2 \frac{dh_2}{dt} = -q_2 = -\frac{\rho_2 g h_2}{R}.
$$
Считая эту скорость приближенно постоянной, находим