› ISS 微小重力 マランゴニ対流・液滴 戻る
宇宙実験・微小重力

ISS 微小重力 マランゴニ対流・液滴シミュレーター

国際宇宙ステーション（ISS）のような微小重力環境では浮力対流が消え、表面張力勾配が駆動する「マランゴニ対流」が主役になります。温度差・流体物性・特性長・重力レベルを変えて、マランゴニ数 Ma と動的 Bond 数を直感的に把握できる無料ツールです。

パラメータ設定
流体
密度・熱膨張係数・表面張力温度勾配 (dσ/dT) を自動設定
温度差 ΔT
K
表面張力 σ
mN/m
参考表示用（駆動力は dσ/dT を使用）
粘度 μ
mPa·s
特性長 L
mm
液柱長または液滴径
重力レベル
微小重力／月／火星／地球を比較
計算結果
—
マランゴニ数 Ma
—
レイリー数 Ra
—
比率 Ma/Ra
—
支配対流
—
表面流速 (mm/s)
—
動的 Bond 数
—
微小重力液滴 — 表面張力流の可視化

温度勾配（青→赤）による表面張力差が液滴表面に流れを生み、内部にロール状の循環ができます。微小重力では液滴は球形のまま浮遊します。

マランゴニ数 Ma と温度差 ΔT
流体別マランゴニ数比較（現在の ΔT・L・μ）
理論・主要公式

$$\mathrm{Ma} = \frac{|\partial\sigma/\partial T|\,\Delta T\,L}{\mu\,\alpha}, \qquad \mathrm{Ra} = \frac{g\,\beta\,\Delta T\,L^{3}}{\alpha\,\nu}$$

マランゴニ数 Ma とレイリー数 Ra。∂σ/∂T：表面張力の温度勾配 [N/m/K]、α：熱拡散率 [m²/s]、β：体積膨張係数 [1/K]、ν=μ/ρ：動粘度 [m²/s]。

$$\mathrm{Bo}_d = \frac{\mathrm{Ra}}{\mathrm{Ma}} = \frac{\rho\,g\,\beta\,L^{2}}{|\partial\sigma/\partial T|}$$

動的 Bond 数。Bo_d ≪ 1 でマランゴニ対流支配、Bo_d ≫ 1 でブイヤンシー対流支配。微小重力では g が約 10⁻⁶ になるため Bo_d が 6 桁減少し、マランゴニ対流が浮上する。

$$u_{\max} \sim \frac{|\partial\sigma/\partial T|\,\Delta T\,L}{\mu}$$

表面最大流速のスケーリング。温度差と特性長に比例し、粘度に反比例する。

ISS 微小重力 マランゴニ対流 — 表面張力駆動

