■Back to TOP-PAGE

 特殊相対論的に正しい「ツィオルコフスキーの公式」。
  ―― ロケットの運動について、特殊相対論を使って、いろいろと調べてみるよ♪

■『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』では、特に、「等加速度運動」する質点について、特殊相対論を使って、詳しく調べてみました。でも、もちろん、「等加速度運動」だけでなく、どんな加速度運動でも、特殊相対論を使って、詳しく調べるコトができます。たとえば、いわゆる「ロケット」の運動は、加速度運動ですけど、一般に、その加速度は時々刻々と変化します。したがって、「等加速度運動」ではありません。でも、特殊相対論を使って、その運動の様子を、詳しく調べるコトができます。

■ところで、改めて、そもそも、「ロケット」とは、ナンでしょうか。人それぞれ、思い浮かべるイメージは、微妙に違うかも知れませんね。ココでは、機体に搭載した「推進剤」を射出し、その反動(反作用)によって、推進力(推力)を得る機関、または、そのような機関を動力源にして運動する機体――と、定義します。推進剤は、あらかじめ機体に搭載されるモノだけではなく、あらかじめ機体に搭載されるモノだけから生成されるモノでも、構わないコトにします(その場合、誤解の恐れがないときには、推進剤の原料とするために、あらかじめ機体に搭載されるモノも、「推進剤の原料」とは言わず、ただ単に、「推進剤」と呼ぶコトが多いです)。


□たとえば、日本が誇る「LE-5B」エンジン(「H-IIA」ロケットの第 2 段エンジンなどに採用)は、搭載した液体水素と液体酸素を化学反応させて噴射し、その反動で推進力を発生します。液体水素と液体酸素の化学反応で生成される物質($\mathrm{H}_{2}\mathrm{O}$=水)が「推進剤」です(ちなみに、厳密な説明をするのは面倒なので、誤解の恐れがないときは、「LE-5B エンジンの推進剤は、液体水素と液体酸素です」みたいな表現をするコトが多いです)。この推進剤を、液体水素と液体酸素の化学反応で生じたエネルギーを使って噴射して、推進力を得ているワケです。その意味で、LE-5B エンジンは、「ロケット」の定義を満たします。また、「イプシロンロケット」の第 1 段モーター(SRB-A)は、搭載した固体燃料を燃焼(化学反応)させて、生じた燃焼ガスを噴射し、その反動で推進力を発生します。したがって、SRB-A も、「ロケット」の定義を満たしています。どちらも、化学反応で発生したエネルギーを使って、推進剤を噴射しているのが、特徴です。そのような特徴を持つロケットを、「化学ロケット」と呼ぶコトがあります。

□他方で、小惑星探査機「はやぶさ 2」のイオンエンジン $\mu 10$ は、搭載しているキセノンガスをイオン化して、機体の外に噴射し、その反動で推進力を発生します。その意味で、コレも「ロケット」と言えます。キセノンのイオン化と、噴射に必要なエネルギーは、太陽電池パネルから得た電気エネルギーを使います。そんなワケで、イオンエンジンは、化学ロケットと区別して、「電気推進ロケット」と呼ばれています。

□なお、ロケットに良く似たモノとして、「ジェットエンジン」があります。ジェットエンジンは、機体の外から空気を取り込んで、機体に搭載した燃料を燃焼させます。そして、その燃焼ガスを噴射する反動で、推進力を発生します。すなわち、噴射する燃焼ガス(=推進剤)は、あらかじめ機体に搭載される燃料と、機体の外から取り込んだ空気の、化学反応で生成されるワケです。推進剤の原料に、あらかじめ機体に搭載されていないモノ(機体の外から取り込んだ空気)が含まれるので、ジェットエンジンは、「ロケット」の定義を満たしません。実は、このコトが、ロケットの運動を計算するときに、1 つのポイントとなります。


■さて、ロケットの運動(具体的には、ロケットの位置や速度や加速度などが、どのように変化するのか)は、ニュートン力学の世界でも、特殊相対論の世界でも、理論的に、厳密に計算するコトができます。計算結果は、「宇宙旅行の父」として知られる コンスタンチン・ツィオルコフスキー にちなんで、「ツィオルコフスキーの公式」と呼ばれています。ツィオルコフスキー は、ライト 兄弟が人類初の有人動力飛行に成功するよりも前の 1897 年に、ロケットの速度変化を支配する公式を発表しました。そして、すでに良く知られ、利用可能だったエネルギー源と材料だけを使うという制約下でも、宇宙旅行が、原理的に可能であるコトを、示したのです。


□藤野は、学生時代に、「軌道エレベーター」について、詳しく計算してみたコトがあります。その結果、軌道エレベーターを建造するのに十分な強度と長さを持つ材料が、この世に存在しない――つまり、未来に開発されるかも知れない夢の新素材が登場するまで、軌道エレベーターを建造するコトは、原理的に不可能である――という結論が導かれて、ガッカリしたコトを覚えています。既存のエネルギー源と材料だけで、宇宙旅行が可能だと知った、ツィオルコフスキー の喜びと興奮は、どれほどのモノだったでしょうね。


■ツィオルコフスキー は、ニュートン力学に基づいて、ロケットの運動を計算しました。では、特殊相対論に基づいて計算した場合は、どうなるのでしょうか。ココで、特殊相対論に基づいて、ロケットの運動を計算し、「特殊相対論的に正しいツィオルコフスキーの公式」を導いて、その結果を、詳しく検討してみましょう。


【印刷用 PDF「特殊相対論的に正しい「ツィオルコフスキーの公式」。」】(約 1.66 MB)

【目 次】


1. 必要な準備(1)
  ―― 物理現象を表す微分方程式を得るための関係式

2. ニュートン力学的に正しいツィオルコフスキーの公式
  ―― 基本事項の復習

3. 必要な準備(2)
  ―― 特殊相対論における基本的な関係式

4. 特殊相対論的に正しいツィオルコフスキーの公式(1)
  ―― 化学ロケット編

5. 特殊相対論的に正しいツィオルコフスキーの公式(2)
  ―― 電気推進ロケット編

6. ツィオルコフスキーの公式を解く
  ―― ロケットの運動と時間の関係


■このページの一番上に戻る。



 ■1. 必要な準備(1) ―― 物理現象を表す微分方程式を得るための関係式

■まず最初に、ニュートン力学の世界でも、特殊相対論の世界でも、ツィオルコフスキーの公式を導くときに、絶大な威力を発揮する、とても重要な関係式を紹介しましょう。ツィオルコフスキーの公式の導出に限ったコトではありませんけど、物理現象の様子を詳しく調べるのに、最も良い方法は、その物理現象を「微分方程式」で表す――そして、その微分方程式を、詳しく調べる――という方法です。よって、ココで紹介する、「物理現象を表す微分方程式を得るための関係式」は、あらゆる物理現象の計算において、最も重要な関係式である――と言うコトができるでしょう。この関係式を用いた、ツィオルコフスキーの公式の導出は、ほぼ、どんな物理現象の計算にも、まったく同じやり方で、応用するコトができます。ぜひ、その方法を、覚えていって下さいね。

■さて、微分(導関数)の定義式を、思い出して下さい。ある変数 $x$ の値に依存して、値が定まる関数 $f(x)$ を考えます。この関数 $f(x)$ の微分(導関数)$df(x)/dx$ は、(高校数学までは)次のように定義されます。

\[ \frac{df(x)}{dx} = \lim_{h \, \to \, 0}{\frac{f(x + h) - f(x)}{h}} \qquad \mbox{または} \qquad \frac{df(x)}{dx} = \frac{f(x + h) - f(x)}{h} \quad \mbox{($h \, \to \, 0$)} \tag{1.1} \]

この定義式は、どんな教科書にも載っていますので、お馴染みのモノだと思います。ところで、この(1.1)式は、次のように変形するコトができます。

\[ f(x + h) = f(x) + \frac{df(x)}{dx} \cdot h \quad \mbox{($h \, \to \, 0$)} \tag{1.2} \]

この(1.2)式こそが、ココで紹介したかった、「物理現象を表す微分方程式を得るための関係式」です。次の第 2 節から、いよいよ、この(1.2)式を駆使して、実際に、ロケットの運動の様子を調べていきます。


□ちなみに、「物理現象を表す微分方程式を得るための関係式」は、藤野の造語ですので、ココ以外では、たぶん、通じません。一般的には、「($h$ に関する)$1$ 次の近似式」と言うコトが、多いようです。近似式ではないのに、なぜ、「近似式」と呼ぶのでしょうか。(1.2)式は、$h$ がゼロ(または、ゼロに限りなく近い値)のときにだけ、成り立ちます。でも、$h$ の値が、ゼロからほんの少しズレただけなら、(1.2)式の両辺の値も、少ししかズレない――つまり、(1.2)式が近似的に成り立つと、みなすコトができます。そこで、$h$ の値が、ゼロの場合だけでなく、ゼロからほんの少しズレただけの場合(ゼロに十分に近い値の場合)も含めて考えるコトが、便宜上、とても多いのです。そのために、「近似式」と呼ぶワケです。

□なお、「近似式」としては、(1.2)式よりも近似の精度を高めた、次の式のほうが有名です。
\[ \qquad f(x + h) \simeq f(x) + \frac{df(x)}{dx} \cdot h + \frac{1}{2!} \frac{d^{2}f(x)}{dx^{2}} \cdot h^{2} + \frac{1}{3!} \frac{d^{3}f(x)}{dx^{3}} \cdot h^{3} + \cdots \quad \mbox{($\lvert h \rvert \ll 1$)} \tag{1.3} \]
(1.3)式を、「テイラー展開(テイラー級数展開)」と呼びます。ココで、「$\lvert h \rvert \ll 1$」は、「$h$ の絶対値は、$1$ よりも十分に小さい」という意味です。コレは、「$h$ は、ゼロか、ゼロに十分に近い値である」と、ほぼ同じ意味です。また、$h$ の値が正であるコトが明らかな場合には、「$h \ll 1$」と書くコトも多いです。ところで、テイラー展開は、物理現象の様子を調べるのに、とても有用なのですけど、それ自体も、数学的に、非常に興味深いモノです。ぜひ、いろいろな側面から、テイラー展開のコトを、調べてみて下さいね。

■このページの一番上に戻る。



 ■2. ニュートン力学的に正しいツィオルコフスキーの公式 ―― 基本事項の復習

■ココでは、「特殊相対論的に正しいツィオルコフスキーの公式」を導く前に、まずは、「ニュートン力学的に正しいツィオルコフスキーの公式」の導き方を、復習しておきましょう。次の図のように、時刻 $t$ のときの、推進剤を含んだロケット全体の質量を $m(t)$、ロケット(rocket)の速さを $u_{\mathrm{R}}(t)$ と表すコトにします。時間が $\varDelta t$ だけ経過すれば、質量と速さは、当然、それぞれ、$m(t+\varDelta t)$ と $u_{\mathrm{R}}(t+\varDelta t)$ になります。ココで、ロケットは、時間 $\varDelta t$ のあいだに、質量 $\varDelta m_{\mathrm{P}}$ の推進剤(propellant)を、ロケットから見て速さ $-\omega_{\mathrm{P}}$ で噴射します。このとき、静止している観測者から見たら、推進剤の噴射速度は $u_{\mathrm{R}}(t)-\omega_{\mathrm{P}}$ になるコトに、注意して下さい。また、ロケットには、ロケット自身の推力(thrust)$F_{\mathrm{T}}$ のほかに、外力 $f$ がはたらくコトも、忘れずに、想定しておきましょう。たとえば、重力に逆らって加速する場合、外力 $f$ として、$f=-m(t)g$($\, g$ は重力加速度)を、計算に入れておく必要があります。どの要素を考えるときも、向きには、特に気をつけて下さいね。
図1:「ニュートン力学的に正しいツィオルコフスキーの公式」の導出(約45KB)

 ■「ツィオルコフスキーの公式」を導く ――  

■以上の設定で、ロケットの運動の様子を、実際に計算してみましょう。ただし、本当に知りたいのは、時刻 $t$ から時刻 $t+\varDelta t$ のあいだ――といった曖昧な状況のハナシではありません。本当に知りたいのは、時刻 $t$ の、その瞬間における、ロケットの運動です。それには、経過時間 $\varDelta t$ が、ゼロ(または、ゼロに限りなく近い値)でなければいけません。このとき、絶大な威力を発揮するのが、第 1 節で紹介した(1.2)式です。(1.2)式より、時刻 $t+\varDelta t$ での、推進剤を含んだロケット全体の質量 $m(t+\varDelta t)$ と、ロケットの速さ $u_{\mathrm{R}}(t+\varDelta t)$ は、

\[ m(t + \varDelta t) = m(t) + \frac{dm(t)}{dt} \cdot \varDelta t \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{2.1a} \] \[ u_{\mathrm{R}}(t + \varDelta t) = u_{\mathrm{R}}(t) + \frac{du_{\mathrm{R}}(t)}{dt} \cdot \varDelta t \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{2.1b} \]

と表せるコトが分かります。また、ロケットの場合は、推進剤が、あらかじめ、ロケットに搭載されています。したがって、噴射した推進剤の質量 $\varDelta m_{\mathrm{P}}$ と同じだけ、ロケット全体の質量は軽くなるハズです。式で表すと、

\[ m(t + \varDelta t) = m(t) - \varDelta m_{\mathrm{P}} \tag{2.2} \]

となります(ちなみに、コレは、$\varDelta t$ の値にかかわらず、ニュートン力学の世界では、必ず成り立っているハズなので、$\varDelta t \, \to \, 0$ の条件は不要です)。この(2.2)式に、(2.1a)式を代入すると、ただちに、

\[ \varDelta m_{\mathrm{P}} = -\frac{dm(t)}{dt} \, \varDelta t \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{2.3} \]

となるコトが分かります。マイナス記号が、少し不自然に感じるかも知れません。でも、推進剤を噴射すると、ロケット全体の質量 $m(t)$ は減少します。つまり、$dm(t)/dt$ は、必ず、負の値になります。よって、(2.3)式は、見た目に反して、噴射した推進剤の質量 $\varDelta m_{\mathrm{P}}$ が、必ず、正の値になるコトを、表しているのです。

■必要な準備が揃ったので、いよいよ、いま考えている状況に、「運動量の保存則」を、適用してみましょう。時間 $\varDelta t$ のあいだ、外力 $f$ がはたらく場合は、力積 $f \cdot \varDelta t$ だけ、全体の運動量が増加するコトに注意すると、この場合の運動量の保存則は、($\varDelta t$ の値にかかわらず)次のように表されます。

\[ \bigl[ m(t) \cdot u_{\mathrm{R}}(t) \bigr] + f \cdot \varDelta t = \bigl[ m(t + \varDelta t) \cdot u_{\mathrm{R}}(t + \varDelta t) \bigr] + \bigl[ \varDelta m_{\mathrm{P}} \cdot \bigl( u_{\mathrm{R}}(t) - \omega_{\mathrm{P}} \bigr) \bigr] \tag{2.4} \]

この(2.4)式に、(2.1a)〜(2.1b)式、および、(2.3)式を代入して、両辺を $\varDelta t$ で割ると、ただちに、

\[ f = m(t) \, \frac{du_{\mathrm{R}}(t)}{dt} + \frac{dm(t)}{dt} \frac{du_{\mathrm{R}}(t)}{dt} \, \varDelta t + \omega_{\mathrm{P}} \, \frac{dm(t)}{dt} \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{2.5} \]

となります。ところで、時刻 $t$ の、その瞬間におけるロケットの運動を知るには、$\varDelta t$ はゼロ(または、ゼロに限りなく近い値)でなければいけませんでした。そこで、(2.5)式に、$\varDelta t = 0$ を代入して、整理すると、

\[ f - \omega_{\mathrm{P}} \, \frac{dm(t)}{dt} = m(t) \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{2.6} \]

という微分方程式が得られます。コレこそが、「ニュートン力学的に正しいツィオルコフスキーの公式」です。この微分方程式を解けば、どんなロケットの運動の様子でも、すべて分かる――というワケです。


□途中までは、「$\varDelta t \, \to \, 0$ ならば、結局、$m(t) = m(t)$ という、当たり前の等式になったり、両辺がともにゼロになるという、つまらない等式になったりするだけなのに、なんで、こんな計算を続けるんだろう?」と、フシギに思ったかも知れませんね。実は、多くの物理現象で、同じような計算をすると、このように、時刻 $t$ の瞬間の状況だけを抜き出した微分方程式が得られるのです。コレこそが、(1.2)式の真価です。

□(2.6)式を、ニュートンの運動方程式 $F = ma$ と見比べてみると、(2.6)式の右辺は、ニュートンの運動方程式の右辺 $ma$(質量×加速度)と、同じであるコトが分かります。よって、(2.6)式の左辺は、ロケットにはたらく力の合計 $F$ を表しているコトになります。実際に、(2.6)式の左辺第 1 項の $f$ は、ロケットにはたらく外力を表していたのでした。そうすると、(2.6)式の左辺第 2 項は、ロケット自身の推力 $F_{\mathrm{T}}$ を表しているに違いありません。すなわち、ロケットの推力 $F_{\mathrm{T}}$ は、
\[ \qquad F_{\mathrm{T}} = -\omega_{\mathrm{P}} \, \frac{dm(t)}{dt} \quad = \mbox{(推進剤の噴射速度)×(推進剤の消費速度)} \tag{2.7} \]
という関係式で与えられる――というコトになります。つまり、もし、より推力の大きいロケットエンジンを作りたいと思ったら、推進剤の噴射速度をより速くするか、推進剤の消費速度(単位時間当たりの推進剤の消費量)をより増やすか、そのどちらか(または、その両方)しかない――というコトが、分かります。


■「ツィオルコフスキーの公式=(2.6)式」を解くコトで、ロケットの運動の様子が、具体的に導かれます。そこで、以下では、いろいろな条件で、実際に(2.6)式を解いてみて、ロケットがどのように運動するのか、詳しく検討してみましょう。

 ■「ツィオルコフスキーの公式」を解く(1):外力がない場合 ――  

■さて、実際に、ツィオルコフスキーの公式=(2.6)式を解いてみる前に、簡単のため、以下では、すべての計算で、「推進剤の消費速度(単位時間当たりの推進剤の消費量)は、一定とする」と仮定します。ココでは、推進剤の消費速度を、$\mu$(定数)と書くコトにしましょう。分かりやすいように、推進剤の消費速度 $\mu$ の値は、必ず正($\, \mu > 0$)であるとします。つまり、$\mu$ の値が大きいほど、単位時間当たりに、たくさんの推進剤を消費するというコトになります。また、時刻 $t = 0$ のときの、推進剤を含んだロケット全体の質量 $m(t)$ を、$m_{0}$ で表すコトにします。まとめると、コレらの記号のあいだには、

\[ \mu \equiv -\frac{dm(t)}{dt} \quad \mbox{($\, \mu$ は定数)} \tag{2.8a} \] \[ m(t) = m_{0} - \mu t \tag{2.8b} \] \[ F_{\mathrm{T}} = \omega_{\mathrm{P}} \mu \tag{2.8c} \]

という関係が成り立っています。最後の(2.8c)式は、先に示した(2.7)式を書き換えたモノで、ロケットの推力 $F_{\mathrm{T}}$ が、推進剤の噴射速度 $\omega_{\mathrm{P}}$ と、推進剤の消費速度 $\mu$ の積で表されるコトを、示しています。


□もちろん、「推進剤の消費速度は、一定とする」という仮定をしなくても、ツィオルコフスキーの公式=(2.6)式を、厳密に解くコトはできます。むしろ、多くの教科書や解説では、このような仮定をせずに、ツィオルコフスキーの公式を解いています。でも、いろいろな設定で解いて、結果を比較する場合に、この仮定は、とても便利なので、ココでは、すべての計算で、この仮定を採用するコトにします。


■さらに、簡単のため、以下のすべての計算で、「推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、一定とする」という仮定も、追加しておきます(コレは、多くの教科書や解説でも、採用されている仮定です)。つまり、ココでは、$\omega_{\mathrm{P}}$ と $\mu$ の積で表される、ロケットの推力 $F_{\mathrm{T}}$ も一定である――というコトになります。


□ちなみに、ロケットの推力 $F_{\mathrm{T}}$ が、一定ではないような場合も、いろいろと、考えられると思いますけど(たとえば、推力が、時間の経過とともに、徐々に大きくなっていく場合とか、ロケットの加速度が、常に一定になるように、推力を調整する場合とか)、ぜひ、いろいろな状況や条件を、自分で想像・設定して、その場合に、ロケットの運動はどうなるのか、計算してみて下さいね。


■それでは、いよいよ、ツィオルコフスキーの公式=(2.6)式を、実際に解いてみましょう。最初は、外力がない場合、つまり、$f = 0$ の場合を考えます。このとき、(2.6)式は、(2.8a)〜(2.8b)式も用いて、

\[ \omega_{\mathrm{P}} \mu = ( m_{0} - \mu t ) \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{2.9} \]

と表されます。この微分方程式は、簡単に解くコトができます。初期条件として、時刻 $t = 0$ のときのロケットの速さ(初速度)$u_{\mathrm{R0}}$ を、$u_{\mathrm{R0}} = 0$ とすると(いろいろな解法や表記法がありますけど、一例として)、

\[ \int_{0}^{u_{\mathrm{R}}(t)}{du_{\mathrm{R}}} = \int_{0}^{t}{\frac{\omega_{\mathrm{P}} \mu}{m_{0} - \mu t} \, dt} \] \[ \quad \mbox{したがって} \quad u_{\mathrm{R}}(t) = -\omega_{\mathrm{P}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \, t \right)} \quad = \omega_{\mathrm{P}} \, \ln{\left( \frac{m_{0}}{m(t)} \right)} \tag{2.10} \]

となります(ほとんどの教科書や解説では、(2.6)式ではなく、(2.10)式のほうを、「ツィオルコフスキーの公式」と呼んでいるようです。ただし、コチラは、あくまでも、「外力がない場合の解」であるコトに、注意して下さい)。(2.10)式で与えられた、時刻 $t$ とロケットの速さ $u_{\mathrm{R}}(t)$ の関係を、次の図に示しています。はじめのうちは、ほぼ、等加速度運動($u_{\mathrm{R}}(t) \simeq (\omega_{\mathrm{P}}\mu/m_{0}) \, t \,$)をしていますけど、ロケット全体の質量 $m(t)$ がゼロに近づくにつれて、ロケットの速さ $u_{\mathrm{R}}(t)$ が、急激に、無限大へと向かって、発散していく様子を、見るコトができるでしょう(もちろん、ロケット全体の質量 $m(t)$ がゼロになるコトは、ありえませんけど……)。
図2:時刻とロケットの速さの関係:ニュートン力学・外力なし(約7KB)


□余談ですけど、実は、外力がない場合は、推進剤の消費速度 $\mu$ が一定でないとしても、推進剤を含んだロケット全体の質量 $m$ が同じならば、そのときのロケットの速さ $u_{\mathrm{R}}(m)$ は、必ず、同じ値になります。


■ところで、最初に搭載していた推進剤を使い切ったら、ロケットは、それ以上、加速できません。外力がない場合、推進剤を使い切ったロケットは、慣性運動=等速度運動をするコトになります。いま、消費する推進剤の総質量を $m_{\mathrm{P} \,}$、推進剤を除いたロケット本体の質量を $m_{\mathrm{R}}$ とします(加速前の、推進剤を含んだロケット全体の質量 $m_{0}$ は、$m_{0} = m_{\mathrm{P}}+m_{\mathrm{R}}$ と表されます)。これらの記号を使うと、推進剤を使い切るときの時刻 $t_{\mathrm{P}}$ が、

\[ t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \tag{2.11} \]

となるコトは、すぐに分かります。つまり、ロケットは、推進剤を消費しながら加速して、時刻 $t_{\mathrm{P}}$ のときに、最大速度 $u_{\mathrm{R \cdot MAX}} \:$ に到達する――というワケです。(2.10)式より、外力がない場合は、時刻 $t \ge t_{\mathrm{P}} = m_{\mathrm{P}}/\mu$ におけるロケットの速さ(最大到達速度)$u_{\mathrm{R \cdot MAX}} \:$ が、

\[ u_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \, \ln{\left( \frac{m_{0}}{m_{\mathrm{R}}} \right)} \quad = \omega_{\mathrm{P}} \, \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} \tag{2.12} \]

となるコトが分かります。(2.12)式より、外力がない場合は、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ が、推進剤の噴射速度 $\omega_{\mathrm{P}}$ と、$R \equiv m_{0}/m_{\mathrm{R}}$(この値は、一般に、「質量比」と呼ばれます)だけで決まる――という結論が、得られました。


□推力の大きさ $F_{\mathrm{T}}$ を決める重要な要素である、推進剤の消費速度 $\mu$ が、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に影響しない――という結果は、ちょっとした驚きかも知れませんね。もちろん、推進剤の消費速度 $\mu$ を大きくすると、推力 $F_{\mathrm{T}}$ が大きくなり、ロケットの加速度は大きくなります。でも、その分、早く、推進剤を使い果たしてしまうので、加速時間が短くなります。2 つの効果が、ちょうど相殺される結果、推進剤の消費速度 $\mu$ が変化しても、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は変化しない――というコトになるのです。もちろん、「とにかく、少しでも早く、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に到達したい」という場合には、推進剤の消費速度 $\mu$ を大きくしなければいけません。目的に応じて、手段を考えるコトが、ナニよりも重要です。

□さて、ココで、実際に、ロケットで宇宙旅行をするコトを、考えてみましょう。ロケットが、地表に落下せず、宇宙に留まるための方法の 1 つは、地球を周回する人工衛星になるコトです。そのような人工衛星になるためには、ロケットの速さが、目的の地球周回軌道の「軌道速度(orbital speed)」に、到達する必要があります。一般に、考える軌道が円軌道の場合、軌道速度 $v_{\mathrm{O}}$ は、重力と遠心力の釣り合いから、
\[ \qquad v_{\mathrm{O}} = \sqrt{\frac{GM}{R_{\mathrm{O}}}} \quad \mbox{($G$:万有引力定数、$M$:考える天体の質量、$R_{\mathrm{O}}$:軌道半径)} \tag{2.13} \]
で与えられます。地球周回軌道の場合には、軌道半径 $R_{\mathrm{O}}$ は、必ず、地球の半径より大きいハズなので、$v_{\mathrm{O}} \le 7.9 \times 10^{3}$ [$\mathrm{m/s}$] となります。$R_{\mathrm{O}}$ が大きくなるほど、$v_{\mathrm{O}}$ は小さくなるワケですけど、地表から、より遠くの軌道へ行く必要が生じます。そこで、なるべく、地表に近い軌道を目指すコトにして、さしあたり、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ として、$7.9 \times 10^{3}$ [$\mathrm{m/s}$] という値を、設定してみます。

□次に、推進剤の噴射速度 $\omega_{\mathrm{P}}$ について、考えてみましょう。ココでは、化学反応で生じた生成物を推進剤として使い、その化学反応で発生するエネルギーによって推進剤を噴射する――いわゆる「化学ロケット」について、ごく簡単に検討してみます。19 世紀当時から現在に至るまで、比較的、利用しやすい化学反応のうち、$1$ [$\mathrm{kg}$] の生成物が生じるときに発生するエネルギーが、最も大きいのは、水素と酸素から水蒸気を得る化学反応(水素の燃焼反応)です。水素の燃焼反応では、生成される水蒸気 $1$ [$\mathrm{kg}$] 当たり、およそ $1.34 \times 10^{7}$ [$\mathrm{J/kg}$] のエネルギー=標準生成熱が発生します。コレが、すべて、運動エネルギー $(1/2)\omega_{\mathrm{P}}^{2}$ に変換されると、推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、およそ $5.18 \times 10^{3}$ [$\mathrm{m/s}$] になります。もちろん、この速さは、原理的に可能な最大速度であって、現実の推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、コレよりも、必ず、小さくなります。実際には、標準生成熱の多くが、温度上昇などに使われて、推進剤の運動エネルギーへと変換できるのは、半分以下だからです。標準生成熱のうち、どの程度を、推進剤の運動エネルギーとして使えるかは、理論的にも、技術的にも、とても難しい問題ですけれど、おおざっぱに、少なめに見積もって、$20 \,$〜$\, 30 \, \%$ とした場合、推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、$3 \times 10^{3}$ [$\mathrm{m/s}$] 前後になるコトが、期待できます(この計算についてや、別の化学反応、化学反応以外のエネルギーを使う場合のコトは、ぜひ、自分で調べてみて下さいね)。

□最後に、消費する推進剤の総質量 $m_{\mathrm{P}}$ を、考えましょう。(2.12)式より、$m_{\mathrm{P}}$ は、ただちに、
\[ \qquad m_{\mathrm{P}} = \left( e^{u_{\mathrm{R \cdot MAX}} \ \ \ / \omega_{\mathrm{P}}} - 1 \right) m_{\mathrm{R}} \tag{2.14} \]
で与えられるコトが分かります。ココまでに想定してきた数値を代入すると、地球周回軌道への投入には、推進剤を除いたロケット本体の質量 $m_{\mathrm{R}}$ の、約 $13$ 倍の推進剤が必要である――という結果になりました。

