Thermal Science & Engineering Vol.15 No.4 (2007)

矩形管内脈動流の熱･物質輸送特性
第1報

脈動周波数が管内に設置された平板面の熱伝達に及ぼす影響齊藤 弘順†

内田 浩二†Heat and Mass Transport Characteristics of Pulsating Flowin a Rectangular Duct
Part 1. Effect of Pulsating Frequency on Heat Transfer of a Flat Plate Installed in the Duct
Hironori SAITOH†

and Kouji UCHIDA†

Abstract
This study deals with flow behavior and associated heat transfer in pulsating duct flow. As well as mechanicalengineers, researchers in bioengineering field also focus on heat and mass transport characteristics of pulsating flowfrom the viewpoint of life science. In spite of more than thirty-year history of pulsating flow and heat transfer studies,
flow structure with periodical flow rate change and its effect on heat transfer are still at question. Researchers showed
conflicting results. Some of them showed heat transfer enhancement due to flow pulsation, the others reported heattransfer prevention or no impact of flow pulsation on heat transfer. In order to make clear the complicated unsteadyflow structure with periodical velocity fluctuation and its effect on heat transfer, it must be necessary to havemacroscopic observation and image processing analysis of both the flow and temperature fields. The objective of thisstudy, therefore, is to make clear the pulsating flow structure and heat transfer characteristics based on the visualizationsof flow and temperature fields. Particle Image Velocimetry (PIV) and color schlieren visualization was applied to flow
and temperature fields, respectively. Results showed that flow pulsation induced complicated spatial velocity
distributions. Series vortices were formed in flow direction and they moved in span-wise direction from the sidewalls to
the center of the duct. Such flow behavior caused locally high and low speed regions. From these results, it is concluded
that heat transfer enhancement or prevention depends on increase or decrease of local flow rate induced by flow
pulsation.
Key Words: Pulsating flow, Convective heat transfer, Temporal and spatial velocity fluctuation, Colorschlieren method, Flow visualization , PIV , Vortices

記 号
A
:
cp
:
de
:
f
:
Henh
:
h
:
Nu
:
q&
:
Ref
:
Rem
:

[m2]
断面積
[kJ/(kg⋅K)]
定圧比熱
[m]
等価直径
[Hz]
周波数
[-]
伝熱促進率
熱伝達率
[W/(m2･K)][-]
ヌッセルト数[W/m2]
熱流束[-]
レイノルズ数(代表流速＝速度振幅)[-]
レイノルズ数(代表流速＝時間平均
流速)

T
uf
um
w
X
Y
δ t′

:
:
:
:
:
:

温度
速度振幅
時間平均流速平板幅
流路方向座標流路幅方向座標: 温度境界層厚さ

[K]
[m/s]
[m/s]
[m]
[m]
[m]

λ
ν
ρ

: 空気の熱伝導率[W/(m･K)][m2/s][kg/m3]: 空気の動粘度: 空気密度

受付日: 2007 年 3 月 12 日, 第 44 回日本伝熱シンポジウムより受付, 担当エディター: 千田
†

[m]

衞

そうじょう

崇 城 大学（旧名：熊本工業大学）工学部機械工学科 (〒860-0082 熊本県熊本市池田 4-22-1)- 211 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

添字
b

:

バルク

エンジンなど往復動内燃機関についても，燃料や作

cross

:

断面

動ガス（空気と燃焼ガス）といったエネルギー変換

in

:

入口

の本質に関わる圧縮性流体の挙動は脈動流が多いこ

m

:

平均

とに改めて気付くと同時に，この往復動という一見

max

:

最大

不恰好な印象を与える機構だからこそ，その潤滑特

min

:

最小

性を生かし“作動ガスの完全密封（漏れ無し）およ

w

:

壁面

びガス流動零（熱伝達損失無し）で静圧膨張”といたきっかけである．ガソリンエンジンやディーゼルうエネルギー変換の理想に近づけることができ，そ１

緒論の高いエネルギー変換効率故に現在でもその優位性を保っているように思われる．脈動流に関して学術

