Fluentの凝固・融解モデルと理論解の比較
By J.Aono
はじめに
Ansys Fluentでは，凝固や融解が（純金属の凝固・融解のように）特定の温度で，発生する流体流れの問題を解析することができます．Ansys Fluentでは，固液の接触領域を明示的に追跡する代わりに，エンタルピー・空隙率法を使用しています．ここでは，氷が融解して水に変化する時間精度を確認します．
Stefan問題
相変化の時間精度を調べるには，Stefan問題がよく知られていて，Alexiadesら1らの理論解があります．この問題は１次元の問題で，相の界面の位置を調べる問題です．
相の界面は次の式で表されます． $$ x_{\Gamma}(t) = 2\gamma\sqrt{\alpha_l t} \tag{1} $$
ここで，$x_{\Gamma}$ は相の界面位置，$\alpha$ は温度拡散率， $t$は時間，$\gamma$ は一定値，下添字$l$は液相です．$\gamma$は次の式の超越方程式を満たします．
$$ \gamma e^{\gamma^2} erf(\gamma) = \frac{C_{p_l} (T_{wall}-T_{melt})}{\hat{H}_{LV} \sqrt{\pi}} \tag{2} $$
ここで，$\hat{H}_{LV}$ は潜熱，$C_p$は定圧比熱，下添字$melt$は凝固点です．
温度分布は次のように表すことができます．
$$ \begin{align} T_l(x,t)=T_{wall} + \frac{T_{melt}-T_{wall}}{erf(\gamma)} erf(\frac{x}{2\sqrt{\alpha_l t}}) \tag{3} \end{align} $$
ここで，$T$ は温度，$\alpha$ は温度拡散率，$\hat{H}_{LV}$ は潜熱，$C_p$は定圧比熱です．この理論解とAnsys Fluentの結果を比較します．
Ansys Fluentの融解・凝固モデル
融解・凝固モデルは弊社記事の鋳造・凝固シミュレーションに詳細がありますので，ここでは割愛します．
Ansys Fluentを使った数値計算
ここでは氷が溶けて水になる数値計算を考えます．固相線温度$T_{solidius}=272.15 K$とし，液相線温度$T_{liquidus}=273.15 K$とします．エンタルピー・空隙率法の半溶融領域の定数はデフォルトの$A_{mush}=100,000$としました．計算モデルは下記のように$ 10 mm$の氷を配置し，左側に氷が融解する温度$T_{wall}=273.15 K$を与えます．
初期温度は$272.15 K$とし，氷とします．
計算格子は100セルの等分割としました．
初期に氷を配置し，左側の壁は$273.15 K$の等温壁を設定します．氷が水になる融解温度は$273.15 K$とします．
次に氷と水の体積分率の計算結果をアニメーションで示します．青色が氷を示し，赤色が水を示します．左側から氷が融解し，水に変化していることが確認できます．
理論解とFluentの計算結果の比較
Alexiadesら1らの理論解とAnsys Fluentの計算結果を比較します．白色から青色のカラーマップは式(3)の結果です．黒い実線が式(1)の結果です．Ansys Fluentの結果は氷と水の体積分率が$0.5$を相界面とし，赤色の丸のマーカーで示しています．理論解と良く一致していることが確認できました．
おわりに
今回は，Ansys Fluentの凝固・融解モデルの時間精度について，理論解と比較しました．半溶融領域の定数や固相・液相線温度のパラメータがありますが，氷と水の場合は今回の設定で理論解と良い一致を得ることができました．
弊社の解析事例
弊社の流体解析事例については，下記のリンクからご覧ください．