• 検索結果がありません。

低いレイノルズ数領域での円柱背面の伝熱流動

N/A
N/A
Protected

Academic year: 2021

シェア "低いレイノルズ数領域での円柱背面の伝熱流動"

Copied!
9
0
0

読み込み中.... (全文を見る)

全文

(1)

低いレイノルズ数領域での円柱背面の伝熱流動

著者

布施 肇, 鳥居 修一, 二見 竜生, 加治屋 重昭

雑誌名

鹿児島大学工学部研究報告

32

ページ

1-7

別言語のタイトル

FLOW AND HEAT TRANSFER CHARACTERISTICS BEHIND

CIRCULAR CYLINDERS IN LOWER REYNOLDS NUMBER

REGION

(2)

低いレイノルズ数領域での円柱背面の伝熱流動

著者

布施 肇, 鳥居 修一, 二見 竜生, 加治屋 重昭

雑誌名

鹿児島大学工学部研究報告

32

ページ

1-7

別言語のタイトル

FLOW AND HEAT TRANSFER CHARACTERISTICS BEHIND

CIRCULAR CYLINDERS IN LOWER REYNOLDS NUMBER

REGION

(3)

低いレイノルズ数領域での円柱背面の伝熱流動

布 施 肇 ・ 鳥 居 修 一

二見竜生*・加治屋重昭**

FLOWANDHEATTRANSFERCHARACTERISTICS mmHINDCIRCULARCYLINDERSINLOWERREYNOLDS NUMBERREGION HajimeFUSE,ShuichiTORII, TatuoFUTAMI,andShigeakiKAJIYA Theaimofthispaperistoclarifyflowandheattransfercharacteristicsbehindcircularcylin-ders・Ourmainattentionisfocusedontheexaminationofthevalidityofheattransfercorrelations, whichhavealreadybeenproposedbyserveralauthors・Heattransferattherearstagnationpointand streamwisevelocityweremeasuredundertheConditionthatfree-streamturbulenceandblockagefac-toraremaintainedconsiderablylowerlevel・TheexaminedrangeoftheReynoldsnumberwasfrom 7000to22000・ IntherelativelyhighReynoldsnumberheattransfercoefficientsarepredictedbymeansofthe existingcorrelations,whiletheirpredictionaccuracyisinsufficientinlowerReynoldsnumbercases・ Acorrespondingexperimentalresultoftheseparatedshearlayershowsthatthelaminar-to-turbu-lenttransitionregionisshiftedintheupstreamwisedirection,astheReynoldsnumberincreases・It wasfound,therefore,thatthetransitionoftheseparatedshearlayerhasasignificanteffectonthe heattransferbehavior. 1 . 緒 目 二次元円柱まわりの流れは熱線風速計や熱交換器な どに見られ,熱設計の見地からその流動特性について の実験的研究が数多く行なわれている。しかしながら、 後流の非定常性や背面近傍の流体力学的挙動が複雑で あるために,理論解析がかなり困難であり,得られた 実験結果から統一した見解は得られていない。これは, それぞれの実験条件すなわち主流乱れ,ブロッケージ 比(円柱直径Dと風洞壁間の距離Lとの比D/L)の 効果に起因すると考えられる。それに関して前者では, 円柱後方岐点の熱伝達が主流の乱れ強さの増加に伴っ て向上することが指摘されており1),2).3),Kestin4) はこれについての研究成果を系統的に纏めている。ま た後者では,プロッケージ比が大きくなる場合,円柱 後方岐点のヌセルト数は風洞壁の干渉効果によって増 加することが報告されており5)・6),この効果を考慮し た熱伝達整理式がある程度整備されている。 Richardoson7)は,従来の実験結果を整理してブロッ ケージ比が零である場合の円柱後方岐点の熱伝達整理 式を求めている一方,五十嵐ら8)は平板背面のはく離 流に適用した体積モデルを円柱に採用して対応する整 理式を導出している。また,プロッケージ比を考慮し た円柱背面熱伝達に関しては桧和田ら5)の研究が挙げ られる。これらの整理式は,円柱まわりの流れの様相 には視点を向けずにレイノルズ数が10000以上を対象 とした伝熱測定から導かれたものであり,レイノルズ 数の10000付近あるいはそれ以下での予測精度につい ては検討されていない。一般に熱伝達整理式は簡便で あるので,適用できる範囲の流れと伝熱の様相が十分 把握できれば,はく離流を伴う機器の熱設計に上記の 整理式を用いる試みは有効であると考えられる。

