数理解析研究所講究録922 巻 1995 年 124-131

124

回転液滴内の流れと変形

農工大・工・応用物理

1.

佐野理

$($

$\mathrm{O}_{\mathrm{S}_{\iota}\mathrm{a}\mathrm{l}}\mathrm{n}n\mathrm{s}_{\mathrm{a}}\mathrm{n}()$

はじめに静止状態で中立浮遊した液滴は表面張力により球形を保っているが、 これに回転を与えていくと次第に偏平になり、 いくつかの形態遷移を経たのちにの実験 1) をはじめとして非常に古い分裂する。 この問題は 19 世紀の歴史をもっているだけでなく、 宇宙での天体の形との類似性から、 Newton,$\mathrm{P}\mathrm{l}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{a}\iota 1$

Maclaurin, Jacobi, Poincare,... など多くの研究者の関心を集めてきた 2)O もっとも後者は自己重力と遠心力などで決まる形であるのに対して、 われわれの系は表面張力と遠心力で決まる形である点が異なっている。 液温の平衡形状についても Chandrasekhar3) による理論計算、 Brown and Scriven4) による数値計算などがあるが、 これらはいずれも液滴の粘性境界条件を無視した取扱である。 最近、 NASA スペースシャトル内での実験結果が Wang5) 達によって報告されている。 これは微小重力下での実験であり、 条件としては望ましいが、液滴内部の流れなどの詳細な測定は行なわれていない。我々は、 地上の実験室ではあるがほぼ中立浮遊の回転液滴を実現させ、形態観測と共に内部の流体の流速分布をレーザー流速計を用いて精密に測定した。 以下ではこれらについて報告し、 流れの機構を流体力学の基礎方程式に基づいて解析する。

2. 実験方法と結果の概要
実験装置および測定の概略を図 1 に示す。 液滴としてはオルトトルイジンを使用した。 これは周囲の液体 (蒸留水) と混じり合わないだけでなく、 密度の温度依存性が高い。 そこで温度を調節することによってほぼ中立浮遊状態の液滴を実現させた。 また、 液滴の上端にモーターで駆動した棒を接触させ、 液滴に任意の回転角速度を与えるように設計した。 モーターの回転数は印加電圧に比例するので、 これをコンピューターで制御する。 実験では到達125

すべき角速度と角加速度の両方を変化させて系統的に形の変化を調べた。
測定には、 液滴の運動に撹乱を与えないようにする必要があるので、 レー
ザー流速計を使用し、 またビデオカメラで形態の変化を記録し画像解析を行なった。

図1

実験装置の概略

以下に主な結果をまとめて述べよう。
(1) 図 2 は液滴の平衡形状の例である。 一般に液滴の形は球\rightarrow 回転楕円体\rightarrow 二分裂形状\rightarrow 西洋梨型と遷移するのが観測された。

図2

液滴の形態遷移126(2) 液滴の回転が比較的遅く、 球または回転楕円体形状が保たれている範囲内で内部の流速を測定した。 図 3 は液滴の赤道面内での周方向速度成分を回転軸からの距離 を変えて測った結果である。 回転角速度が比較的$\mathrm{r}$$\mathrm{v}_{\psi}$小さい場合には流速分布は剛体回転的である (すなわち速度の大きさが中心からの距離に比例して増加する) が、 回転数の増加と共に直線から外れてくる。 このずれは液滴の中心や周辺では$0$であり、 その中間付近で最大である。 さらにずれの最大値は角速度$\Omega$$\mathrm{r}-0.5$の 2 乗に比例して増加している。

$\cup \mathrm{l}\mathrm{s}\mathrm{r}\mathrm{a}\mathrm{n}\mathrm{c}\mathrm{e}r\iota \mathrm{c}\mathrm{m}/\mathrm{S}\rfloor$

図3

赤道面内での周方向速度成分$\mathrm{v}_{\emptyset}$の回転角速度依存性(3) 液滴の回転は–般にモーターの回転より遅い。 これは周囲の流体との粘性摩擦のためである。 しかし球または回転楕円体形状が保たれている範囲では両者はほぼ比例している。 液滴が二分裂形状に遷移すると、 液滴の回転数は著しく減少する。
(4) 液滴の形が回転楕円体から二分裂形状に遷移する臨界角速度を測定した結果を図 4 に示す。 この臨界値 In は液滴の角加速度にも依存する。 そこ$\omega$で角加速度を変化させ、 それぞれの場合に求めた臨界角速度を角加速度$0$とした。 このようにして求の場合に外挿し、 これをもって臨界角速度めた In* と液滴の体積との関係を図 5 に示す。小さな液滴では表面張力の$\omega_{\mathrm{m}}*$$\omega$