1.1 目的
本研究は「“生体の血流による体温調節のメカニズム”の“熱機関･熱機器の脈動流による温度場のア
クティブコントロール”への応用」を視野に入れた生命科学と熱流体工学の複合ともいえる基礎研究である．矩形流路に脈動流を供して時間平均流量，脈動周波数および振幅が流れの状態と壁面熱伝達に及ぼす影響を流れ場と温度場の可視化実験により調べ，結果の画像解析から脈動流の熱･物質伝達および輸送特性を定量評価することが本研究の目的である．

的な知見の積み上げにより現在の内燃機関の姿が確

1.2 着想の経緯と研究意義
脈動流場における熱･物質輸送現象に関する研究に代表されるような可視化画像解析によるデータのは，最近では殊に血流や呼吸など生体内の物質輸送

象解明が可能になってきている．
そこで，流れ場と温度場の可視化実験ならびにその画像解析により脈動流の熱･物質伝達輸送特性に関し，熱･物質伝達率および輸送効率の観点からそのメカニズムを探ろうとするものである．脈動流の輸送特性を系統的に整理できれば，温度場のアクティ
ブコントロールによる熱機関や熱機器の性能向上，医療現場においては例えば人工心肺の性能向上（人
体（細胞）への負担軽減）など応用例は幅広く，本研究の意義は極めて大きいと考えている．

特性として注目されており，動脈･静脈瘤等の発生メ
カニズムおよびその予防･治療法への応用を目指したものとして医工連携テーマの代表例の一つである．
生体内においては体液の流れ（血液，呼気，尿など）はほとんどが脈動流だが，その形態は様々である．
長い進化の過程で現在のような体液輸送システムが出来上がっている事実から，おそらくは脈動流の方が定常流よりも細胞に負担をかけず体温調節の制御性にも優れた特性を有するものと考えられる．しかし，輸送管の多段階（生体工学分野では多世代と表現）にわたる分岐や管自体の柔軟性およびそのぜん立した訳ではない．言わば，先の生体の体液輸送シ
ステムの如く何故かは不明だが，最も効率的な方法になっているのと同様であるように思われる．脈動流に関する流れと熱伝達の研究も多数報告されているが[1-18]，研究者によって相反する結果（伝熱促進／抑制）が示されていることもあり[1-5]，未だ系統的な知見が得られていない．この原因の一つは実験
計測の困難さにあるように思われる．しかし，PIV定量化技術の進歩により非定常流れ場･温度場の現

２

伝熱実験

動運動によってその輸送メカニズムは極めて複雑と

2.1 実験装置および方法
図 1 にテストセクションの概略図を示す．テストなり，「なぜ毛細血管や肺胞など微細な管路内を血セクションは 100 mm×100 mm の断面を有する矩形

液や呼気はスムーズに（低圧力損失で）流れるの管とその中央に設置された加熱平板から成り立ってか？」といった素朴な疑問に対しても学術的に十分いる．加熱平板の両面に厚さ 0.02 mm のステンレスな説明ができるまで現象解明には至っていない．

箔（ヒーター）を貼り付け通電加熱することで一様

かねてより医工連携の研究について興味があった

熱流束条件とした．平板の片面を伝熱面とし，もう事に加え，これまでの研究（往復動内燃機関の燃焼一方の面を平板をよぎる熱の逃げを無視し得るまで技術）を通して「脈動流の積極利用およびその能動小さくすることを目的としたガードヒーターとして的制御によって各種熱機関や熱機器の性能向上においる．伝熱面には流れ方向に 7 点， 0.2 mm の K 型

けるブレークスルー技術へ繋がる可能性があるので熱電対を貼り付け温度計測を行った．また矩形管をはないか」という見解に至ったことが本研究を始め構成する 4 枚の板のうち対向する 2 枚は透明なアク- 212 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

リル製であり，カラーシュリーレン法による温度場（正確には密度場）の可視化も同時に行えるようになっている．
図 2 に実験装置全体図を示す．空気（試験流体）はシロッコファンにより駆動され，流路途中に設けられたバルブの開閉度の時間変化によって脈動を与えられた後，テストセクション内に導かれ，サージ
タンクおよびパイプを通って大気中に放出される．
時間平均流量はサージタンク後のパイプにおいて熱線流速計による速度計測値を基に算出した．流れに脈動を与えるバルブは偏芯カムを利用して往復運動する機構となっている．それに伴いテストセクショ
ン上流の流路断面積が時間とともに正弦波変動することになる．以上の機構により本実験では脈動のモードは連続的正弦波変動としている．図 3 に時間的流量変動のコンセプト図を示す．実験パラメータは
時間平均流量と脈動周波数とし，それぞれシロッコ
ファン出力および偏芯カムの動力源であるモータの回転数により設定した．また，振幅（カムの偏芯量に依存）は一定とした．