(4)

︵蛸︶め 2 ︹四︶ぬ 以上のことから本報は,低いレイノルズ数領域を対 象とした伝熱流動実験を行ない円柱背面の流れと熱伝 達の両特性を明らかにするとともに,従来報告されて いる熱伝達整理式について検討することを目的として いる。 いるのでここでは省略する。円柱まわりの壁面圧力は, 加熱円柱と同直径の銅製円柱に0.2mmの圧力孔を開け て,これを円柱前方岐点から180.まで回転させて求め た。さらに,主流と円柱背面の流れの測定にはI型プ ローブを用い,得られた出力信号を市販のFFT(Fast FourierTransform)解析器で整理してそのパワース ペクトルを求めた。各試験円柱の直径は7mm,8mmお よび10mmであり,これらに対するプロッケージ比はそ れぞれ0.052,0.059,0.063であった。レイノルズ数 Reは7000から22000までの範囲であった。 実験では,円柱上流側の主流方向の速度と乱れ強さ を管断面全域にわたり一様であると仮定しているが, 双方の測定結果から壁面近傍を除く約94%の領域で一 様性が保持されていた。表1は,加熱円柱のかなり上 流域で測定した流路断面中心付近の乱れ強さを示した ものである。Kestinら'0)の結果に従えば,主流乱れ の影響は殆どないものと考えられる。一方,著者ら9) は,主流乱れ強さが影響しない範囲でも高周波成分が 乱れに含まれているような場合は熱伝達に影響するこ とを指摘している。図1は,レイノルズ数が10000に おいて,D=7mm,D=8mmおよびD=10mInの主流パ ワースペクトルである。図から分かるように,高周波 成分は主流乱れに含まれていないのでこの影響はない と判断した。円柱を固定している流路測壁近傍では, 軸方向熱伝導損失によって円柱壁に温度勾配が現われ るが,この領域を除く約86%の部分で一定壁温が得ら れた。一方,円柱まわりの壁温は円柱中央で測定した 記 号 圧力係数=(P−P。。)/(IOU畠/2) 円柱直径 周 波 数 円柱まわりの局所熱伝達係数 長方形流路壁間の高さ ヌ セ ル ト 数 = h D / A 角度8の位置での円柱表面圧力 円柱の影響を受けない主流の圧力 レ イ ノ ル ズ 数 = U D / 〃 修正主流レイノルズ数=U*D/'ノ パ ワ ー ス ペ ク ト ル 時間 主流温度 円柱前方岐点の表面温度 温度比=(T8−T。。)/(To‐Too) 角度βの位置での円柱表面温度 X方向の乱れ強さ=((u'2)1'2/U。。) 主流方向の速度変動成分 主 流 速 度 修正主流速度 =UooI1+0.293(D/L)+0.822(D/L)2 +1.954(D/L)31 流れ方向の座標あるいは円柱中心から流れ方 向への距離 流れと垂直方向の座標あるいは円柱中心から 流れと直角方向への距離 流体の熱伝導率 流体の動粘‘性係数 右まわりを正として,円柱前方岐点から測っ た角度 2.

、“u

函eefl *いり

函orβu

函*

CDfhLNPPRRStTTTTT㎡UU

表 1 各 円 柱 の 主 流 の 乱 れ 強 さ X D(m、) 7 8 10 (%) ,0.350 ,0.350 ・0.350 T u O、180 0.180 0.170 Y 、八〃〃U 鹿 児 島 大 学 工 学 部 研 究 報 告 第 3 2 号 ( 1 9 9 0 ) f H z O 缶 H z O 歪Hz 図1各円柱の主流のパワースペクトル 贈︶の 3 . 実 験 装 置 と 実 験 方 法 2枚の多孔板からなる風洞で整流された空気が、風 洞出口に取り付けた長方形流路に入り、その入口から 300mmの位置に水平に固定された加熱円柱のまわりを 流れ外部へ放出される。使用した加熱円柱の詳細な構 造と円柱まわりの壁温の測定方法は文献(9)で説明して D=8 1m Re=1.0×1〔 I 4 ’ 10Hz −0.0361− 1 0 D=10mm Re雲1.0xlO 1 1 4 1 10Hz -0.0561 − D=7mm Re=1.0×11 ■ AO 10Hz =0.0331 l ▲

(5)

3 2 ので,その精度は上記の熱伝導損失によって左右され ることはない。圧力分布を除く上述の各測定結果の中

で,対応する円柱前方岐点のヌセルト数が理論値'1)

と約3%以内で一致したもののみを実験結果として採 用した。 によれば,一節で述べたように円柱後方岐点のヌセル ト数はブロッケージ比が大きくなるにつれて増加する 傾向にある。これに従えば,図中の実験結果は何れも 整理式より多少高めに位置していなかればならない が,比較的低いレイノルズ数領域の実験点はかなり低 い。一方,レイノルズ数が15000以上では,何れの円 柱の実験結果も上記の整理式と5%以内で一致してい る。各円柱ともレイノルズ数の増加に伴い整理式に漸 近する傾向は同じであるが,D=7mmと8mmの実験点 はレイノルズ数が約13000以上で急激に整理式に近づ くという点で10mmの場合と多少異なっている。従って, 円柱後方岐点のヌセルト数を上記の2つ式の何れかで 整理する際,その予測精度はレイノルズ数が比較的低 い領域では芳しくないことが分かった。そこで,円柱 まわりの壁温と圧力の分布および円柱背面の流れか ら,上述の結果が現われる背景について検討してみよ う。 4.実験結果と考察 4 . 1 円 柱 後 方 岐 点 の ヌ セ ル ト 数 図2は各円柱の後方岐点のヌセルト数Nuをレイノ ルズ数Reとの関係で示したものである。図中の破線 は上記で触れた五十嵐ら8)の整理式であり,次式で表 わされる。 Nu=0.172Re2/3 (1) また、桧和田ら5)も先に述べたように,円柱後方岐点 の熱伝達整理式としてプロッケージ比を考慮した次式 を導出している。 Nu=0.176Re*2/3 (2) この式にD/L=0を代入した結果すなわちプロッ ケージ比が零の場合の熱伝達整理式を図中に実線で併 記している。これらの整理式に対して,得られた実験 結果はブロッケージ比の影響を多少含んである。本実 験のブロッケージ比を式(2)に代入した際,得られたヌ セルト数はブロッケージ比が零の場合に比べて約2% 増加するが,この程度のヌセルト数の増加は無視でき るものと考えられる。円柱後方岐点のヌセルト数とプ ロッケージ比の関係を示した桧和田らの論文(5)の図16 4 . 2 円 柱 ま わ り の 壁 温 と 圧 力 の 分 布 図3(a),(b),(c)は各円柱における円柱まわりの壁温 分布を示している。図中の縦軸は,角度βの円柱表面 と主流の温度差を円柱前方岐点と主流の温度差で除し た無次元温度である。円柱前面の壁温は約90°まで急 激に上昇しているが,これは各円柱ともほぼ同じ曲線 で表わされる。また,円柱前面で見られるこのような 分布形はレイノルズ数が増加しても大きくは変らな い。一方,円柱後方の壁温分布形は円柱径とレイノル ズ数の違いによって異なる。円柱径が7mmと8mmの場 合,Re=7000の壁温は約120.付近までやや増加し, その後円柱後方岐点に向けてゆるやかに減少してい る。Re=15000までの壁温は120。付近で凹部が現われ た後,Re=7000の場合より激しい減少傾向を示して いる。ここで,壁温分布に凹部が見られるのはウニ次 うず'2)が現われその影響が顕著になったためである と考えられる。さらにレイノルズ数が増加した場合, 壁温は後方よどみ点まで単調に減少し,その傾きはこ れまでのレイノルズ数の場合よりも大きい。同様の温 度分布は円柱径が10mmの場合にも見られるが,壁温が 円柱後方岐点まで単調に減少する傾向はレイノルズ数 が13000で既に現われている。こうして見てくると, D=71,,,,=8mmおよびD=10mmの円柱後方岐点の ヌセルト数がレイノルズ数の増加に伴って整理式に漸 近する先の結果は,円柱背面の壁温分布の変化すなわ ち壁温が円柱後方岐点に向けて急激に減少することに 3 UC .△ IgarashietaL O、l72Re2/3 2 6 p△ Hiwad O、176刀 昼 布施・鳥居・二見・加治屋:低いレイノルズ数領域での円柱背面の伝熱流動

