形状係数
形状係数の理論基礎
形状係数（Shape Factor）とは
先生、形状係数って何に使うんですか？
2次元や3次元の定常熱伝導問題を、1次元の熱抵抗に帰着させるための概念だ。形状係数 $S$ を使えば
と表現でき、$\Delta T$ から直接放熱量が分かる。熱抵抗は $R = 1/(kS)$ だ。
複雑な形状を1つの数値に集約できるんですね。
$S$ は形状と境界条件だけで決まる幾何学的な量で、単位は [m]（2D問題では [m/m] = 無次元/単位奥行）だ。
代表的な形状係数
| 形状 | $S$ | 適用条件 |
|---|---|---|
| 無限平板 | $A/L$ | 基本形 |
| 同心円筒 | $2\pi L / \ln(r_2/r_1)$ | 長さ $L$ |
| 同心球 | $4\pi r_1 r_2/(r_2 - r_1)$ | — |
| 埋設球（半無限体表面から深さ$z$） | $4\pi r / (1 - r/(2z))$ | $z > r$ |
| 埋設円筒（半無限体） | $2\pi L / \cosh^{-1}(z/r)$ | $z > r$, $L \gg r$ |
| 2つの平行円筒 | $2\pi L / \cosh^{-1}((d^2-r_1^2-r_2^2)/(2r_1 r_2))$ | 中心間距離 $d$ |
埋設球や埋設円筒の式が実用的そうですね。
地中埋設配管の放熱量計算、建物基礎からの地中伝熱、地熱ヒートポンプの設計などで使われる。手計算で概算が得られるのが形状係数の最大の利点だ。
形状係数の導出
ラプラス方程式 $\nabla^2 T = 0$ の解から導出する。同心円筒の場合
よって $S = 2\pi L / \ln(r_2/r_1)$ が得られる。
解析解から形状係数を逆算するんですね。
解析解が得られない形状ではFEMで数値的に $S = q/(k\Delta T)$ を算出する。
形状係数の定義と物理的意味
形状係数Sは複雑形状での熱流量をq=SkΔTで表す無次元係数だ。1950年代にCarslaw & Jaegerが「固体の熱伝導」で体系化し、埋設配管から球体・円柱まで60種以上の解析解が表として整備された。
形状係数の数値計算手法
数値的な形状係数の算出
複雑な形状の形状係数はどうやって求めるんですか？
FEMで定常熱伝導を解き、以下の手順で $S$ を算出する。
1. 高温面に $T_1$、低温面に $T_2$ のDirichlet条件を設定
2. $k = 1$ W/(m K) に設定（簡略化のため）
3. 解析を実行し、高温面（または低温面）の全熱流量 $q$ を取得
4. $S = q / (k \cdot \Delta T) = q / (T_1 - T_2)$ を算出
$k = 1$ にするのは計算を楽にするためですね。
そう。$S$ は $k$ に依存しない幾何量なので、$k$ の値は結果に影響しない。
メッシュ収束の確認
形状係数は積分量（全熱流量）なので局所的なメッシュ感度は小さいが、形状が複雑な場合は収束確認が必要だ。
| メッシュレベル | 要素数 | $S$ | 誤差 |
|---|---|---|---|
| 粗い | 1,000 | 15.2 m | — |
| 中間 | 10,000 | 15.8 m | 3.9% |
| 細かい | 100,000 | 15.9 m | 0.6% |
| 非常に細かい | 1,000,000 | 15.9 m | 0.0% |
積分量は比較的早く収束しますね。
温度の積分量は応力の局所値より収束が速い。1万要素でも実用精度が得られることが多い。
重ね合わせの原理
形状係数は線形問題なので重ね合わせが可能だ。複数の熱源がある場合、各熱源からの形状係数を足し合わせられる。また対称条件を使えば計算量を削減できる。
対称面で半分のモデルにして $S$ を2倍にすればいいんですね。
その通り。地中埋設円筒の場合、地表面を対称面（断熱面）として扱うと、半無限体から有限領域に問題を縮小できる。
形状係数の解析解一覧の使い方
地中埋設管（直径D、深さz、長さL）の形状係数はS=2πL/ln(4z/D)（z>>D時）。例えばD=100mm、z=1mの東京都水道管20mではS≈27m、k=1.5 W/m·Kの土壌でq≈243WのΔT毎の熱流量が得られる実用計算だ。
形状係数の実務適用
地中埋設配管の放熱計算
形状係数が最も活用される場面は何ですか？
地中埋設配管の放熱量計算だ。地中の温度場は半無限体として扱え、形状係数の公式がそのまま使える。
計算例
蒸気配管（外径114.3mm、断熱材外径214.3mm）を地表から1.5mの深さに埋設。地中温度15℃、断熱材外面80℃。地盤の $k_{\text{soil}} = 1.5$ W/(m K)。
1mあたり184Wの放熱ですか。
100mの配管なら18.4kWの熱損失だ。年間のエネルギーコストに換算して断熱厚の経済最適化を行う。
建物基礎の地中伝熱
スラブオングレード（土間コンクリート床）からの地中伝熱もISO 13370で形状係数ベースの計算が規定されている。
| パラメータ | 影響 |
|---|---|
| 基礎面積/周長比 | 大きいほど地中への放熱が少ない |
| 断熱材配置 | 基礎外周部の断熱が最も効果的 |
| 地盤の $k$ | 砂 1.5、粘土 1.0、泥炭 0.5 W/(m K) |
大きな建物ほど有利なんですね。
面積が $L^2$ に比例し周長が $L$ に比例するため、大きな建物は周長あたりの面積が大きく、地中伝熱の影響が相対的に小さくなる。
形状係数の参考文献
形状係数の網羅的な一覧はIncropera「Fundamentals of Heat and Mass Transfer」のTable 4.1、またはBejan「Heat Transfer」のAppendixに収録されている。特殊形状はHahne & Grigullの文献が充実している。
教科書に載っていない形状はFEMで求めるしかないですね。
その通り。一度FEMで求めた $S$ は設計式としてデータベース化しておけば、以降は手計算で使える。
地中埋設配管の熱損失評価
北海道の地域熱供給では地中80cmに埋設する100A鋼管から年間熱損失が約15 W/mに達する。形状係数を使った解析で断熱被覆の最適厚さを30mmと算出し、ポリウレタンフォーム（k=0.027 W/m·K）の採用で損失を5 W/m以下に低減した事例がある。
形状係数のソフトウェア比較
商用ツールでの形状係数算出
形状係数を商用ツールで求めるにはどうしますか？
どのFEMツールでも基本手順は同じだ。
| ツール | 手順 | 熱流量の取得方法 |
|---|---|---|
| Ansys Mechanical | T固定面のReaction Heat Flowを取得 | FSUM, PRRSOL |
| Abaqus | Output > Node > RFL11（反力熱流量） | History Output |
| COMSOL | Surface Integration > Total Heat Flux | Derived Values |
APDL実装例
埋設円筒の形状係数を数値的に求めるAPDLコード。
```
/PREP7
ET,1,PLANE55,,,1 ! 軸対称ではなく平面
MP,KXX,1,1.0 ! k=1（形状係数算出用）
! 半無限体を有限領域で近似
RECTNG,0,5,0,-3 ! 5m幅 × 3m深さ
CYL4,0,-1.5,0.107,,,, ! 埋設円筒（深さ1.5m, r=107mm）
ASBA,1,2 ! 差分
ESIZE,0.02
AMESH,ALL
/SOL
D,円筒面NODE,,80 ! 管外面80℃
D,地表面NODE,,15 ! 地表15℃
! 側面・底面は断熱（デフォルト）
SOLVE
*GET,Q,FSUM,,HEAT ! 全熱流量
S_factor = Q / (1(80-15)) ! S = q/(kDT)
```
$k = 1$ にしてるから $S = Q/\Delta T$ ですね。
理論値 $S = 1.89$ m/m と比較して検証する。有限領域の大きさが十分か（側面の影響がないか）も確認が必要だ。
地熱設計専用ツール
地熱ヒートポンプの設計では専用ツールが使われる。
| ツール | 用途 |
|---|---|
| GLD (Ground Loop Design) | ボアホール熱交換器の設計 |
| EED (Earth Energy Designer) | 地中熱交換器の長期性能予測 |
| COMSOL Geothermal | 3D地中温度場のFEM解析 |
| FEFLOW | 地下水流を考慮した地中伝熱 |
地下水の流れも影響するんですね。
地下水流がある場合、対流による熱輸送で形状係数が大幅に変わる。純粋な伝導問題ではなくなるため、FEFLOWやCOMSOLでの数値解析が必要になる。
形状係数計算のCAEツール活用
Autodesk CFD 2025は地中埋設構造の形状係数を自動的にポスト処理で出力する機能を持つ。OpenFOAMのlaplacianFoamソルバーで定常熱伝導を解いた後、surfaceIntegrateコマンドで熱流束を積算してS値を算出するワークフローも確立されている。
形状係数の先端研究
複合形状の形状係数
教科書に載っていない複雑な形状はどう扱うんですか？
複合形状を簡単な形状に分解して重ね合わせる方法がある。L字型の壁はストレート部分とコーナー部分に分解し、各部の形状係数を足し合わせる。
コーナーの形状係数はコーナー角度と壁厚で決まる補正値で、教科書に表がある。
パズルみたいに分解するんですね。
分解が難しい場合はFEMで直接算出する。一度求めた $S$ を設計式にフィッティングすれば、パラメトリックな設計計算に使える。
非定常問題への拡張
非定常問題では形状係数が時間の関数 $S(t)$ になる。地中埋設配管の起動時の過渡応答は
$\alpha$ は地盤の温度拡散率、$\gamma = 0.5772$ はオイラー定数。長時間後に定常の $S$ に漸近する。
地中の温度場が定常に達するまでどのくらいかかりますか？
配管サイズと地盤の熱拡散率による。典型的には数日〜数週間。地熱ヒートポンプでは数年かかることもある。
異方性媒体の形状係数
地盤が層状構造（水平方向と鉛直方向で $k$ が異なる）の場合、座標変換で等方性問題に帰着させる。
変換後の等方性媒体（$k_{\text{eff}} = \sqrt{k_x k_y}$）で形状係数を求め、元の座標系に戻す。
座標をスケーリングして等方性に変換するのは賢いですね。
地質の層状構造はまさにこのケースだ。水平方向の透水層は鉛直方向より $k$ が大きいことが多い。
数値形状係数の抽出手法
解析解のない複雑形状の形状係数はFEMで数値的に抽出できる。熱流量qをFEMで計算後、S=q/(kΔT)で逆算する手法をFENICS Projectで実装すると、任意の3D形状に対して1時間以内にSを取得できる。
形状係数のトラブル対応
よくあるトラブルと対策
形状係数の計算で注意すべきことは何ですか？
典型的な問題を整理する。
1. 有限領域サイズの影響
問題: 半無限体をFEMで解くとき、有限領域で打ち切る。領域が小さすぎると $S$ が過小評価される。
対策: 領域の境界を関心領域の10倍以上離す。境界サイズを2倍にして $S$ が1%以内で変わらなければ十分だ。
2. 形状係数の適用条件の逸脱
教科書の公式をそのまま使って間違えるケースですね。
よくある誤用を挙げる。
| 公式 | 適用条件 | よくある誤用 |
|---|---|---|
| 埋設円筒 | $L \gg r$, $z > r$ | 短い管（$L/r < 10$）に適用 |
| 埋設球 | $z > r$ | 地表面に近い球（$z \approx r$）に適用 |
| 2平行円筒 | $d > r_1 + r_2$ | 円筒が接触する場合に適用 |
3. 地盤物性の不確実性
地中伝熱の最大の不確実性は地盤の熱物性だ。
| 地盤 | $k$ [W/(m K)] | 変動幅 |
|---|---|---|
| 乾燥砂 | 0.3 | $\pm$50% |
| 湿潤砂 | 2.0 | $\pm$30% |
| 粘土 | 1.0 | $\pm$40% |
| 岩盤 | 2.5 | $\pm$20% |
地盤の $k$ は含水率で数倍変わるんですね。
TRT（Thermal Response Test）で現地計測するのが最も確実だ。ボアホールに熱負荷を与えて温度応答から $k$ を逆算する。地熱ヒートポンプ設計では必須の工程だ。
4. 重ね合わせの限界
形状係数の重ね合わせは線形問題でのみ有効。放射を含む場合や温度依存 $k$ の場合は重ね合わせが使えない。
高温問題では要注意ですね。
その場合はFEMで直接解くしかない。形状係数は概算ツールであり、限界を理解して使うことが重要だ。
形状係数表の適用条件ミス
Carslaw & Jaegerの形状係数表には各解析解の有効条件（z/D比など）が明記されているが、範囲外適用によって誤差50%超になることがある。地中埋設管でz/D<2の浅い場合は表の式が使えず、Moody補正式か直接FEM解析が必要だ。
なった
詳しく
報告