Fig. 1 Test section for temp. measurement andschlieren visualization.

2.2 データ整理
時間平均流量および振幅に対する無次元パラメー
タとしてそれぞれ Rem および Ref/Rem を定義した．Remおよび Ref はそれぞれテストセクション内の時間平
均流速 um および速度振幅 uf を代表速度としたレイノルズ数である．それぞれの定義式を式(1)および式(2)に示す．

Rem u m d e /

(1)

Re f

(2)

u f de /ここで速度振幅 uf は次式で表される．

uf

u max u m

u m u min1
( u max u min )2

Amax Aminum
Amin Amax(3)

式(3)において A は脈動発生用バルブによって変動

Fig. 2 Apparatus for heat transfer experiments.

する流路断面積を表している．従って速度振幅 uf を
時間平均速度 um で除した値を無次元振幅と定義す
ると無次元振幅は上述のとおり 2 種類のレイノルズ数の比 Ref/Rem で表現される．
定常流における熱伝達特性は壁面温度計測を基に式(4)で示される局所熱伝達率を算出し，局所ヌッセ
ルト数 Nu(X)にて定量評価した．
h( X )

q& /(Tw ( X ) Tb ( X ))

2.000
1.800
1.600
1.400
1.200
1.000
0.800
0.600
0.400
0.200
0.000
0

(4)

90

1.000
0.900
0.800
0.700
0.600
0.500
0.400
0.300
0.200
0.100
0.000
180 270 360 450 540 630 720

Cam angle[deg]

Fig. 3 Design concept of flow rate fluctuation.

- 213 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

Table 1 Experimental conditions.

Y

920, 2000, 3000, 4000

Re f /Re m

0.72

f [Hz]

0(Steady flow), 1.0, 2.0, 3.0, 5.0

δt’ (X )δt(X )Re m

X
TW(X) － Tb(X)Fig. 4 Definition of thermal boundary-layer.
0Hz (Steady flow)

Nu( X ) h( X ) d e / λ(5)

式(4)における Tb は混合平均温度である．ここで熱伝達率に関する基礎概念を式(6)および図 4 で表すと，式(7)で示される t ( X ) を温度境界層厚さとした局所熱伝達率が式(8)によって表現できる．

T
Y Y 0

λ

q&

t (X )h( X )

以上により脈動流においてはカラーシュリーレン法
による温度場の可視化結果を基に式(10)によって伝熱特性を評価した．

Nu( X ) pulsatingh( X ) (Tw ( X ) Tb ( X ))(Tw ( X ) Tb ( X ))
T
Y Y 0

λ/

Fig. 5 Temperature field visualized by color schlieren
method for steady flow case (Rem=920).

t (X )Nu( X ) steady H enh. ( X )

(10)

(6)

(7)

(8)

定常流の場合は式(4)および式(8)より温度境界層厚
さ t ( X ) を算出し，その値を可視化結果に照らし合
わ せ るこ とで カラ ーシ ュリ ー レン 画像 にお いてt ( X ) を示す色の境界線を見出した（カラーシュリ
ーレン法の校正作業）．
脈動流においては混合平均温度を一義的に定義できないため，式(4)による熱伝達率の算出は不可能である．そこで脈動流の場合はカラーシュリーレン法
により温度場を可視化し，可視化画像上で上記手順で求めた色境界の厚さを計測（ t ( X ) を直接計測）

