数理解析研究所講究録 1434 巻 2005 年 3-14

32 円柱を横切る流れの不安定性と遷移
同志社大学工学部
同志社大学工学部
同志社大学工学部

秋永

剛

勇祐
水島 二郎
関(Takeshi Akinaga)(Yusuke Seki)(Jiro Mizushima)

1. はじめに流れの中に円柱がおかれているとき, かなり低いレイノルズ数において定常流から振動流への遷移がみられる. 流れが振動流となると, 円柱が長期にわたり周期的な力を受けて破断に至る現象も報告されており, 流れの振動を制御する技術が研究されてきた. Jackson[l] は 1 本の円柱を過ぎる流れの安定性解析を行い, 円柱を過ぎる流れはレイノルズ数が 454 以上のとき不安定となり対称な定常流から周期流へと遷移することを明らかにした.
流れと直角に並べられた 2 本の円柱を過ぎる流れにおいて, 2 本の円柱間隔がおよそ円柱直径を超える場合, 2 円柱の後流は一対の特徴的な流れとなることが実験的に確認された. また, その流れはレイノルズ数や円柱間隔の大きさによって数種類に分類されることが示された.
円柱間隔が円柱直径程度よりも小さい場合, 高いレイノルズ数 $(Re=10^{3}\sim 10^{4})$において 2 本の円柱間の流れは対称でなく一方に偏ることが確認されている [2].
Xu, Zhou and So[3] はこのような狭い円柱画面の場合について, 熱線流速計や PIVを用いた実験を行い, 円柱間隔の大きさが流れ構造へ与える影響を詳しく調べた.
その結果, レイノルズ数が低いときは, 円柱の下流では 1 つの底地を形成するのに対し, ある値を超えると円柱下流で 2 つの異なる周期の沁沁が発生することを発見し, 臨界レイノルズ数と円柱間隔について整理を行った.
円柱間隔が広い場合についても研究が行われている. Wiiliamson[4] は比較的低いレイノルズ数 $(Re=50\sim 200)$ で円柱間距離がおよそ円柱半径から円柱直径の5 倍の範囲で可視化実験を行い, 円柱間隔が円柱の直径を超えるとき, 2 つの円柱の後流は同期して同位相と逆位相の 2 つの流れパターンが発生することを発見ししかし, 彼はこれらの流れの発生の原因を特定することができなかった.
これに対し, Peschard and Le Ga1[5] は, 2 つの後流の振動を表すモデルとして,た.

連立振幅方程式からなる力学系モデルを提案し, このモデルにより数値シミュレーションを行うと, 2 円柱後方から放出される渦が逆位相モードと同位相モードとなることをに示した. また, これら 2 つのモードの発生は円柱間隔に依存しており, 円柱間隔を大きくしていくことにより, 同位相モードから逆位相モードへと遷移していくことを明らかにした.
最近の研究では, 比較的低いレイノルズ数 $(Re=40\sim 160)$ において, 流れはレイノルズ数や円柱間隔および初期条件などに依存して, 6 つに分類できるという数値シミュレーションの結果も報告されている [6].

4

図 1:$s^{*}:$

2 円柱を横切る流れ (実際の長さに比例していない模式図).

$d^{*}:$

円柱直径,円柱間隔.
しかし, いずれの研究においても,流れの分類は瞬時の流れ場のもつ特徴や形に基づいて行われており, 流れの起こる原因については詳しく調べられていない.
また, それらが遷移するときの臨界レイノルズ数や遷移構造も明らかになっていない. 流れの遷移が起こると流れはより複雑なものとなる. したがって, 遷移したあとの流れがもとの流れからどのように変わったかという観点から分類を行うことが, 最も重要な方法であると考えられる.
本研究では, 一様流流に流れと垂直に並べられた 2 本の円柱を過ぎる流れにつと流れの安定性解析および PIV(Particle Image Veいて, 数値シミュレーシ$\text{ョ}\grave{}$$1\mathrm{o}\mathrm{c}\mathrm{i}_{1}\mathrm{n}\mathrm{e}\mathrm{t}\mathrm{r}\mathrm{y})$装置を用いた実験を行い, 流れの遷移について詳しく調べ, これまで報告されているそれぞれの流れパターンが生じる物理的機構を明らかにする.

