数理解析研究所講究録 1454 巻 2005 年 1-20

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

水島

二郎

末広

典久

(Jiro Mizushima)(Norihisa Suehiro)

1. はじめに
2 円柱を過ぎる流れは工学や土木などのさまざまな分野で調べられてきた. 2円柱の配置はタンデム型と並列型および斜め型の 3 つに分類できる. タンデム型
2 円柱を過ぎる流れの研究は双翼型飛行機の翼を支える柱への応用から始まっているため, 研究者の主な関心は円柱の抗力係数や揚力係数に向けられてきた. 現
在では, 流れと円柱の振動を抑制することを目的として, 海洋工学や土木工学にも応用されている. たとえば, 熱交換器中では高温の流体が流れている管の中に温度計が設置されている. 温度計は流れの中に振動を引き起こし, 流れから繰り返し振動する力を受けることで, 共鳴振動を生じて破断に至る. こうした破断を防ぐためにダミーの円柱を温度計のすぐ前に設置することで, 共鳴振動を抑制するという工夫が考えられている. 温度計の振動を最も効果的に抑制するためのダ

ミーの円柱と温度計との間の最適な距離が存在し得る.
3 つの典型的な配置の中で, 並列型 2 円柱を過ぎる流れは後流の相互作用を考えるための単純なモデルである. 後流の相互作用の結果として同位相と逆位相のは, 22 つの振動流のモードが観察されている. 実際, Peschard and Le
つの振動流のモードを表す連立ランダウ方程式を提案し, 実験で観測される振動流の 2 つのモードと定性的に良く似た結果を再現するモデル方程式を示した.
方, タンデム型 2 円柱を過ぎる流れは並列型配置の 2 円柱後流より複雑な現象を示す. なぜなら, 下流側の円柱は上流側円柱の後流の中に入り込み, 上流側円柱$\mathrm{G}\mathrm{a}\mathrm{I}[1]$によって引き起こされるせん断層や渦などを伴う複雑な流れが出現するからであ
この複雑性は, タンデム型 2 円柱についての研究論文は並列型のものよりかなり少ないという事実にも反映されている. タンデム型 2 円柱を過ぎる流れの最
る.

初のころの研究は主に実験であり, その概要は Zdravkovich[2] と Ohya, Okajima
and Hayashi[3] の 2 っのレビューにまとめられている.

これまでの研究では多くの実験が約 $10^{3}-10^{5}$ の高いレイノルズ数で行われおり[4, 5, 6, 7], 数値シミュレーションは低いレイノルズ数でなされてきた [8,9,10, 11].
タンデム型 2 円柱を過ぎる流れの実験と数値計算の結果で共通しているのは, 円
柱間距離を連続的に変えても, ある臨界距離を越えると流れのパターンが不連続的に変わることである. 2 円柱間隔が臨界値よりも小さなときには, 下流側の円柱は上流側の円柱の後流の中に入っており, 2 円柱の間にはほとんど流れがなく,
2 円柱は一つの大きな物体のようにみなされる. このことは, 下流側の円柱の存在が上流側の円柱からの渦の放出を抑制していることをも示している. しかし, 臨
界値よりも大きな間隔では 2 円柱それぞれの後流がお互いにほとんど独立してい2

るように振舞っている. その結果, 流れを特徴づける抗力係数と揚力係数およびストローハル数などの物理量は間隔比が連続的に増加あるいは減少しても, 臨界値で不連続な変化を示すのである.

この急激な変化の原因を解明するために多くの実験と数値計算がなされてきた.
臨界間隔の存在は Ishigaiet a1.[4] により, レイノルズ数が 1500–15000 の実験によって確認されている. 彼らは臨界間隔を, ストローハル数と表面圧力係数が急激な変化を示す $3.8d$ 付近であるとした, 高いレイノルズ数で行われたほとんどの実験において臨界間隔は $[3.5d-4.0d]$ の範囲にあることが示されている [5, 6, 7].
低いレイノルズ数における, 流れの物理量の不連続な変化についてはまだあまり多くの結果は示されていない. Tanida, Okajima and Watanbe[12] の実験結果では, $Re=80$ での抗力係数が間隔比の関数として, 間隔比 $3.9d$ 付近でつながる 2
つの異なる線分で描かれている. ただし, 彼らは実験結果から判断して, 臨界間隔比は約 $5d$ であると述べている,$Re=1\mathrm{O}\mathrm{O}$

で,

$3d$

から $5d$ の間の間隔で流れパ

ターンに明確な違いがあることが Ohmi and Imaichi[13] によって確認されている.
臨界間隔の存在は 100 から 200 程度の低いレイノルズ数での数値計算でもわかってる [8, 9, 10, 11]. これら $[3.75d-4.0d]$ の間の中でも Sharman et al.[ll] によ$1_{\mathit{1}}\backslash$

る最も最近の結果では臨界間隔が存在することが示されている.
Biermann and Herrnstein[14] は, 下流側円柱の抗力は上流側円柱の存在により減少するが, 上流側円柱の抗力は下流側円柱にはあまり影響されないことを示した. 下流側円柱の抗力が減少する大きさは間隔が $4d$ より大きければ間隔には依存
せず, $4d$ より小さければ下流側円柱の抗力は間隔の減少に伴い急激に減少し, さらに $3d$ より小さければ下流側円柱は前方への力を受け, 両円柱への正味の全抗力
は 1 円柱よりも少なくなることを示した.
タンデム型 2 円柱を過ぎる流れは, 多くの物理的に興味深い現象を含むさまざ
まな流れパターンを示している. 本研究においては, この流れの安定性と遷移について調べ, 特に, 抗力係数と揚力係数およびストローハル数について不連続な変化が起こる原因について詳しく調べる.

2. 基礎方程式
2.1 基礎方程式および境界条件
一様陸中に置かれたタンデム型 2 円柱を過ぎる流れを考える (図 1), 円柱の直径
を , 2 円柱の聞隔ををと定義する. 非圧縮 2 次元として, 間隔比流れ場を仮定し, 流れ関数 $\psi(x, y, t)$ と渦度 $\omega(x, y, t)$ を用いて, 渦度輸送方程式とポアソン方程式を無次元化した形で表すと$d$$\ell$$\Gamma$$\frac{\partial\omega}{\partial t}=$万$[perp]\equiv \mathrm{r}\backslash p/d$$( \psi, \omega)+\frac{1}{Re}\mathcal{M}\omega$$\omega=-\mathcal{M}\psi$

,

,

(1)
(2)

3

図 1: 2 円柱の配置と座標系.

,$\mathcal{M}=\frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}}$$N(f, g)= \frac{\partial f}{\partial x}\frac{\partial g}{\partial y}-\frac{\partial f}{\partial y}\frac{\partial g}{\partial x}$(3)を代表速度, 円柱の全ての変数は上流における一様流の流速
直径 を代表長さとして無次元化している. 時聞は $d/U$ によって無次元化されては動粘性係数である.
いる. レイノルズ数は $Re\equiv Ud/lJ$ で定義する. ただし,無限に広い流れ場を図 1 の ABCD で示される有限の領域で近似する. 一様流が上部境界 (AB) から流入すると仮定し, 流入条件を次のように表す.
となる.

ここで,

$U$

$d$

$lJ$

