ながれ 21（2002） 409409418.

〔特集〕流れを作る：風洞実験

乱流風洞
Turbulence Wind Tunnel∗

豊橋技科大・工学部

蒔田秀治†

Hideharu MAKITA

1

はじめにる乱流発生の最も有効な手段として用いられることとなった．さらに，Dryden2），Batchelor と乱流は流体機械の内部流，自動車や航空機周り

Townsend3），Comte-Bellot と Corrsin4）らによっの流れ，大気，河川，海洋などありとあらゆる流て，主として熱線流速計を用いた格子乱流場の入れ場に現れ，そこで決定的な役割を演じている．

念な計測が行われ，乱流エネルギーの減衰率やエ従って，それらの場にみられる複雑な乱流諸現象

ネルギースペクトル，速度変動の相関など乱流場を忠実に再現し，その観察・計測を通して，乱流の基本特性が求められた．現在の精緻な乱流理論の果たす本質的な役割を理解する必要がある．その体系はそれらの成果を基にして構築されたもののために，目的に応じた様々なタイプの乱流風洞である．図 1 に示すように，乱流格子は丸（角）や，それらに用いる付属装置が開発されてきた．

棒を開口比 55〜60 %程度で格子状に組み合わせ

ここでは，それらのうち代表的なものについて概説する．

2

一様等方性乱流場の形成
一般的に，乱流とは Re 数の大きな流れ場に表Uれる不規則な流体運動であると定義される．乱流場には噴流や境界層，後流などのいくつかの基本M的形態があるが，その中で最も基本的なものは一様等方性乱流場であり，現在の華麗な乱流理論の体系はその概念の上に構築されている．従って，一様等方性乱流の理論体系を論証することは極めて重要であり，風洞中に一様準等方性乱流場を作り，その特性を実験的に検定することが最初に試みられた．風洞内に一様等方性乱流場を作ることは勿論できないが，それに最も近い乱流場として，いわゆる格子乱流が用いられてきた．初めて乱流格子を用いて一様準等方性乱流場を形成し（1934）である．そのたのは，Simons と Salter1）
後，現在に至るまで，乱流格子は風洞実験におけ∗ 〒441-8580 豊橋市天伯町雲雀ヶ丘 1-1†

E-mail: makita@mech.tut.ac.jp

図1

乱流格子とその下流の流れ

410

乱流風洞ただけの極めて簡単な構造であるにもかかわらた．例えば，Gagne ら 6）は，格子乱流場ではなず，その下流に良好な性質を持つ乱流場を誘起でいが，ONERA にある超大型風洞 S1 のリターンきるという特徴を持つ．しかし，この様な格子乱チャンネル（直径 24 m，長さ 150 m，平均流速流場では大きな乱れ強度が得られないこと，下流に向かって乱流エネルギーが減衰すること，その20 m/s）において，Rλ ∼ 2720 に達し，極めて広い慣性小領域（図 12 ）を持つ乱流場を形成し，過程で −5/3 乗則を満たすような明確な慣性小領

その中で速度場の構造関数や確率密度分布等を域が存在する自己保存的なエネルギースペクトル求めている．Kuznetsov7）らも超大型風洞の楕円を持ち得ないことなど，幾つかの本質的な欠点を形ノズル（長軸 24 m，短軸 14 m）からの噴流場持っている．

を熱線流速計で計測し，Rλ ∼ 1660 を得るととも一般的に，格子乱流場の規模は格子間隔 M と流に，乱流場の微細構造の解析を行っている．一方，

速 U に基づく格子レイノルズ数 R M や，マイクロ1/2
スケール λg と特性速度 uc = q2 /3 （q2：乱流エネルギー）に基づく乱流レイノルズ数 Rλ など

Saddoughi ら 8）は NASA Ames Research Centerで表される．例えば，大気乱流等の実験的シミュ実現し，その中で局所等方性の検証を行った．まレーションに用いた場合，少なくとも R M ∼ 107た，大型衝撃波管（直径 30.5 cm）を用いた格子（Rλ ∼ 10 ）が必要であるのに対して，通常の実乱流場 9）（Rλ ∼ 700）の実験なども行われ，格子験室規模の風洞を使用した格子乱流の Rλ 数はた乱流に及ぼす圧縮性（マッハ数 0.3〜0.6）の効果

