§ 1

При некоторых движениях жидкостей, которые мы рассмотрим здесь, полезно знать потенциал эллипсоида, наполненного массой постоянной плотности. Мы не будем выводить выражения этого потенциала. Но приведем его выражение непосредственно и, следуя по весьма простому пути, указанному Дирихле*, докажем его правильность. Пусть уравнение поверхности эллипсоида
$$\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}=1,$$

его плотность равна единице, и его потенциал в точке с координатами $x,y,z$ пусть будет $\mathrm{\Omega }$. Мы сопоставим известные свойства, которыми должно обладать $\mathrm{\Omega }$. Вследствие уравнения (11) шестнадцатой лекции, если точка $(x,y,z)$ лежит внутри эллипсоида, должно быть
$$\mathrm{\Delta }\mathrm{\Omega }=-4\pi ;$$

если она лежит вне, то
$$\mathrm{\Delta }\mathrm{\Omega }=0.$$

Далее $\mathrm{\Omega },\frac{\partial \mathrm{\Omega }}{\partial x},\frac{\partial \mathrm{\Omega }}{\partial y},\frac{\partial \mathrm{\Omega }}{\partial z}$ должны быть однозначны и непрерывны во всем пространстве, и $\mathrm{\Omega }$ должно в бесконечности обращаться в нуль.

Эти свойства однозначно определяют $\mathrm{\Omega }$, как это показывает следующее рассуждение. Положим, что существуют две функции $x,y,z$, которые обладают этими свойствами, и пусть $U$ — их разность. Тогда $U$ имеет такие свойства во всем пространстве: $\mathrm{\Delta }U=0;U,\frac{\partial U}{\partial x},\frac{\partial U}{\partial y},\frac{\partial U}{\partial z}$ — однозначны и непрерывны; в бесконечности $U=0$. Но тогда, согласно разъяснению, сделанному в $\mathrm{\S }7$ шестнадцатой лекции, всюду $U=0$. Пусть будет, как в предыдущей лекции, $u$ — наибольший корень кубического уравнения
$$\frac{x^2}{a^2+u}+\frac{y^2}{b^2+u}+\frac{z^2}{c^2+u}=1$$

если точка ( $x,y,z$ ) лежит вне эллипсоида (1), то $u$ положительно и представляет единственный положительный корень этого уравнения; на поверх-
* Crelle’s Journal, Bd. 32, S. 80.

ности эллипсоида $u=0$. Тогда будем иметь для внешней точки
$$\mathrm{\Omega }=\pi abc{\int }_u^{\mathrm{\infty }}d\lambda \frac{1-\frac{x^2}{a^2+\lambda }-\frac{y^2}{b^2+\lambda }-\frac{z^2}{c^2+\lambda }}{\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}},$$

для внутренней
$$\mathrm{\Omega }=\pi abc{\int }_0^{\mathrm{\infty }}d\lambda \frac{1-\frac{x^2}{a^2+\lambda }-\frac{y^2}{b^2+\lambda }-\frac{z^2}{c^2+\lambda }}{\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}}.$$

Чтобы доказать правильность этого утверждения, нам нужно только показать, что функция, данная уравнениями (3) и (4), обладает свойствами, которые определяют $\mathrm{\Omega }$. Для этого продифференцируем уравнения (3) и (4) сперва по $x$. Для внешней точки мы получим
$$\frac{\partial \mathrm{\Omega }}{\partial x}=-2\pi abcx{\int }_u^{\mathrm{\infty }}\frac{d\lambda }{(a^2+\lambda )\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}},$$

так как член, происходящий от переменности нижнего предела, обращается в нуль вследствие уравнения (2); для внутренней точки
$$\frac{\partial \mathrm{\Omega }}{\partial x}=-2\pi abcx{\int }_0^{\mathrm{\infty }}\frac{d\lambda }{(a^2+\lambda )\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}}.$$

Пусть аналогичные уравнения составлены для $\frac{\partial \mathrm{\Omega }}{\partial y}$ и $\frac{\partial Q}{\partial z}$. При этом заметим, что $\mathrm{\Omega },\frac{\partial \mathrm{\Omega }}{\partial x},\frac{\partial \mathrm{\Omega }}{\partial y},\frac{\partial \mathrm{\Omega }}{\partial z}$ не претерпевают разрыва на поверхности эллипсоида, где $u=0$, и, следовательно, однозначны и непрерывны всюду. Продифференцируем уравнения (5) и (6) еще раз по $x$; тогда для внешней точки получим
$$\frac{{\partial }^2\mathrm{\Omega }}{\partial x^2}=2\pi abc\{\frac{\frac{x}{a^2+u}\frac{\partial u}{\partial x}}{\sqrt{(a^2+u)(b^2+u)(c^2+u)}}-{\int }_u^{\mathrm{\infty }}\frac{d\lambda }{(a^2+\lambda )\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}}\},$$

а для внутренней
$$\frac{{\partial }^2\mathrm{\Omega }}{\partial x^2}=-2\pi abc{\int }_0^{\mathrm{\infty }}\frac{d\lambda }{(a^2+\lambda )\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}}.$$

Складывая с этим уравнением соответствующие уравнения для $\frac{{\partial }^2\mathrm{\Omega }}{\partial y^2}$ и $\frac{{\partial }^2\mathrm{\Omega }}{\partial z^2}$ и пользуясь тем, что
$$\begin{array}{c}\int (\frac{1}{a^2+\lambda }+\frac{1}{b^2+\lambda }+\frac{1}{c^2+\lambda })\frac{d\lambda }{\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}}\cdots \\ -\frac{2}{\sqrt{(a^2+\lambda )(b^2+\lambda )(c^2+\lambda )}},\end{array}$$

для внутренней точки мы получим
$$\mathrm{\Delta }\mathrm{\Omega }=-4\pi ;$$

для внешней, если еще добавим, что по уравнению (17) предыдущей лекции
$$\frac{x}{a^2+u}-\frac{\partial u}{\partial x}+\frac{y}{b^2+u}\frac{\partial u}{\partial y}+\frac{z}{c^2+u}-\frac{\partial u}{\partial z}=2,$$

будем иметь
$$\mathrm{\Delta }\mathrm{\Omega }=0.$$

