フーリエ数（Fo）— 非定常熱伝導の無次元時間

フーリエ数（Fo）の理論基礎

フーリエ数の定義と物理的意味

先生、フーリエ数って何を意味するんですか？ 教科書で Fo って記号が突然出てきて、よく分からなくて…。

ざっくり言うと、フーリエ数は「非定常熱伝導がどれだけ進んだか」を表す無次元の時間だ。定義式はこうなる：

ここで $\alpha$ は熱拡散率 [m²/s]、$t$ は経過時間 [s]、$L$ は代表長さ（板厚の半分や円柱の半径など）[m] だ。

熱拡散率 $\alpha$ って、熱伝導率 $k$ とは違うんですか？

いい質問だ。熱拡散率は熱伝導率を密度と比熱で割ったものだ：

熱伝導率 $k$ が「熱をどれだけ伝えやすいか」なのに対して、熱拡散率 $\alpha$ は「温度変化がどれだけ速く伝わるか」を表す。銅は $k$ が大きいから熱を伝えやすいけど、$\rho c_p$ も結構大きい。一方アルミは $\alpha$ が銅より大きいので、温度変化の伝搬は実はアルミのほうが速いんだ。

なるほど！ じゃあフーリエ数が大きいほど「熱が十分行き渡った」ってことですか？

その通り。例えるなら、お風呂に入った直後は表面だけ温まってるけど（Fo が小さい）、しばらく浸かっていると芯まで温まる（Fo が大きい）。フーリエ数は「物体全体の温度がどれだけ均一化したか」の目安なんだ。

非定常熱伝導の支配方程式と無次元化

フーリエ数は支配方程式のどこから出てくるんですか？

1次元の非定常熱伝導方程式（内部発熱なし、等方性・一定物性）を見てみよう：

ここで無次元変数を導入する。温度を $\theta = (T - T_\infty)/(T_i - T_\infty)$、位置を $\xi = x/L$、時間を $\tau = \alpha t / L^2$ と置く。すると支配方程式はこうなる：

見てのとおり、無次元時間 $\tau$ がまさにフーリエ数 $\mathrm{Fo}$ だ。無次元化すると物性値や寸法のパラメータが全部消えて、同じ Fo と同じ境界条件なら、材料やサイズが違っても同じ無次元温度分布 $\theta(\xi, \mathrm{Fo})$ になる。これが相似則の本質だよ。

え、つまり 10cm のアルミ板と 1m の鉄板でも、Fo が同じなら温度分布のパターンは同じってことですか？

境界条件を同じ Bi（ビオ数）で揃えればね。まさにその考え方で実験の相似模型を作ったり、ハイスラーチャートのような汎用グラフが成立するんだ。

一項近似と Fo > 0.2 の意味

教科書に「Fo > 0.2 なら一項近似が使える」って書いてあったんですけど、これはどういう意味ですか？

非定常熱伝導の解析解は一般に無限級数で表される。例えば、初期温度 $T_i$ の無限平板が急に温度 $T_\infty$ の環境に曝された場合、中心の無次元温度は：

ここで $\zeta_n$ はビオ数から決まる固有値だ。指数関数 $\exp(-\zeta_n^2 \, \mathrm{Fo})$ は $n$ が大きいほど急速に減衰する。Fo が 0.2 を超えると、$n \geq 2$ の項は全て 2% 以下になって無視できる。だから第1項だけ残せば十分な精度が得られる：

Heislerチャートを使うときも Fo > 0.2 が条件って書いてありました。これもそういう理由なんですね！

そのとおり。Heislerチャートは一項近似の結果をグラフ化したものだから、Fo < 0.2 の領域では使えない。初期の短時間領域が必要な場合は、級数を多くの項まで足すか、半無限体の解を使う必要がある。

具体的に、Fo > 0.2 になるのってどれくらいの時間なんですか？

例えば鋼（$\alpha \approx 1.2 \times 10^{-5}$ m²/s）で板厚 20mm（$L = 10$ mm = 0.01 m）の場合：

鋼板 20mm なら約 1.7 秒後には一項近似が使える。一方、コンクリート壁 200mm（$\alpha \approx 5 \times 10^{-7}$ m²/s、$L = 0.1$ m）だと $t \approx 4000$ 秒 = 1 時間以上かかる。材料と寸法で桁違いに変わるんだ。

ビオ数との関係と解析手法の選定

フーリエ数とビオ数って、セットで使うものなんですか？

まさにそう。非定常熱伝導問題の解析手法を決めるとき、ビオ数 Bi とフーリエ数 Fo はペアで考える。ビオ数は：

$h$ は表面の熱伝達係数、$L$ は代表長さ、$k$ は固体の熱伝導率。Bi は「固体内部の熱抵抗」と「表面の対流熱抵抗」の比を表す。

Bi が小さいと集中熱容量法が使えるんですよね？

そう。Bi < 0.1 なら物体内部の温度勾配が無視できるほど小さいので、物体全体を一様温度として扱える——これが集中熱容量法（Lumped Capacitance Method）だ。このとき温度変化は単純な指数減衰になる：