以上述べた評価方法においては最終的なヌッセ
ルト数の主たる不確かさ要因は温度境界層厚さの計
測誤差である．実際の温度境界層は 2～3 mm 程度（後述の図 7 参照）であり，計測に使用した定規の最小目盛りが 0.5 mm，相対正確度が 0.0001 であったので，可視化画像を直接計測したのでは計測誤差が非常に大きくなってしまう．そこで計測誤差を極力小さくするために，可視化画像であるカラーシュ
リーレン像を描画系ソフトウェア上で温度境界層を示す色境界が最も際立つよう色調補正した後，拡大印刷した画像上で温度境界層を計測した．
その結果，温度場の発達領域においては温度境界層厚さの計測
誤差は 8%～12%であった．
表 1 に実験条件範囲を示す．平均流速を代表流速とした Rem および脈動周波数は概ね動物の血流の条件とした．振幅につていは実験装置の機構上の制限によるが，流量の最小値は 0 に近づけた（図 3 参照）．
理由は心臓に近い血管内の血流は間欠流に近いこと
による．

し，式(8)より熱伝達率を算出した．また，流れの脈動による伝熱促進／抑制を定量評価するために式(9)に示す局所伝熱促進率を定義した．

2.3 実験結果
図 5 にカラーシュリーレン法において使用したカ
ラースリットの色相と定常流における温度場の可視

H enh. ( X )t ( X ) steady flowt ( X ) pulsating flow

化結果の一例（Rem=920）を示す．可視化結果にお(9)いて加熱平板表面から平板に垂直方向（Y 方向）に
赤い領域の外縁までの長さが図 4 における t ( X ) に相当する．

- 214 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

Rem=920

Rem=3000

0.00 sec

0.20 sec

0.40 sec

0.80 sec

0.60 sec

1.00 sectime
Fig. 6 Fluctuation of thermal boundary layer of one cycle in the case of pulsating flow with its frequency of 1 Hz.

[mm] Re / Re = 0.72fm
3.50
Rem = 920
3.00[mm] Re / Re = 0.72mf
3.50
Re = 3000
m
3.00

2.50

2.50

2.00

2.00
0 Hz

3 Hz

1.00

1 Hz

5 Hz

0.50

2 Hz

0.50

0.40

0.00
120.000.00

1.50

0.00
120.000.00

0.10

0.20

0.30

0.50

1.50
1.00

0.60

100.00100.00

80.00

80.00

60.00

60.00

40.00

40.00

20.00

20.00

0.00
1.600.00

0.10

0.20

0.30

0.40

0.50

0.60

0.00
3.000.00

1.40

2.50

1.20

2.00

1.00

1.50

0.80
H enh. (X) =

0.60
0.40
0.00

0.10

0.20

0.30

2 Hz

3 Hz

0.20

0.30

0.40

0.50

0.60

0.10

0.20

0.30

0.40

0.50

0.60

H enh. (X) =

0.60

0.00
0.00

5 Hz

0.10

0.50

δt'(X)pulsating
0.50

1 Hz

1.00

δt'(X)steady

0.40

0 Hz

0.10

0.20

0.30

δt'(X)steadyδt'(X)pulsating

0.40

0.50

X / XLX / XL(a)(b)
Fig. 7 Local distributions of t (thickness of thermal boundary layer), Nu (Nusselt number) and
Henh. (index of heat transfer enhancement).

- 215 -C 2007 The Heat Transfer Society of Japan○

0.60

Thermal Science & Engineering Vol.15 No.4 (2007)

図 6 は脈動周波数 1 Hz の脈動流において，ある時（Henh.= 1.5～2.5）が得られた．以上により図 6 にお刻から 1 秒間（つまり 1 周期分）における温度境界

いて予測された脈動流場の伝熱特性が定量的に確認

層の変動の様子を示している．図の上段と下段を比較すると Rem の上昇に伴い総じて t ( X ) が小さいことが認められる．また温度境界層の時間に対する変された．
このような伝熱特性がどのような流れ場に依存するのかを調べるために流れの可視化実験を行った．

動に着目すると，同一周波数であっても Rem=3000

３の場合の方が変動幅が大きい様に見受けられる．このことより伝熱特性の脈動周波数依存度は時間平均流量によって異なることが予測される．この予測を流れの可視化実験

3.1 実験装置および方法
脈動流場における流速変動に対する定量データ検証するために前節で説明した方法により，本実験において定義している温度境界層厚さ t ( X ) を計測の取得を目的として PIV 解析を行うために試験流体して脈動周波数が伝熱特性に及ぼす影響を定量評価用したトレーサー粒子は比重調整されたイオン交換した．