□ところで、いまは、宇宙旅行の実現可能性を考えているので、「推進剤を除いたロケット本体」として、「人間が乗っている宇宙船」を想定する必要があります。目安として、軽トラックでは、乗員・貨物などを満載したときの「車両総重量」が、$1.2 \times 10^{3}$ [$\mathrm{kg}$] 弱のモノが多いみたいです。そこで、軽トラックほどの規模の「宇宙船」で、地球周回軌道へ「宇宙旅行」するには、$1.6 \times 10^{4}$ [$\mathrm{kg}$]($16$ [$\mathrm{t}$])程度の推進剤が必要――という結論になります。コレは、とんでもなく非常識なコトに思えます。そもそも、$16$ [$\mathrm{t}$] もの推進剤を収容するタンクだけで、ロケット本体の質量の大部分を占めてしまいそうです(容量 $200$ [リットル] のドラム缶の質量が、$17$ [$\mathrm{kg}$] 以上あるのです)。ただし、最も大切なのは、「宇宙旅行」は、間違いなく、とてつもない技術的な挑戦にはなるけれども、決して、物理的・原理的には不可能でない――という結果が得られたコトです。宇宙旅行が、夢想ではなく、目指すコトのできる目標になった瞬間です。


■さて、(2.12)式より、外力がない場合、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、推進剤の噴射速度 $\omega_{\mathrm{P}}$ と、質量比 $R = 1+m_{\mathrm{P}}/m_{\mathrm{R}}$($m_{\mathrm{P}}/m_{\mathrm{R}}$ は、ロケット本体の質量 $m_{\mathrm{R}}$ に対する、消費する推進剤の総質量 $m_{\mathrm{P}}$ の比)だけによって決まる――というコトが分かりました。では、もし、あなたが、ロケット工学者ならば、「推進剤の噴射速度 $\omega_{\mathrm{P}}$ を高速化する」、「消費する推進剤の総質量(=搭載可能な推進剤の質量)$m_{\mathrm{P}}$ を増量する」、「ロケット本体の質量 $m_{\mathrm{R}}$ を軽量化する」――のどれを、最も重視して、取り組むべきでしょうか。

■ココで、ちょっとした、数学的テクニックを、紹介しましょう。たとえば、いま、$y = x^{2}$ という関係式を、考えます。この関係式を、$x$ で微分すると、

\[ \frac{dy}{dx} = 2x \qquad \mbox{または} \qquad dy = 2x \, dx \tag{2.15} \]

という新たな関係式が得られます。コレは、「$x$ のを $dx$ だけ変化させると、$y$ のは $2x \, dx$ だけ変化する」というコトを、表しています。この(2.15)式の両辺を $y$ で割ると、次のようなカタチに、変形できます。

\[ \frac{dy}{y} = 2 \, \frac{dx}{x} \tag{2.16} \]

この(2.16)式は、「$x$ の割合を $dx/x$ だけ変化させると、$y$ の割合は $2 \, dx/x$ だけ変化する」というコトを、表しています。ロケット工学者の立場で、「$\omega_{\mathrm{P}}$ を $1$ [$\mathrm{m/s}$] 高速化するのと、$m_{\mathrm{R}}$ を $1$ [$\mathrm{kg}$] 軽量化するのでは、どちらが、$u_{\mathrm{R \cdot MAX}} \:$ の増加量を、より大きくするか」を知りたいときには、(2.15)式のカタチの関係式が便利です。けれど、「$\omega_{\mathrm{P}}$ を $1 \, \%$ 高速化するのと、$m_{\mathrm{R}}$ を $1 \, \%$ 軽量化するのでは、どちらが、$u_{\mathrm{R \cdot MAX}} \:$ の増加率を、より大きくするか」を知りたいときには、(2.16)式のカタチの関係式のほうが便利です。実際の現場では、単位を揃えて(無次元化して)比較できる、(2.16)式のカタチの関係式のほうが、役に立つでしょう。

■そんなワケで、$u_{\mathrm{R \cdot MAX}} \:$ について、(2.16)式のカタチの関係式を、求めてみましょう。(2.12)式より、ただちに、次の 3 つの関係式が得られます。

\[ \frac{du_{\mathrm{R \cdot MAX}}}{u_{\mathrm{R \cdot MAX}}} = \frac{d\omega_{\mathrm{P}}}{\omega_{\mathrm{P}}} \tag{2.17a} \] \[ \frac{du_{\mathrm{R \cdot MAX}}}{u_{\mathrm{R \cdot MAX}}} = \left\{ \frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{( 1 + m_{\mathrm{P}}/m_{\mathrm{R}} ) \, \ln{( 1 + m_{\mathrm{P}}/m_{\mathrm{R}} )}} \right\} \frac{dm_{\mathrm{P}}}{m_{\mathrm{P}}} \tag{2.17b} \] \[ \frac{du_{\mathrm{R \cdot MAX}}}{u_{\mathrm{R \cdot MAX}}} = \left\{ -\frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{( 1 + m_{\mathrm{P}}/m_{\mathrm{R}} ) \, \ln{( 1 + m_{\mathrm{P}}/m_{\mathrm{R}} )}} \right\} \frac{dm_{\mathrm{R}}}{m_{\mathrm{R}}} \tag{2.17c} \]

(2.17a)式から、推進剤の噴射速度 $\omega_{\mathrm{P}}$ を $1 \, \%$ 高速化させると、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は同じ割合($1 \, \%$)だけ速くなる――というコトが、分かります。それに比べて、(2.17b)〜(2.17c)式より、$dm_{\mathrm{P}}/m_{\mathrm{P}}$ と $dm_{\mathrm{R}}/m_{\mathrm{R}}$ の係数の絶対値は、必ず、$1$ より小さくなります(証明は省略しますけど、定義より、$m_{\mathrm{P}}/m_{\mathrm{R}} > 0$ であるコトに注意)。すなわち、搭載可能な推進剤の質量 $m_{\mathrm{P}}$ を $1 \, \%$ 増量した場合、または、ロケット本体の質量 $m_{\mathrm{R}}$ を $1 \, \%$ 軽量化した場合は、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は同じ割合($1 \, \%$)ほどには速くならない――というコトになります。したがって、もし、あなたがロケット工学者なら、まずは、ナニよりも、推進剤の噴射速度 $\omega_{\mathrm{P}}$ の高速化に取り組むべき――という結論が得られました。


□繰り返し、強調しなければいけませんけど、以上の結論は、あくまでも、「外力がない場合」という仮定のもとで、得られたモノです。たとえば、小惑星探査機「はやぶさ 2」などは、この仮定に当てはまると、考えても良いでしょう。したがって、あなたが、「はやぶさ 2」のイオンエンジン開発者ならば、推進剤の消費速度を考慮する必要はなく、さらに、機体の軽量化や、搭載する推進剤の増量よりも、優先すべきは、推進剤の噴射速度の高速化である――というコトになります。とは言え、もちろん、機体の軽量化や、搭載する推進剤の増量も、無視はできません。推進剤の噴射速度の高速化のために、機体の重量化や、搭載する推進剤の減量を強いるコトになる場合には、慎重な検討が必要になります(ちなみに、推進剤の消費速度の低下を強いるだけならば、気にする必要はありません)。

□しつこいようですけど、以上の結論は、あくまでも、「外力がない場合」という仮定のもとで、得られたモノです。外力がある場合には、異なる結論が得られるかも知れません。次に、重力がある場合について、ツィオルコフスキーの公式=(2.6)式を解いてみましょう。


 ■「ツィオルコフスキーの公式」を解く(2):重力に逆らって加速する場合 ――  

■さて、ココでは、外力 $f$ として、ロケットの推力 $F_{\mathrm{T}}$ とは逆向きに、重力 $f = -m(t)g$($\, g$ は重力加速度)がはたらいている場合を、考えます。簡単のために、「重力加速度 $g$ は、一定とする」と仮定しましょう。この仮定は、少し、現実的でない気がするかも知れません。地表から打ち上げられたロケットが、重力とは逆向きに加速していけば、地表から離れるほど、重力は弱くなるハズだからです。でも、心配は要りません。たとえば、国際航空連盟の定款(2022 年版)では、「宇宙飛行(Astronautics)」を、「地表 $100$ [$\mathrm{km}$] 以上で行われる活動」と定義しています。そして、地表 $100$ [$\mathrm{km}$] における重力加速度(およそ $9.5$ [$\mathrm{m/s^{2}}$])は、地表における重力加速度(およそ $9.8$ [$\mathrm{m/s^{2}}$])に比べて、わずか $3 \, \%$ しか違いません。すなわち、地表から打ち上げられたロケットが、宇宙に到達するまでの運動を考えるとき、「重力加速度 $g$ は、一定とする」と仮定するコトは、むしろ、「合理的な近似である」と言うべきなのです。

■そんなワケで、推力とは逆向きに、重力加速度 $g$ が一定の重力がはたらいているとき、ツィオルコフスキーの公式=(2.6)式は、 $f = -m(t)g$、および、(2.8b)式を代入して、次のようになります。

\[ -( m_{0} - \mu t ) \, g + \omega_{\mathrm{P}} \mu = ( m_{0} - \mu t ) \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{2.18} \]
(2.18)式より、ただちに、$\omega_{\mathrm{P}}\mu > m_{0 \,} g$、すなわち、$F_{\mathrm{T}} > m_{0 \,} g$ という条件が満たされていなければ、そもそも、ロケットは、地表から飛び立てない(時刻 $t =0 $ のとき、$du_{\mathrm{R}}(t)/dt$ が正の値にならない)コトが分かります(コレは、ツィオルコフスキーの公式が、常識と一致するコトを示す一例と言えます)。先に想定した数値を代入すると、推進剤の消費速度 $\mu$ は、少なくとも、およそ $56$ [$\mathrm{kg/s}$](工事用の排水ポンプが 1 台あれば、十分可能な値)でなければいけない――というコトになります。また、このとき、推進剤を使い切るまでの時間 $t_{\mathrm{P}}$ は、およそ $285$ [$\mathrm{s}$] です。ただし、コレらは、外力がない場合の想定を用いて、得られた結果です。重力がある場合における見積もりを得るために、計算を続けるとしましょう。


■それでは、(2.18)式を、実際に解いてみましょう。この微分方程式も、外力がない場合と同様に、簡単に解くコトができます。ただし、今回は、推進剤を使い切る前後で、場合分けをして解いてみましょう。先ほどの計算と同じように、初期条件として、時刻 $t = 0$ のときのロケットの速さ(初速度)$u_{\mathrm{R0}}$ を、$u_{\mathrm{R0}} = 0$ とすると(いろいろな解法や表記法がありますけど、一例として)、

\[ t \leq t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ \int_{0}^{u_{\mathrm{R}}(t)}{du_{\mathrm{R}}} = \int_{0}^{t}{\left( -g + \frac{\omega_{\mathrm{P}} \mu}{m_{0} - \mu t} \right) \, dt} \] \[ \qquad \mbox{したがって} \quad u_{\mathrm{R}}(t) = -gt - \omega_{\mathrm{P}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \, t \right)} \quad = -gt + \omega_{\mathrm{P}} \, \ln{\left( \frac{m_{0}}{m(t)} \right)} \tag{2.19a} \] \[ t > t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ \int_{0}^{u_{\mathrm{R}}(t)}{du_{\mathrm{R}}} = \int_{0}^{t_{\mathrm{P}}}{\left( -g + \frac{\omega_{\mathrm{P}} \mu}{m_{0} - \mu t} \right) \, dt} + \int_{t_{\mathrm{P}}}^{t}{\left( -g \right) \, dt} \] \[ \qquad \mbox{したがって} \quad u_{\mathrm{R}}(t) = -gt + \omega_{\mathrm{P}} \, \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} \tag{2.19b} \]

となります。(2.19a)〜(2.19b)式より、ただちに、推進剤を使い切った瞬間(時刻 $t_{\mathrm{P} \,}$)に、ロケットは最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達し、そのときの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、

\[ u_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \left\{ \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} - \frac{m_{\mathrm{P}} g}{\omega_{\mathrm{P}} \mu} \right\} \quad = \omega_{\mathrm{P}} \left\{ \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} - \frac{m_{\mathrm{P}} g}{F_{\mathrm{T}}} \right\} \tag{2.20} \]

であるコトが分かります。外力がない場合の最大到達速度((2.12)式)と比べると、当然、重力に逆らって加速する場合のほうが、$u_{\mathrm{R \cdot MAX}} \:$ は小さくなっています。また、(2.19b)式から、外力がない場合と異なり、時刻 $t_{\mathrm{P}} = m_{\mathrm{P}}/\mu$ を過ぎると(つまり、推進剤を使い切ってしまうと)、ロケットの速さ $u_{\mathrm{R}}(t)$ は、時間の経過とともに、重力加速度 $g$ に比例して、遅くなっていくコトも明らかです。ココまでに得られた結果を、グラフに表したのが、次の図です。
図3:時刻とロケットの速さの関係:ニュートン力学・重力あり(約17KB)


□(2.20)式に、外力がない場合の想定で得ていた数値を代入した場合、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、およそ $5.2 \times 10^{3}$ [$\mathrm{m/s}$] にしかなりません。コレは、地球周回軌道に乗るために必要な $7.9 \times 10^{3}$ [$\mathrm{m/s}$] には足りません。推進剤の消費速度 $\mu$ の増加だけで、不足分を補おうとすると、$\mu = 1.8 \times 10^{3}$ [$\mathrm{kg/s}$] 程度になります。ココで想定している、水素の燃焼反応をエネルギー源とするロケットとして、日本が誇る最強の LE-7A エンジンでさえ、$\mu = 2.6 \times 10^{2}$ [$\mathrm{kg/s}$] 程度であるコトを考えると、コレは、かなり厳しい要求と感じます(ただし、LE-7A エンジンは、推進剤の噴射速度 $\omega_{\mathrm{P}}$ が、$\omega_{\mathrm{P}} = 4.3 \times 10^{3}$ [$\mathrm{m/s}$] 程度であり、ココで想定してきた $\omega_{\mathrm{P}} = 3 \times 10^{3}$ [$\mathrm{m/s}$] を、大きく上回ります)。外力がない場合でさえ、軌道速度に達するコトは、非常に厳しい技術的挑戦なのに、さらに、重力に逆らって宇宙へ行くコトの「無謀さ」は、常軌を逸しているようにさえ思えます。ツィオルコフスキー が指摘したように、「多段式ロケット」のアイデアがなければ、今でも、有人宇宙飛行は、現実していなかったでしょう(多段式ロケットについては、ぜひ、自分で調べてみて下さいね)。現在の宇宙開発の状況は、まさに、「奇跡」と呼ぶべきモノなのです。


■重力がある場合の、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ が分かったので、ロケットの各要素の改善が、$u_{\mathrm{R \cdot MAX}} \:$ に与える影響を、調べてみましょう。(2.20)式より、ただちに、次の 4 つの関係式が得られます。

\[ du_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \left\{ \ln{ \left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} \right\} \frac{d\omega_{\mathrm{P}}}{\omega_{\mathrm{P}}} \tag{2.21a} \] \[ du_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \left\{ \frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{1 + m_{\mathrm{P}}/m_{\mathrm{R}}} \, \frac{m_{0} g}{\omega_{\mathrm{P}} \mu} \right\} \frac{d\mu}{\mu} \tag{2.21b} \] \[ du_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \left\{ \frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{1 + m_{\mathrm{P}}/m_{\mathrm{R}}} \left( 1 - \frac{m_{0} g}{\omega_{\mathrm{P}} \mu} \right) \right\} \frac{dm_{\mathrm{P}}}{m_{\mathrm{P}}} \tag{2.21c} \] \[ du_{\mathrm{R \cdot MAX}} = \omega_{\mathrm{P}} \left\{ -\frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{1 + m_{\mathrm{P}}/m_{\mathrm{R}}} \right\} \frac{dm_{\mathrm{R}}}{m_{\mathrm{R}}} \tag{2.21d} \]

コレらの関係式から、$\omega_{\mathrm{P}}, \ \mu, \ m_{\mathrm{P}}, \ m_{\mathrm{R}}$ を、それぞれ $1 \, \%$ 向上させたときに、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ の増加量がイチバン大きいのは、推進剤の噴射速度 $\omega_{\mathrm{P}}$ の高速化、次いで、ロケット本体の質量 $m_{\mathrm{R}}$ の軽量化であるコトが分かります(証明は省略しますけど、地表から飛び立つために、$\omega_{\mathrm{P}}\mu > m_{0 \,} g$ が成立しているコトに注意)。他方で、搭載する推進剤の質量 $m_{\mathrm{P}}$ を増量すると、ロケットは重くなるので、重力に逆らって加速する場合に、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に与える影響が、気になるトコロです。(2.21c)式より、搭載する推進剤の質量 $m_{\mathrm{P}}$ を増量すると、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、必ず、増加するコトが分かります(ココでも、$\omega_{\mathrm{P}}\mu > m_{0 \,} g$ に注意)。ただし、ロケット本体の質量 $m_{\mathrm{R}}$ の軽量化のほうが、効果は大きいので、搭載する推進剤の質量 $m_{\mathrm{P}}$ の増量のために、推進剤タンクを大型化するなどして、ロケット本体の質量 $m_{\mathrm{R}}$ が増加してしまうと、逆効果になりかねません。もし、あなたがロケット工学者なら、(2.21a)〜(2.21d)式を、注意深く比較・検討して、ロケットの性能向上のための、ベストな方法を探らなければいけない――というコトになりそうです。コレは、とても難しく、悩ましい仕事になるに、違いありません。


□ちなみに、推進剤の消費速度 $\mu$ の増加と、搭載する推進剤の質量 $m_{\mathrm{P}}$ の増量で、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に与える影響の大きさを比較すると、(2.21b)〜(2.21c)式より、推力 $F_{\mathrm{T}}$ が $2m_{0 \,} g$ 以下のときには、推進剤の消費速度 $\mu$ の増加の影響のほうが大きく、推力 $F_{\mathrm{T}}$ が $2m_{0 \,} g$ 以上のときには、搭載する推進剤の質量 $m_{\mathrm{P}}$ の増量の効果のほうが大きくなるコトが、分かります。つまり、推力 $F_{\mathrm{T}}$ の大きさが十分でない場合は、推進剤の消費速度 $\mu$ を増加するコトに努めるべきだし、推力 $F_{\mathrm{T}}$ が十分に大きいなら、搭載する推進剤の質量 $m_{\mathrm{P}}$ を増量するコトに注力したほうが良い――というワケです。とは言え、いずれにせよ、推進剤の噴射速度 $\omega_{\mathrm{P}}$ の高速化の効果のほうが、圧倒的に大きいというコトは、忘れてはいけません。


■ロケットの性能向上について、別の観点から、考えてみましょう。先ほどの「時刻 $t$ とロケットの速さ $u_{\mathrm{R}}(t)$ の関係を表すグラフ」には、外力がない場合と、重力に逆らって加速する場合の、両方が描かれています。図を見ると、重力がなければ到達できたハズの速さに比べて、重力に逆らって加速する場合に、ロケットの速さが、どれだけ「失われるか」が、良く分かります。いま、「重力によって失われる速さ」の基準として、外力がない場合のロケットの最大到達速度と、重力に逆らって加速する場合のロケットの最大到達速度の差を考えて、コレを、「最大到達速度の重力損失 $u_{\mathrm{R \cdot LOSS}} \ \ $」と呼ぶコトにしましょう(注:コレは、藤野の勝手な造語です)。すると、最大到達速度の重力損失 $u_{\mathrm{R \cdot LOSS}} \:$ は、(2.12)式と(2.20)式より、ただちに、

\[ u_{\mathrm{R \cdot LOSS}} = g \, \frac{m_{\mathrm{P}}}{\mu} \tag{2.22} \]

となるコトが分かります。見方を変えれば、重力に逆らって加速するロケットは、最大到達速度に関して、この重力損失の分だけ、「性能が落ちる」と、みなすコトができます。外力がない場合でさえ、軌道速度に到達するコトは、とても困難なのですから、この「性能低下=重力損失」を、できるだけ小さく抑えたいと考えるのは、当然です。(2.22)式より、重力損失を小さくするには、搭載する推進剤の質量 $m_{\mathrm{P}}$ を減量するか、推進剤の消費速度 $\mu$ を増加するかの、どちらかしかありません。ところが、搭載する推進剤の質量 $m_{\mathrm{P}}$ を減量すると、(2.21c)式より、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ が減少してしまいます。コレでは、本末転倒です。結局、重力損失を小さくするには、推進剤の消費速度 $\mu$ を増加させるしかありません。このように、外力がない場合には、ほとんど考慮する必要のなかった、推進剤の消費速度 $\mu$ が、重力に逆らって加速する場合には、とても重要な要素となります。重力中を飛行するロケットも、外力がない宇宙空間を飛行するロケットも、飛行原理は同じなのに、その設計思想は、かなり異なったモノにならざるをえない――というコトが、良く分かります。

■ところで、実際のロケットの運用を考えると、地表から打ち上げられたロケットが、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達した瞬間(つまり、時刻 $t_{\mathrm{P}}$)に、目的の地球周回軌道の高さまで到達していない、または、通過していた――というのでは、非効率です。理想的には、ロケットが最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達した瞬間に、ちょうど、目的の地球周回軌道の高さまで上昇している――という設計が、望まれます。最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達する時刻 $t_{\mathrm{P}}$ におけるロケットの位置 $x_{\mathrm{P}}$ は、簡単に計算できます。初期条件として、時刻 $t = 0$ のときのロケットの位置(初期位置)$x_{0}$ を、$x_{0} = 0$ とすると、時刻 $t_{\mathrm{P}}$ のときのロケットの位置 $x_{\mathrm{P}}$ は、次のようになります。

\[ x_{\mathrm{P}} = \int_{0}^{t_{\mathrm{P}}}{u_{\mathrm{R}}(t) \, dt} = \int_{0}^{t_{\mathrm{P}}}{\left\{ -gt - \omega_{\mathrm{P}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \, t \right)} \right\} \, dt} \] \[ \quad \mbox{したがって} \quad x_{\mathrm{P}} = \omega_{\mathrm{P}} \, \frac{m_{\mathrm{P}}}{\mu} - \omega_{\mathrm{P}} \, \frac{m_{\mathrm{R}}}{\mu} \, \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} - \frac{1}{2} \, g \, \left( \frac{m_{\mathrm{P}}}{\mu} \right)^{2} \tag{2.23} \]

もし、あなたがロケット工学者なら、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を与える(2.20)式と、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達したときのロケットの位置 $x_{\mathrm{P}}$ を与える(2.23)式の両方を睨みながら、具体的に、どのようなロケットを目指し、ナニを重視して開発するのか、そのグランドデザインを検討するコトになるのかも知れませんね。


□もちろん、地表から打ち上げられたロケットが、ずっと、垂直に(重力と逆向きに)上昇している限り、地球周回軌道には入れません。実際の地球周回軌道への投入では、重力の向きと推力の向きの関係は、非常に複雑になります。でも、簡単のため、ココでは、このまま、単純化した設定で、計算を続けましょう。

□「ラグランジュの未定乗数法」を知っている人は、「目指す地球周回軌道の高さ $x_{\mathrm{P}}$ が与えられたなら、(2.23)式を満たし、かつ、(2.20)式の最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を最大化する、最適な $\omega_{\mathrm{P}}, \ \mu, \ m_{\mathrm{P}}, \ m_{\mathrm{R}}$ の組み合わせを、求めるコトができるのでは?」と、考えるかも知れません。けれど、実際に試してみると、条件を満たすような $\omega_{\mathrm{P}}, \ \mu, \ m_{\mathrm{P}}, \ m_{\mathrm{R}}$ の組み合わせは、存在しない(専門用語で言うなら、「少なくとも、$\omega_{\mathrm{P}}, \ \mu, \ m_{\mathrm{P}}, \ m_{\mathrm{R}} > 0$ を満たす場所には、停留点は、存在しない」)というコトが、すぐに分かります。

□興味深いコトに、$m_{\mathrm{P}}$ と $m_{\mathrm{R}}$ の値は固定されたモノとして、$\omega_{\mathrm{P}}, \ \mu$ の組み合わせを調べると、(2.23)式を満たしたまま、いくらでも、$u_{\mathrm{R \cdot MAX}} \:$ を大きくできる(有限の最大値が存在しない)コトが分かります。逆に、(2.23)式を満たし、$u_{\mathrm{R \cdot MAX}} \:$ を最小化するような、(有限の値を持つ)$\omega_{\mathrm{P}}, \ \mu$ の組み合わせなら、存在します(詳しくは、長くなるので割愛します)。そして、$\mu$ の値にかかわらず、$\omega_{\mathrm{P}}$ が、
\[ \qquad \omega_{\mathrm{P}} \ge \frac{m_{\mathrm{P}}/m_{\mathrm{R}}}{m_{\mathrm{P}}/m_{\mathrm{R}} - \ln{( 1 + m_{\mathrm{P}}/m_{\mathrm{R}} )}} \, \sqrt{\mathstrut 2g \, x_{\mathrm{P}}} \tag{2.24} \]
という条件を満たしていないと、ロケットは、推進剤を使い切る時刻 $t_{\mathrm{P}}$ までに、高さ $x_{\mathrm{P}}$ に到達できない――というコトが、分かります。つまり、(2.24)式が、条件を満たす、最小の $\omega_{\mathrm{P}}$ を与えます。

□さらに、今度は、$\omega_{\mathrm{P}}$ と $\mu$ の値が固定されたモノとして、$m_{\mathrm{P}}, \ m_{\mathrm{R}}$ の組み合わせを調べると、(2.23)式を満たしたまま、いくらでも、$u_{\mathrm{R \cdot MAX}} \:$ を大きくできる(有限の最大値が存在しない)コトが分かります。逆に、(2.23)式を満たし、$u_{\mathrm{R \cdot MAX}} \:$ を最小化するような、(有限の値を持つ)$m_{\mathrm{P}}, \ m_{\mathrm{R}}$ の組み合わせは、存在します(その組み合わせは、$\omega_{\mathrm{P}}\mu = m_{0 \,} g$ を満たします)。そして、$m_{\mathrm{P}}$ の値にかかわらず、$m_{\mathrm{R}}$ が、
\[ \qquad m_{\mathrm{R}} \, \left\{ \frac{1}{2} \, \frac{m_{\mathrm{R}} g}{\omega_{\mathrm{P}} \mu} - \frac{1}{2} \, \frac{\omega_{\mathrm{P}} \mu}{m_{\mathrm{R}} g} + \ln{\left( \frac{\omega_{\mathrm{P}} \mu}{m_{\mathrm{R}} g} \right)} \right\} \le \frac{\mu}{\omega_{\mathrm{P}}} \, x_{\mathrm{P}} \tag{2.25} \]
という条件を満たしていないと、ロケットは、推進剤を使い切る時刻 $t_{\mathrm{P}}$ までに、高さ $x_{\mathrm{P}}$ に到達できない――というコトが、分かります。つまり、(2.25)式が、条件を満たす、最大の $m_{\mathrm{R}}$ を与えます。

このように、重力に逆らって加速するロケットでは、目的や、現有技術の水準などによって、最適となる $\omega_{\mathrm{P}}, \ \mu, \ m_{\mathrm{P}}, \ m_{\mathrm{R}}$ の組み合わせが、さまざまに、変わりえます。したがって、設計の方向性や、開発の方針についての、「唯一の正しい答え」といったモノは、存在しないと思われます。実際に、現在、さまざまな目的のために、さまざまな設計思想のもとで開発された、たくさんの種類のロケットが、活躍しています。ココでは、単段式ロケットで、地球周回軌道を目指すコトを、想定してきましたけど(それは、技術的に、あまりに非現実的であるコトが、明らかになりましたけど)、多段式ロケットや、弾道ロケットなどでは、別の設計思想が、必要になってくるに違いありません。ぜひ、具体的に、比較・検討してみて下さいね。


■当然のコトですけれど、推進剤を使い切る時刻(つまり、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達する時刻)$t_{\mathrm{P}}$ におけるロケットの位置 $x_{\mathrm{P}}$ は、ロケットの最大到達高度 $x_{\mathrm{MAX}} \:$ とは、異なります。いま、初期条件として、時刻 $t = 0$ のときのロケットの位置(初期位置)$x_{0}$ を、$x_{0} = 0$ とすると、時刻 $t$ のときのロケットの位置 $x(t)$ は、簡単に計算できます。結果だけを示すと、次のようになります。

\[ t \leq t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ x(t) = -\frac{1}{2}gt^{2} + \omega_{\mathrm{P}} t + \omega_{\mathrm{P}} \frac{m_{0}}{\mu} \left( 1 - \frac{\mu}{m_{0}} \, t \right) \ln{\left( 1 - \frac{\mu}{m_{0}} \, t \right)} \tag{2.26a} \] \[ t > t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ x(t) = -\frac{1}{2}gt^{2} + \omega_{\mathrm{P}} \frac{m_{\mathrm{P}}}{\mu} - \omega_{\mathrm{P}} \frac{m_{0}}{\mu} \left( 1 - \frac{\mu}{m_{0}} \, t \right) \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} \] \[ \hspace{12.4em} = -\frac{1}{2} \, g \, \left\{ t - \frac{\omega_{\mathrm{P}}}{g} \, \ln{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)} \right\}^{2} + x_{\mathrm{P}} + \frac{1}{2g} \, u_{\mathrm{R \cdot MAX}}^{2} \tag{2.26b} \]

