Intel Xeon Phi ∗ ł ¯.´. ˚ 3, ˝ 1,
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org ∂ρ ∂t + ∇ · (ρ~u) = S − D, ∂∂ρtH + ∇ · (ρH~u) = −sH,H2 + S ρρH − DH ρρH , ∂ρH2 + ∇ · (ρH2~u) = sH,H2 + S ρ − D ρ
ρH2 ρH2 , ∂t ∂ρ~u ∂t + ∇ · (ρ~u~u) = −∇p − ρ∇(Φ) + ~vS − ~uD, ∂ρE
+ ∇ · (ρE~u) = −∇ · (p~v) − (ρ∇(Φ), ~u) − Λ + Γ − ε D , ∂t ρ ∂ρε
+ ∇ · (ρǫ~u) = −(γ − 1)ρǫ∇ · ~u − Λ + Γ − ε D , ∂t ρ ρE = 1 ρ~u2 + ρε,
2 p = (γ − 1)ρε, ∂n ∂t + ∇ · (n~v) = D − S, ∂n~v ∂t + ∇ · (n~v~v) = −∇Π − n∇(Φ) + ~uD − ~vS, ∂ρW
+ ∇ · (ρW ~v) = −∇ · (Π~v) − (n∇(Φ), ~v) − Γ + ε D , ∂t ρ ∂Πξξ + ∇ · (Πξξ~v) = −2Π∇ · ~u − Γ + ε D , ∂t 3ρ ρW = 1 ρ~v2 + Πxx + Πyy + Πzz , 2 2 dnH2 = Rgr(T )nH (nH + 2nH2) − (ξH + ξdiss(NH2, AV )) nH2,
dt ßØ Ł º, º ae ε ß, Ł ª Πξξ = (Πxx, Πyy, Πzz) ae Æ aeae ºŒ Ł º Æ aeae ºŒ Ł º ØŒ ß, ρE º Œ º ß, Æ S º º Ø ı Ł aeŒ Ø Œ Ł β = 0.1
M SN
T < 104K ▽ ·~u < 0 ρ > 1.64 pMc−⊙3 S = ΔSN = ddρt = βC τdnyn = βCn3/2s 332πG Γ = 1051 M SN erg
M⊙ Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org ρ3/2 + ρ3/2 n3/2 + n3/2 ρ = √LρL + √ρR n = √LnL + √nR
R R p = pL√ρL + pR√ρR Π = ΠL√nL + ΠR√nR
√ρL + √ρR √nL + √nR U = uL(−λt) + uR(λt) + pL(−λt) − pR(λt)uv (√ρL + √ρR)2
u 2 2 t γL+2γR (ρ3L/2 + ρ3R/2)(pL√ρL + pR√ρR) P = pL(−λt) + pR(λt) + uL(−λt) − uR(λt)uuv γL+γR (ρ3L/2 + ρ3R/2)(pL√ρL + pR√ρR), 2 2 t 2 (√ρL + √ρR)2 V = vL(−µt ) + vR(µt ) + ΠL(−µt ) − ΠR(µt )uuv (√nL + √nR)2 2 2 t3(n3L/2 + n3R/2)(ΠL√nL + ΠR√nR)
v ΠL(−µt ) + ΠR(µt ) + vL(−µt ) − vR(µt )uu3(n3L/2 + n3R/2)(ΠL√nL + ΠR√nR) Π = 2 2 t (√nL + √nR)2 λ = utvu γL+2γR (ρpL3L/√2ρ+Lρ+3R/p2R√ρR) µ = utuv 3(ΠL√n3Ln/2L++nΠ3RR/2√nR) qL(−νt) = qiR − 2νht △qi − qi6 1 − 23νht qR(νt) = qiL + 2νht △qi + qi6 1 − 23νht Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org |~v| = q2(E − ǫ), (E − ~v2/2)/E ≥ 10−3
ρ~v2 ! |ρǫ| = ρE − 2 , (E − ~v2/2)/E < 10−3
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org
I, K, L
~iπij + ~iπkm + ~iπln ! ρi,k,l + ni,k,l = jXmn σjmnexp I K L
,~i φjmn = Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org (80%)
(6%) #pragma omp parallel for ... {} ...
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org P = T otal1
T otalK Sratio = Sreal
Speak 140 120 100 -pU 80 ed 60 e pS 40 20 0 0,8 y it ilb0,6 a l ca0,4 S 0,2 0,0 1 2 4 16 32
64
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org -2-2 -2-2 -2-2 -2-2 T = T otal1
T otalp
0 x 10 kpc c kp 0 0 1 x -1 0 100 200 300 400 500 0 2 4 6 8 10 P. 1017-1024. intera tions // Monthly Noti es of the Royal Astronomi al So iety. 2007. V. 375, I. 3. 4. Sol Alonso M., Lambas D., Tissera P., Coldwell G. A tive gala ti nu lei and galaxy 2011. V. 55, I. 9. P. 770-783. 1. Tutukov A., Lazareva G., Kulikov I. Gas Dynami s of a Central Collision of Two Galaxies: Merger, Disruption, Passage, and the Formation of a New Galaxy // Astronomy Reports.
So iety. 2003. V. 339, I. 2. P. 289-311. 10. Springel V., Hernquist L. Cosmologi al smoothed parti le hydrodynami s simulations: a hybrid multiphase model for star formation // Monthly Noti es of the Royal Astronomi al
Astrophysi al Journal Supplement Series. 1996. V. 105. P. 19-35.
behind Sho k Waves // The Astrophysi al Journal. 2004. V. 612. P. 921-939. frontiers and innovations. 2014. V. 1, I. 3. P. 88-98. of Parallel Numeri al Methods for Plasma Physi s and Astrophysi s // Super omputing 5. Glinskiy B., Kulikov I., Snytnikov A., Romanenko A., Chernykh I., Vshivkov V. Co-design Simulations of Intera ting Galaxies // The Astrophysi al Journal Supplements Series. 2014. V. 214, I. 12. P. 1-12. 7. Kulikov I. GPUPEGAS: A New GPU-a elerated Hydrodynami Code for Numeri al 2005. V. 329. P. 143-152. 3. S hweizer F. Merger-Indu ed Starbursts // Astrophysi s and Spa e S ien e Library.
and Spa e S ien e. 2002. V. 281, I. 1-2. P. 383-387. Alternative to N-body Methods // Monthly Noti es of the Royal Astronomi al So iety. 2013. V. 428. P. 2674-2687. 6. Mit hell N., Vorobyov E., Hensler G. Collisionless Stellar Hydrodynami s as an E ien t IEEE. 2005. V. 93, I. 2. P. 216-231. 19. Frigo M., Johnson S. The Design and Implementation of FFTW3 // Pro eedings of the 20. Kalinkin A., Laevsky Y., Gololobov S. 2D Fast Poisson Solver for High-Performan e Computing // Le ture Notes in Computer S ien e. 2009. V. 5698. P. 112-120. simulations // Computational Mathemati s and Mathemati al Physi s. 2007. V. 47, I. 15. Popov M., Ustyugov S. Pie ewise paraboli method on lo al sten il for gasdynami 12. P. 1970-1989. 11. Sutherland R., Dopita M. Cooling fun tions for low-density astrophysi al plasmas // The Astrophysi al Journal Supplement Series. 1993. V. 88. P. 253-327. 12. Roe P. Approximate Riemann solvers, parameter ve tors, and dieren e solvers // Journal of Computational Physi s. 1997. V. 135. P. 250-258.
Mathemati al Physi s. 2014. V. 54, I. 6. P. 1012-1024.
Equations with Entropy Nonde rease Guarantee // Computational Mathemati s and 18. Godunov S., Kulikov I. Computation of Dis ontinuous Solutions of Fluid Dynami s simulations // Journal of Computational Physi s. 1984. V. 54. P. 174-201. 13. Collela P., Woodward P. The Pie ewise Paraboli Method (PPM) Gas-Dynami al operator-splitting s heme for numeri al simulations of hydrodynami o ws // New 14. Kulikov I., Vorobyov E. Using the PPML approa h for onstru ting a low-dissipation, Astronomy. 2015. (in progress) Igor Kulikov, Igor Chernykh, Vladislav Nenashev and Eugenya Katysheva