08

1 .Z 堂?

七○ ○ 5 4 0 0 0 6 8 1 0 4 2 3 Re 図2円柱後方岐点のヌセノレト数

43

10

(6)

1

i

2

2

o

o

o

1

4 0 対していると考えられる。 図4は,各円柱において,レイノルズ数の違いによ る円柱まわりの圧力分布の変化を示している。はく離 点は,円柱表面上の摩擦応力の測定からあるいは円柱 に塗料を塗りそれがはく離点に吹き寄せられることか ら求められるが,ここではそれらの実験を行なってい ないのではく離点の位置を正確に求めることはできな い。しかしながら,はく離点が圧力の最小値から急激 に上昇している領域にあることは確かである。またこ のことから,実験で現われたはく離が層流はく離であ ることは明らかである。各円柱の円柱まわりの圧力分 布は何れの円柱も同じような形を示しており,円柱後 方の圧力はレイノルズ数が増すにつれて減少してい る。そのような傾向は,レイノルズ数の変化に対して 緩やかなものであり,先の壁温分布の場合のように分 布形が各レイノルズ数によって大きく異なることはな い。 2. 2 占 1 ]Ogg8 1. 0 3 0 6 0 9 0 1 2 0 1 5 0 l 8 0 eo (a) 4.3はく離せん断層の特性 前節から,レイノルズ数の違いによって円柱まわり の壁温分布形に大きな変化が現われることが分かっ た。これは円柱背面の流れの様相がレイノルズ数に よって異なったことが考えられる。また著者ら9)は, -0.5 2 1.0