(2.26b)式より、ロケットの最大到達高度 $x_{\mathrm{MAX}} \:$ は、ただちに、

\[ x_{\mathrm{MAX}} = x_{\mathrm{P}} + \frac{1}{2g} \, u_{\mathrm{R \cdot MAX}}^{2} \tag{2.27} \]

で与えられるコトが分かります。最後に、ロケットが重力に逆らって加速する場合の、打ち上げからの時刻 $t$ と、ロケットの速さ $u_{\mathrm{R}}(t)$、および、位置 $x(t)$ の関係を、グラフにまとめると、次の図のようになります。
図4:時刻とロケットの速さ・位置の関係:ニュートン力学・重力あり(約20KB)


□実際のロケット・宇宙開発の現場では、「比推力」、「LEO(Low Earth Orbit:地球低軌道)打ち上げ能力」、「GTO(Geostationary Transfer Orbit:静止トランスファ軌道)打ち上げ能力」など、ココでは紹介しなかった用語が、たくさん使われています。ココまでのハナシを理解できたなら、コレらの概念が、ココまでのハナシと、どのように関わり、なぜ必要なのか――きっと、より深く、より自然に、納得できるハズです。詳しくは、ぜひ、自分で調べて、コレらの専門用語を、使いこなして下さいね。


 ■「ツィオルコフスキーの公式」を解く(3):空気抵抗がある場合 ――  

■「ニュートン力学的に正しいツィオルコフスキーの公式」の最後に、空気抵抗がある場合の、ロケットの運動について、考えてみましょう。大気中を飛行するロケットにはたらく力(「抗力」と呼ばれます)は、一般に、非常に複雑な問題です。ただし、最も単純化したモデルでは、ロケットの速さ $u_{\mathrm{R}}$ の $2$ 乗に比例する抗力が、ロケットにはたらく――と、みなすコトができます(コレを、「慣性抵抗」と呼びます)。そこで、ココでは、ロケットに、慣性抵抗 $f = -\kappa \, u_{\mathrm{R}}^{2}$($\kappa$ は比例定数)がはたらく場合を、考えてみるコトにします。


□ちなみに、慣性抵抗は、運動量の保存則から、導かれます。いま、断面積 $S$ のロケットが、速さ $u_{\mathrm{R}}$ で飛行しているとしましょう。このロケットが、時刻 $t \,$〜$\, t+\varDelta t$ のあいだに、衝突する空気の体積 $\varDelta V$ は、$\varDelta V = Su_{\mathrm{R}}\varDelta t$ となります。ロケットとの衝突によって、この空気は、ロケットに押しつけられ、ロケットと同じ速さ $u_{\mathrm{R}}$ で運動するコトになる――と仮定します。空気の密度を $\rho$ とすると、空気の運動量変化 $\varDelta p$ は、$\varDelta p = \rho Su_{\mathrm{R}}\varDelta t \cdot u_{\mathrm{R}} = \rho Su_{\mathrm{R}}^{2}\varDelta t$ で表されます。空気の運動量変化 $\varDelta p$ は、空気に加えられた力積 $f_{\mathrm{AIR} \,}\varDelta t$ に等しいコトから、この衝突で空気が受けた力 $f_{\mathrm{AIR} \,}$ は、$f_{\mathrm{AIR} \,} = \rho Su_{\mathrm{R}}^{2}$ となるコトが分かります。「作用・反作用の法則」より、このとき、ロケットにはたらく力 $f$ は、$f = -f_{\mathrm{AIR} \,} = -\rho Su_{\mathrm{R}}^{2}$ となります。ココで、定数 $\kappa \equiv \rho S$ と定義すれば、慣性抵抗の公式が得られます。もちろん、コレは、実際に起きている現象を、あまりに単純化しすぎている――と、思われるかも知れません。でも、この公式は、実測値と、おおむね、一致するようです。したがって、以下では、この慣性抵抗の公式を使いたいと思います。なお、地表では、大気密度 $\rho \simeq 1.23$ [$\mathrm{kg}/\mathrm{m}^{3}$] なので、H-IIA ロケット(断面積 $S \simeq 12.6$ [$\mathrm{m}^{2}$])の場合は、$\kappa \simeq 15.5$ [$\mathrm{kg}/\mathrm{m}$]、イプシロンロケット(断面積 $S \simeq 5.31$ [$\mathrm{m}^{2}$])の場合は、$\kappa \simeq 6.53$ [$\mathrm{kg}/\mathrm{m}$] となります。


■さて、慣性抵抗 $f = -\kappa \, u_{\mathrm{R}}^{2}$ がはたらくとき、ツィオルコフスキーの公式=(2.6)式は、

\[ -\kappa u_{\mathrm{R}}^{2} + \omega_{\mathrm{P}} \mu = ( m_{0} - \mu t ) \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{2.28} \]

となります。以下では、簡単のため、慣性抵抗の公式に出てくる係数 $\kappa$ は、一定とします。もちろん、空気が薄くなれば、慣性抵抗は小さくなります。コレは、係数 $\kappa$ が小さくなるコトを、意味します。でも、そのコトを考慮に入れた計算は、とてもタイヘンです。そこで、ココでは、係数 $\kappa$ が一定の場合だけを、計算します。係数 $\kappa$ が一定のとき、(2.28)式の微分方程式は、簡単に解くコトができます。初期条件として、時刻 $t = 0$ のときのロケットの速さ(初速度)$u_{\mathrm{R0}}$ を、$u_{\mathrm{R0}} = 0$ とすると、推進剤を使い切るまでについては(いろいろな解法や表記法がありますけど、一例として)、

\[ t \leq t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ \int_{0}^{u_{\mathrm{R}}(t)}{\left( \frac{1}{\omega_{\mathrm{P}} \mu - \kappa u_{\mathrm{R}}^{2}} \right) \, du_{\mathrm{R}}} = \int_{0}^{t}{\left( \frac{1}{m_{0} - \mu t} \right) \, dt} \tag{2.29} \]

ココで、「部分分数分解」と呼ばれる数学テクニックを用いると、

\[ \frac{1}{\omega_{\mathrm{P}} \mu - \kappa u_{\mathrm{R}}^{2}} = \frac{1}{2 \sqrt{\mathstrut \omega_{\mathrm{P}} \mu \kappa}} \left( \frac{1}{\sqrt{\omega_{\mathrm{P}} \mu / \kappa} - u_{\mathrm{R}}} + \frac{1}{\sqrt{\omega_{\mathrm{P}} \mu / \kappa} + u_{\mathrm{R}}} \right) \tag{2.30} \]

と変形できるので、(2.29)式の積分は、簡単になります。結果だけを示すと、ロケットの速さ $u_{\mathrm{R}}(t)$ は、

\[ t \leq t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ u_{\mathrm{R}}(t) = \sqrt{\frac{\omega_{\mathrm{P}} \mu}{\kappa}} \, \left\{ \frac{1 - \left( 1 - \frac{\mu}{m_{0}} \, t \right)^{2 \sqrt{\omega_{\mathrm{P}} \, \kappa / \mu}}}{1 + \left( 1 - \frac{\mu}{m_{0}} \, t \right)^{2 \sqrt{\omega_{\mathrm{P}} \, \kappa / \mu}}} \right\} \tag{2.31} \]

と表されます。もちろん、慣性抵抗がある場合も、ロケットが最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達するのは、推進剤を使い切った瞬間(時刻 $t_{\mathrm{P} \,}$)で、そのときの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、次の式で与えられます。

\[ u_{\mathrm{R \cdot MAX}} = \sqrt{\frac{\omega_{\mathrm{P}} \mu}{\kappa}} \, \left\{ \frac{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)^{2 \sqrt{\omega_{\mathrm{P}} \, \kappa / \mu}} - 1}{\left( 1 + \frac{m_{\mathrm{P}}}{m_{\mathrm{R}}} \right)^{2 \sqrt{\omega_{\mathrm{P}} \, \kappa / \mu}} + 1} \right\} \tag{2.32} \]

なお、推進剤を使い切ったあとの、ロケットの速さ $u_{\mathrm{R}}(t)$ は、結果だけを示すと、次のとおりとなります。

\[ t > t_{\mathrm{P}} = \frac{m_{\mathrm{P}}}{\mu} \mbox{ のとき:} \ u_{\mathrm{R}}(t) = \left\{ \frac{\kappa}{m_{\mathrm{R}}} \left( t - \frac{m_{\mathrm{P}}}{\mu} \right) + \frac{1}{u_{\mathrm{R \cdot MAX}}} \right\}^{{\large -1}} \tag{2.33} \]


■ところで、推進剤を使い切るまでのロケットの速さ $u_{\mathrm{R}}(t)$((2.31)式)を眺めていると、興味深いコトに、気がつきます。まず、$1-(\mu/m_{0}) \, t > 0$ は明らかですので((2.8b)式などを参照)、$du_{\mathrm{R}}(t)/dt > 0$ が、常に成り立つコトは、すぐに分かります。コレは、$u_{\mathrm{R}}(t)$ が、単調増加関数であるコトを、表しています。さらに、$d^{2}u_{\mathrm{R}}(t)/dt^{2}$ については、次の関係が成り立っているコトが、(計算は、少しタイヘンですけど)分かります。

\[ \frac{1}{2} \sqrt{\frac{\mu}{\omega_{\mathrm{P}} \kappa}} \ge 1 \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} > 0 \tag{2.34a} \] \[ \frac{1}{2} \sqrt{\frac{\mu}{\omega_{\mathrm{P}} \kappa}} < 1 \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} \ge 0 \quad \left( \ t \le \frac{m_{0}}{\mu} \left\{ 1 - \left( \frac{1 - \frac{1}{2} \sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}}{1 + \frac{1}{2} \sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}} \right)^{\frac{1}{2}\sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}} \ \right\} \ \right) \] \[ \hspace{10.5em} \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} < 0 \quad \left( \ t > \frac{m_{0}}{\mu} \left\{ 1 - \left( \frac{1 - \frac{1}{2} \sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}}{1 + \frac{1}{2} \sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}} \right)^{\frac{1}{2}\sqrt{\mu / \omega_{\mathrm{P}} \, \kappa}} \ \right\} \ \right) \] \[ {} \tag{2.34b} \]

外力がない場合も、重力に逆らって加速する場合も、推進剤を使い切るまで、ロケットの加速度は、時間の経過とともに、大きくなっていきます。でも、慣性抵抗の係数 $\kappa$ が十分に大きい場合には、ある時刻までは、時間の経過とともに、大きくなっていた加速度が、その時刻を境に、逆に、時間の経過とともに、小さくなっていく――そのような特別な時刻、すなわち、「変曲点」が存在するコトになります(ちなみに、変曲点に到達する前に、推進剤を使い切る時刻 $t_{\mathrm{P}}$ に達してしまう場合も、ありえます)。ココで、再び、(2.31)式を見返すと、慣性抵抗がある場合、ロケットの速さ $u_{\mathrm{R}}$ には、いわゆる「終端速度(terminal velocity)」$u_{\mathrm{T}} = \sqrt{\smash[b]{\mathstrut \omega_{\mathrm{P}}\mu/\kappa}}$ が存在するコトが、分かります。慣性抵抗の係数 $\kappa$ が十分に大きい場合は、ロケットの速さ $u_{\mathrm{R}}$ は、漸近的に、終端速度に近づいていきます。でも、慣性抵抗の係数 $\kappa$ が小さい場合には、終端速度に達する前に、推進剤が尽きてしまいます。いずれにせよ、ロケットは、終端速度を超えるコトができません。そんなワケで、ココまで計算してきた結果を、分かりやすく示したのが、下の図です。
図5:時刻とロケットの速さの関係:ニュートン力学・空気抵抗あり(約22KB)

■以上で、「ニュートン力学的に正しいツィオルコフスキーの公式」の復習は、おしまいです。ココで検討したコトの他にも、検討したいコト、検討すべきコトは、たくさん、あると思います。でも、キリがないので、次の節から、いよいよ、本題である「特殊相対論的に正しいツィオルコフスキーの公式」の検討に、進みましょう。


■このページの一番上に戻る。



 ■3. 必要な準備(2) ―― 特殊相対論における基本的な関係式

■さて、ココからは、「特殊相対論的に正しいツィオルコフスキーの公式」について、考えていきます。まず、必要となる、特殊相対論の基本的な関係式を、簡単に、おさらいしておきましょう。ココでは、「特殊相対論的に正しいツィオルコフスキーの公式」を導くのに、最低限、必要となる事柄のみを、復習するコトにします。

■まず最初に、「ロケットの運動を観測する観測者」が誰なのかを、決める必要があります。観測者の運動状態によって、観測されるロケットの運動が、変化してしまうからです(たとえば、地球の運動を観測する場合に、地表に立つ観測者と、月にいる観測者と、太陽系の遥か彼方にいる観測者では、観測結果が、それぞれ異なったモノになるのと、同じです)。ニュートン力学の世界ならば、通常、「静止している観測者」を基準とします。でも、特殊相対論の世界では、「ホントに静止しているのは、誰なのか?」という問題が、生じます。そこで、「少なくとも、加速度運動はしていない観測者」を選びます。加速度運動をしている観測者には、「慣性力」がはたらくので、「ホントに加速度運動をしていないのは、誰なのか?」という問題が、生じないからです。この観測者を、「慣性系にいる観測者」と呼ぶコトにしましょう(なぜ、そう呼ぶのかは、ぜひ、自分で調べてみて下さい)。ちなみに、ある慣性系にいる観測者から見て、さまざまな速さで運動する別の慣性系が、無数に存在します。異なる慣性系にいる観測者が観測すると、同じロケットでも、速さは異なる値になります。でも、コレは、本質的な違いではありません。初期条件で決まる定数を加えるだけで、まったく同じ結果になるからです。

■次に、最も重要なのは、特殊相対論の世界では、ロケットに乗っている観測者が観測する時刻 $\tau$ は、慣性系にいる観測者が観測する時刻 $t$ とは、異なる――というコトです(以後、ロケットに乗っている観測者が観測する物理量は、できるだけギリシャ文字で表し、慣性系にいる観測者が観測する物理量は、できるだけローマ字で表すコトで、区別するようにします)。たとえば、出発時に、時計を合わせておいても、ロケットが目的地に着いたとき、ロケットに乗っている観測者が持つ時計では、$60$ [分間] しか経過していないのに、慣性系にある時計では、$62.6$ [分間] が経過している――というコトが、起こりえます。よって、「いま、考えている物理量は、どの観測者によって観測されるモノか」を、常に、ハッキリと意識しておくコトが、とても大切です。


□ちなみに、上の「ウラシマ効果」の例は、ロケット本体の質量と同じだけの質量の推進剤を、ロケットに乗っている観測者から見て、光速の $0.8$ 倍の速さで噴射し、$60$ [分間] で使い切った、その瞬間における、お互いの時計が示す時刻を、計算してみた結果です。ただし、化学ロケットの場合、推進剤の噴射速度は、光速の $5.8$ 万分の $1$ 倍が、原理的な限界です。実際には、「光速の $0.8$ 倍」という噴射速度を、実現するエネルギー源と推進剤は、なかなか、見つけられないかも知れません。

□ところで、特殊相対論の世界における、時刻の測定について、ココで、改めて、確認しておきましょう。2 つの時計が、同じ場所にあれば、コレらの時計の時刻を比較するのに、ナンの問題もありません。でも、2 つの時計が、離れた場所(たとえば、地球と火星)にある場合、コレらの時計の時刻を比較するコトは、簡単ではありません。光の進む速さが有限なので、「まさに、いま、この瞬間に、時計が示している時刻」の情報が、観測者に届くまでに、時間が経過し、時刻が変化してしまうからです。もし、離れた場所にある時計と観測者の位置関係が、変化しないのであれば――言い換えると、「離れた場所にある時計が、観測者に対して静止しているならば」、この問題を解決するコトができます(詳しくは、『「動いている時計は、遅くなる」の間違いを暴く!』を参照)。

□いま、空間のあらゆる場所に、時計が置いてあるとしましょう。コレらの時計は、すべて、慣性系にいる観測者に対して静止しているとします。お互いに静止しているので、慣性系にいる観測者が持つ時計と、他のすべての時計は、時刻を正しく比較するコトができます。そこで、コレらすべての時計は、慣性系にいる観測者が持つ時計と、正確に、時刻を合わせてあるモノとして、「慣性系にある時計」と名づけましょう。もちろん、慣性系にある、どの 2 つの時計のあいだの距離も、レーザー測距の原理を用いて、正しく測定できます。この距離を、「慣性系で測った距離」と呼び、記号 $x$ で表すコトにしましょう。

□他方で、「加速度運動する観測者に対して、静止している時計」を、特殊相対論の世界で、矛盾なく定義するコトは、恐らく、できないと思われます(詳しくは、『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』を参照)。よって、ロケットに乗っている観測者には、自分が持つ時計の時刻と、離れた場所にある時計の時刻を、正しく比較する手段がありません。さらに、離れた場所にある時計までの距離も、測る手段はありません。つまり、「ロケットに乗っている観測者が測った距離」といったモノは、定義ができないのです。とは言え、いま、自分がいる場所に置いてある(自分と同じ場所にある)「慣性系にある時計」の時刻なら、(たとえ、一瞬で通り過ぎてしまうとしても)自分が持つ時計の時刻と、正しく比較するコトができます。さらに、たくさんの「慣性系にある時計」のうち、自分と同じ場所にあるのが、どの時計なのかを確認するコトで、「慣性系で測った距離」の情報も得られます。このように、特殊相対論の世界では、測定できる物理量と、できない物理量を、判断するコトが大切になります。また、ある観測者によって測定される物理量を、別の観測者が、どのように知り得るのかを考えるコトも、とても有益です。


■観測者によって、観測される時刻が異なるなら、当然、観測される速さや加速度も、観測者によって、異なるコトになります。以下、慣性系にいる観測者が観測するとき、ロケットの速さを $u_{\mathrm{R}}$、加速度を $a_{\mathrm{R}} = du_{\mathrm{R}}/dt$ と書き、推進剤の噴射速度を $u_{\mathrm{P}}$ と書くコトにします。また、ロケットに乗っている観測者が観測するときには、ロケットの速さを $\omega_{\mathrm{R}}$、加速度を $\alpha_{\mathrm{R}} = d\omega_{\mathrm{R}}/d\tau$ と書き、推進剤の噴射速度を $\omega_{\mathrm{P}}$ と書くコトにします。ただし、ロケットに乗っている観測者にとって、自分とロケット本体のあいだの位置関係は変化しない(自分に対して、ロケットは静止している)ので、常に $\omega_{\mathrm{R}} = 0$ であるコトに、注意しましょう。一方で、ロケットに乗っている観測者は、推力が生じているとき、常に「慣性力」を感じているハズなので、$\alpha_{\mathrm{R}} = d\omega_{\mathrm{R}}/d\tau \ne 0$ となります。なお、ロケットに乗っている観測者は、自分と同じ場所にある「慣性系にある時計」の動き(速さ)を、逐次、観測するコトができます。この速さは、「相対性原理」より、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$ と、大きさが同じで、向きが逆になります。したがって、ロケットに乗っている観測者は、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$ を、いつでも、直接、知るコトができます。

■必要な物理量が揃ったので、コレらのあいだに成り立つ関係式を、与えましょう。最初に、慣性系にある時計の時刻が $t$ のときの、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ の値を与える関数を、$\tau = \tau(t)$ と書くコトにします。逆に、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、慣性系にある時計の時刻 $t$ の値を与える関数は、$t = t(\tau)$ と書きます(なお、関数 $t(\tau)$ の具体的なカタチは、第 6 節で、検討します)。このとき、時刻 $t$ と時刻 $\tau$ は、同じ瞬間を表しているコトになります(正確に表現すると、「ロケットがある場所における、同じ瞬間=時間軸上の 1 点」を表しています)。よって、あらゆる物理量の観測値に対して、

\[ m(t) = m\bigl( t(\tau) \bigr) = m(\tau) = m\bigl( \tau(t) \bigr) \ , \qquad u_{\mathrm{R}}(t) = u_{\mathrm{R}}\bigl( t(\tau) \bigr) = u_{\mathrm{R}}(\tau) = u_{\mathrm{R}}\bigl( \tau(t) \bigr) \ , \ \ldots \tag{3.1} \]

が成り立ちます。コレは、「(同じ場所で、)同じ瞬間に、同じ物理量を観測すれば、同じ観測値が得られる」という、当たり前のコトを、厳密に述べただけです。でも、このコトを、最初に確認しておくと、あとで、複雑な計算(特に、変数変換とか)をしていくときに、自信をもって、計算を進めていくコトができるでしょう。


□ちなみに、特殊相対論の世界では、「慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$」と「ロケットに乗っている観測者から見たロケットの速さ $\omega_{\mathrm{R}}$」は、「同じ物理量」ではありません。だから、もちろん、$u_{\mathrm{R}}(t) = \omega_{\mathrm{R}}(t)$ や $u_{\mathrm{R}}(t) = \omega_{\mathrm{R}}(\tau)$ などは、一般に、成り立ちません。他にも、時間と空間(距離)に関わる物理量は、一般に、観測者が異なると、同じ物理量として扱うコトはできませんので、注意が必要です。

□他方で、質量 $m$ や、電気量 $q$ は、時間と空間(距離)に依らない、質点に固有の属性(物理量)です。したがって、コレらは、観測者が異なっても(観測者の運動状態に依らずに)、「同じ物理量」として扱うコトができます。「慣性系にいる観測者から見た(任意の質点の)質量」と「ロケットに乗っている観測者から見た(同じ質点の)質量」は、別々の物理量であると考える必要がない――というワケです。

□ある物理量に対して、観測者が違う場合にも、「同じ物理量」として扱えるのかどうかは、特殊相対論的に正しい運動方程式や、マクスウェルの方程式を、ローレンツ変換した結果に、「相対性原理」を適用するコトで、明確に、判別できます。ただ、たとえ、判別できなくても、心配は要りません。もし、実際には、「同じ物理量」として扱えるのならば、別々の物理量であると仮定して計算したとしても、結局は、「両者は、常に同じ値でなければいけない(同じ物理量でなければいけない)」という結果が得られるからです。

□なお、(3.1)式が成り立つコトを説明したので、この節では、煩雑さを避けるために、時間変数の明示を、省略します(たとえば、$u_{\mathrm{R}}(t)$ や $u_{\mathrm{R}}(\tau)$ のようには書かずに、$u_{\mathrm{R}}$ とだけ書きます)。


■さて、慣性系にある時計の時刻 $t$ と、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ のあいだに、一般に成り立つ関係としては、次の関係式が、特に重要です(以下では、光の速さ(光速)を、記号 $c$ で表します。もちろん、「光速度不変の原理」より、光速 $c$ の値は、観測者に依らず、常に一定です)。

\[ \frac{d\tau}{dt} = \sqrt{1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2}} \qquad \mbox{および} \qquad \frac{dt}{d\tau} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \tag{3.2} \]

次に、慣性系にいる観測者から見たロケットの加速度 $a_{\mathrm{R}} = du_{\mathrm{R}}/dt$ と、ロケットに乗っている観測者が感じる加速度 $\alpha_{\mathrm{R}} = d\omega_{\mathrm{R}}/d\tau$ のあいだに成り立つ関係式は、以下のようになります。

\[ a_{\mathrm{R}} = \left\{ 1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2} \right\}^{3/2} \, \alpha_{\mathrm{R}} \qquad \mbox{すなわち} \qquad \frac{du_{\mathrm{R}}}{dt} = \left\{ 1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2} \right\}^{3/2} \, \frac{d\omega_{\mathrm{R}}}{d\tau} \tag{3.3} \]

ココに現れている $u_{\mathrm{R}}$ は、慣性系にいる観測者から見たロケットの速さですけど、上で説明したように、この値は、ロケットに乗っている観測者も、逐次、知るコトができます。このとき、ロケットに乗っている観測者は、自分が持つ時計の時刻 $\tau$ に基づいて、$u_{\mathrm{R}}$ の時間変化を、($du_{\mathrm{R}}/dt$ ではなくて)$du_{\mathrm{R}}/d\tau$ であると判断するハズです。(3.2)〜(3.3)式より、$du_{\mathrm{R}}/d\tau$ について、次の関係式が成り立つコトが分かります。

\[ \frac{du_{\mathrm{R}}}{d\tau} = \frac{dt}{d\tau} \, \frac{du_{\mathrm{R}}}{dt} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \, \frac{du_{\mathrm{R}}}{dt} \quad = \left\{ 1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2} \right\} \, \frac{d\omega_{\mathrm{R}}}{d\tau} \tag{3.4} \]
□(3.2)式の導出は、『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の(2.9)〜(2.10)式を参照。ただし、ココでは、ロケットに乗っている観測者から見た、ロケットの速さ $\omega_{\mathrm{R}}$(上記記事(2.9)〜(2.10)式の $\boldsymbol{\omega}$ に相当)が、常に $\omega_{\mathrm{R}} = 0$ であるコトに注意。(3.3)式の導出は、『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』の(4.1)式を参照。


■ココで、新たに、質量 $m$ の質点を考えて、この質点が、任意の運動をしているとしましょう。慣性系にいる観測者から見た、この質点の速さを $u$ で表し、ロケットに乗っている観測者から見た、この質点の速さを $\omega$ で表すコトにすると、$u$ と $\omega$ のあいだには、「速度の合成則」の名で良く知られた、次の関係式が成り立ちます(『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』の(2.4)式を参照)。

\[ u = \frac{u_{\mathrm{R}} + \omega}{1 + ( u_{\mathrm{R}} \, \omega/c^{2} )} \qquad \mbox{および} \qquad \omega = \frac{- \ u_{\mathrm{R}} + u}{1 - ( u_{\mathrm{R}} \, u/c^{2} )} \tag{3.5} \]


■また、慣性系にいる観測者から見た、この質点が持つ運動量を $p_{k \,}$、エネルギーを $K_{k}$ で表すと、

\[ \mbox{運動量:} p_{k} = \frac{m u}{\sqrt{1 - ( u/c )^{2}}} \, , \qquad \mbox{エネルギー:} K_{k} = \frac{mc^{2}}{\sqrt{1 - ( u/c )^{2}}} \tag{3.6} \]

となり、ロケットに乗っている観測者から見た、この質点が持つ運動量を $p_{\kappa}$、エネルギーを $K_{\kappa}$ で表すと、

\[ \mbox{運動量:} p_{\kappa} = \frac{m \omega}{\sqrt{1 - ( \omega/c )^{2}}} \, , \qquad \mbox{エネルギー:} K_{\kappa} = \frac{mc^{2}}{\sqrt{1 - ( \omega/c )^{2}}} \tag{3.7} \]

となるコトが知られています(『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の(7.1a)〜(7.1b)式を参照)。


□質点の速さがゼロ($u = 0$ または $\omega = 0$)の場合でも、質点の質量 $m$ がゼロでないならば、質点が持つエネルギーはゼロにならないコトに注意しましょう。すなわち、熱や光と同じように、質量も、それ自体が、エネルギーを持つのです。しかも、(3.6)〜(3.7)式により、そのエネルギーの大きさは、質点の速さとも、密接に関わっているコトが分かります。そこで、質量と速さの両方に依存する、このエネルギー($K_{k}$ および $K_{\kappa \,}$)を、「質量=運動エネルギー」と呼びましょう(注:コレは、藤野の造語です)。


■さらに、慣性系にいる観測者から見た運動量 $p_{k}$ およびエネルギー $K_{k}$ と、ロケットに乗っている観測者から見た運動量 $p_{\kappa}$ およびエネルギー $K_{\kappa}$ は、次のように、お互いに変換できます(『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の(7.8)〜(7.10)式を参照)。

\[ p_{k} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \, \left( p_{\kappa} + \frac{u_{\mathrm{R}}}{c^{2}} K_{\kappa} \right) \, , \qquad K_{k} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \, \bigl( K_{\kappa} + u_{\mathrm{R}} \, p_{\kappa} \bigr) \tag{3.8a} \] \[ p_{\kappa} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \, \left( p_{k} - \frac{u_{\mathrm{R}}}{c^{2}} K_{k} \right) \, , \qquad K_{\kappa} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \, \bigl( K_{k} - u_{\mathrm{R}} \, p_{k} \bigr) \tag{3.8b} \]


■最後に、慣性系にいる観測者から見た質点にはたらく力 $f_{k}$ と、ロケットに乗っている観測者から見た質点にはたらく力 $f_{\kappa}$ の関係は、一般に、とても複雑です。でも、質点の運動の向きと、力の向きが平行の場合には、両者は、同じ値になります(『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の(6.7)式を参照)。つまり、

\[ f_{k} = f_{\kappa} \qquad \mbox{(力の向きが、質点の運動の向きと平行のとき)} \tag{3.9} \]

が成り立ちます。たとえば、1 次元の運動の場合には、(3.9)式が、常に成り立つコトになります。ただし、ニュートン力学の世界と、特殊相対論の世界では、運動方程式のカタチが異なるコトに、注意が必要です。力の定義そのものに、関わるハナシだからです。とは言え、「特殊相対論的に正しいツィオルコフスキーの公式」の計算には、直接、影響するコトはないので、ココでは、深入りせずに、先に進みましょう(興味があるならば、『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の、特に、第 5 節と第 6 節を参照してみて下さい)。ちなみに、次節以降では、実際に、「特殊相対論的に正しいツィオルコフスキーの公式」を導きますけど、第 2 節で、「ニュートン力学的に正しいツィオルコフスキーの公式」を導いたときと同様に、簡単のため、1 次元の運動((3.9)式が成り立つ場合)だけを考えます。