かだか R M ∼ 10 （Rλ ∼ 10 ）程度であり，後述が調べられている．しかし，上述したような大型（図 12 ）するように，そのスペクトル中に乱流場風洞や特殊な装置を使用して実験を行うことは，の速度変動成分が最も本質的な働きをする慣性小経費や使用上の可能性，あるいは，計測上の問題領域を実現できないことから，通常規模の風洞中点等からも種々な困難を伴うことは言うまでもに実現した格子乱流場は必ずしも Re 数が大きいない．

352という乱流場本来の要請に答えるものとは言い難い．の超大型風洞の測定胴天井に乱流境界層を発達させることにより，Rλ ∼ 1450 のせん断乱流場を一方で，通常の実験室規模の風洞を用い，その中に大きな乱れ強度を持つ乱流場を実現する試Kolmogorov の普遍平衡理論や散逸構造に関しみもなされてきた．小型風洞内で強い乱れを持って乱流場の統計的性質の再検討が求められ，そた一様準等方的な乱流場の実現に成功した例とこで，大きな Re 数を有する乱流場を実験的に作して噴流格子 10, 11）（図 2 ）がある．それは，乱り出す研究が行われてきた．例えば，Kistler と流格子を構成するパイプに多数の小さな孔をあ5）

Vrebalovich は大型加圧風洞中（4 atm）に乱流格子を設置し，R M = 2.4 × 106 ，Rλ = 669 の乱け，そこから主流中に噴流を噴出すことにより，流場を作り出した．その結果，スペクトル中に二を誘起する手法である．例えば，Gad-el-Hak と桁以上の波数帯域にわたる十分に広い慣性小領域（図 12 ）を持つ乱流場を風洞中に実現するこCorrsin10）は，回流風洞の測定胴（断面 1×1.3 m2 ，
長さ 10 m）上流に噴流格子を挿入し，その噴きとに初めて成功した．しかし，気流を加圧してい出し条件（上流あるいは下流に噴出）が乱流場

乱流格子と噴流の相乗効果で乱流場に強い乱れるために散逸スケールも小さくなり，この様な乱の特性に与える効果を調べた．その結果，最大流場では，その小規模構造に関する議論をする場
合，流速を検知するプローブの分解能など計測上

Rλ = 160 の準等方性乱流場（u/v ∼ 1.11）が形成された．しかし，彼らによって達成された Rλの問題が生じる．

ではエネルギースペクトル中に明確な慣性小領域

近年，各地で大型風洞が建設され，乱流場の基礎的研究にも，それらの利用が可能となってきは存在せず，十分に大規模な乱流場が実現されているとは言えなかった．

411

蒔田秀治

gridjetU（a）Gad-el Hak & Corrsin（1974）10）（a）Ling & Wan（1972）12）(㬍20Ბ)（b）Tassa & Kamotani（1975）11）

図 2 噴流格子（b）斉藤・佐藤（1974）13）

図3

水流・気流攪拌装置

それとは別に，動的な攪拌装置を用いて乱れ誘起することを試みた．すなわち，気流中の乱れを誘起することも各種提案されてきた．Ling とを取り出し，そこに含まれる超低周波（1 Hz 以Wan12）の実験は水槽中で行われたものであるが，下）のゆらぎを重畳させた準周期的な振動を与えその様な研究例の最初の一つである．彼らは，垂ることによって，約 9 %の乱れ強度を持つ乱流場直に並んだ回転軸に取りつけられた長方形の攪拌を実現した．しかし，エネルギースペクトル中に板（図 3（a））を主流に対して ±19 ˚ の範囲で振は励起周波数成分が残存し，必ずしも良好な結果動させ，開水路（測定胴断面 0.6 × 0.6 m2 ，長さを得ることはできなかった．ただし，この様に流11 m）中に乱流変動を誘起したが，必ずしも大きな乱流レイノルズ数（Rλ ∼ 60）は得られなかっれ場を動的に制御して乱れの強い乱流場を誘起した．同様の手法は風洞中でも試みられている．斎

一様性や等方性を有するなど，良質な大規模乱流

藤と佐藤13）は，図 3（b）に示すような正弦波状の振動板（長さ 1 m, 幅 5 cm）20 段からなる気流
攪拌装置を試作し，それを風洞の測定胴上流に設置した．本装置の特徴は各回転軸端に取り付けたようとする手法は，乱流レイノルズ数が大きく，場を通常規模の風洞中に実現し得る可能性を示したものとして高く評価される．