§

僅 回。 1 1.0

符8

1. 1.0

附 館

0 3 0 6 0 9 0 1 2 0 1 5 0 . 1 8 0 e、 (b) 一ARM山△︹H図﹁1

謎帥.”]叩

0.5

鐘室遮墓

! 】 、=lOnun Re 00.7×10'+ △1.0×10'+ □1.3×lOlO ▽1.5×10'+ ◇1.8×101+ ①2.4×104 鹿 児 島 大 学 工 学 部 研 究 報 告 第 3 2 号 ( 1 9 9 0 ) 1. 0 180

H﹄ 0 3 0 8 -1.5 1. 5

'

L)=lC 6 0 1 2 0 l80 eO 円柱まわりの圧力分布 -1.0 0 ● 《 図4 6 0 9 0 1 2 0 1 5 0 8。 (c) 円柱まわりの壁温分布 0 5

●●

l b1 ,トー 図3

!

、0 、5 0

、=8m、 Re 111111 000000 44.4444 ×××××× 703582 0●●●●● 011112 0△ロ▽◇の

ー § § )◎ ▲ LD 0 。 ◎ O V | ー コム ロ △ 貼 F h − − ムムィ ロロ )。。◎go()o99

9

、クロ の の1 D=7mm e R 444444 000000 111111 ×××××× 703582 ●●●●●● 011112 ○△ロ▽◇①