■以上で、必要な事柄は揃いましたので、次の節から、いよいよ、「特殊相対論的に正しいツィオルコフスキーの公式」の導出に、挑んでみましょう。


■このページの一番上に戻る。



 ■4. 特殊相対論的に正しいツィオルコフスキーの公式(1) ―― 化学ロケット編

■前置きが、随分と長くなりましたけど、ようやく、ココから、「特殊相対論的に正しいツィオルコフスキーの公式」について、本格的な検討を始めるコトができます。まずは、「特殊相対論的に正しいツィオルコフスキーの公式」を、実際に、導いてみましょう。基本的な導き方は、ココまでに復習してきた、「ニュートン力学的に正しいツィオルコフスキーの公式」の導出と、ほとんど同じです。ただし、「エネルギーの保存則」の扱いが、両者で異なるコトが、決定的に重要です。推進剤を噴射するには、エネルギーが必要ですけど、ニュートン力学の世界では、推進剤を噴射するのに、どんなエネルギーを使うのかは考慮せずに、運動量の保存則を考えるコトができました。でも、特殊相対論の世界では、そうはいきません。とても大きな問題が、2 つ、あるからです。それぞれの問題について、詳しく検討してみましょう。

■1 つ目の問題は、特殊相対論の世界では、運動量とエネルギーが、分かちがたく、密接に絡み合っているコトです。(3.8a)〜(3.8b)式で示されたように、慣性系にいる観測者から見て、運動量 $p_{k}$ とエネルギー $K_{k}$ を持つ質点を、ロケットに乗っている観測者が見たときは、運動量 $p_{\kappa}$($p_{\kappa} \ne p_{k}$)とエネルギー $K_{\kappa}$($K_{\kappa} \ne K_{k}$)を持つ質点として、観測されます(実は、(3.8a)〜(3.8b)式は、質点だけでなく、電磁波が持つ運動量とエネルギーについても、成り立ちます)。問題は、それだけではありません。ニュートン力学の世界では、ある観測者から見て、運動量は一定で、エネルギーだけが変化する現象は、別の観測者が見ても、運動量は一定で、エネルギーだけが変化します。でも、特殊相対論の世界では、ロケットに乗っている観測者から見て、運動量は一定で、エネルギーだけが変化する現象を、慣性系にいる観測者が見ると、運動量も、エネルギーも、どちらも変化するのです((3.8a)式を参照)。ニュートン力学の世界とは異なり、特殊相対論の世界では、ロケットの運動を計算する場合に、運動量の保存則だけを考えれば良いというワケにはいかないに違いありません。

■2 つ目の問題は、特殊相対論の世界では、質量そのものが、エネルギーを持っているコトです。推進剤の噴射に使われるエネルギーが、質量=運動エネルギーに由来する場合には、コレは、本質的な問題です。たとえば、化学反応、核分裂、核融合、対消滅などでは、反応前後で、全体の質量が変化します。この質量変化に相当する質量=運動エネルギーが、反応エネルギーとして放出・吸収されるのです(化学反応の場合には、質量変化が、あまりにも小さいので、質量は変わらずに、エネルギー(化学エネルギー)だけが放出・吸収されているように見えますけど、特殊相対論の世界では、そんな、極めてわずかな質量変化さえ、無視するコトはできません)。そのため、推進剤を噴射するために、コレらのエネルギーを使う場合には、推進剤の噴射前と噴射後における、ロケットと推進剤の質量の関係を考えるときに、特に、注意が必要となります。

 ■基本的な仮定 ――  

■コレらの注意事項を踏まえて、第 2 節で、「ニュートン力学的に正しいツィオルコフスキーの公式」を導出したときの、基本的な設定について、特殊相対論の観点から、再検討してみましょう。第 2 節では、「推進剤の消費速度(単位時間当たりの推進剤の消費量)$\mu$ は、一定とする」と仮定しました。でも、特殊相対論の立場では、「誰から見て、一定なのか?」という問題が生じます。ココでは、推進剤の消費速度の制御は、ロケット本体で行う――と考えて、ロケットに乗っている観測者(が持つ時計の時刻 $\tau$)を基準にして、

\[ \mu \equiv -\frac{dm(\tau)}{d\tau} \quad \mbox{($\, \mu$ は定数)} \tag{4.1} \]

と仮定しましょう。つまり、「ロケットに乗っている観測者から見た推進剤の消費速度 $\mu$ は、一定とする」と仮定するワケです。もちろん、この $\mu$ の値は、慣性系にいる観測者から見た推進剤の消費速度 $-dm(t)/dt$ と、同じにはなりません。実際、(3.1)式と(3.2)式より、「合成関数の微分の法則」を用いて、

\[ -\frac{dm(t)}{dt} = -\frac{dm(\tau(t))}{dt} = -\frac{d\tau(t)}{dt}\frac{dm(\tau)}{d\tau} = \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \ \mu \quad \mbox{($\, \mu$ は定数)} \tag{4.2} \]

となるコトが分かります。(4.2)式は、たとえ、ロケットに乗っている観測者から見て、推進剤の消費速度 $\mu$ が一定であったとしても、慣性系にいる観測者から見たら、それは、もはや、一定ではない――というコトを、示しています。つまり、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(t)$ が大きくなるほど、慣性系にいる観測者から見た推進剤の消費速度 $-dm(t)/dt$ は小さくなるのです。さらに、ロケットの速さ $u_{\mathrm{R}}(t)$ が、光速 $c$ に到達すると、慣性系にいる観測者から見た推進剤の消費速度は、ゼロになります。推進剤を消費しなければ、当然、ロケットは、加速できません。よって、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(t)$ は、光速 $c$ を超えられない――というコトになります。


□なぜ、このようなコトが、起きるのでしょうか。実は、コレには、「ウラシマ効果」が影響しています。慣性系にある時計の時刻 $t$ と、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ のあいだには、(3.2)式の関係がありました。(3.2)式は、次のように書き直すコトができます。
\[ \qquad d\tau = \sqrt{1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2}} \, dt \tag{4.3} \]
ココで、$\sqrt{\smash[b]{\mathstrut 1-(u_{\mathrm{R}}/c)^{2}}}$ の値は、必ず、$1$ 以下であるコトに注意すると、ただちに、
\[ \qquad \lvert d\tau \rvert \leq \lvert dt \rvert \tag{4.4} \]
という関係が、必ず、成り立つコトが分かります。たとえば、慣性系にある時計で $1$ 秒が経過したときに($dt = 1$ [秒])、ロケットに乗っている観測者が持つ時計は $1$ 秒未満しか経過していない($d\tau \leq 1$ [秒])――というコトになります。すなわち、慣性系にある時計よりも、ロケットに乗っている観測者が持つ時計のほうが、必ず、進み方が遅くなるのです(「慣性系にいる観測者が持つ時計」ではなく、「慣性系にある時計」であるコトに注意。詳しくは、『「動いている時計は、遅くなる」の間違いを暴く!』を参照)。

□さらに、$u_{\mathrm{R}}= c$ のときには、$\sqrt{\smash[b]{\mathstrut 1-(u_{\mathrm{R}}/c)^{2}}} = 0$ となります。このとき、ロケットに乗っている観測者の立場で見た場合、ロケットに乗っている観測者が持つ時計の時間が少しでも経過すると($d\tau > 0$ [秒])、慣性系にある時計では無限大の時間が経過しなければいけないコトになります($dt = \infty$ [秒])。他方で、慣性系にいる観測者の立場では、慣性系にある時計の時間がどれだけ経過しても($dt > 0$ [秒])、ロケットに乗っている観測者が持つ時計の時間は経過しないコトになります($d\tau = 0$ [秒])。このコトは、当然、慣性系にいる観測者から見た推進剤の消費速度 $-dm(t)/dt$ にも、影響を与える――というワケです。


■同じく、第 2 節では、「推進剤の噴射速度は、一定とする」としました。コレも、推進剤の噴射速度の制御は、ロケット本体で行う――と考えて、「ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、一定とする」と仮定します。もちろん、$\mu$ の場合と同じように、$\omega_{\mathrm{P}}$ の値も、慣性系にいる観測者から見た推進剤の噴射速度 $u_{\mathrm{P}}$ と、同じにはなりません。コチラも、$\mu$ の場合と同じく、注意が必要です。$\omega_{\mathrm{P}}$ と $u_{\mathrm{P}}$ のあいだに成り立つ関係式については、次項以降で、具体的な状況設定を見ながら、検討するコトにします。

■以上の仮定のもとに、次の項では、ロケットが運動している状況を、具体的に描いてみて、「特殊相対論的に正しいツィオルコフスキーの公式」を、実際に、導いてみましょう。

 ■ロケットに乗っている観測者から見た場合の、ツィオルコフスキーの公式の導出 ――  

■まずは、簡単なほう――ロケットに乗っている観測者から見たロケットの運動の様子を、描いてみましょう。次の図のように、時刻 $\tau$ のときの、推進剤を含んだロケット全体の質量を $m(\tau)$、ロケット本体の速さを $\omega_{\mathrm{R}}(\tau)$ と表すコトにします。もちろん、すでに注意しておいたとおり、$\omega_{\mathrm{R}}(\tau) = 0$ です。時間が $\varDelta \tau$ だけ経過すれば、ロケットの質量と速さは、それぞれ、$m(\tau+\varDelta \tau)$、$\omega_{\mathrm{R}}(\tau+\varDelta \tau)$ になります。ココで、ロケットは、時間 $\varDelta \tau$ が経過するあいだに、質量 $\varDelta m_{\mathrm{P}}$ の推進剤を、速さ $-\omega_{\mathrm{P}}$ で噴射します(向きに注意)。ところで、特殊相対論の世界では、推進剤を噴射するためのエネルギーを考えるコトが、特に重要なのでした。ココ(第 4 節)では、質量=運動エネルギーに由来するエネルギーを用いて、推進剤を噴射するロケットを考えましょう。ちなみに、「化学ロケット」などが、コレに該当します。この場合、すでに説明しておいたとおり、反応エネルギーの放出の前後で、全体の質量が減少するコトになります。この反応で失われる質量を、「質量欠損(mass defect)」と呼び、記号 $\varDelta m_{\mathrm{D}}$ で表しましょう。最後に、ロケット本体には、ロケット自身の推力のほかにも、外力 $f$ がはたらくコトを、忘れずに、想定しておきます。以上の設定を図にすると、次のようになるハズです。
図6:「特殊相対論的に正しいツィオルコフスキーの公式」の導出:化学ロケット編(1)(約42KB)


□ちなみに、「ロケットに乗っている観測者から見て、ロケットの速さ $\omega_{\mathrm{R}}$ が、常に $\omega_{\mathrm{R}} = 0$ であるなら、$\omega_{\mathrm{R}}(\tau)$ だけでなく、$\omega_{\mathrm{R}}(\tau+\varDelta \tau)$ も、ゼロになるのでは?」と思ったかも知れません。実は、「どの時刻の観測者から見た速さなのか?」が、その疑問を解く鍵となります。確かに、「時刻 $\tau+\varDelta \tau$ の時点におけるロケットに乗っている観測者」にとっては、$\omega_{\mathrm{R}}(\tau+\varDelta \tau) = 0$ です。でも、いま、知りたいのは、「時刻 $\tau$ の、その瞬間における、ロケットの運動を表す微分方程式」です。だから、あくまでも、「時刻 $\tau$ の時点におけるロケットに乗っている観測者」を基準にして、考えなければいけません。すなわち、「この問題に関しては、$\omega_{\mathrm{R}}(\tau+\varDelta \tau)$ は、時刻 $\tau$ の時点におけるロケットに乗っている観測者から見た、時刻 $\tau+\varDelta \tau$ の時点におけるロケットの速さである」――というコトになります。


■さて、「ニュートン力学的に正しいツィオルコフスキーの公式」を導出したときと同じく、(1.2)式より、時刻 $\tau+\varDelta \tau$ での、推進剤を含んだロケット全体の質量 $m(\tau+\varDelta \tau)$ と、ロケットの速さ $\omega_{\mathrm{R}}(\tau+\varDelta \tau)$ は(ココでは、記号 $\mu$ を(4.1)式のように定義したコト、また、$\omega_{\mathrm{R}}(\tau) = 0$ であるコトも用いて)、

\[ m(\tau + \varDelta \tau) = m(\tau) + \frac{dm(\tau)}{d\tau} \cdot \varDelta \tau \quad = m(\tau) - \mu \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.5a} \] \[ \omega_{\mathrm{R}}(\tau + \varDelta \tau) = \omega_{\mathrm{R}}(\tau) + \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \cdot \varDelta \tau \quad = \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \, \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.5b} \]

と表せます。ただし、推進剤の噴射前後における質量の関係については、質量欠損 $\varDelta m_{\mathrm{D}}$ が存在するために、「ニュートン力学的に正しいツィオルコフスキーの公式」を導出したときとは異なるコトに、注意が必要です。「ニュートン力学的に正しいツィオルコフスキーの公式」における(2.2)式に相当する関係式は、ココでは、

\[ m(\tau + \varDelta \tau) = m(\tau) - \varDelta m_{\mathrm{P}} - \varDelta m_{\mathrm{D}} \tag{4.6} \]

となります。この(4.6)式に、(4.5a)式を代入すると、次の関係が成り立っているコトが、分かります。

\[ \varDelta m_{\mathrm{P}} = \mu \varDelta \tau - \varDelta m_{\mathrm{D}} \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.7} \]


■必要な準備が揃ったので、いま、考えた設定に対して、まず、「運動量の保存則」を、適用してみましょう。ただし、ニュートン力学の世界と異なり、特殊相対論の世界では、ロケットに乗っている観測者から見た運動量が、(3.7)式で与えられます。他方で、特殊相対論の世界でも、外力 $f$ がはたらく場合に、力積(力×時間)の分だけ、全体の運動量が増加するのは、ニュートン力学の世界と同じです(『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の第 6 節(特に、(6.11)式など)を参照)。したがって、運動量の保存則は、$\omega_{\mathrm{R}}(\tau) = 0$ も用いて、

\[ \left[ \frac{m(\tau) \cdot 0}{\sqrt{1 - ( 0/c )^{2}}} \right] + f \cdot \varDelta \tau = \left[ \frac{m(\tau + \varDelta \tau) \cdot \omega_{\mathrm{R}}(\tau + \varDelta \tau)}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \cdot ( -\omega_{\mathrm{P}} )}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right] \tag{4.8} \]

となります。ところで、右辺第 1 項に、(1.2)式を適用すると、$\omega_{\mathrm{R}}(\tau) = 0$ も用いて、

\[ \frac{m(\tau + \varDelta \tau) \cdot \omega_{\mathrm{R}}(\tau + \varDelta \tau)}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} = m(\tau) \, \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \cdot \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.9} \]

となるコトが分かります(ちょっとメンドくさいですけど、計算自体は簡単です)。この(4.9)式、および、(4.7)式を、(4.8)式に代入して、整理すると、

\[ f \, \varDelta \tau = m(\tau) \, \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \, \varDelta \tau - \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta \tau + \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{D}} \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.10} \]

という関係式が得られます。

■(4.10)式には、未知数が 2 つ($d\omega_{\mathrm{R}}(\tau)/d\tau$ と $\varDelta m_{\mathrm{D} \,}$)あります。よって、方程式が 1 つだけでは、完全に解くコトができません。ココで、「エネルギーの保存則」が必要になります。ただし、特殊相対論の世界では、ロケットに乗っている観測者から見たエネルギーが、(3.7)式で与えられます。また、外力 $f$ がはたらく場合は、外力がした仕事(力×移動距離=力×速度×移動時間)の分だけ、全体のエネルギーが増加します(『任意の方向へ運動する慣性系に対するローレンツ変換、そして、特殊相対論的に正しい運動方程式へ。』の第 6 節(特に、(6.38)式など)を参照)。したがって、この場合のエネルギーの保存則は、$\omega_{\mathrm{R}}(\tau) = 0$ も用いて、

\[ \left[ \frac{m(\tau) \, c^{2}}{\sqrt{1 - ( 0/c )^{2}}} \right] + f \cdot 0 \cdot \varDelta \tau = \left[ \frac{m(\tau + \varDelta \tau) \, c^{2}}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \, c^{2}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right] \tag{4.11} \]

となります。ところで、右辺第 1 項に、(1.2)式を適用すると、$\omega_{\mathrm{R}}(\tau) = 0$、および、(4.1)式も用いて、

\[ \frac{m(\tau + \varDelta \tau) \, c^{2}}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} = m(\tau) \, c^{2} - \mu \, c^{2} \cdot \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.12} \]

となるコトが分かります(ちょっとメンドくさいですけど、計算自体は簡単です)。(4.12)式を(4.11)式に代入し、(4.7)式も用いて、整理すると、

\[ 0 = -\mu \, \varDelta \tau + \frac{\mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta \tau - \frac{1}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{D}} \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.13} \]

という関係式が得られます。

■運動量の保存則=(4.10)式と、エネルギーの保存則=(4.13)式から、$\varDelta m_{\mathrm{D}}$ を消去すると、ただちに、

\[ f + \omega_{\mathrm{P}} \mu = m(\tau) \, \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \tag{4.14} \]

という関係式が得られます。コレは、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、ロケットに乗っている観測者が感じる加速度 $d\omega_{\mathrm{R}}(\tau)/d\tau$(および、その加速度に比例する慣性力)を示す、とても重要な方程式です。この(4.14)式が、「ニュートン力学的に正しいツィオルコフスキーの公式」=(2.6)式と、まったく同じカタチをしているコトは、(実は、当然のコトなのですけど)注目に値します。また、(2.7)式と同じ議論によって、この場合にも、ロケットの推力 $F_{\mathrm{T}}$ が、$F_{\mathrm{T}} = \omega_{\mathrm{P}} \mu$ で表されるコトが、分かります。


□ちなみに、(4.14)式は、$\omega_{\mathrm{R}}(\tau)$ に関する微分方程式とみなすコトもできます。この微分方程式を解くコトで得られる $\omega_{\mathrm{R}}(\tau)$ は、「基準となるある時刻の時点におけるロケットに乗っている観測者から見た、時刻 $\tau$ の時点における(計算上の)ロケットの速さ」を表します。実は、この $\omega_{\mathrm{R}}(\tau)$ は、光速 $c$ を超えるコトができます。でも、「ロケットに乗っている観測者が測った距離」は、定義ができません(第 3 節を参照)。そのため、ロケットに乗っている観測者は、基準となるある時刻から、時刻 $\tau$ までのあいだに、ロケットがどれだけの距離を進んだのか――すなわち、$\omega_{\mathrm{R}}(\tau)$ の値を、実際に測定する手段がありません。ロケットに乗っている観測者は、「これまでに、十分な大きさの加速度(慣性力)を感じ続けてきたから、このロケットは、光速を超えているハズだ。直接、測ったワケではないけど」と、思うに違いありません。


■(4.14)式によって、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、ロケットに乗っている観測者の感じる加速度 $d\omega_{\mathrm{R}}(\tau)/d\tau$ を、(あらゆる時刻 $\tau$ において)知るコトができます。でも、この加速度を積分して、速さ $\omega_{\mathrm{R}}(\tau)$ を計算しても 、物理的な意味は、あまり、ありません。ところが、ココで、ロケットに乗っている観測者から見た慣性系の速さ(慣性系にある時計がすれ違う速さ)$-u_{\mathrm{R}}$ と、ロケットに乗っている観測者の感じる加速度 $d\omega_{\mathrm{R}}(\tau)/d\tau$ のあいだに、(3.4)式で与えられる関係があるコトを、思い出しましょう。(4.14)式に、(3.4)式を代入すると、ただちに、

\[ f + \omega_{\mathrm{P}} \mu = \frac{m(\tau)}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{4.15} \]

という微分方程式が得られます。この(4.15)式こそ、「特殊相対論的に正しいツィオルコフスキーの公式」と呼ぶにふさわしい微分方程式です。この(4.15)式を解くコトで、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$ が、得られるからです。


□この議論については、次のように考えてみると、分かりやすいかも知れません。いま、窓のないロケットのなかに、観測者が閉じ込められているトコロを、イメージしてみて下さい。ロケットが加速すると、この観測者には、慣性力がはたらきます。そのコトで、窓のないロケットのなかに閉じ込められていても、この観測者は、自分が加速度運動をしているコトが、分かります。加速度の大きさと時間から、いまのロケットの速さを、計算するコトもできるでしょう。ただし、窓のないロケットのなかに閉じ込められているので、それを検証する手段は、ありません。ところが、このロケットに、窓が見つかったとします。すると、窓の外の景色(の動き)を観測するコトで、ロケットの速さが分かります。もちろん、特殊相対論の世界では、体感した加速度から求めた計算上の速さと、窓の外の景色の動きから分かる速さが、等しくはなりません。でも、両者は、まったくの無関係というワケではありません。この 2 つの異なる「速さ」の関係を与える公式こそが、(3.4)式であり、 2 つの微分方程式=(4.14)式と(4.15)式をつなぐカギなのです。


■ところで、いま、時刻 $\tau = 0$ のときの、推進剤を含んだロケット全体の質量 $m(\tau)$ を、$m(0) \equiv m_{0}$ と表すコトにすると、(4.1)式で定義した推進剤の消費速度 $\mu \equiv -dm(\tau)/d\tau$($\, \mu$ は定数)も用いて、$m(\tau)$ は、

\[ m(\tau) = m_{0} - \mu \tau \tag{4.16} \]

と表されるコトは、明らかです。したがって、(4.15)式は、

\[ f + \omega_{\mathrm{P}} \mu = \frac{m_{0} - \mu \tau}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{4.17} \]

と書くコトもできます。外力 $f$ が定数である場合には、この微分方程式を解くのは、とても簡単です。でも、それは、第 6 節で、じっくりと楽しむコトにして、いまは、先に進みましょう。


□ちなみに、運動量の保存則=(4.10)式と、エネルギーの保存則=(4.13)式から、$\varDelta m_{\mathrm{D}}$ ではなく、$d\omega_{\mathrm{R}}(\tau)/d\tau$ のほうを消去すると、(4.7)式も用いて、ただちに、
\[ \qquad \varDelta m_{\mathrm{D}} = \left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, \mu \, \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.18a} \] \[ \qquad \varDelta m_{\mathrm{P}} = \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \mu \, \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{4.18b} \]
という関係式が得られます。(4.18a)〜(4.18b)式より、ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ が大きくなるほど、時間 $\varDelta \tau$ のあいだに消費される推進剤の質量 $\mu \varDelta \tau$ のうち、実際に噴射される推進剤の質量 $\varDelta m_{\mathrm{P}}$ の割合が小さくなるコトが分かります。また、ロケットに乗っている観測者から見た推進剤の噴射速度を、光速に等しくするコト($\omega_{\mathrm{P}} = c$)は不可能であるコトも、分かります。ただし、コレは、あくまでも、質量=運動エネルギーに由来するエネルギーを用いて、質量を持つ推進剤を噴射するロケットについての計算結果であるコトに、注意が必要です。たとえば、「光子ロケット」のように、質量を持たない電磁波を発振する反動で推進する場合(コレも、「ロケット」の定義を満たすコトに注意)は、もちろん、推進剤の噴射速度が、光速になります。ココでは扱いませんけれど、このようなロケットの運動の計算も、とても興味深い結果が得られるので、ぜひ、挑戦してみて下さいね。


 ■慣性系にいる観測者から見た場合の、ツィオルコフスキーの公式の導出 ――  

■ココまでは、ロケットに乗っている観測者から見たロケットの運動を、計算してきましたけど、まったく同じロケットの運動を、慣性系にいる観測者が見た場合は、どうなるでしょうか。その場合には、次の図のように、時刻 $t$ のときの、推進剤を含んだロケット全体の質量が $m(t)$、ロケットの速さが $u_{\mathrm{R}}(t)$ となります。時間が $\varDelta t$ だけ経過すれば、ロケットの質量と速さは、それぞれ、$m(t+\varDelta t)$、$u_{\mathrm{R}}(t+\varDelta t)$ になります。ココで、ロケットは、時間 $\varDelta t$ が経過するあいだに、質量 $\varDelta m_{\mathrm{P}}$ の推進剤を、速さ $u_{\mathrm{P}}$ で噴射します(あとで計算しやすいよう、慣性系にいる観測者から見て、右向きのとき $u_{\mathrm{P}} > 0$、左向きのとき $u_{\mathrm{P}} < 0$ となるように、$u_{\mathrm{P}}$ を定めます)。また、このときに生じる質量欠損を、$\varDelta m_{\mathrm{D}}$ と書きましょう。最後に、ロケット本体には、ロケット自身の推力のほかに、外力 $f$ がはたらくコトが、想定されています。ちなみに、いまは、直線上(1 次元)の運動だけを考えているので、(3.9)式で説明したように、ロケットに乗っている観測者にとっても、慣性系にいる観測者にとっても、外力 $f$ の値は、同じになります。以上の設定を図にすると、次のようになるハズです。
図7:「特殊相対論的に正しいツィオルコフスキーの公式」の導出:化学ロケット編(2)(約45KB)

■さて、この場合、推進剤の噴射前後での、全体の質量の関係は、質量欠損 $\varDelta m_{\mathrm{D}}$ があるコトを考慮して、

\[ m(t + \varDelta t) = m(t) - \varDelta m_{\mathrm{P}} - \varDelta m_{\mathrm{D}} \tag{4.19} \]

となります。この(4.19)式に、(1.2)式を適用すると、ただちに、

\[ \varDelta m_{\mathrm{P}} = -\frac{dm(t)}{dt} \, \varDelta t - \varDelta m_{\mathrm{D}} \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{4.20} \]

となるコトが分かります。(4.20)式に現れる $-dm(t)/dt$ は、慣性系にいる観測者から見た推進剤の消費速度を表します。慣性系にいる観測者から見た推進剤の消費速度 $-dm(t)/dt$ と、ロケットに乗っている観測者から見た推進剤の消費速度 $\mu = -dm(\tau)/d\tau$(ココでは、$\mu$ は定数と仮定していたコトに注意)の関係は、(4.2)式で与えられているので、(4.20)式は、

\[ \varDelta m_{\mathrm{P}} = \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, \mu \, \varDelta t - \varDelta m_{\mathrm{D}} \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{4.21} \]

と表すコトもできます。

■必要な準備ができたので、早速、いま、考えた設定に対して、「運動量の保存則」を、適用してみましょう。この場合の、運動量の保存則は、

\[ \left[ \frac{m(t) \cdot u_{\mathrm{R}}(t)}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \right] + f \cdot \varDelta t = \left[ \frac{m(t + \varDelta t) \cdot u_{\mathrm{R}}(t + \varDelta t)}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \cdot u_{\mathrm{P}}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} \right] \tag{4.22} \]

となります。この(4.22)式の右辺第 1 項に、(1.2)式を適用すると、

\[ \begin{align} \frac{m(t + \varDelta t) \cdot u_{\mathrm{R}}(t + \varDelta t)}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} &= \frac{m(t) \cdot u_{\mathrm{R}}(t)}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, u_{\mathrm{R}}(t) \, \frac{dm(t)}{dt} \cdot \varDelta t \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \cdot \varDelta t \nonumber \\[0em] &\hspace{20em} \mbox{($\varDelta t \, \to \, 0$)} \tag{4.23} \end{align} \]

となるコトが分かります(ちょっとメンドくさいですけど、計算自体は簡単です)。(4.22)式の右辺第 2 項は、「速度の合成則」=(3.5)式を用いて、慣性系にいる観測者から見た推進剤の噴射速度 $u_{\mathrm{P}}$ を、ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$(ココでは、$\omega_{\mathrm{P}}$ は定数と仮定していたコトに注意)に変換できるので、(ロケットの速さ $u_{\mathrm{R}}(t)$ によって値が変化する $u_{\mathrm{P}}$ ではなく)定数である $\omega_{\mathrm{P}}$ を使って、