3 せん断乱流場の形成

ステッピングモータにより振動板の駆動モードを現在，気象や大気拡散，建物周りの流れなど風制御可能としたことである．彼らは大気乱流と同環境工学分野における実験的研究が盛んに行われ様なスペクトル分布を有する乱流場を風洞中にるようになった．この種の実験を行うため環境・

412

乱流風洞

図 4 気象風洞（コロラド州立大学）14）䉴䊌䉟䊟䊷plates forming channel☻ᐲscreen slots
generator
図6

速度成層装置（Champagne el at.，197018））

上述した様な風環境問題など，せん断乱流場中
図 5 風洞中に設置されたスパイヤーと粗度の一例 16）で行われるシミュレーション実験でまず第一に重要なことは，主流に正確な平均速度分布を与え気象風洞と呼ばれる大型の装置が建設されていることである．そのためには，簡単な構造を持つる．図 4 に示すコロラド州立大学の大型気象風

不等間隔格子 17）が最初に用いられた．これは水（1971 年）はその最初の例である．この気洞 14）

平に並んだ二次元的な格子の間隔や太さを高さ

象風洞は，その中で接地境界層を模擬するために方向に不均一に変えることによって局所的な抵長大な測定胴（長さ 29 m）を持ち，流れ方向の圧抗係数を変化させて平均流に速度分布を与える力勾配を調整するための可変天井（2.0 ± 0.5 m）ものである．それを一歩進めたものに速度成層装の採用や，温度成層流を形成するために測定胴

置 18）（図 6 ）がある．速度成層装置は流路を数段床面の温度制御装置（加熱・冷却）等を備えていに分割し，それらに開口比の異なる金網や多孔板る．以来，この様な環境・気象風洞を用いて，都を挿入することによって各段を通る流れに異なっ市の風環境や汚染物質の拡散現象等の解明を目的た抵抗値を与え，その下流に任意の速度分布を形とした多くの実験が行われてきた．しかし，この成する装置である．このような装置を用いて目的様な大型風洞を用いて大気境界層を模擬するためとする平均速度分布を正確に形成するためには，には，単に平均速度分布のみならず，乱流強度や各段毎に様々な開口比の金網や多孔板を用意し，レイノルズ応力等の乱流特性の分布を実際の接地速度分布を何度も確認しながら，その抵抗値を繰境界層と相似なものにする必要がある．実際の風り返し調整し直す必要がある．それを改良したも環境実験では，より精密な大規模せん断乱流場をに示す可変格子を用いた速度成のとして，（図 7 ）実現するために，様々な装置上の工夫が凝らされ層装置 19）がある．それは 4 枚の素子を動かすこている．例えば，図 5 にみられるように，測定胴とにより各段の開口比を微妙に調整することがで上流にスパイヤー 15）と呼ばれる特殊な三角錐のき，種々の速度分布を容易に実現することが可能突起を設置したり，種々な粗度15）を測定胴床面上に敷き詰めることによってその実現が図られている．となっている．
しかし，これらの手法では平均流速分布を正確に模擬することはできても，大気乱流場の実験的

413

蒔田秀治

シミュレーション等を行うために要求される強いが可能な多数のファンを用いた乱流風洞を開発し乱流強度や大きな渦スケールを伴った大規模せんた（図 8 ）．この乱流風洞には直径 0.18 m の AC

断乱流場を風洞中に実現することは極めて困難でファンを縦に 11 個，横に 9 列配置してあり，さらある．すなわち，平均速度分布に加えて乱れ強度に，測定胴下流に 21 枚の振動翼列（NACA0010）や渦スケール，レイノルズ応力など乱流諸量の鉛が設置されている．これらのファンと振動翼の直方向分布を制御し得るものでなければならない

動作モードをコンピュータ制御することによっが，実際，これらの静的な装置のみでは不可能でて，図 9 に示すように，平均速度 U ，乱れ強度あった．

Iu ，インテグラルスケール Lu 等の鉛直方向分布

このような問題を解決するために，流れ場にを風洞中に実現している．特に注目すべきは，実能動的な制御を加えて大規模せん断乱流場を実地観測した流速の瞬間波形（図 10 上段）を基に現し，同時にその中における乱流特性量の分布を制御しようとする試みがなされている．例えば，Nishi ら 20）はサーボモータにより風量制御