Наконец, чтобы убедиться в том, что определяемая уравнением (3) функция $\mathrm{\Omega }$ в бесконечности обращается в нуль, надо только заметить, что вследствие уравнения (2) $u$ бесконечно велико на бесконечности.

Добавим к этому еще определение потенциала цилиндра, наполненного массой с плотностью, равной единице, уравнение боковой поверхности которого есть
$$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1\text{. }$$

и уравнения оснований суть
$$z=-\gamma \text{ и }z=-\gamma ,$$

где $\gamma$ должно быть рассматриваемо как величина бесконечно большая по сравнению с $a,b$ и с координатами точек, к которым потенциал относится. Обозначим последний через $\mathrm{\Omega }$; тогда по уравнению (30) шестнадцатой лекции будем иметь
$$\mathrm{\Omega }=2\pi ab\lg 2\gamma +U,U=-2\int df\lg \rho ,$$

где $df$ — элемент площади, ограниченной эллипсом $\frac{x^2}{a^2}+\frac{y^2}{b^2}=1,z:=0$, $\rho$ — расстояние этого элемента от точки $(x,y,0)$, если $(x,y,z)$ есть точка, к которой относятся $U$ и $\mathrm{\Omega }$. По проведенному в шестнадцатой лекции исследованию $U$ определено с точностью до произвольной постоянной, если оно со своими первыми производными однозначно и непрерывно на всей плоскости $xOy$ и удовлетворяет уравненио
$$\frac{{\partial }^{2}U}{\partial x^2}+\frac{{\partial }^{2}U}{\partial y^2}\cdots -4\pi \text{ или }\frac{{\partial }^{2}U}{\partial x^2}-\frac{{\partial }^{2}U}{\partial y^2}=0,$$

что зависит от того, лежит ли точка $(x,y)$ внутри или вне указанного эллипса, а $\frac{\partial U}{\partial x}$ и $\frac{\partial U}{\partial y}$ на бесконечности обращаются в нуль. $U$ будет обладать этими свойствами, если положим для внутренней точки
$$U=\pi ab({\int }_0^{u^{\prime }}\frac{1-\frac{x^2}{a^2+\lambda }-\frac{y^2}{b^2+\lambda }}{\sqrt{(a^2+\lambda )(b^2-\lambda )}}d\lambda +1-\lg u^{\prime })\text{, }$$

а для внешней точки
$$U=\pi ab({\int }_u^{u^{\prime }}\frac{1-\frac{x^2}{a^2+\lambda -\lambda -y^2}-\lambda }{\sqrt{(a^2+\lambda )(b^2+\overline{\lambda })}}d\lambda +1-\lg u^{\prime }),$$

при этом под $u^{\prime }$ понимаем бесконечно большое постоянное, а под $u$ — положительный корень уравнения
$$\frac{x^2}{a^2+u}+\frac{y^2}{b^2+u}=1.$$

Это вытекает из рассуждений, подобных рассуждениям первой части этого параграфа.

Для доказательства того, что эти выражения для $U$ не вводят иного постоянного, кроме содержащегося для $U$ в формуле (7), надо только показать, что если мы положим $x^2+y^2=p^2$, где ${\rho }^2$ бесконечно велико, но бесконечно мало по сравнению с $u^{\prime }$, то получим из (9)
$$U=-2\pi ab\lg \rho .$$

Это следует из того, что при сделанном предположении $t={\rho }^2$, и
$${\int }_{{\rho }^2}^{u^{\prime }}\frac{1}{\lambda }(1-\frac{{\rho }^2}{\lambda })d\lambda =\lg u^{\prime }-2\lg \rho -1.$$

Нетрудно произвести интегрирование, указанное в выражениях (8) и (9). Тогда для внутренней точки найдем
$$U=\pi ab(2\lg \frac{2}{a+b}+1)-2\pi \frac{b{x}^2+a{y}^2}{a+b}.$$

Если эллипс превращается в круг радиуса $a$, то это уравнение дает
$$U=\pi a^2(1-2\lg a)-\pi {\rho }^2,$$

где по-прежнему $x^2+y^2={\rho }^2$; для внешней точки при этом будем иметь
$$U=-2\pi a^2\lg \rho .$$
§ 2
Положим ${}^{32}$
$$\phi =M\frac{\partial \mathrm{\Omega }}{\partial z}-z,$$

где $M$ означает постоянное, а $\mathrm{\Omega }$ имеет данное уравнением (3) значение; тогда в пространстве, внешнем относительно эллипсоида (1), ч удовлетворяет уравнению $\mathrm{\Delta }\phi =0$.

Поэтому уравнение (11) представляет возможное движение жидкости в этом пространстве. На бесконечности имеем
$$\frac{\partial \phi }{\partial x}=0,\frac{\partial \phi }{\partial y}=0,\frac{\partial \phi }{\partial z}=\cdots 1$$

следовательно, в бесконечности жидкость течет со скоростью, равной единице, в направлении, противоположном оси $z$. Мы покажем, что постоянное $M$ можно определить так, чтобы на поверхности эллипсоида (1) было
$$\frac{\partial \phi }{\partial n}=0\text{. }$$

Если постоянное определено именно так, то уравнение (11) пригодно для случая, когда жидкость движется определенным образом в бесконечность и в ней покоится эллипсоид (1). Обозначим через $n_a$ нормаль к элементу поверхности эллипсоида, направленную внутрь жидкости, нормаль же, направленную внутрь эллипсоида, обозначим через $n_i$; тогда получим условие, которое должно быть выполнено при надлежащем выборе $M$ :
$$M\frac{{\partial }^2\mathrm{\Omega }}{\partial n_a\partial z}=\cos (n_a,z),$$

или, так как по уравнению (10) шестнадцатой лекции
$$\frac{{\partial }^2\mathrm{\Omega }}{\partial n_i\partial z}+\frac{{\partial }^2\mathrm{\Omega }}{\partial n_a\partial z}=4\pi \cos (n_a,z)$$

то
$$M\frac{{\partial }^2\mathrm{\Omega }}{\partial n_i\partial z}=(4\pi M-1)\cos \left(n_{a}z\right).$$

Мы ввели вместо производной по $n_a$ производную по $n_i$, так как определяемое из уравнения (4) значение $\mathrm{\Omega }$ для внутренней точки проще, чем определяемое уравнением (3) для внешней точки. Действительно, для внутренней точки по уравнению (4) имеем выражение
$$\mathrm{\Omega }=\text{ const }-\pi (A{x}^2+B{y}^2+C{z}^2),$$

