test Browse by Author Names Browse by Titles of Works Browse by Subjects of Works Browse by Issue Dates of Works
       

Advanced Search
Home   
 
Browse   
Communities
& Collections
  
Issue Date   
Author   
Title   
Subject   
 
Sign on to:   
Receive email
updates
  
My Account
authorized users
  
Edit Profile   
 
Help   
About T-Space   

T-Space at The University of Toronto Libraries >
School of Graduate Studies - Theses >
Doctoral >

Please use this identifier to cite or link to this item: http://hdl.handle.net/1807/29832

Title: Magneto-hydrodynamics Simulation in Astrophysics
Authors: Pang, Bijia
Advisor: Pen, Ue-Li
Department: Physics
Keywords: Magneto-hydrodynamics
simulation
Issue Date: 31-Aug-2011
Abstract: Magnetohydrodynamics (MHD) studies the dynamics of an electrically conducting fluid under the influence of a magnetic field. Many astrophysical phenomena are related to MHD, and computer simulations are used to model these dynamics. In this thesis, we conduct MHD simulations of non-radiative black hole accretion as well as fast magnetic reconnection. By performing large scale three dimensional parallel MHD simulations on supercomputers and using a deformed-mesh algorithm, we were able to conduct very high dynamical range simulations of black hole accretion of Sgr A* at the Galactic Center. We find a generic set of solutions, and make specific predictions for currently feasible observations of rotation measure (RM). The magnetized accretion flow is subsonic and lacks outward convection flux, making the accretion rate very small and having a density slope of around $-1$. There is no tendency for the flows to become rotationally supported, and the slow time variability of the RM is a key quantitative signature of this accretion flow. We also provide a constructive numerical example of fast magnetic reconnection in a three-dimensional periodic box. Reconnection is initiated by a strong, localized perturbation to the field lines and the solution is intrinsically three-dimensional. Approximately $30\%$ of the magnetic energy is released in an event which lasts about one Alfv\'en time, but only after a delay during which the field lines evolve into a critical configuration. In the co-moving frame of the reconnection regions, reconnection occurs through an X-like point, analogous to the Petschek reconnection. The dynamics appear to be driven by global flows rather than local processes. In addition to issues pertaining to physics, we present results on the acceleration of MHD simulations using heterogeneous computing systems \cite{shan2006heterogeneous}. We have implemented the MHD code on a variety of heterogeneous and multi-core architectures (multi-core x86, Cell, Nvidia and ATI GPU) using different languages (FORTRAN, C, Cell, CUDA and OpenCL). Initial performance results for these systems are presented, and we conclude that substantial gains in performance over traditional systems are possible. In particular, it is possible to extract a greater percentage of peak theoretical performance from some heterogeneous systems when compared to x86 architectures.
URI: http://hdl.handle.net/1807/29832
Appears in Collections:Doctoral

Files in This Item:

File Description SizeFormat
Pang_Bijia_201106_PhD_video.mp4Video20.88 MBMPEG
View/Open
Pang_Bijia_201106_PhD_thesis.pdfMain article6.59 MBAdobe PDF
View/Open

Items in T-Space are protected by copyright, with all rights reserved, unless otherwise indicated.

uoft