2. 基礎方程式
2.1 基礎方程式および境界条件
流れと垂直に並べておかれた 2 本の円柱を過ぎる流れを考える (図 1). 2 本の円軸をとる. 上流に軸をとり,軸と垂直に柱間の中心を通り流れと平行にであおいて流れは速度の一様流であると仮定する. 2 円柱の直径は等しくるとし, 円柱表面間の距離をとする. 流れを特徴づける無次元パラメータとしは流体を定義する. ただしおよびて, レイノルズ数の動粘性係数である. ここでは, $Re=10\sim 200$ および $\Gamma=0.3\sim 5.0$ の範囲で調べることとする.
流れは 2 次元非圧縮であると仮定し, 流れ関数 と渦度 を導入する. 流れの$x^{*}$$x^{*}$$y^{*}$$d^{*}$$U^{*}$$s^{*}$$Re\equiv U^{*}d^{*}/\nu^{*}$$\iota/^{*}$$\Gamma\equiv s^{*}/d^{*}$$\psi$$\omega$

5
基礎方程式は, 渦度輸送方程式(1)$\frac{\partial\omega}{\partial t}=\frac{1}{Re}\lambda 4\omega+N(\psi, \omega)$と流れ関数についてのポアソン方程式$\psi$$\mathcal{M}\psi=-\omega$$\mathcal{M}\equiv\frac{\partial^{2}}{\partial’x^{2}}+\frac{\partial^{2}}{\partial y^{2}},$

(2)

,(3)

$N( \psi,\omega)\equiv\frac{\partial\psi}{\partial x}\frac{\partial\omega}{\partial y}-\frac{\partial\psi}{\partial y}\frac{\partial\omega}{\partial x}$

である.

および
流れは計算領域の上流 (図 1, AB) と側面 (図 1,) において一様
ではゾンマーフェであるとする. 円柱上で滑りなし条件を課し, 流出口 (図 1,)
ルト放射条件を用いる.
$\mathrm{A}\mathrm{D}$$\mathrm{B}\mathrm{C}$$\mathrm{C}\mathrm{D}$2.2 定常解
流れは低いレイノルズ数のとき対称である. 定常解方程式および定常ポアソン方程式を満たす.
$0= \frac{1}{Re}\mathcal{M}\overline{\omega}+N(\overline{\psi},\overline{\omega})$$\mathcal{M}\psi^{\overline{/}}=$

式 (4) と (5) を境界条件のもとに, 定常解$\overline{\omega}$

一$(\overline{\psi},\overline{\omega})$は次の定常渦度輸送(4),(5).

$(\overline{\psi},\overline{\omega})$を数値的に求める.

3. 数値計算法
3.1 一般座標
複数物体を過ぎる流れの数値シミュレーションと定常解の数値計算を行うため,物体に沿った格子を数値的に形成する。 このため直交座標系から一般曲線座標系への座標変換 $x=x(\xi, \eta),$ $y=y(\xi, \eta)$ を行う. このとき空間微分は次のように変換する.
$\frac{\partial}{\partial x}=$$\frac{1}{J}(\frac{\partial y}{\partial\eta}\frac{\partial}{\partial\xi}-\frac{\partial y}{\partial\xi}\frac{\partial}{\partial\eta})$,$\frac{\partial}{\partial y}=\frac{1}{J}(\frac{\partial x}{\partial\xi}\frac{\partial}{\partial\eta}-\frac{\partial x}{\partial\eta}$

釘

たとえば, 渦度輸送方程式 (1) およびポアソン方程式 (2) に現れるオペレーターを一般座標に変換すると, 次のようになる.と亙(6)$\mathcal{M}$$(\psi, \omega)$,$\mathcal{M}=\frac{1}{J^{2}}(\alpha\frac{\partial^{2}}{\partial\xi^{2}}-2\beta\frac{\partial^{2}}{\partial\xi\eta}+\gamma\frac{\partial^{I2}}{\partial\eta^{2}}+\sigma\frac{\partial}{\partial\xi}+\tau\frac{\partial^{r}}{\partial\eta})$(7)

$\epsilon$

図 2: 数値計算に用いた計算格子例. $\Gamma=1.0$ .

$N(’ \emptyset, \omega)=-\frac{1}{J}(\frac{\partial\psi}{\partial\eta}\frac{\partial\omega}{\partial\xi}-\frac{\partial\psi}{\partial\xi}\frac{\partial\omega}{\partial\eta})$

(8)

.

ここで,
$\alpha=(\frac{\partial x}{\partial\eta})^{z}+(\frac{\partial y}{\partial\eta})^{z}$$J= \frac{\partial x}{\partial\xi}\frac{\partial y}{\partial\eta}-\frac{\partial x}{\partial\eta}\frac{\partial y}{\partial\xi},$$\gamma=(\frac{\partial x}{\partial\xi})^{2}+(\frac{\partial y}{\partial\xi})^{2}$$\beta=\frac{\partial x}{\partial\xi}\frac{\partial x}{\partial\eta}+\frac{\partial y}{\partial\xi}\frac{\partial y}{\partial\eta},$,$Dx= \alpha\frac{\partial^{2}}{\partial\xi^{2}}-2\beta\frac{\partial^{2}x}{\partial\xi\partial\eta}+\gamma\frac{\partial^{2}xx)}{\partial\eta^{2}’}\sigma=\frac{1}{J,x}(\frac{\partial x}{\partial\eta}Dy-\frac{\partial y}{\partial\eta}D)\tau=\frac{1}{J}(\frac Dx-y)Dy=\alpha\frac{\partial y\partial\xi\partial^{2}y}{\partial\xi^{2}}-2\beta\frac{\frac\partial xD\partial\xi\partial^{2}y}{\partial\xi\partial\eta}+’\gamma\frac{\partial^{2}y}{\partial\eta^{2}}$’$\mathit{1}\mathit{1}X=\alpha-\cdot\angle lJ+\gamma\overline{\partial\xi^{2}}\overline{\partial\xi\partial\eta}\overline{\partial\eta^{2}}$$\mathit{1}/y=\alpha_{\overline{\partial\xi^{2}}}-\overline{\angle}\mathrm{P}_{\overline{\partial\xi\partial\eta}}+\gamma_{\overline{\partial\eta^{2}}}$である.

