ラバランプシミュレーター 戻る
熱解析

ラバランプ シミュレーター

熱膨張・アルキメデス浮力・粘性抵抗を使ったラバランプの物理シミュレーション。ブロブが上昇・下降するサイクルで自然対流の本質を体験しよう。

熱膨張 + 浮力 クリックで熱源追加 アルキメデスの原理
ラバランプ流動
パラメータ
加熱強度
ブロブ数
粘度(抵抗)
ブロブサイズ
px
速度倍率
x
カラースキーム
コントロール
統計
計算結果
0
高温ブロブ
0
低温ブロブ
0.5
平均温度
60
FPS
理論・主要公式

$$F_b = \rho_{fluid} g V$$

浮力(アルキメデスの原理):\(\rho_{fluid}\) 流体密度、\(g\) 重力加速度、\(V\) 排除体積

$$\rho(T) = \rho_0 [1 - \beta(T - T_0)]$$

熱膨張による密度変化:\(\beta\) 熱膨張係数、\(T_0\) 基準温度

$$F_{net} = (\rho_{fluid} - \rho_{wax}) g V$$

ワックスの浮沈は密度差で決定。\(F_{net}>0\) で上昇、\(<0\) で下降

ラバランプの物理とは

🙋
ラバランプって、あのゆらゆら動くインテリアですよね?あれがどうして動き続けるんですか?
🎓
大まかに言うと、温まると膨らんで軽くなり、冷えると縮んで重くなる、その繰り返しだよ。このシミュレーターで「加熱強度」のスライダーを上げてみて。ブロブが温まるスピードが上がって、もっと活発に動き出すのがわかるよ。
🙋
「軽くなる」って、油の中なのに浮くんですか?水に氷が浮くのと同じ原理ですか?
🎓
その通り!アルキメデスの原理、つまり浮力が働いているんだ。ブロブの密度が周りの流体より小さくなると、浮力で上に押し上げられる。逆に冷えて密度が大きくなると沈む。「熱膨張係数」を大きくすると、温度変化で密度が激しく変わるから、浮き沈みがよりドラマチックになるんだ。確認してみて。
🙋
でも、上まで行ったブロブがすぐに落ちてこないのはなぜ?なんかゆっくりしてます。
🎓
良いところに気づいたね。それは「粘性」のせいだ。油の中を動くのは水の中を動くよりずっと抵抗が大きいんだ。このシミュレーターの「流体粘性」パラメータをいじるとその影響がはっきり分かるよ。粘性を高くすると動きが鈍重になり、低くするとサッと動くようになる。現実のラバランプのゆったりした動きは、高い粘性と浮力のバランスで生まれているんだ。

よくある質問

ブロブの温度が十分に上がっていない可能性があります。ヒーター出力を上げるか、シミュレーション速度を遅くして加熱時間を確保してください。また、粘性抵抗が高すぎると上昇しにくくなるため、流体の粘度パラメータを下げてみてください。
はい、影響します。浮力はブロブの体積に比例するため、大きなブロブほど強い浮力を受けます。ただし、体積が大きいと粘性抵抗も増えるため、上昇速度は単純に比例しません。サイズと温度のバランスが動きの鍵です。
αを大きくすると、同じ温度変化でもブロブの密度が大きく減少し、浮力が強くなります。結果としてブロブが早く上昇し、サイクルが短くなります。逆にαを小さくすると変化が鈍くなり、ラバランプの動きが緩慢になります。
本シミュレーターは熱膨張・浮力・粘性抵抗という主要な物理原理をモデル化していますが、表面張力や化学的な界面効果などは省略しています。そのため定性的な動きの再現は可能ですが、実際のランプの複雑な挙動を完全に再現するものではありません。

実世界での応用

建築・環境工学(室内気流解析):ラバランプと同じ自然対流の原理は、室内の暖気や空調による気流のシミュレーションに応用されます。暖房器具の配置による温度ムラや、換気効率を予測するのに役立ちます。

電子機器の冷却設計:発熱する電子部品(CPUなど)の上に生じる上昇気流(熱プルーム)は、自然対流冷却の基本です。ヒートシンクのフィン間の空気流動を最適化する際の基礎物理となります。

