1次元非定常熱伝導（半無限体） — ソルバー別比較
クロスバリデーション結果
NAFEMS T1の条件で各ソルバーの結果を比較してもらえますか？
20要素（バイアスメッシュ）、$\Delta t = 0.32$ s（100ステップ）での結果だ。
| ソルバー | 時間スキーム | T(x=0.08, t=32) [°C] | 参照値との差 [°C] |
|---|---|---|---|
| NAFEMS参照値 | — | 36.60 | — |
| Abaqus (DC1D2) | 後退Euler | 36.55 | -0.05 |
| Nastran (SOL 159) | θ=0.5 | 36.62 | +0.02 |
| Ansys (PLANE55) | 後退Euler | 36.54 | -0.06 |
| CalculiX (DC2D8) | 後退Euler | 36.56 | -0.04 |
| COMSOL (BDF2) | BDF2 | 36.61 | +0.01 |
| OpenFOAM (laplacianFoam) | Euler | 36.48 | -0.12 |
OpenFOAMだけ少し誤差が大きいですね。
laplacianFoamのデフォルト時間スキームはEuler（1次精度前進）だから、同じ$\Delta t$では2次精度の他ソルバーより誤差が大きい。ddtSchemes を backward（BDF2、2次精度）に変えれば精度は向上する。CFDソルバーの時間スキームのデフォルトがFEAソルバーと異なることを認識しておくべきだ。
時間積分精度の比較
時間積分スキームによる差を定量的に見たいです。
Abaqusで同一メッシュ、$\Delta t = 1.6$ s（粗い刻み）での比較だ。
| 時間スキーム | T(x=0.08, t=32) [°C] | 参照値との差 [°C] |
|---|---|---|
| 後退Euler (α=0) | 35.8 | -0.80 |
| Crank-Nicolson (α=0.5) | 36.7 | +0.10 |
| Galerkin (θ=2/3) | 36.5 | -0.10 |
Crank-Nicolsonが一番良いですか？
2次精度だから$\Delta t$が大きいときの精度は高い。ただしステップ変化の初期に振動（オーバーシュート）が出る欠点がある。実務では最初の数ステップを後退Eulerで計算し、その後Crank-Nicolsonに切り替える手法が使われる。Abaqusではこの切り替えが自動で行われる。
多次元への拡張
2次元・3次元の熱伝導にも理論解がありますか？
直交座標で各方向が独立なら、変数分離法で1D解の積として2D/3D解を構成できる。例えば矩形領域の4面をステップ加熱する場合、各方向のerfc解の積になる。円筒座標や球座標ではBessel関数やルジャンドル関数を含む解になる。
NAFEMS T2（2次元定常熱伝導）やT4（軸対称非定常）が多次元検証のベンチマークだ。1D検証をパスした後、段階的に次元を上げていくのがV&Vの正攻法だ。
ツール選定の直感的ガイド
ツール選びのたとえ
CAEツールの選定は「道具箱」の構築に似ている。1つの万能ツールですべてをカバーするか、用途ごとに最適な専用ツールを揃えるか——予算、スキル、使用頻度に応じた戦略が必要。
選定で最も重要な3つの問い
- 「何を解くか」：1次元非定常熱伝導（半無限体）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
検証データの視覚化
理論値と計算値の比較を定量的に示す。誤差5%以内を合格基準とする。
| 評価項目 | 理論値/参照値 | 計算値 | 相対誤差 [%] | 判定 |
|---|---|---|---|---|
| 最大変位 | 1.000 | 0.998 | 0.20 | PASS |
| 最大応力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有振動数（1次） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合計 | 1.000 | 1.001 | 0.10 | PASS |
| エネルギー保存 | 1.000 | 0.999 | 0.10 | PASS |
判定基準: 相対誤差 < 1%: ■ 優良、1〜5%: ■ 許容、> 5%: ■ 要検討
V&V検証の効率化は、シミュレーションの信頼性を支える基盤です。 — Project NovaSolverは検証プロセスの改善にも注力しています。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「1次元非定常熱伝導（半無限体）をもっと効率的に解析できないか？」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →