スライド 1

Download Report

Transcript スライド 1

太陽・恒星フレアにおける輻射
流体シミュレーション
京大・宇宙物理
川道 俊見
概要
• イントロダクション
– 輻射流体の背景
– 現在の輻射流体シミュレーション
• 世界における輻射流体シミュレーション
– Ding さん(Nanjing University)グループ
– Carlsson さん(Oslo University)グループ
• 京大グループの今後(太陽系外フレア)
イントロダクション
• 輻射流体の背景
– フレアに伴う大気変動
• Machado (1980) など
– 非熱的電子の効果
• Ding & Fang (1997) など
• いずれもダイナミクス・エネルギーの取り扱い
は単純
Semi-empirical models
• Basic equations
radiative transfer:
statistical equilibrium:

dI 
d
 I  S
N
ni

N
j  1, j  i
hydrostatic equilibrium:
I  (   0 ,  ) 
 

h
0
  dh


0
S  (t  )e
n
( R ij  C ij ) 
j
(R
ji
 C ji )
j  1, j  i
mg  (1 . 1n H  n e )kT
 t / 
dt  / 
S 
   (n i  n j g i / g j )
2h 
c
e
2
3
1
nig
2
m ec
f
j
/n jgi 1
1
  D
H (a , p )
Semi-empirical models
• Method
MODEL
T (m ), v t (m )
compute
c
I
compare
adjust the model
with
c
I
o
I
Comparison of some models

Quiet-Sun:
VAL-C

Small flare:
F1

Big flare:
F2

White-light
flare:
F1*,F3
Flare models including a nonthermal
electron beam
radiative transfer:
statistical equilibrium:

dI 
d
 I  S
N
ni

N
B
( R ij  C ij  C
C ijBij ) 
j  1, j  i
hydrostatic equilibrium: mg  (1 .1n H  n e )kT

n j (R
j  1, j  i
ji
C
B
ji
C
C Bjiji )
Effects of the nonthermal electron beam and velocity
fileds on the H profile (Ding & Fang 1997)
Both blue and red asymmetries can be produced by downward
motions, depending on the strength of the electron beam
現在の輻射流体シミュレーション
• プラズマのダイナミクス・エネルギー収支の取
り扱い
– 非熱的電子の効果
– Radiative backwarming の効果
• より深い大気の加熱
– White light flares (WLFs) や Black light flares
(BLFs)、恒星フレアの再現
White light flare (太陽)
Liu et al. (2001)
M型星(AD Leo)
Haisch et al. (1991)
Black light flare (太陽)
Ding et al. (2003)
Ding et al. (2003)
• エネルギーバランスを考えた数値計算
– 非熱的電子が Coulomb 衝突で落とすエネル
ギーと放射冷却がつりあうように温度発展
– Collisional excitaion & ionization : H 4レベル
+連続
– Radiative cooling : H, H
• 結果
-
– H による Cooling で WLFs, BLFs を説明
– TMR (Temperature minimum region) の大気
Ding et al. (2003)
• 非熱的電子が落とすエネルギー
• Radiative cooling rate
• Energy balance
Continuum(λ=8500Å) の増し分
Continuum contrast
Flare maximum phase
(しだいに増光: WLFs)
Early impulsive
phase
(最初、減光: BLFs)
非熱的電子のエネルギーフラックス
TMR 付近の温度
1011 erg cm-2 s-1
1010 erg cm-2 s-1
初期大気
Allred et al. (2006)
• M 型星フレアにおける輻射流体
– ダイナミクス、エネルギーをきちんと解く
– XEUV backwarming
– H, He, CaII
• 結果
– 観測と合う Doppler shift と Line broadening
(Stark broadening)
– 一方で、Continuum の特徴は観測と合わず
基礎方程式(Abbett & Hawley 1999)
Cooling (H) et al. (2006)
Allred
Completely
Ionized (H)
Cooling (He)
Allred et al. (2006)
Completely
Ionized (He)
Allred et al. (2006)
Line broadening
Continuum dimming
Balmer & Paschen jump
京大グループの今後
• 可視光
– 太陽:エラーマンボム(中村くん)
– 恒星フレア:観測と協力(野上さん、栗山くん)
• 興味深いラインの挙動等が発見 (V711 Tau)
• 輻射流体シミュレーションで再現(川道)
• 電波とX線(川道)
– Benz 則
• 褐色矮星フレア
• 太陽系外惑星フレア
Benz 則
Absorption and Emission Coefficients
Absorption coefficient
 

ji
i


g
 n  i n  ( ) 
j
 i g
 ij
j


 (n
 ni e
*
i
 h  / kT
) ik ( )  n k n e kk ( )( 1  e
 h  / kT
)  n e e
i
Emission coefficient
j 
2h
c
2
3


 i

j i
 g 
n j  i  ij ( ) 
g 
 j
n
*
i
i
ik
( ) e
 h  / kT

  n k n e kk ( )  n e e J 