地学・海洋学(マントル対流):地球のマントルは固体ですが、長い時間スケールで見ると粘性の高い流体のように振る舞い、内部の温度差によって対流しています。プレートテクトニクスを駆動する大規模な「地球のラバランプ」モデルです。

化学プロセス工学:反応槽内で温度差によって生じる自然対流は、物質の混合や熱伝達に影響を与えます。撹拌機を使わない静かな混合プロセスを設計する際の検討要素になります。

よくある誤解と注意点

このシミュレーターを使い始めるときに、いくつかつまずきやすいポイントがあるから気をつけてね。まず一つ目は、「加熱強度を最大にすれば一番激しく動く」と思いがちなこと。確かにブロブの加熱は速まるけど、上がりすぎるとブロブがすぐに天井まで到達してしまい、逆に下降して冷却されるまでの「ゆったりしたサイクル」が崩れてしまうんだ。現実のラバランプの落ち着いた動きを再現したいなら、加熱強度は中程度から始めて、他のパラメータとバランスを取るのがコツだよ。

二つ目は、「熱膨張係数」と「流体の粘性」の影響を混同してしまうこと。熱膨張係数(α)は、温度変化に対する密度の変化の激しさを決めるパラメータだ。例えばαを2倍にすると、同じ温度上昇で密度が2倍の勢いで下がるから、浮力が強く、急激に上昇するようになる。一方、粘性はその動きに対するブレーキの強さだ。αを大きくして浮力を強くしても、粘性も同時に大きくしすぎると、ブロブは重たい油の中を這うようにしか動かなくなる。パラメータをいじる時は、「どっちの効果を変えているのか」を意識してみよう。

最後に、これは実務でもありがちな落とし穴だけど、シミュレーションの「見た目」だけを追いかけないこと。例えばブロブの形が完全な球から少し崩れても、本質である「浮力と粘性抵抗のバランスによる対流」の原理は変わらない。数値シミュレーションでは、このような単純化されたモデルで現象の本質を抽出し、理解することが第一歩なんだ。細部のリアルさにこだわりすぎて、核心を見失わないようにね。

使い方ガイド

  1. 加熱強度スライダ(Heat)を0~100の範囲で調整し、ラバランプ内部の温度差を設定します。値が大きいほど底部ヒーターの熱出力が増加します
  2. 熱膨張係数スライダ(Thermal Expansion)で流体の密度変化率を0.0001~0.01/℃の範囲で設定し、温度1℃当たりの体積膨張を制御します
  3. 粘性抵抗スライダ(Drag)で流体の粘度係数を調整(0.1~10 Pa·s相当)し、ブロブの上昇・下降速度を変更します
  4. ブロブサイズスライダ(Size)で直径5~30mm範囲のパラメータを設定し、浮力と抗力のバランスに影響させます
  5. 再生ボタンを押してサイクルを開始し、温度成層化によるブロブの周期的な上下運動を観察します

具体的な計算例

標準的なシリコンオイル(ρ=960 kg/m³、α=0.0007/℃)を使用した場合、加熱強度60、膨張係数0.0005/℃に設定すると、底部が60℃上昇した時点でブロブ密度は0.968 kg/m³に低下します。直径15mmの球形ブロブ(体積1.77×10⁻⁶ m³)の浮力は約1.52×10⁻⁵ Nとなり、粘性抵抗3 Pa·s下ではおよそ8秒で上昇が完了。冷却開始から約45秒でサイクルが一周します

実務での注意点

  1. 実際のラバランプ(Lava Lamp)ではワックス密度が0.9~0.95 g/cm³で、液体石油が0.82 g/cm³のため、温度差10℃程度で浮力が反転します。シミュレータの膨張係数を0.0005~0.001/℃に設定して現物に近づけてください
  2. 粘性抵抗が高すぎるとブロブが停滞し、低すぎると急速に上下して観察困難になります。粘度2~5 Pa·s相当が安定した周期運動を生成します
  3. 容器底部の加熱均一性が重要で、実機では温度勾配により複数のブロブが異なるタイミングで上昇します。シミュレータの加熱強度を40~80の間で変動させて、このばらつきを再現できます