。△

。△ロ

○△ロ 。△ロ

。△ロ

ロI 61△lpI 。△ロ 0△□ ◎ 00 日 ▽

① 、〃 の ▼ ▼ ▼ ▼

◇砂

| hr ① 4 の ‐ 心 ▽ ◇ ▽▽,

I

b

①の − ①両

(7)

X/、 5 0 . 5 1 . 0 1 . 円柱背面の熱伝達ははく離せん断層内の特性と密接な 関係があることを明らかにしているので,以下ではこ れについて検討を試みる。 先ずはく離せん断層の形状を求めるには円柱後方の 速度分布が必要である。この測定はプローブ先端を上 流側に向けこれをY方向に移動して求めた。図5は, D=7mmでレイノルズ数が10000の場合に,後流の速 度分布が下流方向に変化する過程を表わした一例であ る。特に円柱背面近傍では流れが逆流しているので, この領域の測定値を主流方向のものと看倣すには多少

問題がある。一方,安達ら'3)の定義に従って,せん

断層の外側,内側境界を速度が極大,極小値をとる位 置とし,この領域をせん断層の厚さと定めた。従って, 先に触れたように主流速度の測定精度は円柱背面近傍 で多少劣ってはいるが,これによって求めたせん断層 の形状に大きな違いが現われるとは思われない。 図6(a),(b),(c)は,3つの異なったレイノルズ数に ついて得られた各円柱のはく離せん断層の形状を比較 したものである。レイノルズ数が比較的小さい場合, せん断層は下流に向けて多少厚くなってはいるもの の,何れの円柱のものもかなり薄い。一方レイノルズ 数が13000の場合のせん断層は各位置で対応する先の 結果より厚く,後方に向けてのせん断層の広がり方は 円柱径が大きいほど激しい。ざらにレイノルズ数が増 加した場合(Re=22000),せん断層は各円柱ともほ ぼ同じ形状であり,その厚さはレイノルズ数が13000 でD=10mmの場合に近い。すなわち,Re=13000以上 では,D=71nmとD=8mmのせん断層の形状はD= 10mmの結果に急激に近づいていることが窺える。以上 の結果から,レイノルズ数の違いによる各円柱のせん 断層の変化が明らかになったので,次にその内部の様 相を検討する。 図7(a),(b)は,D=7mmでレイノルズ数が10000の 場合に,円柱中心から0.1直径だけ下流で測定したせ ん断層中心付近の速度変動波形とそのパワースペクト ルである。ここで,図7(b)の縦軸の目盛は任意である。 0.8 布施・鳥居・二見・加治屋:低いレイノルズ数領域での円柱背面の伝熱流動 ロヘン 0.5 0 ,0.5 (a) 1.0

/

0. 00OO、 ロヘン 1.5 0 Re=1.0×101+ D=7mm 06 8コヘヨ 5 8

