同じダクトから、レイノルズ数が2つ出てきた
400×300mm の空調ダクトに空気を 8m/s で流す。レイノルズ数を出そうとして資料をめくると、「相当直径」という言葉が2か所に出てくる。圧損線図の脇に載っている式で計算すると 377.709mm。非円形断面の代表長さとして教科書に載っている 4A/P で計算すると 342.857mm。どちらも日本語では同じ「相当直径」と呼ばれている。
前者を入れると Re=199,928、後者を入れると Re=181,480。差は10.2%。桁が変わるわけではないので、どちらを使っても「乱流」という結論は動かない。だから気づかないまま次の計算へ進める。そして Re は圧力損失にも熱伝達にも入っていくので、この10%は最後まで消えずに設計値へ乗る。
この2つの「相当直径」は、名前が同じだけで別の量だ。片方はレイノルズ数に入れるための長さ、もう片方は圧損線図から丸ダクトの径を読むための長さ。このツールは両方を並べて出し、Re に使ったのがどちらかを結果に明記する。以下はその食い違いの正体と、円管以外の形状で「長さに何を入れるのか」という、もっと手前にある問題の話。
なぜ作ったのか — 見つかるのは円管専用の計算機ばかりだった
日本語で「レイノルズ数 計算」を探すと、無料の計算ツールがいくつも出てくる。調べた範囲で4つ見つけたが、いずれも円管専用だった。流速・内径・動粘度の3欄に数字を入れると Re が出る。ところが実務でいちばん間違えるのは、その「内径」の欄に何を入れるかのほうだ。矩形ダクトなら。平板なら。球なら。撹拌槽なら。そこは利用者に丸投げされている。
動粘度が手入力なのも気になった。cSt と cP は名前が似ているが、前者は動粘度、後者は粘度で別の量だ。cP の数字を cSt の欄に入れても計算機は何も言わない。単位換算を1回飛ばすだけで Re の桁が6つ動き、層流だったものが乱流に化ける。
決定打は自分のサイトの中にあった。配管圧力損失、ポンプ揚程、オリフィス、強制対流、撹拌動力 —— 6本のツールが内部でレイノルズ数を計算しているのに、Re そのものを出すページが1枚も無かった。検索コンソールで確認すると、「レイノルズ」を含むクエリはサイト全体で90日にたった1回しか表示されていなかった。受け皿が無いのだから当然だ。
だから独自性は3つに絞った。(1) 円管・矩形ダクト・平板・球・撹拌槽の5形状を切り替え、代表長さの取り方と臨界レイノルズ数を同時に入れ替えること。(2) 矩形ダクトでは水力直径 Dh と円形相当直径 De を区別して両方出し、Re に使ったほうを明示すること。(3) 動粘度を物性ライブラリと規格の式から機械生成して、流体と温度を選ぶだけで引けるようにすること。
レイノルズ数とは何か — 慣性力と粘性力の比を第一原理から
レイノルズ数 とは — 力の大小ではなく「比」
水を張った洗面器で手を振ると、指の後ろに小さな渦ができて、手を止めても水はしばらく回り続ける。同じことを蜂蜜でやると渦はできない。手の形にまとわりついて流れ、手を止めた瞬間に蜂蜜も止まる。手の速さも手の大きさも同じなのに、起きていることがまるで違う。
この違いを1つの数にしたのがレイノルズ数だ。
Re = ρ v L / μ = v L / ν
ρ: 密度[kg/m³] / v: 代表速度[m/s] / L: 代表長さ[m]
μ: 粘度[Pa·s] / ν: 動粘度[m²/s](ν = μ / ρ)
分子の ρvL は流体が動き続けようとする勢い(慣性力)、分母の μ は周りに引きずられて止まろうとする効果(粘性力)を表す。洗面器の水は慣性が勝つので渦が残り、蜂蜜は粘性が勝つので渦ができない。つまり Re は「流れが速いか」でも「どろどろか」でもなく、その2つのどちらが優勢かという比を見ている。
無次元だから相似則が成り立つ — レイノルズ数 求め方の要点
比なので単位が消える。[m/s] × [m] ÷ [m²/s] を約分すると何も残らない。この無次元性こそがレイノルズ数の効きどころで、実物と模型が同じ Re なら同じ流れになる(相似則)。船や航空機の風洞試験、化学プラントのスケールアップが成り立つのはこのためだ。
計算する側から見ると、無次元であることは「単位さえ揃えれば何を使ってもいい」という意味になる。よくある事故は動粘度の単位で、カタログの 46 mm²/s をそのまま 46 として式に入れると答えが10⁶倍ずれる。m²/s への換算を先に済ませてから代入する。
代表長さ 決め方 — ここだけは第一原理から一意に決まらない
ρ・v・μ は測れば決まる。決まらないのは L だ。「流れが感じている大きさ」という意味なので、形状ごとに約束が要る。
| 形状 | 代表長さ L | 取り方 |
|---|---|---|
| 円管 | 内径 D | 内径をそのまま使う |
| 矩形ダクト | 水力直径 Dh | Dh = 4A/P = 2ab/(a+b) |
| 平板 | 先端からの距離 x | 位置ごとに Re が変わる |
| 球 | 直径 d | 球・円柱まわりは直径 |
| 撹拌槽 | 翼径 d | 代表速度も n·d に変わる |
目を引くのは平板だ。同じ流れの中で、場所によって Re が変わる。先端では 0、下流へ行くほど大きくなり、どこかで乱流境界層に切り替わる。「この流れの Re は」という言い方そのものが、平板では成り立たない。
撹拌槽も変則で、代表速度に流速ではなく翼まわりの速さ n·d(n は毎秒回転数)を使う。結果として Re = n·d²/ν になり、長さが2乗で入る。rpm を毎秒回転数に直し忘れると60倍ずれる。
臨界レイノルズ数 2300 は何を指しているのか
慣性が粘性に勝ち始めると、流れは層状のまま進めなくなって乱れる。円管でその境目にあたるのが Re=2,300 で、そこから Re=4,000 あたりまでが遷移域、それ以上が乱流というのが日本の実務書と日本語版 Wikipedia の扱いだ。このツールも同じ 2,300 / 4,000 を使う。
ただしこの数字は形状ごとに違う。平板の境界層は Re_x ≈ 5×10⁵、撹拌槽は 10 と 10⁴、球は 1(ストークス域の上限)・1,000(ニュートン域の入口)・3×10⁵(抗力危機)と、バンドの数まで違う。形状を変えたら、長さの取り方と一緒に境目の数字も入れ替わる。ツールで形状を切り替えると数直線のバンド境界が動くのは、そのためだ。
層流か乱流かで何が変わるのか — 圧力損失と熱伝達で安全側が逆になる
管摩擦係数の式が丸ごと入れ替わる
層流の管摩擦係数は f = 64/Re の一本で決まる。管の内面が磨かれていようが錆びていようが関係ない —— 粗さは式に入っていない。ところが乱流になると Colebrook-White の陰関数式になり、相対粗さが効いてくる。境目をまたぐと、使う式そのものが変わる。圧損を出す前に Re を確かめるのは、電卓を叩く順序の問題ではなく、どの式を開くかという分岐の問題だ。
熱伝達係数は数倍変わる
管内の熱伝達はもっと極端だ。層流では Nu が一定値に張りつく(等熱流束の円管で 4.36)。流速を上げても熱伝達係数はほとんど増えない。乱流では Nu が Re の0.8乗で伸びる。冷却が足りないという不具合の原因が「思ったより Re が低くて層流のままだった」であることは珍しくない。流量を増やして解決しようとしたのに効かない、という展開になる。
遷移域は「どちらでもない」ではなく「設計してはいけない」帯
Re=2,300〜4,000 の遷移域は、同じ配管でも日によって層流に寄ったり乱流に寄ったりする。入口の擾乱、曲がりの位置、温度の振れで簡単に行き来する。ここで最も重要な帰結を書く —— 圧力損失を見るなら乱流側、熱伝達を見るなら層流側で評価するのが安全側になる。乱流のほうが損失は大きく、層流のほうが熱は伝わりにくいからだ。つまり目的によって安全側が逆を向く。「遷移域なら乱流側で見ておけばいい」という定番のアドバイスは、圧損では正しく、熱設計では危険側になる。
冷間始動 — 温度だけで判定が変わる
実害が出やすいのは油だ。油圧配管 25A を 4m/s で流すと、ISO VG46 の油温 40℃ では Re=2,400 で遷移域。同じ配管・同じ流速のまま油温を 20℃ に落とすと Re=821 で完全な層流に落ちる(使用例のケース3・4)。流速は1mmも変えていないのに判定が変わる。しかも層流では f = 64/Re なので、Re が下がるほど管摩擦係数は大きくなる。冬の朝だけポンプの起動トルクが足りない、リリーフ弁が吹く、という故障はここから来ている。40℃の物性値だけで検討していると、この経路は最後まで見えない。
こんな場面で使う
圧損を出す前の分岐確認。 配管の圧力損失を計算するとき、最初にやるのは管摩擦係数をどの式で出すかの決定だ。Re を先に出して層流か乱流かを確かめてから、配管圧力損失計算へ渡す。
ダクトの相当直径で迷ったとき。 圧損線図から読んだ相当直径をレイノルズ数に入れてよいのか、という迷いは設備設計でくり返し出てくる。Dh と De が並んで出れば、どちらが何のための長さか一度で分かる。丸ダクト径の選定そのものはダクトサイズ計算ツールが担当する。
撹拌槽のスケールアップ。 試験槽で乱流だったものが、実機で翼径と回転数を変えた結果、遷移域に落ちていることがある。翼径と回転数を入れて Re を比べれば、流動状態が保たれているかどうかが先に分かる。動力そのものは撹拌機 動力計算シミュレーターへ。
伝熱面のどこから乱流になるか。 ヒートシンクのフィンや熱交換器の平板では、先端から何mmで乱流境界層に切り替わるかで熱伝達が変わる。遷移位置を出しておくと、強制対流熱伝達係数 計算ツールでどの相関式を選ぶかの判断がつく。
基本の使い方
① 形状を選ぶ。 円管・矩形ダクト・平板・球・撹拌槽の5つ。ここを選んだ時点で、代表長さの取り方と臨界レイノルズ数が同時に決まる。選択に応じて入力欄そのものが入れ替わるので、使わない寸法を埋める必要はない。
② 流体と温度を入れる。 水・海水(塩分3.5%)・空気・エチレングリコール水溶液30%/50%・ISO VG32/46/68 から選んで温度を入れると、動粘度と密度が入る。実測値やカタログ値があるなら直接入力に切り替えればよく、そのときは温度欄は使わない。
③ 寸法と速度を入れる。 円管なら内径と流速、矩形ダクトなら幅・高さと風速、平板なら先端からの距離と流速、球なら直径と流速。撹拌槽だけは流速ではなく回転数[rpm]と翼径を入れる。撹拌槽で流速欄が出てこないのは仕様で、式に流速が入らないためだ。
結果にはレイノルズ数と流動状態のほか、使った代表長さ・動粘度・密度、そして**臨界レイノルズ数になる流速(撹拌槽は回転数)**が出る。「いくつまで落とせば層流か」が数字で出るので、判定にどれだけ余裕があるかを掴める。あわせて Re の対数数直線に形状ごとのバンドが描かれ、いまの値がどのバンドのどこにいるかを目で追える。
具体的な使用例(7ケース)
以下はすべて、実装した計算エンジンに同じ入力を与えて出力を確認した値。
ケース1: 給水管50A・水20℃・流速2m/s
入力 → 形状=円管 / 流体=水 20℃ / 内径 52.9mm / 流速 2m/s
結果 → Re = 105,441(乱流)。動粘度 ν=1.0034 mm²/s、体積流量 15.82 m³/h。Re=2,300 になる流速は 0.0436 m/s、Re=4,000 になる流速は 0.0759 m/s
解釈 → いちばんありふれた条件で、余裕をもって乱流。層流に落とすには流速を 0.0436 m/s まで、設計流速 2m/s から桁を2つ近く下げなければならない。水の配管を扱っている限り層流を気にする場面はほぼ無い、ということが数字で確認できる。
ケース2: 空調ダクト400×300・空気20℃・風速8m/s
入力 → 形状=矩形ダクト / 流体=空気 20℃ / 400×300mm / 風速 8m/s
結果 → 水力直径 Dh=342.857mm で Re = 181,480(乱流)。円形相当直径 De=377.709mm(Dh より +10.2%)、風量 3,456 m³/h、ν=15.1138 mm²/s
解釈 → この記事の山場。同じダクトを De=377.709mm で計算すると Re=199,928 になり、Dh 基準より10.2%大きい。Re に入れるのは Dh のほうで、De は圧損線図から等摩擦の丸ダクト径を読むための長さだ。どちらも「相当直径」と訳されるため、画面に両方を出して取り違えを防いでいる。
ケース3: 油圧配管25A・ISO VG46 40℃・流速4m/s
入力 → 形状=円管 / 流体=ISO VG46 40℃ / 内径 27.6mm / 流速 4m/s
結果 → Re = 2,400(遷移域)。ν=46.0 mm²/s、流量 8.62 m³/h。Re=2,300 になる流速は 3.83 m/s、Re=4,000 になる流速は 6.67 m/s
解釈 → 層流の上限 2,300 をわずかに超えただけの遷移域。流速を 3.83 m/s まで落とせば層流、6.67 m/s まで上げれば乱流で、設計流速 4m/s はいちばん決まらない帯のど真ん中にいる。圧損を見るなら乱流側、熱を見るなら層流側で評価する場面だ。
ケース4: 同じ配管の冷間始動(ISO VG46 20℃)
入力 → ケース3と同じ配管・同じ流速で、油温だけ 20℃
結果 → Re = 821(層流)。ν=134.49 mm²/s(40℃の2.92倍)。Re=2,300 になる流速は 11.21 m/s、Re=4,000 になる流速は 19.49 m/s
解釈 → 流速も配管も変えていないのに、油温が 20℃ というだけで遷移域から層流へ落ちる。この状態で遷移域を抜けるには 11.21 m/s が要る。動粘度が2.92倍になった分だけ、同じ Re を保つのに流速を上げなければならないからだ。冬の朝だけ圧損が跳ね上がる理由がここにある。
ケース5: 平板・空気20℃・流速10m/s・先端から500mm
入力 → 形状=平板 / 流体=空気 20℃ / 流速 10m/s / 先端からの距離 500mm
結果 → Re_x = 330,823(層流境界層)。遷移位置 x_cr = 0.756 m。ν=15.1138 mm²/s
解釈 → 先端から500mm の点は、まだ層流境界層の中。ただし遷移位置は先端から 0.756m なので、板が1m あるなら後半は乱流境界層になる。同じ板の上で、前半と後半で使う熱伝達の相関式が変わるということで、平板だけが持つ性質がそのまま出ている。
ケース6: 硬式野球ボールφ73mm・空気20℃・40m/s
入力 → 形状=球 / 流体=空気 20℃ / 直径 73mm / 流速 40m/s
結果 → Re = 193,201(ニュートン域)。抗力危機(Re=3×10⁵)に達する流速は 62.11 m/s
解釈 → 速球の域でも Re はニュートン域にとどまり、抗力係数が急落する抗力危機には届かない。届くのは 62.11 m/s、時速224km 相当の速度だ。プロの投球でも抗力危機は起きていない、というのが数字の答えになる。
ケース7: 撹拌槽・ISO VG68 20℃・翼径500mm・60rpm
入力 → 形状=撹拌槽 / 流体=ISO VG68 20℃ / 翼径 500mm / 回転数 60rpm
結果 → Re = 1,168(遷移域)。ν=214.12 mm²/s。乱流(Re=10⁴)になる回転数は 513.88 rpm
解釈 → 撹拌槽の乱流下限は 10⁴ なので、60rpm では遷移域に留まる。乱流にするには 513.88rpm 必要だが、翼径500mm の槽では現実的な回転数ではない。高粘度の槽は乱流にできないという前提で、動力数を遷移域の値で取ることになる。
7ケースのうち、流速も寸法も変えずに判定が動いたのがケース3と4。温度だけで遷移域と層流を行き来している。動粘度を温度から引けるようにした理由がここにある。
仕組み・実装の詳細
4つの設計判断
Re = vL/ν という定義式そのものに選択の余地は無い。実装で決めることは4つあった。
- 代表長さを形状ごとに切り替えるか、円管専用にするか — 切り替える。既存の日本語ツールが円管専用で、実務でいちばん間違うのが円管以外の L だからだ。ここを削るなら、このページを作る意味が無い。
- 矩形ダクトの「相当直径」をどちらの定義にするか — 両方出す。Re に使うのは Dh、De は圧損線図用だと結果に明記する。
- 動粘度をどう用意するか — 温度から引けるようにする。手入力だけでは cSt と cP の取り違えを誰も止められない。
- 臨界Reをどこに持つか — 形状ごとのバンドを1つのテーブルに集約し、判定を1関数に閉じ込める。各所に
Re < 2300を書き散らすと形状を足すたびに漏れる。警告の出し分け・数直線のバンド塗り・臨界流速の逆算は、すべてこの1テーブルを読んでいる。
水力直径 Dh と円形相当直径 De は別の量
水力直径(レイノルズ数に入れる長さ)
Dh = 4A / P = 4ab / 2(a+b) = 2ab / (a+b)
円形相当直径(圧損線図で等摩擦の丸ダクト径を読む長さ・Huebscher)
De = 1.30 (ab)^0.625 / (a+b)^0.25
Dh は「断面積を濡れ縁長さで割って4倍したもの」で、非円形断面の流れを円管に読み替えるための幾何学的な長さ。Hydraulic diameter の定義そのもので、この項目は円形相当直径に一切触れていない。一方 De は ASHRAE の等摩擦法で使われる Huebscher の実験式で、「同じ風量で単位長さあたりの圧力損失が等しくなる丸ダクトの径」を与える。目的の違う2つの長さに「相当直径」という同じ訳語が当たっているのが混同の原因で、英語なら hydraulic diameter と circular equivalent diameter なので取り違えにくい。
自社ツールではダクトサイズ計算ツールが De を使う側、このツールは Re に Dh を使う側。両方が同じ画面に出るようにしたのは、対になる関係を一目で見せるためだ。
計算例: 400×300mm・空気20℃・8m/s を1行ずつ
Dh = 2 × 400 × 300 / (400 + 300) = 240000 / 700 = 342.857 mm
De = 1.30 × (400 × 300)^0.625 / (400 + 300)^0.25 = 377.709 mm
ν(空気 20℃) = 15.1138 mm²/s = 15.1138×10⁻⁶ m²/s
Re(Dh) = 8 × 0.342857 / 15.1138×10⁻⁶ = 181,480 ← Re に使うのはこちら
Re(De) = 8 × 0.377709 / 15.1138×10⁻⁶ = 199,928 ← +10.2%
正方形ダクトだと Dh は辺長そのものになるので直感と合うが、扁平になるほど Dh は辺長から離れていく。ここを De で代用すると差は開く一方だ。
動粘度はどこから来ているか
水・海水(塩分3.5%)・空気・エチレングリコール水溶液30%/50% は、流体物性ライブラリ CoolProp 8.0.0 から1気圧の密度と粘度を取り、ν=μ/ρ を計算して表にした。孫引きの物性表は使っていない。
表は温度の刻みで持つが、動粘度は温度に対して指数的に変わる(水は0→100℃ で約6倍、ISO VG46 の油は0→100℃ で86倍)。線形に内挿すると刻みの中央で誤差が跳ねるので、ν は exp(lerp(ln ν)) の対数線形、ρ は線形で内挿している。この方式での再現誤差は全流体・全温度で1%以内。線形内挿だと1.4%まで悪化する。
ISO VG 油は表を持たない。ASTM D341(Walther)の式を2点から解いて、温度から直接計算している。
log10(log10(ν + 0.7)) = A − B × log10(T[K])
(40℃, ν40) と (100℃, ν100) の2点から A, B を解く
ν40 は ISO 3448 の中心値(VG32=32.0 / VG46=46.0 / VG68=68.0)
ν100 は粘度指数 VI=100 を満たす値を ASTM D2270 の表から逆算
40℃の値は ISO 3448(JIS K 2001)の中心値そのままなので、潤滑油の粘度グレード換算ブリッジと完全に同じ数字が出る。同じ量を2本のツールが別々の値で持つと、利用者から見れば単にサイトが矛盾しているだけなので、共有できる定数は共有した。
定数の根拠と、文献の幅
円管の 2,300 / 4,000 は日本語版 Wikipedia の記述と、自社の配管圧力損失計算・強制対流熱伝達係数の実装値を揃えた値だ。ただし英語版 Wikipedia は乱流側を 2,900 としていて、どこから乱流と呼ぶかの合意は文献で割れている。平板の 5×10⁵ も、乱れ強さや表面粗さによって 3×10⁵〜3×10⁶ まで動く。球の抗力危機は 3×10⁵ を採ったが、これは平滑球の値で、表面が粗ければもっと低い Re で起きる。ツール側では、この幅を注記として出している。
ISO VG 油の100℃側に VI=100 を仮定したのは、ISO 3448 が規定しているのが40℃の値だけで、100℃側はグレードだけからは決まらないからだ。VI=100 は鉱油系のパラフィン基油の基準値で、市販の VG46 作動油のカタログ値(100℃で 6.8 前後)とも整合する。VI が大きく異なる油(高VIの合成油など)を扱うなら、40℃以外の温度では実測値を直接入力に入れたほうが確実だ。
境界ちょうどの入力をどう扱うか
臨界流速を表示して、利用者がその値を流速欄に戻す —— という往復は、このツールの設計上ふつうに起きる。戻ってきた値が境界のわずかに下へ落ちれば、さっき「Re=2,300 になる流速」として出したはずの入力が「層流」と表示される。素直に Re < 2300 と書いていると、この食い違いを止める手段が無い。
実際に踏むのは、むしろ丸い入力のほうだ。内径52.9mm の管に ISO VG46 の40℃油を 2m/s で流すと、実数なら 2 × 0.0529 ÷ 46×10⁻⁶ は 2,300 ちょうどになる。ところが動粘度は ASTM D341 の式から復元した値なので 46.00000000000006 mm²/s であり、Re は 2,299.9999999999973 に着地する。素の比較では「層流」、正しくは遷移域だ。そこで判定は Re < 上限 × (1 − 1e-9) という相対許容で行っている。Re は 1 から 10⁷ まで7桁の幅があるので、絶対値の許容誤差では効かない。
既存のレイノルズ数計算ページと何が違うのか
既存の無料計算ページが円管専用で動粘度は手入力、という前提は前半で触れたとおり。足しているのは次の5点だ。
1. 形状5択で代表長さと臨界Reが同時に切り替わる。 円管・矩形ダクト・平板・球・撹拌槽。選ぶと、長さに何を入れるか(内径/水力直径/先端からの距離/球の直径/翼径)と、どの境界で判定するか(2,300と4,000/5×10⁵/1・1,000・3×10⁵/10と10⁴)が一緒に決まる。
2. 水力直径 Dh と円形相当直径 De を両方出す。 そのうえで Re に入れたのがどちらかを明示する。
3. 流体8種の動粘度と密度を温度から引く。 水・海水・空気・エチレングリコール30%/50%・ISO VG32/46/68。実測値があるなら直接入力もできる。
4. 臨界Reになる流速(撹拌槽は回転数)を逆算する。 「いまは乱流」だけでなく「どこまで絞れば層流に落ちるか」が同時に出る。
5. 平板は遷移位置を長さで出す。 板のどこから乱流境界層に変わるかが先端からの距離で返る。
関連ツールとの役割分担も分けてある。このツールは Re と流動状態を出すところで止まり、その先は目的で分かれる——圧力損失なら配管圧力損失計算、熱伝達係数なら強制対流熱伝達係数 計算ツール、所要動力なら撹拌機 動力計算シミュレーター、絞り流量計の流出係数ならオリフィス板流量計設計。臨界Reはこれらの実装と1つずつ突き合わせて揃えてあるので、「同じ配管なのにページによって判定が違う」ということは当サイト内では起きない。ダクトだけは対の関係で、ダクトサイズ計算ツールが De(等摩擦の丸ダクト径)を使う側、このツールが Dh(Re に入れる側)を使う側になる。
染料の一筋がほどける瞬間——レイノルズ数にまつわる話
ガラス管と染料。 この数が生まれたのは1883年、オズボーン・レイノルズの実験だ。水で満たしたガラス管の中心に細い染料の筋を流し込み、弁を少しずつ開いて流速を上げていく。ある速さを越えると、それまでまっすぐ伸びていた筋が突然ほどけて管いっぱいに散る。彼が論文で使った語は「direct(まっすぐな)」と「sinuous(くねった)」で、いまの層流・乱流という呼び名はここから来ている。染料を一筋流すという素朴きわまりない可視化から、無次元数という抽象が出てきたのが面白い(Osborne Reynolds)。
2,300 という半端な数字の正体。 なぜ2,000でも2,500でもないのか。臨界値が実験の丁寧さに強く依存するからだ。入口の擾乱を徹底的に抑えた装置では、はるかに高い Re でも層流が保たれることが後の実験で分かっている。つまり 2,300 は「これを超えたら必ず乱流」ではなく、「これ以下なら擾乱を与えても層流に戻ってくる」という下限側の値なのだ。上側は「どこから乱流と呼ぶか」の合意そのものが無いので、文献ごとに数字が割れる。英語版 Wikipedia の Reynolds number は乱流側を 2,900 とし、日本の実務書は 4,000 と書く。どちらかが誤りなのではなく、遷移という連続的な現象のどこに線を引くかの違いにすぎない。
表面を汚くしたほうが速い。 球まわりの流れには、Re が 3×10⁵ あたりで抗力係数が 0.5 から 0.2 付近へ急落するという有名な逆転がある(抗力危機)。境界層が層流のまま剥がれていたのが、乱流化して剥離点が後ろへずれ、後流が細くなるためだ。表面を粗くすると境界層が早く乱れるので、この急落がもっと低い Re で起こる——ゴルフボールのディンプルは、それを意図的に前倒しする仕掛けだ。滑らかに磨き上げるより凹ませたほうが飛ぶという、直感に逆らう工学の代表例になっている(Drag crisis)。
使いこなしの小技
- 水20℃なら「Re ≈ 1000 × 流速[m/s] × 内径[mm]」で暗算できる。 水20℃の動粘度が 1.0034 mm²/s とほぼ 1 だからで、この近似と実際の値は1%も違わない。打ち合わせ中に桁だけ合わせたいときに効く。
- 空気の動粘度は水の約15倍ある。 20℃で水が 1.0034、空気が 15.1138 mm²/s。同じ寸法・同じ速度なら空気のほうが Re は15分の1になる。「空気はサラサラだから乱れやすい」という感覚は、動粘度で見ると逆向きだ。
- 数直線は対数目盛だと意識して読む。 流速を10倍してもマーカーは1目盛ぶんしか動かない。Re は桁で捉える量であって、有効数字3桁目を気にする量ではない。
- 判定が境界の前後1割に入ったら断定しない。 その帯では入口形状や振動といった、式に入っていない条件で実物の挙動が動く。両側の式で計算して差の大きさを見ておくほうが、片方を信じるより安全だ。
- 矩形ダクトの幅と高さは内寸で測る。 保温材込みの外寸や板厚を含んだ寸法を入れると Dh がそのまま外れ、Re も同じ割合でずれる。図面の呼称寸法が内寸か外寸かは会社ごとに違うので、一度確認しておく価値がある。
よくある質問
矩形ダクトの「相当直径」は、Dh と De のどちらを入れればいい?
レイノルズ数に入れるのは水力直径 Dh のほうだ。De(円形相当直径)は摩擦損失線図やダクトサイズ表で「同じ圧損・同じ風量になる丸ダクト径」を引くための量で、Re の定義に代入するための量ではない。同じ訳語で呼ばれる別物が2つある事情と差の出方は、前半の仕組みの節で計算例つきに説明している。使う側としては、矩形ダクトを選ぶと Dh と De の両方が結果に並び、Re に使ったのがどちらかも書かれるので、手元の資料の「相当直径」がどちらの値に近いかを突き合わせれば取り違えにその場で気づける。
動粘度の単位は cSt と mm²/s と m²/s のどれで入れるのか
入力欄の単位は mm²/s で、cSt(センチストークス)と mm²/s はまったく同じ量だ。1 cSt = 1 mm²/s = 1×10⁻⁶ m²/s なので、cSt 表記のカタログ値はそのまま打ち込める。注意が要るのは cP(センチポアズ=mPa·s)で、こちらは動粘度ではなく粘度。密度で割らないと使えない(ν[mm²/s] = μ[mPa·s] ÷ ρ[g/cm³])。換算を忘れて cP の数字を直接入れると Re が桁違いにずれるので、直接入力に切り替えたときは欄の下に同じ注意が出る。グレードを跨いだ粘度換算そのものは潤滑油の粘度グレード換算ブリッジ側の担当だ。
撹拌槽だけ計算式が違うのはなぜか
撹拌槽には「管の内径」にあたる寸法も「平均流速」にあたる速度も無いので、代表長さに翼径 d、代表速度に n·d(n は毎秒回転数)を使うという別の約束が置かれている。結果として Re = n·d²/ν という、他の形状とは形の違う式になる。実務で多い事故は n を毎分回転数のまま入れることで、そのままだと Re が60倍になる。このツールは回転数を rpm で受けて内部で毎秒に直すので、カタログの rpm をそのまま入れてよい。もう一つ、長さに入れるのは翼径であって槽径ではない。槽径を入れると Re が数倍大きく出る。
臨界レイノルズ数が資料によって違うのはなぜか。このツールはどの値を使っている?
どの値を採ったか、なぜその値にしたかは前半の仕組みの節で出典つきに説明している。使う側にとって大事なのは、臨界Reが当サイトの配管・熱・撹拌のツールと同じ値に揃えてあることだ。円管と矩形ダクトの 2,300 / 4,000 は配管圧力損失計算・強制対流熱伝達係数 計算ツールと、平板の 5×10⁵ は強制対流熱伝達係数 計算ツールと、撹拌槽の 10 / 10⁴ は撹拌機 動力計算シミュレーターと一致する。Re をここで出して各ツールへ渡しても、途中で判定が入れ替わることはない。適用した臨界Reは結果に書き出しているので、社内基準が別の値なら差を見たうえで読み替えてほしい。
流体の有効温度範囲を外れた温度を入れるとどうなるのか
結果を消さずに、範囲の端の物性で計算を続けて「範囲外なので何℃の値で計算した」という注記を出す。ここで計算を止める作りにすると、空気で150℃を入れたあと水に切り替えた瞬間に結果が真っ白になり、流体どうしの比較ができなくなる。有効範囲は流体ごとに違い(水と油は0〜100℃、海水は0〜40℃、空気がいちばん広い)、選択中の流体の範囲は温度欄の下に出る。範囲外の値は外挿ではなく端で頭打ちになっているだけなので、設計の数字にするなら実測値を直接入力に入れたほうがいい。
入力した寸法や温度はどこかに送信されるのか
されない。計算はすべてブラウザの中で完結していて、入力値がサーバーへ送られることも保存されることもない。動粘度の表も換算式もページに同梱してあり、計算のたびに外部へ問い合わせる作りにはなっていない。社外に出せない配管径や運転条件を入れても手元から出ないので、その点は気にせず使ってほしい。結果はコピーボタンでテキストとして取り出せるので、必要なら自分で計算書や議事録に貼り付ければいい。
まとめ
レイノルズ数の式そのものは vL/ν の一行で、難しいのは v でも ν でもなく L の決め方のほうだ。形状ごとに長さの取り方も判定の境界も変わり、矩形ダクトに至っては「相当直径」という同じ訳語が2つの別の量を指している。形状を選べば両方が同時に切り替わるようにしたのは、そこを間違えたまま先へ進むのがいちばん高くつくからだ。Re が出たら、圧力損失は配管圧力損失計算、熱伝達係数は強制対流熱伝達係数 計算ツール、撹拌の所要動力は撹拌機 動力計算シミュレーターへ渡してほしい。この形状も入れてほしい、判定が手元の資料と合わない、といった指摘はお問い合わせページから歓迎する。
Mahiro
Mahiro Appの開発者。空調ダクトの相当直径を圧損線図の De のままレイノルズ数に入れて、10%ずれた値で熱交換器を選定しかけた設備設計者。以来、Re を出すときは「長さに何を入れたか」を必ず書き残すようにしている。
運営者情報を見る