ステファン問題 — CAE用語解説

厳密には、移動境界を伴う相変化問題の古典的な数学モデルを指します。例えば、氷が溶けて水になる際、固相と液相の境界（相界面）が時間とともに移動します。この移動境界の位置自体が未知数であり、その位置で潜熱の吸収・放出が起こるという条件（ステファン条件）が課されるのが特徴です。

支配方程式自体は各相内で通常の非定常熱伝導方程式、例えば

良い指摘です。古典的なステファン問題は、純物質の等温相変化を仮定します。しかし実務では、合金の凝固のように「凝固温度範囲」を持つ問題がほとんどです。この場合、固相率を導入する「エンタルピー法」や「温度回復法」といった近似解法が用いられます。例えば、アルミニウム合金A356の凝固範囲は約575°Cから615°Cです。この温度域では固相と液相が混在する「マッシューゾーン」が存在し、境界がシャープではなくなります。

主に2つのアプローチがあります。1つは「界面追跡法」で、メッシュを界面に合わせて変形・再構成します。もう1つは「界面捕捉法」で、固定メッシュ上で界面の位置を補助変数（例えばレベルセット関数やVOF関数）で表現します。界面追跡法は精度が高い反面、大変形時のメッシュ品質悪化が課題です。鋳造シミュレーションでは、後者のアプローチをベースにしたソリューションが多く見られます。

そこで「エンタルピー法」がよく使われます。温度

します。これが大きな課題です。温度-エンタルピー関係が不連続に近いため、反復計算の中で「オーバーシュート」が起こり、1つのタイムステップでセルが過剰に融解・凝固して発振することがあります。これを防ぐために、「エンタルピーソース項法」や「温度回復法」といった手法が開発されています。実装の際は、潜熱の処理を「陰的」に行い、反復収束を慎重にモニタリングする必要があります。多くの商用ソルバーは、この安定化処理をデフォルトで組み込んでいます。

3点が重要です。第一に「潜熱」の値。JISやASMハンドブックなど信頼できるデータベースから取得してください。第二に「固相率-温度関係」。純物質でない限り、これは線形ではなく、Scheil式やレバーロール則に基づくデータが必要です。第三に「凝固・融解温度範囲での熱伝導率と比熱」。これらは固相率の関数として変化することが多く、定数だと誤差の原因になります。例えば、マグネシウム合金AZ91の凝固範囲は約470°Cから595°Cで、この間の物性は大きく変わります。

界面の位置が事前にわからないので、相変化が起こり得る領域全体に対して、ある程度細かいメッシュが必要です。経験則として、少なくとも凝固層の厚さ方向に5〜10要素は確保したい。例えば、10mmの凝固層を予想するなら、メッシュサイズは1〜2mm以下にします。また、温度勾配が急峻な金型近くや冷却水管周りは、さらに細かくする必要があります。メッシュ依存性の検証は必須です。粗いメッシュでは潜熱解放が早すぎて、凝固完了時間が実際より早く算出される傾向があります。

これが最も不確実性の高いパラメータの一つです。鋳造では、コーティングの有無、金型温度、圧力によって大きく変わります。例えば、ダイカストでアルミニウムを鋳造する場合、金型との界面熱伝達率は初期は非常に高く（10,000 W/m²K以上）、凝固収縮で空隙ができると急激に低下（1,000 W/m²K以下）します。実測データがない場合は、文献値や類似プロセスからの推定値を使い、感度解析を行って結果への影響を評価するのが現実的です。JFS規格（日本鋳造工学会）などにも参考データがあります。

いいえ、通常はより高レベルな機能を使います。Ansys Fluentでは「Solidification/Melting」モデルを有効にするだけで、背後でエンタルピー法が動きます。ユーザーは材料に潜熱と固液相線温度を定義するだけです。COMSOL Multiphysicsでは、「Heat Transfer Module」の「Phase Change」ノードを追加し、同様に物性を設定します。どちらも古典的なステファン条件を直接入力するインターフェースはなく、実用上の近似解法が実装されています。

専用ソフトは、鋳造プロセスに特化した物理モデルと前後処理が充実しています。例えば、ProCAST（ESIグループ）は、収縮孔の予測に必要な「二相域でのダルシー流れモデル」や、微細組織予測の「CAFE法」を標準で備えています。また、金型の加熱サイクルや鋳湯の注入プロセスを考慮した連続鋳造のシミュレーションも得意です。一方、Fluentなどの汎用ソフトは、凝固・融解そのものの基礎プロセスは計算できますが、鋳造特有の現象を再現するにはユーザーが追加モデルを構築する必要があり、難易度が高くなります。

Abaqus/Standardでは、「熱電気」ステップで温度場を計算した後、「構造」ステップで熱歪みを負荷する順連成解析が一般的です。相変化に伴う体積変化（例えば、鋳鉄のグラファイト化膨張）は、ユーザーサブルーチン「UEXPAN」や「UMAT」を用いて、温度と相（固相率）の関数として熱歪みを定義することで表現できます。あるいは、Abaqusの「Model Change」機能を使って、凝固した要素の物性（熱膨張係数）を動的に変化させる方法もあります。いずれにせよ、純粋なステファン問題を超えた、材料工学の知識が必要になります。

最も多い原因は「タイムステップが大きすぎる」ことです。潜熱の解放・吸収は非常に大きな熱源項であり、これを陽的に扱うと不安定になります。対策としては、まずタイムステップを小さく（例えば0.01秒から0.001秒に）して試してください。第二に、ソルバーの「緩和係数」を確認します。温度やエンタルピーの緩和係数をデフォルト（1.0）から0.5や0.7などに下げると安定化することがあります。Ansys Fluentでは、「Solidification/Melting Parameters」で「Mushy Zone Parameter」を調整（デフォルト10^5を10^6〜10^7に増加）することで、凝固領域の擬似的な粘性を上げ、発散を防ぐ方法もあります。

第一に「潜熱の値」と「固液相線温度」を再確認してください。データシートの値は純度99.99%のような高純度物質のもので、実際の合金とは異なる場合があります。第二に「界面熱伝達率」。これが小さすぎると熱が逃げにくく凝固が遅く、大きすぎると早くなります。第三に「メッシュサイズ」。先述のように粗すぎると熱拡散が早く計算され、凝固が早く完了します。メッシュ依存性の検証を行い、実測値と合うように界面熱伝達率を逆算して調整（キャリブレーション）する作業が、実務ではよく行われます。

物性、特に「固相率-温度関係」の設定が不適切な可能性が高いです。純物質の設定（固相率が融点で0から1に瞬間変化）をしているのに、数値拡散の影響でマッシューゾーンが広がって見えることもあります。合金の場合は、平衡状態図に基づく正確な固相率曲線を入力してください。また、熱伝導率を固相率の関数として定義していないと、マッシューゾーン内の熱伝達が不正確になり、ゾーンが不自然に広がる原因になります。ソフトウェアのマニュアルで、固相率と物性を連動させる設定方法を確認しましょう。