実務での判断フローをまとめるとこうなる：

Bi と Fo の2つをチェックするだけで、どの解法を使うか決められるんですね。すごくスッキリしました！

フーリエ数（Fo）の数値計算手法

要素フーリエ数と陽解法の安定条件

先生、FEM の非定常熱解析でもフーリエ数って出てきますよね？

出てくる。数値解析では要素フーリエ数（メッシュフーリエ数）が重要になる。代表長さ $L$ の代わりに要素サイズ $\Delta x$ を、時間 $t$ の代わりに時間刻み $\Delta t$ を使う：

陽解法（前進Euler法、前進差分法）を使う場合、$\mathrm{Fo}_e \leq 0.5$（1次元）でないと解が発散する。これが有名な CFL 条件の熱伝導版だ。2次元なら $\mathrm{Fo}_e \leq 0.25$、3次元なら $\mathrm{Fo}_e \leq 1/6$ になる。

メッシュを細かくすると時間刻みもすごく小さくしないといけないんですか？

そう、陽解法では $\Delta t \leq \mathrm{Fo}_{e,\max} \cdot \Delta x^2 / \alpha$ だから、メッシュサイズ $\Delta x$ を半分にすると許容時間刻みは4分の1になる。計算ステップ数は2倍になるから、合計の計算量は8倍に膨らむ。これが陽解法の最大の弱点だ。

時間刻みの設定指針

じゃあ実務では時間刻みをどうやって決めるんですか？

フーリエ数ベースの考え方が便利だ。まず解析全体の時間スケールを見積もる：

• 特性時間 $t_c = L^2 / \alpha$ ——温度変化が物体全体に伝わるまでの目安時間
• 解析終了時間——定常に近づく（Fo $\approx$ 1〜2）か、関心のある時間まで

次に時間刻みの目安：

陰解法なら Fo_e をいくら大きくしても大丈夫なんですか？

安定性の面では大丈夫だけど、精度は落ちる。時間刻みが大きすぎると温度の過渡的な変化を捉え損なう。特に熱衝撃の初期応答を見たい場合は、最初の数ステップは小さい $\Delta t$ を使って、温度変化が緩やかになったら大きくする「自動時間刻み制御」が実務の定番だ。

陰解法と陽解法の選択

結局、どっちを使えばいいんですか？

一般的な非定常熱伝導では陰解法が圧倒的に有利だ。理由は単純で、熱伝導は拡散型の方程式だから時間スケールが長い。陽解法だと安定性のために非常に小さい $\Delta t$ を強いられるけど、物理的には温度変化はゆっくりなので、陰解法で大きな $\Delta t$ を使ったほうが効率的だ。

ただし例外もある：

• レーザー加熱やプラズマ照射のような超短時間の局所加熱 → 陽解法のほうが適切な場合がある
• 構造-熱の陽的連成（LS-DYNA等）→ 構造側が陽解法なので熱側も陽解法に合わせることがある

フーリエ数（Fo）の実務適用

フーリエ数を使った解析設計フロー

先生、実務でフーリエ数を使う場面を具体的に教えてもらえますか？

CAE で非定常熱解析をセットアップするとき、フーリエ数は以下の4つの場面で使うよ：

1. 解析終了時間の見積もり：物体が定常に近づく時間 $t_c = L^2/\alpha$ を計算して、解析時間の目安にする。Fo $\approx$ 1〜2 で概ね定常。
2. 時間刻みの妥当性チェック：要素フーリエ数 Fo_e を計算して、安定性条件や精度の確認。
3. メッシュと時間刻みのバランス確認：メッシュを細分化したら Fo_e が大きくなりすぎていないかチェック。
4. 結果の検証：中心温度の時刻歴が一項近似の理論値と合っているか、Fo > 0.2 の範囲で比較。

解析を始める前にフーリエ数を電卓で計算しておくだけで、かなり見通しが良くなりそうですね。

その通り。CAE はブラックボックスになりがちだけど、フーリエ数とビオ数の2つだけでも事前に計算しておくと、「この解析結果はおかしい」と気づける力がつく。結果を鵜呑みにしないためのリテラシーだよ。

代表的な材料のフーリエ数計算例

いくつか具体的な計算例を見せてもらえますか？

よくある設定で計算してみよう：

食品の殺菌って、15分も経たないと一項近似すら使えないんですね…。レトルトカレーの加工って大変なんだ。

食品メーカーでは中心温度が規定温度に達する時間を正確に制御するのが品質管理の要だからね。フーリエ数はそういう現場でも使われている概念なんだ。

よくある失敗と対策

フーリエ数がらみで初心者がやりがちな失敗ってありますか？

3つの典型的な失敗パターンを紹介しよう：

ソフトウェアでの扱い

各ソルバーでのFo関連設定

商用CAEソフトでフーリエ数って意識するところはありますか？