物体に沿った格子を数値的に形成には, 次のポアソン方程式の解,$\frac{\partial^{2}\xi}{\partial x^{2}}+\frac{\partial^{2}\xi}{\partial y^{2}}=P(\xi, \eta)$$\frac{\partial^{2}\eta}{\partial x^{2}}+\frac{\partial^{2}\eta}{\partial y^{2}}=Q(\xi, \eta)$.

$(\xi, \eta)$を用いる.

(9)

ここで,上式の右辺の関数 $P(\xi, \eta)$ と $Q(\xi, \eta)$ を適当に選ぶことにより格子点分布を制御する. この選び方については Steger and Sorensen[7] によって提案された方線を円柱表面に沿って描き,法を用いる.
線を 線に直交させる. 今回の計算で主に用いた格子は円柱に近づくほど小さく, 最小格子幅は円柱表面近くで 002,最大格子幅は 0.1 である. 格子点数は 方向に 310 点であり,方向には 229\sim 399
点の範囲で に応じて適当な数を選んぶ. このようにして作った計算格子の例を
図 2 に示す.
$\xi$$\eta$

$\xi$

$\xi$

$\eta$$\Gamma$

3.2 数値シミュレーションと定常解の数値計算法
数値計算においては, 計算領域を $L_{1}^{*}=5.5d^{*},$ $L_{2}^{*}=20.5d^{*},$ $L_{3}^{*}=9d^{*}+0.5s^{*}$とする. 数値シミュレーションでは, 空間微分項を 2 次精度の中心差分法により

7

図 3: 実験装置: 水槽, : ハイスピードビデオカメラ, : パソコン,
モーター.
ター, : レーザー, : ブラシ,$\mathrm{b}$$\mathrm{a}$$\mathrm{c}$$\mathrm{d}$: ブイル$\mathrm{f}$$\mathrm{g}:\mathrm{D}\mathrm{C}$$\mathrm{e}$差分近似し, 時間積分には 4 次精度の Runge-Kutta 法を用いる.

このとき, 時間とする.
刻み幅を同様に定常解の数値計算においても, 空間微分項を 2 次精度の中心差分法を用いて差分近似し, SOR 法を用いて計算を行う. SOR 法における緩和係数 vH よ,の計算において,軸に対して0.1\sim 15 の範囲で適当な値を選ぶ. 定常解対称な流れ場を求めるためには反対称条件$\triangle t=0.005$$(\overline{\psi},\overline{\omega})$$\overline{\psi}(x, y)=-/\overline{\psi}(x, -y)$

,$x$

