Sink particle accretion test
|
|
- Cory Green
- 5 years ago
- Views:
Transcription
1 Sink particle accretion test David A. Hubber & Stefanie Walch 1 Objectives Simulate spherically-symmetric Bondi accretion onto a sink particle for an isothermal gas. Calculate the accretion rate onto a sink particle for cases where the sink radius is greater than, and less than, the sonic radius. Determine how quickly steady state accretion can be achieved, if at all. For steady-state accretion, compare the accretion rate to the expected Bondi accretion rate for all sink sizes. 2 Introduction Sink particles (Bate et al., 1995) have become a vital algorithmic tool for simulators when performing gravitational collapse problems such as star and planet formation. In the simplest case, sink particles model accretion onto a protostar by defining an accretion radius (often simply called the sink radius) where any matter that falls within the sink radius that is bound, is instantaneously accreted to the sink by adding its mass, momentum and angular momentum to that of the sink. However, there are differences between the various implementations of sink particles in both Smoothed Particle Hydrodynamics (SPH) and mesh codes. SPH codes discretise the fluid into discrete masses. Therefore, the simplest approach to modelling accretion is to simply accrete whole SPH particles when they enter the sink radius. Various accretion criteria can be employed (e.g. Bate et al., 1995) to control how particles accrete. However, this method usually results in a void of particles inside the sink leading to inaccurate hydrodynamical forces near the sink. Hubber et al. (2013) attempted to improve the situation by allowing partial accretion of an SPH particle. Instead of accreting the whole particle, a simple sub-grid accretion model was developed using the two limiting cases of purely spherical accretion, and rotational (i.e. disc) accretion and interpolating between the two cases. Mesh codes use a different approach to sink accretion. In the original mesh sink implementation of Krumholz et al. (2004), they calculate any excess mass inside the sink that would lead to the gravitational instability and then accrete only that excess to the sink. Therefore there is always a non-zero density of gas inside the sink radius to ensure continuity of the density field. Other mesh codes (e.g. Federrath et al., 2010) are based on this approach, although contain some variations. 1
2 2.1 Aim of the test Although simple to implement, sink particles can potentially introduce large errors into simulations if their effects are not treated properly. The biggest problem with sink particles is that removing gas from a simulation also removes its hydrostatic support from the surrounding gas. Instead of modelling a contracting protostar with an inwards pressure gradient (and therefore outwards pressure force), simulators could be modelling a density trough and hence an outwards pressure gradient (and therefore inwards pressure force) which could lead to an artificially enhanced accretion rate. For sink particle prescriptions that do not accrete all gas inside them and use some other presciption for removing gas, then the opposite effect is also possible, i.e. too little gas is removed leading to artificially lower accretion rates. The main aim of this test is to determine (i) if sink particles introduce any unwanted side effects, and (ii) if the effect is to artificially increase or decrease the accretion rate. Since disc accretion is highly dependent on the effective numerical viscosity, we model simple spherically symmetric accretion following Bondi (1952). This will help determine if various sink particle implementations are able to consistently model simplistic steady-state accretion. More sophisticated tests modelling disc accretion may be developed in the future. 3 Test The Bondi accretion problem is described in many articles and text books, so we only summarise the main practical results here which are important in setting up the test. Bondi accretion models the steady-state spherically-symmetric accretion of a non-self gravitating gas onto a central star of mass M from a uniform background gas of density ρ 0. For simplicity, we will model an isothermal gas of pure hydrogen with (dimensionless) temperature T = 1 and dimensionless isothermal sound speed a = 1. The accretion rate onto the star is given by Ṁ = e3/2 πg 2 M 2 ρ 0 a 3, (1) where is the constant e = At large distances from the central star, the gas moves subsonically; near to the star, the gas moves super-sonically such that the gas is effectively in freefall towards the star. The transition point between subsonic and supersonic is called the sonic radius, R SONIC = GM 2a 2. (2) Hubber et al. (2013) found that if the sink particle size, R S, was smaller than the sonic radius (R S R SONIC ), then the lack of hydrodynamical forces was not important and accretion was correctly modelled. If the sink particle was larger than the sonic radius (R S R SONIC ), then any problems with the accretion are apparent and the accretion rate may be incorrectly modelled. Figure 1 shows the behaviour of the original sink particle implementation (black dashed line) and the newer improved sinks (red solid line) for SPH. The computational domain will be of size 8 (in dimensionless units). This can be either a spherical region of radius 8 (if spherical coordinates are used) or a cubic region of range 8 < x < +8 in each dimension. 2
3 Figure 1: Accretion rates for isothermal Bondi accretion onto a sink particle for various sink particle sizes using the Hubber, Walch & Whitworth (2013) sink particle algorithm. The black dashed line shows the case for simple accretion (where particles are accreted immediately as they enter the sink); in this case, accretion is accurate for supersonic infall (i.e. R < R SONIC ) but breaks down for subsonic infall (i.e. R > R SONIC ) where the inaccurate hydrodynamical forces cause large errors. For the new accretion model, the accretion rates are accurate for all sink radii, both for supersonic and subsonic infall. We place a sink particle of mass m SINK at the origin which has a sink radius r SINK. The sink radius will take values either side of the sonic radius. Following Hubber et al. (2013), we take the values R S /R SONIC = 1/8, 1/4, 1/2, 1, 2, 4.8. We will then measure the steady-state accretion rate for each sink particle implementation. We note that the sink radius will be the only variable quantity in this test. The simulation run-time will be measured in units of the quantity t 0 = R SONIC /a which is the time for a sound-wave to cross the distance of a single sonic radius. In order to ensure a steady-state is achieved for each simulation, we will run the simulation for t END = 4t 0. In our dimensionless units, R SONIC = 1/2; therefore t 0 = 1/2 so t = 1/10 and t END = Bondi solution The main initial conditions of the test are a spherically symmetric density and velocity field setup to model the steady-state trans-sonic solution of the Bondi accretion problem (Bondi, 1952). The Bondi equations were solved numerically and outputted to a file BONDI SOLUTION.dat which is available on the Starbench web-page. The file is a simple four-column data file with the columns recording (1) radial position, (2) integrated mass (i.e. all mass between origin and current radius), (3) gas density, and (4) inward velocity (N.B. no negative sign in table). Each solution is in dimensionless form with the radial position being in units of the sonic radius, 3
4 R SONIC, masses in units of 4πρ 0 R 3 SONIC /3, densities in units of the background density, ρ 0, and velocities in units of the isothermal sound speed, a. 3.2 Additional caveats The following caveats must be considered when preparing and running the simulations The gas is not self-gravitating and therefore only gravity due to the sink itself should be considered. The gravitational mass of the sink must not change, even when accreting more mass. This is to ensure the validity of the steady-state solution. This might require the user to modify some lines of code in the respective accretion or gravitational routines. Due to inaccuracies in the hydrodynamics method and the initial accretion of mass at the start of the simulation, it may take a short amount of time for the simulation to settle into a steady-state. This time will be of order the sound-crossing time of the sink, R S /a. Therefore the simulation must be run long enough to allow the numerical steady-state to be achieved. 3.3 SPH setup For the SPH simulations, we have prepared files containing particle positions, masses, velocities and internal energies in order to simplify the process of setting up initial conditions. The files were created by first creating a closed hexagonal-packed array of particles, ignoring all particles greater than unity distance from the origin, hence creating a unit sphere. This lattice is then stretched radially in order to reproduce the density distribution from the tabulated numerical solution. For authors who wish to set-up their own initial conditions, or who are required to set-up them up internally in their own code, then the file BONDI SOLUTION.dat contains the numerical solution. Since the SPH particle distribution is finite, the outer edge of the particle distribution will start to leak away due to a rarefaction wave that propagates inwards. Since the gas is isothermal and moving subsonically at large radii, we can estimate the time it takes for the wave to reach the central sink and therefore ensure the simulation time is short enough to prevent this artifact from corrupting the solution. 3.4 AMR setup For setting-up the mesh simulations, the properties of each cell of the density and velocity should be interpolated from the BONDI SOLUTION.dat file that is provided. Note that the radial positions are logarithmically spaced and therefore must be accounted for when interpolating the table properties to the mesh initial conditions. If the mesh code allows continuous/free-flowing boundary conditions, then these should be employed to minimise any boundary effects in the code. 4
5 For the resolution, it is advised to properly resolve the sonic radius, and/or the sink radius, either with the full grid resolution (if using a uniform grid) or with the use of AMR. To resolve the smallest sink particle (R S = 1 8 R SONIC = 1/16) in a region of size 16 across requires an effective resolution of at least If the mesh code allows mirror symmetry such that only one octant needs to be simulated, then this can reduce the run-time and memory constraints. 4 Output Format Each simulation should generate three separate output files. The first should record the mass of the central sink at regular time intervals t = 1 5 t 0 = 1/10 until the simulation end time (t END = 2). This output filename should have the form SINKACC RSXXX CODENAME.dat. RSXXX is the sink radius (e.g. RS0.25) For mesh simulations, the resolution should also be appended to the filename description (e.g. SINKACC RS0.125 RES256 CODENAME.dat). The format of these files should be a simple 2-column format, the first column with the time and the second column with the mass of the sink. The second and third files should contain the radial profiles of the gas properties at the beginning and end of the simulation, i.e. either grid cells or SPH particles. These files should contain the radial positions, densities and radial velocity component (i.e. v ˆr) for each gas element in a simple 3-column format. The filenames should have the form BONDI TFIRST RSXXX CODENAME.dat and BONDI TLAST RSXXX CODENAME.dat for the first and last snapshots respectively. For mesh codes, this should also include the resolution as described above. References Bate, M. R., Bonnell, I. A., and Price, N. M. (1995). Modelling accretion in protobinary systems. MNRAS, 277: Bondi, H. (1952). On spherically symmetrical accretion. MNRAS, 112: Federrath, C., Banerjee, R., Clark, P. C., and Klessen, R. S. (2010). Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH. ApJ, 713: Hubber, D. A., Walch, S., and Whitworth, A. P. (2013). An improved sink particle algorithm for SPH simulations. MNRAS, 430: Krumholz, M. R., McKee, C. F., and Klein, R. I. (2004). Embedding Lagrangian Sink Particles in Eulerian Grids. ApJ, 611:
Collapse 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 informationThe Initial Mass Function Elisa Chisari
The Initial Mass Function AST 541 Dec 4 2012 Outline The form of the IMF Summary of observations Ingredients of a complete model A Press Schechter model Numerical simulations Conclusions The form of the
More informationmc 2, (8.1) = R Sch 2R
Chapter 8 Spherical Accretion Accretion may be defined as the gravitational attraction of material onto a compact object. The compact object may be a black hole with a Schwarzschild radius R = 2GM /c 2
More informationAccretion phase of star formation in clouds with different metallicities
doi:10.1093/mnras/stu2633 Accretion phase of star formation in clouds with different metallicities Masahiro N. Machida and Teppei Nakamura Department of Earth and Planetary Sciences, Faculty of Sciences,
More informationThe Effects of Radiative Transfer on Low-Mass Star Formation
The Effects of Radiative Transfer on Low-Mass Star Formation Stella Offner NSF Fellow, ITC Dense Cores in Dark Clouds Oct 23 2009 Collaborators: Chris McKee (UC Berkeley), Richard Klein (UC Berkeley; LLNL),
More informationastro-ph/ Oct 1995
Mon. Not. R. Astron. Soc. 000, 000{000 (1995) Modelling accretion in protobinary systems Matthew R. Bate, Ian A. Bonnell and Nigel M. Price Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA Accepted
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 informationFLASH Code Tutorial. part III sink particles & feedback. Robi Banerjee Hamburger Sternwarte
FLASH Code Tutorial part III sink particles & feedback Robi Banerjee Hamburger Sternwarte banerjee@hs.uni-hamburg.de Motivation: modelling of dense regions in collapse simulations, e.g. star formation
More informationStellar structure and evolution. Pierre Hily-Blant April 25, IPAG
Stellar structure and evolution Pierre Hily-Blant 2017-18 April 25, 2018 IPAG pierre.hily-blant@univ-grenoble-alpes.fr, OSUG-D/306 10 Protostars and Pre-Main-Sequence Stars 10.1. Introduction 10 Protostars
More informationLecture 22 Stability of Molecular Clouds
Lecture 22 Stability of Molecular Clouds 1. Stability of Cloud Cores 2. Collapse and Fragmentation of Clouds 3. Applying the Virial Theorem References Myers, Physical Conditions in Molecular Clouds in
More informationAccretion Mechanisms
Massive Protostars Accretion Mechanism Debate Protostellar Evolution: - Radiative stability - Deuterium shell burning - Contraction and Hydrogen Ignition Stahler & Palla (2004): Section 11.4 Accretion
More informationThe Fragmentation of expanding shells. Kazunari Iwasaki (Nagoya Univ.) Collaborate with S. Inutsuka (Nagoya Univ.) T. Tsuribe (Osaka Univ.
.... The Fragmentation of expanding shells Kazunari Iwasaki (Nagoya Univ.) Collaborate with S. Inutsuka (Nagoya Univ.) T. Tsuribe (Osaka Univ.) 1 / 23 Outline Introduction Three Dimensional SPH Simulations
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 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 informationVariability of accreting black holes induced by shocks in low angular momentum flows
induced by shocks in low angular momentum flows Astronomical Institute of the CAS Prague, Czech Republic Cooperation with Agnieszka Janiuk, CFT PAN, Vladimír Karas ASU CAS 23.10.2017 Low angular momentum
More informationBinary star formation
Binary star formation So far we have ignored binary stars. But, most stars are part of binary systems: Solar mass stars: about 2 / 3 are part of binaries Separations from: < 0.1 au > 10 3 au Wide range
More informationThe effect of magnetic fields on the formation of circumstellar discs around young stars
Astrophysics and Space Science DOI 10.1007/sXXXXX-XXX-XXXX-X The effect of magnetic fields on the formation of circumstellar discs around young stars Daniel J. Price and Matthew R. Bate c Springer-Verlag
More informationASP2062 Introduction to Astrophysics
School of Physics and Astronomy ASP2062 2015 ASP2062 Introduction to Astrophysics Star formation II Daniel Price Key revision points 1. Star formation is a competition between gravity and pressure 2. Timescale
More informationGravitational fragmentation of discs can form stars with masses
Gravitational fragmentation of discs can form stars with masses from ~3 M J to ~200 M J (0.2M ) Defining stars, brown dwarfs Stars and planets Objects formed by gravitational instability on a dynamical
More informationStar formation Part III
Lecture 4 Star formation Part III Lecture Universität Heidelberg WS 11/12 Dr. C. Mordasini Based partially on script of Prof. W. Benz Matthew R. Bate Mentor Prof. T. Henning Lecture 4 overview 1. Heating
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 informationTowards Understanding Simulations of Galaxy Formation. Nigel Mitchell. On the Origin of Cores in Simulated Galaxy Clusters
Towards Understanding Simulations of Galaxy Formation Nigel Mitchell On the Origin of Cores in Simulated Galaxy Clusters Work published in the Monthly Notices of the Royal Astronomy Society Journal, 2009,
More informationThe Migration of Giant Planets in Massive Protoplanetary Discs
The Migration of Giant Planets in Massive Protoplanetary Discs Kate E. M. Robinson A thesis submitted in partial fulfilment for the requirements for the degree of Master of Science at the University of
More informationThe dependence of star cluster formation on initial conditions. Matthew Bate University of Exeter
The dependence of star cluster formation on initial conditions Matthew Bate University of Exeter Stellar properties do not greatly depend on initial conditions constellation Little evidence for variation
More informationAn overview of star formation
An overview of star formation Paul Clark ITA: Ralf Klessen Robi Banerjee Simon Glover Ian Bonnell Clare Dobbs Jim Dale Why study star formation? Stars chemically the enrich the Universe, so star formation
More informationIV From cores to stars
IV From cores to stars 4.0 The Jeans condition When the supporting pressure in a region is not able to hold that region up against gravitational collapse it is said to be Jeans unstable. The Jeans length
More informationThe accuracy of the gravitational potential was determined by comparison of our numerical solution
the results of three key accuracy tests. 48 A6.1. Gravitational Potential The accuracy of the gravitational potential was determined by comparison of our numerical solution with the analytical result for
More informationON THE EVOLUTION OF ULTRACOMPACT H ii REGIONS Eric Keto
Astrophysical Journal, 580:980 986, 2002 December # 2002. The American Astronomical Society. All rights reserved. Printed in U.S.A. ON THE EVOLUTION OF ULTRACOMPACT H ii REGIONS Eric Keto Harvard-Smithsonian
More informationN-body Dynamics in Stellar Clusters Embedded in Gas
N-body Dynamics in Stellar Clusters Embedded in Gas Matthew Bate, University of Exeter Kurosawa, Harries, Bate & Symington (2004) Collaborators: Ian Bonnell, St Andrews Volker Bromm, Texas Star Formation
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 27 May, 2004 1.30 to 3.30 PAPER 64 ASTROPHYSICAL FLUID DYNAMICS Attempt THREE questions. There are four questions in total. The questions carry equal weight. Candidates
More informationThe Competitive Accretion Debate
The Competitive Accretion Debate 1,2 Paul C. Clark 2 Ralf S. Klessen 3 Ian A. Bonnell 3 Rowan J. Smith 1 KITP 2 University of Heidelberg 3 University of St Andrews What is CA and how does it work? Theory
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 informationarxiv: v1 [astro-ph.sr] 29 Jul 2011
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 2 August 2011 (MN LATEX style file v2.2) arxiv:1108.0009v1 [astro-ph.sr] 29 Jul 2011 Collapse of a molecular cloud core to stellar densities: the formation
More informationFormation Mechanisms of Brown Dwarfs: Observations & Theories. Dan Li April 2009
Formation Mechanisms of Brown Dwarfs: Observations & Theories Dan Li April 2009 What is brown dwarf (BD)? BD Mass : upper-limit ~ 0.075 M lower-limit ~ 0.013 M (?) Differences between BD and giant planet:
More informationEnrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México
Enrique Vázquez-Semadeni Centro de Radioastronomía y Astrofísica, UNAM, México 1 Javier Ballesteros-Paredes Centro de Radioastronomía y Astrofísica, UNAM, México 2 Collaborators: Javier Ballesteros-Paredes
More informationNumerical Simulations of Star Cluster Formation. Matthew Bate
Numerical Simulations of Star Cluster Formation Matthew Bate http://www.astro.ex.ac.uk/people/mbate Typical molecular cloud Jeans mass 1 M, Opacity limit 3 MJ, P(k) k -4 Denser cloud Jeans mass 1/3 M Lower
More informationCosmic Structure Formation on Supercomputers (and laptops)
Cosmic Structure Formation on Supercomputers (and laptops) Lecture 4: Smoothed particle hydrodynamics and baryonic sub-grid models Benjamin Moster! Ewald Puchwein 1 Outline of the lecture course Lecture
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 informationStellar Winds. Star. v w
Stellar Winds Star v w Stellar Winds Geoffrey V. Bicknell 1 Characteristics of stellar winds Solar wind Velocity at earth s orbit: Density: Temperature: Speed of sound: v 400 km/s n 10 7 m 3 c s T 10 5
More informationCollapse of magnetized dense cores. Is there a fragmentation crisis?
Collapse of magnetized dense cores Is there a fragmentation crisis? Patrick Hennebelle (ENS-Observatoire de Paris) Collaborators: Benoît Commerçon, Andréa Ciardi, Sébastien Fromang, Romain Teyssier, Philippe
More informationThe Formation of Population III stars in a ΛCDM universe
Chapter 4 The Formation of Population III stars in a ΛCDM universe 4.1 Summary In this chapter I discuss aspects of primordial star formation in a ΛCDM universe. The collapse of gas in a representative
More informationarxiv: v1 [astro-ph.sr] 28 Aug 2009
Outflow Feedback Regulated Massive Star Formation in Parsec-Scale Cluster Forming Clumps Peng Wang 1, Zhi-Yun Li 2, Tom Abel 1 and Fumitaka Nakamura 3 arxiv:0908.4129v1 [astro-ph.sr] 28 Aug 2009 ABSTRACT
More informationDisk Formation and Jet Driving in Collapsing Cloud Cores
Disk Formation and Jet Driving in Collapsing Cloud Cores Masahiro Machida (Kyushu University) Star Formation Process Observations have shown that Low-velocity outflows and high-velocity jets are ubiquitous
More informationarxiv: v1 [astro-ph.sr] 23 Sep 2014
Mon. Not. R. Astron. Soc. 000,?? (04) Printed 4 September 04 (MN LATEX style file v.) arxiv:409.658v [astro-ph.sr] 3 Sep 04 Towards a more realistic sink particle algorithm for the RAMSES code Andreas
More informationStellar, brown dwarf and multiple star properties from a radiation hydrodynamical simulation of star cluster formation
Mon. Not. R. Astron. Soc. 419, 3115 3146 (2012) doi:10.1111/j.1365-2966.2011.19955.x Stellar, brown dwarf and multiple star properties from a radiation hydrodynamical simulation of star cluster formation
More informationEnrique Vázquez-Semadeni. Centro de Radioastronomía y Astrofísica, UNAM, México
Enrique Vázquez-Semadeni Centro de Radioastronomía y Astrofísica, UNAM, México 1 Collaborators: CRyA UNAM: Abroad: Javier Ballesteros-Paredes Pedro Colín Gilberto Gómez Recent PhDs: Alejandro González
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 informationBROWN DWARF FORMATION BY DISC FRAGMENTATION
BROWN DWARF FORMATION BY DISC FRAGMENTATION Ant Whitworth, Cardiff Dimitri Stamatellos, Cardiff plus Steffi Walch, Cardiff Murat Kaplan, Cardiff Simon Goodwin, Sheffield David Hubber, Sheffield Richard
More informationSAM GEEN (ITA/ZAH HEIDELBERG)
SAM GEEN (ITA/ZAH HEIDELBERG) WITH PATRICK HENNEBELLE JUAN SOLER AND RALF KLESSEN Credit: Lost Valley Observatory Star formation is self regulating HII HII regions, regions, supernovae supernovae Molecular
More informationImpact of Protostellar Outflow on Star Formation: Effects of Initial Cloud Mass
Impact of Protostellar Outflow on Star Formation: Effects of Initial Cloud Mass Masahiro N. Machida 1 and Tomoaki Matsumoto 2 ABSTRACT Star formation efficiency controlled by the protostellar outflow in
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 informationDoes a prestellar core always become protostellar? Tracing the evolution of cores from the prestellar to protostellar phase
MNRAS 430, 1854 1866 (2013) doi:10.1093/mnras/stt006 Does a prestellar core always become protostellar? Tracing the evolution of cores from the prestellar to protostellar phase S. Anathpindika 1 and James
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 informationTheory of star formation
Theory of star formation Monday 8th 17.15 18.00 Molecular clouds and star formation: Introduction Tuesday 9th 13.15 14.00 Molecular clouds: structure, physics, and chemistry 16.00 16.45 Cloud cores: statistics
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 informationSPH simulations of star/planet formation triggered by cloud-cloud collisions
Exoplanets: Detection, Formation and Dynamics Proceedings IAU Symposium No. 249, 2008 Y.-S. Sun, S. Ferraz-Mello & J.-L. Zhou, eds. c 2008 International Astronomical Union DOI: 00.0000/X000000000000000X
More informationPAPER 68 ACCRETION DISCS
MATHEMATICAL TRIPOS Part III Tuesday, 2 June, 2009 9:00 am to 11:00 am PAPER 68 ACCRETION DISCS There are THREE questions in total Full marks can be obtained by completing TWO questions The questions carry
More informationStar formation. Protostellar accretion disks
Star formation Protostellar accretion disks Summary of previous lectures and goal for today Collapse Protostars - main accretion phase - not visible in optical (dust envelope) Pre-main-sequence phase -
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 informationEquations of Stellar Structure
Equations of Stellar Structure Stellar structure and evolution can be calculated via a series of differential equations involving mass, pressure, temperature, and density. For simplicity, we will assume
More informationHigh Energy Astrophysics
High Energy Astrophysics Accretion Giampaolo Pisano Jodrell Bank Centre for Astrophysics - University of Manchester giampaolo.pisano@manchester.ac.uk April 01 Accretion - Accretion efficiency - Eddington
More informationThe Origin of the Initial Mass Function
The Origin of the Initial Mass Function Ian A. Bonnell University of St Andrews Richard B. Larson Yale University Hans Zinnecker Astrophysikalisches Institut Potsdam Bonnell et al.: Origin of the Initial
More informationarxiv:astro-ph/ v1 17 Mar 2006
The Origin of the Initial Mass Function Ian A. Bonnell University of St Andrews Richard B. Larson Yale University Hans Zinnecker Astrophysikalisches Institut Potsdam arxiv:astro-ph/0603447v1 17 Mar 2006
More informationA hybrid SPH/N-body method for star cluster simulations
MNRAS 430, 1599 1616 (2013) doi:10.1093/mnras/sts694 A hybrid SPH/N-body method for star cluster simulations D. A. Hubber, 1,2 R. J. Allison, 1,3 R. Smith 4 and S. P. Goodwin 1 1 Department of Physics
More information3D SPH Simulations of Shocks in Accretion Flows around black holes
Mon. Not. R. Astron. Soc., () Printed 8 December 218 (MN LATEX style file v2.2) 3D SPH Simulations of Shocks in Accretion Flows around black holes G. Gerardi 1, D. Molteni 1 and V. Teresi 1 1 Dipartimento
More informationAdvection Dominated Accretion Flows. A Toy Disk Model. Bohdan P a c z y ń s k i
ACTA ASTRONOMICA Vol. 48 (1998) pp. 667 676 Advection Dominated Accretion Flows. A Toy Disk Model by Bohdan P a c z y ń s k i Princeton University Observatory, Princeton, NJ 8544-11, USA e-mail: bp@astro.princeton.edu
More informationCollapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores
Zurich Open Repository and Archive University of Zurich Main Library Strickhofstrasse 39 CH-8057 Zurich www.zora.uzh.ch Year: 2011 Collapse, outflows and fragmentation of massive, turbulent and magnetized
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 informationFrom Massive Cores to Massive Stars
From Massive Cores to Massive Stars Mark Krumholz Princeton University / UC Santa Cruz Collaborators: Richard Klein, Christopher McKee (UC Berkeley) Kaitlin Kratter, Christopher Matzner (U. Toronto) Jonathan
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationPoS(ICRC2017)1086. Accretion onto Black Holes. Myeong-Gu Park. Kyungpook National University, Daegu 41566, KOREA
Accretion onto Black Holes Kyungpook National University, Daegu 41566, KOREA E-mail: mgp@knu.ac.kr Most energetic events or objects in the universe are explained by accretion onto black holes. However,
More information[2 marks] Show that derivative of the angular velocity. What is the specific angular momentum j as a function of M and R in this Keplerian case?
k!! Queen Mary University of London M. Sc. EXAM I N AT1 0 N ASTMOOS Angular Momentum and Accretion in Astrophysics Fkiday, 26th May, 2006 18:15-19:45 Time Allowed: lh 30m This paper has two Sections and
More informationarxiv:astro-ph/ v1 11 Aug 2004
Astronomy & Astrophysics manuscript no. February 2, 2008 (DOI: will be inserted by hand later) Planets opening dust gaps in gas disks Sijme-Jan Paardekooper 1 and Garrelt Mellema 2,1 arxiv:astro-ph/0408202v1
More informationHigh-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion. Overview
High-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion Robert Laing Overview Evidence for black holes in galaxies and techniques for estimating their mass Simple
More informationLecture 6: Protostellar Accretion
Lecture 6: Protostellar Accretion Richard Alexander, 23rd March 2009 The angular momentum problem When we try to understand how young stars form, the obvious question to ask is how do young stars accrete
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 informationDr G. I. Ogilvie Lent Term 2005 INTRODUCTION
Accretion Discs Mathematical Tripos, Part III Dr G. I. Ogilvie Lent Term 2005 INTRODUCTION 0.1. Accretion If a particle of mass m falls from infinity and comes to rest on the surface of a star of mass
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 informationEvolution of protoplanetary discs
Evolution of protoplanetary discs and why it is important for planet formation Bertram Bitsch Lund Observatory April 2015 Bertram Bitsch (Lund) Evolution of protoplanetary discs April 2015 1 / 41 Observations
More informationWhere do Stars Form?
Where do Stars Form? Coldest spots in the galaxy: T ~ 10 K Composition: Mainly molecular hydrogen 1% dust EGGs = Evaporating Gaseous Globules ftp://ftp.hq.nasa.gov/pub/pao/pressrel/1995/95-190.txt Slide
More informationThe Jeans mass and the origin of the knee in the IMF
Mon. Not. R. Astron. Soc. 368, 1296 1300 (2006) doi:10.1111/j.1365-2966.2006.10214.x The Jeans mass and the origin of the knee in the IMF I. A. Bonnell, 1 C. J. Clarke 2 and M. R. Bate 3 1 School of Physics
More informationastro-ph/ Nov 94
Mon. Not. R. Astron. Soc. 000, 000{000 (1994) The formation of close binary systems Ian A. Bonnell and Matthew R. Bate Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA Accepted 1994 August 31.
More informationStellar Atmospheres. University of Denver, Department of Physics and Astronomy. Physics 2052 Stellar Physics, Winter 2008.
Stellar Atmospheres University of Denver, Department of Physics and Astronomy Physics 2052 Stellar Physics, Winter 2008 By Kathy Geise Introduction A star does not have a solid surface, so the definition
More informationThe Stellar Black Hole
The Stellar Black Hole Kenneth Dalton e-mail: kxdalton@yahoo.com Abstract A black hole model is proposed in which a neutron star is surrounded by a neutral gas of electrons and positrons. The gas is in
More informationBONDI-HOYLE ACCRETION IN A TURBULENT MEDIUM
The Astrophysical Journal, 638:369 381, 2006 February 10 # 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. A BONDI-HOYLE ACCRETION IN A TURBULENT MEDIUM Mark R. Krumholz
More informationarxiv:astro-ph/ v1 17 Mar 2006
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 21 January 2014 (MN LaT E X style file v2.2) The Jeans mass and the origin of the knee in the IMF arxiv:astro-ph/0603444v1 17 Mar 2006 I. A. Bonnell
More informationwhile the Planck mean opacity is defined by
PtII Astrophysics Lent, 2016 Physics of Astrophysics Example sheet 4 Radiation physics and feedback 1. Show that the recombination timescale for an ionised plasma of number density n is t rec 1/αn where
More informationM J v2 s. Gm (r) P dm = 3. For an isothermal gas, using the fact that m r, this is GM 2 R = 3N ot M 4πR 3 P o, P o = 3N ot M
1 Stellar Birth Criterion for stability neglect magnetic fields. Critical mass is called the Jean s mass M J, where a cloud s potential energy equals its kinetic energy. If spherical with uniform temperature,
More informationStar Cluster Formation and the Origin of Stellar Properties. Matthew Bate University of Exeter
Star Cluster Formation and the Origin of Stellar Properties Matthew Bate University of Exeter Typical molecular cloud (Bate et al. 2003) Denser cloud (Bate & Bonnell 2005) Jeans mass 1 M, Opacity limit
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 informationThree Dimensional Models of RR Lyrae Pulsation
Regional Variable Star Conference: Physics & Astronomy Department, Michigan State University: 40 Years of Variable Stars: A Celebration of Contributions by Horace A. Smith ed. K. Kinemuchi (Sunspot, NM:
More informationThe Formation of Star Clusters
The Formation of Star Clusters Orion Nebula Cluster (JHK) - McCaughrean Jonathan Tan University of Florida & KITP In collaboration with: Brent Buckalew (ERAU), Michael Butler (UF u-grad), Jayce Dowell
More informationThe effect of magnetic fields on star cluster formation
Mon. Not. R. Astron. Soc. 000, 1 16 (2007) Printed 15 January 2008 (MN LATEX style file v2.2) The effect of magnetic fields on star cluster formation Daniel J. Price and Matthew R. Bate School of Physics,
More informationPOPULATION III STAR FORMATION IN A CDM UNIVERSE. I. THE EFFECT OF FORMATION REDSHIFT AND ENVIRONMENT ON PROTOSTELLAR ACCRETION RATE
The Astrophysical Journal, 654:66Y92, 2007 January 1 # 2007. The American Astronomical Society. All rights reserved. Printed in U.S.A. POPULATION III STAR FORMATION IN A CDM UNIVERSE. I. THE EFFECT OF
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 informationAstro Instructors: Jim Cordes & Shami Chatterjee.
Astro 2299 The Search for Life in the Universe Lecture 8 Last time: Formation and function of stars This time (and probably next): The Sun, hydrogen fusion Virial theorem and internal temperatures of stars
More informationMulti-physics CFD simulation of three-phase flow with MPS method
APCOM & ISCM 11-14 th December, 2013, Singapore Abstract Multi-physics CFD simulation of three-phase flow with MPS method *Ryouhei Takahashi¹, Makoto Yamamoto 2 and Hiroshi Kitada 1 1 CMS Corporation,
More informationObservational signatures of proto brown dwarf formation in protostellar disks
Mem. S.A.It. Vol. 84, 88 c SAIt 213 Memorie della Observational signatures of proto brown dwarf formation in protostellar disks O. V. Zakhozhay 1, E. I. Vorobyov 2,3, and M. M. Dunham 4 1 Main Astronomical
More informationOutflows from young stellar objects and their impact on star formation
Outflows from young stellar objects and their impact on star formation Robi Banerjee 1, @ 1 : Hamburg Sternewarte, University of Hamburg Gojenbergsweg 112 21209 Hamburg - Germany Jets and outflows are
More information