影響が大きいので回転角速度をかなり大きくしなければ遷移が起こらないが、 大きな液滴では遠心力の影響が大きいので回転角速度をかなり小さく127しても遷移が起こる。 ちなみに実験的には$\omega_{\mathrm{m}}*$は体積 V の -0.16 乗に比例しているようである。

図4

臨界角速度の角加速度依存性

図5

臨界角速度の義民体積依存性

3. 流れの理論解析
液面内部の流れを Navier-Stokes 方程式と連続の方程式に基づいて解析してみよう。 ただし液滴内の流体は粘性率
液体外部の流体は粘性率,t’の非圧縮性流体とし、 球座標系 $(r, \Theta, \phi)$ で考える。 いま液滴は外部からのトルクによって–定の角速度 で回転していると仮定すると、 基礎方程’$\mu_{\text{、}}$$\omega$

式および境界条件は)$\frac{1}{\mathrm{r}^{2}}\frac{\partial}{\partial \mathrm{r}}\langle_{\mathrm{r}^{2}\mathrm{V}}\mathrm{r}$$+ \frac{1}{\mathrm{r}\sin\theta}\frac{\partial}{\partial\theta}(_{\mathrm{s}\mathrm{i}\theta \mathrm{v}_{\ominus}}\mathrm{n})=0$

,

(1a)

$\Delta \mathrm{v}_{\mathrm{r}^{-}\mathrm{r}}\frac{2}{\mathrm{r}^{2}}\mathrm{v}-\frac{2}{\mathrm{r}^{2}}\frac{\partial \mathrm{v}_{\Theta}}{\partial\theta}-\frac{2\mathrm{v}_{\Theta}}{\mathrm{r}^{2}}\mathrm{c}\mathrm{o}\mathrm{t}\mathrm{e}-\frac{\partial \mathrm{p}}{\partial \mathrm{r}}={\rm Re}(\mathrm{v}_{\mathrm{r}^{\frac{\partial \mathrm{v}_{\mathrm{r}}}{\partial \mathrm{r}}}}+\frac{\mathrm{v}_{\Theta}}{\mathrm{r}}\frac{\partial \mathrm{v}_{\mathrm{r}}}{\partial\theta}-\frac{\mathrm{v}_{\Theta}^{2_{+\mathrm{V}}2}\varphi}{\mathrm{r}})$

(1b)

,

$\Delta \mathrm{v}_{\Theta^{-}}\frac{\mathrm{v}_{\Theta}}{\mathrm{r}^{2}\sin^{2}\theta}+\frac{2}{\mathrm{r}^{2}}\frac{\partial \mathrm{v}_{\mathrm{r}}}{\partial\theta}-\frac{1}{\mathrm{r}}\frac{\partial \mathrm{p}}{\partial\theta}=\mathrm{R}\ovalbox{\tt\small REJECT} _{\mathrm{f}}\frac{\partial \mathrm{v}_{\Theta}}{\partial \mathrm{r}}+\frac{\mathrm{v}_{\Theta}}{\mathrm{r}}\frac{\partial \mathrm{v}_{\Theta}}{\partial\theta}+\frac{\mathrm{v}_{\mathrm{r}}\mathrm{v}_{\Theta}}{\mathrm{r}}-\frac{\mathrm{v}^{2}(1)\cot\theta}{\mathrm{r}})$,$(1\mathrm{c})$,$\Delta \mathrm{v}_{\phi)}-\frac{\mathrm{v}_{\phi}}{\mathrm{r}^{2}\sin^{2}\theta}={\rm Re}(\mathrm{V}_{\mathrm{r}}\frac{\partial \mathrm{v}_{\phi}}{\partial \mathrm{r}}+\frac{\mathrm{v}_{\Theta}}{\mathrm{r}}\frac{\partial \mathrm{v}_{\phi}}{\partial\theta}+\frac{\mathrm{v}_{\mathrm{r}}\mathrm{v}_{\phi}}{\mathrm{r}}+\frac{\mathrm{v}_{\Theta}\mathrm{V}_{\phi}\mathrm{c}\mathrm{o}\mathrm{t}\mathrm{e}}{\mathrm{r}}+\frac{2}{\mathrm{R}\mathrm{o}}(\mathrm{v}_{\mathrm{r}}\sin\theta+\mathrm{v}_{\Theta}\cos\theta)1$$(1\mathrm{d})$128となる$0$

ただし\Deltaはラプラシアン、 $Re=Ua/\mathrm{v}$$= \frac{1}{\mathrm{r}^{2}}\frac{\partial}{\partial \mathrm{r}}(\mathrm{r}2_{\frac{\partial}{\partial \mathrm{r}}})+\frac{\mathrm{l}}{\mathrm{r}\sin\theta}2\frac{\partial}{\partial\theta}(\sin \mathrm{e}\frac{\partial}{\partial\theta})$と $Ro=U/\mathrm{e}\mathrm{o}a$ はそれぞれ Reynolds 数および ROssby 数であり、は液滴の半径、 は動粘性率である。ここでは定常軸対称流を考えているので、 時間および方位角 に関する微分は省いてある。 式 (ld) の右辺括弧内第 2 項はCoriolis 力を表わす。 液滴は自由表面であるが表面張力が充分大きくて球形が保たれていると仮定し、 境界条件は液滴表面で速度および応力が連続であるとする。 以下では $Re$ が小さいと仮定し、 解を$\mathrm{a}$$fJ$$\mathrm{t}$$\phi$$\mathrm{r}^{=}1$$\mathrm{v}=\mathrm{v}^{(\{))}+{\rm Re} \mathrm{v}^{(1)}+\ldots.$.,$\mathrm{p}=\mathrm{p}^{(0)}+{\rm Re} \mathrm{p}^{(1)}+..,.$.

$(2 \mathrm{a}, \mathrm{b})$と展開する。

3.1$\mathrm{o}({\rm Re})0$の解

このオーダーの解は$\mathrm{v}_{\mathrm{r}}^{(0)}=\mathrm{v}_{\Theta}^{(0)}=0$,$\mathrm{p}^{(0)}=2\mathrm{e})(=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{S}\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{t})$$\mathrm{v}_{\phi}^{(0)}=\frac{1}{\mathrm{R}\mathrm{o}}\mathrm{r}\sin\theta$$= \frac{1}{\mathrm{R}\mathrm{o}}\frac{1}{\mathrm{r}^{2}}\sin\theta$for,,$(3\mathrm{a},\mathrm{b}_{\mathrm{C})}$$\mathrm{r}\leq 1$