図 9 乱流風洞（宮崎大）によって実現された平均
図7

速度，乱れ強度，インテグラルスケールの可変格子を用いた速度成層装置鉛直分布 20）（佐藤・恩田，197319））㩿㪸㪀ታ࿾᷹ⷰߦࠃࠆᵄᒻ
㩿㪹㪀C ࠍၮߦౣ⃻ߐࠇߚᵄᒻ
図 10
図 8 多段ファンを用いた乱流風洞
（宮崎大）

乱流風洞
（宮崎大）
によって再現された速度変動の波形 20）

414

乱流風洞

して，AC ファンと振動翼の動作をフィードバッ全断面積（0.7 × 0.7 m2 ）にわたって縦横各 20 本ク制御することにより風洞中にそれを実現した結の回転軸（格子間隔 M = 35 mm）を格子状に組果（図 10 下段）であり，本風洞の性能の高さがみ，それぞれの回転軸に取り付けた多数の小さなうかがえる．本風洞の構成上，高い周波数の変動矩形攪拌翼を反転振動させることにより乱流場の（数十 Hz 以上）まで厳密に再生することは難し励起を行ったものである．各回転軸端には外部のいが，突風の影響等，風環境に関する実験的な研ステッピングモータが直結され，それらは個別に究分野では高い利便性を発揮する装置であると思駆動される．各回転軸の動作モードはワンボードわれる．

コンピュータ，および，それらに制御条件を与えは，気流をランダムに攪拌して，スるパーソナルコンピュータにより平均迎角 α とペクトルの低周波領域にエネルギーを与えるこ

回転範囲 β が実験目的に従って設定され，その条とにより，実験室規模の風洞中に大きな速度変動件下で攪拌翼をランダムに反転振動させている．

を持つ発達した乱流場を作るという発想から，多本装置では回転速度や各軸の α と β を変化させ数の攪拌翼を有する回転軸を格子状に組んだ乱ることにより，一様準等方性乱流，一様せん断乱流発生装置（図 11 ）を開発した．それは流路の流，乱流境界層など任意の流速分布や乱流特性を21）

Makita㫊㫋㪼㫇㫇㫀㫅㪾㩷㫄㫆㫋㫆㫉㫅㫆㫑㫑㫃㪼㪽㫃㫆㫎

有する乱流場を下流の測定胴内に形成することがd㪌 㪉㪅㪌㪻㫉㫀㫍㪼㫉㪸㪾㫀㫋㪸㫋㫆㫉㩷㫎㫀㫅㪾できる．
本装置の特色は，スペクトルの低周波側（約㪊㪌㪟 㪔㪎㪇㪇20 Hz 以下）に直接ランダムな三次元的撹乱を与㪊㪌㪽㫃㫆㫎㪉㪇㩷㫃㪸㫐㪼㫉㫊㪺㫆㫄㫇㫌㫋㪼㫉㫇㪼㫉㪽㫆㫉㪸㫋㪼㪻㩷㫇㫃㪸㫋㪼㫌㫅㫀㫋㩷㪑㩷㫄㫄 㫊㪿㪼㪸㫉㩷㪽㫃㫆㫎㩷㪾㪼㫅㪼㫉㪸㫋㫆㫉㪊㪣㪔㪍㬍㪈㪇㫋㪼㫊㫋㩷㫊㪼㪺㫋㫀㫆㫅

図 12エネルギースペクトルの比較●：通常の乱流格子 21）
R M = 1.56 × 104 ，Rλ = 25▲：大型加圧風洞（Kistler et al 5））
R M = 2.4 × 106 ，Rλ = 669○：動的乱流格子（Makita 21））R M = 1.56 × 104 ，Rλ = 390▽：超大型風洞（Saddoughi et al 8）），Rλ = 1450■：超大型風洞（Gagne et al 6）），Rλ = 2720
図 11 大気乱流風洞およびせん断乱流発生装置（豊橋技術科学大学）25）△：潮流（Grant et al 22））
，Re = 108 ，Rλ > 2000□：大気乱流（佐々ら 23）），Rλ = 2010

415