которым пользовались уже в двенадцатой лекции и где $A,B,C$ суть постоянные, значения которых были определены уравнением (4) той же лекции.
Отсюда следует, что
$$\frac{{\partial }^2\mathrm{\Omega }}{\partial n_i\partial z}=-2\pi C\frac{\partial z}{\partial n_i}=2\pi C\cos \left(n_{a}z\right).$$

Поэтому уравнение (12) будет выполнено, если мы положим
$$M=\frac{1}{2\pi (2-C)}.$$

Если для установившегося движения жидкости с потенциалом скоростей этот потенциал $\phi$ известен, то определение линий тока потребует еще интегрирования дифференциальных уравнений
$$dx:dy:dz=\frac{\partial \phi }{\partial x}:\frac{\partial {\phi }_{-}}{\partial y}:\frac{\partial \phi }{\partial z},$$

интегралы которых, содержащие две произвольные постоянные, дают уравнения линий тока.

Если $\phi$ есть функция двух аргументов: $z$ и $\sqrt{x^2+y^2}$ [условие, которое в рассматриваемом здесь случае будет выполнено, если эллипсоид (1) есть эллипсоид вращения и $a=b$ ], то интегрирование уравнений (14) приводится к квадратурам. В самом деле, положим
$$\sqrt{x^2+y^2}=\rho$$

и будем рассматривать $\phi$ как функцию $z$ и $\rho$; тогда получим
$$\frac{\partial \phi }{\partial x}=\frac{\partial \phi }{\partial \rho }\cdot \frac{x}{\rho },\frac{\partial \phi }{\partial y}=\frac{\partial \phi }{\partial \rho }\frac{y}{\rho }.$$

Поэтому одно из уравнений (14) будет
$$xdy-ydx=0,$$

откуда следует
$$\frac{x}{y}=\text{ const, }$$
т. е. каждая линия тока лежит в плоскости, проходящей через ось (z). Из (15) в связи с уравнением
$$xdx+ydy=\rho d\rho$$

получается далее
$$dx=\frac{xd\rho }{\rho },$$

вследствие чего второе из уравнений (14) примет вид
$$\frac{\partial \phi }{\partial z}d\rho -\frac{\partial \phi }{\partial \rho }dz=0.$$

Но из уравнения $\mathrm{\Delta }\phi =0$ можно убедиться, что интегрирующий множитель этого уравнения есть $\rho$, т. е. левая часть уравнения
$$\rho \frac{\partial \phi }{\partial z}d\rho -\rho \frac{\partial \phi }{\partial \rho }dz=0$$

есть полный дифференциал. Действительно, мы имеем
$$\begin{array}{l}\frac{{\partial }^2\phi }{\partial x^2}=\frac{{\partial }^2\phi }{\partial {\rho }^2}\frac{x^2}{{\rho }^2}+\frac{\partial \phi }{\partial \rho }\frac{y^2}{{\rho }^3},\\ \frac{{\partial }^2\phi }{\partial y^2}=\frac{{\partial }^2\phi }{\partial {\rho }^2}\frac{y^2}{{\rho }^2}+\frac{\partial \phi }{\partial \rho }\frac{x^2}{{\rho }^3},\end{array}$$

откуда следует
$$\mathrm{\Delta }\phi =\frac{{\partial }^2\phi }{\partial z^2}+\frac{{\partial }^2\phi }{\partial {\rho }^2}+\frac{1}{\rho }\cdot \frac{\partial \phi }{\partial \rho }=0,$$

или
$$\frac{\partial }{\partial z}\left(\rho \frac{\partial \phi }{\partial z}\right)=-\frac{\partial }{\partial \rho }\left(\rho \frac{\partial \phi }{\partial \rho }\right),$$

чем и доказано высказанное утверждение. Таким образом, существует функция $U$ переменных $z$ и $\rho$, имеющая то свойство, что
$$\frac{\partial U}{\partial z}=\rho \frac{\partial \phi }{\partial \rho },\frac{\partial U}{\partial \rho }=-\rho \frac{\partial \phi }{\partial z};$$

если она найдена, то мы получим
$$U=\text{ const }$$

как интеграл уравнения (16) и как уравнение линий тока.
Пусть будет известна функция $V$ переменных $z$ и $\rho$, для которой
$$\frac{\partial V}{\partial z}=\phi \text{ и }\mathrm{\Delta }V=\frac{{\partial }^{2}V}{\partial z^2}+\frac{{\partial }^{2}V}{\partial {\rho }^2}+\frac{1}{\rho }\frac{\partial V}{\partial \rho }=0;$$

гогда можно будет положить
$$U=\rho \frac{\partial V}{\partial \rho },$$

так как при этом будут удовлетворены оба уравнения (17). Следовательно, уравнение линий тока будет
$$\rho \frac{\partial V}{\partial \rho }=\text{ const. }$$

В случае движения жидкости, представляемого уравнениями (11) в предположении, что $a=b$, условиям (18) можно удовлетворить, положив
$$V=-M\mathrm{\Omega }-\frac{z^2}{2}+\frac{{\rho }^2}{4};$$

тогда уравнение линий тока будет
$$\rho (M\frac{\partial \mathrm{\Omega }}{\partial \rho }+\frac{\rho }{2})=\text{ const. }$$

Если эллипсоид вращения превращается в шар радиуса $R$, то для внешней точки будем иметь
$$\mathrm{\Omega }=\frac{4\pi }{3}{R}^3\frac{1}{r}\text{. }$$

Из этого, пользуясь предложением, что $\frac{\partial \mathrm{\Omega }}{\partial z}$ не претерпевает разрыва на поверхности шара, можно вывести, что $C=\frac{2}{3}$, и, следовательно, по (13), $M=-\frac{3}{8\pi }$. Отсюда следует, что
$$\phi =\frac{R^3}{2}\frac{\partial \frac{1}{r}}{\partial z}-z$$

и что уравнение линий тока есть
$${\rho }^2(1-\frac{R^3}{r^3})=\text{ const. }$$