樹脂であり平均粒径は 200 m である．図 8 に実験

図 7 に温度境界層厚さ, ヌッセルト数および伝熱促進率の流れ方向局所分布（ t ( X ) , Nu(X)および装置概略図を示す．水（試験流体）はヘッドタンク

Henh.(X)）を示す．局所温度境界層厚さ t ( X ) につい

時間平均流量はヘッドタンク後およびロータフローては時間平均値である（1 秒間に撮影された 16 枚の画像に対してそれぞれ t ( X ) を計測し，その算術平メーター前のバルブの開度によって設定できるよう均値を時間平均値とした）．同図(a)および(b)はそれものと基本構造は同じである．テストセクションはぞれ Rem=920 および Rem=3000 の場合に対応してい寸法も伝熱実験のものと同一であるが，流れの可視を水として粒子懸濁法により流れを可視化した．使とオーバーフロータンクのヘッド差のみで駆動され，になっている．脈動発生装置は伝熱実験で使用した

る．横軸はいずれも流れ方向位置 X を伝熱面長さ XLで除した無次元流路位置を表している（原点 X=0 は伝熱面開始位置）．Rem=920 の場合（図 7(a)），いず
れの周波数に対してもヌッセルト数は下流にいくにしたがい単調減少（温度境界層厚さは単調増加）している．また，流路入口近傍（X/XL<0.25）で流れの脈動によって 20%～40％の伝熱促進（Henh.= 1.2～
1.4）の傾向が認められるが，流路中間より下流の
X/XL>0.5 の領域では伝熱促進の効果は数％のオーダ
ーにまで低下している．しかし，脈動周波数の変化による t ( X ) および Nu(X)の変化はほとんど認められず，伝熱特性の脈動周波数依存度は小さいといえる．
それに対して Rem=3000 の場合（図 7(b)）は，伝熱特性が Rem=920 の場合とは全く異なっている．
t ( X ) および Nu(X)は X/XL に対して単調な増加／減

Fig. 8 Experimental setup for flow visualization.
45

少傾向を示さず，流れ方向領域 0.15<X/XL<0.25 付近
でそれぞれ極大値および極小値をとる傾向を示して

0.208

15

て伝熱が抑制されていることが確認できる．これら

5の特徴は脈動周波数に関わらず認められるが，周波-5数が大きい程顕著である．また伝熱促進率の脈動周-15波数依存度が大きく，周波数が大きい程，伝熱促進

-25

0.375
0.542
X
YW Flat Plate
Side WallY
0の傾向が顕著に示されており，周波数が 3Hz, 5Hz の場合は X/XL>0.3 の領域で 50％～150%の伝熱促進

W

0.042

um = 10.9 mm/s X/XL= 0.3

25いる．更に X/XL<0.10 の領域では流れの脈動によっY/YRef / Rem = 0.72 Rem = 920, f =2Hz

35

0.5

1

1.5

0.708
0.875

2

time [s]Fig. 9 Temporal fluctuations of local velocities at X/XL=0.3.

- 216 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

化用に内部の平板を含め全て透明アクリル製である．

してはレーザーライトシートを流路奥行き方向にト

光源には LD 励起 YAG/YVO4 固体レーザー（出力：

ラバースし，可視化画像に有意差がないことを確認

30 mW, 波長：532 nm, 緑色）を用い，シリンドリカすることで流れ場の二次元性（流路奥行き方向の一ルレンズによってシート状に拡げ，テストセクショ次元性）を確認している．

ンの中央断面における 2 次元速度場を可視化し，デPIV による画像処理には市販ソフトウェア(DIPPィジタルビデオカメラにて撮影した．尚，実験に際

Flow : ㈱ディテクト社製)を使用した．可視化領域にMean flow rate

Flow ratetime
15.0 [mm/s]

7.50

0.00

X

Y

X

Y

Flat plate

Side wall

Flat plate

Side wall

Flat plate

Side wall

(c)
(a)
(b)
Fig. 10 Temporal and spatial periodic fluctuation of flow field in forms of velocity vector and scalar(Rem=920, Ref /Rem=0.72, 2 Hz).
- 217 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