1.0 3 図6 (b) 0.8 0.5 DO246802 ノ...。... ロヘン 1.3 Q・△□①▽O︾小 0.5 0.5 0 0.5 1.0 (c) X/, 0 1.0 各円柱のせん断層の形状 Y/, 後流の速度分布 図5 Re=0 D(四) ○ 7 □ 8 △ 10 Re=2● 言 詔 、 D(−1 ○ 7 ロ 8 △ 10

臣 :

Re=1 一一 、(−1 ○ 7 □ 8 △ 10 司 出 ー 垂

(8)

図7(a)の正弦波形から求めた渦放出周波数は,図7(b) のスペクトルの最大値のそれと対応している。こうし て求めたパワースペクトルの最大値すなわち渦放出周 波数を持つエネルギの流れ方向の変化を,図6の結果 と対応させて図8(a),(b),(c)に纏める。前で述べたよ うに、円柱まわりの境界層は層流はく離しているので, はく離直後のせん断層は層流である。周知のとおり, 層流から乱流への流れの遷移は一点で起こるのではな く流れ方向にある幅を持つ領域内で現われるので,乱 流遷移点を厳密に定めることはできない。一方,層流 はく離した境界層が層流から乱流に遷移する際, Bloor'4)は乱流への遷移に先駆けて流れの中に遷移波 と呼ばれる比較的規則正しい速度変動が現われること を報告している。そこでこれに従って,図7(a)の波形 に高周波の波形が現われ始めた位置を乱流遷移の開始 点とし,それが波形全体に現われた位置を乱流終了点 として,両点を半黒と全黒で図中に併記し遷移領域を 表わす。レイノルズ数が7000の場合,何れの円柱とも 乱流遷移はかなり下流位置で起こっている。一方, Re=13000の場合の乱流遷移は,各円柱ともRe= 7000の場合より上流側へ移動しており,特にこの傾向 はD=10mmで顕著である。更にレイノルズ数が増加し た場合,何れの乱流遷移も円柱にかなり接近している。 せん断層が乱流に遷移すれば,主流とせん断層,せん 断層と死水域の間の運動量輸送が活発に行なわれてせ ん断層の速度勾配が早期に平滑化され,せん断層の渦 度もこれに伴って拡散されるので,せん断層は急激に 厚くなると考えられる。またこれに伴って,円柱背面 からの熱の拡散も活発になるであろう。従って,先に 述べたように,高レイノルズ数域での円柱背面の壁温 が後方岐点に向けて急激に減少したのは,乱流遷移が 早期に起こったことに起因していると考えられる。さ らに,D=7mmとD=8mmの遷移領域がレイノルズ数 の増加によって急激に円柱に近づいていることは,実 験点がRe=1,000以上で整理式に漸近することに対 応している。一方,D=10mmのせん断層は他の円柱よ 6 0 1.5 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 X/、 (a) 6.0V ● 0 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 X/、 (c) 図8せん断層内パワースペクトルの主流方向の変化 ︵蛸︶② ︵単︶の . O −6.0V 1.0 0 . 0 1 0 . 0 m S O f H z t ( a ) 渦 放 出 周 波 数 ( b ) パ ワ ー ス ペ ク ト ル 図7はく離せん断層内の速度変動 0.5 色彦三.8 鹿 児 島 大 学 工 学 部 研 究 報 告 第 3 2 号 ( 1 9 9 0 ) 0.5 2.0 X/、 ︵蛸︶の 1.5 1.0 0.5 0 1.5 (b) 2.0 0 . 2 0 . 4 0 . 6 0 . 8 1 ︵串︶の 1.0 / 、 /h 〃 、 へ /( ハ 、 ヅ 、ダ 、ノ 、=7m、 Re=1.0×10 X/D=0.1 4 20Hz 0.07 ヘーーI 75 三 590Hz 0.0679 l L 一Re= 一 一 一 U 0.7×101+ ' ○ □ D(函) 7 8 一一一一 一一一一 ﹃﹃﹃一 Re=1.U 3×101+

'

│ ロ DIm1 7 8 ∼

(9)

布施・鳥居・二見・加治屋:低いレイノルズ数領域での円柱背面の伝熱流動 7 りも早い段階で乱流に遷移しているので,ヌセルト数 もこれに対応してかなり低いレイノルズ数領域から熱 伝達整理式に近づいている。何れにしても,円柱後方 岐点の実験値が熱伝達整理式で纏められる結果は,せ ん断層の乱流遷移がかなり円柱に近づいて起こったも のであると考えられる。 以上の状況から判断して,円柱後方岐点の熱伝達は せん断層の乱流遷移と関係があり,用いた熱伝達整理 式で纏められるものとそうでないものとは円柱背面の 流れの様相がかなり異なっていることが分かった。 5 . ま と め レイノルズ数が比較的低い領域において円柱を加熱 した伝達実験を行ない,以下のような知見を得た。 (1)レイノルズ数が10000以下の領域では,円柱後 方岐点のヌセルト数は式(1),(2)では整理しにくい。 (2)円柱まわりの圧力分布の変化はレイノルズ数の 変化に対して壁温分布ほど敏感ではない。 (3)式(1)あるいは式(2)で整理できる結果はせん断層 の乱流遷移が円柱に接近しているのに対して,そ のでないものは乱流遷移がかなり下流域で行なわ れることから,円柱背面熱伝達はせん断層の特性 に左右される。 文 献 1)Petrie,A・MandSimpson,HC.,Int.J・Heat MassTransfer,15−8(1972),1497. 2)Boulos,M1・andPei,,.C、T、,Int、JHeatMass Transfer,17−7(1974),767. 3)Yardi,N・RandSukhtame,S.P.,Proc、6thInt, HeatTransferConf.,FC(b)−29(1978),347.