Положим, что величина, обозначенная в этом уравнении через const, равна нулю; тогда мы удовлетворим этому уравнению, полагая $\rho =0$ или $r=R$. В этом случае линия тока состоит из отрезка оси $z$, лежащего вне шара, и из полукруга, по которому поверхность шара пересекается с полуплоскостью, ограниченной осью $z$ и проходящей через ось $\rho$. В точках, в которых части оси $z$ встречаются с этим кругом, т. е. в точках $\rho 0$ и $z=\pm R$, скорость равна нулю, как это следует из уравнений
$$\frac{\partial \phi }{\partial z}=\frac{R^3}{3}(\frac{3z^2}{r^5}-\frac{1}{r^3})-1$$

и
$$\frac{\partial p}{\partial \rho }\cdot \frac{R^3}{2}\cdot \frac{3zp}{r^5}.$$

которые, вообще, определяют компоненты скорости по осям $z$ и $\rho$. Эти точки обладают тем свойством, что никакая частица жидкости, лежащая в данный момент на конечном от них расстоянии, никогда не может их достигнуть. Вычисление показывает, что необходимое для этого время будет бесконечно велико. Мы покажем это для частиц, лежащих на названном полукруге. Для такой частицы $r=R$; следовательно, $\frac{\partial \phi }{\partial z}=\frac{3}{2}(\frac{z^2}{R^2}-1)$,
и потому $\frac{dz}{dt}:=\frac{3}{2}(\frac{z^2}{R^2}-1)$. Обозначим значение $t$, для которого $z:=0$, через $t_0$; тогда интегрированием получим
$$\frac{3}{4}(t-t_0)-\lg \frac{R-z}{R+z},$$

откуда действительно получаем
$$t=\mathrm{\infty }\text{ для }z\cdots -R.$$
§ 3
Выведенным уравнениям мы можем придать еще другое, отличное от прежнего, значение. Согласно разъяснению, данному в § 4 четвертой лекции, для системы координат, оси которой двикутся поступательно с постоянной скоростью, пригодны те же дифференциальные уравнения движения, как и для неподвижной. Представим себе, что оси $x,y,z$ неизменно связаны с эллипсоидом и движутся с ним в некотором направлении с постоянной скоростью; тогда выведенные в предыдущем параграфе формулы будут пригодны для движения жидкости относительно эллипсоида. Допустим, что движение происходит в направлении оси $z$ со скоростью, равной единице, так что при этом абсолютная скорость частиц жидкости в бесконечности равна нулю.

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

Воспользуемся обозначениями пятой лекции и введем две прямоугольные системы координат, одна из которых ( $\xi ,\eta ,\zeta$ ) неподвижна в пространстве, а другая ( $x,y,z$ ) неподвижна в теле; напишем уравнения между координатами одной и той же точки в обеих системах:
$$\begin{array}{l}\xi =\alpha +{\alpha }_{1}x+{\alpha }_{2}y+{\alpha }_{3}z,\\ \eta =\beta +{\beta }_{1}x+{\beta }_{2}y+{\beta }_{3}z,\\ \zeta =\gamma +{\gamma }_{1}x+{\gamma }_{2}y+{\gamma }_{3}z.\end{array}$$

Обозначим далее через $u$, $v$, w компоненты скорости начала координат системы $x,y,z$ по осям $x,y,z$ и через $p,q,r$ — компоненты угловой скорости тела относительно тех же осей. Тогда выражения (1) шестой лекции, именно выражения
$$\begin{array}{c}u+zq-yr,\\ v-xr-zp,\\ w+yp-xq,\end{array}$$

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

Из последнего можно заключить, что если мы отнесем $x,y,z$ к этой же жидкой частице, то получим
$$\begin{array}{l}\frac{dx}{dt}=\frac{\partial \phi }{\partial x}-u-zq+yr,\\ \frac{dy}{dt}=\frac{\partial \phi }{\partial y}-v-xr+zp,\\ \frac{dz}{dt}=\frac{\partial \phi }{\partial z}-w-yp+xq.\end{array}$$

Интегрирование этих уравнений дает относительное движение всех жидких частиц по отношению к телу, если $u,v,w,p,q,r$ — известные функции $t$ и $\phi$. Если бы мы пожелали найти их абсолютное значение, то надо было бы вычислить $\xi ,\eta ,\zeta$ из уравнений (20), пользуясь определенными значениями $x,y,z$.

Посмотрим теперь, как определить $\phi$. Обозначим через $n$ направленную внутрь жидкости нормаль к элементу поверхности тела; тогда выражения (21), по умножении их на $\cos (nx),\cos (ny),\cos (nz)$ и после сложения, дадут компоненту скорости рассматриваемого элемента на направление $n$, но эта компонента должна быть равна $\frac{\partial \phi }{\partial n}$, т. е.
$$\frac{\partial \phi }{\partial n}=(u+zq-yr)\cos (nx)+(v+xr-zp)\cos (ny)+(w+yp-xq)\cos (nz).$$

K этому добавляются следующие условия: в бесконечности $\frac{\partial \phi }{\partial x},\frac{\partial \phi }{\partial y}$, $\frac{\partial \phi }{\partial z}$ обращаются в нуль; в объеме, наполненном жидкостью, $\mathrm{\Delta }\phi =0$ и $\frac{\partial \phi }{\partial x},\frac{\partial \phi }{\partial y},\frac{\partial \phi }{\partial z}$ однозначны и непрерывны. Что $\phi$ однозначно, мы уже допустили; не ограничивая общности, мы можем предположить, что $\phi$ также непрерывно, так как пока оно определено только через свои производные. Согласно разъяснению, сделанному в § 7 шестнадцатой лекции, $\phi$ определено здесь до добавочной постоянной; оно будет определено вполне, если мы примем, а это допустимо, что оно обращается в нуль в бесконечности. Все эти требования будут удовлетворены, если мы положим
$$\phi =u{\phi }_1+v{\phi }_2+w{\phi }_3+p{\phi }_4+q{\phi }_5+r{\phi }_6$$

и определим шесть функций ${\phi }_1,{\phi }_2,\dots$ так, чтобы каждая из них удовлетворяла уравнению $\mathrm{\Delta }\phi =0$, заданным для $\phi$ условиям непрерывности и обращалась в нуль в бесконечности, а на поверхности движущегося тела удовлетворяла уравнениям
$$\begin{array}{rlrl}\frac{\partial {\phi }_1}{\partial n}& =\cos (nx),& \frac{\partial {\phi }_4}{\partial n}& =y\phi \cos (nz)-z\cos (ny),\\ \frac{\partial {\phi }_2}{\partial n}& =\cos (ny),& \frac{\partial {\phi }_5}{\partial n}& =z\cos (nx)-x\cos (nz),\\ \frac{\partial {\phi }_3}{\partial n}& =\cos (nz),& \frac{\partial {\phi }_6}{\partial n}& =x\cos (ny)-y\cos (nx).\end{array}$$