$\psi=y$

側面 (AD と$\mathrm{B}\mathrm{C}$,$\omega=0$$\psi=\psi_{2}$$\psi_{1}$

と

(4)

AB.

) でも一様流を仮定し, 次のようにおく.
$\omega=0$$\psi=\psi_{1\}}$

ただし,

on

$\psi_{2}$,$\omega=0$

onon$\mathrm{A}\mathrm{D}$$\mathrm{B}\mathrm{C}$

,

.

(5)

は定数である. 流出口 (CD) での流出条件はゾンマーフェルト

放射条件$\frac{\partial\psi}{\partial t}+c\frac{\partial\psi}{\partial x}=0$を用いる. ただし, 位相速度$c$,$\frac{\partial\omega}{\partial t}+c\frac{\partial\omega}{\partial x}=0$on$\mathrm{C}\mathrm{D}$(6)の与え方については後で説明する. 円柱表面では滑りなし条件が適用される. これは,$\psi=\psi_{\epsilon}$,$\omega=-\frac{\partial^{2}\psi}{\partial n^{2}}$(7)は円柱表面に垂直方向の座標である. これらの条件では,が与えられない限り, (1) と (2) の解を一意に定義するには不充分である. 円柱表を見積もるために, 流れ場全体で圧力を 1 価関数として面上での流れ関数の値のように表される.

$\psi_{S}$$n$$\psi_{\mathit{8}}$

条件を課す. これは,$\oint_{C}\nabla p\cdot tds=0$

(8)

4

はそれぞれの円柱を任意に囲む $C$ に沿った, 接線方は表面に滑向の座標と単位ベクトルを表す. ここでは, 円柱表面の圧力勾配りなし条件を適用することで次のように単純に表せるため, 円柱表面に沿った $C$と表される. ただし,

$s$

と

$t$

$\nabla p$を用いる.
$\nabla p=-\frac{1}{Re}\nabla \mathrm{x}\omega$(9).

2.2 定常解低いレイノルズ数においては, 流れは 2 円柱の中心を通る$x$軸に対して対称で定常である. 対称定常流に対応する解はレイノルズ数の大きさに関わらず境界条
件 (4) $-(7)$ のもとで基礎方程式 (1) と (2) を満たすが, この流れは臨界値を越えが線形安定性解析の主流でると不安定となる. この対称定常流, つまり,$(\overline{\psi}\mathrm{J}\overline{\omega})$

ある. 定常解は,次の定常渦度輸送方程式およびポアソン方程式を満たす.
$N( \psi,\overline{\omega})+\frac{1}{Re}\mathcal{M}\overline{\omega}=0$$\overline{\omega}=-\mathcal{M}\overline{\psi}$

(10)

,

(11)

.

方程式 (10) と (11) を適切な境界条件のもとに解くことで定常解は数値的に得られる. 対称定常解の境界条件は式 (6) の時間微分項がないことと円柱表面が 0 であることを除いて, 式 (4) $-(7)$ と同様である.
での$(\overline{\psi},\overline{\omega})$$\psi_{S}$2.3 定常解の線形安定性解析
定常解の安定性を解析するために, 定常解

渦度および流れ関数を式 (10) を引くことで,$\omega=\overline{\omega}+\omega’,$$\omega’$を付け加えて,に撹乱と表す. これらを式 (1) に代入し,$(\psi’,\omega’)$$(\overline{\psi},\overline{\omega})$$\psi=\overline{\psi}+\psi’$の撹乱方程式$\frac{\partial\omega’}{\partial t}=N(\psi’,\overline{\omega})+N(\overline{\psi}, \omega’)+N(\psi’\}\omega’)+\frac{1}{Re}\mathcal{M}\omega’$が得られる. 式 (12) において, 撹乱の非線形項
存性を指数関数と仮定して, $\psi’=\hat{\psi}(x, y)\exp(\lambda t),$$(\psi’)\omega’)$(12)を無視し, 撹乱の時間依と表すと,$\omega’=\hat{\omega}(x, y)\exp(\lambda t)$

線形化した撹乱方程式を得る.
$\lambda\hat{\omega}=N(\hat{\psi},\overline{\omega})+N(\overline{\psi},\hat{\omega})+\frac{1}{Re}\mathcal{M}\hat{\omega}$

.

(13)

この式を, 撹乱のポアソン方程式$\hat{\omega}=-\mathcal{M}\hat{\psi}$と連立して解く. ここで, 係数 は撹乱の複素線形増幅率である. その実部
虚部
は撹乱の増幅率ど角速度 (振動数) を表している.
$\lambda$$\lambda_{\mathrm{i}}$(14)$\lambda_{\mathrm{r}}$

と

5

撹乱

$(\hat{\psi})\hat{\omega})$の上流および側面の境界条件は,$\hat{\psi}=0$,

on$\hat{\omega}=0$(AB,$\mathrm{A}\mathrm{D},$$\mathrm{B}\mathrm{C}$

)

(15)

で与えられる. 下流の流出条件は,$\lambda\ovalbox{\tt\small REJECT}+c\frac{\partial\hat{\psi}}{\partial x}=0$$\lambda\hat{\omega}+c\frac{\partial\hat{\omega}}{\partial x}=0$on

CD(16)と表される. 円柱表面では滑りなし条件が適用される.
$\hat{\psi}=0$,$\hat{\omega}=\frac{\partial^{2}\hat{\psi}}{\partial n^{2}}$

.

(17)

式 (13) と式 (14) は境界条件 (15)-(17) のもとで, 固有値問題として解と表し, その固有関数
このとき固有値 は無限個あるが, これを固有値

ここで,

く.
を

$\lambda$$(\hat{\psi}_{k},\hat{\omega}_{k}),$$\lambda_{k}$$(k=1,2, \cdots)$と表して, それぞれをモードと呼ぶ. 固有値$\lambda_{h}$の実とと順序をつけて, 最大増幅モード
部の大きい順にならばとおく. このとき対称定常流はに着目し, これらを単にのときのレイノルズ数が主流の臨界レならば安定である.
不安定,$\lambda_{1}$$\lambda_{1\mathrm{r}}\geq\lambda_{2\mathrm{r}}\geq\cdots$$\lambda,$$(\hat{\psi}_{1},\hat{\omega}_{1})$$\lambda_{\mathrm{r}}>0$$(\hat{\psi},\hat{\omega})$$\lambda_{\mathrm{r}}=0$$\lambda_{\mathrm{r}}<0$

イノルズ数 Re。である.

2.4 非線形周期解
レイノルズ数が臨界値より大きくなると, 線形不安定性により, 撹乱は増幅す
るが充分時間が経過した後では, 非線形効果により, 平衡振幅をもつ振動流となをフーリエ級数展
る. 数値的に平衡状態の振動流を求めるために, 撹乱$(\psi’, \omega’)$

開する.
$\psi’=\sum_{m=-M}^{M}\tilde{\psi}_{m}\exp(\mathrm{i}m\Omega t)$,$\omega’=\sum_{m=-M}^{M}\tilde{\omega}_{m}\exp(\mathrm{i}m\Omega t)$

.

(18)