for$(3\mathrm{d},\mathrm{e})$$\mathrm{r}\geq 1$で与えられる。 式 (3d,e) はそれぞれ剛体回転による内外の流速を表わしている。 ここで、$\sigma$は表面張力であり、 これが充分大きくて液滴は球形が保たれていると仮定している。 また、 上記の運動が保たれるためには外部からのトルクが与えられている必要がある。 以下では、 液滴内部の解のみ示す。

$8\pi \mathrm{a}^{3}\mu\omega \mathrm{I}$3.2$\mathrm{o}({\rm Re}^{1})$の解

このオーダーの方程式は$\frac{1}{\mathrm{r}^{2}}\frac{\partial}{\partial \mathrm{r}}(_{\mathrm{r}^{2}\mathrm{v}_{\mathrm{r}}^{()}}\iota)+\frac{1}{\mathrm{r}\sin\theta}\frac{\partial}{\partial\theta}\{\sin\theta \mathrm{V})(\Theta=01)$,$(- 4\mathrm{a})$$\Delta \mathrm{V}_{\mathrm{r}}^{(1)}-\frac{2}{\mathrm{r}^{2}}\mathrm{V}_{\mathrm{r}}^{(})-1\frac{2}{\mathrm{r}^{2}}\frac{\partial \mathrm{v}_{\Theta}^{(1)}}{\partial\theta}-\frac{2\mathrm{v}_{\Theta}^{(1)}}{\mathrm{r}^{\mathrm{A}}\circ}\cot\theta-\frac{\partial \mathrm{p}^{(1)}}{\partial \mathrm{r}}=-\frac{1}{\mathrm{r}}(\mathrm{v}_{\phi}^{(0)})^{2}$,$\Delta \mathrm{v}_{\Theta}^{(1\rangle}-\frac{\mathrm{v}_{\Theta}^{(1\rangle}}{\mathrm{r}^{2}\sin^{2}\theta}+\frac{2}{\mathrm{r}^{2}}\frac{\partial \mathrm{v}_{\mathrm{r}}^{()}1}{\partial\theta}-\frac{1}{\mathrm{r}}\frac{\partial \mathrm{p}^{(1)}}{\partial\theta}=-\frac{\cot\theta}{\mathrm{r}}(\mathrm{v}_{\phi}^{(})^{2}0)$

,

(4b)

(4c)

129

$\Delta \mathrm{v}_{\phi}^{(1)}-\frac{\mathrm{v}_{\phi}^{(1)}}{\mathrm{r}^{2}\mathrm{s}\mathrm{i}\mathrm{I}1^{\gamma}\sim\theta}=0$.

$(4\mathrm{d}\rangle$である。 式 (4d) の解のうち、 液滴中心で正則かつ境界で(1)$\mathrm{v}_{\phi}$$=0$$0$となるものはだけである。 また、 式 (4a-4c) の非同次の特解は$\mathrm{v}_{\mathrm{r}}^{(1)}\mathrm{p}=\mathrm{v}_{\Theta}^{(1\mathrm{p})}=0,$$\mathrm{p}^{(1)}\mathrm{P}=\frac{1}{2\mathrm{R}\mathrm{o}^{2}}\mathrm{r}^{2}\sin^{2\theta}$.

(5)である。 ここで上付き添字の第 1 番目は Re の次数を、 第 2 番目は解の性格を表わすものとする。式 (4a-4c) の同次の–般解のうち液滴内部で有界であり境界で(i)$0$となるものを、 解の対称性を考慮していくつか示す。

1 セル型 .
$\cdot$A$\mathrm{v}_{\mathrm{r}}^{(11)}=2$$0^{(11)}=- 20$

ここで $A$$(1 - \mathrm{r}^{2})\cos\theta$A,$\mathrm{v}_{\Theta}^{(11)}=- 2$

A(1-2r ) sin6,$\mathrm{v}_{\phi}^{(11)}=0$,$(6\mathrm{a},\mathrm{b},\mathrm{c}\rangle$.

$\mathrm{r}\cos\theta$(6d)は任意定数である。 また$\mathrm{v}_{\mathrm{r}}=\frac{\mathrm{l}}{\mathrm{r}\sin\theta}\frac{\partial\psi}{\partial\theta}2$,$\mathrm{v}_{\Theta}=-\frac{1}{\mathrm{r}\sin\theta}\frac{\partial\psi}{\partial \mathrm{r}}$で定義される流れの関数$\eta^{(11)}’=\mathrm{A}\mathrm{r}^{2}(1.

$\psi$-\mathrm{r}^{2})\sin^{2}\theta$(7)を用いると、 式 (6) は.

(8)と表わされる。 この流れは Hill の球形渦と同じ表現を持っており、 その概略は図 6 に示したようなものである。

図61 セル型の流れ

図72 セル型の流れ

130(i) 2$\cdot \mathrm{b}[]\triangleright$

型..

$\mathrm{v}_{\mathrm{r}}^{(12)}=4\mathrm{B}\mathrm{r}(1- \mathrm{r}^{2})\mathrm{P}_{2}(\cos\theta)$,$\mathrm{p}^{(12)}=$

ここで$\mathrm{P}_{2}(\cos\theta)=\frac{1}{4}(1+3\cos 2\theta)$$\mathrm{v}_{\mathrm{e}^{12}}^{()}=- \mathrm{B}\mathrm{r}(3- 5\mathrm{r}^{2})\sin 2\theta$

-28B

,,$B$

は任意定数である。

$-$.となる。 図 7 にその概形を示す。

(1 )$\psi^{(12)}=2\mathrm{B}\mathrm{r}^{3}-(1- \mathrm{r}^{2})\sin^{2}\theta\cos\theta$$\mathrm{o}({\rm Re}^{2})$$\mathrm{v}_{\phi}^{(\iota 2)}=0$(.9 a-d)$\mathrm{r}^{2}\mathrm{p}2(\cos\theta)$

は Legendre 関数であり、流れの関数を用いるとこれは3.3,$()$の解このオーダーの解は)$\mathrm{v}_{\mathrm{r}}^{(2)}=\mathrm{v}_{\Theta}=\langle\angle$$\mathrm{P}^{(}2$) $=0$ ,$(11 \mathrm{a}, \mathrm{b}, \mathrm{c})$および$\mathrm{v}^{(2)}$$\partial \mathrm{v}^{(^{(})})$$\mathrm{v}^{(1)\partial}\mathrm{v}^{(}\mathrm{t}))$$\mathrm{v}^{(1)}\mathrm{v}(0_{)}$$\mathrm{v}\mathrm{c}(1)(())\mathrm{o}\mathrm{t}\theta$$\Delta \mathrm{v}_{\phi}^{(_{\sim}^{\gamma_{)}}}-\frac{\phi}{\mathrm{r}^{2}\sin^{2}\theta}=\mathrm{v}_{\mathrm{r}^{1}}^{()_{\frac{\phi}{\partial \mathrm{r}}}}+\frac{\Theta}{\mathrm{r}}\frac{\phi}{\partial\theta}+\frac{\mathrm{r}\phi}{\mathrm{r}}+\frac{\Theta\phi}{\mathrm{r}}+-2\mathrm{R}\overline{\mathrm{o}}(\mathrm{v}_{\mathrm{r}}^{()}\sin\theta+\cos\theta\iota(1))\Theta$(1(i) 式 (3d) および$\mathrm{v}_{\phi}^{(2)}=()$$(6\mathrm{a},\mathrm{b})$$2\rangle$

を式 (12) に代入し r $=()$ で正則かつ境界条件 $r=1$ でを満たすものを求めると$\mathrm{v}_{\phi}^{()}21=-\frac{4A}{7R_{\mathit{0}}}\mathrm{r}^{2}(1- \mathrm{r}^{2})\sin\theta\cos\theta$(ii) 式 (3d) および$\mathrm{v}_{\phi}^{(2)}=(\mathrm{J}$’$(9\mathrm{a},\mathrm{b})$

.

(13)

を式 (12) に代入し r $=()$ で正則かつ境界条件 $r=1$ でを満たすものを求めると「( l-r2)$\mathrm{v}_{\phi}^{\langle 22)}=\frac{B}{R_{\mathit{0}}}\mathrm{r}\sin\theta$ $\frac{2}{5}-\frac{16}{15}\mathrm{r}\mathrm{L}2+\frac{2}{3_{-}}\mathrm{r}^{\mathfrak{l}_{-}}49\Gamma^{2}$(14)$\cos 2\theta\rfloor$を得る。

液滴の赤道面内での速度分布は(15)$\mathrm{v}_{\phi}=\mathrm{v}_{\varphi}^{(0)}+{\rm Re}(_{\phi}^{(}221\rangle+\mathrm{V}\mathrm{v}_{\mathfrak{p}(}^{(_{\grave{\mathrm{A}}}2)})+\ldots$.

で与えられる。 –例として分布を図 8 に示す。

$R_{\mathit{0}}=2,$,$R\epsilon\supset=1()0,$$B=2\mathrm{x}10^{-4}$の場合の速度1314議論

我々の理論で予想される速度分布の剛体回転からのずれは$\mathfrak{d}\mathrm{v}_{\phi}={\rm Re}^{22)}\mathrm{v}_{\phi}^{(2}+\ldots=\frac{\mathrm{B}{\rm Re}^{2}}{\mathrm{R}\mathrm{o}}\mathrm{F}(\mathrm{r})+\ldots$

ただし$\mathrm{F}(\mathrm{r})=\mathrm{r}$(1,$6\mathrm{a}\rangle$$(16\mathrm{b}\rangle$$\llcorner||^{\frac{2}{-5}-}|\frac{28}{4_{-}5}\mathrm{r}^{2}+\frac{2}{9}\mu]$である。 これを実験結果と比較してみよう。 関数 $F(r)$ は $r=0.5()212$ で最大値 0.12917 をとる。 このピーク値とモーター回転数の関係を図 9 に示す。
ただし、 実験条件を鑑みてら明らかに$\Omega_{\mathrm{m}_{\sim}^{\langle}}25$

[$)$

[

$\mathrm{r}\mathrm{p}_{\mathrm{l}}\mathrm{n}|$$\Omega_{\mathrm{m}}$で\mbox{\boldmath $\delta$}V\mbox{\boldmath $\phi$}は$R\mathrm{e}$に比例し $Ro$ は–定であるとした。 図か$\max/\Omega_{\mathrm{m}}^{2}=$一定となっているいることがわかる。 すなわち、 液滴内部の流れは COriolis 力に起因すると考えてよさそうである ($25()\iota_{\mathrm{r}\mathrm{D}\ln}|$

図8以上の回転数でこの理論からはずれるのは液滴が非対称な赤道面内の流れ

図9

ずれの大きさとモーター回転数

参考文献
1) Hildebrandt,S. and Trolmba,A. : Mathematics$\partial I\mathit{1}d$Optim$\mathrm{a}\mathit{1}$Form,

Sci. Am Library, 14 (Freemann, New York, 1985).
2) 例えば 蜂巣 泉: 京都大学数理解析研究所講究録 510(1984)
3) Chandrasekhar,S. :Proc $.R(yy.soc$ $.LoIld.$ , A286 (1965) 1-26.
$.Lond$ .
(1980) 331-357.
:Proc5) Wang,T.G., Trinh,E.H., Croonquist,A.P. and Elleman,D.D. : Phys.Rev.Lett.

4) Brown .R.A. and Scriven,

56 (1986)452-455.

$\mathrm{L}.\mathrm{E}.$$.R\mathrm{o}_{\nu}V.soc$$,\mathrm{A}371$