マランゴニ対流 — CAE用語解説

良い質問だ。浮力対流が体積力（重力）に起因するのに対し、マランゴニ対流は界面に働く表面張力の不均一が駆動力だ。具体的には、溶融プールの溶接シミュレーションで顕著だな。溶接部の中心は温度が高く、周辺は低い。表面張力は温度が高いほど小さくなる（多くの金属や溶融ガラスで

その通り。無重力環境では浮力対流が抑制されるため、マランゴニ対流が支配的になる。実際、国際宇宙ステーション（ISS）での材料実験（例えばJAXAの「MELT」実験）では、溶融材料中のマランゴニ対流の挙動が重点的に研究されている。地上では両者が競合するが、微小重力下では純粋なマランゴニ効果を観測できる。支配方程式では、境界条件として表面張力の効果を考慮する必要がある。

体積力項ではなく、自由表面における応力境界条件として組み込まれる。具体的には、自由表面におけるせん断応力の連続条件が、表面張力勾配と釣り合う形で記述される。温度勾配が原因の場合、その境界条件は以下のように書ける。

その通りだ。少なくとも3点が重要だ。第一に、表面張力係数を「定数」ではなく「温度の関数」として定義する必要がある。多くのソフトウェアでは、線形関数

マランゴニ数（Ma）やカプリ数（Ca）が指標になるが、実用的には境界層を捉えることが重要だ。表面張力勾配によって生じる速度境界層の厚さ

温度場と流れ場が強く結合するため、分離解法より完全陰的連成解法（Coupled solver）の収束性が良い場合が多い。Ansys Fluentで言えば「Pressure-Based Coupled Algorithm」を使う。また、表面張力の効果は明示的に扱うと時間ステップ制限（CFL条件）が非常に厳しくなるため、CSF（Continuum Surface Force）モデルなどを用いた陰的な取り扱いが標準だ。収束が難しい場合は、まずマランゴニ効果なしで流れ・温度場を初期化し、その後マランゴニ効果を徐々に立ち上げる（ランプする）手法も有効だ。

まず、材料物性の入力が全ての基礎だ。最低限必要なのは、密度、粘度、熱伝導率、比熱に加え、表面張力の温度係数

無次元数の比、ボンド数（Bo）を見るのが一般的だ。浮力と表面張力勾配力の比で定義される。

直接流速を測るのは難しいので、マランゴニ対流が及ぼす二次的な効果を測定値と比較する。最も一般的なのは「溶融プール形状」だ。マランゴニ対流はプール内の熱輸送を促進するため、プールがより深く、狭くなる傾向がある。また、添加材が分散する「マクロ偏析」のパターンも重要な比較対象だ。シミュレーションでは、温度分布と流速ベクトルを可視化し、表面から中心に向かう流れ（温度係数が負の場合）が形成されているかを確認せよ。定量的には、プール最深部の温度や、特定点の冷却曲線を実験データと比較する。

コアとなる物理モデルは同じだが、インターフェースと設定の柔軟性が大きく異なる。Ansys Fluentでは、「Multiphase」モデル（VOF）を有効にした上で、「Surface Tension」モデルを有効にし、「Marangoni Stress」にチェックを入れる。表面張力係数は「constant」または「piecewise-linear」で温度関数を定義できる。ユーザー定義関数（UDF）を使えばより複雑な関数も可能だ。ただし、浮力とマランゴニの連成を安定して解くには、かなり経験を要する。

「CFD Module」と「Microfluidics Module」のいずれかが必要だ。設定は直感的で、「Two-Phase Flow, Level Set」インターフェースなどを追加し、「Surface Tension」設定で「Include Marangoni effect」をオンにする。ここで表面張力係数を温度の関数として直接入力できる。COMSOLの強みは、熱伝導や化学種輸送との連成がGUI上でシームレスに行える点だ。例えば、溶質濃度勾配によるマランゴニ効果（溶質マランゴニ効果）も比較的簡単に設定できる。

標準の多相流ソルバー（例えば`interFoam`）には表面張力モデル（CSF）は実装されているが、マランゴニ効果はデフォルトでは含まれていない。しかし、`interFoam`の派生ソルバーである`interThermalFoam`や、コミュニティで開発されたソルバー（`reactingInterFoam`を改造したもの等）を使うか、自分で境界条件を実装する必要がある。具体的には、`alphaEqn.H`や`UEqn.H`を修正して、表面張力係数を温度の関数とし、その勾配を表面せん断応力として境界条件に加えるコードを追加する。柔軟性は最高だが、ハードルも高い。

確かにある。例えば、結晶成長シミュレーションの分野では、CGSim (STR Group) というソフトウェアが広く使われている。これはシリコン単結晶引き上げ法（Cz法）などのプロセス最適化が目的で、浮力、マランゴニ、遠心力などの効果を高度にモデル化している。また、ガラス溶融炉の設計には、Glass Service 社の `GlasSim` などがある。これらの専用ソフトは、産業特有の物理現象と材料物性データベースが最初から組み込まれており、汎用CFDソフトより設定が容易で信頼性が高いが、当然ながら適用範囲は限定される。

最も多い原因は3つだ。第一に、表面張力係数の温度関数が極端な値を持っている。例えば、計算領域内の温度範囲で表面張力が負の値になっていないか確認せよ。第二に、初期条件が悪い。静止状態から突然マランゴニ効果をオンにすると、大きなせん断応力が発生して不安定になる。まずはマランゴニ効果なしで熱対流をある程度発達させ、その後マランゴニ効果をランプ関数でゆっくりと立ち上げる。第三に、メッシュが粗すぎる。先述した表面境界層を解像できておらず、勾配の計算が暴れる。

それは「マランゴニ・ベナール対流」の可能性が高い。均一な加熱ではなく、下面加熱のような一様な温度勾配条件下では、マランゴニ効果は表面張力の不安定性を引き起こし、規則的なセル状の対流パターン（ベナールセル）を形成する。しかし、溶接のような局所加熱でそれが起きるなら、注意が必要だ。まず、境界条件や初期条件の対称性が保たれているか確認せよ。数値的な丸め誤差が非対称性の種になることがある。また、メッシュが非対称だと、それだけで非対称な解にバイアスがかかる。物理的に非定常・乱流に遷移しているのかもしれない。時間平均を取ると対称性が回復するか確認すること。

非常に感度が高いパラメータだ。符号が逆転すれば流れの方向が逆になる。絶対値も流速にほぼ比例して影響する。例えば、溶融スズの溶接シミュレーションで、

界面の振動（パラサイトカレント）はVOF法の宿敵で、マランゴニ計算では特に致命的だ。対策はいくつかある。第一に、界面再構築スキームを「Geo-Reconstruct」や「HRIC」のような高精度なものにすること。第二に、表面張力のモデル化を見直す。Fluentでは「Surface Tension Model」で「Implicit」を選択し、「Curvature Specification」を「Curvature from Cell Gradient」よりも安定性の高い「Brackbill et al.」などに変えてみる。第三に、根本的な解決策として、界面追跡法をLevel Set法に切り替えることを検討する。COMSOLではLevel Setがデフォルトのことも多く、界面の滑らかさという点ではVOFより優れている場合がある。