\[ \frac{\varDelta m_{\mathrm{P}} \cdot u_{\mathrm{P}}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{u_{\mathrm{R}}(t) - \omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{P}} \tag{4.24} \]

と表せるコトが分かります(かなりメンドくさいですけど、計算自体は簡単です。ただし、$u_{\mathrm{P}}$ と $\omega_{\mathrm{P}}$ の向きの決め方に、注意して下さい)。(4.23)〜(4.24)式を、(4.22)式に代入して、整理すると、

\[ \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, f \, \varDelta t = u_{\mathrm{R}}(t) \, \frac{dm(t)}{dt} \, \varDelta t + \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \, \varDelta t + \frac{u_{\mathrm{R}}(t) - \omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{P}} \] \[ \hspace{31em} \mbox{($\varDelta t \, \to \, 0$)} \tag{4.25} \]

という関係式が得られます。

■次に、「エネルギーの保存則」を、適用してみましょう。この場合の、エネルギーの保存則は、

\[ \left[ \frac{m(t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \right] + f \cdot u_{\mathrm{R}}(t) \cdot \varDelta t = \left[ \frac{m(t + \varDelta t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \, c^{2}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} \right] \tag{4.26} \]

となります。先ほどと同じように、この(4.26)式の右辺第 1 項に、(1.2)式を適用すると、

\[ \begin{align} \frac{m(t + \varDelta t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} &= \frac{m(t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, c^{2} \, \frac{dm(t)}{dt} \cdot \varDelta t \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{m(t) \, u_{\mathrm{R}}(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \cdot \varDelta t \nonumber \\[0em] &\hspace{20em} \mbox{($\varDelta t \, \to \, 0$)} \tag{4.27} \end{align} \]

となるコトが分かります(ちょっとメンドくさいですけど、計算自体は簡単です)。(4.26)式の右辺第 2 項も、先ほどと同じように、ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$(=定数)を用いて、

\[ \frac{\varDelta m_{\mathrm{P}} \, c^{2}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{c^{2} - \omega_{\mathrm{P}} \, u_{\mathrm{R}}(t)}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{P}} \tag{4.28} \]

と表せるコトが分かります(かなりメンドくさいですけど、計算自体は簡単です。ただし、$u_{\mathrm{P}}$ と $\omega_{\mathrm{P}}$ の向きの決め方に、注意して下さい)。(4.27)〜(4.28)式を、(4.26)式に代入して、整理すると、

\[ \begin{align} \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, f \, u_{\mathrm{R}}(t) \, \varDelta t &= c^{2} \, \frac{dm(t)}{dt} \, \varDelta t + \frac{m(t) \, u_{\mathrm{R}}(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \, \varDelta t \nonumber \\[0.5em] &\qquad + \frac{c^{2} - \omega_{\mathrm{P}} \, u_{\mathrm{R}}(t)}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \varDelta m_{\mathrm{P}} \hspace{4em} \mbox{($\varDelta t \, \to \, 0$)} \tag{4.29} \end{align} \]

という関係式が得られます。

■運動量の保存則=(4.25)式と、エネルギーの保存則=(4.29)式から、$\varDelta m_{\mathrm{P}}$ を消去すると、ただちに、

\[ f - \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{dm(t)}{dt} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{4.30} \]

という微分方程式が得られます。この(4.30)式も、「特殊相対論的に正しいツィオルコフスキーの公式」と呼ぶべき、もう 1 つの微分方程式です。(4.30)式を解けば、慣性系にある時計の時刻が $t$ のときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(t)$ が得られるからです。


□実は、一般に、(4.30)式を解くのは、極めて難しいと思われます。ただし、外力 $f$ がゼロの場合に、推進剤を含んだロケット全体の質量 $m$ の値が与えられたときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(m)$ を求めるコトならば、とても簡単です。この計算は、第 6 節で、詳しく検討しましょう。

□ちなみに、運動量の保存則=(4.25)式と、エネルギーの保存則=(4.29)式から、$\varDelta m_{\mathrm{P}}$ ではなく、$du_{\mathrm{R}}(t)/dt$ のほうを消去すると、(4.2)式、および、(4.21)式も用いて、ただちに、
\[ \begin{align} \qquad \varDelta m_{\mathrm{D}} &= -\left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, \frac{dm(t)}{dt} \, \varDelta t \nonumber \\[0.5em] &= \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, \left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, \mu \, \varDelta t \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{4.31a} \end{align} \] \[ \begin{align} &\qquad \varDelta m_{\mathrm{P}} = -\sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \frac{dm(t)}{dt} \, \varDelta t = \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \mu \, \varDelta t \nonumber \\[0em] &\hspace{31em} \mbox{($\varDelta t \, \to \, 0$)} \tag{4.31b} \end{align} \]
という関係式が得られます。ただし、一般に、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$ を、慣性系にある時計の時刻 $t$ を用いて表す――すなわち、$u_{\mathrm{R}}(t)$ の具体的なカタチを求めるコトは、高校までに習う知識だけでは、極めて難しいと思われます。このコトについても、第 6 節で、詳しく検討します。

□なお、慣性系にいる観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(4.30)式は、慣性系にいる観測者から見た物理量を、ロケットに乗っている観測者から見た物理量に変換する関係式((4.2)式、および、(3.4)式)を用いると、ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(4.15)式に、変換できます。観測者が違うだけで、同じロケットの運動を観測しているので、当然と言えば当然です。


■さて、ココまでの議論で、質量=運動エネルギー(質量欠損)に由来するエネルギーを用いて、推進剤を噴射するロケットについて、「特殊相対論的に正しいツィオルコフスキーの公式」を、導くコトができました。ココまでの議論で得られた微分方程式を、実際に解いて、ロケットの運動の様子を具体的に検討するのは、第 6 節のお楽しみとしましょう。次の第 5 節では、推進剤を噴射するためのエネルギーが、外部から供給される場合について、「特殊相対論的に正しいツィオルコフスキーの公式」を導いてみます。


□ところで、この第 4 節での議論には、大きな問題が、1 つ、あります。それは、「質量欠損に由来する質量=運動エネルギーが、すべて、推進剤の噴射に使われる」――と、(暗黙のうちに)仮定しているコトです。日本が誇る LE-7A エンジンでさえも、質量欠損に由来する質量=運動エネルギーのうち、実際に、推進剤の噴射に使えるのは $70 \, \%$ に満たない――というコトを考えると、エネルギーロスがまったくない、効率 $100 \, \%$ のエンジンを仮定するのは、非現実的と言わなければいけないでしょう。エネルギー変換効率を考慮に入れた計算は、ココでは割愛しますけど、興味があれば、ぜひ、挑戦してみて下さいね。

■このページの一番上に戻る。



 ■5. 特殊相対論的に正しいツィオルコフスキーの公式(2) ―― 電気推進ロケット編

■前節(第 4 節)では、推進剤そのものだけでなく、推進剤を噴射するために必要なエネルギーについても、あらかじめ、ロケットに搭載されている――という設定で、検討を進めました。つまり、完全に、ロケット単独で推進するコトが、想定されていました。ココ(第 5 節)では、推進剤を噴射するのに必要なエネルギーが、外部から供給される場合を、考えてみましょう。もちろん、ロケットなので、推進剤そのものは、あらかじめ、ロケットに搭載されていなければなりません。でも、推進剤を噴射するためのエネルギーは、航行しながら調達する――というワケです。よって、外部からのエネルギー供給がなくなると、このタイプのロケットは、推進剤を噴射できません。すなわち、ロケット単独では、推力を得るコトができない――というコトです。たとえば、小惑星探査機「はやぶさ 2」のイオンエンジン $\mu 10$ は、太陽電池パネルで発生する電気エネルギーを使って、推進剤(キセノンイオン)を噴射します。太陽電池パネルに、十分なエネルギーを持つ光が当たらないと(外部からのエネルギー供給がないと)、推進剤を噴射するコトができません。そういう視点で見ると、小惑星探査機「はやぶさ 2」のイオンエンジン $\mu 10$ は、このタイプのロケットに、分類するコトができるでしょう。そんなワケで、ココ(第 5 節)では、そのようなタイプのロケット――推進剤を噴射するのに必要なエネルギーが、外部から供給される場合の、「特殊相対論的に正しいツィオルコフスキーの公式」を、導いてみましょう。

■ところで、ロケットに、外部からエネルギーを供給する方法は、いろいろと考えられます。けれども、実際に試してみると、ロケットの運動を、簡単に計算できるようなエネルギー供給方法は、ほとんど、ないようです。もちろん、「極めて難しいけど、解こうと思えば解ける」というモノも、ないわけではありません。そもそも、どんな方法も、数値計算によって、ロケットの運動を調べるコトはできます。ただ、ココでは、イチバン、計算が簡単なエネルギー供給方法を、考えるコトにしましょう。それは、次のような方法です。

■まず、ココでも、ロケットに乗っている観測者から見た推進剤の消費速度 $\mu = -dm(\tau)/d\tau$、および、推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、一定とします。のちの計算からも明らかになるように、コレを実現するためには、ロケットに乗っている観測者から見て、単位時間当たり一定のエネルギーが、外部から供給され続ける必要があります。また、一般に、外部からエネルギーを受け取るときには、同時に、運動量も受け取ります。たとえば、太陽電池パネルに光を当てるコトで、外部からエネルギーを受け取るとき、同時に、光圧に由来する力積がはたらくコトで、運動量も受け取るコトになります。でも、簡単のため、ココでは、ロケットに乗っている観測者から見て、外部から受け取る運動量は、常にゼロであると仮定しましょう。コレを実現するためには、たとえば、ロケットの側面全体に太陽電池パネルを配置するコトが考えられるかも知れません。反対側の太陽電池パネルが受け取る光圧による力積が、お互いに打ち消し合い、全体では、外部から受け取る運動量は、常にゼロになるというワケです。または、余分な運動量が打ち消されるように、受け取ったエネルギー=運動量の一部を、放出するという方法もあります。現実的ではありませんけど、とりあえず、叩き台として、この設定で、計算してみましょう。

 ■ロケットに乗っている観測者から見た場合の、ツィオルコフスキーの公式の導出 ――  

■早速、上で説明した設定で、ロケットに乗っている観測者から見たロケットの運動の様子を、描いてみます。時刻 $\tau$ のときの、推進剤を含んだロケット全体の質量を $m(\tau)$、ロケット本体の速さを $\omega_{\mathrm{R}}(\tau)$ と表しましょう。ココでも、$\omega_{\mathrm{R}}(\tau) = 0$ です。時間が $\varDelta \tau$ だけ経過したとき、ロケットの質量と速さは、それぞれ、$m(\tau+\varDelta \tau)$、$\omega_{\mathrm{R}}(\tau+\varDelta \tau)$ になります。ロケットは、時間 $\varDelta \tau$ が経過するあいだに、外部からエネルギー $\varepsilon \varDelta \tau$ を受け取り、質量 $\varDelta m_{\mathrm{P}}$ の推進剤を、速さ $-\omega_{\mathrm{P}}$ で噴射します(向きに注意)。なお、$\varepsilon$ は、「ロケットに乗っている観測者から見た、単位時間当たりに外部から供給されるエネルギー」を表します。すでに説明したとおり、ロケットに乗っている観測者から見て、外部から受け取る運動量は、完全に打ち消し合って、常にゼロになる(そのように設計・調整する)と仮定されています。最後に、ロケット本体には、ロケット自身の推力のほかに、外力 $f$ がはたらくコトを、忘れずに、想定しておきます。以上の設定を図にすると、次のようになるハズです。
図8:「特殊相対論的に正しいツィオルコフスキーの公式」の導出:電気推進ロケット編(1)(約43KB)

■さて、今回の設定では、質量欠損は生じません。したがって、推進剤の噴射前と噴射後における、全体の質量の関係から、(1.2)式と、推進剤の消費速度 $\mu$ の定義($\mu \equiv -dm(\tau)/d\tau$)を用いて、ただちに、

\[ \begin{align} m(\tau) &= m(\tau + \varDelta \tau) + \varDelta m_{\mathrm{P}} \nonumber \\[0.5em] &= m(\tau) + \frac{dm(\tau)}{d\tau} \cdot \varDelta \tau + \varDelta m_{\mathrm{P}} = m(\tau) - \mu \varDelta \tau + \varDelta m_{\mathrm{P}} \quad \mbox{($\varDelta \tau \, \to \, 0$)} \nonumber \\[0.7em] &\hspace{-1em} \mbox{したがって} \quad \varDelta m_{\mathrm{P}} = \mu \varDelta \tau \quad \mbox{($\varDelta \tau \, \to \, 0$)} \tag{5.1} \end{align} \]

となるコトが分かります。

■それでは、ココで、「運動量の保存則」を、適用してみましょう。この場合の、運動量の保存則は、外部から供給される運動量がゼロであるコトに注意して、

\[ \left[ \frac{m(\tau) \cdot 0}{\sqrt{1 - ( 0/c )^{2}}} \right] + f \cdot \varDelta \tau + 0 = \left[ \frac{m(\tau + \varDelta \tau) \cdot \omega_{\mathrm{R}}(\tau + \varDelta \tau)}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \cdot (-\omega_{\mathrm{P}})}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right] \tag{5.2} \]

となります。第 4 節と同じ要領で、(1.2)式を適用して計算すると、(5.1)式も用いて、さらに、その結果に、$d\omega_{\mathrm{R}}(\tau)/d\tau$ と $du_{\mathrm{R}}(\tau)/d\tau$ のあいだの関係を表す(3.4)式を適用するコトで、最終的に、

\[ f + \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} = m(\tau) \, \frac{d\omega_{\mathrm{R}}(\tau)}{d\tau} \quad = \frac{m(\tau)}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{5.3} \]

という関係式が得られます。さて、いま、時刻 $\tau = 0$ のときの、推進剤を含んだロケット全体の質量 $m(\tau)$ を、$m(0) \equiv m_{0}$ と表すコトにすると、$m(\tau)$ は、推進剤の消費速度 $\mu$ を用いて、$m(\tau) = m_{0}-\mu \tau$ と表されるコトは明らかです((4.16)式を参照)。したがって、(5.3)式は、ただちに、

\[ f + \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} = \frac{m_{0} - \mu \tau}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{5.4} \]

と書き直されます。(5.4)式が、この場合の「特殊相対論的に正しいツィオルコフスキーの公式」と呼ぶべき微分方程式です。この(5.4)式を解くコトで、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$ が、得られるからです。外力 $f$ が定数である場合には、この微分方程式を解くのは、とても簡単です(第 6 節で、実際に、解いてみます)。なお、外部からのエネルギー供給がない場合(推進剤を噴射するためのエネルギーが、推進剤自身の質量欠損に由来する場合)の「特殊相対論的に正しいツィオルコフスキーの公式」=(4.17)式と、ココで求めた(5.4)式が、似たカタチをしているコトは、注目に値します。でも、いまは、そのコトを指摘するにとどめて、先へ進みましょう。

■次に、「エネルギーの保存則」を、適用してみましょう。この場合の、エネルギーの保存則は、微小時間 $\varDelta \tau$ のあいだに、外部から供給されるエネルギーが $\varepsilon \varDelta \tau$ であるコトに注意して、

\[ \left[ \frac{m(\tau) \, c^{2}}{\sqrt{1 - ( 0/c )^{2}}} \right] + f \cdot 0 \cdot \varDelta \tau + \varepsilon \varDelta \tau = \left[ \frac{m(\tau + \varDelta \tau) \, c^{2}}{\sqrt{1 - ( \omega_{\mathrm{R}}(\tau + \varDelta \tau)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \, c^{2}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right] \tag{5.5} \]

となります。先ほどと同じように計算すると、最終的に、

\[ \varepsilon = \left( \frac{1}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} - 1 \right) \, \mu c^{2} \tag{5.6} \]

という関係式が得られます。エネルギーが外部から供給される場合、単位時間当たりに消費できるエネルギーが制限されるので、推進剤の消費速度 $\mu$ と噴射速度 $\omega_{\mathrm{P}}$ は、この関係式によって規定されます。なお、(5.6)式からは、$\mu$ と $\omega_{\mathrm{P}}$ をともに一定にするためには、ロケットに乗っている観測者から見た、単位時間当たりに外部から供給されるエネルギー $\varepsilon$ が、一定でなければいけないコトが分かります。


□実は、この場合には、慣性系にいる観測者から見た、単位時間当たりに外部から供給されるエネルギー $e$ も、一定になるコトが、証明できます。いま、慣性系にいる観測者から見て、単位時間当たりにロケットに供給されるエネルギーを $e$、単位時間当たりにロケットに供給される運動量を $p$ で表します。このとき、慣性系にいる観測者から見ると、時間 $\varDelta t$ のあいだに、エネルギー $e \varDelta t$、運動量 $p \varDelta t$ が、ロケットに供給されるコトになります。他方で、ロケットに乗っている観測者から見たときには、時間 $\varDelta \tau$ のあいだに、ロケットに供給されるエネルギーは $\varepsilon \varDelta \tau$、運動量は $0$ でしたので、(3.8a)式により、
\[ \begin{align} \qquad e \varDelta t &= \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \bigl( \varepsilon \varDelta \tau + u_{\mathrm{R}}(t) \cdot 0 \bigr) = \frac{\varepsilon \varDelta \tau}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} = \varepsilon \varDelta t \tag{5.7a} \\[0em] &\hspace{26em} \mbox{($\varDelta t, \ \varDelta \tau \, \to \, 0$)} \nonumber \end{align} \] \[ \begin{align} \qquad p \varDelta t &= \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \left( 0 + \frac{u_{\mathrm{R}}(t)}{c^{2}} \, \varepsilon \varDelta \tau \right) = \frac{(u_{\mathrm{R}}(t)/c^{2}) \, \varepsilon \varDelta \tau}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} = \frac{u_{\mathrm{R}}(t)}{c^{2}} \, \varepsilon \varDelta t \tag{5.7b} \\[0em] &\hspace{26em} \mbox{($\varDelta t, \ \varDelta \tau \, \to \, 0$)} \nonumber \end{align} \]
となります(ココでは、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、慣性系にある時計の時刻 $t$ のあいだに成り立つ関係を与える(3.2)式も使いました)。(5.7a)式より、この場合に、$e = \varepsilon$ となるコトが分かります。強調しておきたいのですけど、$e = \varepsilon$ となるのは、あくまでも、ココでの設定に限ったハナシですので、注意が必要です(もちろん、一般には、$e = \varepsilon$ という関係は、成り立ちません)。

□なお、(5.7b)式では、ロケットに乗っている観測者から見た、単位時間当たりに供給される運動量がゼロなのに、慣性系にいる観測者から見た、単位時間当たりに供給される運動量 $p$ はゼロにならないコトが、示されています。コレには、「光行差現象」が影響しています。雨が降っている状況を、想像してみて下さい。立ち止まっているときには、雨粒は、上から落ちてきます。でも、走っているときには、雨粒は、前方からぶつかってきます。このように、当事者の運動状態によって、雨粒のおよぼす影響が変わります。コレが、光行差現象(正しくは、光行差現象がロケットへのエネルギー供給に与える影響の喩え話)です。実際、ココで検討したエネルギー供給方法は、具体的には、たとえば、次の図のようにして実現できます。この図から、(5.7a)〜(5.7b)式の成り立つ状況が、理解できるでしょう。
図9:「特殊相対論的に正しいツィオルコフスキーの公式」:光行差現象(約68KB)

□ちなみに、光=電磁波のカタチで、エネルギーを供給する場合、(3.8a)式を使わずに、マクスウェルの方程式から導かれる、特殊相対論的に正しい光行差現象の公式だけから、(5.7a)〜(5.7b)式を導くコトができます。興味があれば、ぜひ、挑戦してみて下さいね。


 ■慣性系にいる観測者から見た場合の、ツィオルコフスキーの公式の導出 ――  

■次に、まったく同じロケットの運動を、慣性系にいる観測者が見た場合について、計算してみましょう。その場合には、次の図のように、時刻 $t$ のときの、推進剤を含んだロケット全体の質量が $m(t)$、ロケットの速さが $u_{\mathrm{R}}(t)$ となります。時間が $\varDelta t$ だけ経過すれば、ロケットの質量と速さは、それぞれ、$m(t+\varDelta t)$、$u_{\mathrm{R}}(t+\varDelta t)$ になります。ココで、ロケットは、時間 $\varDelta t$ が経過するあいだに、外部からエネルギー $e \varDelta t$ と、運動量 $p \varDelta t$ を受け取り、質量 $\varDelta m_{\mathrm{P}}$ の推進剤を、速さ $u_{\mathrm{P}}$ で噴射します(あとで計算しやすいよう、慣性系にいる観測者から見て、右向きのとき $u_{\mathrm{P}} > 0$、左向きのとき $u_{\mathrm{P}} < 0$ となるように、$u_{\mathrm{P}}$ を定めます。また、$e$ は「慣性系にいる観測者から見た、単位時間当たりに外部から供給されるエネルギー」を、$p$ は「慣性系にいる観測者から見た、単位時間当たりに外部から供給される運動量」を表します)。そして、最後に、ロケットには、ロケット自身の推力のほかに、外力 $f$ がはたらくと、想定されています。ちなみに、いまは、直線上(1 次元)の運動だけを考えているので、(3.9)式で説明したように、ロケットに乗っている観測者にとっても、慣性系にいる観測者にとっても、外力 $f$ の値は、同じになります。以上の設定を図にすると、次のようになるハズです。
図10:「特殊相対論的に正しいツィオルコフスキーの公式」の導出:電気推進ロケット編(2)(約47KB)

■さて、推進剤の噴射前と噴射後における、全体の質量の関係から、(1.2)式と(4.2)式を用いて、

\[ m(t) = m(t + \varDelta t) + \varDelta m_{\mathrm{P}} \quad = m(t) + \frac{dm(t)}{dt} \cdot \varDelta t + \varDelta m_{\mathrm{P}} \quad \mbox{($\varDelta t \, \to \, 0$)} \] \[ \quad \mbox{したがって} \quad \varDelta m_{\mathrm{P}} = -\frac{dm(t)}{dt} \, \varDelta t = \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, \mu \, \varDelta t \quad \mbox{($\varDelta t \, \to \, 0$)} \tag{5.8} \]

となるコトが分かります(いま、$\mu$ は定数と仮定していたコトに注意)。

■ココで、「運動量の保存則」を、適用してみましょう。この場合の、運動量の保存則は、考える微小時間 $\varDelta t$ のあいだに、外部から $p \varDelta t$ の運動量が供給されるコトに注意して、

\[ \left[ \frac{m(t) \cdot u_{\mathrm{R}}(t)}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \right] + f \cdot \varDelta t + p \varDelta t = \left[ \frac{m(t + \varDelta t) \cdot u_{\mathrm{R}}(t + \varDelta t)}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \cdot u_{\mathrm{P}}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} \right] \tag{5.9} \]

となります。この式の右辺は、(4.23)〜(4.24)式で、すでに、与えられています。よって、(5.8)式も用いて、整理すると、最終的に、

\[ \begin{align} f + p &= \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \left( u_{\mathrm{R}}(t) - \frac{u_{\mathrm{R}}(t) - \omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right) \, \frac{dm(t)}{dt} \tag{5.10} \end{align} \]

という関係式が得られます。

■引き続いて、「エネルギーの保存則」を、適用してみましょう。この場合の、エネルギーの保存則は、考える微小時間 $\varDelta t$ のあいだに、外部から $e \varDelta t$ のエネルギーが供給されるコトに注意して、

\[ \left[ \frac{m(t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \right] + f \cdot u_{\mathrm{R}}(t) \cdot \varDelta t + e \varDelta t = \left[ \frac{m(t + \varDelta t) \, c^{2}}{\sqrt{1 - ( u_{\mathrm{R}}(t + \varDelta t)/c )^{2}}} \right] + \left[ \frac{\varDelta m_{\mathrm{P}} \, c^{2}}{\sqrt{1 - ( u_{\mathrm{P}}/c )^{2}}} \right] \] \[ {} \tag{5.11} \]

となります。この式の右辺は、(4.27)〜(4.28)式で、すでに、与えられています。よって、(5.8)式も用いて、整理すると、最終的に、

\[ \begin{align} f \cdot u_{\mathrm{R}}(t) + e &= \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{ m(t) \, u_{\mathrm{R}}(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \nonumber \\[0.5em] &\qquad + \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \left( c^{2} - \frac{c^{2} - \omega_{\mathrm{P}} \, u_{\mathrm{R}}(t)}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \right) \, \frac{dm(t)}{dt} \tag{5.12} \end{align} \]

という関係式が得られます。

■ココで、$p$ と $e$ のあいだに成り立つ関係について、考えてみましょう。(3.6)式より、一般に、特殊相対論の世界では、慣性系にいる観測者から見た、任意の質点のエネルギー $K_{k}$ と運動量 $p_{k}$ のあいだに、

\[ p_{k} = \frac{u}{c^{2}} \, K_{k} \qquad \mbox{($u$ は、慣性系にいる観測者から見た、考える質点の速さ)} \tag{5.13} \]

という関係が、成り立っているコトが分かります。ココで想定している、ロケットに供給される外部エネルギー(エネルギー塊)に当てはめると、$p_{k} = p \varDelta t, \ K_{k} = e \varDelta t$ となります。また、$u$ については、次のように考えるコトができます。まず、1 次元の運動だけを仮定しているので、エネルギー塊の速さ $u$ については、ロケットの進行方向に対して平行な成分だけを考えれば、十分です。ところで、ロケットに乗っている観測者から見て、外部から受け取る運動量は、ゼロであると仮定していました。コレを実現するには、エネルギー塊は、ロケットに対して、相対速度がゼロでなければいけません。空中給油機で、飛行中の航空機に燃料を供給するトコロを、想像してみて下さい。航空機の前、または、後ろから、勢い良く燃料を注いだ場合、その勢いによって、航空機の運動は、乱されてしまうでしょう。コレは、航空機が、供給された燃料から運動量を受け取ったコトを、意味しています。コレを防ぐには、供給される燃料と航空機の相対速度を、常に、ゼロにしておく必要があります。同じように、ロケットに乗っている観測者から見て、外部から受け取る運動量をゼロにするためには、ロケットに供給されるエネルギー塊は、ロケットに対して、相対速度がゼロでなければいけません。よって、この場合、(5.13)式の $u$ は、$u = u_{\mathrm{R}}(t)$ となります。以上の考察によって、ただちに、$p$ と $e$ のあいだに、

\[ p = \frac{u_{\mathrm{R}}(t)}{c^{2}} \, e \tag{5.14} \]

という関係が、成り立っていなければならないコトが分かります。


□(5.13)式は、「質点」のエネルギーと運動量を想定していますけど、もちろん、「電磁波(光子)」のエネルギーと運動量についても、(5.13)式が成り立ちます。電磁波のカタチでエネルギーを供給する場合、電磁波とロケットの進行方向のあいだの角度 $\varphi$ が、$u_{\mathrm{R}}(t) = c \cdot \cos{\varphi}$ という条件から、得られます。


■(5.10)式、(5.12)式、(5.14)式から、$p$ および $e$ を消去すると、(5.8)式も用いて、最終的に、

\[ f + \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(t)/c )^{2}}} \, \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{5.15a} \] \[ \quad \mbox{または} \quad f \, \sqrt{1 - \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} - \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \frac{dm(t)}{dt} = \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{5.15b} \]

という関係式が得られます。この(5.15a)〜(5.15b)式が、「特殊相対論的に正しいツィオルコフスキーの公式」と呼ぶべき、もう 1 つの微分方程式になります。(5.15a)〜(5.15b)式を解くコトで、慣性系にある時計の時刻が $t$ のときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(t)$ が得られるからです。ただし、一般に、(5.15a)〜(5.15b)式を解くのは、極めて難しいと思われます。他方で、(5.15b)式のほうは、外力 $f$ がゼロの場合に、推進剤を含んだロケット全体の質量 $m$ が与えられたときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(m)$ を求めるコトならば、とても簡単です。詳しくは、第 6 節で、検討しましょう。


□なお、慣性系にいる観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(5.15a)式は、慣性系にいる観測者から見た物理量を、ロケットに乗っている観測者から見た物理量に変換する関係式である(3.4)式、および、(4.16)式を用いると、ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(5.4)式に、変換できます。観測者が違うだけで、同じロケットの運動を観測しているので、当然と言えば当然です。


■ちなみに、(5.10)式、(5.12)式、(5.14)式から、$p$ および $e$ ではなく、$f$ のほうを消去すると、

\[ e = \left( \frac{1}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} - 1 \right) \, \mu c^{2} \tag{5.16} \]

という関係式が得られます。この(5.16)式の右辺は、(5.6)式の右辺と、まったく同じです。したがって、このコトからも、この場合に、(5.7a)式の関係($e = \varepsilon$)が成り立つコトが、分かります。


□繰り返しになりますけれど、$e = \varepsilon$ が成り立つのは、「ロケットに乗っている観測者から見て、外部から受け取る運動量は、常にゼロである」と仮定した場合だけですので、注意が必要です。

□ところで、実は、このような仮定をしない、一般の場合でも、「ロケットに乗っている観測者から見て、単位時間当たりに外部から供給されるエネルギーと運動量は、常に一定である」という条件を満たすなら、ロケットの運動は、比較的、簡単に計算できます(いずれにせよ、極めて非現実的な設定ですけれども)。その計算は、ココでは割愛しますけど、興味があれば、ぜひ、挑戦してみて下さいね。


■なお、(5.6)式、ないし、(5.16)式によって、この場合の「特殊相対論的に正しいツィオルコフスキーの公式」=(5.4)式、および、(5.15a)式の左辺第 2 項が($\omega_{\mathrm{P}} > 0$ としたとき)、

\[ \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} = c \mu \, \sqrt{\left( \frac{e}{\mu c^{2}} + 1 \right)^{2} - 1} \quad = c \mu \, \sqrt{\left( \frac{\varepsilon}{\mu c^{2}} + 1 \right)^{2} - 1} \tag{5.17} \]

と表せるコトが分かります。実は、(2.7)式と同じ議論が、ココでも成り立ちます。つまり、(5.17)式が、この場合における、ロケットの推力 $F_{\mathrm{T}}$ を表しています。$\sqrt{\smash[b]{\mathstrut 1 - ( \omega_{\mathrm{P}}/c )^{2}}} \le 1$ ですので、ロケットに乗っている観測者から見た推進剤の消費速度 $\mu$ と噴射速度 $\omega_{\mathrm{P}}$ の条件が同じならば、推進剤の質量欠損からエネルギーが供給される場合のロケットの推力 $F_{\mathrm{T}} = \omega_{\mathrm{P}} \mu$((4.14)式以下の議論を参照)より、外部からエネルギーが供給される場合のロケットの推力 $F_{\mathrm{T}} = \omega_{\mathrm{P}} \mu/\sqrt{\smash[b]{\mathstrut 1 - ( \omega_{\mathrm{P}}/c )^{2}}}$ のほうが、必ず、大きくなります。同じ質量の推進剤を消費しても、質量欠損によって失われる推進剤がない分だけ、外部からエネルギーが供給される場合のロケットのほうが、実際に噴射される推進剤の質量が大きくなるからです。

 ■現実的なエネルギー供給方法を考えた場合 ――  

■最初に述べたように、ココで考えた、ロケットへのエネルギー供給方法は、まったく現実的ではありません。それでは、もう少し、現実的なエネルギー供給方法で計算した場合は、どうなるでしょうか。たとえば、次の図のように、慣性系にある発振器から、ロケットに向かって、光ビーム(レーザーやマイクロ波など)を送って、その光ビームが持つエネルギーを、ロケットが受け取る――という方法ならば、実現しても、おかしくはない気がします。ちなみに、この設定は、太陽から放射される光(光エネルギー)を、太陽電池パネルで受け取って、電力(電気エネルギー)に変換して、その電気エネルギーを使って推進剤を加速・噴射する、「はやぶさ 2」のイオンエンジン $\mu 10$ など、現実の「電気推進ロケット」のエネルギー供給方法と、ほとんど同じです。でも、その場合、太陽までの距離によって、ロケットが太陽光から受け取れる光エネルギーの量が変化してしまうので(太陽から離れるほど、太陽光は弱くなってしまうので)、計算がタイヘンです。そこで、ココでは、光ビームは減衰せず、慣性系にある発振器から出力された光ビームのエネルギーのすべてを、ロケットは完全に受け取るコトができる――と、仮定しましょう。
図11:「特殊相対論的に正しいツィオルコフスキーの公式」の導出:電気推進ロケット編(3)(約45KB)

■さて、この場合、慣性系にある発振器は、単位時間当たり一定の大きさのエネルギー $E$ を出力していても、ロケットの速さ $u_{\mathrm{R}}$ が変化すると、「ドップラー効果」の影響によって、ロケットが受け取る単位時間当たりのエネルギー $\varepsilon$ は、変化してしまいます。よって、ロケットに乗っている観測者から見た推進剤の消費速度 $\mu$ と噴射速度 $\omega_{\mathrm{P}}$ の両方を、一定に保つコトはできません。技術的には、推進剤の消費速度 $\mu$ のほうを一定にして(制御して)、供給されるエネルギーに応じて決まる噴射速度 $\omega_{\mathrm{P}}$ で、推進剤を噴射する仕組みを考えるのが、イチバン、簡単そうです。このとき、ロケットの運動を与える方程式は、ロケットに乗っている観測者から見た場合と、慣性系にいる観測者から見た場合で、結果だけを示すと、それぞれ、

\[ f + \left( \frac{1 - u_{\mathrm{R}}(\tau)/c}{1 + u_{\mathrm{R}}(\tau)/c} \right) \frac{E}{c} + \mu c \, \sqrt{\left\{ 1 + \left( \frac{1 - u_{\mathrm{R}}(\tau)/c}{1 + u_{\mathrm{R}}(\tau)/c} \right) \frac{E}{\mu c^{2}} \right\}^{2} - 1} = \frac{m_{0} - \mu \tau}{1 - (u_{\mathrm{R}}(\tau)/c)^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \] \[ {} \tag{5.18a} \] \[ \begin{align} &f + \left( \frac{1 - u_{\mathrm{R}}(t)/c}{1 + u_{\mathrm{R}}(t)/c} \right) \frac{E}{c} + \mu c \, \sqrt{\left\{ 1 + \left( \frac{1 - u_{\mathrm{R}}(t)/c}{1 + u_{\mathrm{R}}(t)/c} \right) \frac{E}{\mu c^{2}} \right\}^{2} - 1} \nonumber \\[0.5em] &\qquad = \frac{1}{1 - (u_{\mathrm{R}}(t)/c)^{2}} \, \frac{1}{\sqrt{1 - (u_{\mathrm{R}}(t)/c)^{2}}} \left( m_{0} - \mu \int_{0}^{t}{\sqrt{1- \left( \frac{u_{\mathrm{R}}(t)}{c} \right)^{2}} \, dt} \right) \frac{du_{\mathrm{R}}(t)}{dt} \tag{5.18b} \end{align} \]

となります。この方程式を解くのは、外力 $f$ がゼロの場合でも、かなり難しそうです。それでも、コレよりも計算が簡単になる、現実的な外部からのエネルギー供給方法は、なかなか、見つけられないと思います。


□余談ですけど、慣性系にある発振器から届く光ビームのエネルギーを、推進剤を噴射するためではなく、そのまま反射して推力を得るために使う、いわゆる「宇宙ヨット」を、考えるコトができます。すなわち、慣性系にある発振器から、大きな「帆」を持つ宇宙ヨットに向かって、光ビームを送ります。宇宙ヨットの帆が、その光ビームを反射するコトで生じる反作用によって、推力を得る仕組みです。地球上のヨットが、その帆に風を受けて、自由自在に進むのと、同じ原理です(風圧の代わりに、「光圧」を使うワケです)。実際に、小型ソーラー電力セイル実証機「IKAROS」が、この原理で、宇宙航行に成功しています。実は、慣性系にある発振器から届く光ビームを受けて航行する宇宙ヨットの運動は、比較的、簡単に計算するコトができます。宇宙ヨットは、定義上、「ロケット」には当てはまらないので、ココでは扱いませんけれど、現実的というだけでなく、とてもオモシロいので、興味があれば、ぜひ、挑戦してみて下さいね。


 ■ロケットに搭載されたエネルギー源から、エネルギーが供給される場合 ――  

■惑星探査機「ボイジャー」は、原子力電池を搭載しています。もちろん、この原子力電池は、電子機器に電力を供給するためのモノで、「ボイジャー」は、「電気推進ロケット」ではありません。けれども、ロケット本体に搭載したエネルギー源(原子力電池でも、核融合炉でも、なんでも OK)から得られるエネルギーを使って、推進剤を加速・噴射するような「電気推進ロケット」を、考えるコトができます。このような電気推進ロケットは、外部から十分なエネルギーが得られない場合(たとえば、太陽光がほとんど届かないトコロを目指すとか、大きな太陽電池パネルを装備できないとか)に、宇宙空間を航行するための有力な手段となるかも知れません。ただし、このような「電気推進ロケット」は、ココ(第 5 節)で検討した、「推進剤を噴射するために必要なエネルギーが、外部から供給される場合」には、当然のコトながら、当てはまりません。


□その意味で、この第 5 節のサブタイトルが、「電気推進ロケット編」であるのは、極めて不適切です。でも、文字数などの制限上、適当なサブタイトルを決めるコトができませんでした……許して下さいね。


■ところで、ロケット本体に搭載したエネルギー源から得られるエネルギーは、どのようなエネルギー源でも、究極的には、質量欠損に由来すると、みなすコトができます。たとえば、原子力電池なら、得られるエネルギーは、原子核崩壊、すなわち、原子核の質量欠損に由来します。したがって、このような「電気推進ロケット」の運動は、第 4 節で説明した、「質量=運動エネルギー(質量欠損)に由来するエネルギーを用いて、推進剤を噴射するロケット」における「特殊相対論的に正しいツィオルコフスキーの公式」に従うコトになるハズです。ただし、この場合、推進剤には質量欠損が生じないので、消費した推進剤は、その質量のすべてが、実際に噴射されます。そのため、第 4 節での議論よりも、推進剤の消費速度 $\mu$ の影響が、大きくなります。具体的には、第 4 節で得られた「特殊相対論的に正しいツィオルコフスキーの公式」の、推進剤の消費速度 $\mu$ を、すべて、

\[ \mu \quad \longrightarrow \quad \mu' = \frac{\mu}{\sqrt{1 - (\omega_{\mathrm{P}}/c)^{2}}} > \mu \tag{5.19} \]

で与えられる $\mu'$ で置き換えると、「ロケット本体に搭載したエネルギー源から得られるエネルギーを用いて、推進剤を噴射するロケット」における「特殊相対論的に正しいツィオルコフスキーの公式」になります。


□もちろん、以上の議論は、ココまでのすべての議論と同じく、「推進剤の噴射速度 $\omega_{\mathrm{P}}$、および、推進剤の消費速度 $\mu$ が、ともに一定である」という前提のもとでの結果であるコトに、注意しましょう。また、ココまでのすべての議論と同じく、この場合も、ロケットに乗っている観測者から見た、単位時間当たりに消費されるエネルギー(質量欠損)は、一定でなければいけないコトが、簡単な計算から分かります。

□ちなみに、ロケット本体に搭載したエネルギー源が、静電エネルギーなど、質量以外のエネルギーに由来する場合でも、$E = mc^{2}$ によって、そのエネルギーが、どれだけの質量の質量欠損に相当するのか、換算してから計算するコトで、まったく同じ議論が成り立ちます。


■さて、推進剤を噴射するためのエネルギーが、推進剤自身の質量欠損に由来する場合と、ロケット本体に搭載されたエネルギー源に由来する場合の違いは、数式上は、「特殊相対論的に正しいツィオルコフスキーの公式」に含まれる $\mu$ が、(5.19)式で与えられる $\mu'$ に、置き換わるだけです。けれども、その帰結として、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に関する議論に、大きな違いが生じます。まず、推進剤を噴射するためのエネルギーが、推進剤自身の質量欠損に由来する場合は、原理的には、最初に搭載する推進剤の質量を、無限に増加させる(正確には、「質量比 $R = m_{0}/m_{\mathrm{R}}$($m_{\mathrm{R}}$:推進剤を除いたロケット本体の質量。$m_{0}$:推進剤を含めた加速前のロケット全体の質量)」を、無限に増加させる)コトで、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を、いくらでも、光速 $c$ に近づけるコトができます(詳しくは、第 6 節で計算します)。


□もちろん、搭載する推進剤を増やすには、より大きな推進剤タンクが必要になり、その分、ロケット本体の質量が増加してしまうでしょう。重いロケットを加速するために、その分、余計に推進剤が必要になり、さらに大きな推進剤タンクが必要になり、ロケット本体の質量が増加して、……というコトになって、実際には、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、光速 $c$ よりも遥かに遅い速さで、限界を迎えるハズです。


■他方で、推進剤を噴射するためのエネルギーが、ロケット本体に搭載されたエネルギー源に由来する場合は、搭載する推進剤を増やしても、それだけでは、推進剤を噴射するためのエネルギーは増えません。したがって、推進剤に十分な噴射速度を与えるコトができず、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、期待するほどには向上しません(推進剤を噴射するためのエネルギーが、推進剤自身の質量欠損に由来する場合は、推進剤を増やした分だけ、同時に、推進剤を噴射するためのエネルギーも増えるコトになる――という点に、注意しましょう)。推進剤を噴射するためのエネルギーを増やすには、ロケット本体に搭載されるエネルギー源に蓄えるエネルギー(質量欠損)を増やす必要があります。でも、そうすると、ロケット全体の質量も増えてしまうので、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ への影響は、複雑です。実は、推進剤を噴射するためのエネルギーが、ロケット本体に搭載されたエネルギー源に由来する場合には、そのエネルギー源のエネルギー密度 $E_{\rho}$ に応じて、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を最大にするような、最適な推進剤の質量(ロケット本体の質量に対する推進剤の質量の割合)が存在します。詳しい説明は割愛しますけど、計算自体は簡単ですので、興味があれば、ぜひ、最適な推進剤の質量と、そのときのロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ の計算にも、挑戦してみて下さいね。


□ちなみに、ロケット本体に搭載するエネルギー源が、核融合炉の場合には、そのエネルギー密度 $E_{\rho}$ を、$6.3 \times 10^{14}$ [$\mathrm{J/kg}$](pp 連鎖反応。核燃料以外の部分の質量をゼロとした場合)以上にするコトは、原理的に不可能です。このとき、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、$2.9 \times 10^{4}$ [$\mathrm{km/s}$]($0.095 \, c$)が、理論的な上限値となります。太陽系内を航行するだけならば、十分かも知れませんけれど、星間旅行は難しそうです(最も近い恒星でも、到着まで、片道で 40 年以上もかかるコトになります)。ましてや、ウラシマ効果やスターボウといった、特殊相対論の世界に特有の現象を体感するコトは、ほとんど、できないでしょう。

□さらに、現実的なハナシとして、ロケット本体に搭載するエネルギー源に、原子力電池を採用する場合、そのエネルギー密度 $E_{\rho}$ を、(核燃料以外の部分の質量は無視して)$8.2 \times 10^{13}$ [$\mathrm{J/kg}$](ウランの核分裂)と考えると、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、$1.0 \times 10^{4}$ [$\mathrm{km/s}$]($0.034 \, c$)が、理論的な上限値です。

□もっと身近なハナシになると、スマホなどに使われているリチウムイオン電池のエネルギー密度 $E_{\rho}$ は、最も高性能なモノでも、$1.0 \times 10^{6}$ [$\mathrm{J/kg}$] 程度のようです。すると、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、$1.1$ [$\mathrm{km/s}$] が、理論的な上限値となります。小惑星探査機「はやぶさ(初号機)」のイオンエンジン $\mu 10$ が、地球−小惑星イトカワの往復において、延べ $39600$ [時間](4 基合計)の運転で、$2.2$ [$\mathrm{km/s}$] の増速を果たしたのと比較すると、少し物足りなく感じるかも知れません。

□なお、以上の計算例は、「ロケット本体の質量が、限りなくゼロに近い」と仮定して、求めたモノです。もちろん、ロケット本体の質量をゼロにするのは不可能なので、実際のロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ は、コレらの計算例よりも、遥かに遅いモノになるハズです。

■このページの一番上に戻る。



 ■6. ツィオルコフスキーの公式を解く ―― ロケットの運動と時間の関係

■ココまでに、いろいろな設定における、「特殊相対論的に正しいツィオルコフスキーの公式」を、微分方程式のカタチで、導いてきました。得られた微分方程式を解くコトができれば、特殊相対論の世界での、ロケットの運動について、さらに、いろいろなコトが分かります。最後に、ココまでに得られた、「特殊相対論的に正しいツィオルコフスキーの公式」を、実際に解いてみて、特殊相対論の世界での、ロケットの運動について、詳しく調べてみましょう。

■まず、外力 $f$ がゼロの場合の「特殊相対論的に正しいツィオルコフスキーの公式」について、検討します。この場合には、ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ が一定ならば、観測者がロケットに乗っているのか、慣性系にいるのかにかかわらず、ロケットの質量が $m$ のときの、ロケットの速さ $u_{\mathrm{R}}(m)$ を求めるコトは、とても簡単です(この場合は、推進剤の消費速度には、特に、条件を課す必要はありません)。次に、外力 $f$ がゼロでない場合を、考えます。この場合には、ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ だけでなく、外力 $f$ と、ロケットに乗っている観測者から見た推進剤の消費速度 $\mu$ も、それぞれ一定ならば、時刻が $\tau$ のときの、ロケットの速さ $u_{\mathrm{R}}(\tau)$ のほか、時刻 $\tau$ と時刻 $t$ の関係(「ウラシマ効果」と呼ばれる現象)など、さらに、さまざまなコトが、簡単に計算できます。早速、計算してみましょう。

 ■外力がゼロの場合:ロケットの質量と速さの関係 ――  

まずは、外力 $f$ がゼロの場合を、考えましょう。第 4 節では、推進剤を噴射するエネルギーが、推進剤自身の質量欠損に由来する場合の、「特殊相対論的に正しいツィオルコフスキーの公式」を導きました。ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(4.15)式、慣性系にいる観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(4.30)式は、推進剤の消費速度 $\mu$ が(4.1)式で定義されていたので、外力 $f$ がゼロのとき、それぞれ、次のように、書き直すコトができます。

\[ - \omega_{\mathrm{P}} \, \frac{dm(\tau)}{d\tau} = \frac{m(\tau)}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.1a} \] \[ - \omega_{\mathrm{P}} \, \frac{dm(t)}{dt} = \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{6.1b} \]

また、第 5 節で考えたような、(ちょっと特殊な)外部からのエネルギー供給方法で、推進剤を噴射する場合は、ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(5.3)式、慣性系にいる観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(5.15b)式は、同様に、外力 $f$ がゼロのとき、それぞれ、次のように、書き直すコトができます。

\[ - \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \frac{dm(\tau)}{d\tau} = \frac{m(\tau)}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.2a} \] \[ - \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \, \frac{dm(t)}{dt} = \frac{m(t)}{1 - ( u_{\mathrm{R}}(t)/c )^{2}} \, \frac{du_{\mathrm{R}}(t)}{dt} \tag{6.2b} \]


■ココで、「ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、一定とする」と、仮定しましょう(推進剤の消費速度については、ココでは、特に、条件を課しません)。このとき、(6.1a)〜(6.2b)式の左辺の $dm(\tau)/d\tau, \ dm(t)/dt$ の係数は、いずれも、時刻に依らない定数となります。そこで、便宜的に、新たな記号 $\omega_{\mathrm{P}}'$ を、

\[ \mbox{第 4 節の場合:} \ \omega_{\mathrm{P}}' \equiv \omega_{\mathrm{P}} \, , \qquad \mbox{第 5 節の場合:} \ \omega_{\mathrm{P}}' \equiv \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \tag{6.3} \]

のように定義します(もちろん、$\omega_{\mathrm{P}}'$ は、時刻に依らない定数です)。ちなみに、$\omega_{\mathrm{P}}$ と $\omega_{\mathrm{P}}/(\sqrt{\smash[b]{\mathstrut 1-(\omega_{\mathrm{P}}/c)^{2}}}\, )$ の関係は、次の表のようになります。

\[ \hspace{-2em} \begin{array}{c|cccccccccc} \omega_{\mathrm{P}} & 0.1 \, c & 0.2 \, c & 0.3 \, c & 0.4 \, c & 0.5 \, c & 0.6 \, c & 0.7 \, c & 0.8 \, c & 0.9 \, c & c \\ \hline \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} & 0.101 \, c & 0.204 \, c & 0.315 \, c & 0.436 \, c & 0.577 \, c & 0.750 \, c & 0.980 \, c & 1.333 \, c & 2.065 \, c & \infty \end{array} \]

さらに、便宜的に、どの観測者が持つ時計の時刻を考えるのかに応じて、新たな記号 $\sigma$ で $\sigma \equiv \tau$ または $\sigma \equiv t$ を表すとすると、(6.1a)〜(6.2b)式は、まとめて、次のように表すコトができます。

\[ - \omega_{\mathrm{P}}' \, \frac{dm(\sigma)}{d\sigma} = \frac{m(\sigma)}{1 - ( u_{\mathrm{R}}(\sigma)/c )^{2}} \, \frac{du_{\mathrm{R}}(\sigma)}{d\sigma} \qquad \mbox{($\sigma$ は、$\sigma \equiv \tau$ または $\sigma \equiv t$ を表す)} \tag{6.4} \]

したがって、(6.4)式を解けば、(6.1a)〜(6.2b)式の、すべての微分方程式が、解けたコトになります。

■ところで、(6.4)式は、$m(\sigma)$ の具体的なカタチが分からなくても(たとえ、$m(\sigma)$ が、$m(\sigma) = m_{0}-\mu \sigma$ で表されないとしても)、ロケット全体の質量が $m$ のときの、ロケットの速さ $u_{\mathrm{R}}$(すなわち、$u_{\mathrm{R}} = u_{\mathrm{R}}(m)$)を与える微分方程式とみなした場合には、簡単に解くコトができます。いろいろな解法や表記法がありますけど、一例として、ココでは、次のように解いてみましょう。まず、(6.4)式を、次のように変形します。

\[ - \frac{\omega_{\mathrm{P}}'}{m(\sigma)} \, \frac{dm(\sigma)}{d\sigma} = \frac{1}{1 - ( u_{\mathrm{R}}(\sigma)/c )^{2}} \, \frac{du_{\mathrm{R}}(\sigma)}{d\sigma} \tag{6.5} \]

次に、(6.5)式の両辺を被積分関数として、時刻 $\sigma = 0$ から時刻 $\sigma$ までの区間で積分をしてみます。すなわち($\omega_{\mathrm{P}}'$ が、時刻 $\sigma$ に依らない定数であったコトに注意して)、

\[ -\omega_{\mathrm{P}}' \int_{0}^{\sigma}{\frac{1}{m(\sigma)} \, \frac{dm(\sigma)}{d\sigma}} \, d\sigma = \int_{0}^{\sigma}{\frac{1}{1 - ( u_{\mathrm{R}}(\sigma)/c )^{2}} \, \frac{du_{\mathrm{R}}(\sigma)}{d\sigma}} \, d\sigma \tag{6.6} \]

さらに、左辺は $m = m(\sigma)$、右辺は $u_{\mathrm{R}} = u_{\mathrm{R}}(\sigma)$ で、それぞれ、置換積分してみましょう。ただし、時刻 $\sigma = 0$ のとき、それぞれ、$m_{0} \equiv m(0), \ u_{\mathrm{R0}} \equiv u_{\mathrm{R}}(0)$ と書くコトにします(それぞれ、時刻 $\sigma = 0$ のときの、推進剤を含んだロケット全体の質量、および、ロケットの速さ(初速度)を表しています)。すると、置換積分の公式を用いて、(6.6)式を、次のように書き直すコトができます。

\[ -\omega_{\mathrm{P}}' \int_{m_{0}}^{m}{\frac{1}{m}} \, dm = \int_{u_{\mathrm{R0}}}^{u_{\mathrm{R}}}{\frac{1}{1 - ( u_{\mathrm{R}}/c )^{2}}} \, du_{\mathrm{R}} \tag{6.7} \]

(6.7)式の左辺の積分は、簡単です。右辺の積分も、(2.30)式と同じ要領で、「部分分数分解」と呼ばれる数学テクニックを用いると、簡単に計算できます。結果だけを示すと、

\[ \omega_{\mathrm{P}}' \, \ln{\left( \frac{m_{0}}{m} \right)} = \frac{c}{2} \left\{ \ln{\left( \frac{1 + u_{\mathrm{R}}/c}{1 - u_{\mathrm{R}}/c} \right)} - \ln{\left( \frac{1 + u_{\mathrm{R0}}/c}{1 - u_{\mathrm{R0}}/c} \right)} \right\} \tag{6.8a} \] \[ \begin{align} \quad \mbox{したがって} \quad u_{\mathrm{R}} &= c \ \frac{U_{0} \, \left( m_{0}/m \right)^{2\omega_{\mathrm{P}}'/c} - 1}{U_{0} \, \left( m_{0}/m \right)^{2\omega_{\mathrm{P}}'/c} + 1} \quad = c \ \frac{U_{0} \, \left( m_{0} \right)^{2\omega_{\mathrm{P}}'/c} - \left( m \right)^{2\omega_{\mathrm{P}}'/c}}{U_{0} \, \left( m_{0} \right)^{2\omega_{\mathrm{P}}'/c} + \left( m \right)^{2\omega_{\mathrm{P}}'/c}} \tag{6.8b} \\[0.5em] &\hspace{15em} \left( \, \mbox{ただし、} U_{0} \equiv \frac{1 + u_{\mathrm{R0}}/c}{1 - u_{\mathrm{R0}}/c} \ \right) \nonumber \end{align} \]

となります(なお、$U_{0}$ は、初速度 $u_{\mathrm{R0}}$ だけに依存する定数で、$u_{\mathrm{R0}} = 0$ のときに、$U_{0} = 1$ となるコトに注意。また、初速度の絶対値は、光速 $c$ を超えられないので、$0 \le U_{0} \le \infty$ であるコトにも注意)。(6.8b)式より、推進剤の噴射速度 $\omega_{\mathrm{P}}$ が一定のとき、(推進剤の消費速度が、一定かどうかにかかわらず)ロケットの速さ $u_{\mathrm{R}}$ は、そのときの、推進剤を含んだロケット全体の質量 $m$ の値だけで決まる――というコトが分かります。


□(6.8b)式より、$\omega_{\mathrm{P}}' > 0$ の場合、$du_{\mathrm{R}}/dm < 0$、および、$u_{\mathrm{R}} \, \to \, c$($m \, \to \, 0$)であるコトが分かります。つまり、推進剤を含んだロケット全体の質量 $m$ が減少するにしたがって、ロケットの速さ $u_{\mathrm{R}}$ は増加し、$m \, \to \, 0$ の極限で、光速 $c$ に達するコトになります。もちろん、ロケットの質量 $m$ がゼロになってしまうコトはありえないので、ロケットは、光速 $c$ に到達する前に、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達してしまうコトになります。(6.8b)式は、その最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を与える公式であると、みなすコトもできます。


 ■外力、および、推進剤の消費速度が、ともに一定の場合(1):ロケットの時刻と速さの関係 ――  

次に、外力 $f$ が一定の場合を、考えましょう。なお、「外力が一定の場合」というのは、常に同じ大きさの外力が、ロケットにはたらく場合だけでなく、ロケットにはたらく外力が、常にゼロである場合も含まれます。


□むしろ、特殊相対論の世界で、ロケットの運動を考える場合には、「ロケットにはたらく外力はゼロ」、または、「ロケットにはたらく外力は(十分に小さいので)無視する」として、計算すべき状況のほうが、圧倒的に多いかも知れません。そもそも、「ロケットに一定の外力がはたらく」という状況は、ほとんど、考えられません。でも、「外力がゼロ」の場合は、「一定の外力がはたらく」という条件での計算結果に、$f = 0$ を代入すれば良いだけなので、ココでは、あえて、この設定で、計算を続けてみます。


■さらに、「ロケットに乗っている観測者から見た推進剤の消費速度 $\mu$ は、一定である」と仮定しましょう。もちろん、これまでと同じく、「ロケットに乗っている観測者から見た推進剤の噴射速度 $\omega_{\mathrm{P}}$ は、一定である」と仮定します。このように、$f, \ \mu, \ \omega_{\mathrm{P}}$ のすべてが常に一定の場合には、コレまでに得られた「特殊相対論的に正しいツィオルコフスキーの公式」から、より多くの、とても興味深い結果が導かれます。最初に、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$ の関係(すなわち、時刻 $\tau$ のときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(\tau)$)を、求めてみましょう。

■さて、第 4 節で考えた、推進剤を噴射するエネルギーが、推進剤自身の質量欠損に由来する場合における、ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(4.17)式は、

\[ f + \omega_{\mathrm{P}} \mu = \frac{m_{0} - \mu \tau}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.9} \]

また、第 5 節で考えた、(ちょっと特殊な)外部からのエネルギー供給方法で、推進剤を噴射する場合には、ロケットに乗っている観測者から見た「特殊相対論的に正しいツィオルコフスキーの公式」=(5.4)式は、

\[ f + \frac{\omega_{\mathrm{P}} \mu}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} = \frac{m_{0} - \mu \tau}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.10} \]

でした。いま、$f, \ \mu, \ \omega_{\mathrm{P}}$ は、すべて、常に一定と仮定していたので、ココで、便宜的に、新たな記号 $\varOmega_{\mathrm{P}}$ を、

\[ \mbox{第 4 節の場合:} \ \varOmega_{\mathrm{P}} \equiv \frac{f}{\mu} + \omega_{\mathrm{P}} \, , \qquad \mbox{第 5 節の場合:} \ \varOmega_{\mathrm{P}} \equiv \frac{f}{\mu} + \frac{\omega_{\mathrm{P}}}{\sqrt{1 - ( \omega_{\mathrm{P}}/c )^{2}}} \tag{6.11} \]

のように定義します(もちろん、$\varOmega_{\mathrm{P}}$ は、時刻に依らない定数です)。そうすると、(6.9)〜(6.10)式は、まとめて、次のように表すコトができます。

\[ \varOmega_{\mathrm{P}} \mu = \frac{m_{0} - \mu \tau}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.12} \]

したがって、(6.12)式を解けば、(6.9)〜(6.10)式の微分方程式が、どちらも解けたコトになります。

■ところで、実際に、(6.12)式は、簡単に解くコトができます。いろいろな解法や表記法がありますけど、一例として、ココでは、次のように解いてみましょう。まず、(6.12)式を、次のように変形します。

\[ \frac{\varOmega_{\mathrm{P}} \mu}{m_{0} - \mu \tau} \, d\tau = \frac{1}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, du_{\mathrm{R}}(\tau) \tag{6.13} \]

次に、この式の両辺を、時刻 $\tau = 0$ から時刻 $\tau$ の区間で積分します。ただし、時刻 $\tau = 0$ のときの、ロケットの速さ(初速度)を、$u_{\mathrm{R0}} \equiv u_{\mathrm{R}}(0)$ と書くコトにしましょう。すると、考える積分は、

\[ \varOmega_{\mathrm{P}} \mu \int_{0}^{\tau}{\frac{1}{m_{0} - \mu \tau} \, d\tau} = \int_{u_{\mathrm{R0}}}^{u_{\mathrm{R}}(\tau)}{\frac{1}{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}} \, du_{\mathrm{R}}(\tau)} \tag{6.14} \]

と表されます($\varOmega_{\mathrm{P}}, \ \mu, \ m_{0}$ は、時刻 $\tau$ に依らない定数であったコトに注意)。この左辺の積分は、簡単です。さらに、右辺の積分も、(2.30)式と同じ要領で、「部分分数分解」と呼ばれる数学テクニックを用いると、簡単に計算できます。結果だけを示すと、

\[ \varOmega_{\mathrm{P}} \, \ln{\left( \frac{m_{0}}{m_{0} - \mu \tau} \right)} = \frac{c}{2} \left\{ \ln{\left( \frac{1 + u_{\mathrm{R}}(\tau)/c}{1 - u_{\mathrm{R}}(\tau)/c} \right)} - \ln{\left( \frac{1 + u_{\mathrm{R0}}/c}{1 - u_{\mathrm{R0}}/c} \right)} \right\} \tag{6.15a} \] \[ \begin{align} \quad \mbox{したがって} \quad u_{\mathrm{R}}(\tau) &= c \ \frac{U_{0} \, \left( \frac{m_{0}}{m_{0} \ - \, \mu \tau} \right)^{2\varOmega_{\mathrm{P}}/c} - 1}{U_{0} \, \left( \frac{m_{0}}{m_{0} \ - \, \mu \tau} \right)^{2\varOmega_{\mathrm{P}}/c} + 1} \quad = c \ \frac{U_{0} - \left(1 - \frac{\mu}{m_{0}} \tau \right)^{2\varOmega_{\mathrm{P}}/c}}{U_{0} + \left(1 - \frac{\mu}{m_{0}} \tau \right)^{2\varOmega_{\mathrm{P}}/c}} \tag{6.15b} \\[0.5em] &\hspace{15em} \left( \, \mbox{ただし、} U_{0} \equiv \frac{1 + u_{\mathrm{R0}}/c}{1 - u_{\mathrm{R0}}/c} \ \right) \nonumber \end{align} \]

となります(なお、$U_{0}$ は、初速度 $u_{\mathrm{R0}}$ だけに依存する定数で、$u_{\mathrm{R0}} = 0$ のときに、$U_{0} = 1$ となるコトに注意。また、初速度の絶対値は、光速 $c$ を超えられないので、$0 \le U_{0} \le \infty$ であるコトにも注意)。この(6.15b)式によって、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$(すなわち、時刻 $\tau$ のときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(\tau)$)の関係が、与えられる――というワケです(ただし、$f, \ \omega_{\mathrm{P}}, \ \mu$ は、すべて、時刻 $\tau$ に依らない定数としたコトに注意)。


□ちなみに、(6.15a)式の右辺に、良く知られた近似式 $\ln{(1+X)} \simeq X$($\lvert X \rvert \ll 1$)を適用すると、
\[ \qquad u_{\mathrm{R}}(\tau) \simeq u_{\mathrm{R0}} + \varOmega_{\mathrm{P}} \, \ln{\left( \frac{m_{0}}{m_{0} - \mu \tau} \right)} \quad \mbox{($\, \lvert u_{\mathrm{R}}(\tau) \rvert, \ \lvert u_{\mathrm{R0}} \rvert \ll c \,$)} \tag{6.16} \]
という関係式が得られます。実は、この(6.16)式は、「ニュートン力学的に正しいツィオルコフスキーの公式」=(2.6)式を、$f, \ \omega_{\mathrm{P}}, \ \mu$ が一定という条件で解いた結果と、まったく同じカタチをしています。つまり、$\lvert u_{\mathrm{R}}(\tau) \rvert, \ \lvert u_{\mathrm{R0} \,} \rvert \ll c$ のとき、特殊相対論の世界とニュートン力学の世界は、一致するのです。


■いま導いた(6.15b)式、つまり、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$(時刻 $\tau$ のときの、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}(\tau)$)を、グラフに表すと、次の図のようになります(ただし、この図は、$\varOmega_{\mathrm{P}} > 0$、および、ロケットの初速度がゼロ、すなわち、$u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合を、描いているコトに、注意して下さい)。この図には、さまざまな値の $\varOmega_{\mathrm{P}}$ について、特殊相対論の世界での計算結果(青いグラフ)と一緒に、同じ条件のもとで、ニュートン力学の世界で計算した結果(黒いグラフ)も、描かれていますので、ぜひ、比較してみて下さいね。
図12:時刻とロケットの速さの関係(1):特殊相対論の世界(約67KB)

■ココで、(6.15b)式を、さらに詳しく検討してみましょう。まず、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、推進剤を含んだロケット全体の質量 $m(\tau)$ について、$m(\tau) = m_{0}-\mu\tau > 0$ であるコトは明らかですので、$\varOmega_{\mathrm{P}} > 0$ のとき、$du_{\mathrm{R}}(\tau)/d\tau > 0$ が、常に成り立つコトは、すぐに分かります。コレは、$u_{\mathrm{R}}(\tau)$ が、単調増加関数であるコトを、表しています。さらに、$d^{2}u_{\mathrm{R}}(\tau)/d\tau^{2}$ については、$\varOmega_{\mathrm{P}} > 0$ のとき、次の関係が成り立っているコトが、分かります(ちょっとメンドくさいですけど、計算自体は簡単です)。

\[ 0 < \varOmega_{\mathrm{P}} \le \frac{c}{2} \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(\tau)}{d\tau^{2}} > 0 \tag{6.17a} \] \[ \varOmega_{\mathrm{P}} > \frac{c}{2} \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(\tau)}{d\tau^{2}} \ge 0 \quad \left( \ \tau \le \frac{m_{0}}{\mu} \left\{ 1 - \left( U_{0} \, \frac{2\varOmega_{\mathrm{P}} - c}{2\varOmega_{\mathrm{P}} + c} \right)^{c/(2\varOmega_{\mathrm{P}})} \ \right\} \ \right) \] \[ \hspace{8.3em} \frac{d^{2} u_{\mathrm{R}}(\tau)}{d\tau^{2}} < 0 \quad \left( \ \tau > \frac{m_{0}}{\mu} \left\{ 1 - \left( U_{0} \, \frac{2\varOmega_{\mathrm{P}} - c}{2\varOmega_{\mathrm{P}} + c} \right)^{c/(2\varOmega_{\mathrm{P}})} \ \right\} \ \right) \tag{6.17b} \]

つまり、$0 < \varOmega_{\mathrm{P}} \le c/2$ の場合、ロケットの加速度は、時間が経過するとともに、常に大きくなっていきます。でも、$\varOmega_{\mathrm{P}} > c/2$ の場合には、ある時刻 $\tau_{\mathrm{C}}$ までは、時間の経過とともに、加速度が大きくなっていくものの、その時刻 $\tau_{\mathrm{C}}$ を境に、逆に、時間の経過とともに、加速度は小さくなっていく――そのような特別な時刻 $\tau_{\mathrm{C}}$、すなわち、「($u_{\mathrm{R}}(\tau)$ のグラフにおける)変曲点」が存在します。(6.17b)式より、この変曲点 $\tau_{\mathrm{C}}$ は、

\[ \varOmega_{\mathrm{P}} > \frac{c}{2} \mbox{ の場合:} \ \tau_{\mathrm{C}} = \frac{m_{0}}{\mu} \left\{ 1 - \left( U_{0} \, \frac{2\varOmega_{\mathrm{P}} - c}{2\varOmega_{\mathrm{P}} + c} \right)^{c/(2\varOmega_{\mathrm{P}})} \ \right\} \, , \quad u_{\mathrm{R}}(\tau_{\mathrm{C}}) = c \, \frac{1}{2(\varOmega_{\mathrm{P}} / c)} \tag{6.18} \]

であるコトが分かります。$u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合、$\tau_{\mathrm{C}}$ の具体的な値は、次の表のようになります。

\[ \hspace{2em} \begin{array}{c|ccccccc} \varOmega_{\mathrm{P}} & 0.2 \, c & 0.4 \, c & 0.6 \, c & 0.8 \, c & c & 1.5 \, c & 2.0 \, c \\ \hline \tau_{\mathrm{C}} & \mbox{なし} & \mbox{なし} & 0.864 \, \frac{m_{0}}{\mu} & 0.600 \, \frac{m_{0}}{\mu} & 0.423 \, \frac{m_{0}}{\mu} & 0.206 \, \frac{m_{0}}{\mu} & 0.120 \, \frac{m_{0}}{\mu} \end{array} \]
□なお、ココに出てくる $m_{0}/\mu$ は、定義より明らかに、$m(\tau) = 0$ となる時刻です((4.16)式を参照)。この表のとおりに、変曲点(すなわち、$u_{\mathrm{R}}(\tau)$ のグラフが、下に凸の曲線から、上に凸の曲線に、変化する点)が存在するコトを、先に示したグラフから、確認してみて下さいね。

□ところで、このような $u_{\mathrm{R}}(\tau)$ のふるまいが生じるのは、なぜでしょうか。それには、2 つの要素が影響しています。まず、推進剤の消費に伴って、ロケットの質量が減少していきます。推力は変わらなくても、質量が減少していくと、加速度は、どんどん増大していきます。ところが、ロケットの速さ $u_{\mathrm{R}}$ が光速 $c$ に近づいていくと、今度は、特殊相対論的効果によって、ロケットの速さは増加しにくくなります。もしも、ロケットの推力が小さいと(つまり、$\varOmega_{\mathrm{P}}\mu$ が小さいと)、特殊相対論的効果が顕著になるほどの速さに、なかなか達しません。その前に、質量減少による加速度増大の効果が大きくなりすぎているので、ロケットの質量がゼロになる直前に、一気に、ロケットの速さが光速に向かってしまう――というワケです。


 ■外力、および、推進剤の消費速度が、ともに一定の場合(2):ロケットの最大到達速度 ――  

■先に示した($\varOmega_{\mathrm{P}} > 0, \ u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合における)$u_{\mathrm{R}}(\tau)$ のグラフには、ロケットの質量 $m(\tau)$ がゼロになる時刻(ロケットの速さが光速に達する時刻)までの、$u_{\mathrm{R}}(\tau)$ のふるまいが描かれていました。もちろん、実際には、ロケットの質量がゼロになるコトはありえません。したがって、ロケットの速さは、光速 $c$ に到達する前に、最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ に達してしまうコトになります。そこで、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ について、少し詳しく、検討してみましょう。

■いま、ロケットが推進剤の噴射を始める前、つまり、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ で、$\tau = 0$ のときの、ロケットに搭載されている推進剤の総質量を $M_{\mathrm{P} \,}$、推進剤を除いたロケット本体の質量を $M_{\mathrm{R}}$ と書くコトにしましょう。ちなみに、これまでと同じく、時刻 $\tau = 0$ のときの、推進剤を含んだロケット全体の質量は、$m_{0}$ と書くコトにします。このとき、$m_{0 \,}$、および、推進剤を使い切るときの時刻 $\tau_{\mathrm{P}}$ が、それぞれ、

\[ m_{0} = M_{\mathrm{P}} + M_{\mathrm{R}} \tag{6.19a} \] \[ \tau_{\mathrm{P}} = \frac{M_{\mathrm{P}}}{\mu} \tag{6.19b} \]

と表されるコトは、明らかです。もちろん、推進剤を使い切ると、ロケットは、それ以上は加速ができません。したがって、(6.15b)式より、ただちに、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ が、次のように与えられます。

\[ u_{\mathrm{R \cdot MAX}} = c \ \frac{U_{0} \, \left( \frac{m_{0}}{M_{\mathrm{R}}} \right)^{2\varOmega_{\mathrm{P}}/c} - 1}{U_{0} \, \left( \frac{m_{0}}{M_{\mathrm{R}}} \right)^{2\varOmega_{\mathrm{P}}/c} + 1} \quad = c \ \frac{U_{0} \, \left( 1 + \frac{M_{\mathrm{P}}}{M_{\mathrm{R}}} \right)^{2\varOmega_{\mathrm{P}}/c} - 1}{U_{0} \, \left( 1 + \frac{M_{\mathrm{P}}}{M_{\mathrm{R}}} \right)^{2\varOmega_{\mathrm{P}}/c} + 1} \tag{6.20} \]


■ところで、第 2 節「ニュートン力学的に正しいツィオルコフスキーの公式」では、消費する推進剤の総質量を $m_{\mathrm{P} \,}$、推進剤を除いたロケット本体の質量を $m_{\mathrm{R}}$ と書きました。でも、第 4 節の設定では、最初に搭載した推進剤が、すべて、実際に噴射されるワケではありません。推進剤の質量の一部(質量欠損)は、推進剤を噴射するためのエネルギーに使われるからです。そこで、ココでは、第 2 節での設定と区別するために、時刻 $\tau = 0$ のときの、ロケットに搭載されている推進剤の総質量を $M_{\mathrm{P} \,}$、推進剤を除いたロケット本体の質量を $M_{\mathrm{R}}$ と、第 2 節とは違う記号で、書いていたワケです。以下では、しばらく、第 4 節の場合(推進剤を噴射するためのエネルギーが、推進剤の質量欠損に由来する場合)について、より詳しく、検討してみます。

■第 4 節の設定(推進剤を噴射するためのエネルギーが、推進剤の質量欠損に由来する場合)では、時刻 $\tau = 0$ のときの、ロケットに搭載されている推進剤の総質量 $M_{\mathrm{P}}$ のうち、質量欠損となる質量の合計 $m_{\mathrm{D}}$ と、実際に噴射される質量の合計 $m_{\mathrm{P}}$ は、それぞれ、(4.18a)〜(4.18b)式と、(6.19b)式より、ただちに、

\[ m_{\mathrm{D}} = \int_{0}^{\tau_{\mathrm{P}}}{\left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, \mu \, d\tau} \quad = \left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, M_{\mathrm{P}} \tag{6.21a} \] \[ m_{\mathrm{P}} = \int_{0}^{\tau_{\mathrm{P}}}{\sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \mu \, d\tau} \quad = \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, M_{\mathrm{P}} \tag{6.21b} \] \[ M_{\mathrm{P}} = m_{\mathrm{D}} + m_{\mathrm{P}} \tag{6.21c} \]

となるコトが分かります(いま、$\mu, \ \omega_{\mathrm{P}}$ は一定と仮定していたコトに注意)。(3.7)式より、質量欠損となる質量の合計 $m_{\mathrm{D}}$ は、ロケットに乗っている観測者から見て、$K_{\kappa} = m_{\mathrm{D}} \, c^{2}$ の質量=運動エネルギーに相当します(ロケットに乗っている観測者から見て、ロケットに搭載されている状態の推進剤は、常に静止しているコトに注意。詳しくは、第 3 節を参照)。すなわち、総質量 $M_{\mathrm{P}}$ の推進剤ならば、推進剤の噴射のために、$m_{\mathrm{D}} \, c^{2}$ のエネルギーに相当する質量欠損がなければいけない――というワケです。したがって、この推進剤に求められるエネルギー密度(つまり、単位質量当たりに得られるエネルギー)$E_{\rho}$ は、(6.21a)式も用いて、

\[ E_{\rho} = \frac{m_{\mathrm{D}} c^{2}}{M_{\mathrm{P}}} = \left( 1 - \sqrt{1 - \left( \frac{\omega_{\mathrm{P}}}{c} \right)^{2}} \, \right) \, c^{2} \tag{6.22} \]

と表されます。(6.22)式を書き直すと、ただちに($\omega_{\mathrm{P}} > 0$ としたとき)、

\[ \omega_{\mathrm{P}} = c \, \sqrt{1 - \left( 1 - \frac{E_{\rho}}{c^{2}} \right)^{2}} \quad = \sqrt{2E_{\rho} - \left( \frac{E_{\rho}}{c} \right)^{2}} \tag{6.23} \]

となります。(6.22)〜(6.23)式は、「推進剤の噴射速度を $\omega_{\mathrm{P}}$ にするには、エネルギー密度 $E_{\rho}$ の推進剤が必要」、または、「エネルギー密度 $E_{\rho}$ の推進剤を使うと、推進剤の噴射速度は $\omega_{\mathrm{P}}$ になる」というコトを、表しています。たとえば、軽水素の核融合(pp 連鎖反応)の場合には、エネルギー密度 $E_{\rho} \simeq 6.3 \times 10^{14}$ [$\mathrm{J/kg}$] より、推進剤の噴射速度 $\omega_{\mathrm{P}} \simeq 3.5 \times 10^{4}$ [$\mathrm{km/s}$]($0.118 \, c$)となります。先に示した、ロケットの速さ $u_{\mathrm{R}}(\tau)$ のグラフを見ると、この噴射速度は、特殊相対論の世界を堪能するには、少し物足りなく感じるかも知れません。他にも、ウランの核分裂の場合は、$E_{\rho} \simeq 8.2 \times 10^{13}$ [$\mathrm{J/kg}$] なので、$\omega_{\mathrm{P}} \simeq 1.3 \times 10^{4}$ [$\mathrm{km/s}$]($0.043 \, c$)となり、水素と酸素の化学反応では、$E_{\rho} \simeq 1.3 \times 10^{7}$ [$\mathrm{J/kg}$] なので、$\omega_{\mathrm{P}} \simeq 5.1$ [$\mathrm{km/s}$] となります。


□もちろん、コレらの計算は、エネルギーロスがまったくない、効率 $100 \, \%$ の場合を、想定しています。実際の推進剤の噴射速度は、コレらの計算結果よりも、遥かに遅くなるハズです。

□定義より、明らかに、$M_{\mathrm{P}} \ge m_{\mathrm{D}}$ ですので、(6.22)式より、$E_{\rho} \le c^{2}$ となります。エネルギー密度 $E_{\rho}$ が最大($E_{\rho} = c^{2}$)となるのは、$M_{\mathrm{P}} = m_{\mathrm{D} \,}$ の場合、すなわち、推進剤の質量が、すべて、質量欠損になる場合です。藤野の知る限り、コレを実現するのは、「対消滅(物質と反物質が合体・消滅して、光子を生成する反応)」だけです。なお、藤野の知る限り、エネルギー密度 $E_{\rho}$ が核融合よりも大きくなりうるのは、対消滅と、「ホーキング放射(ブラックホールの蒸発)」だけです。物質・反物質やブラックホールを自在に操れるようになるまでは、人類が、特殊相対論の世界を堪能するのは、難しいのかも知れません。

□逆に、エネルギー密度 $E_{\rho}$ が $c^{2}$ よりも十分に小さい場合には、(6.23)式より、ただちに、$\omega_{\mathrm{P}} \simeq \sqrt{\smash[b]{\mathstrut 2E_{\rho}}}$($E_{\rho} \ll c^{2}$)となるコトが分かります。コレは、ニュートン力学の世界で、「エネルギーの保存則」から、推進剤の噴射速度 $\omega_{\mathrm{P}}$ を計算した結果と、一致しています。

□ココで、「質量比 $R$」を、$R \equiv m_{0}/M_{\mathrm{R}}$ で定義すると、(6.20)式、および、(6.11)式と(6.23)式より、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ が、$R, \ E_{\rho}$ で決まるコトが分かります。さらに、ただちに、
\[ \qquad \mbox{$\omega_{\mathrm{P}} > 0$ の場合:} \ \frac{d \, u_{\mathrm{R \cdot MAX}}}{dE_{\rho}} \ge 0 \quad \mbox{($E_{\rho} \le c^{2}$)} \tag{6.24a} \]\[ \qquad \mbox{$\varOmega_{\mathrm{P}} > 0$ の場合:} \ \frac{d \, u_{\mathrm{R \cdot MAX}}}{dR} > 0 \, , \qquad u_{\mathrm{R \cdot MAX}} \, \to \, c \quad \mbox{($R \, \to \, \infty$)} \tag{6.24b} \]
となるコトが分かります。特に(6.24b)式は、質量比 $R$ を、いくらでも大きくするコトで、ロケットの最大到達速度 $u_{\mathrm{R \cdot MAX}} \:$ を、光速 $c$ に、いくらでも近づけられるコトを、示しています。


 ■外力、および、推進剤の消費速度が、ともに一定の場合(3):2 つの時刻のあいだの関係 ――  

■ココからは、再び、(6.15b)式に戻りましょう。(6.15b)式で与えられる $u_{\mathrm{R}}(\tau)$ を、慣性系にある時計の時刻 $t$ と、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ のあいだの関係を与える(4.3)式に代入すると、この 2 つの時刻のあいだの関係――つまり、「ウラシマ効果」の具体的なカタチが得られます。すなわち、

\[ \begin{align} dt &= \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}(\tau)/c )^{2}}} \, d\tau \nonumber \\[0.5em] &= \frac{1}{2} \, \left\{ \sqrt{U_{0}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{-\varOmega_{\mathrm{P}}/c} + \frac{1}{\sqrt{U_{0}}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{\varOmega_{\mathrm{P}}/c} \ \right\} \, d\tau \tag{6.25} \end{align} \]

という、微分形式の関係式が得られます。さらに、(6.25)式の両辺を、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ で、時刻 $\tau = 0$ から時刻 $\tau$ までの区間で、積分してみましょう。ただし、簡単のため、ロケットに乗っている観測者が持つ時計の時刻が $\tau = 0$ となるときに、慣性系にある時計の時刻が $t = 0$ になるように、お互いの時計の時刻を合わせておくコトにします。すると、(6.25)式の、考える区間での積分は、

\[ \int_{0}^{t(\tau)}{dt} = \int_{0}^{\tau}{\frac{1}{2} \, \left\{ \sqrt{U_{0}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{-\varOmega_{\mathrm{P}}/c} + \frac{1}{\sqrt{U_{0}}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{\varOmega_{\mathrm{P}}/c} \ \right\} \, d\tau} \tag{6.26} \]

と表されるコトになります。ココで、$t(\tau)$ は、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、慣性系にある時計の時刻を表します。(6.26)式の両辺の積分は簡単です(ただし、右辺の積分については、場合分けが必要なので、注意しましょう)。結果だけを示すと、ただちに、

\[ \begin{align} &\mbox{$\varOmega_{\mathrm{P}} \ne \pm c$ の場合:} \ t(\tau) = \frac{m_{0}}{2\mu} \left[ \ \frac{\sqrt{U_{0}}}{1 - \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{1 - \varOmega_{\mathrm{P}}/c} \ \right\} \right. \nonumber \\[0.5em] &\hspace{15em} \left. + \frac{1}{\sqrt{U_{0}}} \, \frac{1}{1 + \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{1 + \varOmega_{\mathrm{P}}/c} \ \right\} \ \right] \tag{6.27a} \end{align} \] \[ \mbox{$\varOmega_{\mathrm{P}} = c$ の場合:} \ t(\tau) = \frac{m_{0}}{2\mu} \left[ \ -\sqrt{U_{0}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \tau \right)} + \frac{1}{2\sqrt{U_{0}}} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \ \right\} \ \right] \] \[ {} \tag{6.27b} \] \[ \mbox{$\varOmega_{\mathrm{P}} = -c$ の場合:} \ t(\tau) = \frac{m_{0}}{2\mu} \left[ \ -\frac{1}{\sqrt{U_{0}}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \tau \right)} + \frac{\sqrt{U_{0}}}{2} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \ \right\} \ \right] \] \[ {} \tag{6.27c} \]

となります。(6.27a)〜(6.27c)式こそが、ロケットに乗っている観測者が持つ時計の時刻が $\tau$ のときの、慣性系にある時計の時刻 $t(\tau)$ を与える公式です(ただし、ココでは、$f, \ \mu, \ \omega_{\mathrm{P}}$ が常に一定と仮定していたコトに注意)。ロケットの初速度がゼロ、すなわち、$u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合の、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、慣性系にある時計の時刻 $t$ の関係を、グラフに表したのが次の図です。ニュートン力学の世界では、もちろん、この 2 つの時刻に、違いはありません。でも、特殊相対論の世界では、慣性系にある時計の時刻 $t$ が経過=増加しても、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ は、それほど経過=増加しなくなる様子が、この図から、とても良く、見てとれるでしょう。
図13:2 つの時刻のあいだの関係:特殊相対論の世界(約38KB)

■ココで、(6.27a)〜(6.27c)式を、さらに詳しく検討してみましょう。まず、最初に、ロケットの質量がゼロになる時刻(ロケットの速さが光速に達する時刻)は、ロケットに乗っている観測者から見て、$\tau = m_{0}/\mu$ であるのは、明らかです((4.16)式を参照)。それでは、慣性系にいる観測者が観測したとき、ロケットの質量がゼロになる時刻(ロケットの速さが光速に達する時刻)は、どうなるでしょうか。コレは、(6.27a)〜(6.27c)式に、$\tau = m_{0}/\mu$ を代入するコトで、求められます。ロケットの初速度がゼロ($u_{\mathrm{R0}} = 0, \ U_{0} = 1$)の場合の計算結果は、次の表のようになります($\lvert \varOmega_{\mathrm{P}} \rvert > c$ の場合には、$\tau = m_{0}/\mu$ を代入すると、(6.27a)式の右辺に、必ず、$1/0 \to \infty$ となる項が生じるコトに注意)。

\[ \hspace{2em} \begin{array}{c|ccccccc} \varOmega_{\mathrm{P}} & 0.2 \, c & 0.4 \, c & 0.6 \, c & 0.8 \, c & c & 1.5 \, c & 2.0 \, c \\ \hline t \left( \tau = \frac{m_{0}}{\mu} \right) & 1.042 \, \frac{m_{0}}{\mu} & 1.190 \, \frac{m_{0}}{\mu} & 1.563 \, \frac{m_{0}}{\mu} & 2.778 \, \frac{m_{0}}{\mu} & \infty & \infty & \infty \end{array} \]

この結果より、$\lvert \varOmega_{\mathrm{P}} \rvert > c$ のときには、慣性系にいる観測者から見て、ロケットの速さは、決して、光速には到達しない(光速に達するのに、無限の時間がかかる)コトが分かります。

 ■外力、および、推進剤の消費速度が、ともに一定の場合(4):慣性系の時刻とロケットの速さの関係 ――  

■(6.27a)〜(6.27c)式は、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ が与えられたときの、慣性系にある時計の時刻 $t = t(\tau)$ を示しています。もしも、逆に、慣性系にある時計の時刻 $t$ が与えられたときの、ロケットに乗っている観測者が持つ時計の時刻 $\tau = \tau(t)$ が得られたならば、それを、(6.15b)式に代入するコトで、慣性系にある時計の時刻が $t$ のときのロケットの速さ $u_{\mathrm{R}}(t)$ が得られるハズです。とても残念なコトに、(6.27a)〜(6.27c)式を変形して、 $\tau = \tau(t)$ のカタチにするのは、極めて難しそうです(少なくとも、初等関数だけで表すのは、無理だと思われます)。


□実は、$\varOmega_{\mathrm{P}} = \pm c$、かつ、ロケットの初速度がゼロ($u_{\mathrm{R0}} = 0, \ U_{0} = 1$)の場合に限るならば、「ランベルトの W 関数」を用いて、$u_{\mathrm{R}}(t)$ を、簡単なカタチで表すコトができます。ランベルトの W 関数 $W(y)$ は、次のように、$y = x \cdot e^{x}$ の逆関数として定義されます。
\[ \qquad y = x \cdot e^{x} \quad \Longrightarrow \quad x = W(y) \qquad \left( y \ge -\frac{1}{e} \, , \quad W(y) \ge -1 \right) \tag{6.28} \]
さて、$\varOmega_{\mathrm{P}} = \pm c$ の場合の $t(\tau)$ を表す(6.27b)〜(6.27c)式は、$u_{\mathrm{R0}} = 0$($U_{0} = 1$)のとき、どちらも、
\[ \begin{align} &\qquad t(\tau) = \frac{m_{0}}{2\mu} \left[ \, -\ln{\left( 1 - \frac{\mu}{m_{0}} \tau \right)} + \frac{1}{2} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \ \right\} \, \right] \nonumber \\[0em] &\hspace{27em} \mbox{($\, \varOmega_{\mathrm{P}} = \pm c \, , \ u_{\mathrm{R0}} = 0 \,$)} \tag{6.29} \end{align} \]
という、まったく同じ式になります。さらに、ランベルトの W 関数を用いると、
\[ \qquad e^{1 - 4 \, \frac{\mu}{m_{0}} t} = \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \cdot e^{\left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2}} \quad \Longrightarrow \quad \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} = W \left( e^{1 - 4 \, \frac{\mu}{m_{0}} \, t} \right) \] \[ \qquad \quad \mbox{したがって} \quad \tau(t) = \frac{m_{0}}{\mu} \, \left\{ \ 1 - \sqrt{W \left( e^{1 - 4 \, \frac{\mu}{m_{0}} \, t} \right)} \ \right\} \quad \mbox{($\, \varOmega_{\mathrm{P}} = \pm c \, , \ u_{\mathrm{R0}} = 0 \,$)} \tag{6.30} \]
のように、$\tau = \tau(t)$ のカタチに変形できます。コレを(6.15b)式に代入すると、慣性系にある時計の時刻が $t$ のときのロケットの速さ $u_{\mathrm{R}}(t)$ は、
\[ \qquad u_{\mathrm{R}}(t) = c \ \frac{1 - \left\{ W \left( e^{1 - 4 \, \frac{\mu}{m_{0}} \, t} \right) \right\}^{\pm 1}}{1 + \left\{ W \left( e^{1 - 4 \, \frac{\mu}{m_{0}} \, t} \right) \right\}^{\pm 1}} \quad \mbox{($\, \varOmega_{\mathrm{P}} = \pm c \, , \ u_{\mathrm{R0}} = 0 \,$)} \tag{6.31} \]
で与えられます。ランベルトの W 関数 $W(y)$ の性質は、良く知られているので、(6.30)〜(6.31)式によって、$\varOmega_{\mathrm{P}} = \pm c, \ u_{\mathrm{R0}} = 0$ のときの、$\tau = \tau(t)$ や $u_{\mathrm{R}}(t)$ のふるまいを、詳しく調べるコトができます。ココでは、具体的な検討は割愛しますけど、興味があれば、ぜひ、挑戦してみて下さいね。