ソルバー自体が「Fo を入力してください」と聞いてくることはないけど、時間刻みの設定と安定性の判断の裏では常にフーリエ数が関わっている。具体的に見ていこう：

ソフトが自動で $\Delta t$ を調整してくれるなら、手動でフーリエ数を計算する意味ってあるんですか？

大あり。自動制御は万能じゃない。特に以下のケースでは手動チェックが重要だ：

• 初期時間刻みの設定ミス：自動制御の「初期値」が大きすぎると、最初のステップで非物理的な温度分布が出て、以降のステップも引きずる
• 局所的に細かいメッシュがある場合：全体の平均メッシュサイズでは問題なくても、最小要素で Fo_e が大きすぎることがある
• 結果の妥当性検証：「この時間で Fo はいくつか？」を計算して、理論的に温度がどの程度均一化しているはずかを事前に見積もれる

発展的トピック

多次元問題でのフーリエ数

ここまで1次元の話でしたけど、2次元・3次元ではどうなるんですか？

有限の長方形体（$2L_1 \times 2L_2 \times 2L_3$）の過渡温度分布は、各方向の1次元解の積で表せる（変数分離法による積の解）：

ここで $\mathrm{Fo}_i = \alpha t / L_i^2$ だ。つまり各方向でフーリエ数を独立に計算して、最も小さい方向（最も厚い方向）が温度均一化のボトルネックになる。例えば薄い板なら板厚方向が最初に均一化して、面内方向は後から追いつく。

数値解析の安定条件も3次元だと変わるんですよね？

そう。等方的なメッシュ（$\Delta x = \Delta y = \Delta z$）の場合、3次元の陽解法安定条件は：

1D なら $\leq 1/2$、2D なら $\leq 1/4$、3D なら $\leq 1/6$。次元が増えるほど安定条件が厳しくなるから、3次元の陽解法は本当に時間刻みが小さくなる。だからこそ実務では陰解法が主流なんだ。

非線形問題への拡張

温度依存の物性値（$k$ が温度で変わる場合）だとフーリエ数はどう扱うんですか？

厳密には、$\alpha$ が温度の関数 $\alpha(T)$ になるから、フーリエ数は一定値にならない。このときは：

• 参照温度での $\alpha$ を使って概算値を出す（初期温度や平均温度で評価）
• Kirchhoff変換：$U = \int_0^T k(T')\,dT'$ で変数変換すると、新しい変数について線形方程式に帰着できる場合がある
• 相変化（融解・凝固）がある場合は潜熱の影響で $c_p$ が見かけ上巨大になるため、Fo が極端に小さくなる時間帯がある。この領域ではフーリエ数による見積もりは使えないので、数値解析に頼るしかない

非線形でも、概算としてフーリエ数を使う価値はあるってことですね。

そう。完璧な値は出せなくても、オーダー推定には十分使える。「この問題の特性時間は数秒オーダーなのか数時間オーダーなのか」——それだけでも解析の設計が全然違ってくるからね。

フーリエ数（Fo）のトラブル対応

Foに関連する典型的なトラブル

先生、非定常熱解析でよくあるトラブルを教えてください。

フーリエ数に関連するトラブルと対策をまとめておこう：

1. 温度が振動する（非物理的な温度の上下動）

• 原因：Crank-Nicolson 法で Fo_e が大き過ぎる場合に起きやすい。完全陰解法（後退Euler）では振動しないが、Crank-Nicolson は2次精度の代償として短時間の振動モードを持つ
• 対策：$\Delta t$ を小さくするか、完全陰解法に切り替える。ANSYS の場合は「Theta = 1.0」に設定

2. 解が発散する（温度が無限大に飛ぶ）

• 原因：陽解法で安定条件 Fo_e ≤ 1/(2d) を超えている。特に局所的に細かいメッシュ部分が原因になりやすい
• 対策：最小メッシュサイズで Fo_e を確認。$\Delta t \leq \Delta x_{\min}^2 / (2d \cdot \alpha)$ を満たすようにする

3. 定常に達するまでの時間が予想と違う

• 原因：$\alpha$ の値が桁違い（物性データの単位ミス、温度依存性の未考慮）、または $L$ の定義が不適切
• 対策：$t_c = L^2/\alpha$ で特性時間を手計算して比較。2倍程度のずれなら形状効果、10倍以上なら物性値か $L$ の定義を再確認

4. Heislerチャートと数値解析の結果が合わない

• 原因：Heislerチャートは一項近似（Fo > 0.2）の無限平板/円柱/球が前提。有限体、複合材、非均一境界条件では適用外
• 対策：前提条件を確認。Fo > 0.2 かつ形状が単純（1D 近似可能）な場合のみ比較に使う

全部フーリエ数を事前にチェックしていれば防げるトラブルですね。手計算の大切さがよく分かりました！

「解析が合わない」と思ったら、まず Fo と Bi を計算するところから始めるのが鉄則だ。この2つの無次元数は、非定常熱伝導問題の羅針盤みたいなものだからね。