蒔田秀治
表 1 乱流特性量の比較U0MX/MRMu /U u /v
Luxλg
RL
Experiment［m/s］ ［mm］
（103 ） ［%］［mm］ ［mm］
Turbulence grid 4）
1050.8
4234.02.22
1.024
4.9355
Turbulence grid 5）
61
171.0
45
2400
1.97
1.23
48
2.4
13300
Jet gird 10）
4.2
102.0
46
41.8
3.27
1.11
134
9.45
2280
Active gird 12）
0.003
64.0
45
0.94
4.24
1.05
59
26.7
75
Active gird 21）
5
46.7
50
15.6
16.4
1.22
197
8.29126
Atmosperic turbulence 24）
4.3
12.0
144.0

1.2

ケードを念頭において，十分に広い慣性小領域を1Y/ δ持つ自己相似的な乱流場の実現を直接的に目指したものである．結果として，表 1 に示すように，0.80.6

乱れ強度 u /U > 35〜16 %，インテグラルスケー

0.4

ル Lux = 190 mm に達し，Rλ = 390 の大規模な0.2

一様準等方性乱流場が形成された（図 12 ）．
上述した各種装置によって形成された乱流場の特性を比較したものを表 1 に示す．Kistler とVrebalovich5）では，インテグラルスケール Lux に基づいた乱流レイノルズ数 RL や，Rλ は他の乱流場よりも大きいものの，乱れ強度は u /U ∼ 2 %72669
160343904830

X/M=1001.4えたことである．これは乱流場のエネルギーカス

Rλ00

図 130.5U/U

1.000.05 0.1u'/U

00.002 0.004
2-uv/U

せん断乱流発生装置によって形成された
平均速度，乱れ強度，レイノルズ応力の鉛直方向分布 25）
×：Case 1，○：Case 2，−：U ∝ Y 1/7：大気乱流 27）と通常の格子乱流場と同程度であり，マイクロス：発達した乱流境界層 28）

ケール λg も小さい．また，Jet grid10）では，乱さらに，蒔田ら 25）は流路を 20 段に分割したれ強度や Lux ，λg は通常の格子乱流よりは大きい

速度成層装置（図 11 ）を乱流発生装置の上流にものの，高レイノルズ数乱流の条件を十分に満併設し，せん断乱流発生装置とした．その各段にたすものであるとは言い難い．Ling と Wan12）の挿入する金網や多孔板の開口比を変えることに

Active gird では通常の格子乱流場よりも 2 倍程

よって，一様せん断乱流場や境界層厚さ，δ，が大度の乱れ強度を得ているが，水槽実験であること

きい乱流境界層（図 13 ）などをより正確に形成から，十分な乱流レイノルズ数を得られていない．

できる．その結果，δ = 400〜500 mm，Rλ = 80021）

Makita の乱流発生装置で得られた諸特性では，
λg は Jet grid10）や Ling と Wan12）の Active girdに達する大規模な乱流境界層が実現されている．

よりも小さいが，大気乱流に匹敵する大きな乱す．Case 1 は乱流発生装置を静止させ，通常の乱れ強度を実現しており，その結果，Rλ ∼ 390 と流格子として用い，それと速度成層装置を併用しいう比較的大きな値が得られている．また，エネて形成した乱流場，Case 2 は攪拌翼を各段毎に異ルギーを含む渦のスケールとされる Lux も十分になった動作モードで駆動させて形成した乱流場で大きく，図 12 に示すように，エネルギースペクある．両乱流場とも平均速度分布は平坦な地形上トル中に −5/3 乗の勾配を有する慣性小領域が二に発達する 1/7 乗則分布 26）を与えている．両者桁以上にわたって明確に存在している．これはの結果を比較すると，Case 2 では，実際の接地境

Kistler と Vrevalovich の結果である R M > 10

6界層と同様に，外層部にも大きな乱れ強度やレイに匹敵する大規模な乱流場が実現されていること

ノルズ応力が分布しており，壁面近傍におけるレを示している．

イノルズ応力分布を除いて，エラーバーで示した

5）

本乱流境界層の測定結果の一例 25）を図 13 に示

416

乱流風洞

ば，ビル風や煙突からの排煙の拡散等の問題を正Lwx =0.2d

確にシミュレーションするためには平均速度分
Flow

布，乱れ強度やレイノルズ応力分布等を実際の接地境界層内の特性と一致させることに加えて，対象となる流体現象と風洞気流中の乱流渦のスケールの比まで等価的に合致させることが必要であることを示している．

d/㱐 =0.15

Lux =0.7d