おけるセル毎のベクトル算出は，一般的な時系列画に流量が時間と共に連続的に正弦波変動すると，速像の相互相関関数に基づく確率統計処理によるもの

度場は時間のみならず空間的にも複雑に変動するこである．図 9 に局所流速（流路位置 X/XL=0.3 においとが明らかとなった．この流れ構造は実験したいずて，平板から側壁までの流路幅方向 6 点の流れ方向

れの周波数においても確認された．しかし，Rem=920速度）の時間変化を示す．図 9 と図 3 を照らし合わの場合，図 9 および図 10 の速度スカラー3D 表示図

せると設計コンセプトどおり，流量は連続的正弦波

において平板近傍の流速の大きさに着目するとほぼ変動していることが認められる（2 Hz の条件で 2 秒

流量変動と同期した変動の様相を示しており，時間間に 4 サイクル分の波形が得られている）．

平均すると定常流における速度と大差のない状態である．図 7(a)と照らし合わせると僅かな伝熱促進な

3.2 PIV 解析結果および伝熱特性に関する考察
図 10 および図 11 にそれぞれ Rem=920, 2Hz およびらびに伝熱特性の周波数依存度が小さいという結果

Rem=3000, 3Hz の実験条件における PIV 解析結果をる．

は上記のような流れ構造に起因するものと考えられ示す．両図とも上段は主流方向速度スカラーの 3D

それに対し Rem=3000 の場合（図 11），減速期に側

表示図，下段は主流方向速度ベクトル図である．ベ壁近傍から流れに乱れが生じ始めるという特徴はクトルに関しては，その長さまたは色でスカラー情

Rem=920 の場合同様に認められるものの，流れ方向

報も表現するのが一般的であるが，非定常現象でもに連なる渦列の形成は明確には確認できない（左かあり解析画像一枚あたりの情報量があまりにも多い

ら 2 番目の図，その拡大図は図 11(a)）．また，流量ため，流れ場をより明確に表現するためにベクトル

変動１周期を通して流れに乱れが生じる領域は側壁

図は速度の方向のみを示し，その大きさであるスカ近傍に限定されており，平板近傍の流路幅中央付近

ラー情報は別に 3D 表示とした．

は常に下流へ向かう整然とした流れが存在している．

図 10 および図 11 は共に１周期分の解析結果であ側壁付近においては乱れから逆流へと移行することり，それぞれ左端の図より中央の図に向かって流量が中央の図（拡大図は図 11(b)）より確認できる．時（平均流速）減少期，中央の図より右端の図に向か間平均的には一定の流量がテストセクション内に供って流量（平均流速）増加期に対応している．

給されているため，このような逆流域の存在によりRem=920 の場合（図 10），流れは時間的に変動す

流路中央付近（つまり平板周り）の局所的流量は定ると同時に空間的にも変動し，極めて複雑な構造を常流の場合よりも増加しているはずである．速度ス

とることが確認できる．流路下流に向かう整然としカラーの 3D 表示においては横断面の形状がその Xた流れ（左端の図）が，減速期において側壁近傍で座 標 にお ける 瞬間 の速 度分 布 を示 すわ けだ が，

不安定となり，乱れが生じ始め，その乱れは流れ方

Rem=3000 の場合は流量の減少期においても上記の向全域にわたっている（左より 2 番目の図）．図 10(a)ような平板周りの局所的増速が得られるため，平板はその拡大図である．時間的に流速が変動すること

表面においては常に大きな速度勾配が生じているこ

で流線間にせん断力が生じ渦が形成されるものと考とが確認できる．このような流れ構造により，図 7(b)えられる．更にその渦列は時間とともに流路幅方向で示したような流れの脈動による 50％～150%ものに両側壁から中央の平板に向かって移動する．渦列

伝熱促進が得られたものと考えられる．図 7(b)のが平板にまで移動した時点（流量が最も小さくなる

3Hz の条件（◆）で X/XL=0.3 で約 50%の伝熱促進（Nu

時期とほぼ一致）では側壁近傍で逆流が認められるが定常流の場合の 1.5 倍）を示しているが，図 11 の（中央の図）．図 10(b)にその拡大図を示す．その後，