Этими условиями функции ${\phi }_1,{\phi }_2,\dots$ вполне определены; они не зависят от движения тела, но зависят исключительно от его формы. Если тело представляет эллипсоид, уравнение поверхности которого есть уравнение (1), то функции ${\phi }_1,{\phi }_2,\dots$ легко определить на основании результатов $\mathrm{\S }2$. Именно, пользуясь принятыми там обозначениями и принимая для нормали $n$ данное здесь направление, мы получим
$$\begin{array}{l}\frac{{\partial }^2\mathrm{\Omega }}{\partial n\partial x}=2\pi (2-A)\cos (nx),\\ \frac{{\partial }^{2}Q}{\partial n\partial y}=2\pi (2-B)\cos (ny),\\ \frac{{\partial }^2\mathrm{\Omega }}{\partial n\partial z}=2\pi (2-C)\cos (nz).\end{array}$$

Отсюда прежде всего будем иметь
$${\phi }_1=\frac{1}{2\pi (2-A)}\frac{\partial \mathrm{\Omega }}{\partial x},{\phi }_2=\frac{1}{2\pi (2-B)}\frac{\partial \mathrm{\Omega }}{\partial y},{\phi }_3=\frac{1}{2\pi (2-C)}\frac{\partial Q}{\partial z}.$$

Далее, мы можем заметить, что если $N$ обозначает определенное постоянное, то
$${\phi }_6=N(x\frac{\partial \mathrm{\Omega }}{\partial y}-y\frac{\partial \mathrm{\Omega }}{\partial x}),\frac{3}{2}$$

или, что то же самое,
$${\phi }_6=N-\frac{\partial \mathrm{\Omega }}{\partial \vartheta },$$

если положить
$$x=\rho \cos \vartheta ,y=\rho \sin \vartheta .$$

Все поставленные для ${\phi }_6$ условия будут удовлетворять определенным выражениям при любом значении $N$, за исключением последнего из уравнений (23), которое выполняется только при некотором определенном значении $N$. Именно, вследствие уравнений (24) и уравнений
$$\begin{array}{ll}\frac{\partial \mathrm{\Omega }}{\partial x}=-2\pi Ax,& \frac{\partial \mathrm{\Omega }}{\partial y}=-2\pi By,\\ \frac{\partial x}{\partial n}=\cos (nx),& \frac{\partial y}{-}=\cos (ny),\end{array}$$

это уравнение примет вид
$$x\cos (ny)-y\cos (nx)=2\pi N[x\cos (ny)(2-B+A)-y\cos (nx)(2-A+\phi )]\text{. }$$

Оно линейно и однородно относительно $x$ и $y$; на этом основании и так как
$$\cos (nx):\cos (ny)=\frac{x}{a^2}:\frac{y}{b^2},$$

можно положить, что $x$ и $y$ соответственно равны $a^2\cos (nx)$ и $b^2\cos (ny)$; тогда в обе части уравнения $\cos (nx)\cos (ny)$ войдет общим множителем. По сокращении на этот множитель мы получим
$$N=\frac{a^2-b^2}{2\pi [2(a^2-b^2)+(A-B)(a^2+b^2)]}.$$

Выражения, подобные (25), которым определено ${\phi }_6$, можно составить и для ${\phi }_4$ и ${\phi }_5$.

Если $a^2-b^2$ обращается в нуль, т. е. если эллипсоид переходит в эллипсоид вращения, ось которого совпадает с осью $z$, то отношение $(A-B):(a^2-b^2)$, а также и $N$ получают легко определяемые значения. Но множитель при $N$ в уравнении (25), а следовательно, и ${\phi }_6$ обращается в нуль.
Если эллипсоид превращается в шар радиуса $R$, то
$${\phi }_4=0,{\phi }_5=0,{\phi }_6=0,$$
т. е. вращение шара вокруг его центра не оказывает никакого влияния на движение жидкости. Далее, согласно сделанному в § 2 определению, будем иметь
$${\phi }_1=\frac{R^3}{2}\frac{\partial -\frac{1}{r}}{\partial x},{\phi }_2=\frac{R^3}{2}\frac{\partial \frac{1}{r}}{\partial y},{\phi }_3=\frac{R^3}{2}\frac{\partial \frac{1}{r}}{\partial z},$$

где
$$r=\sqrt{x^2+y^2+z^2}\text{. }$$

Поэтому для шара мы имеем
$$\phi =\frac{R^3}{2}(u\frac{\partial \frac{1}{r}}{\partial x}+v\frac{\partial \frac{1}{r}}{\partial y}+w\frac{\partial \frac{1}{r}}{\partial z}).$$

Это $\phi$ удовлетворяет на поверхности шара условию
$$\frac{\partial \phi }{\partial n}=u\cos (nx)+v\cos (ny)+w\cos (nz).$$

Здесь уместно упомянуть, что найденное для $\phi$ выражение согласуется с потенциалом молекулярного магнита, находящегося в центре шара, магнитная ось которого имеет направление движения, и магнитный момент которой равен скорости центра, умноженной на $\frac{R^3}{2}$, относительно полюса, лежащего в точке $(x,y,z)$ и содержащего единицу магнетизма. Это следует из того, что скорость в точке $(x,y,z)$ по величине и направлению равна силе, с которой этот молекулярный магнит действует на полюс.

Движение шара в жидкости впервые было исследовано Дирихле *, эллипсоида — Клебшем **.
$\mathrm{\$}4$
В предыдущих параграфах мы предполагали, что на конечном расстоянии жидкость ограничена только поверхностью движущегося твердого тела. В этом случае следует отнести потенциал скоростей $\phi$ к системе координат, неподвижно связанной с телом, потому что тогда он зависит исключительно от формы тела и его движения в рассматриваемый момент. Если кроме данного тела на конечном расстоянии находятся еще другие твердые тела, которые движутся или покоятся, то потенциал скоростей всегда будет зависеть от относительного положения всех тел. Тогда целесообразно отнести его прямо к неподвижной в пространстве системе координат. Мы будем теперь представлять себе, что в бесконечной жидкости на конечном расстоянии движутся два твердых тела и что система осей $x,y,z$ неподвижна в пространстве. Пусть $u,v,w$-компоненты скорости точки первого тела, $u^{\prime },v^{\prime }{w}^{\prime }$ — точки второго, $p,q,r$ — компоненты угловой ско-
* Monatsberichte der Berliner Academie, 1852, S. 12.
** Crelle’s Journal, Bd. 52, S. 103 u. Bd. 53, S. 287.