また，本大規模乱流場は乱流現象の素過程等の（a）乱流格子+速度成層装置解明にも用いられている．例えば，乱流発生装置の下流にハニカムやスリットを設置して，鉛直Lwx =1d方向あるいは水平方向の速度変動を抑制するこFlowとにより，強い異方性を有する軸対称異方性乱流場 30）や 3 次元異方性乱流場 31）も形成することができる．著者らの研究室では，この様な異方性乱流場の等方化プロセスの解明などを通じて

Rotta のモデルや SSG モデルの検証 30, 31）など乱Lux =27d

流場の基礎過程に関する研究も行っている．さら（b）せん断乱流発生装置

図 14に，風洞乱流における本質的欠陥の一つである乱建物周りのながれ 25）

流エネルギーの減衰が極めて小さい乱流場の実実地観測結果 27）とよく近似した分布が得られて現等の研究へと，テーマを発展させている．さらいる．壁面近傍の乱流特性をより正確に合致させに，Sassa と Makita32）はこの乱流発生装置と温るためには，粗度等を併用する必要があるが，せ度成層装置を組み合わせることにより，高レイノん断乱流発生装置を用いることにより実際の接地

ルズ数安定成層流中の熱輸送について研究を行

境界層の乱流特性がほぼ正確に再現されているこい，Mydlarski と Warhaft33）も本装置の手法を取とがわかる．

り入れて Rλ ∼ 730 の乱流場を実現し，パッシブ

このような乱流場に建物模型を挿入して，水平スカラーの統計的性質を解析している．現在，こに設置（Y/δ = 0.12）したスモークワイヤーによれまでに述べた様な実験的研究を通して，種々のり，周囲流の様子を観察した結果

25）を図 14 に示す．Case 1 の乱流場では，建物前面の角部で流れがはく離し，背後に明確な死水領域が形成され大規模乱流場の構造が詳細に解明されつつある．

4 おわりにている．一方，Case 2 のせん断乱流場では，はく現在，理論的研究や CFD の発達によって得ら離渦が建物背後まで大きく回り込んでいる様子がれた成果を検証するために，理想的な乱流場中で観察できる．これは建物の幅（d = 60 mm）と主乱流の基礎構造をより厳密に計測することが望ま

流中の乱流渦のスパン方向インテグラルスケールれている．風洞実験によって与えられるものは，Lwx が同程度であり，それが建物模型から発生す常に実在する流体現象である．しかし，我々実験る渦構造と相互に干渉し，後流中の乱流構造に強者が風洞実験で直接観察し得るものは，あくまい影響を与えているためである．また，本乱流場でも複数の流体力学的事象が複合されたものでで得られる大きな拡散係数（C ∼ 147 cm /s）をあり，常にあいまいさを伴う．さらに，そこから229）なども行われてい抽出する個々の計測結果には，誤差も含まれる．

るが，これらの結果は大気乱流場の諸現象，例え風洞実験を試みるとき，まず第一に重要なことは利用して，大気拡散の実験

417

蒔田秀治

「良質な流れ場を得ること」である．そのためには，今後も風洞やその付属装置に一層の改良に加generated turbulence, J. Fluid Mech. 432（2001）219.

える必要がある．さらに，計測や可視化の手法，

10） Gad-el-Hak, M. & Corrsin, S. : Measurements of

データ処理の手段等に多くの工夫を加え，乱流場

the nearly isotropic turbulence behind a uniformを記述するのに必要な高次の相関量や統計量につjet grid, J. Fluid Mech. 62（1974）115–143.

いても信頼性の高い実験データを得て，そこから

11） Tassa, Y. & Kamotani, Y. : Experiments on tur-

有意義な情報を正確に抽出しなければならない．

bulence behind a grid with jet injection in down-

その上に立って，理論的研究や数値実験の結果とstream and upstreamdirection, Phys. of Fluids 18

互いに検証し合いながら，複雑な流体現象の真の姿に少しでも接近することが風洞実験の永遠のテーマであることは言うまでもない．

（1975）411.

12） Ling, S. C. & Wan, C. A. : Decay of isotropic turbulence generated by mechanically agitated gird,
Phys. of Fluids 15（1972）1363.

引用文13）斉藤・佐藤 : 大気の乱れを風洞の中に作る一方法,献1） Simons, L. F. & Salter, C. : Experimental investigation and analysis of the velocity variations, Proc.
Roy. Soc. Ser. A 415（1934）212.