主流方向速度スカラー値で解析したところ，同位置

増速期に入ると流路中央の平板付近から下流に向かでの局所的流速増加率は 70%（局所流速が時間平均

う整然とした流れが再度形成され始めると同時に両流速の 1.7 倍）であった．局所流速の増加が伝熱促

側壁近傍では逆流から渦列が形成され（右から 2 番進に寄与したことは定量的にも説明されたわけであ目の図），
見かけ上は減速期における側壁近傍の渦列

るが，Rem=3000 という条件は管内流としては乱流遷

形成時と類似した流速分布になっている．図 10(c)移域に属するので，局所流速の増加により境界層がにその拡大図を示す．その後，側壁近傍における渦乱流遷移したとも考えられる（円管内乱流における列が崩壊し，再度流路全域にわたって下流に向かう

Colburn の経験式によれば Nu は Re の 0.8 乗に比例す整然とした流れとなって元の状態に戻る．このようるため流速（Re）が 1.7 倍になれば Nu は 1.70.8=1.53- 218 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

Mean flow rate

Flow ratetime
80.0 [mm/s]

40.0

0.00

X

Y

X

Y

Flat plate

Side wall

Flat plate

Side wall

(a)

(b)

Fig. 11 Temporal and spatial periodic fluctuation of flow field in forms of velocity vector and scalar(Rem=3000, Ref /Rem=0.72, 3 Hz).
倍となり得られた結果と一致する）．しかし，今回得また，側壁近傍の逆流と平板近傍の順流により流られた実験結果から流れが層流であるか乱流である

路内には平板を境として大規模な一対の渦が形成さかを議論するにはデータ不足であるので，境界層のれ，それにより流路入口付近では主流方向流速の Y

乱流遷移という観点からの考察は今後の課題とした方向（流路幅方向）分布が定常流の場合と大きく異い．

なり，その結果として入口付近では平板近傍の主流- 219 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

速度成分は総じて定常流の場合よりも小さくなるものと予想される．Rem=3000 の場合に流路入口付近で伝熱抑制の傾向が得られたのは，このような局所的減速によるものと考えられる．
４

結

参考文献

論

以上の実験結果ならびに解析結果より，矩形管内
脈動流の流れ構造とそれに伴う管内に設置された平板面の伝熱特性に関し以下に示す結論を得た．
伝熱特性
総じて流れの脈動によって平板面の伝熱は促進されるが，局所的には伝熱が抑制される場
合もある．
時間平均流量（Rem）がある値を超えると局所
熱伝達率（ヌッセルト数）は流れ方向に単調減少の傾向を示さず，流路入口近傍で一旦極小値を示した後増加する傾向を示す．
時間平均局所熱伝達率（ヌッセルト数）の脈動周波数依存度は時間平均流量によって異なり，時間平均流量が大きいほど脈動周波数によって伝熱促進率は大きく異なる．
上記流れの脈動による伝熱促進/抑制の特徴は以下に示す複雑な流れ構造に伴う局所的増速/減速に起因する．
脈動流の流れ構造
流量が時間的に正弦波変動すると速度場は時間的に変動するだけでなく，局所的に渦列形成と崩壊を繰り返し，空間的にも複雑に変動する．
流量の時間的変動に関し，流量減少期に側壁近傍で乱れが生じ流れ方向に連なる渦列が形成される．
ある時間平均流量以下（本研究ではRem<1000）では，側壁近傍で形成された渦列は時間の経過とともに流路幅方向内側に向かって移動する．時間平均流量の増加に伴い，渦を含め流れに乱れが生じる領域は側壁近傍に限定され中央の平板面上では常に大きな速度勾配が生じる．
謝によるものである．
ここに記すと共に謝意を表する．
また実験実施およびデータ整理に協力してくれた学部 4 年生大麻 慎哉君ならびに德桝 亮君の両名に感謝申し上げる．

