軸対称-北海道大学.ppt

３次元 MHD Nested Grid シミュレーションによる星間磁気雲の重力収縮 町田正博(北海道大学)　富阪幸治(国立天文台)　　　松本倫明(法政大人間環境) この研究での目的 以下の研究をする事を今回の目的とした 星間磁気雲から星形成までの過程 等温磁気雲の成長過程 軸対称(m=0)、非軸対称(m=2)揺らぎによる成長の違い 中心密度の時間進化 非軸対称の場合の軸の成長 将来、アウトフロ－現象、星間雲の分裂などを研究するための足がかり 研究手法 3次元 MHD nested grid シミュレーション 流体部分 数値流速補正?マルチタイムステップ　–Roeの方法 時間積分　– fractional timestep 空間２次精度　MUSCL法 勾配制限関数　–TVD 等温?断熱?ポリとロープ 自己重力 Multigrid法 Nested grid 最細グリッドでJeans conditionが破れると新しいグリッドを生成 初期条件 グリッド数 ： 64*64*64 Scale : λ/2 (λ:most unstable wavelength) Unit : r0c=1 　tff=1 X-Y平面の密度、磁場分布 : cylinder の平衡解 磁場、回転の強さ : a=B2/8prcs2=1.0 w=0.1 Perturbation case A： z方向にz=0を中心とする密度揺らぎ case B： z方向の密度揺らぎ + X-Y平面にm=2の密度揺らぎ 終了条件 : 15段目の格子のJeans 条件が破れた時点で計算終了 まとめ?今後の課題 平衡状態にあるCylinder状の等温ガス雲に摂動（m=0,m=2,z軸上）を与えた場合の収縮過程についてのシミュレーションを行った Z軸対称(m=0)、Z軸非対称の場合(m=2)で収縮に違いを求める事が出来た 軸対称、非軸対称共に自己相似的に成長する事が確かめられた 今後、星形成機構の解明のため、アウトフローの発生機構を研究するためにガス雲が光学的に厚くなった後(n1010)のシミュレーションをする必要がある * * L=1 L=2 L=3 Nested grid Z X Y 回転 磁場の向き λ/2 0 　　t = 0.0 ρc = 0.0 　　t = 4.79 ρc = 1.69 　　t = 6.36 ρc = 3.35 　　t = 7.16 ρc = 7.93 　　t = 7.29 ρc = 10.8 　　t = 7.65 ρc = 39.5 ★Y-Z方向の進化（横から見た図） 計算結果 L = 12,13,14 ρc = 2.2×108 L = 9,10,11 ρc = 5.2×106 L = 7,8,9 ρc = 1.8×105 L = 12,13,14 ρc = 2.4×106 L = 9,10,11 ρc = 5.2×106 L = 6,7,8 ρc = 7.2×105 ★X-Y方向の進化（横から見た図） [上段 m=0(軸対称)　、下段m=2(非軸対称)] ★m=0モード（軸対称）の計算の最終段階 L=13,14,15t=8.02rc=1.91x109fig.a:y-z図、密度コントア、速度ベクトル、磁力線fig.b:x-y図、密度コントア、速度ベクトル、密度　　　等高線fig.c:y-z図、磁場コントア、磁力線、磁場等高線 Fig.a Fig.b Fig.c ★m=2モード（非軸対称）の計算の最終段階 L=13,14,15t=6.83rc =2.01x109fig.a:y-z図、密度コントア、速度ベクトル、磁力線fig.b:x-y図、密度コントア、速度ベクトル、密度　　　等高線fig.c:y-z図、磁場コントア、磁力線、磁場等高線 Fig.a Fig.b Fig.c ★3次元で見た図(最終段階、15段目)　① 軸対称 非軸対称 磁力線 密度コントアと速度ベクトル Z=0面 ★3次元で見た図(最終段階、15段目)　② 軸対称 非軸対称 Y=0面 ★3次元で見た図(最終段階、15段目)　③ ほぼ真上(Z軸)からみた図 軸対称 非軸対称 軸の成長　① 上図：m=0の軸の成長下図：m=2の軸の成長 y軸－軸の長さ x軸－密度 軸の定義：ρ=0.1ρmaxの範囲にあるガスλ:Jeans length (λ=2p/kj　　kj=(4pGr/cs2)) Z軸の長さ XY平面の軸の長さ(長軸、短軸) 　軸対称 非軸対称 Jeans length Iij=∫r r max rirj r(r) dr M= ∫r 0.1 r max