凝固・鋳造シミュレーション
理論と物理
概要
先生、凝固のCFDって何を計算するんですか？
溶融金属が冷却されて固化する過程を予測する。鋳造の充填・凝固解析、連続鋳造のシェル成長、溶接のビード形成、金属3Dプリンティング（PBF/DED）の溶融池挙動など、固液相変化を伴う流動・伝熱問題だ。
液体が固体になるのを流体力学で扱うんですか？
エンタルピー法（Enthalpy-Porosity法）が標準的な手法だ。固液界面を明示的に追跡せず、各セルの液相分率 $f_L$（0〜1の連続変数）で凝固の進行度を表現する。
支配方程式
エンタルピー法の方程式を教えてください。
エネルギー方程式はエンタルピー $H$ で記述する。
エンタルピーは顕熱と潜熱の和だ。
$h = \int_{T_{ref}}^{T} c_p dT$ は顕熱、$f_L$ は液相分率、$L$ は凝固潜熱だ。液相分率はリキダス温度 $T_L$ とソリダス温度 $T_S$ から決まる。
固化した部分の流れはどうなるんですか？
Enthalpy-Porosity法では、固化した領域を高抵抗の多孔質体として扱う。運動量方程式にDarcyの抗力ソース項を追加する。
$C$ はMushy zone定数（典型値$10^5$〜$10^8$）、$\epsilon$ はゼロ除算回避の小定数、$\mathbf{u}_{pull}$ は引き抜き速度（連続鋳造の場合）だ。$f_L \to 0$（完全固化）でソース項が無限大に近づき、速度がゼロに減衰する。
マランゴニ対流
溶融池の流れにはどんな駆動力がありますか？
表面張力の温度依存性（$d\sigma/dT$）によるマランゴニ対流が溶融池の流れパターンを支配する。
$$ \tau_{Ma} = \frac{d\sigma}{dT} \nabla_s T $$
多くの金属で $d\sigma/dT < 0$（高温ほど表面張力が低い）なので、溶融池中央から外向きの流れが発生する。ただし酸素やイオウの含有量で符号が反転し、流れパターンが劇的に変わることがある。
Coffee Break よもやま話
凝固の微視的世界——デンドライト成長とStephan問題
金属凝固の物理は「固液界面がどのように移動するか」というStephan問題として定式化されます。純金属の場合、界面は融点等温面に沿って移動し、その速度は潜熱と熱流束のバランスで決まります。しかし実際の合金では局所的に組成過冷が生じ、樹枝状（デンドライト）の複雑な固液界面が自発的に形成されます。このデンドライト成長をCFDで直接計算するPhase-Field法は1990年代以降急速に発展し、2020年代には10^6メッシュのデンドライト計算がGPUで実時間に近い速度で実行可能になりました。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
変数 SI単位 注意点・換算メモ 速度 $u$ m/s 入口条件で体積流量から換算する際、断面積の単位に注意 圧力 $p$ Pa ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 密度 $\rho$ kg/m³ 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C 粘性係数 $\mu$ Pa·s 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 レイノルズ数 $Re$ 無次元 $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 CFL数 無次元 $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結
溶融池の流れにはどんな駆動力がありますか？
表面張力の温度依存性（$d\sigma/dT$）によるマランゴニ対流が溶融池の流れパターンを支配する。
多くの金属で $d\sigma/dT < 0$（高温ほど表面張力が低い）なので、溶融池中央から外向きの流れが発生する。ただし酸素やイオウの含有量で符号が反転し、流れパターンが劇的に変わることがある。
凝固の微視的世界——デンドライト成長とStephan問題
金属凝固の物理は「固液界面がどのように移動するか」というStephan問題として定式化されます。純金属の場合、界面は融点等温面に沿って移動し、その速度は潜熱と熱流束のバランスで決まります。しかし実際の合金では局所的に組成過冷が生じ、樹枝状（デンドライト）の複雑な固液界面が自発的に形成されます。このデンドライト成長をCFDで直接計算するPhase-Field法は1990年代以降急速に発展し、2020年代には10^6メッシュのデンドライト計算がGPUで実時間に近い速度で実行可能になりました。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
凝固シミュレーションの数値的なポイントを教えてください。
最大の課題は凝固潜熱の取り扱いだ。潜熱の放出はエンタルピーの非線形性を生み、温度場の急激な変化を引き起こす。タイムステップの管理と反復計算の収束が鍵になる。
Mushy zone定数の影響
Mushy zone定数 $C$ はどう設定するんですか？
$C$ の値は結果に大きな影響を与える。
| $C$ の値 | 効果 | 用途 |
|---|---|---|
| $10^4$〜$10^5$ | 緩やかなダンピング | 低速凝固、連鋳 |
| $10^5$〜$10^6$ | 標準的 | 一般的な鋳造 |
| $10^7$〜$10^8$ | 急速なダンピング | 急速凝固、溶接 |
理想的にはDarcy透過率の実験データから決めるべきだが、実務的には $10^5$〜$10^6$ から始めて感度解析を行うことが多い。
VOF法との連成
鋳造の充填過程はどうやって計算するんですか？
充填過程はVOF法で溶湯の自由表面を追跡し、充填完了後（または充填と同時に）凝固を計算する。FluentやFlow-3DではVOF + Solidification/Meltingモデルの同時計算が可能だ。
充填中の湯流れパターン、酸化膜の巻き込み、ガス巻き込みが鋳造欠陥に直結するため、充填解析は鋳造シミュレーションの重要な要素だ。
ツール別の実装
| ツール | 凝固モデル | VOF連成 | 応力解析連携 |
|---|---|---|---|
| Ansys Fluent | Solidification/Melting | 対応 | Mechanical連携 |
| Flow-3D | TruVOF + Solidification | ネイティブ | 限定的 |
| STAR-CCM+ | Solidification | 対応 | 構造連成 |
| ProCAST (ESI) | FEM凝固 + 充填 | 専用充填 | 応力・変形解析統合 |
| MAGMASOFT | 独自FVM | 専用モデル | 残留応力予測 |
ProCASTとMAGMASOFTは鋳造専用ソフトですよね？
その通り。鋳造プロセスに特化した専用ツールで、充填・凝固・収縮巣予測・残留応力・変形まで一貫して計算できる。汎用CFDより鋳造固有の物理モデル（例：フィーディング、Niyamaクライテリア）が充実している。
Enthalpy-Porosity法——実用凝固CFDの標準アプローチ
工業的な凝固CFD（連続鋳造、鋳造充填）で最も広く使われるのがEnthalpy-Porosity（エンタルピー-多孔質）法です。固液共存域（mushy zone）を「浸透率が小さい多孔質」として扱い、Kozeny-Carman式で局所速度を抑制することで固化の進行を再現します。このアプローチはVoller & Prakkashが1987年に提案し、ANSYS FluentやOpenFOAMのsolidificationMeltingSourceに実装されています。モデルのキャラクタライゼーション定数（Amush）はデフォルト値1.6x10^5が使われることが多いですが、この値が不適切だと凝固速度が30〜50%ずれます。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
凝固シミュレーションの手順を教えてください。
アルミニウム重力鋳造の解析を例に説明しよう。
1. 形状作成: 製品 + 湯道系 + 押湯 + 鋳型
2. メッシュ: 製品部は2〜5 mm、薄肉部は3層以上確保
3. 物性値: 温度依存の密度・比熱・熱伝導率・粘度、凝固潜熱
4. 充填計算: VOF法で溶湯注入（注湯温度700℃程度）
5. 凝固計算: Solidification/Meltingモデルで冷却・凝固
6. 界面熱抵抗: 鋳型-溶湯間のGap熱伝達係数を設定
7. 後処理: 凝固時間分布、最終凝固位置、Niyama値
鋳造欠陥の予測
鋳造欠陥もCFDで予測できるんですか？
主な欠陥予測指標を示そう。
| 欠陥 | 予測指標 | 方法 |
|---|---|---|
| 引け巣（Shrinkage porosity） | Niyama criterion $Ny = G/\sqrt{\dot{T}}$ | $Ny < 1$ で引け巣リスク |
| ガスポロシティ | 水素溶解度 vs 局所圧力 | Sievert's法則 |
| 湯境 / コールドシャット | 湯先温度分布 | 温度 < リキダスの領域 |
| 酸化膜巻き込み | VOFの液面折れ畳み | 自由表面の追跡 |
| 凝固収縮 | $\Delta V / V = (\rho_s - \rho_l)/\rho_s$ | 体積変化の計算 |
Niyamaクライテリアって何ですか？
Niyama et al.（1982）が提案した引け巣予測のための無次元基準だ。温度勾配 $G$（K/m）と冷却速度 $\dot{T}$（K/s）の比で定義される。鋼鋳物では $Ny < 1$ (K·s)$^{1/2}$/mm で引け巣リスクが高い。鋳造専用ソフト（ProCAST, MAGMASOFT）にはNiyama値の自動計算機能が組み込まれている。
物性値の重要性
凝固シミュレーションで特に重要な物性は何ですか？
凝固潜熱 $L$ と液相分率-温度関係 $f_L(T)$ が最重要だ。合金の場合はScheil式やLever Ruleで $f_L(T)$ を計算する。JMatProやThermo-Calcなどの計算状態図ソフトウェア（CALPHAD）から物性データを取得するのが理想的だ。
連続鋳造のCFD——鉄鋼生産の「心臓部」を守るシミュレーション
連続鋳造（Continuous Casting）は世界の粗鋼生産の95%以上が通る工程で、溶鋼が鋳型内で急速冷却されながら連続的に引き出されます。鋳型内の流動と凝固のCFD解析は、表面欠陥（スカム巻き込み、割れ）の発生機構の理解と工程最適化に不可欠です。日本製鉄は1990年代からCFDを用いた鋳型内流動制御（電磁攪拌）の設計を行っており、インクルージョン（酸化物介在物）の捕捉率を25%低減した事例を学術論文で公表しています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
凝固・鋳造シミュレーションのツールを比較してください。
| ツール | 充填 | 凝固 | 応力 | 欠陥予測 | 特徴 |
|---|---|---|---|---|---|
| MAGMASOFT | 独自FVM | 対応 | 対応 | Niyama, Porosity | 鋳造専用、最も広く使用 |
| ProCAST (ESI) | FEM | 対応 | 対応 | 多数の基準 | 合金物性統合(CALPHAD) |
| Flow-3D CAST | TruVOF | 対応 | 限定的 | Defect tracking | 自由表面精度が高い |
| Ansys Fluent | VOF | Enthalpy-Porosity | Mechanical連携 | UDFで拡張可能 | 汎用CFDの柔軟性 |
| STAR-CCM+ | VOF | 対応 | 構造連成 | 基本指標 | 汎用CFD |
| SOLIDCast | 独自 | 対応 | 限定的 | FCC/Gating | 低コスト・簡易 |
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 砂型鋳造（鋳鉄、鋼） | MAGMASOFT, ProCAST | 鋳造専用機能の充実 |
| ダイカスト | MAGMASOFT, Flow-3D CAST | 高速充填の精度 |
| 連続鋳造 | Fluent, STAR-CCM+ | カスタマイズの柔軟性 |
| 溶接 | Fluent, STAR-CCM+ | マランゴニ対流、レーザー熱源 |
| 金属3Dプリンティング | Fluent, Flow-3D AM | 溶融池 + 凝固 + 残留応力 |
| 学術研究 | OpenFOAM | ソースコードアクセス |
鋳造専用ソフトと汎用CFDはどう使い分けるんですか？
鋳造プロセスの最適化（湯道設計、押湯配置、冷し金位置）なら専用ソフトが圧倒的に効率的だ。鋳造固有のワークフロー、材料データベース、欠陥予測基準が統合されている。一方、新しい物理モデルの開発や非定型的な問題（溶接、AM、結晶成長など）では汎用CFDの柔軟性が必要だ。
Flow-3D AM（Additive Manufacturing）は金属3Dプリンティング専用モジュールで、レーザー出力・走査速度と溶融池形状・凝固組織の関係を予測できる。AM分野で急速に採用が進んでいる。
MAGMAsoft vs ProCAST——鋳造CFDの二大巨頭
鋳造プロセスCFDの分野ではMAGMAsoft（MAGMA社）とProCAST（ESI Group）が世界シェアの大半を占めます。MAGMAsoftは工程エンジニアが使いやすいGUIと豊富な合金データベース（1000種以上）が強みで、自動車アルミダイカスト業界では事実上の標準です。ProCASTはFEM（有限要素法）ベースで熱-構造連成（残留応力・変形予測）に強く、航空機ターボブレードの精密鋳造（investment casting）で多く採用されています。両ツールとも年間ライセンスが数百万円規模で、専用UIと合金データベースの充実度が選ばれる理由です。
選定で最も重要な3つの問い
- 「何を解くか」：凝固・鋳造シミュレーションに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
凝固シミュレーションの最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
Cellular Automaton-Finite Element（CAFE）法
凝固組織（デンドライト構造、柱状晶/等軸晶の遷移CET）を予測するCAFE法が発展している。FEMで温度場を解き、Cellular Automatonで結晶核生成と成長を追跡する。ProCASTに統合されており、鋳造組織の予測が可能だ。
Phase-Field法による凝固組織シミュレーション
Phase-Field法って何ですか？
秩序パラメータ $\phi$（固相=1、液相=0）の時間発展方程式で界面を暗示的に追跡する手法だ。デンドライトの樹枝状成長を直接シミュレーションできる。
デンドライトアーム間隔やマイクロ偏析を予測できるが、計算コストが膨大で、現在は数百μm〜mm スケールに限定されている。
金属3Dプリンティングのマルチスケールシミュレーション
金属AM（Additive Manufacturing）では、溶融池スケール（CFD）→凝固組織（Phase-Field/CA）→部品スケール（FEM構造解析）のマルチスケール連成が重要な研究テーマだ。レーザーのパラメータ（出力、速度、ハッチ間隔）と組織・残留応力の関係を予測する。
3DプリンティングにCFDが使われているんですね。
PBF（粉末床溶融）やDED（指向性エネルギー堆積）の溶融池内のマランゴニ対流、キーホール形成、スパッタ飛散をVOF法で直接計算する研究が急増している。Khairallah et al.（LLNL, 2016）のALE3D計算が先駆的だ。
機械学習による鋳造条件最適化
CFDの凝固シミュレーション結果をデータとして、湯道設計や鋳造条件の最適化を機械学習で行う研究が進んでいる。数百〜数千ケースのCFD結果からサロゲートモデルを構築し、リアルタイムの条件探索に活用する。
Phase-Field法——凝固CFDのフロンティア
Phase-Field法は固相・液相・共存域を「相場φ（0〜1）」という連続的なスカラー場で表現し、界面を暗的に追跡する手法です。界面位置を陽的に追跡する必要がなく、トポロジー変化（デンドライトの枝分かれ）を自然に扱えます。金属合金凝固でのPhase-Field法は、Karma & Rappelの薄界面漸近展開（1996年）により定量的予測が可能になりました。東北大学と大阪大学を中心とした日本のグループは、Al-Si合金のデンドライト組織をPhase-Field-CFD連成で予測し、放射光X線その場観察との直接比較で形状を再現することに成功しています。
トラブルシューティング
トラブルシューティング
凝固シミュレーションでよくあるトラブルを教えてください。
順番に見ていこう。
1. 凝固が始まらない / 進まない
対策:
- リキダス温度 $T_L$ とソリダス温度 $T_S$ が正しく設定されているか確認
- 凝固潜熱 $L$ の値と単位系を確認
- 鋳型の境界条件（熱伝達係数、鋳型温度）が適切か確認
- Mushy zone定数が小さすぎないか確認
2. Mushy zoneでの数値振動
凝固前線付近で温度が振動します…
対策:
- タイムステップを小さくする（凝固の潜熱放出を十分に解像）
- エネルギー方程式の under-relaxation を下げる（0.7〜0.9）
- メッシュをMushy zone領域で細かくする
- Mushy zone定数 $C$ を調整（大きすぎると速度の急変で不安定化）
3. 凝固後に非物理的な流れが残る
症状: 完全固化（$f_L = 0$）した領域でまだ速度が残る。
対策:
- Mushy zone定数 $C$ を大きくする（$10^7$以上）
- $\epsilon$ の値が大きすぎないか確認（$10^{-3}$程度が標準）
- 固相の粘度が十分に高くなっているか確認
4. 充填+凝固の同時計算で発散
対策:
- 充填計算と凝固計算を分離する（まず充填、完了後に凝固）
- VOF + Solidificationの同時計算では、タイムステップを十分小さく
- 充填初期に溶湯温度を高めに設定して凝固の開始を遅らせる
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Solidification/Melting ModelのPull Velocity設定（連鋳の場合のみ） |
| Flow-3D | Solidification Drag Coefficientの調整がキー |
| ProCAST | CALPHAD物性データの温度範囲がシミュレーション範囲をカバーしているか確認 |
| MAGMASOFT | メッシュ自動生成の品質確認。薄肉部の層数が不足していないか |
ショートサーキット凝固——湯廻り不良の診断とCFD
鋳造充填CFDで発生する「ショートサーキット（短絡凝固）」は、溶湯の流路が早期凝固によって閉塞し、鋳型の一部に湯が届かない現象です。CFDで見逃しやすいのは薄肉部への流動末端部で、温度降下が速すぎてシミュレーション初期から固化している場合があります。対策として初期壁温を実測値に設定し、Heat Transfer Coefficient（HTC）の感度解析を行うことが有効です。HTCは金型材料・離型剤種・接触圧力によって10〜1000 W/m^2Kの幅があり、この不確かさを±50%の範囲で変えてシミュレーションした場合の充填率変化を確認することが、設計余裕評価の実務的アプローチです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——凝固・鋳造シミュレーションの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告