рости по осям, параллельным осям координат для первого тела, $p^{\prime },q^{\prime },r^{\prime }$ соответствующие величины для второго. Тогда, соответственно уравнению (22), можно положить
$$\begin{array}{l}\phi =u{\phi }_1+v{\phi }_2+w{\phi }_3+p{\phi }_4+q{\phi }_5+r{\phi }_6+\\ +u^{\prime }{\phi }_1^{\prime }+v^{\prime }{\phi }_2^{\prime }+w^{\prime }{\phi }_3^{\prime }+p^{\prime }{\phi }_4^{\prime }+q^{\prime }{\phi }_5^{\prime }+r^{\prime }{\phi }_6^{\prime },\end{array}$$

где ${\phi }_1,{\phi }_2,\dots ,{\phi }_1^{\prime },{\phi }_2^{\prime },\dots -$ функции $x,y,z$, которые зависят не только от движения обоих тел, но и от их мгновенного положения. Қаждая из них должна в наполненном жидкостью пространстве удовлетворять уравнению $\mathrm{\Delta }\phi =0$, а также быть однозначной и непрерывной вместе со своими первыми произвоннми, обращаться в нуль на бесконечности и на поверхности обоих тел удовлетворять некоторым двум уравнениям. А именно, пусть $a,b,c$ и $a^{\prime },b^{\prime },c^{\prime }$ — координаты двух точек, компоненты скоростей которых обозначены через $u,v$, w и $u^{\prime }$. $u^{\prime },{\omega }^{\prime };$ тогда на пэверхности первого тела должно быть
$$\begin{array}{ll}\frac{\partial {\phi }_1}{\partial n}=\cos (nx),& \frac{\partial {\phi }_1^{\prime }}{\partial n}=0,\\ \frac{\partial {\phi }_2}{\partial n}=\cos (ny),& \frac{\partial {\phi }_2^{\prime }}{\partial n}=0,\\ \frac{\partial {\phi }_3}{\partial n}=\cos (nz),& \frac{\partial {\phi }_3^{\prime }}{\partial n}=0,\\ \frac{\partial {\phi }_4}{\partial n}=(y-b)\cos (nz)-(z-c)\cos (ny),& \frac{\partial {\phi }_4^{\prime }}{\partial n}=0,\\ \frac{\partial {\phi }_5}{\partial n}=(z-c)\cos (nx)-(x-a)\cos (nz),& \frac{\partial {\phi }_5^{\prime }}{\partial n}=0,\\ \frac{\partial {\phi }_6}{\partial n}=(x-a)\cos (ny)-(y-b)\cos (nx),& \frac{\partial {\phi }_3^{\prime }}{\partial n}=0.\end{array}$$

н на поверхности второго
$$\begin{array}{l}\frac{\partial {\phi }_1}{\partial n}=0,\frac{\partial {\phi }_1^{\prime }}{\partial n}=\cos (nx),\\ \frac{\partial {\phi }_2}{\partial n}=0,\frac{\partial {\phi }_2^{\prime }}{\partial n}=\cos (ny),\\ \frac{\partial {\phi }_3}{\partial n}=0,\frac{\partial {\phi }_3^{\prime }}{\partial n}=\cos (nz),\\ \frac{\partial {\phi }_4}{\partial n}=0,\frac{\partial {\phi }_4^{\prime }}{\partial n}=(y-b^{\prime })\cos (nz)-(z-c^{\prime })\cos (ny),\\ \frac{\partial {\phi }_5}{\partial n}=0,\frac{\partial {\phi }_5^{\prime }}{\partial n}=(z-c^{\prime })\cos (nx)-(x-a^{\prime })\cos (nz),\\ \frac{\partial {\phi }_6}{\partial n}=0,\frac{\partial {\phi }_6^{\prime }}{\partial n}=(x-a^{\prime })\cos (ny)-(y-b^{\prime })\cos (nx).\end{array}$$

Эти уравнения, которые соответствуют уравнениям (20), выведены таким же путем, как и полледние. Указанные условия определяют вполне двенадцать функций ${\phi }_1,{\phi }_2,\dots$

Простейшим случаем будет тот, когда оба тела — шары, и точки $(a,b,c)$ и $(a^{\prime },b^{\prime },c^{\prime })$ — их центры. Тогда условия, определяющие ${\phi }_4,{\phi }_5,{\phi }_6$, ${\phi }_4^{\prime },{\phi }_5^{\prime },{\phi }_6^{\prime }$, приведут к тому, что эти шесть функций обращаются в нуль, и, следовательно, мы получим
$$\phi =u{\phi }_1+v{\phi }_2+w{\phi }_3+u^{\prime }{\phi }_1^{\prime }+u^{\prime }{\phi }_2^{\prime }+w^{\prime }{\phi }_3^{\prime }.$$

В этом случае с помощью так называемых шаровых функций можно найти выражение $\phi$ в виде бесконечного ряда, всегда сходящегося и тем быстрее сходящегося, чем больше расстояние между шарами по сравнению с их радиусами. Мы не будем касаться здесь теории шаровых функций, и потому наметим только в общих чертах способ решения названной задачи и определим результат лишь постольку, поскольку он будет необходим для дальнейших исследований.

Шаровые функции непосредственно позволяют найти функцию $U$, которая вне одного шара удовлетворяет дифференциальному уравнению $\mathrm{\Delta }U=0$, однозначна и непрерывна вместе со своими первыми производными, обращается в нуль на бесконечности и на поверхности шара $\frac{\partial U}{\partial n}$ получает произвольно заданное непрерывно изменяющееся значение. Из бесконечного числа подобных функций можно составить сходящийся ряд для потенциала скоростей жидкости, в которой таким образом движутся два шара. Чтобы это показать, назовем шар, центр которого — $a,b,c$, первым, а другой вторым.

