运动方程式と正则化-东海大学.pptVIP

* ヘリウム-３輸送用， 月-地球周回軌道の計算 Calculation of a periodical orbit　 between the earth and the moon for 3He transportation 九州東海大学大学院 工学研究科　生産工学 御手洗研究室　 増田　孝弘 １．まえがき ? 　中性子も発生せず，放射性物質であるトリチウム（T）も 使用しないクリーンな発電方式であるD-3He核融合発電は人 類にとって究極のエネルギーである． 　この発電炉の燃料である重水素（D）は海水中に無尽蔵に あるが，ヘリウム-３（3He）は地球上にほとんど無い．しか し，月から持ち帰られたレゴリスと呼ばれる小さな砕石の中 に太陽の核融合で作られ，太陽風によって運ばれた3Heが多量に蓄積されている． 　 　本研究ではこの3Heを月から地球へ効率的かつ低コストで輸送する方法を考案するためのもので，少ない燃料で月?地球間を往復できるような輸送機の周回軌道を数値計算によって見い出し，その軌道特性を明らかにした． ２．運動方程式と正則化 ２．１．輸送機（ペイロード）の運動方程式 　地球?月周回軌道を見出すために，地球と月の重心を中心として月の 公転周期Ωで回転する座標系を考える．（いわゆる円制限３体問題） 地球?月重心系回転座標 月の地球と 月重心からの座標 （X，Y） ペイロードの地球と月重心からの座標 （Xe，Ye ） 地球の地球と 月重心からの座標 （Xm，Ym） 地球?月重心系回転座標系におけるペイロードの運動方程式 ただし, 　Gは万有引力定数，m，Me，Mmはペイロード，地球，月の質量． （１），（２）式を４次のルンゲ?クッタ?ギル法で解く． （１） （２） 二つの天体が大接近すると重力が大きくなり特異点に近くなるため，一般に数値解はその誤差が急に大きくなりその後の軌道を正確に計算することが困難になる．よって，ペイロードが地球や月に接近した時には計算に誤差が生じないようにレビ?チビタの変換を利用し，この特異点を取り除くいわゆる正則化を行わなければならない． ［１］　レビ?チビタ変換とは 　このレビ?チビタ変換は数学的には２つの複素平面において，Z平面から W平面に等角写像することに等しい． より、 ２．２．レビ?チビタ変換による正則化 レビ?チビタ変換 ［２］正則化方程式の求め方 ペイロード,地球，月の位置と速度 ペイロードと地球，月の距離 （４） （３） ペイロードが地球に近づく場合 （X，Y）空間から（U，V）空間へのレビ?チビタ変換 （５） （７） （６） （８） （９） 　（10） （11） （12） （13） （15） （14） 以上から ペイロードの全エネルギー これらから特異点のない（いわゆる正則化された）運動方程式が得られる． （16） （ U，V）空間の正則化した運動方程式 （18） ペイロードと地球の距離がある距離以上に離れた場合 　　逆に，（u，v），（u’，v’）の最後の値を　　　，　　　に変換し初期値とし，（X，Y）空間における運動方程式を用いて軌道を計算する． 図の表示のために（X，Y）空間に逐次変換する． （19） （20） ペイロードと地球の距離がある距離に近づいた場合，それまで計算した，　　　，　　　の最後の値を初期値として，(u，v),(u’，v’)に変換し，正則化運動方程式を解く． 　３．正則化計算プログラムのチェック 　 ３．１．簡単化のため月の回転速度Ωを0とした場合 　　　　　　　　ペイロードの軌道が急激に変化する場合 　　　　　　（A）正則化なし　　　計算が不安定になる 　　　初期条件：座標X = ‐386,238 km，Y = 0 km，初速度2.3034 km/s （A） 正則化なし ＊　この初期条件ではペイロードが地球に近づいた時に計算が不安定　　になり，右上方に急速に飛び去ってしまう． ＊　これは２つの物体が大接近して起こる特異点のためである． （B）正則化あり　　　計算が安定になる. 同じ初期条件でペイロードと天体の距離が30,000 km 以下で正則化を行った結果，月まで戻ってくる軌道を計算することができる．全運動エネルギーが一定で，軌道が正しく計算されていることがわかる． （B） 正則化あり ３．２．月が公転している実際の場合 （月の公転速度Ω＝2.649×10-6 km/s） 　地球上の地点：X =‐17,122 km，Y =‐21,316 kmで，初速度：VX = 4.7245 km/s，VY =‐2.7597 km/sで月をめがけてペイロードを打ち上げる． ＊正則化をしていない場合