□ちなみに、第 4 節の設定(推進剤を噴射するためのエネルギーが、推進剤の質量欠損に由来する場合)の「光子ロケット」は、外力 $f$ がゼロのときに、実際に、$\tau = \tau(t)$ や $u_{\mathrm{R}}(t)$ が、(6.30)〜(6.31)式で与えられます。光子ロケットとは、光子を推進剤に使うロケット――つまり、光(電磁波)を発振する反動で推進するロケットのコトです。もちろん、光子ロケットでは、どの観測者から見ても、推進剤の噴射速度が、光速 $c$ になります(すなわち、$\omega_{\mathrm{P}} = u_{\mathrm{P}} = \pm c$ となります)。よって、(6.22)式より、光子ロケットの推進剤のエネルギー密度 $E_{\rho}$ は、$E_{\rho} = c^{2}$ でなければいけないコトが分かります。藤野の知る限り、この条件を満たしうるのは、「物質・反物質」を推進剤とする「対消滅エンジン」だけです。詳しい説明は割愛しますけど、興味があれば、ぜひ、光子ロケットや、対消滅エンジンについて、自分で調べて、いろいろと計算・検討してみて下さいね。


■さて、慣性系にある時計の時刻が $t$ のときのロケットの速さ $u_{\mathrm{R}}(t)$ は、簡単なカタチで表せない――というコトが、明らかになりました。でも、$u_{\mathrm{R}}(\tau)$ と $t = t(\tau)$ ならば、(6.15b)式と(6.27a)〜(6.27c)式で、厳密に与えられています。したがって、数値計算によって(具体的には、(6.27a)〜(6.27c)式の $t = t(\tau)$ に、$t$ の値を代入して、数値計算で、等式を満たす $\tau$ を求めたあとに、その $\tau$ を(6.15b)式の $u_{\mathrm{R}}(\tau)$ に代入するコトで)、慣性系にある時計の時刻が $t$ のときのロケットの速さ $u_{\mathrm{R}}(t)$ の、厳密な値を知るコトができるハズです。実際に、そのような数値計算を用いて、慣性系にある時計の時刻 $t$ とロケットの速さ $u_{\mathrm{R}}(t)$ の関係を、グラフに表したのが次の図です(ただし、この図は、ロケットの初速度をゼロ($u_{\mathrm{R0}} = 0, \ U_{0} = 1$)として、描かれているコトに注意。また、ココでも、もちろん、$f, \ \omega_{\mathrm{P}}, \ \mu$ が常に一定であると仮定されています)。
図14:時刻とロケットの速さの関係(2):特殊相対論の世界(約65KB)


