JGeek Log

虫記事 科学記事 SF About 物置 Amazon.co.jp
ディアスポラ数理研 /白熱光メモ 読み中

2014-10-10

[] 地球に穴をあけたらどうなる?

この質問。http://q.hatena.ne.jp/1412645084

たかだか標高8kmの山に登った程度で命に関わるくらい空気が薄くなるわけで、じゃあ反対に6000kmくらい下にいったら超はげしく空気が濃くなるよね、って話。

理想気体を仮定して計算すると指数関数的に濃くなっていき、数百kmで液体の濃度に近くなり理想気体として扱えなくなります。

ちなみに、圧力の増加=高度差×モル濃度×重力加速度×モル質量、これを解くとモル濃度=地表でのモル濃度×exp(深度×モル質量×重力加速度)

理想気体で計算できる範囲では、指数的増加のスピードはモル質量が重いほど速いので、地表では少ないキセノンなんかが途中で酸素や窒素を追い抜いて最初に液体の濃度に達します。

そうなるともう酸素や窒素は共存できなくなって途中からはほぼ100%キセノンのみという状態になるはずです。

ちょいと面倒ですが実際に計算してみました。

混合気体の状態方程式は、ファンデルワールス方程式を使います。理想気体からのズレを表す係数は各気体成分の値から分率を使って計算されます。

圧力の増加=高度差×モル濃度×重力加速度×モル質量、で圧力の増加を計算し、状態方程式からモル濃度の増加を計算。

次に各成分について自由エネルギー= 高度×モル濃度×重力加速度×モル質量- RT モル濃度 Log(モル濃度)、の変分を使って、これと全体の増分が上の計算と一致する拘束条件でラグランジュの未定乗数で各成分の増分を計算することを繰り返します。

100km付近でキセノンの一人勝ちになり他の気体が減っていきます。

f:id:ita:20141010221555p:image

割合はこちら。最初に重い酸素やアルゴンが増えていくけど、途中からキセノンに押され150kmからはほぼキセノンが100%に。

f:id:ita:20141010221554p:image

そこから下のXeの圧力変化はこんなかんじ。周囲の温度はどっかから適当な値を拾ってきて線形補完しました。右の相図はシミュレーションの論文から引っ張ってきました。圧力が上がって1000km付近で固化。地球中心付近では重力が0に近づくので圧力はサチるけど温度が上がっていくならまた液化するかも、というところ。ちょうど相境界のへんをうろうろしてるんでどっちになるか微妙ですが。

ちなみにXeの最大圧力は80GPa程度に対し、周囲の地球の圧力は300GPa以上。個体Xeよりはるかに重い鉄とかいろいろあるんでトンネルの壁がなければさらに潰されてしまいます。圧力が周囲と同じ、という条件でやればもっと早く固化するでしょう。

http://www.nature.com/nature/journal/v401/n6752/fig_tab/401432a0_F1.html#figure-title

f:id:ita:20141010221556p:image

Contact: Mitsuhiro Itakura/板倉充洋
ita.mitsu spam @ gmail, sausage, spam, egg, and com
最近のコメント

#. DATE NAME

1. 09/21 ita
2. 09/21 たけぞー
3. 09/14 ita
4. 09/14 syayuuzin
5. 08/31 ita
6. 08/31 星谷 仁
7. 07/15 ita
8. 07/15 客員研究員
9. 08/21 客員研究員
最近のTB

#. DATE  NAME

CALENDAR
<< 2014/10 >>
1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30 31