🙋
先生、ISS の宇宙飛行士が水を球形のまま空中に浮かべている映像を見ました。あの「無重力の水」って中で何か動いているんですか？それとも完全に静止しているんでしょうか？
🎓
ぱっと見は静かそうだけど、実は中で結構ダイナミックな流れがある場合が多いんだ。微小重力では「下から温められて上に上がる」という普通の対流＝ブイヤンシー（浮力）対流はほぼ消える。でも液滴の表面に少しでも温度ムラがあると、表面張力の差で表面が動き出す。これが「マランゴニ対流」。地上では浮力対流が圧倒的に強くて気づかないけど、宇宙ではこれが主役になるんだ。
🙋
なぜ温度で表面張力が変わると流れができるんですか？表面張力って「水たまりが丸くなる力」というイメージしかなくて…
🎓
表面張力は「表面を縮ませようとする引っ張り合い」のことだ。水だと温度が上がるほど表面張力は下がる（dσ/dT が負）。たとえば液滴の片側が熱くて反対側が冷たいと、冷たい側の方が強く引っ張る。すると表面の水分子は「熱い側→冷たい側」へ引きずられる。表面が動けば内部も粘性で連れていかれて、ロール状の循環ができる。これがマランゴニ対流の正体だよ。
🙋
左の「重力レベル」を微小重力から地球に変えてみると、Ma/Ra 比がガクッと小さくなって「ブイヤンシー支配」に切り替わりますね。これって地上では本当にマランゴニは無視していいんですか？
🎓
いい質問！動的 Bond 数 Bo_d = Ra/Ma が決め手で、L²·g に比例する。だから「g が小さい（宇宙）」か「L が小さい（薄膜・小さな液滴）」と Ma が勝つ。地上でも、はんだ付けの溶融はんだボール、薄液膜の蒸発、結晶成長のメルトゾーンなんかは L が小さいからマランゴニが効く。「地上では完全に無視」ではなくて「大きなスケールでだけ無視できる」が正確だね。
🙋
流体プリセットを「溶融金属」にすると Ma がすごく大きくなって表面流速も速くなりました。これは何の研究と関係があるんですか？
🎓
溶融金属は dσ/dT が水の 2 倍以上ある上に、半導体結晶（Si、GaAs、Ge）を宇宙で育てる「浮遊帯溶融法」では液柱を電気炉で加熱する。地上だと自重で潰れて液柱を作れないけど、微小重力なら長い液柱が保てる。問題は、ここでマランゴニ対流が振動的になると結晶に縞模様（ストリエーション）が入って欠陥になる。Ma を 500〜1000 以下に抑えるための温度プロファイル設計が ISS 実験の主要テーマの一つなんだ。
🙋
なるほど、つまり「マランゴニ対流が強すぎる ＝ 結晶が荒れる」ということですね。じゃあ宇宙実験で得た知見は地上にも応用できるんでしょうか？
🎓
そう、ISS のマランゴニ実験（日本の「きぼう」モジュールの MEIS や FPEF、欧州の Geoflow）で得た「対流の遷移臨界条件」は地上のミニチュアスケール（マイクロ流体、インクジェット、レーザー溶接の微小溶融池）に直接活きる。Ma の閾値を超えるとロールが定常→振動→乱流と遷移する。これを CFD で予測する際の検証データを宇宙で取って、地上のプロセス制御に使う、というのが宇宙微小重力流体研究の本流だよ。

よくある質問

液体の表面張力が温度や濃度で変わるとき、表面に沿って張力が高い側へ流体が引っ張られて生じる対流のことです。一般的な液体では温度が上がると表面張力が下がるので、熱い領域から冷たい領域へ向けて表面流が走り、その下にロール状の循環ができます。地上では浮力（ブイヤンシー）対流の方が圧倒的に強いため目立ちませんが、ISS の微小重力では浮力が消えるためマランゴニ対流が支配的になります。本シミュレーターは温度差・物性・特性長から Ma を計算し、Ra との比 Ma/Ra で支配機構を判定します。
Ma = (∂σ/∂T)·ΔT·L / (μ·α) で定義されます。∂σ/∂T は表面張力の温度勾配（水で 0.15 mN/m/K 程度、負の値ですが絶対値を使うことが多い）、ΔT は温度差、L は特性長（液柱長や液滴径）、μ は粘度、α は熱拡散率です。Ma が大きいほど熱拡散より表面張力対流が速いことを意味し、おおむね Ma > 80 で対流が立ち上がり、Ma > 500〜1000 で振動的・乱流的になります。
本シミュレーターのデフォルト（水、ΔT=20K、L=20mm、微小重力）で Ma≈400,000 となります。これは地上の浮力対流が完全に支配する条件（Ra≈10〜100）と比べて 4 桁以上大きく、マランゴニ対流が完全に支配します。地上でも特に薄い液層や小さな液滴（L が小さい）では Ma が Ra を上回ることがあり、はんだ付けや結晶成長で重要になります。
動的 Bond 数 Bo_d = Ra/Ma は、浮力対流と表面張力対流の強さの比を表す無次元数です。Bo_d ≫ 1 ならブイヤンシー対流（地上の通常の対流）、Bo_d ≪ 1 ならマランゴニ対流（微小重力や微小スケール）が支配します。ISS のような微小重力では Bo_d が地上の 10⁻⁶ 倍になるので、地上ではマランゴニ対流が無視できる系でも宇宙では主役になります。

実世界での応用

宇宙での半導体結晶成長：シリコンやガリウム砒素の単結晶を微小重力で育てる「浮遊帯溶融法（フローティングゾーン）」では、液柱内のマランゴニ対流が結晶品質を直接左右します。Ma が大きいと対流が振動的になり、結晶に縞模様（ストリエーション）として欠陥が残ります。本ツールの Ma/Ra 比で、地上では浮力対流に隠れていた表面張力対流の寄与を切り分け、宇宙実験条件の事前検討に使えます。

