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

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

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

熱膨張 + 浮力 クリックで熱源追加 アルキメデスの原理
ラバランプ流動
パラメータ
加熱強度
ブロブ数
粘度(抵抗)
ブロブサイズ
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\) で下降

ラバランプの物理とは

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

よくある質問

ブロブの温度が十分に上がっていない可能性があります。ヒーター出力を上げるか、シミュレーション速度を遅くして加熱時間を確保してください。また、粘性抵抗が高すぎると上昇しにくくなるため、流体の粘度パラメータを下げてみてください。
はい、影響します。浮力はブロブの体積に比例するため、大きなブロブほど強い浮力を受けます。ただし、体積が大きいと粘性抵抗も増えるため、上昇速度は単純に比例しません。サイズと温度のバランスが動きの鍵です。
熱膨張係数 α は密度式 ρ = ρ₀[1−α(T−T₀)] の係数で、α が大きいほど同じ温度変化でも密度が大きく減少し浮力が強まります。本ツールでは α は内部で一定値(0.4)に固定されており調整スライダはありません。動きの活発さは「加熱強度」と「粘度(抵抗)」で調整してください。
本シミュレーターは熱膨張・浮力・粘性抵抗という主要な物理原理をモデル化していますが、表面張力や化学的な界面効果などは省略しています。そのため定性的な動きの再現は可能ですが、実際のランプの複雑な挙動を完全に再現するものではありません。

実世界での応用

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

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

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

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

よくある誤解と注意点

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

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

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

使い方ガイド

  1. 加熱強度スライダ(0.2〜3.0)で底部ヒーターの強さを設定します。値が大きいほどブロブが速く温まり活発に動きます
  2. ブロブ数(5〜25)・ブロブサイズ(10〜40px)でブロブの個数と大きさを調整します
  3. 粘度(抵抗)スライダ(0.005〜0.08、無次元の減衰係数)で動きの粘っこさを、速度倍率(0.5〜3.0)でシミュレーションの時間スケールを変えます
  4. カラースキームを選び、キャンバスをクリックすると局所熱源が追加されます。一時停止/リセットで再生を制御します
  5. ※ 熱膨張係数 α は内部で一定値に固定されており(参考・本計算では未使用の調整UIはありません)、密度変化 ρ = ρ₀[1−α(T−T₀)] は温度 T のみで決まります

具体的な計算例

本ツールは定性的な可視化モデルです。密度は ρ = 1 − 0.4·(T−0.5)(T は 0〜1 の無次元温度)で表され、浮力 ∝ (ρ_fluid − ρ) でブロブが上下します。例として加熱強度を 1.0 から 2.5 に上げると、底部のブロブが速く高温(T→1)になり密度が下がって上昇が活発になります。粘度(抵抗)を 0.08 側へ上げると上下動がゆっくりになり、現実のラバランプに近い落ち着いたサイクルになります。表示される数値は高温/低温ブロブ数・平均温度・FPS で、SI 単位の力や秒数は算出しません。

実務での注意点

  1. 実際のラバランプではワックス密度が0.9〜0.95 g/cm³、液体石油が0.82 g/cm³で、温度差10℃程度で浮力が反転します。本ツールの密度モデルは無次元の定性表現であり、SI密度には対応しません
  2. 粘度(抵抗)が高すぎるとブロブが停滞し、低すぎると急速に上下して観察困難になります。0.02〜0.04 付近が安定した周期運動を生みます
  3. 実機では底部加熱の温度勾配で複数ブロブが異なるタイミングで上昇します。本ツールでも加熱強度を中程度(1.0〜2.0)にし、クリック熱源を併用するとそのばらつきを観察できます