1 / 27

GRMHD Astrophysics Simulations using Cosmos++

GRMHD Astrophysics Simulations using Cosmos++. Joseph Niehaus , Chris Lindner, Chris Fragile. Why do Computational Astrophysics?. Tests the extremes of space that cannot be simulated by conventional means Many vital parameters cannot be observed Many problems have no exploitable symmetry.

javen
Download Presentation

GRMHD Astrophysics Simulations using Cosmos++

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. GRMHD Astrophysics Simulations using Cosmos++ Joseph Niehaus, Chris Lindner, Chris Fragile

  2. Why do Computational Astrophysics? • Tests the extremes of space that cannot be simulated by conventional means • Many vital parameters cannot be observed • Many problems have no exploitable symmetry

  3. Finite Volume Simulations • Divide the computational area into zones • Each zone contains essential data about the material contained inside • The simulation is evolved in time through a series of time steps • As the simulation progresses, cells communicate with each other

  4. Highlights of Cosmos++ • Developers: P. Anninos, P. C. Fragile, J. Salmonson, & S. Murray • Anninos & Fragile (2003) ApJS, 144, 243 • Anninos, Fragile, & Murray (2003) ApJS, 147, 177 • Anninos, Fragile & Salmonson (2005) ApJ, 635, 723 • Multi-dimensional Arbitrary-Lagrange-Eulerian (ALE) fluid dynamics code • 1, 2, or 3D unstructured mesh • Local Adaptive Mesh Refinement (Khokhlov 1998)

  5. Highlights of Cosmos++ • Multi-physics code for Astrophysics/Cosmology • Newtonian & GR MHD • Arbitrary spacetime curvature (K. Camarda -> Evolving GRMHD) • Relativistic scalar fields • Radiation transport (Flux-limited diffusion -> Monte Carlo) • Equilibrium & Non-Equilibrium Chemistry (30+ reactions) • Radiative Cooling • Newtonian external & Self-gravity • Developed for large parallel computation • LLNL Thunder, NCSA Teragrid, NASA Columbia, JPL Cosmos, BSC MareNostrum

  6. Local Adaptive Mesh Refinement

  7. GRMHD Equations in Cosmos++Extended Artificial Viscosity (eAV) mass conservation momentum conservation induction “divergence cleanser”

  8. Active Galaxy Centaurus A

  9. Describing a Black Hole • Three possible intrinsic properties: • Mass • Angular momentum (spin) • Electric charge • Nothing else can be known about a black hole • “No hair” theorem Astrophysically unlikely

  10. Black Hole Accretion Disks • Often formed from binary star systems • Black hole accretes matter from donor star • Disk of plasma forms around black hole • Angular momentum is exchanged through • Magnetic fields • Magnetically dominated flux points away from black hole’s poles, forming jets

  11. Accretion Disks: What we don’t know • Jets • What powers jets? • What sets their orientation? • How is the black hole oriented? • Cooling and Heating • What type of radiative transport occurs in the disk? • How does this effect disk structure? • How does this effect what we observe? • QPOs • What is the source of these phenomena? Total intensity image at 4.85 GHz of SS433 Blundell, K. M. & Bowler, M. G., 2004, ApJ, 616, L159

  12. What determines jet orientation in accretion disk systems? We can answer this question by simulating systems where the angular momentum of the disk is not aligned with the angular momentum of The black hole • “Tilted accretion disks” • (Fragile, Mathews, & Wilson, 2001, Astrophys. J., 553, 955) • Can arise from asymmetric binary systems • Breaks the main degeneracy in the problem

  13. Most commonly used type of grid for accretion disk simulations good angular momentum conservation easy to accommodate event horizon Not very good for simulating jets in 3D zones get very small along pole forcing a very small integration timestep pole is a coordinate singularity creates problems, particularly for transport of fluid across the pole Spherical-Polar Grid

  14. Common in atmospheric codes Not seen as often in astrophysics Adequate for simulating disks good angular momentum conservation easily accommodates event horizon Advantages for simulating jets nearly uniform zone sizing over entire grid no coordinate singularities (except origin) Cubed-Sphere Grid

  15. The Cubed Sphere Six cubes are projected into segments of a sphere Each block has its own coordinate system

  16. Jet Orientation

  17. Energy Equations in Cosmos++Extended Artificial Viscosity (eAV) internal energy total energy conservation

  18. Why Two Energy Equations? • Tracked Simultaneously through code • Attempt to recapture as much heat as possible • Attempting to counteract numerical diffusion • Used when total energy below error • Both energies compared if both below error • Higher energy chosen

  19. Heating Processes • Magnetic • Magnetic Reconnection • Recaptured through total energy equation • No explicit term • Hydrodynamic • Shockwaves & Gas Compression • Handled directly by both energy equations • Viscous • Internal heating due to fluid dynamics • Recaptured through total energy

  20. Radiative Cooling Processes • Bremsstrahlung • “Braking” cooling, emits radiation when decelerating • Synchrotron • Relativistic electrons & positrons • Inverse Compton • Electrons colliding with photons • Becomes prevalent as optical depth increases

  21. Radiative Cooling Processes

  22. 2.5D Simulations • Initial stable solution for rotating torus • Set up for MRI growth • Poloidal fields • No mass or energy transported azimuthally • Vectors tracked numerically

  23. 2.5D Simulations • 3 Scenarios for Comparison • M • Similar to past runs • No heating or cooling • Physical assumption • TM • Heating included • Total energy & Internal energy equations • TMC • Heating and Cooling Processes • Total energy & Internal energy

  24. 2D Simulations - Results torus2d.m.h torus2d.tm.h torus2d.tmc.h

  25. Conclusions • Cosmos++ • GR MHD • AMR • Radiative cooling • Accretion Disks • Cooling/Heating • Jets/Tilted Disks • QPO’s

  26. Untilted Disk Jets Unbound Material Magnetic Field Lines

  27. 15 Degree Tilt Jets Unbound Material Magnetic Field Lines

More Related