(10)

$\overline{\omega}(x, y)=-\overline{\omega}(x, -y)$を課す 44. 実験装置および方法
4.1 実験装置
実験に用いた循環式水槽の概略図を図 3 に示す。 アクリル製の循環式水槽は水平に設置され, その中に水が満たされている. 水温は, 室温を調節することでほモーター (オリエンタ
ぼ 20 に保つ. 水槽内の流れは,ルモータ製, AXU590-GU) でロールブラシを回転することにより駆動し制御する.
流れをテストセクション断面内で一様にするため, 測定部の前後にステンレス製の格子フィルターを設置する.
であり,, 高さ 100, 長さ 450
テストセクション (図 3) は幅 200) を垂直に設置する. こ
その中央部に 2 本の円柱 (真鍮, 直径 2Inm, 長さ 100
れらの円柱の一端 (上端) をアクリル製の角棒に適当な間隔で固定する. 流れを特徴づけるパラメータとして, 間隔比(=s*/d りとレイノルズ数 $Re(=U^{*}d^{*}/\nu^{*})$
は円柱表面間の距離,は流入口での流は円柱直径,を定義する. ただし$\Gamma=1.0\sim5.0$ の範囲では水の動粘性係数である. 実験は $Re=30\sim 200$ ,速,$\zeta \mathrm{C}$$(\nu=1.010\mathrm{m}\mathrm{m}^{2}/\mathrm{s})$$\mathrm{D}\mathrm{C}$$[\mathrm{m}\mathrm{m}]$$[\mathrm{m}\mathrm{m}]$$[\mathrm{m}\mathrm{m}]$$\mathrm{m}\mathrm{m}$$\Gamma$$d^{*}$$\nu^{*}$

行う.

$s^{*}$$U^{*}$

$\epsilon$(b)$(_{\backslash }\mathrm{a})$図 4: 1 円柱の流れパターン. (a) $Re=40,$ $(\mathrm{b})Re=50$ .

4.2 PIV 測定
速度場を得るために PIV による測定を行う. YAG レーザー (Photonics Industries) を 0.5の厚-スイッチレーザー DMIO-527, 出力 15製,の高さの位置に水平に入射する. レーザーに照らされたさで水槽底から 60) を水槽の真下か, 直径 48\sim 52
トレーサー粒子 (オルガソル, 密度 103$\mathrm{m}\mathrm{J}/\mathrm{p}\mathrm{u}\mathrm{l}\mathrm{s}\mathrm{e}$$\mathrm{m}\ln$$\mathrm{N}\mathrm{d}:\mathrm{Y}\mathrm{L}\mathrm{F}\mathrm{Q}$$\mathrm{m}\mathrm{m}$$\mathrm{g}/\mathrm{c}\mathrm{c}$$\mu \mathrm{m}$らハイスピードビデ三郎メラ (Vision Research 製, Phantom V50, ピクセ) 数(0.0251632 1024 pixels) で撮影する. 測定領域の実寸法は 408mm$\mathrm{X}30.0\mathrm{m}\mathrm{m}$$\mathrm{x}$に移動する時間をパルス間隔
トレーサーが下流に 0.1とり, その 5 分の 1 をパルス幅とする. たとえば $Re=200$ のとき, パルス間隔は$\Delta t=1000\mu \mathrm{s}$ である.

$\mathrm{m}\mathrm{m}/\mathrm{p}\mathrm{i}\mathrm{x}\mathrm{e}\mathrm{l})$である.

$\Delta t$$\mathrm{m}\mathrm{m}$040020$v_{1}$$Re_{\mathrm{c}}=45.5[searrow]$

0-0.20-040
4055

45$Re$

図 5: 1 円柱を過ぎる流れの分岐図.

6095. 結果および考察
5.1 1 本円柱を過ぎる流れ
1 本円柱を過ぎる流れは低いレイノルズ数において対称定常流であり, ある臨界レイノルズ数 Re。を超えると振動流へ遷移する. ここで用いる数値計算法の正しさを確認するため, 1 本の円柱を過ぎる流れについて数値シミュレーションおよび線形安定性解析を行い,これまで報告されている結果と比較する.

レイノルズ数が $Re=40$ および $Re=50$ のときの 1 本円柱を過ぎる流れの流れ
でにおいて流れは対称定常流であるのに対して, 図場を図 4 に示す. 図は非対称な振動流となっている. 円柱から下流へ $1d$ 離れた点 (図 1) を代表点にをいろいろなレイノルズ数について計算した.
とり, この点における 方向速度
その結果を図 5 に示す. この図から, 流れは臨界レイノルズ数 $Re_{\mathrm{c}}=45.5$ を超えると定常対称流から振動流へと遷移し, そのときの分岐構造はホップ分岐であることがわかる. 対称定常流の線形安定性解析を行ったところ, 臨界レイノルズ数は $Re_{\mathrm{c}}=45.420$ であることがわかった. この値は数値シミュレーションで得られた値とよく一致をしている. また, Jackson[l] の安定性解析によれば, 臨界レイノルズ数の値は $Re_{c}=45.406$ であり, ここで得られた値とも良く一致している. これより, ここで用いる数値計算法は正しく, かつ精度良い結果を得ることができると結論づける.
$4(\mathrm{b})$$4(\mathrm{a})$$\mathrm{P}_{1}$$y$$v_{1}$

5.22 本円柱を過ぎる流れ
速度場測定

2 本円柱を過ぎる流れの速度場を, PIV 装置を用いて測定した,

測定は間隔比 $\Gamma=[1,5]$ の範囲において行い, 流れ場を分類した. 分類にあたって

各瞬間における速度場を数多く測定してそれらの平均速度場を対称定常速度場と仮定し, 各瞬聞の速度場から対称定常速度場を差し引いたものを撹乱速度場と呼ぶ. こうして得られた撹乱速度場の性質から流れ場の分類を行う. なお, ここで定義した撹乱速度場は対称定常流の線形安定性解析から得られる撹乱の速度場とほぼ同じであることが示される. PIV 測定により得られた瞬間速度場のベクはそれぞれレイノルズ数と間隔比がと
トル図は図 6 のようになった. 図は,$6(\mathrm{a})$$(Re, \Gamma)=(130.0,2.0)$$6(\mathrm{b})$と (135, 4.0) のときの速度場であり, 矢印はそれぞれの位置での速度ベクトルを表している. また色の濃淡は流れの速さを表しており, 色が黒
では円柱背後の 2 つの後流が同に近づくほど遅く, 白に近づくほど速い. 図
では逆位相で振動している流れのパターンとなって位相で振動しており, 図とそれぞれの場合の対称定常流を求めると, 図いることが確認できる.のような速度場が得られる. ここで, 対称定常速度場の計算には, 時間的に変動する流れ場よりランダムに 1000 の瞬間流れ場を取り出したものを平均した. それぞ$6(\mathrm{a})$$6(\mathrm{b})$$6(\mathrm{c})$$6(\mathrm{d})$れの場合について各瞬間の速度場から対称定常流を差し引いた撹乱速度場を計算

10(a)

(b)

(c)

(d)

(e)

(f)

図 6: PIV 測定により得られた速度場. (a), (b) 瞬間速度場. (c), (d) 平均速度場.
逆(e), (f) 撹乱速度場. (a), (c), (e) 同位相周期流 $Re=130,$ $\Gamma=2.0$ .
$(\mathrm{b}),$位相周期流. $Re=135,$ $\Gamma=4.0$ .

$(\mathrm{d}),$$(\mathrm{f})$11

図 7: 撹乱の速度場の分布図.

すると図$6(\mathrm{e})$と$6(\mathrm{f})$のようになる. 図$\nabla$$6(\mathrm{e})$: 同位相型,$\mathrm{v}$: 逆位相型.

では 2 本の円柱後方のベクトルの向き

ではその向きが逆 (逆位相) である. 得が同じ (同位相) であるのに対して, 図$(\Gamma,Re)$平面で整理すると図 7 のようになる. ここで,られた撹乱のベクトル場をはそれぞれ同位相と逆位相の撹乱が観測されたことを表す, この図から, 不と安定な撹乱のもつ空間的な構造は円柱間隔によって異なり, その構造は $\Gamma=2\sim 3$$6(\mathrm{f})$$\nabla$$\mathrm{V}$の間で入れ替わることがわかる.
対称定常流の線形安定性解析は Akinaga and Mizushima [8] によって既に行われでは最も不安定な撹乱の構造は同位相ており, その結果によれば,ではその構造は逆位相の流れパターンをもつことの流れパターンをもち,がわかっている. 今回得られた実験結果は Akinaga and Mizushima の結果と定性的に一致している. ここで, 強調しておく必要がある事柄は, $\Gamma=1$ のとき 2 つの振動パターンのどちらにも属さない複雑な流れを観測したことである. しかし,$\Gamma<1$ の範囲については円柱間隔が狭く実験が困難であるため,数値シミュレーションを行い詳しく調べる. その結果を次節で説明する.
数値シミュレ– ション 間隔比 $\Gamma=[0.3,0.9]$ の範囲について数値シミュレーションを行った. 出現する流れ場はレイノルズ数 $Re$ と間隔比 によって異なる.
これらの流れ場を図 8 の 4 つに分類した.
),
問隔比 $\Gamma=0.5$ の場合には, $Re=40$ のとき流れは対称な定常流であり (図). この$Re=50$ で流れが 軸に対して同位相で振動する振動流となった (図ような対称定常流から振動流への遷移について調べるために, 円柱から lld 下流とレイノルズ数 $Re$ との関係を調べた. この点 (図 1 戸こおける 方向速度をの最大値と最小値をレイの値を, 振動流の場合にはこで, 定常流の場合にはである.
ノルズ数の関数として分岐図を描いた. こうして得られた分岐図が図$\Gamma<\Gamma_{\mathrm{c}}=2.34$$\Gamma<\Gamma_{\mathrm{c}}$$\Gamma$$8(\mathrm{a})$$8(\mathrm{b})$$x$$\mathrm{P}_{2}$$y$$v_{2}$$v_{2}$$v_{2}$$9(\mathrm{a})$この分岐図より, 流れはある臨界レイノルズ数 Re。を超えると対称定常流から振動流へと遷移し, そのときの分岐構造はホップ分岐であることがわかる.

またそ12(a)

(b)

(c)

(d)

図 8: 流れパターン. (a) 対称定常流$Re=50,$$(\mathrm{c})$$\Gamma=0.5,$

振動偏向流. $\Gamma=0.5,$ $Re=60,$のときの臨界レイノルズ数は$Re_{\mathrm{c}}=4\ell \mathrm{d}.1$$(\mathrm{d})$$Re=45,$振動流$\Gamma=0.5$,非対称定常流. $\Gamma=0.6,$ $Re=55$ .

であった.

構造にもつ対称定常流から振動流の遷移は,$(\mathrm{b})$このようなホップ分岐を分岐$\Gamma=[0.3,0.5]$

および $\Gamma=[0.65,0.9]$ の範囲で生じることを確認した.
これまでよく知られていなかった流れが聞隔比 $\Gamma=0.5$ の場合に $Re=60$ において見いだされた. 流れは 2 本の円柱問で片方の円柱の方向に片寄り, 円柱の下流で振動している (図). このような流れはこれまで知られている力学系の基本的な分岐構造では理解することができない. この振動流から振動偏向流への遷移は,$8(\mathrm{c})$流れが対称定常流から非対称定常流へ遷移するピッチフォーク分岐とも関係していると考えられるので, 定常解探索を行い定常解の遷移について調べる. その結果を 2 本円柱下流の点まとめると, 図$9(\mathrm{b})$$\mathrm{P}_{1}$(図 1 戸こおける速度のようになる.

$v_{1}$とレイノルズ数 $Re$ との関係として

ここで点線は定常解探索により得た速度をプから, 対称定常流は臨界レイノルズ数 $Re_{\mathrm{c}}=54.5$ロットしたものである. 図を超えると非対称定常流へと遷移し, そのときの分岐構造はピッチホーク分岐であることがわかる. つまり $\Gamma=0.5$ の場合には, $Re_{\mathrm{c}}=45.6$ でホップ分岐が起き,さらにレイノルズ数を上げると, $Re_{\mathrm{c}}=54.5$ でピッチホーク分岐が生じている.
$9(\mathrm{b})$ここで述べた結果については説明が必要である. ホップ分岐が生じた後の流れになぜピッチフォーク分岐が起こる可能性があるのかという点については, (‘分岐の局所性” あるいは “不安定性の局所性” という新しい概念を導入しなければならない. すなわち, ある広い流れ場のある 1 点ではホップ分岐が生じ, そこから離13(a)(b)003

0OSC102006$e_{c}=43.1$$e_{c}=54.5$$.\ldots..-\cdot$.

$\ldots\ldots\ldots\ldots\ldots\ldots$$e_{c}-54.5[searrow]^{-}$$e_{c}-43.1[searrow]^{-}$$[searrow]$$[searrow]$$v_{1}^{0.02}0.040$$v_{2}0.010$

-0.01
-0.02
-0.03
40

45

50

-0.02
-0.04
-0.06
-0.08
50

$\ldots\ldots$

60

55

65

$\ldots\ldots\ldots..-\cdots\ldots$

70

$Re$

$Re$

図 9: 分岐図, $\Gamma=0.5$ . (a) ホップ分岐,
フォーク分岐,: 点 での 方向速度.
$\mathrm{P}_{1}$$v_{1}$.

$v_{2}$: 点 P2 での$y$

方向速度. (b) ピッチ$y$

れた別の点ではより大きなレイノルズ数でピッチフォーク分岐が生じると考えるのである. このとき, それぞれの分岐に対応する不安定モードはそれぞれの場所に局所化している. ただ, 今回の計算結果ではもう少し複雑であり, レイノルズ数が大きくなるに従って, 局所化していた撹乱が相互作用するようになる. このことを詳しく調べるために,$Re=55.0$ のときの流れを数値シミュレー). す
ションと定常解探索で求めると, 流れは非対称定常流となっている (図.
$=0.6$ ではホップ分岐が生じるレイノルズ数よりも小さいレイノルズ数
なわち,$\Gamma=0.6$ で$8(\mathrm{d})$$I^{\neg}$でピッチフォーク分岐が生じており, これら 2 つの分岐に対応する臨界レイノルズ数の大小関係が, $\Gamma=0.5$ と 06 の間で入れ替わるのである.
このようだして, 流れと垂直に並べられた 2 本の円柱を過ぎる流れの不安定性と分岐を $\Gamma=[0.3,5.0]$ の範囲で調べた結果をまとめると図 10 のような遷移図が得の結果も示されている. この図で,られた. この図には Akinaga and Mizushim$\mathrm{n}\mathrm{a}[8]$および・は, それぞれ同位相判の撹乱および逆位型の撹乱対して不安定となる臨
は数値シミュレーションにより得られた結果でと界点を表している. また,あり, それぞれピッチホーク分岐点とホップ分岐点を表している. 今回の研究で$\mathrm{o}$$\rangle\langle$$\triangle$は非対称定常流が見いだされたが, この流れが観測されるのは $\Gamma=[0.54,0.62]$ の極めて狭い領域である. ただし, そのときのピッチフォーク分岐の構造はもっと広い範囲まで及んでおり, その影響で偏流した振動流が観測されるのである.

参考文献[1] C. P. Jackson, A finite-element study of the onset of vortex shedding in flow
past variously shaped bodies”, J. Fluid Mech., 182 (1987) 23-45.
“

14: 同位二型撹乱に対する臨界レイノルズ数,型撹乱対する臨界レイノルズ数, , : Akinaga and Mizushima ([S]),フォーク分岐点, : ホップ分岐点.
図 10: ながれの遷移図.

$\bullet$$\circ$$0$$\bullet$$\mathrm{X}$: 同位置: ピッチ$\triangle$

The interaction between a pair of[2] P. W. Bearm an and A. J. Wadcock,astream”,
Fluid Mech., 61 (1973) 499-511.
J.
circular cylinders normal to“[3] S. J. Xu, Y. Zhou and R. M. C. So, Reynolds number effects on the flow
structure behind two side-by-side cylinders”, Phys. Fluids, 15 (2003) 12141219.
“Evolution of a single wake behind a pair of bluff[4] C. H. K. Williamson,bodies”, J. Fluil Mech., 159 (1985) 1-18.
“[5] I. Peschard and P. Le Gal,77 (1996) 3122-3125.

“

Coupled wakes of cylinder”, Phys. Rev. Lett.,[6] S. Kang, Characteristics of flow over two circular cylindrs in a side-by-sidearrangement at low Reynolds numbers’) , Phys. Fluids, 15 (2003) 2468-2498.
“

Automatic mesh-point clustering near
[7] L. J. Steger and L. R. Sorenson,
a boundary in grid generation with elliptic partial differential equation”,
J. Cornp. Phys., 33 (1979) 405-410.
“[8] T. Akinaga and M. Mizushima, Linear stability of flow past two cylinders
in a side-by-side arrangement” , to appear in J. Phys. Soc. $Jpn$ . (2005).
“