4

K

e

s

t

i

n

,

,

A

d

v

a

n

c

e

s

i

n

H

e

a

t

T

r

a

n

s

f

e

r

,

v

o

l

AcademicPress,NewYork,1966.

5)桧和田・ほか3名,機論,42−360,2部(1967),

2481.

6)桧和田・揚淵,機論,46-409,B(1980),1750.

7

R

i

c

h

a

r

d

s

o

n

,

.

.

,

C

h

e

m

E

n

g

r

S

c

i

.

1

9

6

3

)

149.

8)五十嵐・平田,機論,39−322,2部(1973),1890.

9)布施・ほか2名,機論,50-453,B(1984),1302.

1

0

K

e

s

t

i

n

,

a

n

d

W

o

o

d

,

,

,

T

r

a

n

s

A

S

M

E

,

S

e

r

C,93−4(1971),321.

1

1

G

o

l

d

s

t

e

i

n

M

o

d

e

r

n

D

e

v

e

l

o

p

m

e

n

t

s

i

n

F

l

u

i

d

Dynamics,(1964),632,Dover、

1

2

S

o

n

,

.

a

n

d

H

a

n

r

a

t

t

y

,

,

F

l

u

i

d

M

e

c

h

.

,

(1969),35(1969),353.

13)安達・加藤,日本航空宇宙学会誌,23-256,

(1975),45.

1

4

B

l

o

o

r

,

,

J

F

l

u

i

d

M

e

c

h

.

,

(

1

9

6

4

)

,

参照

関連したドキュメント

振動流中および一様 流中に没水 した小口径の直立 円柱周辺の3次 元流体場 に関する数値解析 を行った.円 柱高 さの違いに よる流況および底面せん断力

Due to the monotone increase with time and monotone decrease with radial dis- tance of the temperature (Figure 4.1), the gas is evaporated from the inner cylinder and in the course

In this paper, we discuss the nature of incompressible conducting fluid flow around a circular cylinder in the presence of an external magnetic field for a range of Reynolds

Max-flow min-cut theorem and faster algorithms in a circular disk failure model, INFOCOM 2014...

So, the aim of this study is to analyze, numerically, the combined effect of thermal radiation and viscous dissipation on steady MHD flow and heat transfer of an upper-convected

On the other hand, the magnitude of the cross-flow velocity increases with the increase in either suction pa- rameter or frequency parameter, while it increases near the

In this paper, we discuss the nature of incompressible conducting fluid flow around a circular cylinder in the presence of an external magnetic field for a range of Reynolds

27 found that multiple solutions exist for a certain range of ratio of the shrinking velocity to the free stream velocity which again depends on the unsteadiness parameter for