は平衡状態での撹乱の振動数であり, 打ち切り波数 $M$ は充分大きくとという関係が得られる.
, 飢の実数条件よりる. 撹乱の係数を等しいとお式 (18) を式 (12) に代入して, 同じフーリエ成分
ここで,$\Omega$$\tilde{\psi}_{-m}=\tilde{\psi}_{m}^{*}$$(\psi’, \omega’)$$m=\tilde{\omega}_{m}^{*}$$\exp(\mathrm{i}m\Omega t)$

くと,$(\tilde{\psi}_{m},\tilde{\omega}_{m})$の式$\mathrm{i}m\Omega\tilde{\omega}_{m}=N(\tilde{\psi}_{m)}\overline{\omega})$$+N( \overline{\psi},\tilde{\omega}_{m})+\sum_{p+q=m}N(\psi_{p}^{\tilde{/}},\tilde{\omega}_{q})+\frac{1}{Re}\mathcal{M}\tilde{\omega}_{m}$(19)を得る. ポアソン方程式 (11) より.
$\tilde{\omega}_{m}=-\mathcal{M}\tilde{\psi}_{m}$が導かれる.

(20)

$\mathrm{G}$

図 2: 計算格子の例. $\Gamma=2$ .

3. 数値計算法
3.1 格子形成
発展方程式の数値シミュレーションや定常方程式および線型安定性解析の計算など全ての数値計算において, 時間微分や空間微分を有限差分で近似する, 格子形成のために, 次のポアソン方程式による格子形成法を用いる [17].
$\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)$

.

(21)

は変換後の座標であり, 上式の右辺の関数 $P(\xi, \eta)$ と $Q(\xi, \eta)$ は円柱の表面付近で格子点が集中するように定める. ここで, $P(\xi, \eta),$ $Q(\xi, \eta)$ の形はここで,$(\xi, \eta)$

Steger and Sorenson によって提案された方法で定義する [18]. 座標

は円柱表面をに沿うようにとり,に直交するようにとる. また, $\eta=0$ は円柱表面であ$\Gamma=2$
る. 計算格子の例として,
における計算格子を図 2 に示す. この図では格$260$ 点であり, 最小格子幅は円柱表面付近で 002, 最大格子
子点は全部で 371$\eta$

$\xi$

$\xi$

$\rangle\langle$

幅は流出および上下境界で 02 である,

3.2 時間発展法
数値シミュレーションにおいて, 基礎方程式 (1) と (2) を境界条件 (4) $-(8)$ と適

切な初期条件のもとに時間発展法を用いて解く. 空間微分には 2 次精度の中心差
分を, 時間微分には 4 次精度の Runge-Kutta 法を用いた, ポアソン方程式 (2) の解法には SOR (Suecessive Over Relaxation) 法を用いた. 緩勅係数 については,
レイノルズ数や間隔比によって適切に決めるが, 多くの場合 $\epsilon=1.5$ とした. ポ
アソン方程式の収束条件としては,について, 次のステップとの相対誤差が計算領域の全格子点で 10 以下になったとき十分収束したとみなした. 隣接した 2 っの時問ステップ聞のおよびの相対誤差の合計がそれぞれ $10^{-10}$ より小さくなったとき, 流れが定常状態であると判断した.
$\epsilon$$(\psi, \omega)$$\psi$$\omega$7

3.3 SOR 法
定常流 (主流) を得るためと, 線形安定性解析の固有値問題を解くのには SOR法を用いた. 式 (10) と (11) を境界条件 (4) $-(7)$ のもとで数値的に解く . 空間微と $Re$ の値によって,分を 2 次精度の差分で近似し, SOR 法の緩和係数 は$0.3<\epsilon<0.6$ の範囲で定める, 計算時間を節約するために
が
軸に対
および
して反対称であるという条件を用いる. すなわち,として定常解の計算を行った.
式 (13) と (14) および境界条件 (15)–(17) からなる撹乱の固有値問題を定常解の$\Gamma$$\epsilon$$(\overline{\psi},\overline{\omega})$$x$$\overline{\psi}(x, -y)=-\overline{\psi}(x, y)$$\overline{\omega}(x, -y)=-\overline{\omega}(x, y)$も反復によって求める従属変数
数値計算と同様の手順で数値的に解く. ここで,であるとみなす. 解の収束については, 定常解の判定法と同様の方法で判定した 8最大増幅モードの固有関数は対称性 $\hat{\psi}(x, -y)=\hat{\psi}(x, y)$ およびがあるので, この対称性を数値計算に用いた.
非線形周期解は適切な初期条件があれば SOR 法を用いて計算できる. 非線形周
を非線形周期解を得るための従属変数とする.
期解の数値計算において, 振動数とは異なるた
なぜなら,は非線形効果により線形安定性解析で見積もられためである. 計算時間と容量を節約するために, 偶数次のフーリエモードには反対
称条件を, 奇数次のフーリエモードには対称条件を課す. すなわち,$\lambda$$\hat{\omega}(x, -y)=\hat{\omega}(x, y)$$\Omega$$\Omega$$\omega_{\mathrm{i}}$$\tilde{\psi}_{m}(x, -y)=$$(-1)^{m+1}\tilde{\psi}_{m}(x, y)$

および $\tilde{\omega}_{m}(x, -y)=(-1)^{m+1}\tilde{\omega}_{m}(xy,)$ とした.

4. 計算結果
4.1 計算方法の精度見積もりと有用性
計算結果を紹介する前に, 今回用いた計算方法の精度について見積りを行う. 全ての数値計算において, 計算領域を $L_{1}=5d,$ $L_{2}=20d$ および $W=8d$ とした.
および $W$ を順に変えて計算を行い,計算領域の大きさを検定するために,と揚力
ストローハル数 $St$ および上流側円柱と下流側円柱の抗力係数 $C_{D1},$$L_{1},$$L_{2}$$C_{D2}$

係数$C_{L1},$$C_{L2}$の数値誤差を見積った. このようにして得られた $\Gamma=3,$$Re=1\mathrm{O}\mathrm{O}$における数値シミュレーションの計算結果を表 1 に示す. $L_{1}=5d$ と $7.5d$ を用いの絶対値は非常て計算したストローハル数 $St$ の相違は 07% 以内であった.
の誤差は 45% 程度となった. 計算領域を $L_{1}=5d,$ $L_{2}=20d$に小さいので,および $W=8d$ としたとき, 計算結果に含まれる数値誤差は 45%以内であると判断し, この計算領域で全ての数値計算を行う.
で 4 つの異なる格子点数を選んで,格子点数を決めるために $\Gamma=3.0,$
数値シミュレーションを行った. この結果を表 2 に示す. この表より数値誤差は1,8% 以内であることがわかる. よって, 数値シミュレーションを含む数値計算で$C_{D2}$$C_{D2}$$Re=1\mathrm{O}\mathrm{O}$は $380\cross 200$ の格子点を用いる.

8

表 1: 計算領域の大きさの検証.
$St$(a) 上流長さ$C_{D1}$$C_{D2}$$C_{L2}$$C_{L1}$$L_{1}$$(L_{2}=20d, W=8d)$0.12191.181 -0.0421 0.1951 0.0643
0.12281.177 -0.0441 0.1905 0.0633$L_{1}=5d$$L_{1}=7.5d$(b) 下流長さ$L_{2}$$(L_{1}=5d, W=8d)$$L_{2}=20d$$L_{2}=25d$0.1219 1.181 -0.0421 0.1951 0.0643
0.1226 1.179 -0.0433 0.1889 $0.0639$(c) 流域幅 $W$$(L_{1}=5d, L_{2}=20d)$$W=8d$$W=1\mathrm{O}d$0.1219 1.181 -0.0421 $0.1951$ 0.0643
0.1226 1.178 -0.04310.1885 $0.0638$

