連結階層シミュレーションアルゴリズムの開発

Similar documents
連結階層シミュレーションアルゴリズムの開発

07_Shiota.pptx

資料

連結階層シミュレーションアルゴリズムの開発

連結階層シミュレーションアルゴリズムの開発

s

名古屋大学宇宙地球環境研究所 ( 所長 : 町田忍 ) の塩田大幸 ( しおただいこう ) 特任助教と国立極地研究所の片岡龍峰 ( かたおかりゅうほう ) 准教授からなる研究チームは 太陽で起きる爆発の結果として巨大なプラズマ ( 注 1) と磁場が宇宙空間を通して地球に到達する過程について これま

1 a b cc b * 1 Helioseismology * * r/r r/r a 1.3 FTD 9 11 Ω B ϕ α B p FTD 2 b Ω * 1 r, θ, ϕ ϕ * 2 *

Spacecraft Propulsion Using Solar Energy Spacecraft with Magnetic Field Light from the Sun Solar Wind Thrust Mirror Solar Sail Thrust production by li

(MHD) ( ) MHD 2

パネルディスカッション (太陽研究長期計画)

PowerPoint プレゼンテーション

anan.pptx

PowerPoint Presentation

SFN

PowerPoint プレゼンテーション

Microsoft PowerPoint - okamura.ppt[読み取り専用]

宇宙天気予報と磁気流体力学

Gravothermal Catastrophe & Quasi-equilibrium Structure in N-body Systems

Annual Report of the Earth Simulator Center April March 2014 using the CReSS-SDM. The simulations are carried out on an density of aerosols. As

28 Horizontal angle correction using straight line detection in an equirectangular image

BH BH BH BH Typeset by FoilTEX 2

今後 UDAS に解析ソフトが追加される予定のデータの紹介 京都大学 飛騨天文台所有の太陽画像データ 京都大学大学院理学研究科附属天文台 飛騨天文台 上野悟 1

GJG160842_O.QXD

磯部委員提出資料

太陽活動領域における電流層の構造と上安定性

200708_LesHouches_02.ppt

スケーリング理論とはなにか? - --尺度を変えて見えること--

2012専門分科会_new_4.pptx

超新星残骸Cassiopeia a と 非球対称爆発

Solar Flare neutrino for Super Novae Conference

本講演の内容 1. はじめに 2. 天体 MHD 現象のおもしろさ 3. 天体 MHD シミュレーションは超困難 4. 天体 MHD シミュレーションの魅力 5. 天体 MHD シミュレーションの魔力 はまると危険 ( 落とし穴の数々 ) 6. むすび : ノーベル賞課題 ( 超難問 ) に挑戦せよ

5 時 10 分 22 秒 5 時 15 分 25 秒 5 時 20 分 27 秒 図 年 4 月 10 日のフレアの Hα 単色像 ( 京都大学飛騨天文台ドームレス太陽望遠鏡による ) フレアの典型的な大きさは 1 万km~10 万km 寿命は数分 ~ 数時間 解放される全エネルギー

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.

Akira MIZUTA(KEK) AM, Nagataki, Aoi (ApJ, , 2011) AM + (in prep)

観測的宇宙論WS2013.pptx

PSTEP サマースクール陸別 2017 シラバス 講義タイトル : 太陽活動現象 時 :2017 年 7 31 ( )10:30-12:00 講師 : 清 敏 ( 宇宙航空研究開発機構宇宙科学研究所 准教授 ) 連絡先メールアドレス :shimizu<at>solar.isas.jaxa.jp 概

B

磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論

18 2 F 12 r 2 r 1 (3) Coulomb km Coulomb M = kg F G = ( ) ( ) ( ) 2 = [N]. Coulomb

alternating current component and two transient components. Both transient components are direct currents at starting of the motor and are sinusoidal

q quark L left-handed lepton. λ Gell-Mann SU(3), a = 8 σ Pauli, i =, 2, 3 U() T a T i 2 Ỹ = 60 traceless tr Ỹ 2 = 2 notation. 2 off-diagonal matrices

2. 2 P M A 2 F = mmg AP AP 2 AP (G > : ) AP/ AP A P P j M j F = n j=1 mm j G AP j AP j 2 AP j 3 P ψ(p) j ψ(p j ) j (P j j ) A F = n j=1 mgψ(p j ) j AP

Studies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth

Microsoft PowerPoint - hoshino_part1.pptx

観測的宇宙論workshop.pptx

PowerPoint プレゼンテーション

プラズマ核融合学会誌11月【81‐11】/小特集5

OHO.dvi

untitled