Составим функцию $U_1$, для которой на поверхности первого шара будет $\frac{\partial U_1}{\partial n}=\frac{\partial \phi }{\partial n}$, т. е. равно $u\cos (nx)+v\cos (ny)+w\cos (nz)$, и функцию $U_1^{\prime }$, для которой на поверхности второго шара будет $\frac{\partial U_1^{\prime }}{\partial n}=\frac{\partial \phi }{\partial n}$, т. е. равно $u^{\prime }\cos (nx)+v^{\prime }\cos (ny)+w^{\prime }\cos (nz)$, и положим
$$U_1+U_1^{\prime }=V_1\text{. }$$

Тогда это $V_1$ будет первым членом ряда, представляющего ч. Функция $\phi -V_1$ должна для обоих шаров удовлетворять условиям, что на первом
$$\frac{\partial (\phi -V_1)}{\partial n}=-\frac{\partial U_1^{\prime }}{\partial n},$$

и на втором
$$\frac{\partial (\phi -V_1)}{\partial n}=-\frac{\partial U_1}{\partial n}.$$

Составим функции $U_2$ и $U_2^{\prime }$, для которых на первой шаровой поверхности
$$\frac{\partial U_2}{\partial n}=-\frac{\partial U_1^{\prime }}{\partial n}$$

и на второй
$$\frac{\partial U_2^{\prime }}{\partial n}=-\frac{\partial U_1}{\partial n},$$

и положим
$$U_2+U_2^{\prime }=V_2$$

тогда $V_2$ будет вторым членом ряда для $\phi$. Теперь для первой шаровой поверхности должно быть
$$\frac{\partial (\phi -V_1-V_2)}{\partial n}=-\frac{\partial U_2^{\prime }}{\partial n},$$

для второй
$$\frac{\partial (\phi -V_1-V_2)}{\partial n}=-\frac{\partial U_2}{\partial n}.$$

Положим
$$U_3+U_3^{\prime }=V_3,$$

причем $U_3$ и $U_3^{\prime }$ должны быть определены так, чтобы на первой шаровой поверхности
$$\frac{\partial U_3}{\partial n}=-\frac{\partial U_2^{\prime }}{\partial n}$$

и на второй
$$\frac{\partial U_3^{\prime }}{\partial n}=-\frac{\partial U_2}{\partial n}.$$

Продолжая таким образом дальше, мы получим
$$\phi =V_1+V_2+V_3+\dots$$

Этот ряд везде сходится, чего, однако, мы здесь не будем доказывать. Если радиусы обохх шаров рассматривать как бесконечно малые сравни. тельно с расстоянием между шарами, то при этом каждый последующий член ряда будет бесконечно мал сравнительно с предыдущим. Қаждая из величин $V_2,V_3,\dots$ сама выражается в шаровых функциях бесконечным рядом, который обладает также свойством, что каждый следующий член бесконечно мал сравнительно с предыдущим, если радиусы шаров бесконечно малы сравнительно с расстоянием между ними. Если при таком предположении желательно вычислить $\phi$ с точностью только до величин известного порядка, то можно принимать в расчет лишь ограниченное число величин $V$ и для каждой из них ограниченное число членов.

Величины $U_1$ и $U_1^{\prime }$, а также и величина $V_1$ могут быть найдены непосредственно из уравнения (27). Обозначим расстояния точки ( $x,y,z)$ от центров обоих шаров через $r$ и $r^{\prime }$, радиусы их через $R$ и $R^{\prime }$ и воспользуемся тем, что $r$ есть функция $(x-a),(y-b),(z-c),r^{\prime }$ — функция $(x-a^{\prime }),(y-b^{\prime }),(z-c^{\prime })$; тогда получим
$$\begin{array}{l}{U}_1=-\frac{R^3}{2}(u\frac{\partial \frac{1}{r}}{\partial a}+v\frac{\partial \frac{1}{r}}{\partial b}+w\frac{\partial \frac{1}{r}}{\partial c})\\ U_1^{\prime }=-\frac{R^{\prime }}{2}(u^{\prime }\frac{\partial \frac{1}{r^{\prime }}}{\partial a^{\prime }}+v^{\prime }\frac{\partial \frac{1}{r^{\prime }}}{\partial b^{\prime }}+w^{\prime }\frac{\partial \frac{1}{r^{\prime }}}{\partial c^{\prime }}).\end{array}$$

Сумма этих двух выражений, т. е. $V_1$, определяет значение $\phi$ в первом приближении. Примем $R$ и $R^{\prime }$ за бесконечно малые первого порядка и расстояние между системами за конечную величину; тогда для того, чтобы при этом потенциал скоростей и скорость были, воощще говоря, конечными, $u,v,w,u^{\prime },v^{\prime },w^{\prime }$ должны быть бесконечно большими третьего порядка. Тогда на шаровых поверхностях $\phi$ со своими первыми производными должны быть бесконечно велики.

На конечном расстоянии от шаровых поверхностей уравнение $\phi =V_1$ дает скорости до бесконечно малых величин, но не на самих поверхностях. Чтобы получить их на этих поверхностях, надо составить $U_2$, но при этом надо принять в расчет члены высшего порядка. Это пригодит к тому, что $U_2$ и $U_2^{\prime }$ опять надо будет вы числять из уравнения (27). Чтобы найти $U_2$, следует составить для первой шаровой поверхности $\frac{\partial U_1^{\prime }}{\partial n}$, т. е.
$$\frac{\partial U_1^{\prime }}{\partial n}=\frac{\partial U_1^{\prime }}{\partial x}\cos (nx)+\frac{\partial U_1^{\prime }}{\partial y}\cos (ny)+\frac{\partial U_1^{\prime }}{\partial z}\cos (nz).$$

Подставим сюда вместо $U_1^{\prime }$ его значение (28), введем вместо производных по $x,y,z$, проязво ұные по $a^{\prime },b^{\prime },c^{\prime }$, взятые с обратным знаком. Заменим, что допустимо пра требуемой точности, $r^{\prime }$ на $r_0$, где $r_0$ — pacстояние между центрами шаров, следозательно, получим
$$r_0=\sqrt{{(a-a^{\prime })}^2+{(b-b^{\prime })}^2+{(c-c^{\prime })}^2}$$