計算結果の誤差を見積もり, 計算方法の正しさを確認するために, 一様流中に置いた 1 円柱を過ぎる流れの臨界レイノルズ数を求めた. 1 円柱を過ぎる流れは,本研究では $\Gamma=0$ の場合に相当する. この流れは安定性も含め理論的および実験的に詳細に調べられてきた最も単純な流れである. 低いレイノルズ数において, 1

円柱を過ぎる流れはどのような初期条件であっても十分時聞が経過すると定常状態となる. 定常状態では流れ場は円柱の中心を通る流線に対して対称である. 数
値シミュレーションによって得られた分岐図を解析した結果,

臨界レイノルズ数
は $Re_{\text{。}}=46.731$ であった, また, 定常流の線形安定性解析から得られた臨界値は$Re_{\text{。}}=46.088$ であった.
数値シミュレーションと線形安定性解析における Re。の違いは約 1.4% であり, どちらの方法を用いても正しい結果が得られている.
Mathis, Provansal and Boyer[19] が実験によって見積もった臨界レイノルズ数は$Re_{c}=47.\mathrm{O}$ であり,Jackson [20] が線形安定性解析によって見積もった理論値は$Re_{\text{。}}=46.184$ であった.
これらの値と今回計算で得られた値の違いは 2% 以内である.

4.2 流れパターン
抗力・揚力. ストローハル数などの物理量の不連続な変化は間隔比 $\Gamma=.3.5-3.9$で起こっていると報告されているので, 本研究では, $\Gamma=[1,6]$ の範囲で, さまざ
まな値のについて計算を行った. ここでは, 典型的なケースとして $\Gamma=2$ と 4の 2 つの間隔比について, その不安定性と遷移について説明する. 間隔比の値にかかわらず, 低いレイノルズ数においては, 円柱の中心を通り 軸に対して対$\Gamma$$\Gamma$$x$

9

表 2: 計算格子戸の検証.

Grid number$380$ $\mathrm{x}200$

760
380
760

$\rangle \mathrm{e}200$$\mathrm{x}400$$\mathrm{x}400$$St$

0.12190.122401224
0.1227$C_{D2}$$C_{D1}$1.1810-0.0421
1.1842 -0.0415
1.1831-0.0419
1.1865-0.0411$C_{L1}$$C_{L2}$

0.19590.1986
0.1966
0.19940.064000632
0.06350.0631(b)(a)

図 3: 流れ場 (流線). 定常対称流 $Re=60$ .

$(\mathrm{a})\Gamma=2$.

$(\mathrm{b})\Gamma=4$.

称で定常な流れとなる. このような対称定常流は定常方程式 (10) と (11) を解くことで得られ, その線形安定性を数値的に調べる. $Re=60$ における $\Gamma=2$ と 4 のとに示す. 下流側円柱直後の双子渦の長さは
流れパターンをそれぞれ図
どちらの流れパターンもほとんど同$\Gamma=4$ の方が $\Gamma=2$ より短いことを除けば,$3(\mathrm{a})$様である. また,

$3(\mathrm{b})$上流側円柱からはがれた流れが下流側円柱に再付着し閉じた流線が両方の場合に見られる. これらは,

2 円柱聞の閉じた流線を外側輪郭とする表面をもつ一つの物体を過ぎる流れであるかのように見える. $\Gamma=4$ は非線形性を考慮に入れると不安定となるが, どちらの流れも線形安定性解析によって安定であることが示されている.
$\Gamma=2$ と 4 のどちらの場合においても,
流れは(a)$Re=1\mathrm{O}\mathrm{O}$で対称性を失ない, 振(b)

図 4: 流れ場 (流線). 振動流 $Re=100$ .

$(\mathrm{a})\Gamma=2$.

$(\mathrm{b})\Gamma=4$.

$\mathrm{J}0$

(b)

(a)

$\mathrm{x}10^{\sim 2}$$\mathrm{x}\mathrm{I}\mathrm{O}^{-2}$

60

1.5

4.0

1.0

2. 0

0.5

夏

$’<^{\mathrm{h}}$0$\nwarrow Re_{c}=78.5$

-0.5
-1.0

0
-2.0

$\nwarrow Re_{c}=68.0$

-4.0
-6.0

85

80

75

66

68

70

72

$Re$
$Re$

図 5: 線形増幅率$\lambda_{\mathrm{r}}$.

$(\mathrm{a})\Gamma=2$.

$(\mathrm{b})\Gamma=4$.

十分時間が経過したあとで平衡状態になることが数値計算によって明に示す. どちらとらかになった. $\Gamma=2$ と 4 の瞬間の流れパターンを図の場合も 2 円柱後流で振動流が見られ, その振幅については $\Gamma=4$ の方が $\Gamma=2$の振幅よりも大きい. $\Gamma=2$ のとき, 2 円柱の間の領域では振動していないことがわかる (図). 上流側円柱からはがれた流線が下流側円柱で再付着し, 2 円

動し始め,$4(\mathrm{a})$$4(\mathrm{b})$$4(\mathrm{a})$で流れは閉じた流線
柱を囲う閉じた流線が存在している. それゆえに,の形をもつ 1 つの物体を過ぎる流れであるかのようにみなされる. 一方, 図 4(b)$Re=1\mathrm{O}\mathrm{O}$で示される $\Gamma=4$ については 2 円柱の間では流れが振動流となっていることが観察できる. この $\Gamma=2$ と $\Gamma=4$ の流れパターンの違いが臨界間隔付近での急激な物理量の変化を示すことを表している.

4.3 定常流の線形安定性解析
数値シミュレーションを行った結果,
れは $Re=60$ で対称定常流$\Gamma=2$$Re=1\mathrm{O}\mathrm{O}$と 4 のどちらの場合についても, 流で振動流が現れることを確認した. これより, 対称定常流はその安定性を臨界レイノルズ数で失っていることになる. この流れの線形安定性を調べるために, 数値的に方程式 (10) と (11) を解くことで,を計算し, 方程式 (13) と (14) からなる固有値問題を境界条件
対称定常解(15) -(17) のもとで解く.
とに図示する. 計算により線を図
固有値の実部である線形増幅率
は $\Gamma=4$ のときの結果で, これら
は $\Gamma=2$ , 図
形増幅率 を求めた. 図$\Gamma=4$ では $Re_{c}=68.0$の図より臨界レイノルズ数 Re。は, $\Gamma=2$ でが得られた. これらの臨界レイノルズ数より判断すると, 流れは無限小撹乱では$\Gamma=2$ の方が $\Gamma=4$ よりも安定であることがわかる.
これより, 1 円柱を過ぎる
流れの臨界レイノルズ数は $Re_{\text{。}}\sim 46-47$ であることを考慮すると, ある物体の近くにもう一つ別の物体が置かれているとき, 流れは安定化するといえる.
$(\overline{\psi},\overline{\omega})$$\lambda_{\mathrm{r}}$$\lambda$$5(\mathrm{a})$$5(\mathrm{a})$$5(\mathrm{b})$$5(\mathrm{b})$$Re_{\mathrm{c}}=78.5,$