ISS「きぼう」での流体実験（MEIS／FPEF）：JAXA の流体物理実験装置 FPEF は、シリコーン油の液柱を加熱して Ma の臨界値を測定する代表的実験です。本シミュレーターでシリコーン油を選び、ΔT を 10〜50K で振ると、Ma が 10⁴ オーダーに達し、地上ではあり得ない強いマランゴニ対流が発生することが分かります。

地上のはんだ付け・レーザー溶接：溶融はんだボールや溶接プールでは特性長 L が 1〜3mm と小さく、地上でも Ma が Ra を上回ります。「金属」プリセット＋L=2mm で計算すると、地上重力でもマランゴニ対流が支配的になることが確認できます。これがビードの形状や深さ（キーホール）に影響するため、実プロセス設計に直結する知見です。

マイクロ流体・インクジェット：µL オーダーの液滴では L=0.1〜1mm と非常に小さく、表面張力が支配します。インクジェット噴射、PCR の油中液滴、ラブオンチップの混合促進など、温度勾配で意図的にマランゴニ流を作って混合・輸送する設計が進んでいます。

よくある誤解と注意点

まず最大の落とし穴が、「微小重力＝完全な静止流体」と思い込むこと。ISS 内でも残留加速度（g-jitter、宇宙飛行士の運動、推進器噴射）は 10⁻⁶〜10⁻³ g 程度あり、低周波の擾乱として加わります。さらに本ツールで示すように、わずか 1K の温度差でもマランゴニ対流の Ma は数千〜数万に達し、流体は地上の自然対流と同等以上のスピードで循環します。「無重力 ＝ 何も動かない」は誤りで、「無重力 ＝ 表面張力と濃度差が主役になる」が正しい理解です。実験設計でも、温度均一性を 0.1K オーダーまで管理しないとマランゴニ流が乱す事例が多くあります。

次に、「表面張力 σ そのものが Ma を決める」という誤解。本ツールでは入力としてσ（mN/m）を表示用に用意していますが、Ma の駆動力は σ ではなく 温度勾配 ∂σ/∂T です。水の σ は 72 mN/m ありますが、∂σ/∂T はわずか 0.15 mN/m/K。一方、溶融金属は σ こそ 500 mN/m 程度ですが ∂σ/∂T は 0.36 mN/m/K と大きく、温度差があれば強い駆動力になります。教科書で Ma の式を見るとき、「σ ではなく dσ/dT」を必ず確認してください。本ツールの内部計算も dσ/dT（fluidProps.sigma_T）を用いています。

最後に、「Ma が大きい＝対流が良いこと」ではない。結晶成長やコーティングでは、Ma が臨界値（流体形状や Pr 数で変わるが概ね 500〜2000）を超えると対流が定常状態から「振動マランゴニ対流」に遷移し、表面波が立ち、結晶や薄膜に周期的な欠陥として刻まれます。ISS 実験の主目的の一つはこの遷移臨界の精密測定で、地上では浮力対流に隠れて測れない値です。本ツールで Ma > 500 になると warn / ok の verdict が切り替わるのは、この遷移目安を反映しています。

使い方ガイド

  1. 温度差(K)、表面張力勾配(mN/m/K)、動粘度(mPa·s)、特性長(mm)を入力します
  2. シミュレーターがマランゴニ数Maとレイリー数Raを計算し、微小重力環境(g=9.81×10⁻⁶ m/s²)での対流挙動を判定します
  3. Ma/Ra比率から支配的な対流モード(マランゴニ優位/ブイヤンシー優位/混合)を特定し、表面流速と動的Bond数を読み取ります

具体的な計算例

シリコーンオイル(μ=20 mPa·s、dσ/dT=-0.16 mN/m/K)の液滴、特性長3mmで温度差5K印加時：Ma=[(0.16×5×3)/(20×10⁻⁶)]≈1,200。ISS微小重力(g=1.0×10⁻⁶ m/s²)でRa≈0.08と極小化するため、マランゴニ対流が圧倒的に優位となり、赤道付近での表面流速は約8.5 mm/sに達します。動的Bond数Bo*=0.15以下で液滴は球形維持。

実務での注意点