□このグラフより、慣性系にいる観測者が観測する場合も、ロケットに乗っている観測者が観測する場合と同じく、ロケットの質量がゼロになる時刻($m(t) = 0$ となる時刻)に、ロケットの速さ $u_{\mathrm{R}}(\tau)$ が光速 $c$ に達するコトが分かります。ただし、(6.27a)〜(6.27c)式から示されたように、慣性系にいる観測者が観測する場合、$\lvert \varOmega_{\mathrm{P}} \rvert \ge c$ ならば、ロケットの質量がゼロになるには、無限の時間がかかるので、ロケットの速さは、いくらでも、光速に近づいていくものの、いつまでも、光速に達するコトができません。

□もちろん、このコトは、$\lvert \varOmega_{\mathrm{P}} \rvert$(ないし、推進剤の噴射速度 $\omega_{\mathrm{P}}$)が小さいほど、ロケットの速さが光速に達しやすい――というコトではありません。グラフから明らかなように、推進剤の噴射速度 $\omega_{\mathrm{P}}$ が遅い場合は、ロケットの質量がゼロになる直前に、ロケットの速さが、急激に、光速に向かいます。逆に、推進剤の噴射速度 $\omega_{\mathrm{P}}$ が速い場合は、ロケットの質量がゼロに近づく遥か前に、ロケットの速さが、すでに、光速に近くなっているのです。実際、時刻 $t = 0$ のときの、推進剤を含んだロケット全体の質量を $m_{0}$ とすると、ロケットの速さが $0.9 \, c$(光速の $90 \, \%$)に達するときのロケットの質量 $m$ は 、$\varOmega_{\mathrm{P}} = c$ では、$m = 0.23 \, m_{0}$ ですけど、$\varOmega_{\mathrm{P}} = 0.2 \, c$ の場合には、$m = 0.00064 \, m_{0}$ です。コレは、時刻 $t = 0$ のときのロケット全体の質量のうち、$99.9 \, \%$ 以上が推進剤でなければいけないというコトです。$\lvert \varOmega_{\mathrm{P}} \rvert$(ないし、推進剤の噴射速度 $\omega_{\mathrm{P}}$)が十分に大きくないと、光速 $c$ に近い速さまで、ロケットを加速するのは、極めて難しいのです。