11

(b)

(a)

図 6: 流れ場 (流線). 最も不安定な撹乱の野曝$Re=6\mathrm{S}.0$$\hat{\psi}_{\mathrm{r}}$.

$(\mathrm{a})1^{\neg}=2,$$Re=78.2$ .

$(\mathrm{b})\Gamma=4$,.

は撹乱の角速度であり, 流れにおける振動の無次元振動数であるスト
二部のとのように計算できる. $St$ の値は図
ローハル数 $St$ はレイノルズ数の範囲ではほぼ一定で, $\Gamma=2$ で 0.122, $\Gamma=4$ では 0.160 であり,$\lambda_{\mathrm{i}}$$5(\mathrm{a})$$St=\lambda_{\mathrm{i}}/(2\pi)$$5(\mathrm{b})$この値は $Re$ にほとんど依存しない 3の実部は, $\Gamma=2$ と 4 のそれぞれの
撹乱の流れパターン, つまり, 固有関数のようになる. ここで, 固有関数 は下流直円柱から約と臨界値で図として正規化している.
において$6.3d$ 後方, 中心線 $(y=0)$ 上の代表点
ただし, 採用した点と下流側円柱の距離は正確に $6.3d$ ではない. なぜなら, 本研$\hat{\psi}_{\mathrm{r}}$$\hat{\psi}$$6(\mathrm{a})$$6(\mathrm{b})$$\hat{\psi}=1$$\mathrm{P}_{1}$

究で用いた計算格子は代表点の位置での正確な格子点をもっておらず, 代表点にで2 円
最も近い点を採用しているからである. 2 っの流れパターンの図は図柱の間に見られる閉じた流線を除き大きな違いはない. しかし, 固有関数のわず$6(\mathrm{b})$かな違いが流れの物理量に起こる急激な変化の説明に対する重大な役割を果たしている. 固有関数の乙部の流れパターンはこれらの図とほぼ同様なので省略する.

4.4 分岐図
対称定常流 (主流) が不安定になると, 振動流を生じる. 流れの遷移を調べるた
めに, 数値シミュレーションおよび周期解の直接的な数値計算の結果を用いて, 流
れの分岐図を描いた, 流れの中の振動の大きさを表す代表的な物理量として, 下の振動振幅を用いる.
流側円柱から約 $6.3d$ 後方で 方向速度である. これらの図では, 流れが振動
こうして得られた分岐図が図の振幅 (最大値, 最小値) が $Re$ の関数として描かれている. 分岐のするときの$v_{1}$$y$$7(\mathrm{a})-7(\mathrm{e})$$v_{1}$

