11 iteration ALI Accelerated Lambda Iteration Variable Eddington

11 数値解法 反復解法が必要 ・Λ-iteration ・ALI (Accelerated Lambda Iteration) ・Variable Eddington Factor 法 ・Feautrier法 ・Rybicki法 ・Complete Linearization 法 ・. . . ・Monte Carlo法 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

13 数学的構造 Feautrier 法 (Rybicki法, Complete Linearization法, . . . ) Λ-iteration法 Accelerated Lambda Iteration (ALI) (Operator Perturbation Method) 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

20 Comoving Frame O(v/c) O(1) の式を得る static diffusion dynanic diffusion free streaming relativistic flow すべての場合で無視できるもの以外は残す Lab Frame 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

22 解法 1. Flux Limited Diffusion (FLD) 法 2. Simple Method 3. Variable Eddington Factor (VEF) 法 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

1. Flux Limited Diffusion (FLD) 法 23 Diffusion 近似 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

26 3. Variable Eddington Factor (VEF) 法 Eddington Factor 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

36 原始星の構造推定 Ø 2 -D axisymmetric Ø three components; (Nakazato, TN, & Umemura 2003) Observer i Envelope • Central star: Luminosity Temperature Mass • Circumstellar disk: Surface density at 1 AU Power law index • Envelope: Density at 1 AU Circumstellar disk AU 0 0 10 100 AU Power law index Central star Opening angle ØTemperature: Radiative equilibrium Outflow region 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

37 parameters Model SED νLν [erg s-1] 1036 観測角度依存性 大 i の推定 原始星構造の推定 (νLν)max 1034 1032 1030 1028 1011 1012 1013 1014 1015 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

38 Ø IRAS 04365+2535（TMC 1 A） νLν [erg s-1] Class I 天体 (Nakazato, TN, Umemura 2003) 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

3 D 輻射輸送 (TN, Umemura, & Susa 2001) N 3＝ 1283 3, 1282 in (8 Mpc) Nangle = Isotropic background UV: I 21 z =0. 1 Zel’dovich approximation: = 15 41 Radiative Transfer Ionization Equilibrium Neutral Fraction: CP-PACS 256 PU (理論ピーク ~80 GFlops) ~ 100 hour 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

42 Reionization History of an Inhomogeneous Universe I 21=0. 1 Z=15 Z=9 Z=7 Z=5 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

44 Shadowing Effect (1) 3 D RT Raises Neutral Fraction in Optically Thick Region (2) Inhomogeneity Raises Neutral Fraction in All Region Shadowing Effect is prominent at (TN, Umemura, & Susa 2001) 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

45 3 D RHD: Cosmic Reionization (Gnedin 2000) XHI z=9 Cosmological HD ~1010 -12 Msun RT (Local Optical Depth 近似) Star Formation H, He H 2 form/dest ngas T 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

Optically Thin Variable Eddington Tensor 46 (Gnedin & Abel 2001) XHI LOD OTVET 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

47 3 D RHD (SPH): 銀河形成 (Susa & Umemura, in prep. ) Initial Condition by GRAFIC in COSMICS + Rigid rotaion with λ=0. 05 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史

50 Heterogeneous Multi Computer System GRAPE 6 8 boards~8 Tflops PCI bus ~200 MB/s × 8 Alpha Cluster (16×alpha 21264) CP‐PACS 2048 PU~600 Gflops (peak) Internal Communication 1： 1→ 280 MB/s 16 parallel 40 MB/s (peak) 第 15回 理論懇シンポジウム 国立天文台 2002年 12月26日 中本泰史