辞本研究遂行に要した経費（設備備品および消耗品の購入費用）の一部は独立行政法人日本学術振興会
科学研究費補助金（若手研究(B)課題番号：15760139）[1] Niida, T., Yoshida, T., Yamashita, R. and Nakayama,
S. “The influence of Pulsation on Laminar Heat
Transfer in Pipes”, Heat Transfer Japan. Res. 3-3(1974), 19-28.
[2] Bayley, F. J., Edwards, P. A. and Singh, P. P. “The
Effect of Flow Pulsation on Heat Transfer by Forced
Convection from a Flat Plate”, Proc. of the FirstInternational Heat Transfer Conference, Boulder(1961), 499-509.
[3] Fujita, N. and Tsubouchi, T., “An Experimental
Study of Unsteady Heat Transfer from a Flat Plate toan Oscillating Air Flow”, Heat Transfer Japan. Res.
11 (1982), 31-43.
[4] Miller, J. A., “Heat Transfer in the Oscillating
Turbulent Boundary Layer”, Journal of Eng. Power,
91 (1969), 239-244.
[5] Jackson, T. W. and Purdy, K. R., “Resonant
Pulsating Flow and Convective Heat Transfer”,
ASME. Journal of Heat Transfer, 87 (1965),507-512.
[6] Kearney, S. P., Jacobi, A. M. and Lucht, R. P.,“Time-Resolved Thermal Boundary-Layer Structure
in a Pulsatile Reversing Channel Flow”, Trans.
ASME. Journal of Heat Transfer, 123-4 (2001),
655-664.
[7] Hessami, M. A., Berryman, A. and Bandopdhayay,P., “HEAT TRANSFER ENHANCEMENT IN AN
ELECTRICALLY HEATED HORIZONTAL PIPE
DUE TO FLOW PULSATION”, ASME. HTD.,
374-3 (2003), 49-56.
[8] Rocha, L. A. O. and Bejan, A., “Geometric
Optimization of Periodic Flow and Heat Transfer ina Volume Cooled by Parallel Tubes”, Trans. ASME.
Journal of Heat Transfer, 123-2 (2001), 233-239.
[9] Li, P., and Yang, K. T., “Fluid mechanics and heattransfer in zero-mean oscillatory flows in shortchannels”, ASME. HTD., 361-1 (1998), 69-76.
[10] Myrum, T. A., Acharya, S. and Inamdar, S.,“Laminar Pulsed Forced and Mixed Convection in a
Vertical Isothermal Tube”, Journal of Thermophys.
Heat Transfer, 11-3 (1997), 423-428.
[11] Mladin, E. C. and Zumbrunnen, D. A., “Dependenceof convective heat transfer to a pulsating stagnation
flow on incident velocity pulse characteristics”,
ASME. HTD., 271 (1994), 11-21.
[12] Zhang, J. G. and Kurzweg, U. H., “Numerical
Simulation of Time-Dependent Heat Transfer inOscillating Pipe Flow”, Journal of Thermophys.

- 220 -C 2007 The Heat Transfer Society of Japan○

Thermal Science & Engineering Vol.15 No.4 (2007)

Heat Transfer, 5-3 (1991), 401-406.
[13] Iguchi, M. and Miura, H., “Radial Propagation of
Turbulence in Pulsatile Pipe Flow”, JSMEInternational Journal, Ser. II, 32-2 (1989), 173-179.
[14] Siegel, R., “Influence of oscillation-induced
diffusion on heat transfer in a uniformly heatedchannel”, Trans. ASME. Journal of Heat Transfer,109-1 (1987), 244-247.
[15] Liao, N. S., Wang, C. C. and Hong, J. T., “Aninvestigation of heat transfer in pulsating turbulent
pipe flow”, ASME. HTD., 42 (1985), 53-60.

[16] Ohmi, M., Usui, T., Tanaka, O. and Toyama, M.,“Pressure and Velocity Distributions in Pulsating
Turbulent Pipe Flow, Part 2 Experimental
Investigations”, Bulletin of the JSME, 19-134 (1976),
951-957.
[17] Mizushina, T., Maruyama, T. and Hirasawa, H.,“STRUCTURE OF THE TURBULENCE IN
PULSATING PIPE FLOWS”, Journal of Chemical
Engineering of Japan, 8-3 (1974)[18] Gerrard, J. H., “An experimental investigation ofpulsating turbulent water flow in a tube”, Journal of
Fluid Mech., 46 (1970), 43-64.

- 221 -C 2007 The Heat Transfer Society of Japan○