形状の変化を詳細に追うために $\Gamma=2$ と 4 に加えて, $\Gamma=3,3.1$ および 33 の間隔比を採用した. $\Gamma=2$ (図 7(a)) では, 分岐図はホップ分岐を構成しており, 1
円柱を過ぎる流れの分岐図とよく似ている. これは, 流れが振動して 1 るときに2 円柱が流れ場の中で単独の障害物と見なされるという 42 節の結果を裏づける) は $Re<Re_{c}=78.2$ (AB) の範囲で安定,ものである. 対称定常流 (この図で$\mathrm{A}\mathrm{C}$

12

(b)

(a)

o.os

0.3
$\mathrm{D}$$\mathrm{D}$

0.2

0.04

$\mathrm{T}$

0.1
A

\^a

———————–$\mathrm{c}$0$\underline{\mathrm{B}}_{---------------------------}\mathrm{c}$$\mathrm{A}$$\mathrm{B}\nwarrow_{Re_{L}=78.2}--------------------$

-0.1
-0.2

$\mathrm{T}’$$\mathrm{D}’$$\mathrm{D}’$$-0.0870$$-0.3$

100

90

80

90

$\mathrm{S}0$

$Re$

120

110

100
$Re$

(d)

(c)

06

06

$\mathrm{D}$$\mathrm{E}$T2$\mathrm{D}$$\mathrm{T}_{2}$

0.4

0.4

$^{\sim_{------\sim}}.j_{\mathrm{T}_{1}}$$r\mathrm{r}$$\sim\sim\sim-----------^{\mathrm{T}}----\sim$$.\mathrm{F}\sim\sim---_{--\sim}j\mathrm{T}_{1}$0202$—–\mathrm{F}_{\dot{i}}.\underline{\swarrow^{Re_{t2}=86}}\dot{} jRe_{t1}=11\mathrm{a}_{}j$$R\dot{\text{\’{e}}_{1}}\mathrm{F},=11_{\delta}^{r}[searrow]^{}i$$\mathrm{B}_{----------arrow-}^{}R\dot{\text{\’{e}},1}=11_{\delta}^{r}\underline{[searrow]}i$$\tilde{\mathrm{p}}$$\mathrm{c}$$\mathrm{A}$0$\mathrm{B}_{-----------*-,\dot{}}.\cdot\backslash$$0$$\tilde{\mathrm{i}\supset}$$-arrow_{Re_{\mathrm{c}}=83.3\overline{\nwarrow}\dot{_{\grave{\dot{}}}}}^{\overline{}}---------=-$$arrow Re_{\mathrm{c}}=83.\dot{3}\nwarrow$

:

-0.2

$\dot{}$-0.2$Ret\mathit{2}=10_{\sigma}0\dot{}$$Re_{t2}=10_{\sigma}0\dot{}$$-\yen---\neg\dot{}\mathrm{T}_{\acute{2}1}\mathrm{F}’\dot{}3_{\acute{}}^{i}\tilde{^{-------- \mathrm{R}--_{\overline{\mathrm{T}}}}\mathrm{V}}He_{c}=\mathrm{S}0.7j\underline{\mathrm{q}}_{\underline{2}}----*./\dot{^{}}$$_{_{\wedge}-\sim--\sim--\mathrm{T}_{\acute{1}}}^{}.\mathrm{F}’\dot{^{}}$$_{_{\wedge-\sim-----}}^{}.\mathrm{F}’\dot{^{}}\mathrm{T}_{1}’$

-0.4

-0.4

$\mathrm{T}’$$\mathrm{T}_{2}’$$[perp]_{2}$$\mathrm{D}’$$\mathrm{D}’$$\mathrm{E}’$

-0.6

70

80

100

90

110

-0.6

120

90

80

100

110

$Re$
$Re$

(e)
06
$\mathrm{E}$

0.4

T2

$\mathrm{D}$$\mathrm{T}_{2}$$i\mathrm{E}$.

$Re_{\{}=\mathrm{S}8_{\delta}^{r}$$\mathrm{A}j\.H^{Re_{\{}=\mathrm{S}8.\mathit{0}}^{\backslash }i\backslash \backslash \backslash _{\backslash \underline{\mathrm{B}\mathrm{C}}}ji.r$$jneq^{\backslash }i\backslash \backslash \backslash _{\backslash \underline{\mathrm{B}\mathrm{C}}}\dot{j}$$0.20$$\mathrm{A}$歌$j\mathrm{G}$-0.2$/\dot{i}//$$\prime i\nearrow/!.\nwarrow_{Re_{c}=78.2}$$i_{}$$\mathrm{T}_{2}^{l\dot{}}_{}i_{/}/\nearrow i_{\mathrm{G}^{J!}}/’i_{}_{17’}.\overline{.\nwarrow}Re_{c}=78.2$.

$\dot{}$$-0.4$$\mathrm{T}_{2}^{l}$.

$\mathrm{D}’$$\mathrm{E}’$

-0.6

80

60

100

$Re$

図 7: 分岐図. (a) $\Gamma=2.0$ .

$(\mathrm{b})\Gamma=3.0$.

$(\mathrm{c})\Gamma=3.1$. $(\mathrm{d})\Gamma=3.3$ . $(\mathrm{e})\Gamma=4.0$ .

13では不安定である. Re=Re。で超臨界ホップ分岐が起こり, 分岐
では振動流でのの最大値と最小値を表している.
した枝の BD とにおける点) が現れ間隔比が 3 まで大きくなると分岐図に変曲点 (図る. この変曲点は $\Gamma=3.1$ では 2 つの折りたたみ構造になり, その転換点は図のと T2 で示されている. このときには, 3 重解が存在する. つまり,
この
$Re_{t2}\leq Re\leq Re_{t1}$ の範囲において 2 つの安定な解と 1 っの不安定解があり,
結果, ヒステリシスが生じる. すなわち, レイノルズ数が増加するとき, A からかを通る. 一方, レイノルズ数が減少するとき,への解の道筋は$Re>Re_{c}(\mathrm{B}\mathrm{C})$$\mathrm{B}\mathrm{D}’$$v_{1}$$\mathrm{T}$$7(\mathrm{b})$$7(\mathrm{c})$$\mathrm{T}_{1}$$\mathrm{D}$$\mathrm{B},$ら A への道筋は$\mathrm{T}_{1},$$\mathrm{D}$$\mathrm{E}$を通る. 以上から,の範囲において, レイノルズ数が変化する向きに応じて 2 つの異なる解が生じることが明らかにされた.
$\mathrm{T}_{2},$$\mathrm{F},$$\mathrm{B}$$Re_{t2}\leq Re\leq Re_{t1}$

(b)

(a)

(c)

. $\Gamma=3.3$ . (a): 定常状態 (図図 8: 多重解.
態(点), (c): 定常状態 (Q3 点).
$Re=1\mathrm{O}\mathrm{O}$$7(\mathrm{d})$の$\mathrm{Q}_{1}$点), (b):不安定状$\mathrm{Q}_{2}$に示されている $\Gamma=3.3$ の分岐図は $\Gamma=3.1$ の分岐図と同様の形であのる. 唯一の違いはホップ分岐と転換点のレイノルズ数 Re。と Retl およびとの値は Re。よりも小さ値だけである. 間隔比がより大きくなると,と合体する. この結果 亜は $\Gamma=4.0$ でホップ分岐点とくなり転換点図$7(\mathrm{d})$$Re_{t2}$$Re_{t1}$$Re_{t2}$$\mathrm{B}$$\mathrm{T}_{1}$$\mathrm{T}_{1}’$). この分岐図でもヒステリシスが存在する. レ
臨界ホップ分岐が生じる (図への解の道筋は から へとジャンプし,イノルズ数が増加するとき, A から
へとジャン
から
へは
から
他方, レイノルズ数が減少する場合には,$7(\mathrm{e})$$\mathrm{E}$$\mathrm{B}$$\mathrm{D}$$\mathrm{D}$$\mathrm{A}$$\mathrm{G}$$\mathrm{T}_{2}$プする.

2 っの安定解と 1 っの不安定解が聞丁零 $\Gamma=3.1$ と 3.3 で,$Ret2\leq Re\leq Re_{t1}$

14

120

Periodic flow

$1\backslash b_{\iota}$

100$\backslash \backslash \backslash \backslash \backslash \backslash$80

避

2L‘\sim へ\sim \--60

オ

Steady flow

20

01

3

2

5

4

6

$\Gamma$

図 9: 遷移図.

$0$: ボツプ分岐点$(Re_{\text{。}}),$$\triangle$: 転換点$(Re_{t})$. $arrow:1$ 円柱の場合のホップ

分岐点 4の範囲に存在する.

$\Gamma=3.3$

図 8 に示す ただし, 図

での,$Re=1\mathrm{O}\mathrm{O}$とおよび$8(\mathrm{a})$$8(\mathrm{b})$のときの 3 つの解の流れパターンを$8(\mathrm{c})$は$\mathrm{Q}_{1}$と Q2 および$\mathrm{Q}_{3}$で示される解である. いうまでもなく, 流れパターンだけではこれらの安定性をはっきり
と
に描かれ
より, 図
見分けることはできない. しかしながら, 図ている振動流の振動振幅が最大のものと最小のものが安定であり, 中程度の振幅(図$8(\mathrm{b})$$8(\mathrm{c})$$8(\mathrm{a})$$7(\mathrm{d})$) をもつ振動流が不安定であることがわかる.

4.5 遷移図で亜臨界ホップ分岐が生じ,

この 2つの間隔比の間では分岐図の形は複雑に変化することが前節で示された.

ここで間隔比 $\Gamma=2$ で超臨界ホップ分岐,
は,$1\leq\Gamma\leq 6$$\Gamma=4$の範囲で分岐図を求めて, 超臨界ホップ分岐点, 亜臨界ホップ分と岐点 Re。と転換点を評価した. これを図 9 に示す. この図において
によって示されている, 矢印はホップ分岐の臨界レイノルズ数は , 転換点は
1 円柱を過ぎる流れのホップ分岐による臨界レイノルズ数 $Re_{\text{。}}=46.7$ を表す. この図より, 分岐は $\Gamma<3.4$ あるいは $\Gamma>5.5$ で超臨界ホップ分岐, $3.4<\Gamma<5.5$で亜臨界ホップ分岐である. また, 1 組の転換点が $\Gamma\sim 3.05$ で現れ, レイノルズ数が小さい側の転換点は $3.05<\Gamma<5.5$ の範囲で存在する. 一方, 大きい側の転で消失する.
換点はホップ分岐が超臨界から亜臨界へと変化する
狭い間隔比では, 分岐図は 1 円柱を過ぎる流れのものと似ている. このことは, 間隔比がの範囲では臨界レイノルズ数 Re。は の値に依存するが, 流れパターンは 2 円柱が 1 つの結合した物体であるかのように考えられるというこれまでの考察を支持するものである. 一方で, $\Gamma>5.5$ の広い問隔比については 2 つの後流はお互いにほとんど独立しており, 臨界レイノルズ数は 1 円柱の値に近づいている. これらの間の間隔比では, 上流側円柱からはがれた流れが$Re_{t1}$$Re_{t2}$$\circ$$\triangle$$\Gamma\sim 3.8$$\Gamma\leq 3.1$$\Gamma<3.1\sim$$\Gamma$15

し$\Gamma$

図 10: 平均抗力係数$C_{D}$$C_{D}$.

$\circ$: 上流側円柱の$C_{D}$.

$\bullet$;

下流側円柱の$C_{D}$. $arrow:1$ 円柱の.

下流側円柱に対し, 非線形効果によって分岐図は転換点で折り返し, 結果的に解が多数生じるという複雑な相互作用をしていることになる.

4.6 揚力と抗力
本研究における最大の関心事は抗力係数と揚力係数における不連続な変化である. これらの物理量はこれまで行われてきたほとんどの研究における共通の関心事であり, 広範囲にわたって調べられており,$\Gamma$

が連続して変化すでこれらの物理量の急激な変化が起こることがわかって
本節では, 実験結果と比較するために, 1 つのレイノルズ数の値について

るとき,
いる,

これにより,$\Gamma=3.5\sim 3.9$$(Re=1\mathrm{O}\mathrm{O})$,

さまざまな値の$\Gamma$で,計算によって両方の円柱の抗力係数

数値シミュレーションと非線形周期解の数値$C_{D}$と揚力係数$C_{L}$を求めた.

流れによって物体に作用する力は物体周りの圧力分布および流れと物体表面の摩擦によって生じている. ここで, 物体に働く力の流れ方向の成分は抗力 $D$ , 垂直な成分を揚力 である. 抗力係数と揚力係数は無次元化し, 次の式で求められる.
$L$,$C_{D}=2D=-2 \oint p_{b}\frac{\partial x}{\partial\xi}d\xi+\frac{2}{Re}\oint\omega_{b}\frac{\partial y}{\partial\xi}d\xi$$C_{L}=2L=-2 \oint p_{b}\frac{\partial y}{\partial\xi}d\xi+\frac{2}{Re}\oint\omega_{b}\frac{\partial x}{\partial\xi}d\xi$

ただし,$p_{b}$と$\omega_{b}$は物体表面の圧力と渦度であり,

$\xi$

.

(22)
(23)

は物体表面に沿った座標である.

抗力は流れが振動しているときには時間的に変化するので, 1 周期分の時間で平均化した値を求めた. 抗力係数は間隔比が変わると図 10 のように変化する. ここと( ) はそれぞれ上流側円柱と下流側円柱に作用する平均抗力
で,$C_{D1}(\circ)$$C_{D2}$$\bullet$1Bの値はおよそ 12\sim 14 で, 最も信頼で係数である. 一般に 1 円柱の場合,
は
きる値は 1.33[11] である. これは図 10 に矢印で示されている. 抗力係数
これらの値は, ほとんど$\Gamma<3.5$ でおよそ 12, $3.2\leq\Gamma\leq 6.0$ では 1.34 である.
$C_{D}$$C_{D1}$の値より判断すると, 実際には上流側円柱は下
1 円柱の場合と同じである.
の流側円柱の影響を受けていないように見えるが, 図 10 を注意深く見るととの領域でで 2 つの転換点が見られる.
と 35
線上にの 3 重関数である. このことは, この範囲において解が多数存在しは,は$C_{D1}$$C_{D1}$$\Gamma\sim 3.1(\mathrm{D}_{1})$$C_{D1}$$\Gamma$ていることを示している.

の値は$Re=1\mathrm{O}\mathrm{O}$が増加するときへはジャンプをする. 他方,からと至る. このときには,$C_{D1}$$\mathrm{C}_{1}$$\mathrm{D}_{1}$$(\mathrm{C}_{1})$$\Gamma$において,$\mathrm{A}_{1}\mathrm{B}_{1}\mathrm{C}_{1}$$\mathrm{E}_{1}$$\Gamma$$\mathrm{D}_{1}$から$C_{D1}$$\mathrm{E}_{1}\mathrm{F}_{1}$の測定を実験により行うと,

へと至る.

このとき,

は
が減少するときへとジャンプする.
にも現れている. 抗力係数
$C_{D1}$$\mathrm{F}_{1}\mathrm{C}_{1}\mathrm{D}_{1}$から$\mathrm{C}_{1}$から$\mathrm{B}_{1}\mathrm{A}_{1}$へ$\mathrm{B}_{1}$よりは多重解の存在は下流側円柱ので下流側円柱に及ぼされる力は抗力ではなく推も大きく変化をしている.
の値は $3.1<\Gamma<3.5$ で 3 つのが増加するにつれて,進力である. 間隔比
で 1.1 まで増加する. $\Gamma>4.0$ では $C_{D2}\sim 1.0-1.1$ となりほ
値をもつが,
ぼ一定である. このことは, 2 円柱の後流がお互いにほとんど独立していることが負になることは, 下流側円柱が上流側円柱の後流の証明である. 一方で,に完全に包まれ,
では上流側円柱から大きく影響を受けていることを表しが 3 つの値をもっている. 2 つの後流の非線形相互作用は $3.1<\Gamma<3.5$ でことを見ても大変複雑になっている. ここでの結論としてはこれまでの研究結果と矛盾していない.
も臨界レイノルズ数よりも大きいレイノルズ数では周期円柱に作用する揚力的に変化する. しかし, 計算領域の対称性から判断して, 1 周期の時間で平均化した値は 0 となる. よって, ここでは, 上流側円柱と下流側円柱における揚力係
数の振動振幅( ) を図 11 に $1.5\leq\Gamma\leq 6.0$ の範囲で示す. い
うまでもなく,の両方の線は $3.1<\Gamma<3.5$ で多重解の存在を示す 3 重解となる.
での 1 円柱における揚力の振幅は $C_{L}\sim 0.08$ と報告されている [12].
これを図 11 で矢印により示す. $\Gamma<2.5$ では上流側円柱, 下流側円柱のどちらも$C_{L1}\sim C_{L2}\sim \mathrm{O}.0-0.07$ となり 1 円柱の値に良く似ている.
しかしながら, これらの値は $3.1<\Gamma<3.5$ で突然に変化し $\Gamma>4.0$ でかなり大きくなる.
$\Gamma>4.5$ ではの値はおよそ 16 であり, これは下流側円柱に働く平均抗力と同じかそれよりも大きい. 下流側円柱については広い間隔比 $(\Gamma>4.0)$ では抗力と同様に大きな揚力も受けているということが結論付けられる. 一方で, 上流
側円柱については下流側円柱よりも大きな抗力を受けるが, 揚力 $(C_{L1})$ はそれほど大きくはない.
一般に, 1 円柱を過ぎる流れのストローハル数 $St$ #までおよそ 0.164である [11]. これは図 12 の矢印で示されている. 図 12 からわかるように 1 円柱の場合とタンデム型 2 円柱の場合どちらも 0.12-0.16 でほぼ同じ値である. $\Gamma<3.5$の狭い間隔比の場合はストローハル数は約 0.12 であり, 広い間隔比 $(\Gamma>3.5)$ に$C_{D2}$$C_{D2}$$\Gamma\sim<3$$C_{D2}$$\Gamma$$\Gamma\geq 4$$C_{D2}$$\Gamma\leq 3.0$

$C_{D2}$$L$$C_{L1}(\circ),$$C_{L}$$C_{L2}$$\bullet$$C_{L}$$Re=1\mathrm{O}\mathrm{O}$$C_{L2}$$Re=1\mathrm{O}\mathrm{O}$$C_{D1}$

17

果$\mathrm{F}$

図 11: 揚力係数の振動 Q.

$Re=1\mathrm{O}\mathrm{O}$.

$\circ$:$C_{L1}$(上流側円柱).

$\bullet$:$C_{L2}$(下流側円柱). $arrow:C_{L}$ ( 円柱).
$1$$\ddagger\tilde{\mathrm{o}}$$\Gamma$

図 12: ストローハル数.
なると 0.16 まで大きくなる.

$Re=1\mathrm{O}\mathrm{O}$.

$arrow:1$$3.1<\Gamma<3.5$ の範囲では,円柱の St.

ストローハル数は急激に変化し, 結果的にこの範囲で多数の解を生じる.
および $St$ は同様に解の多重性を示している. この解の多重と物理量の振動振幅を(下流側円柱後方約 $3.6d$ ) での流速性は図 13 のように代表点の振動振幅は $\Gamma<3.0$
代表物理量にとることにより明確に描くことができる,$C_{D}$$C_{L}$$v_{1}$$\mathrm{P}_{1}$$v_{1}$

では小さく,$\Gamma\sim 2.5$で最小になる. また,$\Gamma>4.0$

では大きくなる.

5. 結言
および $St$ の物理量に起こる不連続な変化は基礎方程式の解が多数存在することが原因であることを示した. 一方, 実験において, この本研究では,$C_{D}$と$C_{L}$

不連続な変化はもっと高いレイノルズ数 $Re\sim 10^{3}-10^{4}$ で指摘されてきた. 本研

究で得た結論は高いレイノルズ数の流れにおいて直接には適用できないが, 高い18$\dot{\mathrm{i}\not\supset}$$\Gamma$

図 13:$v_{1}$の振幅.

$Re=1\mathrm{O}\mathrm{O}$.

レイノルズ数の流れで小さい規模の渦運動を乱流粘性として取り込むことにより,低いレイノルズ数での分岐構造が高いレイノルズ数でも再現されるという再帰構
造をしていることを期待して, ここでの議論が高レイノルズ数においても適用できる可能性を指摘しておく. ただし, このことは今後の研究課題である.

謝辞
この研究を進めるにあたり, 同志社大学大学院生除補義信君と稲生洋一君の協力を得た. 両君に感謝をする. また, この研究は文部科学省および日本学術振興会の科学研究費補助金 (基盤研究 (C) (2)) の支援を受けて行った.

参考文献[1] I. Peschard and P. Le Gal, “Coupled wakes of cylinder,” Phys. Rev. Lett,
3122 (1996).

$77_{7}$[2] M. M. Zdravkovich, “Review of flow interference between two circular cylinders In various arrangement,” J. Fhtids Eng., 99, 618 (1977).
[3] Y. Ohya, A. Okajima and M. Hayashi, “Wake interference and vortex shedding,” Encyclopedia of Fluid Mechanics. Edited by N. P. Cheremisinoff (Gulf
Pub., Houston, 1988), vol. 8, pp. 323-389.

[4] S. Ishigai, E. Nishikawa, K. Nishimura and K. Cho, “Experimental study
on structure of gas flow in tube banks with tatbe axes normal to flow,”
Bull. JSME, 15, 949 (1972).

19[5] A. Okajima, “Flows around two tandem circuiar cylinders at very highReynolds number,” Bull. JSME, 22, 504 (1979).

[6] M. Arie, M. Kiya, M. Moriya and H. Mori, “Pressure fluctuations of thesurface of two circular cylinders in tandem arrangement}’’ J. Ffuid Eng., 105,
161 (1983).
[7] G. Xu and Y. Zhou,$Exp$ . Fluids, 37, 248 (2004).

$\zeta\zeta \mathrm{S}\mathrm{t}\mathrm{r}\mathrm{o}\mathrm{u}\mathrm{h}\mathrm{a}\mathrm{l}$numbers in the wake of two inline cylinders,”[8] W. Eissler, P. Drtina and A. Frohn, “Cellular automata simuiation of ftow
around chains of cylinders,” Int. J. $Num$ . Meth. Fluids, 34, 773 (1992).
incompress[9] S. Mittal, V. Kumar and A. Raghuvanshi,
and staggered arrangement,”
ible flows past two cylinders in tandInt. J. Nurn. Meth. Fluids, 25, 1315 (1997).
$‘\zeta \mathrm{U}\mathrm{n}\mathrm{s}\mathrm{t}\mathrm{e}\mathrm{a}\mathrm{d}\mathrm{y}$$\mathrm{e}\mathrm{m}$[10] J. R. Meneghini, F. Saltara, C. L. R. Siqueira and J. A. Ferrari, . , “Numerical simulation of flow interference between two circular cylinders in tandem
and side-by-side arrangements,” J. Fluids Struct., 15, 327 (2001).
$\mathrm{J}\mathrm{r}$[11] B. Sharman, F. S. Lien, L. Davidson and C. Norberg, “Numerical predictions of low Reynoids number flows over two tandem circular cylinders,”
Int. J. $Num$ . Meth. Fluids, 47, 423 (2005).
of circular cylimder os[12] Y. Tanida, A. Okajima and Y. Watanabe,
cillating in uniform or in wake,” J. Fluil Mech., 61, 769 (1973).
($‘ \mathrm{S}\mathrm{t}\mathrm{a}\mathrm{b}\mathrm{i}\mathrm{l}\mathrm{i}\mathrm{t}\mathrm{y}$[13] K. Ohmi and K. Imaichi, “Vortex wake visualization of two circular cylinProc. 6th Int. Symp. , Yokohama, Japan, 5-9
ders,” in Florn Visualization
October 1992 (Pacific Convention Plaza, Yokohama, 1992), pp. 322-326.
$IV.\cdot$[14] D. Biermann and W. H. Herrnstein Jr., “The interference between struts
in various combination,” National Advisory Commitee for Aeronautics,
Tech. Rep., 468, 515 (1933).

[15] P. Le Gal and M. P. Chauve, ‘ Coupled wakes betfind two circular cylinders,”
Phys. Rev., A 41, 4566 (1990).
$\langle$of the flow around two circular cylinders ar[16] T. Igarashi,
ranged in tandem (1st report),” Bull. JSME, 24, 323 (1981).
$‘\zeta \mathrm{C}\mathrm{h}\mathrm{a}\mathrm{r}\mathrm{a}\mathrm{c}\mathrm{t}\mathrm{a}\mathrm{r}\mathrm{i}\mathrm{s}\mathrm{t}\mathrm{i}\mathrm{c}\mathrm{s}$[17] J. F. Thompson, F. C. Tham es and C. W. Mastin, “Automatic numericalgeneration of body-fitted curvilinear coordinate system for field containing

20

any number of arbitrary two-dimensionai bodies,” J. Comp. Phys., 15, 299(1974).

[18] J. L. Steger and R. L. Sorenson, “Automotic mesh-point clustering neara boundary in grid generation with elliptic partial differential equations,”
J. Comp. Phys., 33, 405 (1979).
[19] C. Mathis, M. Provansal and L. Boyer, “The Benard-von Karman instability:
an experimental study near the threshold,” J. Phys. Left., 45, L-483 (1984).
[20] C. P. Jackson, “A finite-element study of onset of vortex shedding in flow pastvariously shaped bodies,” J. Fluid. Mech., 182, 23 (1987).
[21] S. Taneda, “Experimental investigation of the wake behind cylinders and platsat low Reynolds number,” J. Phys. Soc. Japan, 11, 302 (1956).