14） Cermak, J. E. : Laboratory simulation of the atmospheric boundary layer, AIAA J. 9（1971）1746.

2） Dryden, H. L. : Isotropic turbulence in theory andexperiment, T. von Karman Anniv.（1941）85.
3） Batchelor, G. K. & Townsend, A. A. : Decay of
isotropic turbulence in the initial period,

乱流シンポジウム講演論文集（1974）103.

Proc.

Roy. Soc. A 193（1948）539.

15） Rae, W. H. & Pope, A. : Low-speed wind tunnel
Testing（John Wiley & Sons., 1984）486.
16）村上 : 高層建築物による風害の予測・対策と風洞実（1983）
験手法について, 東京大学生産技術研究所

77.

4） Comte-Bellot, G. & Corrsin, S. : The use of a con-17） Rose, W.G. : Results of an attempt to generate a

traction to improve the isotropy of grid-generated

homogeneous turbulent shear flow, J. Fluid Mech.

turbulence, J. Fluid Mech. 25（1963）657.

25（1966）97.

5） Kistler, A. L. & Vrevalovich, T. : Grid turbulence18） Champagne, F. H.,

Harris, V. G. & Corrsin,

at large Reynolds numbers, J. Fluid Mech. 26

S. : Experiments on nearly homogeneous turbulent（1966）37.

shear flow, J. Fluid Mech. 41（1970）81–139.

6） Gagne, Y., Hopfinger, E. J. & Frish, U. : A New
Universal Scaling for Fully Developed Turbulence,NATO ASI 237（1990）315.

19）佐藤・恩田・斉藤 : 東京大学宇宙航空研究所報告9（1973）675.
20） Nishi, A., Kikugawa, H., Matsuda, Y. & Tashiro,7） Kuznetsov, V. R., Praskovsky, A. A. & Sabelnikov,

D. : Turbulence control in multiple-fun wind tun-V. A. : Fine-Scale Turbulence Structure of Inter-nels, J. Wind Engi and Industrial Aero 67（1997）

mittent Shear Flows, J. Fluid Mech. 243（1992）

861.

595.
8） Saddoughi, S. G. & Veeravalli, S. V. : Local
Isotropy in Turbulence Boundary Layers at HighReynolds Number, J. Fluid Mech. 268（1994）
333.
9） Briassulis, G., Agui, J. H. & Andreopoulos,
Y. : The structure of weakly compressible grid-21） Makita, H. : Realization of a large-scale turbulence
field in a small wind tunnel, Fluid Dyna. Res. 8（1991）53.

22） Grant, H.L., Stewart, R.W. & Moilliet, A. : Turbulence spectra from a tidal channel, J. Fluid Mech.
12（1962）241.
23）佐々・千葉・蒔田 : 接地境界層中の乱流微細構造

418

乱流風洞に及ぼす安定度の影響, 日本流体力学会第 27 回

29）蒔田・佐々・飯田・茂山 : 大気乱流拡散の実験的乱流シンポジウム講演論文集（1995）203.

シミュレーション, 機論 522-B（1990）388.

24） Sheih, C.M., Tennekes, H. & Lumley, J.L. : Air-30）蒔田・岩崎・飯田 : 軸対称異方性乱流場の実験的borne hot-wire measurements of the small-scale
structure of atmospheric turbulence, Phys. of Fluids 14（1971）201.

解析, B54-505（1988）408.

31）蒔田・南 : 三次元異方性乱流場の減衰過程に関する実験的研究, B61-581（1995）26.

25）関下・蒔田・一郷・藤田 : 大気乱流境界層の実験

32） Sassa, K. & Makita, H. : Heat transfer in stable

的シミュレーション, 機論 B 68-665（2002）55.

stratified and high Re turbulence, Internal work-26） Davenport, A.G. : Proc. ASCE 86, Spr 5（1969）

shop on physical modeling of flow and dispersion39.
27） Hanafusa, T. & Fujitani, T. : Meteorology andGeophysics 32（1981）19.
28） Klebanoﬀ, P.S. : NACA TN 3178（1954）.

phenomena（2001）.
33） Mydlarski, L.& Warhaft, Z. : Passive scalar statistics in high-Peclet-number gird turbulence, J. Fluid
Mech. 358（1998）135.