и заменим производные по $a^{\prime },b^{\prime },c^{\prime }$ производными по $a,b,c$; тогда
$$\begin{array}{l}-\frac{R^{\prime }}{2}(u^{\prime }\frac{{\partial }^2\frac{1}{r_n}}{\partial b\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b^2}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial c})\cos (ny)-\\ -\frac{R^{\prime 3}}{2}(u^{\prime }\frac{{\partial }^2-\frac{1}{r_0}}{\partial c\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_n}}{\partial c^2})\cos (nz).\end{array}$$

Вместе с тем по (27) будем иметь
$$\begin{array}{l}{U}_2=\frac{R^3{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2-\frac{1}{r_0}}{\partial a^2}+v^{\prime }\frac{{\partial }^2-\frac{1}{r_0}}{\partial a\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial z\partial c})\frac{\partial \frac{1}{r}}{\partial a}+\\ +\frac{R^3{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b}+v^{\prime }\frac{{\partial }^2\frac{1}{r_n}}{\partial b^2}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial c})\frac{\partial \frac{1}{r}}{\partial b}+\\ +\frac{R{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_{\mathrm{n}}}}{\partial c\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial {\iota }^2})\frac{\partial \frac{1}{r^2}}{\partial c}.\end{array}$$

Соответственно будет
$$\begin{array}{rl}{U}_2^{\prime }& =\frac{R^3{R}^3}{4}(u\frac{{\partial }^2\frac{1}{r_n}}{\partial a^2}+v\frac{{\partial }^2\frac{1}{r_0}}{\partial a\partial b}+w\frac{{\partial }^2\frac{1}{r_n}}{\partial a\partial c})\frac{1}{\frac{1}{r^{\prime }}}+\\ & +\frac{R^3{R}^{\prime 3}}{4}(u\frac{{\partial }^2\frac{1}{r_n}}{\partial b\partial a}+v\frac{{\partial }^2\frac{1}{r_0}}{\partial b^2}+w\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial c})\frac{\partial \frac{1}{r^{\prime }}}{\partial b^{\prime }}+\\ & +\frac{R^3{R}^{\prime 3}}{4}(u\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial a}+v\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial b}+w\frac{{\partial }^2\frac{1}{r}}{\partial c^2})\frac{\partial \frac{1}{r}}{\partial c^{\prime }}.\end{array}$$

Составим с помощью этих значений $U_2$ и $U_2^{\prime }$ уравнение $\phi =V_1+V_2$, тогда оно даст также и на шаровой поверхности скорости с точностью до бесконечно малых величин и значение $\phi$ со включением бесконечно малых величин первого порядка. Мы будем выполнять вычисление так. чтобы $\phi$ на каждой шаровой поверхности было определено с такой же точностью. Чтобы найти $\phi$, надо еще преоб разовать значения $U_1$ и $U_1^{\prime }$, определяемые (28). Займемся определением $\phi$ для первой шаровой поверхности. Величину $\frac{1}{r^{\prime }}$, входящую в выражение для $U_1^{\prime }$, разложим по степеням $x-a,y-b,z-c$, которые являются бесконечно малыми первого порядка, причем необходимо принять во внимание только первые степени. Так кая
$$x-a=R\cos (nx),y-b=R\cos (ny),z-c=R\cos (nz),$$

то поэтому можно положить
$$\frac{1}{r^{\prime }}=\frac{1}{r_0}+R[\frac{\partial \frac{1}{r_0}}{\partial a}\cos (nx)+\frac{\partial \frac{1}{r_0}}{\partial b}\cos (ny)+\frac{\partial -\frac{1}{r_0}}{\partial c}\cos (nz)]$$

отсюда получим
$$\begin{array}{c}{U}_1^{\prime }=\frac{R^{\prime 3}}{2}(u^{\prime }\frac{\partial -\frac{1}{r_0}}{\partial a}+v^{\prime }\frac{\partial \frac{1}{r_0}}{\partial b}+w^{\prime }\frac{\partial \frac{1}{r_0}}{\partial c})+\\ +\frac{R{R}^{\prime 3}}{2}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a^2}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a\partial c})\cos (nx)+\\ +\frac{R{R}^{\prime 3}}{2}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b^2}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial c})\cos (ny)+\\ +\frac{R{R}^{\prime 3}}{2}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c^2})\cos (nz).\end{array}$$

Далее из (28) следует
$$U_1=-\frac{R}{2}[u\cos (nx)+v\cos (ny)+w\cos (nz)],$$

так как
$$\frac{\partial \frac{1}{r}}{\partial a}=\frac{1}{R^2}\cos (nx),\frac{\partial \frac{1}{r}}{\partial b}=\frac{1}{R^2}\cos (ny),\frac{\partial \frac{1}{r}}{\partial c}=\frac{1}{R^2}\cos (nz).$$

Воспользуемся этими соотношениями, чтобы составить также определен. ное (29) выражение для $U_2$, и заметим, что $U_2^{\prime }$ высшего порядка малости, чем $U_2$; тогда для первой шаровой поверхности с требуемой степеньк точности мы получим
$$\begin{array}{c}\phi =-\frac{R}{2}[u\cos (nx)+v\cos (ny)+w\cos (nz)]+\\ +\frac{R^{\prime 3}}{2}(u^{\prime }\frac{\partial \frac{1}{r_0}}{\partial a}+v^{\prime }\frac{\partial \frac{1}{r_0}}{\partial b}+w^{\prime }\frac{\partial \frac{1}{r_0}}{\partial c})+\\ +\frac{3R{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a^2}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial a\partial c})\cos (nx)+\\ +\frac{3R{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b^2}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial b\frac{\partial c}{}})\cos (ny)-\\ +\frac{3R{R}^{\prime 3}}{4}(u^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial a}+v^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c\partial b}+w^{\prime }\frac{{\partial }^2\frac{1}{r_0}}{\partial c^2})\cos (nz).\end{array}$$

Значение $\phi$ на поверхности второго шара мы найдем, если переставим в этом выражении буквы со штрихами и буквы без штрихов.

Более общая задача, чем излсженная здесь, разрешена Бьеркнесом (Bjerknes) в его мемуаре: «Sur le mouvement simultané de corps sphériques variables dans un fluide indéfini et incompressible, presenté à la société des Sciences de Christiania le 15 sept. 1871».