G (n) (x 1, x 2,..., x n ) = 1 Dφe is φ(x 1 )φ(x 2 ) φ(x n ) (5) N N = Dφe is (6) G (n) (generating functional) 1 Z[J] d 4 x 1 d 4 x n G (n) (x 1, x 2

( )


Optical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd:

July 28, H H 0 H int = H H 0 H int = H int (x)d 3 x Schrödinger Picture Ψ(t) S =e iht Ψ H O S Heisenberg Picture Ψ H O H (t) =e iht O S e i


研究成果報告書

4/15 No.

発表資料

Microsoft PowerPoint - Fukazawa pptx

太陽からの電波

untitled

Kerr 時空における球対称流に対するコリメーション効果 ( CQG, 26, , 2009 ) 髙見健太郎 ( 広島大学 / Albert-Einstein-Institute) 共同研究者 : 小嶌康史 ( 広島大学 ) 2009 年 10 月 01 日駒場宇宙コロキウム

スライド 1

Microsoft PowerPoint - 島田美帆.ppt

Hanbury-Brown Twiss (ver. 2.0) van Cittert - Zernike mutual coherence

- 3 Solar-B -

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1)

Vol.-ICS-6 No.3 /3/8 Input.8.6 y.4 Fig....5 receptive field x 3 w x y Machband w(x =

OPA277/2277/4277 (2000.1)

chap1_MDpotentials.ppt

日立金属技報 Vol.34

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

余剰次元のモデルとLHC

Microsoft PowerPoint Aug30-Sept1基研研究会熱場の量子論.ppt

untitled

September 9, 2002 ( ) [1] K. Hukushima and Y. Iba, cond-mat/ [2] H. Takayama and K. Hukushima, cond-mat/020


iBookBob:Users:bob:Documents:CurrentData:flMŠÍ…e…L…X…g:Statistics.dvi

08 p Boltzmann I P ( ) principle of equal probability P ( ) g ( )g ( 0 ) (4 89) (4 88) eq II 0 g ( 0 ) 0 eq Taylor eq (4 90) g P ( ) g ( ) g ( 0

nsg02-13/ky045059301600033210

0946: : : :30 UT 1004:30 UT global expansion global expansion pseudo breakup pseudo breakup 3 1: Kotzebue UT [Sh

輻射の量子論、選択則、禁制線、許容線

tsunoda _1

A5 PDF.pwd


SiC SiC QMAS(Quantum MAterials Simulator) VASP(Vienna Ab-initio Simulation Package) SiC 3C, 4H, 6H-SiC EV VASP VASP 3C, 4H, 6H-SiC (0001) (11 20) (1 1

ニュートン重力理論.pptx

PowerPoint Presentation

A

(77) F Moreton wave EIT wave, pp (78) K ( ) ( ) X, pp (79) K Tongjiang Wang,Jia-long Wang,Yihua Yan,H.Kurokawa,and K.Shibata The Large-S

Microsoft PowerPoint - 太陽黒点と磁場.pptx

建築構造力学 I ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.

MHD) MHD MHD 1977 MHD ApJ MHD simulation Abstract/Keywords ADS simulation ApJ MHD simulation 1982

Introduction SFT Tachyon condensation in SFT SFT ( ) at 1 / 38

JKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n =

PowerPoint プレゼンテーション

1 ( ) Einstein, Robertson-Walker metric, R µν R 2 g µν + Λg µν = 8πG c 4 T µν, (1) ( ds 2 = c 2 dt 2 + a(t) 2 dr 2 ) + 1 Kr 2 r2 dω 2, (2) (ȧ ) 2 H 2

J-PARC E15 K K-pp Missing mass Invariant mass K - 3 He Formation K - pp cluster neutron Mode to decay charged particles p Λ π - Decay p Decay E15 dete

Transcription:

2013 年天文夏の学校 2013.7.30 太陽フレアはいつ起きるか? 草野完也 名古屋大学太陽地球環境研究所

太陽フレア フレアの際の soft X-ray 流束の時間変化 2013/7/30

3

太陽嵐の影響と被害 宇宙天気短期的な宇宙環境変動 太陽フレア コロナ質量放出 (CME) 磁気嵐 電離圏嵐 環境 社会への甚大な影響 Coronal Mass Ejection 太陽フレア爆発の予測 宇宙プラズマ爆発現象の理解 衛星通信システム被曝影響送電網 パイプライン Quebec blackout A single serious event might cause $2 trillion in damages. (NRC)

磁気リコネクション フレア モデル Shiota et al. 2005 Flare ribbons corresponds to the feet of reconnected field lines. 2013/7/30

フレア爆発のトリガ問題いつ どこで なぜ どうやって発生するか? 光球面 彩層 Hinode/SOT

爆発的 ( 突発 ) 現象 地震 雪崩 噴火 株暴落 革命 破局 http://en.wikipedia.org/wiki/black_monday_%281987%29 2013/7/30

爆発的現象の特徴 自由エネルギー の蓄積と解放 (storage-and-release processes) stable 自由エネルギー の蓄積 自由エネルギー の解放 unstable 不安定性 ( フィードバック ) low price 例 : 大暴落のメカニズム fear of investor sell order 2013/7/30

ヘリシティの過飽和による不安定化 磁気ヘリシティ H = A B dv double mode state single mode state H=0 ( ポテンシャル磁場 ) H>0 single mode state 不安定 double mode state 安定 (Kusano et al. 1995, Kusano & Nishikawa 1996)

ひので (2006~) 2013/7/30

Nonlinear Force-Free Field Model Inoue, Kusano et al. ApJ 2011 B = αb Force-free field 太陽表面観測磁場を境界条件とした実データフォースフリー磁場モデル Twist T n = (4π ) 1 field line [ B] B z z dl Smaller than the critical

フレア発生にとって重要な条件は? Large-scale Component 磁気シア ( 磁気ヘリシティ ) 自由エネルギー Savcheva et al. 2012 + 構造の複雑さ? 不安定化 ( トリガ ) Zirin and Wang 1993 Small-scale Component 問題 : どのような複雑構造がフレアのトリガとなるか?

フレア発生条件として磁場構造 Long PIL of strong magnetic shear (Hagyard, et al. 1984) sigmoidal structure (Rust & Kumar 1996; Canfield et al. 1999) flux rope (Forbes & Priest 1995; Torok & Kliem 2005) magnetic helicity injection (Kusano et al. 1995, 1996) double loop structure (Hanaoka 1997) flux cancellation (van Ballegooijen & Martens 1989) converging foot point motion (Inhester et al. 1992) narrow magnetic lanes between major sunspots (Zirin & Wang 1993) emerging magnetic fluxes (Heyvaerts, Priest & Rust 1977; Moore & Roumeliotis 1992; Feynman & Martin 1995; Chen & Shibata 2000) the sharp gradient of magnetic field (Schrijver 2007) reversed magnetic shear (Kusano et al. 2004) topological complexity (Schmieder et al. 1994) intermittency and multifractality (Abramenko & Yurchyshyn 2010) multipolar topologies (Antiochos et al. 1999)

太陽面爆発のトリガとなる磁場構造は何か? 複雑すぎてデータを見ているだけでは何が重要な構造なのか分からない 人間は注目したいものしか注目しない それ故 シミュレーションを利用すべき

シミュレーションによる 2 つの戦略 戦略 1 Data-driven Simulation 戦略 2 Ensemble Simulation 様々な初期境界条件 パラメタ 初期境界条件パラメタ 仮説提示 観測による検証

Parameters in Ensemble Simulation Large field (free energy) & Small field (trigger) azimuth: ϕ e shear angle: θ 0 Box:Rectangle including PIL Initial condition:lfff offset δ flux Φ 161 cases 3D MHD 256x1024x512 grids output: 800 GB/run Earth Simulator (JAMSTEC)

Parameter Space: θ 0 vs. ϕ e strong shear weak shear potential field Right Polarity Opposite Polarity Right Polarity Normal Shear Reversed Shear

シミュレーション結果 strong shear no flare flare no flare weak shear potential field Right Polarity Opposite Polarity Normal Shear Reversed Shear

Flare Phase Diagram strong shear no flare OP flare RS no flare weak shear potential field Right Polarity Opposite Polarity Normal Shear Reversed Shear

反極性型 (OP-type)

逆磁場型 (RS-type)

Magnetic structure just prior to flare OP type Bp current sheet chromosphere Flux tube is formed by reconnection in chromosphere.

Magnetic structure just prior to flare OP type flux tube Bp current sheet chromosphere Flux tube is formed by reconnection in chromosphere.

F=r c I 2 B p I トーラス不安定性 Ideal MHD F 電流 I リコネクション

2 つのフレアトリガ シナリオ Eruption-induced Reconnection OP type MHD Instability Magnetic Reconnection RS type Reconnection-induced Eruption

観測的検証 OP RS 2006-12-13 02:22 UT RS OP OP RS Zirin and Wang 1993 Kurokawa, Wang & Ishii 2002 Green, Kliem & Wallace 2011

フレアトリガ領域の形成過程 Bamba, et al 2013 (submitted) Toriumi, Iida et al 2013 (in press) magnetic flux トリガ磁場の磁束 Ca line emission プリフレア発光 (Ca II H) flare

非線形不安定性の分岐遷移 Critical Phenomena Triggered Phenomena How large free energy has been loaded? free energy When will he come? Linear instability Nonlinear instability

まとめ 太陽フレアは大きなスケール ( 数 10Mm 以上 ) におけるエネルギー蓄積と小さなスケール ( 数 Mm 以下 ) の擾乱の相互作用によって発生する非線形不安定化現象である 逆極性 (OP) 型磁場と逆シア型 (RS) 型磁場は 効果的なフレアトリガとして働く それゆえ 大きなスケールと小さなスケール双方の磁場観測によって フレア発生予測は原理的に可能である

News: Solar Flares are Predictable! 朝日新聞 2012 年 11 月 3 日

予測と近代科学 ハレー彗星 Wikipedia より エドモンド ハレー (1656 年 10 月 29 日 - 1742 年 1 月 14 日 ) 1682 年に出現した彗星の観測データとニュートン力学から この彗星が 76 年の周期を持つ楕円軌道を持つと結論 ( プリンキピア出版は1687 年 ) 過去の記録から 1531 年 1607 年に出現した彗星が同一のものと推測 次回の回帰が1758 年であると予測 1758 年 12 月 25 日 予測通り彗星が出現