■ところで、慣性系にある時計の時刻が $t$ のときのロケットの速さ $u_{\mathrm{R}}(t)$ が、簡単なカタチで表せなくても、$u_{\mathrm{R}}(\tau)$ と $t = t(\tau)$ が、簡単なカタチ((6.15b)式と(6.27a)〜(6.27c)式)で与えられているコトから、$u_{\mathrm{R}}(t)$ のふるまいについて、かなり詳しく、調べるコトができます。たとえば、「合成関数の微分の法則」と、「逆関数の微分公式」より、

\[ \frac{du_{\mathrm{R}}(t)}{dt} = \frac{d\tau}{dt} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} = \frac{1}{dt/d\tau} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \tag{6.32} \]

となりますけど、(6.25)式より $dt/d\tau \ge 1$ ですし、(6.15b)式より $\varOmega_{\mathrm{P}} > 0$ ならば $du_{\mathrm{R}}(\tau)/d\tau > 0$ なので、(6.32)式は、 $\varOmega_{\mathrm{P}} > 0$ ならば $du_{\mathrm{R}}(t)/dt > 0$ であるコトを、示しています。コレは、$u_{\mathrm{R}}(t)$ が、単調増加関数であるコトを、表しています。さらに、$d^{2}u_{\mathrm{R}}(t)/dt^{2}$ については、

\[ \begin{align} \frac{d^{2}u_{\mathrm{R}}(t)}{dt^{2}} &= \frac{d}{dt} \, \left\{ \frac{1}{dt/d\tau} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \right\} = \frac{d\tau}{dt} \, \frac{d}{d\tau} \, \left\{ \frac{1}{dt/d\tau} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \right\} \nonumber \\[0.5em] &= \frac{1}{( dt/d\tau )^{3}} \, \left\{ \frac{dt}{d\tau} \, \frac{d^{2}u_{\mathrm{R}}(\tau)}{d\tau^{2}} - \frac{d^{2}t}{d\tau^{2}} \, \frac{du_{\mathrm{R}}(\tau)}{d\tau} \right\} \tag{6.33} \end{align} \]

となるので、$\varOmega_{\mathrm{P}} > 0$ のとき、次の関係が成り立っているコトが、分かります(かなりメンドくさいですけど、計算自体は簡単です。まず、(6.33)式より、$d^{2}u_{\mathrm{R}}(t)/dt^{2} = 0$ を満たす $\tau$ を求めるコトができます。そして、(6.27a)〜(6.27c)式を使って、その $\tau$ を、$t$ に変換すると、次の関係が得られるハズです)。

\[ 0 < \varOmega_{\mathrm{P}} \le \frac{c}{3} \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} > 0 \tag{6.34a} \] \[ \varOmega_{\mathrm{P}} > \frac{c}{3} \mbox{ の場合:} \ \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} \ge 0 \quad \mbox{($\ t \le t_{\mathrm{C}} \ $)} \, , \quad \frac{d^{2} u_{\mathrm{R}}(t)}{dt^{2}} < 0 \quad \mbox{($\ t > t_{\mathrm{C}} \ $)} \tag{6.34b} \]

ココで、(6.34b)式における $t_{\mathrm{C}}$(「変曲点」)は、次のように与えられます。

\[ \begin{align} &\mbox{$\varOmega_{\mathrm{P}} \ne \pm c$ の場合:} \ t_{\mathrm{C}} \equiv \frac{m_{0}}{2\mu} \left[ \ \frac{\sqrt{U_{0}}}{1 - \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( U_{0} \, \frac{3\varOmega_{\mathrm{P}} - c}{3\varOmega_{\mathrm{P}} + c} \right)^{\frac{1}{2} \, ( c / \varOmega_{\mathrm{P}} \ - \ 1 )} \ \right\} \right. \nonumber \\[0.5em] &\hspace{15em} \left. + \frac{1}{\sqrt{U_{0}}} \, \frac{1}{1 + \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( U_{0} \, \frac{3\varOmega_{\mathrm{P}} - c}{3\varOmega_{\mathrm{P}} + c} \right)^{\frac{1}{2} \, ( c / \varOmega_{\mathrm{P}} \ + \ 1 )} \ \right\} \ \right] \nonumber \end{align} \] \[ {} \tag{6.35a} \] \[ \mbox{$\varOmega_{\mathrm{P}} = c$ の場合:} \ t_{\mathrm{C}} \equiv \frac{m_{0}}{4\mu} \left[ \ \frac{1}{\sqrt{U_{0}}} - \sqrt{U_{0}} \, \left\{ \frac{1}{2} + \ln{\left( \frac{U_{0}}{2} \right)} \right\} \ \right] \tag{6.35b} \]

このように、慣性系にいる観測者が観測するとき、$0 < \varOmega_{\mathrm{P}} \le c/3$ の場合に、ロケットの加速度は、時間が経過するとともに、常に大きくなっていきます。けれども、$\varOmega_{\mathrm{P}} > c/3$ の場合には、変曲点 $t_{\mathrm{C}}$ までは、時間の経過とともに、加速度が大きくなっていくものの、変曲点 $t_{\mathrm{C}}$ を過ぎたとたん、逆に、時間の経過とともに、加速度は小さくなっていきます。質量減少による加速度増大の効果よりも、特殊相対論的効果による加速度減少の効果のほうが、大きくなってしまうからです。ちなみに、$u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合には、$t_{\mathrm{C}}$ の具体的な値は、次の表のようになります。

\[ \hspace{2em} \begin{array}{c|ccccccc} \varOmega_{\mathrm{P}} & 0.2 \, c & 0.4 \, c & 0.6 \, c & 0.8 \, c & c & 1.5 \, c & 2.0 \, c \\ \hline t_{\mathrm{C}} & \mbox{なし} & 1.047 \, \frac{m_{0}}{\mu} & 0.680 \, \frac{m_{0}}{\mu} & 0.438 \, \frac{m_{0}}{\mu} & 0.298 \, \frac{m_{0}}{\mu} & 0.141 \, \frac{m_{0}}{\mu} & 0.081 \, \frac{m_{0}}{\mu} \end{array} \]
□ロケットに乗っている観測者が観測する場合((6.17a)〜(6.17b)式)は、変曲点が存在するのは、$\varOmega_{\mathrm{P}} > c/2$ のときであるのに対して、慣性系にいる観測者が観測する場合には、変曲点が存在する条件が、$\varOmega_{\mathrm{P}} > c/3$ となるのは、とても興味深いコトです。さらに、$\tau_{\mathrm{C}}$ と $t_{\mathrm{C}}$ の表を見比べると、同じ $\varOmega_{\mathrm{P}}$ のとき、必ず、$t_{\mathrm{C}} < \tau_{\mathrm{C}}$ が成り立っているコトに、気がつきます。このコトは、$dt > d\tau$、つまり、$t, \ \tau > 0$ のとき、必ず、$t > \tau$ が成り立つという事実を思い返すと((6.25)式を参照)、不思議に感じるかも知れません。でも、加速度が、「速度の変化量÷時間」で与えられるコトを思い出せば、ナニも不思議ではありません。$dt > d\tau$ なので、ロケットに乗っている観測者が持つ時計で $1$ 秒が経過したときに($d\tau = 1$ [秒])、慣性系にある時計では $1$ 秒以上が経過している($dt > 1$ [秒])――というコトになります。すると、慣性系から見たロケットの速さがある値だけ変化したときに、慣性系にある時計では、ロケットに乗っている観測者が持つ時計よりも、長い時間が経過しているコトになります。「加速度=速度の変化量÷時間」なので、より長い時間が経過する慣性系では、その分だけ、加速度が小さくなります($du_{\mathrm{R}}/dt < du_{\mathrm{R}}/d\tau$ となります)。その結果、より早く、変曲点(加速度が減少し始める時刻)を迎えるコトになる――というワケです。


 ■外力、および、推進剤の消費速度が、ともに一定の場合(5):ロケットの移動距離 ――  

最後に、ロケットの移動距離を、考えてみましょう。ココで、特殊相対論の世界における「距離」について、第 3 節の説明を、思い出す必要があります。なによりも重要なのは、ロケットに乗っている観測者にとって、「距離」を知る手段が、極めて限られる――というコトです。すなわち、ロケットに乗っている観測者自身は、慣性系にいる観測者を頼らずには、距離を測れないのです(詳しくは、『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』を参照)。そんなワケで、ココで言う「ロケットの移動距離」は、必然的に、「慣性系で測った距離 $x$」だけを指すコトになります。


□ちなみに、慣性系にいる観測者を頼れるならば、ロケットに乗っている観測者が、「距離」を知るのに、たとえば、次のような方法があります。まず、慣性系にいる観測者に、慣性系で測った距離で、等間隔に、標識(距離標)を並べておいてもらいます。ロケットに乗っている観測者は、ロケットが、いつ、どの標識を横切ったのかを観測するコトで、ロケットの移動距離(慣性系で測った距離)を、知るコトができます。


■さて、時刻 $t = 0$ から時刻 $t$ までのあいだに、ロケットが進む距離(慣性系で測った距離)$x$ は、定義より、次の式で表されます(正しくは、「変位」と言うべきですけど、分かりやすさを優先して、ココでは、すべて、「移動距離」と書くコトにします)。

\[ x = \int_{0}^{t}{u_{\mathrm{R}}(t) \, dt} \tag{6.36} \]

(6.36)式は、$u_{\mathrm{R}}(t)$ が簡単なカタチで表せないので、このままでは、これ以上は計算できません。けれども、(6.15b)式で $u_{\mathrm{R}}(\tau)$ が、(6.25)式で $dt$ と $d\tau$ の関係が、それぞれ、与えられているので、ただちに、

\[ x = {\Large \int}_{0}^{\tau}{c \ \frac{U_{0} - \left(1 - \frac{\mu}{m_{0}} \tau \right)^{2\varOmega_{\mathrm{P}}/c}}{U_{0} + \left(1 - \frac{\mu}{m_{0}} \tau \right)^{2\varOmega_{\mathrm{P}}/c}} \cdot \frac{1}{2} \, \left\{ \sqrt{U_{0}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{-\varOmega_{\mathrm{P}}/c} + \frac{1}{\sqrt{U_{0}}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{\varOmega_{\mathrm{P}}/c} \ \right\} \, d\tau} \] \[ \quad = \frac{c}{2} \, {\large \int}_{0}^{\tau}{\left\{ \sqrt{U_{0}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{-\varOmega_{\mathrm{P}}/c} - \frac{1}{\sqrt{U_{0}}} \, \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{\varOmega_{\mathrm{P}}/c} \right\} \, d\tau} \tag{6.37} \]

と書き直すコトができます。この積分は、とても簡単です。結果だけを示すと、次のようになります。

\[ \begin{align} &\mbox{$\varOmega_{\mathrm{P}} \ne \pm c$ の場合:} \ x = \frac{c \, m_{0}}{2\mu} \left[ \ \frac{\sqrt{U_{0}}}{1 - \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{1 - \varOmega_{\mathrm{P}}/c} \ \right\} \right. \nonumber \\[0.5em] &\hspace{15em} \left. - \frac{1}{\sqrt{U_{0}}} \, \frac{1}{1 + \varOmega_{\mathrm{P}}/c} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{1 + \varOmega_{\mathrm{P}}/c} \ \right\} \ \right] \tag{6.38a} \end{align} \] \[ \mbox{$\varOmega_{\mathrm{P}} = c$ の場合:} \ x = \frac{c \, m_{0}}{2\mu} \left[ \ -\sqrt{U_{0}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \tau \right)} - \frac{1}{2\sqrt{U_{0}}} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \ \right\} \ \right] \] \[ {} \tag{6.38b} \] \[ \mbox{$\varOmega_{\mathrm{P}} = -c$ の場合:} \ x = \frac{c \, m_{0}}{2\mu} \left[ \ \frac{1}{\sqrt{U_{0}}} \, \ln{\left( 1 - \frac{\mu}{m_{0}} \tau \right)} + \frac{\sqrt{U_{0}}}{2} \, \left\{ \, 1- \left( 1 - \frac{\mu}{m_{0}} \tau \right)^{2} \ \right\} \ \right] \] \[ {} \tag{6.38c} \]

(6.38a)〜(6.38b)式を、グラフに表したのが次の図です(ただし、この図は、$\varOmega_{\mathrm{P}} > 0$、および、ロケットの初速度がゼロ、すなわち、$u_{\mathrm{R0}} = 0$($U_{0} = 1$)の場合を、描いているコトに、注意して下さい)。
図15:時刻とロケットの移動距離の関係:特殊相対論の世界(約67KB)


□このグラフには、(6.38a)〜(6.38b)式で与えられた、「ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、慣性系で測ったロケットの移動距離 $x$ の関係 $x = x(\tau)$(青線)」だけでなく、「慣性系にある時計の時刻 $t$ と、慣性系で測ったロケットの移動距離 $x$ の関係 $x = x(t)$(緑線)」も、描かれています。ちなみに、(6.38a)〜(6.38c)式で与えられる $x = x(\tau)$ と異なり、$x = x(t)$ のほうは、簡単なカタチで表すコトが、できそうにありません。上の $x = x(t)$ のグラフは、(6.27a)〜(6.27b)式の $t = t(\tau)$ に、$t$ の値を代入して、数値計算で、等式を満たす $\tau$ を求めたあとに、その $\tau$ を、(6.38a)〜(6.38b)式の $x(\tau)$ に代入する――という作業を(関数グラフ描画ソフトが自動的に)繰り返すコトで、描いています。なお、ココでも、「ロケットに乗っている観測者が持つ時計の時刻が $\tau = 0$ となるときに、慣性系にある時計の時刻が $t = 0$ になるように、お互いの時計の時刻を合わせておく」という条件(初期条件)が、採用されています((6.27a)〜(6.27c)式の説明を参照)。


■さて、グラフを見ると、慣性系で測ったロケットの移動距離 $x$ が、(どんな値でも良いですけど)ある値に達したとき、慣性系にある時計の時刻 $t$ よりも、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ のほうが、小さい($\tau < t$ である)コトに気づきます。コレが、いわゆる「ウラシマ効果」と呼ばれる現象です。具体的に説明してみましょう。いま、出発地(たとえば、地球とか)から、慣性系で測った距離で、距離 $x_{\mathrm{d}}$ だけ離れた目的地(たとえば、他の恒星系とか銀河とか)を目指して、ロケットが発進するとします(もちろん、出発地と目的地は、お互いに静止している――すなわち、同じ慣性系にあるとします)。また、ロケットの初速度はゼロとして、出発の瞬間に、慣性系にある時計の時刻 $t$ と、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ を、$t = \tau =0$ に合わせます(もちろん、ココで使う時計は、正確に、まったく同じ性質・性能であるとします)。ところが、先ほどの $x = x(t)$ と $x = x(\tau)$ のグラフによると、このロケットが目的地に到着したとき(慣性系で測ったロケットの移動距離が、$x_{\mathrm{d}}$ になったとき)には、目的地にある時計(慣性系にある時計)と、ロケットに乗っている観測者が持つ時計は、異なる時刻を指している――というコトになります。いま説明した状況を、分かりやすく示したのが、下の図です。
図16:ウラシマ効果:特殊相対論の世界(約18KB)


残念ながら、$t = t(x)$ や $\tau = \tau(x)$ を、簡単なカタチで表すコトは、できそうにありません。そのため、上の図における「ウラシマ効果による時間のズレ」を、具体的に計算するのは、ちょっと、タイヘンです。さらに、実際には、ロケットが出発してから、目的地に到着するまで、推進剤の噴射速度も消費速度も一定である――という設定は、不自然でしょう。現実的には、目的地に到着するときに、ロケットの速さがゼロになるように(目的地で停止するように)、推進剤の噴射速度や消費速度を、コントロールするハズです。そんなワケで、ロケットによる宇宙航行においては、「ウラシマ効果による時間のズレ」を、具体的に計算するコトは、かなり、骨の折れる作業になりそうです。ぜひ、自分で、いろいろな目的地や、その目的地を目指すための航行計画を考えて、実際に、ウラシマ効果の計算を楽しんでみて下さいね。


■ところで、先ほど示した図を見ると、ロケットに乗っている観測者が持つ時計の時刻 $\tau$ と、慣性系で測ったロケットの移動距離 $x$ の関係($x = x(\tau)$)を表したグラフの傾き(すなわち、$dx/d\tau$)が、時間の経過とともに大きくなり、無限大に向かって、発散していくように感じられます。実際に、(6.38a)〜(6.38c)式より、$dx/d\tau \to \pm \infty$($\tau \to m_{0}/\mu$)となるコトは、すぐに、確かめられます。「速さ=移動距離÷時間」であるので、コレは、ロケットの速さが、光速 $c$ を超えて、無限大へと発散してしまうコトを、示しているのでしょうか。もちろん、そうではありません。定義に戻って考えてみましょう。$x = x(\tau)$ のグラフの傾き $dx/d\tau$ は、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$ の定義($u_{\mathrm{R}} \equiv dx/dt$)、および、(3.2)式を用いて、

\[ \frac{dx}{d\tau} = \frac{dt}{d\tau} \cdot \frac{dx}{dt} = \frac{1}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}}} \cdot u_{\mathrm{R}} \tag{6.39} \]

と表されるコトになります。(6.39)式により、$u_{\mathrm{R}} \to c$ のときに、$dx/d\tau \to \infty$ となるコトは、明らかです。つまり、$dx/d\tau \to \infty$ となるコトと、慣性系にいる観測者から見たロケットの速さ $u_{\mathrm{R}}$(すなわち、ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R} \,}$)が、光速 $c$ を超えられないというコトのあいだには、ナニも、矛盾はないのです。


□実は、このコトには、いわゆる「ローレンツ収縮」と呼ばれる現象が、関係しています。ローレンツ収縮とは、お互いに静止している 2 点のあいだの距離を測るときに、その 2 点に対して静止している観測者が測定する距離 $\varDelta x$ よりも、その 2 点に対して運動している観測者が測定する距離 $\varDelta \xi$ のほうが、短くなる――という現象のコトです。よって、慣性系で測った距離で、等間隔に並べられた標識のあいだの距離 $\varDelta x$ は、ロケットに乗っている観測者から見ると、$\varDelta x$ よりも狭くなります。このため、コレらの標識のあいだ(距離 $\varDelta x$)を通過するのに要する時間(ロケットに乗っている観測者が持つ時計で測った時間)$\varDelta \tau$ も、短くなります。速さ $u_{\mathrm{R}}$ が増加すると、ロケットに乗っている観測者から見た標識のあいだの距離は、より狭くなり、通過時間 $\varDelta \tau$ は、速さ $u_{\mathrm{R}}$ の増加以上に短くなります。こうして、$\varDelta x/\varDelta \tau \to \infty$ となるのです。

□上の議論は、意外と奥が深いので、もう少しだけ、このハナシを続けましょう(詳しくは、『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』を、読んで下さい)。そもそものハナシでは、「加速度運動をしているロケットに乗っている観測者が測定する距離」を、正しく定義する手段はない――と、説明されていました(第 3 節を参照)。すると、「速さ」の定義(速さ=移動距離÷時間)により、「加速度運動をしているロケットに乗っている観測者が測定する速さ」も、正しく定義する手段はない――というコトになってしまうハズです。それならば、ココまでの議論で、たびたび登場してきた、「ロケットに乗っている観測者から見た慣性系の速さ $-u_{\mathrm{R}}(\tau)$」とは、厳密には、ナニを意味しているのでしょうか。

ココで、ロケットに乗っている観測者が持つ時計で測る、極めて短い時間 $\varDelta \tau$ に起こるコトについて、考えてみましょう。実際には、ロケットは加速度運動をしていますけど、微小時間 $\varDelta \tau$ は、あまりに短い時間なので、このあいだは、ロケットの速度変化が、ほとんどゼロであると、みなせるとします。つまり、微小時間 $\varDelta \tau$ のあいだは、ロケットは等速度運動をしている、すなわち、ロケットに乗っている観測者は慣性系にいる――と、みなすコトができるワケです。観測者が慣性系にいるなら、レーザー測距の原理で、(微小時間 $\varDelta \tau$ のあいだに、レーザーが往復できる距離までに限れば)微小距離 $\varDelta \xi$ が、正しく測定できるハズです。こうして、ロケットに乗っている観測者から見た速さ $\omega \equiv \varDelta \xi/\varDelta \tau$ が、定義されるのです。

こうして定義される、ロケットに乗っている観測者から見た速さ $\omega \equiv \varDelta \xi/\varDelta \tau$ は、もちろん、極めて短い時間間隔、極めて狭い距離空間のなかだけで、成り立つモノです。けれども、慣性系に並べられた標識を、すれ違う瞬間に観測して、すれ違う速さ(すなわち、ロケットに乗っている観測者から見た慣性系の速さ)を測るのには、十分です。さらに、ロケットに乗っている観測者が持つ時計で測った時間 $\varDelta \tau$ と、慣性系にある時計で測った時間 $\varDelta t$ のあいだには、(3.2)式で与えられる関係が、成り立っています。そして、慣性系で測った距離 $\varDelta x$ を、ロケットに乗っている観測者が測ったときの距離(慣性系にいるとみなせる観測者が、慣性系とみなせる距離空間のなかで測った距離)$\varDelta \xi$ は、ローレンツ収縮の効果によって、
\[ \qquad \varDelta \xi = \sqrt{1 - \left( \frac{u_{\mathrm{R}}}{c} \right)^{2}} \ \varDelta x \tag{6.40} \]
となるコトが知られています(詳しくは、多くの特殊相対論の教科書や解説などで説明されていますので、ぜひ、調べてみて下さいね)。したがって、慣性系に並べられた標識を、ロケットに乗っている観測者から見たときに、標識がすれ違う速さ(すなわち、ロケットに乗っている観測者から見た慣性系の速さ)$\omega$ は、次のように表されます(マイナス符号は、ロケットに乗っている観測者から見て、慣性系にある標識が $\varDelta \xi$ だけ移動したとき、慣性系にいる観測者から見て、ロケットは $-\varDelta x$ だけ移動するコトを、表します)。
\[ \qquad \omega \equiv \frac{\varDelta \xi}{\varDelta \tau} \ = \frac{-\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}} \ \varDelta x}{\sqrt{1 - ( u_{\mathrm{R}}/c )^{2}} \ \varDelta t} = -\frac{\varDelta x}{\varDelta t} \equiv -u_{\mathrm{R}} \tag{6.41} \]
(6.41)式は、慣性系にいる観測者から見たロケットの速さが $u_{\mathrm{R}}$ のとき、ロケットに乗っている観測者から見た慣性系の速さが $-u_{\mathrm{R}}$ となるコトを、示しています。つまり、「ロケットに乗っている観測者から見た慣性系の速さは、慣性系にいる観測者から見たロケットの速さと、大きさが同じで、向きが逆になる」という「相対性原理」が、この場合にも、確かに、成り立っている――というワケです。

□ところで、ココまでの説明のなかで、当たり前のコトであるかのように、「微小時間 $\varDelta \tau$ のあいだは、ロケットは等速度運動をしている、すなわち、ロケットに乗っている観測者は慣性系にいる――と、みなすコトができる」と述べていますけど、この説明は、どこまで、本当に正しいのでしょうか。実は、この疑問を、トコトン、突き詰めるコトこそが、「等価原理」から、「一般相対論」へと続く道を切り拓く、大きな手がかりの 1 つとなるのです。詳しいハナシは、『特殊相対論で、加速度運動。そして、ウラシマ効果 & 等価原理へ。』のなか(特に『おまけ 2: 「等価原理」から一般相対論へ』)に、書いてありますので、興味があれば、ぜひ、読んでみて下さいね。

■このページの一番上に戻る。


■Back to TOP-PAGE