Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence
|
|
- Jemima Hines
- 5 years ago
- Views:
Transcription
1 Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence Junseong Park, Dongsu Ryu Dept. of Physics, Ulsan National Institute of Science and Technology. Ulsan, Korea 2016 KNAG meeting KASI, Daejeon, Korea
2 Purpose of Study Smith et al studied the distribution of Mach number in driven supersonic HD and MHD turbulence. However, they did not consider fast and slow Mach numbers distribution in molecular cloud conditions. In this study we will identify differences between the two kinds of shocks, because density and velocity structures in turbulence are associated with the fast and slow shocks. 1. The main goal is to identify the formulas for fast and slow Mach numbers and energy dissipation rate at shocks and understand the characteristic properties of fast and slow MHD shocks. 2. We studied MHD turbulence in a variety of astrophysical environments by performing numerical simulations of isothermal, compressible MHD turbulence with resolution up to Contents MHD equations of compressible, isothermal gas - simulation initial conditions Formula for Fast and Slow Mach numbers - fast and slow Mach number D structures - fast and slow Mach number distribution Formula for Energy dissipation rate at shocks - fast and slow energy dissipation rate Conclusion and Future work
3 MHD equations of compressible, isothermal gas ρ + ρv = 0 t v t + v v + c 2 s ρ ρ 1 ρ B v B = 0 t B B = 0 CONTINUITY MOMENTUM INDUCTION with the additional constraint for the absence of magnetic monopoles, B = 0 Simulation initial conditions Isothermal TVD Code (Kim, et al. 1999) Turbulence was driven by solenoidal ( v = 0) forcing (Stone et al. 1999; Mac Low 1999) Initial conditions: - computational domain L 0 = L x = L y = L z = 10 with periodic boundaries - uniform density ρ 0, magnetic field B 0 = 2/β, 0,0 where β = P g /P B = 2c s 2 /V A 2 = 0.1, 1, 10 - velocity forcing drawn from Gaussian random field, P k k 6 exp( 8k/k peak ), where k peak = 2k 0 (k 0 = 2π/L 0 ) - M s = v rms /c s = 0.5, 1,, 8 at saturation with time intervals t = 0.001L 0 /c s - flow length scale L d = 2π/k peak = L 0 /2, dynamical time t d = L d /v rms = L 0 / 2M s c s - resolution: 128, 256, grid cells - self-gravity was ignored.
4 MHD equations of compressible, isothermal gas ρ + ρv = 0 t v t + v v + c 2 s ρ ρ 1 ρ B v B = 0 t B B = 0 CONTINUITY MOMENTUM INDUCTION with the additional constraint for the absence of magnetic monopoles, B = 0 Simulation initial conditions Saturation was reached around ~2t d Comparison of observation & simulation : - cloud clump of size L 0 = 2pc, T = 10K, c s 0.2km/s, n H2 = 10 cm (observed MCs) - sound crossing time: t s = L 0 /c s ~10Myr, dynamical time t d = L d /v rms = L 0 / 2M s c s = 5Myr/M s - B 0 = 1.4μGβ 1/2 T/10K 1/2 n H2 /10 2 cm 1/2 ~ 14, 4.4, 1.4 μg for β = 0.1, 1, 10
5 Comparison of sliced density images (n = 512, t/t d = 9) M s = 0. 5 (subsonic) M s = 1 (transonic) M s = (supersonic) M s = 8 (supersonic) β = 0. 1 (sub-alfvénic) β = 1 (trans-alfvénic) B 0 β = 10 (super-alfvénic) Filaments and sheets with high density are formed in a flow with higher M s Supersonic turbulence Dense core formation
6 Formula for Fast and Slow Mach numbers Shock jump conditions of both sides of the shock front are as follows: ρ 1 v 1 = ρ 2 v 2 ρ 1 v c 2 s ρ B 1 2 = ρ 2 v c 2 s ρ B 2 2 (1) (2) plane-parallel shock pre-shock Post-shock ρ 1, v 1, v 1 ρ 2, v 2, v 2 B, B 1, c s B, B 2, c s ρ 1 v 1 v 1 B B 1 = ρ 2 v 2 v 2 B B 2 v 1 B 1 v 1 B = v 2 B 2 v 2 B () (4) n sh : parallel to shock normal : perpendicular to shock normal B, c s : same in pre-shock and Post-shock χ = ρ 2 /ρ 1 = v 1 /v 2 > 1 preferred frame for shock surface B 1 /B = v 1 /v 1, B 2 /B = v 2 /v 2 = v 2 χ/v 1 From 2, equations, we can rewrite as follows: ρ 1 v c 2 s ρ B 1 2 = ρ 2 1v 1 χ + χc s 2 ρ B 2 2 ρ 1 v 2 1 B 2 B 1 = ρ 1v 1 χ B 2 B 2 2 (5) (6) B 2 substituted into the equation 5, then 2ρ 1 v 6 1 4B B 1 ρ c 2 s ρ 1 χ + B 2 1 ρ 2 1 χ 2 v B 2 c 2 s ρ B 2 B B 4 ρ 1 χ 2 v 2 1 2B 4 c 2 s ρ 1 χ = 0
7 Formula for Fast and Slow Mach numbers Shock jump conditions of both sides of the shock front are as follows: ρ 1 v 1 = ρ 2 v 2 ρ 1 v c 2 s ρ B 1 2 = ρ 2 v c 2 s ρ B 2 2 ρ 1 v 1 v 1 B B 1 = ρ 2 v 2 v 2 B B 2 v 1 B 1 v 1 B = v 2 B 2 v 2 B (1) (2) () (4) divide by 2ρ 1 c s 6, and set M s = v 1 /c s, c A = B / ρ 1, c A = B 1 / ρ 1, α = c A /c s, α = c A /c s M s α α 2 χ α 2 χ 2 M s 4 + α α 2 + α 2 χ 2 M s 2 α 4 χ = 0 (7) the phase velocity of fast and slow modes is derived from equation 7 with χ = 1 as follows: 2 c fa,sl 2 = 1 c s α 2 + α 2 ± 1 + α α 2 2 4α Fast and Slow Mach numbers M fa = M s c s /c fa, M sl = M s c s /c sl (8)
8 Formula for Fast and Slow Mach numbers How to find M s from 6th order equation? let x = M s 2 x + C a x 2 + C b x 2 + C c = 0 where C a = 1 + 2α α 2 χ α 2 χ 2 C b = α α 2 + α 2 χ 2 C c = α 4 χ Compute Q = C a 2 C b /9, R = 2C a 9C a C b + 27C C /54 If Q 2 < R, then the cubic equation has three real roots (Press et al. 1986) as follows: x 1 = 2 Qcos θ C a x 2 = 2 Qcos θ + 2π C a x = 2 Qcos θ 2π C a where θ = arccos R/ Q Otherwise Q 2 > R, the cubic equation has only one real root. Compute A = sgn R R + R 2 Q 1, B = Q/A A 0 or 0 A = 0 and calculate the real root as follows: x 1 = A + B a/
9 Formula for Fast and Slow Mach numbers How to find M s from 6th order equation? let x = M s 2 x + C a x 2 + C b x 2 + C c = 0 where C a = 1 + 2α α 2 χ α 2 χ 2 C b = α α 2 + α 2 χ 2 C c = α 4 χ Compute Q = C a 2 C b /9, R = 2C a 9C a C b + 27C C /54 If Q 2 < R, then the cubic equation has three real roots (Press et al. 1986) as follows: x 1 = 2 Qcos θ slow mode C a x 2 = 2 Qcos θ + 2π fast mode C a x = 2 Qcos θ 2π C a where θ = arccos R/ Q Otherwise Q 2 > R, the cubic equation has only one real root. Compute A = sgn R R + R 2 Q 1, B = Q/A A 0 or 0 A = 0 and calculate the real root as follows: x 1 = A + B a/ fast mode 0 slow mode
10 Find shocks - Shock cells: ρ i+1 ρ i 1 ρ i+1 +ρ i 1 = χ 1 χ+1 > (i.e. χ = ) and v < 0 - Shock center: v is local minimum and labeled this center as part of a shock surface. - Shock quantities: max(ρ i+2, ρ i 2 ) is post-shock and min(ρ i+2, ρ i 2 ) is pre-shock Find Fast and Slow shocks - Fast shock B 2 /B 1 1, M fa 1.0, χ Slow shock 0 < B 2 /B 1 < 1, 1 < M sl < c A /c sl, 1 < χ < 2α 2 + α α 2 + α α 2 /4 ρ 1 v c 2 s ρ B 1 2 = ρ 2 1v 1 χ + χc s 2 ρ B 2 2 ρ 1 v 2 1 B 2 B 1 = ρ 2 1v 1 χ B 2 B 2 (5) (6) with B 2 = 0
11 Fast Mach number D structures (n = 512, t/t d = 9) M s = 0. 5 (subsonic) M s = 1 (transonic) M s = (supersonic) M s = 8 (supersonic) β = 0. 1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) x y B 0 z The locations of identified fast shocks along with the gas density, velocity and magnetic field distributions
12 Fast Mach number D structures (n = 512, t/t d = 9) M s = 0. 5 (subsonic) M s = 1 (transonic) M s = (supersonic) M s = 8 (supersonic) β = 0. 1 (sub-alfvénic) anisotropic β = 1 (trans-alfvénic) β = 10 (super-alfvénic) isotropic z x y B 0 The turbulence remains isotropic for M s = 8 and β = 10, but when it becomes M s = 0.5 and β = 0.1, the turbulence becomes highly anisotropic.
13 Slow Mach number D structures (n = 512, t/t d = 9) M s = 0. 5 (subsonic) M s = 1 (transonic) M s = (supersonic) M s = 8 (supersonic) β = 0. 1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) x y B 0 z The locations of identified slow shocks along with the gas density, velocity and magnetic field distributions
14 Slow Mach number D structures (n = 512, t/t d = 9) M s = 0. 5 (subsonic) M s = 1 (transonic) M s = (supersonic) M s = 8 (supersonic) β = 0. 1 (sub-alfvénic) anisotropic β = 1 (trans-alfvénic) β = 10 (super-alfvénic) isotropic z x y B 0 The turbulence remains isotropic for M s = 8 and β = 10, but when it becomes M s = 0.5 and β = 0.1, the turbulence becomes highly anisotropic.
15 Fast Mach number distribution with the number of shocks (t/t d = 4~9) β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) The start in the power law is found to occur at M fa min given by M fa min ~2.5M s 0.2 /n 0.1 dn = na b M dm fa M fa where a~2.7 ± 0.2 s
16 Fast Mach number distribution with the number of shocks (t/t d = 4~9) β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) The start in the power law is found to occur at M fa min given by M fa min ~2.5M s 0.2 /n 0.1 dn = na b M dm fa M fa where a~2.7 ± 0.2 s Power-law index b with M s b~m s 0.6 M s b As M s increases, b becomes smaller, but it becomes larger with increasing resolution Power law index is not significantly affected by plasma beta, and it follows 0.6 a power law b~m s
17 Slow Mach number distribution with the number of shocks (t/t d = 4~9) β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) The start in the power law is found to occur at M sl min given by M sl min ~2.5M s 0.2 /n 0.1 dn = na b M dm sl M sl where a~2.7 ± 0.2 s
18 Slow Mach number distribution with the number of shocks (t/t d = 4~9) β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) The start in the power law is found to occur at M sl min given by M sl min ~2.5M s 0.2 /n 0.1 dn = na b M dm sl M sl where a~2.7 ± 0.2 s Power-law index b with M s b~m s 0.49 M s b As M s increases, b becomes smaller, but it becomes larger with increasing resolution Power law index is not significantly affected by plasma beta, and it follows 0.49 a power law b~m s This relationship work very well for super-alfvénic case, but breaks down at sub/trans-alfvénic cases with M s =, 8
19 Slow Mach number distribution with the number of shocks (t/t d = 4~9) β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) The start in the power law is found to occur at M sl min given by M sl min ~2.5M s 0.2 /n 0.1 dn = na b M dm sl M sl where a~2.7 ± 0.2 s Upstream Alfvénic Mach number M A (= ρ 1 v 1 /B 1 ) with M s fast shock slow shock When the turbulence becomes sub-alfvénic or trans-alfvénic, M A 0.2, anisotropies arise in the gas, because the back reaction of the magnetic field onto the flow is extremely strong for flows perpendicular to the magnetic field lines (Cho & Lazarian, 200; Esquivel & Lazarian, 2011) magnetic field and anisotropic become increasingly important
20 Formula for energy dissipation rate at shock When we adopt an isothermal equation of state (Mouschovias et al. 1974) that mimics the effect of radiative losses, the total energy is not conserved if shock waves are present in the system. E = 1 2 ρv2 + PlnP B2 dv From the conservation equation of energy, we obtain t 1 2 ρv2 + PlnP B ρv2 + PlnP + P v + B v B = 0 where v = v e + v e, B = B e + B e t 1 2 ρv2 + PlnP B2 + x 1 2 ρv2 + PlnP + P + B 2 v B B v + B v = 0 Some of the energy is dissipated at shocks with the isothermal conditions.
21 Formula for energy dissipation rate at shock When we adopt an isothermal equation of state (Mouschovias et al. 1974) that mimics the effect of radiative losses, the total energy is not conserved if shock waves are present in the system. E = 1 2 ρv2 + PlnP B2 dv From the conservation equation of energy, we obtain t 1 2 ρv2 + PlnP B ρv2 + PlnP + P v + B v B = 0 where v = v e + v e, B = B e + B e t 1 2 ρv2 + PlnP B2 + 1 x 2 ρv2 + PlnP + P + B 2 v B B v + B v = 0 The energy dissipation rate Q at shock discontinuity Q = 1 2 ρv2 + PlnP + P + B 2 v B B v + B v ρv2 + PlnP + P + B 2 v B B v + B v 2
22 Formula for energy dissipation rate at shock The energy dissipation rate Q at shock discontinuity Q = 1 2 ρv2 + PlnP + P + B 2 v B B v + B v ρv2 + PlnP + P + B 2 v B B v + B v 2 In the preferred frame for shock surface B /B = v /v then B 2 v B B v + B v = 0 and with ρ 1 v 1 = ρ 2 v 2 or P 1 v 1 = P 2 v 2 energy dissipation rate Q becomes Q = 1 2 ρ 1v ρ 1 lnρ 1 c s 2 v ρ 2v ρ 2 lnρ 2 c s 2 v 2 After some arithmetic manipulation Q = M s ρ 1 c 1 s 2 M s χ M s 2 α 2 2 α 1 M s 2 2 α 2 M 2 s α 2 χ 2 lnχ,for M s 1, most of the upstream flow energy is dissipated Q ~ M s
23 Fast and Slow Energy dissipation rate Q at shock (t/t d = 4~9) Q = 1 2 : 1 2 ρv2 + PlnP + P + B 2 v B B v + B v β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) ~M fa Fast shock ~M sl Slow shock Q fa,sl increases with fast and slow Mach numbers, asymptotically follows a power law Q fa,sl ~M fa,sl The energy is dissipated by a small range of stronger shocks
24 Fast and Slow Energy dissipation rate Q at shock (t/t d = 4~9) Q = 1 2 : 1 2 ρv2 + PlnP + P + B 2 v B B v + B v β = 0.1 (sub-alfvénic) β = 1 (trans-alfvénic) β = 10 (super-alfvénic) ~M fa Fast shock ~M sl Slow shock How much percent of the turbulence energy is dissipated due to fast and slow shocks?
25 Conclusion We studied compressible isothermal MHD turbulence in a variety of M s and β, with resolution up to 1024 The power law for fast and slow Mach numbers distribution with the number of shocks - dn = na b M dm fa,sl M fa,sl s where a~2.7 ± as M s increases, b becomes smaller, but it becomes larger with increasing resolution Fast and Slow Energy dissipation rate at shock - as the plasma beta β increases, Q fa,sl becomes smaller. - Q fa,sl increases with fast and slow Mach numbers, asymptotically follows a power law Q fa,sl ~M fa,sl - energy is dissipated by a small range of stronger shocks Future work No physical dissipation terms are modelled, but numerical diffusion is unavoidable. In saturation Energy injection = Energy dissipation at shocks + Energy dissipation due to viscosity & resistivity how to quantify this value is the next project!!
Intracluster Shock Waves
ρ Intracluster Shock Waves T R 200 Dongsu Ryu (UNIST, Ulsan National Institute of Science and Technology, Korea) Hyesung Kang (Pusan National U, Korea) Sungwook E. Hong (Korea Institute for Advanced Study)
More informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
More informationThe Magnetic Field of GMCs. Paolo Padoan (ICREA & Institute of Cosmos Sciences - University of Barcelona)
The Magnetic Field of GMCs Paolo Padoan (ICREA & Institute of Cosmos Sciences - University of Barcelona) Two different views of the magnetic field in MCs: 1. The old view (Shu et al. 1987) Strong mean
More informationMagnetic Fields (and Turbulence) in Galaxy Clusters
Magnetic Fields (and Turbulence) in Galaxy Clusters Dongsu Ryu (UNIST, Ulsan National Institute of Science and Technology, Korea) with Hyesung Kang (Pusan Nat. U, Korea), Jungyeon Cho (Chungnam Nat. U.),
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationCompressible MHD Turbulence in Strongly Radiating Molecular Clouds in Astrophysics*
Poster T25 Compressible MHD Turbulence in Strongly Radiating Molecular Clouds in Astrophysics* D.D. Ryutov LLNL, Livermore, CA 94551, USA Presented at 8 th International Workshop on the Physics of Compressible
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationDirect Evidence for Two Fluid Effects in Molecular Clouds. Dinshaw Balsara & David Tilley University of Notre Dame
Direct Evidence for Two Fluid Effects in Molecular Clouds Dinshaw Balsara & David Tilley University of Notre Dame 1 Outline Introduction earliest stages of star formation Theoretical background Magnetically
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationMixing in Highly Compressible Turbulence
Mixing in Highly Compressible Turbulence Liubin Pan Evan Scannapieco School of Earth and Space Exploration Arizona State University Astrophysical motivation: Heavy elements from supernova and stellar winds
More informationAlexei Kritsuk. CalSpace IGPP/UCR Conference Palm Springs, March 27-30, 2006 c Alexei Kritsuk
Supersonic Turbulence in Molecular Clouds 1 Alexei Kritsuk University of California, San Diego In collaboration with Rick Wagner, Mike Norman, & Paolo Padoan, UCSD http://akpc.ucsd.edu/isothermal CalSpace
More informationASTRONOMY AND ASTROPHYSICS The distribution of shock waves in driven supersonic turbulence
Astron. Astrophys. 362, 333 341 (2000) ASTRONOMY AND ASTROPHYSICS The distribution of shock waves in driven supersonic turbulence M.D. Smith 1, M.-M. Mac Low 2, and F. Heitsch 3 1 Armagh Observatory, College
More informationThe density variance Mach number relation in supersonic turbulence I. Isothermal, magnetized gas
Mon. Not. R. Astron. Soc. 423, 2680 2689 2012 doi:10.1111/j.1365-2966.2012.21075.x The density variance Mach number relation in supersonic turbulence I. Isothermal, magnetized gas F. Z. Molina, 1 S. C.
More informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More information26. Non-linear effects in plasma
Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless
More informationShock Waves. = 0 (momentum conservation)
PH27: Aug-Dec 2003 Shock Waves A shock wave is a surface of discontinuity moving through a medium at a speed larger than the speed of sound upstream. The change in the fluid properties upon passing the
More informationFUNDAMENTALS OF MAGNETOHYDRODYNAMICS (MHD)
FUNDAMENTALS OF MAGNETOHYDRODYNAMICS (MHD) Dana-Camelia Talpeanu KU Leuven, Royal Observatory of Belgium Basic SIDC seminar ROB, 7 March 2018 CONTENTS 1. Ideal MHD 2. Ideal MHD equations (nooooooo.) 2.1
More informationCosmological Shocks and Their Signatures
Cosmological Shocks and Their Signatures Dongsu Ryu (Chungnam National U, Korea) Hyesung Kang (PNU, Korea), Renyi Ma (CNU, Korea), Jungyeon Cho (CNU, Korea), David Porter (U of Minnesota), T. W. Jones
More informationTurbulence and Magnetic Field in High-β Plasma of Intracluster Medium
Turbulence and Magnetic Field in High-β Plasma of Intracluster Medium Coma cluster of galaxies Dongsu Ryu UNIST (Ulsan National Institute of Science and Technology, Korea) KASI (Korea Astronomy and Space
More informationTHE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK
THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK G.P. Zank, I. Kryukov, N. Pogorelov, S. Borovikov, Dastgeer Shaikh, and X. Ao CSPAR, University of Alabama in Huntsville Heliospheric observations
More informationGrowth of Magnetic Fields by Turbulent Motions and Characteristic Scales of MHD Turbulence
Growth of Magnetic Fields by Turbulent Motions and Characteristic Scales of MHD Turbulence Jungyeon Cho Chungnam National University, Korea Cho, Vishniac, Beresnyak, Lazarian & Ryu 2009, ApJ Cho & Ryu,
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationVII. Hydrodynamic theory of stellar winds
VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the
More informationAmplification of magnetic fields in core collapse
Amplification of magnetic fields in core collapse Miguel Àngel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz Universitat de València; Max-Planck-Institut
More informationThe Janus Face of Turbulent Pressure
CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Christoph Federrath, Monash University Patrick Hennebelle, CEA/Saclay Alexei Kritsuk, UCSD yt-project.org Seminar über Astrophysik,
More informationExploring Astrophysical Magnetohydrodynamics Using High-power Laser Facilities
Exploring Astrophysical Magnetohydrodynamics Using High-power Laser Facilities Mario Manuel Einstein Fellows Symposium Harvard-Smithsonian Center for Astrophysics October 28 th, 2014 Ø Collimation and
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationSolutions to Merav Opher (2010) Problems
Solutions to Merav Opher 00 Problems. The normal of the shock is Since from the plot you can obtain all the three components of Bu and Bd, the normal can be easily found. The shock speed is: With the shock
More informationICMs and the IPM: Birds of a Feather?
ICMs and the IPM: Birds of a Feather? Tom Jones University of Minnesota 11 November, 2014 KAW8: Astrophysics of High-Beta Plasma in the Universe 1 Outline: ICM plasma is the dominant baryon component in
More informationTurbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields
Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Jungyeon Cho Dmitry Shapovalov MWMF Madison, Wisconsin April 2012 The Large Scale Dynamo The accumulation of magnetic
More informationMagnetospheric Physics - Final Exam - Solutions 05/07/2008
Magnetospheric Physics - Final Exam - Solutions 5/7/8. Dipole magnetic field a Assume the magnetic field of the Earth to be dipolar. Consider a flux tube with a small quadratic cross section in the equatorial
More information! what determines the initial mass function of stars (IMF)? dn /dm
! what determines the initial mass function of stars (IMF)? dn /dm M ! what determines the initial mass function of stars (IMF)?! What determines the total mass of stars that can form in the cloud? dn
More informationMHD Flow Field and Momentum Transfer Process of Magneto-Plasma Sail
J. Plasma Fusion Res. SERIES, Vol. 8 (2009) MHD Flow Field and Momentum Transfer Process of Magneto-Plasma Sail Hiroyuki NISHIDA, Ikkoh FUNAKI, Yoshifumi INATANI 1) and Kanya KUSANO 2) University of Tokyo,
More informationShock Waves: II. HII Regions + Planetary Nebulae
Dejan Urošević Department of Astronomy, Faculty of Mathematics, University of Belgrade Shock Waves: II. HII Regions + Planetary Nebulae Ionization Fronts and Shocks in HII Regions two stages of evolution
More information7 Stellar Winds & Supernovae Remnants
36 Physics 426 Notes Spring 2009 7 Stellar Winds & Supernovae Remnants Now let s look at some supersonic flow situations which involve shocks. You remember that we worked with smooth transonic flows last
More informationA Multi-ion Model of the Heliosphere with Secondary Charge Exchange
A Multi-ion Model of the Heliosphere with Secondary Charge Exchange Matthew Bedford, University of Alabama in Huntsville, Department of Space Science Nikolai Pogorelov, faculty advisor The heliosphere
More informationarxiv: v1 [astro-ph.ga] 19 Mar 2010
DRAFT VERSION OCTOBER 29, 2018 Preprint typeset using LATEX style emulateapj v. 11/10/09 arxiv:1003.3697v1 [astro-ph.ga] 19 Mar 2010 VELOCITY FIELD OF COMPRESSIBLE MHD TURBULENCE: WAVELET DECOMPOSITION
More informationThe Turbulent WIM: Distribution and MHD Simulations
The Turbulent WIM: Distribution and MHD Simulations Alex S. Hill University of Wisconsin-Madison Collaborators: Robert A. Benjamin, Grzegorz Kowal, Ronald J. Reynolds, L. Matthew Haffner, Alex Lazarian
More informationLES Simulations of Quiet Sun Magnetism
LES Simulations of Quiet Sun Magnetism Matthias Rempel HAO/NCAR Quiet sun magnetism Origin and spatial distribution of quiet sun field Small scale dynamo? Remnant field from large scale dynamo? Vögler,
More informationLinear and non-linear evolution of the gyroresonance instability in Cosmic Rays
Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationDoes magnetic-field-angular-momentum misalignment strengthens or weakens magnetic braking? (Is magnetic braking dynamically important?
Does magnetic-field-angular-momentum misalignment strengthens or weakens magnetic braking? (Is magnetic braking dynamically important?) Yusuke Tsukamoto Kagoshima University S. Okuzumi, K. Iwasaki, M.
More informationarxiv: v2 [astro-ph.he] 4 Mar 2013
DRAFT VERSION MARCH 5, 2013 Preprint typeset using L A TEX style emulateapj v. 5/2/11 UNIVERSALITY AND INTERMITTENCY IN RELATIVISTIC TURBULENT FLOWS OF A HOT PLASMA DAVID RADICE Max-Planck-Institut für
More informationCosmic-ray Acceleration and Current-Driven Instabilities
Cosmic-ray Acceleration and Current-Driven Instabilities B. Reville Max-Planck-Institut für Kernphysik, Heidelberg Sep 17 2009, KITP J.G. Kirk, P. Duffy, S.O Sullivan, Y. Ohira, F. Takahara Outline Analysis
More informationShocks. Ellen Zweibel.
Shocks Ellen Zweibel zweibel@astro.wisc.edu Departments of Astronomy & Physics University of Wisconsin, Madison and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas Shocks
More informationIntroduction and Basic Concepts
Topic 1 Introduction and Basic Concepts 1 Flow Past a Circular Cylinder Re = 10,000 and Mach approximately zero Mach = 0.45 Mach = 0.64 Pictures are from An Album of Fluid Motion by Van Dyke Flow Past
More informationMHD SIMULATIONS IN PLASMA PHYSICS
MHD SIMULATIONS IN PLASMA PHYSICS P. Jelínek 1,2, M. Bárta 3 1 University of South Bohemia, Department of Physics, Jeronýmova 10, 371 15 České Budějovice 2 Charles University, Faculty of Mathematics and
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationSimulating magnetic fields within large scale structure an the propagation of UHECRs
Simulating magnetic fields within large scale structure an the propagation of UHECRs Klaus Dolag ( ) Max-Planck-Institut für Astrophysik ( ) Introduction Evolution of the structures in the Universe t =
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationDissipation Scales & Small Scale Structure
Dissipation Scales & Small Scale Structure Ellen Zweibel zweibel@astro.wisc.edu Departments of Astronomy & Physics University of Wisconsin, Madison and Center for Magnetic Self-Organization in Laboratory
More informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
More informationarxiv: v4 [astro-ph.im] 27 Apr 2011
A ROBUST NUMERICAL SCHEME FOR HIGHLY COMPRESSIBLE MAGNETOHYDRODYNAMICS: NONLINEAR STABILITY, IMPLEMENTATION AND TESTS K. WAAGAN, C. FEDERRATH, AND C. KLINGENBERG arxiv:1101.3007v4 [astro-ph.im] 27 Apr
More informationTurbulence, kinematics & galaxy structure in star formation in dwarfs. Mordecai-Mark Mac Low Department of Astrophysics
Turbulence, kinematics & galaxy structure in star formation in dwarfs Mordecai-Mark Mac Low Department of Astrophysics Outline Turbulence inhibits star formation, but slowly Interplay between turbulence
More informationUNIVERSITY of CALIFORNIA SANTA CRUZ
UNIVERSITY of CALIFORNIA SANTA CRUZ EVOLUTION OF BLISTER-TYPE HII REGIONS IN A MAGNETIZED MEDIUM A thesis submitted in partial satisfaction of the requirements for the degree of BACHELOR OF SCIENCE in
More informationStructure Formation and Particle Mixing in a Shear Flow Boundary Layer
Structure Formation and Particle Mixing in a Shear Flow Boundary Layer Matthew Palotti palotti@astro.wisc.edu University of Wisconsin Center for Magnetic Self Organization Ellen Zweibel University of Wisconsin
More informationarxiv: v1 [astro-ph] 1 Jan 2008
Studies of regular and random magnetic fields in the ISM: statistics of polarization vectors and the Chandrasekhar-Fermi technique Diego Falceta-Gonçalves 1 arxiv:0801.0279v1 [astro-ph] 1 Jan 2008 Núcleo
More informationDiffusive shock acceleration: a first order Fermi process. jan.-fév NPAC, rayons cosmiques E. Parizot (APC)
1 Diffusive shock acceleration: a first order Fermi process 2 Shock waves Discontinuity in physical parameters shock front n 2, p 2, T 2 n 1, p 1, T 1 v 2 v 1 downstream medium (immaterial surface) upstream
More informationA Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows 4
draft of August 0, 998 A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows 4 Dongsu Ryu, Francesco Miniati,T.W.Jones, and Adam Frank 3 ABSTRACT A description is given for preserving
More informationA shocking group! Turbulent dissipation and star formation in Stephan s Quintet
Physics of Groups, IAP, Dec 12-15 A shocking group! Turbulent dissipation and star formation in Stephan s Quintet Pierre Guillard, P. Appleton, F. Boulanger, M. Cluver, P. Lesaffre, G. Pineau des Forêts,
More informationarxiv:astro-ph/ v1 14 Feb 2000
A&A manuscript no. (will be inserted by hand later) Your thesaurus codes are: 08 (02.08.1; 02.19.1; 02.20.1; 09.03.1; 09.11.1 ) ASTRONOMY AND ASTROPHYSICS 1.2.2008 arxiv:astro-ph/0002291v1 14 Feb 2000
More informationSimulations of Supersonic Turbulence in Molecular Clouds: Evidence for a New Universality
NUMERICAL MODELING OF SPACE PLASMA FLOWS: ASTRONUM-28 ASP Conference Series, Vol. 46, 29 Nikolai V. Pogorelov, Edouard Audit, Phillip Colella, and Gary P. Zank, eds. Simulations of Supersonic Turbulence
More informationarxiv: v2 [astro-ph.ga] 6 May 2015
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 26 February 2018 (MN LATEX style file v2.2) The density variance Mach number relation in isothermal and non-isothermal adiabatic turbulence arxiv:1504.04370v2
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationShocks in Heliophysics. Dana Longcope LWS Heliophysics Summer School & Montana State University
Shocks in Heliophysics Dana Longcope LWS Heliophysics Summer School & Montana State University Shocks in Heliophysics Vasyliunas, v.1 ch.10 Shocks in Heliophysics Tsuneta and Naito 1998 Forbes, T.G., and
More informationp = nkt p ~! " !" /!t + # "u = 0 Assumptions for MHD Fluid picture !du/dt = nq(e + uxb) " #p + other Newton s 2nd law Maxwell s equations
Intro to MHD Newton s 2nd law Maxwell s equations Plasmas as fluids Role of magnetic field and MHD Ideal MHD What do we need to know to understand the sun, solar wind&shocks, magnetospheres? Some material
More informationVELOCITY ANISOTROPY IN SELF-GRAVITATING MOLECULAR CLOUDS. I: SIMULATION
Draft version January 1, 217 Preprint typeset using L A TEX style AASTeX6 v. 1. VELOCITY ANISOTROPY IN SELF-GRAVITATING MOLECULAR CLOUDS. I: SIMULATION Frank Otto Department of Physics, The Chinese University
More information1 Energy dissipation in astrophysical plasmas
1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationSimulating supersonic turbulence in magnetized molecular clouds
Journal of Physics: Conference Series Simulating supersonic turbulence in magnetized molecular clouds To cite this article: Alexei G Kritsuk et al 9 J. Phys.: Conf. Ser. 8 View the article online for updates
More informationCollapse of Low-Mass Protostellar Cores: Part I
Collapse of Low-Mass Protostellar Cores: Part I Isothermal Unmagnetized Solutions and Observational Diagnostics Andrea Kulier AST 541 October 9, 2012 Outline Models of Isothermal Unmagnetized Collapse
More informationarxiv:astro-ph/ v2 23 Oct 2001
Kolmogorov Burgers Model for Star Forming Turbulence Stanislav Boldyrev 1 Institute for Theoretical Physics, Santa Barbara, CA 93106 arxiv:astro-ph/0108300v2 23 Oct 2001 ABSTRACT The process of star formation
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Shocks occur in supersonic flows; Shocks are sudden jumps in velocity, density and pressure; Shocks satisfy flux in = flux out principle
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationShocks in the ICM and the IPM
Shocks in the ICM and the IPM Tom Jones (University of Minnesota) 1 Outline Setting the stage for following talks The Interplanetary and Intracluster Media as Collisionless Plasmas Basic Introduction to
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Shocks occur in supersonic flows; Shocks are sudden jumps in velocity, density and pressure; Shocks satisfy flux in = flux out principle
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Central concepts: Phase velocity: velocity with which surfaces of constant phase move Group velocity: velocity with which slow
More information13. ASTROPHYSICAL GAS DYNAMICS AND MHD Hydrodynamics
1 13. ASTROPHYSICAL GAS DYNAMICS AND MHD 13.1. Hydrodynamics Astrophysical fluids are complex, with a number of different components: neutral atoms and molecules, ions, dust grains (often charged), and
More informationSuperdiffusive and subdiffusive transport of energetic particles in astrophysical plasmas: numerical simulations and experimental evidence
Superdiffusive and subdiffusive transport of energetic particles in astrophysical plasmas: numerical simulations and experimental evidence Gaetano Zimbardo S. Perri, P. Pommois, and P. Veltri Universita
More informationarxiv:astro-ph/ v1 14 Jan 2002
Accretion onto a Black Hole in the Presence of Bremsstrahlung Cooling - Parameter Space Study Sivakumar G. Manickam and Sandip K. Chakrabarti arxiv:astro-ph/020203v 4 Jan 2002 S. N. Bose National Centre
More informationTurbulence simulations with ENZO and FLASH3
Turbulence simulations with ENZO and FLASH3... not yet Adaptive-mesh simulations with FLASH Christoph Federrath Institute for Theoretical Astrophysics Heidelberg Collaborators: Ralf Klessen, Robi Banerjee,
More informationTurbulence and Magnetic Fields in the Large Scale Structure of the Universe
arxiv:0805.2466v1 [astro-ph] 16 May 2008 Turbulence and Magnetic Fields in the Large Scale Structure of the Universe Dongsu Ryu, 1 Hyesung Kang, 2 Jungyeon Cho, 1 and Santabrata Das 3 1 Department of Astronomy
More informationA super-alfvénic model of dark clouds
A super-alfvénic model of dark clouds Paolo Padoan Instituto Nacional de Astrofísica, Óptica y Electrónica, Apartado Postal 216, 72000 Puebla, México Åke Nordlund Astronomical Observatory and Theoretical
More informationDynamo on the OMEGA laser and kinetic problems of proton radiography
Dynamo on the OMEGA laser and kinetic problems of proton radiography Archie Bott, Alex Rigby, Tom White, Petros Tzeferacos, Don Lamb, Gianluca Gregori, Alex Schekochihin (and many others ) 2 th August
More information3D Reconnection of Weakly Stochastic Magnetic Field and its Implications
3D Reconnection of Weakly Stochastic Magnetic Field and its Implications Alex Lazarian Astronomy Department and Center for Magnetic Self- Organization in Astrophysical and Laboratory Plasmas Collaboration:
More informationSpace Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering
Space Physics ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen Aalto University School of Electrical Engineering The 6 th week: topics Last week: Examples of waves MHD: Examples
More informationCosmic-ray acceleration by compressive plasma fluctuations in supernova shells
Cosmic-ray acceleration by compressive plasma fluctuations in supernova shells Ming Zhang Florida Institute of Technology E-mail: mzhang@fit.edu We suggest that the production of Galactic cosmic rays in
More informationStudies of dissipative standing shock waves around black holes
Mon. Not. R. Astron. Soc. 401, 2053 2058 (2010) doi:10.1111/j.1365-2966.2009.15793.x Studies of dissipative sting shock waves around black holes Santabrata Das, 1,2 Sip K. Chakrabarti 3,4 Soumen Mondal
More informationCollapse of Massive Cloud Cores
Collapse of Massive Cloud Cores Robi Banerjee ITA, University of Heidelberg Based on 3D MHD, AMR* Simulations * Adaptive Mesh Refinement Collapse of Hydrostatic Cores Cloud formation by TI Slowly rotating
More informationThree-Dimensional Evolution of the Parker Instability under a Uniform Gravity
Three-Dimensional Evolution of the Parker Instability under a Uniform Gravity Jongsoo Kim 1,2,S.S.Hong 2, Dongsu Ryu 3, and T. W. Jones 4 ABSTRACT Using an isothermal MHD code, we have performed three-dimensional,
More informationSAM GEEN (ITA, HEIDELBERG)
Episode IV: THE RETURN TO SAM GEEN (ITA, HEIDELBERG) WITH PATRICK HENNEBELLE PASCAL TREMBLIN AND JOAKIM ROSDAHL UV FEEDBACK IN CLOUDS HII HII regions, regions, supernovae, supernovae, Molecular Molecular
More informationarxiv: v1 [astro-ph] 12 Jul 2007
SUBMITTED TO THE ASTROPHYSICAL JOURNAL (LETTERS), MAY 6 Preprint typeset using L A TEX style emulateapj v. // THE GENERATION AND DISSIPATION OF INTERSTELLAR TURBULENCE - RESULTS FROM LARGE SCALE HIGH RESOLUTION
More informationMagnetic field structure from Planck polarization observations of the diffuse Galactic ISM
Magnetic field structure from Planck polarization observations of the diffuse Galactic ISM François Boulanger Institut d Astrophysique Spatiale on behalf of the Planck Consortium Outline The Planck data
More informationCreation and destruction of magnetic fields
HAO/NCAR July 30 2007 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More information