Radiation Transfer and Radiation Hydrodynamics
|
|
- Matthew Chase
- 5 years ago
- Views:
Transcription
1 2018/02/02 Winter SHAO Radiation Transfer and Radiation Hydrodynamics Kengo TOMIDA ( 富田賢吾, Osaka University) Some Slides are provided by Yan-Fei Jiang (UCSB)
2 Textbooks Rybicki & Lightman, Radiative Process in Astrophysics Mihalas & Mihalas, Foundations of Radiation Hydrodynamics Shu, The Physics of Astrophysics, Volume I. Radiation Castor, Radiation Hydrodynamics
3 Heat Transport Mechanisms There are three different modes in heat transport Conduction - diffusion of thermal energy not always efficient, sometimes anisotropic relatively easy to implement Convection advection (+mixing) of thermal energy automatically included in the hydrodynamics although high resolution may be required for accuracy Radiation transport by photons (neutrinos) works in many different applications in astrophysics can be diffusion like (optically thick) or streaming (thin) can transport not only energy but also momentum can have much faster time scale than dynamics frequency dependent complicated behaviors
4 Radiation Transfer PROBLEM 1: Photons are collisionless. PROBLEM 2: Photons travel (always) at light speed. Collision less = photons properties are not determined locally. (c.f. hydrodynamics = gas properties are fully set by ρ & T) We have to solve photons properties non-locally, including their distribution not only in space but also in directions. Radiation Transfer equation (Boltzmann equation for photons) Intensity I(r, n, ν, t) [erg s -1 cm -2 Hz -1 sr -1 ] : power carried by photons per unit area, unit solid angle, unit frequency, at position r, time t, frequency ν, and traveling to the direction of n 6 (3 positions, 2 directions, 1 frequency) +1 (time) dimensions!
5 Radiation Transfer Equation Radiation transfer equation: evolution of I along each ray The simplest case free streaming: I ν t + cn I ν = 0 or: I ν + c I ν = 0 (s: distance along the ray) t s Interaction with material source and sink terms: 1 I ν c t + n I ν = j ν k ν I ν j: emissivity, k: absorption (note this term is proportional to I) Note: the source term can depend not only on the gas properties, but also radiation field (I in different directions) e.g. scattering 1 I ν c t + n I ν = j ν k ν + σ ν I ν + σ ν J ν (J ν : mean intensity)
6 Formal Solutions If the radiation time scale is much shorter than the dynamical one, we can omit the time derivative term. Then, along the ray, di ν ds = (k ν+σ ν )I ν + j ν + σ ν J ν Using α ν = k ν + σ ν, dτ ν = α ν ds, S ν = j ν+σ ν J ν α ν di ν dτ ν = I ν + S ν Formal solution: I ν (τ ν ) = I ν (0)e τ ν + 0 τ ν dτ ν S ν (τ ν)e (τ ν τ ν ) Constant source function: I ν τ ν = I ν 0 e τ ν + S ν 1 e τ ν τ ν : optical depth, mean free path: l ν = 1 Τα ν, τ ν ~l ν α ν ~1 Optically thick (τ ν 1): I ν τ ν Optically thin (τ ν 1): I ν τ ν S ν = I ν 0 (1 τ ν ) + S ν τ ν
7 Radiation Moments It is convenient to calculate some quantities integrated (averaged) over directions (angle) J ν = 1 4π න I νdω H ν = 1 4π න I νndω K ν = 1 4π න I νnndω Note: n is the normal unit vector in which the ray is traveling. Note: you can define higher order moments.
8 Radiation Moment Equations It is more convenient to modify them into physically relevant units E ν = 1 c I ν dω Radiation energy density F ν = I ν ndω P ν = 1 c I ν nndω Radiation energy flux Radiation pressure tensor Using these quantities, you can rewrite the transfer equations into moment equations (assuming sources and sinks are isotropic): E ν t + F ν= c(s ν σ ν E ν ) F ν t + c2 P ν = cσ ν F ν Note: the moment equations continue forever
9 Gray Approximation When the radiation field of your interest has a broad spectrum (continuum radiation), it is appropriate to take frequency average. This is often used when we consider thermal radiation (I ν = B ν ). 1 I t + n I = σ a B T I c E t + F = c(σ aat 4 σ E E) F t + c2 P = cσ F F Note: at 4 Stefan-Boltzmann, a= x erg cm 3 K 4 Note: opacities σ must be averaged with an adequate weighting The moment equations are similar to hydrodynamic equations (only) if we can truncate it at a small order of moments
10 Comoving frame Define the radiation quantities on the comoving frame with fluid Relatively simple, but non-conservative terms appear. Laboratory frame Define the radiation quantities on the laboratory (static) frame Transport part remains simple but the RHS becomes a mess Mixed frame Define radiation on the laboratory frame and interaction on the commoving frame, and connect accordingly getting popular Coupling to Moving Fluid When we want to couple radiation with moving fluids, we need to deal with frame transformation, since interaction between fluid and radiation is usually defined on the fluid frame. As radiation is intrinsically relativistic, we need Lorentz transform.
11 Comoving Frame The Transport equation becomes (Lowrie et al. 2000, JQSRT): And the (frequency integrated) moment equations: E t + (ue) + F + P: u = c(σ aat 4 σ E E) F t + uf + c2 P + F u = cσ F F The additional terms appear due to the motion of the frame (conservative) and velocity gradient (non-conservative). Note: The comoving frame equations are not exact, but the approximation is acceptable in non-relativistic regimes.
12 Source Terms in Hydro Equations The gas interacts with the radiation field in two ways: Absorption/emission and radiation force Mass conservation: ρ t + (ρu) = 0 Eq. of motion: ρu t + ρuu + pi = σ F c F Gas energy: e t + e + p u = c σ aat 4 σ E E + σ F c F u Radiation energy: E t + (ue) + F + P: u = c(σ aat 4 σ E E) Radiation flux: F t + uf + c2 P + F u = cσ F F Note: if frequency-dependency is included, we need to incorporate the Doppler shift effect in the frequency space.
13 e.g. Jiang et al. 2012, 2014 Mixed Frame Eq. of motion: ρu + ρuu + pi = S t m e Gas energy: + e + p u = cs t e Radiation energy: E + F = cs t e Radiation flux: 1 c 2 F t + P += S m Mom. coupling: S m = σ a+σ s c Ene. coupling: S e =σ a at 4 E + (σ a σ s ) u c [F (ue + u P)] + u c σ a at 4 E 2 [F (ue + u P)]
14 Different Levels of RHD Level 0: only local source terms of radiation heating/cooling Level 0.5: local source terms with simple optical-depth estimates --- radiation hydro in a broad sense --- Level 1: solve radiation fields but in simplified situations e.g. plane-parallel or direct radiation from a star, ionization only --- radiation hydro with some approximation --- Level 2: solve radiation fields but in an approximated way e.g. flux limited diffusion and other moment methods --- high-level radiation hydro --- Level 3a: solve radiation fields using ray-tracing e.g. variable Eddington factor, direct method (e.g. Jiang+ 2014) Level 3b: solve frequency-dependent radiation fields Level 3c: solve radiation and chemical reactions consistently Level 3d: solve atomic/molecular line transfers
15 Moment Methods The transport equation is a 6 (gray: 5) + 1(time) dimensional problem computationally expensive 1 I c t + n I = σ a B T I The moment equations are 3D (if gray) +1 (time). E t + F = c(σ aat 4 σ E E) F t + c2 P = cσ F F In general the moment equations continue forever. However, IF we can find a physically reasonable closure relation to truncate the moment equations, it can greatly reduce the costs.
16 Simplest Closure: Diffusion Approx. Closure: a relation between higher-moment and lower-moment. The simplest closure: in the optically-thick limit where radiation fields become almost isotropic. The radiation tensor becomes: P = 1 3 EI (I: unit matrix) And you can ignore the time dependency in the flux equation the radiation flux becomes a simple diffusion flux: F = c 3σ E Finally you obtain the radiation diffusion equation: E t c 3σ E = c σ aat 4 σ E E The system is reduced into a single 3D equation!
17 Levermore & Pomraning 1981 Levermore 1984 Flux Limited Diffusion The diffusion approximation is invalid in the optically thin region: it can violate the causality. The flux must be limited by light speed: F = cen, n = E E Flux limited diffusion forces this limit using a flux limiter function F = cλ σ E 2 + R E λ R =, R = 6 + 2R + R2 σe The radiation pressure tensor is modified to connect both limits: P = DE, D = 1 χ 3χ 1 I + nn, χ = λ + λ 2 R Note: although FLD satisfies causality, it is merely a diffusion approximation anyway and not accurate in optically thin regions. We will see this later. Note2: Do not use it with relativity because it still violates causality.
18 Dubroca & Feugeas 1999 Chan et al M1 Closure Scheme Is there any better option? Obviously, it should be better to solve higher order moment, while FLD solves only the zeroth moment. For this, a closure relation for higher-order moment is needed: Eddington Tensor: P=DI D = 1 χ 3χ 1 I + nn, n = F F 2 2 How to calculate χ? Assume either (note: they are equivalent) Radiation field must be isotropic on a certain frame, and is represented by its Lorentz boosting Entropy of radiation should be maximized χ = f 2 F, f = f 2 ce Note: this is Lorentz covariant and is compatible with relativity
19 1D Test FLD M1 A very simple test: radiation flux propagation FLD requires a gradient to produce energy flow. Also FLD requires short timestep because of the non-linearity in the limiter
20 Gonzalez et al D Shadow Test FLD M1 An optically thick obstacle placed in optically thick gas. Radiation comes from left, and should produce a shadow. Obviously FLD fails this test, because it cannot produce a shadow because the scheme is intrinsically isotropic. M1 closure, on the other hand, produces a reasonably good shadow, although there is non-zero diffusion.
21 2D Crossing Beam Test However, M1 closure fails when there are two or more directions in the radiation field, because it assumes the radiation field has only one directionality. We need to calculate the closure directly, or go to the direct method that solves the intensity equations directly.
22 Ray Tracing Long characteristic: cast rays between all the cells Short Characteristic: trace rays cell by cell ART: cast long rays across the grid Monte-Carlo: samples photon tracks randomly Long and Monte-Carlo: O(N 6 ), Short and ART: O(N 5 ) Note: Monte-Carlo require many photons to converge, but it is relatively easy to implement and parallelize, including complicated physics like lines and scattering
23 Variable Eddington Tensor Usually, the time scale of radiation is much faster than dynamics. So we can omit the time-dependency in the transport equations. Then we calculate the Eddington tensor D using ray tracing. Ray tracing itself is relatively cheap as it is just 1D transport equations along rays. (note: scattering makes this complicated.) Once D is given, we update the radiation quantities using the moment equation. Wait we already calculate the radiation intensities. Why do we need to use the moment equations? Can t we directly calculate the radiation moments from the intensities?
24 Advantages of Moment Methods Moment equations are 3D (with a good closure), while the transfer equations are 5D (with the gray approximation). This greatly reduces the computational costs. This is especially beneficial if we use implicit time integrators that requires inversion of a huge matrix to get a consistent solution between all the variables. Also, the moment equations can be treated similarly to the hydrodynamics, and conservation laws can be explicitly applied. On the other hand, if the hydrodynamics is fast enough (e.g. black hole accretion disks), or the reduced speed of light approximation is applicable, the direct method works.
25 (Based on Yan-Fei s slide) Shadow Test with Two Beams Sadowski et al McKinney t al Jiang et al VET can solve difficult situations with more than two directions where M1 closure fails.
26 (Based on Yan-Fei s slide) Disk Test Radiation fields from optically thick disks. Here, M1 closure produces solutions too-strongly directed, while FLD cannot reproduces the shadow. Which scheme do you think is better?
27 (Based on Yan-Fei s slide) Binary Test Radiation fields from two points sources. Again, M1 closure produces structures like a shock in radiation, which is not realistic but an artifact due to treatment like fluid. Interestingly, higher order moment is not necessarily better! Choose your scheme wisely based on their behaviors.
28 Direct Method Instead of solving radiation intensity along rays, it is possible to consider the radiation intensities as volume-averaged quantities within a cell, keeping the time derivative time. This is analogy to the Boltzmann equation. (Jiang et al. 2014) (Ohsuga & Takahashi 2016) Note: this method is recently developed, and is good especially for high-energy system because they both use explicit time integrators. If an implicit integrator is needed, it must be much more expensive.
29 Time Integrators Usually we utilize an explicit time integrator for hydrodynamic simulations and most of additional physics, because it is much simpler to implement, and yet more accurate. U n+1 U n = f(u t n ) However, we have to use small time steps for this method t < x x2 hyperbolic, t < v 2D parabolic or t < τ (where τ is typical time scale) An implicit method is unconditionally stable: U n+1 U n = f(u t n+1 ) But we need to solve the matrix to calculate the unknow U n+1, and often iterations are required if the system is non-linear.
30 Matrix Inversion (FLD case) To solve the non-linear system of radiation transfer (with FLD approx.) we use Newton-Raphson iterations. (I m showing them just to tell you how depressing they are.) The matrix is often ill behaving (not diagonally dominated) A robust sparse matrix solver is needed, but it is difficult (robustness vs scalability).
31 Caution on Implicit Methods Stable does not mean accurate It will approach the steady state in the limit of Δt But the time scale is different (exponential vs polynominal) You often need to use simple, crude, less accurate schemes Large time step does not mean fast or scalable Global matrix inversion is memory and communication intensive Some times the matrix is ill conditioned and difficult to converge It is not trivial to get good performance and scalability Use it only if The accuracy is tolerable for your problem You can get reasonable speed up with it You can develop it within a reasonable time frame
32 Some Useful Techniques Splitting the transport/diffusion term and the source term You can split them and use an implicit integrator for the source. If the stiffness comes from the source term this will work well. Applying the reduced speed of light (signal) approximation When the transport/diffusion speed is very fast, the field should reach a quasi steady state quickly. In this case, it is OK to reduce the transport/diffusion speed as long as the relation between different terms is unchanged. (e.g. Skinner & Ostriker 2013) (Caution: this does not work when there are regions with different time scales) Super Time Stepping (Alexiades et al. 1996) Some schemes special integrators allow larger time steps while keeping the explicit formulation. There are certain conditions (e.g. not work well with hyperbolic systems) but easy to implement.
33 S. JAMSTEC Note on Operator Splitting We cannot split the terms if they are both faster than dynamics!
34 Star Formation Applications Note: these are not exhaustive at all in any sense. 1D radiation hydrodynamics (thermal radiation) Larson 1969, MNRAS, 145, 271 The first simulation of star formation was already radiation hydrodynamics in 1D, implicit, with diffusion approximation Masunaga et al. 1998, Masunaga & Intusuka D RHD (VEF) simulations on low-mass star formation Vaytet et al. 2013, 2017 Systematic survey of 1D star formation models with M1 closure
35 Star Formation Applications Multi-dimensional R(M)HD simulations (thermal, diffuse) Krumholz et al. 2009, Myers et al D R(M)HD (FLD) simulations of massive star formation Whitehouse & Bate 2006, Bate 2010, D SPH RHD (FLD) simulations on low-mass star formation Commercon et al. 2010, Tomida et al. 2010, D RMHD (FLD) simulations of low-mass star formation Multi-dimensional RHD simulations (thermal, direct + diffuse) Yorke & Sonnhalter 2002, Kuiper et al. 2011, D RHD (direct+fld) of massive star formation with feedback. Rosen et al D RHD (direct+fld) of massive star formation with feedback. Harries et al D RHD (Monte-Carlo) of massive star formation with feedback.
36 Radiation Effects in Star Formation Tomida et al The temperature starts to rise when the radiation cooling becomes inefficient The outer region is heated up by radiation from the central hot core.
37 Disk Structure No Radiation With Radiation Radiation heats up the disk and stabilizes against gravity
38 Rosen et al Massive Star Formation Radiation feedback reduces the star formation efficiency 150Ms in the initial cloud ~50Ms after one free-fall time
39 Summary Radiation hydrodynamics is still a developing field. As radiation is the fundamental physics in astrophysics, it will get more and more important (and common). Although the current version only supports source functions, in the near future Athena++ will offer : Direct method for high energy systems Flux Limited Diffusion with an implicit integrator Variable Eddington Factor Method Stay tuned, and have a happy simulation life!
(Numerical) study of the collapse and of the fragmentation of prestellar dense core
(Numerical) study of the collapse and of the fragmentation of prestellar dense core Benoît Commerçon Supervisors: E. Audit, G. Chabrier and P. Hennebelle Collaborator: R. Teyssier (3D - AMR) CEA/DSM/IRFU/SAp
More informationOutline of lectures. Adding more physics. Doing Better: source terms in CTU unsplit integrator
Grid-based methods for hydrodynamics, MHD, and radiation hydrodynamics. Outline of lectures (Four lectures) Lecture 1. Introduction to physics and numerics Lecture 2. Operator split (ZEUS-like) methods
More informationFormation of massive stars: a review
Formation of massive stars: a review Patrick Hennebelle Thanks to: Benoît Commerçon, Marc Joos, Andrea Ciardi Gilles Chabrier, Romain Teyssier How massive stars form? -Can we form massive stars in spite
More informationRadiative transfer equation in spherically symmetric NLTE model stellar atmospheres
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov 20.03.2008 p. Outline 1.
More informationBlack Hole Formation in the Early Universe
1 Black Hole Formation in the Early Universe Luo Yang (Osaka University) XC30 2 1. Introduction Observations of high-redshift quasars reveal that some supermassive black holes (SMBH) with masses exceeding
More informationNew Results from 3-D supernova models with spectral neutrino diffusion
New Results from 3-D supernova models with spectral neutrino diffusion Stuart C. Whitehouse and Matthias Liebendörfer Department of Physics, University of Basel, Switzerland Overview Introduction The Isotropic
More informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More informationFormation of massive stars : a review
Formation of massive stars : a review Patrick Hennebelle (a former «star formation ignorant» Ant s postdoc) Benoit Commerçon, Marc Joos, Andrea Ciardi, Gilles Chabrier One of the Constellation network
More informationAbstract. Department of Astronomy. radiation hydrodynamics (RHD), based on the first two moments of the radiative
Abstract Title of Dissertation: AN EFFICIENT METHOD FOR RADIATION HYDRODYNAMICS IN MODELS OF FEEDBACK-REGULATED STAR FORMATION Michael Aaron Reed Skinner, Doctor of Philosophy, 213 Dissertation directed
More informationRay-Tracing and Flux-Limited-Diffusion for simulating Stellar Radiation Feedback
Ray-Tracing and Flux-Limited-Diffusion for simulating Stellar Radiation Feedback Rolf Kuiper 1,2 H. Klahr 1, H. Beuther 1, Th. Henning 1, C. Dullemond 3, W. Kley 2, R. Klessen 3 1 - Max Planck Institute
More informationRadiation-hydrodynamic Models for ULXs and ULX-pulsars
Radiation-hydrodynamic Models for ULXs and ULX-pulsars Tomohisa KAWASHIMA Division of Theoretical Astrophysics, NAOJ in collaboration with Ken OHSUGA, Hiroyuki TAKAHASHI (NAOJ) Shin MINESHIGE, Takumi OGAWA
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationAstro 305 Lecture Notes Wayne Hu
Astro 305 Lecture Notes Wayne Hu Set 1: Radiative Transfer Radiation Observables From an empiricist s point of view there are 4 observables for radiation Energy Flux Direction Color Polarization Energy
More informationOpacity and Optical Depth
Opacity and Optical Depth Absorption dominated intensity change can be written as di λ = κ λ ρ I λ ds with κ λ the absorption coefficient, or opacity The initial intensity I λ 0 of a light beam will be
More informationModel Atmospheres. Model Atmosphere Assumptions
Model Atmospheres Problem: Construct a numerical model of the atmosphere to estimate (a) Variation of physical variables (T, P) with depth (b) Emergent spectrum in continuum and lines Compare calculated
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationRadiative Transfer Plane-Parallel Frequency-Dependent
4 Radiative Transfer Plane-Parallel Frequency-Dependent variables I ν J ν H ν K ν in astronomy 4. Basic Equations After Fukue, J. 2, PASJ, 63, in press We here assume the followings: i) The disk is steady
More informationIf light travels past a system faster than the time scale for which the system evolves then t I ν = 0 and we have then
6 LECTURE 2 Equation of Radiative Transfer Condition that I ν is constant along rays means that di ν /dt = 0 = t I ν + ck I ν, (29) where ck = di ν /ds is the ray-path derivative. This is equation is the
More informationStars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11)
Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Kenneth Wood, Room 316 kw25@st-andrews.ac.uk http://www-star.st-and.ac.uk/~kw25 What is a Stellar Atmosphere? Transition from dense
More informationElectrodynamics of Radiation Processes
Electrodynamics of Radiation Processes 7. Emission from relativistic particles (contd) & Bremsstrahlung http://www.astro.rug.nl/~etolstoy/radproc/ Chapter 4: Rybicki&Lightman Sections 4.8, 4.9 Chapter
More informationRadiation-hydrodynamics from Mpc to sub-pc scales with RAMSES-RT
Radiation-hydrodynamics from Mpc to sub-pc scales with RAMSES-RT Joki Rosdahl Centre de Recherche Astrophysique de Lyon With Aubert, Blaizot, Bieri, Biernacki, Commercon, Costa, Dubois, Geen, Katz, Kimm,
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 informationOutline. Today we will learn what is thermal radiation
Thermal Radiation & Outline Today we will learn what is thermal radiation Laws Laws of of themodynamics themodynamics Radiative Radiative Diffusion Diffusion Equation Equation Thermal Thermal Equilibrium
More informationNeutrinos, nonzero rest mass particles, and production of high energy photons Particle interactions
Neutrinos, nonzero rest mass particles, and production of high energy photons Particle interactions Previously we considered interactions from the standpoint of photons: a photon travels along, what happens
More information2. NOTES ON RADIATIVE TRANSFER The specific intensity I ν
1 2. NOTES ON RADIATIVE TRANSFER 2.1. The specific intensity I ν Let f(x, p) be the photon distribution function in phase space, summed over the two polarization states. Then fdxdp is the number of photons
More informationde = j ν dvdωdtdν. (1)
Transfer Equation and Blackbodies Initial questions: There are sources in the centers of some galaxies that are extraordinarily bright in microwaves. What s going on? The brightest galaxies in the universe
More informationTheory of optically thin emission line spectroscopy
Theory of optically thin emission line spectroscopy 1 Important definitions In general the spectrum of a source consists of a continuum and several line components. Processes which give raise to the continuous
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationarxiv: v2 [astro-ph.im] 28 Jan 2015
Accepted to ApJS Preprint typeset using L A TEX style emulateapj v. 5/2/11 MONTE CARLO RADIATION-HYDRODYNAMICS WITH IMPLICIT METHODS Nathaniel Roth 1 and Daniel Kasen 1,2,3 Accepted to ApJS arxiv:144.4652v2
More informationCHAPTER 26. Radiative Transfer
CHAPTER 26 Radiative Transfer Consider an incoming signal of specific intensity I ν,0 passing through a cloud (i.e., any gaseous region). As the radiation transits a small path length dr through the cloud,
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More information[Toward] Simulating Cosmological Reionization with Enzo
[Toward] Simulating Cosmological Reionization with Enzo Texas Cosmology Network Meeting 2009 29 October 2009 Daniel R. Reynolds M.L. Norman, P. Paschos, G. So [UCSD] J.C. Hayes [LLNL] Sponsor: NSF Astronomy
More informationAccretion disks. AGN-7:HR-2007 p. 1. AGN-7:HR-2007 p. 2
Accretion disks AGN-7:HR-2007 p. 1 AGN-7:HR-2007 p. 2 1 Quantitative overview Gas orbits in nearly circular fashion Each gas element has a small inward motion due to viscous torques, resulting in an outward
More informationRadiative Processes in Astrophysics
Radiative Processes in Astrophysics 6. Relativistic Covariance & Kinematics Eline Tolstoy http://www.astro.rug.nl/~etolstoy/astroa07/ Practise, practise, practise... mid-term, 31st may, 9.15-11am As we
More informationRADIATIVE TRANSFER IN AXIAL SYMMETRY
Title : will be set by the publisher Editors : will be set by the publisher EAS Publications Series, Vol.?, 26 RADIATIVE TRANSFER IN AXIAL SYMMETRY Daniela Korčáková and Jiří Kubát Abstract. We present
More informationAy Fall 2004 Lecture 6 (given by Tony Travouillon)
Ay 122 - Fall 2004 Lecture 6 (given by Tony Travouillon) Stellar atmospheres, classification of stellar spectra (Many slides c/o Phil Armitage) Formation of spectral lines: 1.excitation Two key questions:
More informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
More informationSIMPLE RADIATIVE TRANSFER
ASTR 511/O Connell Lec 4 1 SIMPLE RADIATIVE TRANSFER The theory of radiative transfer provides the means for determining the emergent EM spectrum of a cosmic source and also for describing the effects
More informationSimulations of massive magnetized dense core collapse
Simulations of massive magnetized dense core collapse Matthias González Benoît Commerçon, Neil Vaytet Raphaël Mignon-Risse Laboratoire AIM, Université Paris Diderot-CEA-CNRS, CEA Saclay Outline 1 Context
More informationAdvanced Cosmological Simulations
Advanced Cosmological Simulations John Wise (Georgia Tech) Enzo Workshop 19 Oct 2013 1 Outline We will consider additional physics in Thursday s AMR (no nested grids) cosmology simulation. Refresher on
More informationKengo TOMIDA Kohji Tomisaka, Masahiro N. Machida, Tomoaki Matsumoto, Satoshi Okuzumi, Yasunori Hori
09/19/2014 Department of Astrophysical Sciences, Princeton University Department of Physics, University of Tokyo JSPS Research Fellow Kengo TOMIDA Kohji Tomisaka, Masahiro N. Machida, Tomoaki Matsumoto,
More informationI ν. di ν. = α ν. = (ndads) σ ν da α ν. = nσ ν = ρκ ν
Absorption Consider a beam passing through an absorbing medium. Define the absorption coefficient, α ν, by ie the fractional loss in intensity in travelling a distance ds is α ν ds (convention: positive
More informationTRANSFER OF RADIATION
TRANSFER OF RADIATION Under LTE Local Thermodynamic Equilibrium) condition radiation has a Planck black body) distribution. Radiation energy density is given as U r,ν = 8πh c 3 ν 3, LTE), tr.1) e hν/kt
More informationEinstein Toolkit Workshop. Joshua Faber Apr
Einstein Toolkit Workshop Joshua Faber Apr 05 2012 Outline Space, time, and special relativity The metric tensor and geometry Curvature Geodesics Einstein s equations The Stress-energy tensor 3+1 formalisms
More informationElectromagnetic Counterparts to Gravitational Wave Detections: Bridging the Gap between Theory and Observation
Electromagnetic Counterparts to Gravitational Wave Detections: Bridging the Gap between Theory and Observation Prof. Zach Etienne, West Virginia University 4 km General Relativity, Briefly Special Relativity:
More informationLecture 3: Emission and absorption
Lecture 3: Emission and absorption Senior Astrophysics 2017-03-10 Senior Astrophysics Lecture 3: Emission and absorption 2017-03-10 1 / 35 Outline 1 Optical depth 2 Sources of radiation 3 Blackbody radiation
More information3D Simulations of Core-collapse Supernovae. Tomoya Takiwaki(NAOJ) Kei Kotake(Fukuoka U) Yudai Suwa(YITP) Tomohide Wada(vis) And many collaborators
2013/12/3 MMCOCOS@Fukuoka University 3D Simulations of Core-collapse Supernovae Tomoya Takiwaki(NAOJ) Kei Kotake(Fukuoka U) Yudai Suwa(YITP) Tomohide Wada(vis) And many collaborators Plan 1. Brief summary
More informationAccretion Disks. 1. Accretion Efficiency. 2. Eddington Luminosity. 3. Bondi-Hoyle Accretion. 4. Temperature profile and spectrum of accretion disk
Accretion Disks Accretion Disks 1. Accretion Efficiency 2. Eddington Luminosity 3. Bondi-Hoyle Accretion 4. Temperature profile and spectrum of accretion disk 5. Spectra of AGN 5.1 Continuum 5.2 Line Emission
More informationˆd = 1 2π. d(t)e iωt dt. (1)
Bremsstrahlung Initial questions: How does the hot gas in galaxy clusters cool? What should we see in their inner portions, where the density is high? As in the last lecture, we re going to take a more
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More informationObservational Appearance of Black Hole Wind Effect of Electron Scattering
Observational Appearance of Black Hole Wind Effect of Electron Scattering Kazuyuki OGURA Astronomical Institute Osaka Kyoiku Univ. 29 Jun 2013 Meeting of BH Horizon Project @Nagoya Univ. Contents Introduction
More informationEquilibrium Structure of Radiation-dominated Disk Segments
Equilibrium Structure of Radiation-dominated Disk Segments Shigenobu Hirose The Earth Simulator Center, JAMSTEC, Japan collaborators Julian Krolik (JHU) Omer Blaes (UCSB) Reconstruction of standard disk
More informationSolution of the radiative transfer equation in NLTE stellar atmospheres
Solution of the radiative transfer equation in NLTE stellar atmospheres Jiří Kubát kubat@sunstel.asu.cas.cz Astronomický ústav AV ČR Ondřejov Non-LTE Line Formation for Trace Elements in Stellar Atmospheres,
More informationIntroduction to the School
Lucio Crivellari Instituto de Astrofísica de Canarias D.pto de Astrofísica, Universidad de La Laguna & INAF Osservatorio Astronomico di Trieste (Italy) Introduction to the School 10/11/17 1 Setting the
More informationRadiative MHD. in Massive Star Formation and Accretion Disks. Rolf Kuiper, Hubert Klahr, Mario Flock, Henrik Beuther, Thomas Henning
Radiative MHD in Massive Star Formation and Accretion Disks, Hubert Klahr, Mario Flock, Henrik Beuther, Thomas Henning, Radiative MHD with Makemake and Pluto : We developed a fast 3D frequency-dependent
More informationPHAS3135 The Physics of Stars
PHAS3135 The Physics of Stars Exam 2013 (Zane/Howarth) Answer ALL SIX questions from Section A, and ANY TWO questions from Section B The numbers in square brackets in the right-hand margin indicate the
More information2. Basic Assumptions for Stellar Atmospheres
2. Basic Assumptions for Stellar Atmospheres 1. geometry, stationarity 2. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres!
More informationElectromagnetic Spectra. AST443, Lecture 13 Stanimir Metchev
Electromagnetic Spectra AST443, Lecture 13 Stanimir Metchev Administrative Homework 2: problem 5.4 extension: until Mon, Nov 2 Reading: Bradt, chapter 11 Howell, chapter 6 Tenagra data: see bottom of Assignments
More informationLet s consider nonrelativistic electrons. A given electron follows Newton s law. m v = ee. (2)
Plasma Processes Initial questions: We see all objects through a medium, which could be interplanetary, interstellar, or intergalactic. How does this medium affect photons? What information can we obtain?
More informationCompschool, Copenhagen Core-Collapse Supernovae. Large cancellation effects in the total energy budget:
Compschool, Copenhagen 2009 Core-Collapse Supernovae M. Liebendörfer University of Basel Collapse phase: Dynamics & ν-interactions Postbounce phase: ν-transport & explosion mechanisms Models: Approximations
More informationDisk modelling by global radiation-mhd simulations
Disk modelling by global radiation-mhd simulations ~Confrontation of inflow & outflow~ Shin Mineshige (Kyoto) & Ken Ohsuga (NAOJ) Magnetic tower jet by RMHD simulation (Takeuchi+11) Outline Introduction
More informationParticle acceleration and generation of high-energy photons
Particle acceleration and generation of high-energy photons For acceleration, see Chapter 21 of Longair Ask class: suppose we observe a photon with an energy of 1 TeV. How could it have been produced?
More informationarxiv:astro-ph/ v1 7 Jul 1999
Gamma-ray Burst Energetics Pawan Kumar Institute for Advanced Study, Princeton, NJ 08540 Abstract arxiv:astro-ph/9907096v1 7 Jul 1999 We estimate the fraction of the total energy in a Gamma-Ray Burst (GRB)
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
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 informationAtomic Physics 3 ASTR 2110 Sarazin
Atomic Physics 3 ASTR 2110 Sarazin Homework #5 Due Wednesday, October 4 due to fall break Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationarxiv: v1 [astro-ph.sr] 23 Nov 2011
Astronomy & Astrophysics manuscript no. 17808-eps c ESO 2011 November 28, 2011 On the stability of radiation-pressure-dominated cavities R. Kuiper 1,2, H. Klahr 2, H. Beuther 2, and Th. Henning 2 1 Argelander
More informationRadiative transfer in planetary atmospheres
Chapter 9 Radiative transfer in planetary atmospheres Warning: This chapter is just a very rough introduction to planetary atmospheres. Since the discovery of planets around other stars than the sun, the
More informationACTIVE GALACTIC NUCLEI: FROM THE CENTRAL BLACK HOLE TO THE GALACTIC ENVIRONMENT
Julian H. Krolik ACTIVE GALACTIC NUCLEI: FROM THE CENTRAL BLACK HOLE TO THE GALACTIC ENVIRONMENT PRINCETON UNIVERSITY PRESS Princeton, New Jersey Preface Guide for Readers xv xix 1. What Are Active Galactic
More information1. Why photons? 2. Photons in a vacuum
Photons and Other Messengers 1. Why photons? Ask class: most of our information about the universe comes from photons. What are the reasons for this? Let s compare them with other possible messengers,
More informationInterstellar Medium Physics
Physics of gas in galaxies. Two main parts: atomic processes & hydrodynamic processes. Atomic processes deal mainly with radiation Hydrodynamics is large scale dynamics of gas. Start small Radiative transfer
More informationHIGH ENERGY ASTROPHYSICS - Lecture 7. PD Frank Rieger ITA & MPIK Heidelberg Wednesday
HIGH ENERGY ASTROPHYSICS - Lecture 7 PD Frank Rieger ITA & MPIK Heidelberg Wednesday 1 (Inverse) Compton Scattering 1 Overview Compton Scattering, polarised and unpolarised light Di erential cross-section
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationZEUS-2D: A Radiation Magnetohydrodynamics Code for. Astrophysical Flows in Two Space Dimensions: III. The. Radiation Hydrodynamic Algorithms and Tests
ZEUS-2D: A Radiation Magnetohydrodynamics Code for Astrophysical Flows in Two Space Dimensions: III. The Radiation Hydrodynamic Algorithms and Tests James M. Stone 1 National Center for Supercomputing
More informationProbing Neutron Star Physics using Thermonuclear X-ray Bursts
Probing Neutron Star Physics using Thermonuclear X-ray Bursts Sudip Bhattacharyya University of Maryland (CRESST) NASA s Goddard Space Flight Center Outline Neutron Stars: why do we care? Thermonuclear
More informationAGN Feedback at the Parsec Scale
AGN Feedback at the Parsec Scale Feng Yuan Shanghai Astronomical Observatory, CAS with: F. G. Xie (SHAO) J. P. Ostriker (Princeton University) M. Li (SHAO) OUTLINE Intermittent activity of compact radio
More informationTHE ROLE OF RADIATION PRESSURE IN HIGH-Z DWARF GALAXIES
THE ROLE OF RADIATION PRESSURE IN HIGH-Z DWARF GALAXIES John Wise (Georgia Tech) Tom Abel (Stanford), Michael Norman (UC San Diego), Britton Smith (Michigan State), Matthew Turk (Columbia) 14 Dec 2012
More informationBremsstrahlung Radiation
Bremsstrahlung Radiation Wise (IR) An Example in Everyday Life X-Rays used in medicine (radiographics) are generated via Bremsstrahlung process. In a nutshell: Bremsstrahlung radiation is emitted when
More informationGravitational Waves from Supernova Core Collapse: Current state and future prospects
Gravitational Waves from Core Collapse Harald Dimmelmeier harrydee@mpa-garching.mpg.de Gravitational Waves from Supernova Core Collapse: Current state and future prospects Work done with E. Müller (MPA)
More informationBremsstrahlung. Rybicki & Lightman Chapter 5. Free-free Emission Braking Radiation
Bremsstrahlung Rybicki & Lightman Chapter 5 Bremsstrahlung Free-free Emission Braking Radiation Radiation due to acceleration of charged particle by the Coulomb field of another charge. Relevant for (i)
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationUNDERSTANDING POWERFUL JETS AT DIFFERENT SCALES IN THE FERMI ERA. B. THE KPC SCALE: X-RAY EMISSION MECHANISM AND JET SPEED
UNDERSTANDING POWERFUL JETS AT DIFFERENT SCALES IN THE FERMI ERA. A. THE PC SCALE: LOCATING THE Γ-RAY EMISSION SITE B. THE KPC SCALE: X-RAY EMISSION MECHANISM AND JET SPEED Markos Georganopoulos 1,2 Eileen
More informationCOX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE
COX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE Extended Second Edition A. Weiss, W. Hillebrandt, H.-C. Thomas and H. Ritter Max-Planck-lnstitut fur Astrophysik, Garching, Germany C S P CONTENTS PREFACE
More informationCore-Collapse Supernovae and Neutrino Transport
Core-Collapse Supernovae and Neutrino Transport SFB TR7 Gravitational Wave Astronomy Video Seminar B.Müller 5.7.2010 Core Collapse Supernovae, General Relativity & Gravitational Waves Core collapse supernova
More information3 Some Radiation Basics
12 Physics 426 Notes Spring 29 3 Some Radiation Basics In this chapter I ll store some basic tools we need for working with radiation astrophysically. This material comes directly from Rybicki & Lightman
More informationSelected Topics in Plasma Astrophysics
Selected Topics in Plasma Astrophysics Range of Astrophysical Plasmas and Relevant Techniques Stellar Winds (Lecture I) Thermal, Radiation, and Magneto-Rotational Driven Winds Connections to Other Areas
More informationRadiative Processes in Astrophysics
Radiative Processes in Astrophysics 11. Synchrotron Radiation & Compton Scattering Eline Tolstoy http://www.astro.rug.nl/~etolstoy/astroa07/ Synchrotron Self-Absorption synchrotron emission is accompanied
More informationarxiv:astro-ph/ v1 1 Mar 1999
A Possible Explanation for the Radio Flare in the Early Afterglow of GRB990123 Xiangdong Shi and Geza Gyuk arxiv:astro-ph/9903023v1 1 Mar 1999 Department of Physics, University of California, San Diego,
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationCMFGEN. Associate every atomic species with a given identification number
CMFGEN Associate every atomic species with a given identification number For each species, set: atomic mass atomic number solar abundance species abbreviations species identification Read in MODEL_SPEC
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationCompton Scattering. hω 1 = hω 0 / [ 1 + ( hω 0 /mc 2 )(1 cos θ) ]. (1) In terms of wavelength it s even easier: λ 1 λ 0 = λ c (1 cos θ) (2)
Compton Scattering Last time we talked about scattering in the limit where the photon energy is much smaller than the mass-energy of an electron. However, when X-rays and gamma-rays are considered, this
More informationStellar Atmospheres: Basic Processes and Equations
Stellar Atmospheres: Basic Processes and Equations Giovanni Catanzaro Abstract The content of this chapter is a very quick summary of key concepts that concern the interaction between photons created in
More informationRadiation Processes. Black Body Radiation. Heino Falcke Radboud Universiteit Nijmegen. Contents:
Radiation Processes Black Body Radiation Heino Falcke Radboud Universiteit Nijmegen Contents: Planck Spectrum Kirchoff & Stefan-Boltzmann Rayleigh-Jeans & Wien Einstein Coefficients Literature: Based heavily
More informationThe First Stars. Simone Ferraro Princeton University. Sept 25, 2012
The First Stars Simone Ferraro Princeton University Sept 25, 2012 Outline Star forming minihalos at high z Cooling physics and chemistry Gravitational Collapse and formation of protostar Magnetic fields
More informationProtostars 1. Early growth and collapse. First core and main accretion phase
Protostars 1. First core and main accretion phase Stahler & Palla: Chapter 11.1 & 8.4.1 & Appendices F & G Early growth and collapse In a magnetized cloud undergoing contraction, the density gradually
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
More informationPhysical Processes in Astrophysics
Physical Processes in Astrophysics Huirong Yan Uni Potsdam & Desy Email: hyan@mail.desy.de 1 Reference Books: Plasma Physics for Astrophysics, Russell M. Kulsrud (2005) The Physics of Astrophysics, Frank
More information