PhD research at UC Riverside

My area of research at UCR is galactic formation and evolution using numerical simulations. This is a very rich topic that I have been interested in since my first astrophysics course as an undergraduate.  My advisor is Gabriela Canalizo.

In order to simulate a galaxy, many phenomena must be included because galaxies are much more than simply gravitationally-bound systems of stars.  For instance, a good simulation must include a consistent simulation of the gas and dust in the galaxy so that star formation can be properly taken into account.  Stellar population evolution must be included because different types of stars interact with the interstellar medium differently at various stages in their evolution.  The interaction between AGN and the host galaxies must be included because AGN perturb the galactic dust and gas and thus affect star formation.  Features of the galactic environment must also be included because gravitational interactions, in-flowing and out-flowing gas, and mergers between galaxies are important in the evolution process.  Finally, cosmological considerations must be included in a good model; in particular, the average density, opacity, temperature, and rate of expansion of the universe can influence galactic evolution.

At the moment, I am interested in learning how different theories of gravitation impact the simulations.  All of the simulations, of which I am aware, use Newtonian gravitation, which is a causality-violating theory because it assumes that gravitational effects propagate infinitely fast.  Galaxies are sufficiently large and stars move sufficiently fast that the causality violation may be an important issue. Theories of gravitation that are consistent with causality contain extra effects (such as velocity-dependent forces) which are not reproduced by Newtonian gravity and thus are not included in the current simulations of galaxies.  Although there are theoretical arguments which show that the influence of these post-Newtonian effects should be small, I am interested in comparing simulations with and without causality violation to see whether the differences are significant in comparison to the typical numerical errors in the gravity solvers used by state-of-the-art simulation codes. I'm particularly interested in performing this comparison using cluster-scale simulations. Before tackling that problem though, I am gaining proficiency running and analyzing GADGET-2 simulations.  My doctoral research focuses on more mundane, but important, details of high-resolution simulations using a version of GADGET coupled with the visualization code, Sunrise, and observational astronomy tools.

Master's research at George Mason University

During my time at GMU, I worked in the Space Weather group. My advisor was Robert Weigel. About half of my time was spent learning to do research and acquiring the background knowledge required for the project that I would eventually work on.  I learned UNIX/Linux system administration, shell scripting, Apache web server configuration, $\LaTeX$, basic fluid dynamics, plasma physics, and MHD.  I also took Rainald Löhner's computational fluid dynamics course where I learned about unstructured mesh generation, methods of solution for PDEs on unstructured meshes, and techniques of high-performance computing.

The focus of my research was inner-magnetosphere modeling.  The original intent was to acquire the source code to the Rice Convection Model (RCM), which is the standard code for modeling the conditions in the inner magnetosphere.  The plan was to modify the code such that it would be able handle a more realistic (less symmetric) geometric configuration so that the seasonal dependence of the model could be tested.  My advisor and I were unable to acquire a copy of the RCM.  It eventually became clear that the Authors at Rice University were not anxious to share, so I had to write the entire code myself.  I read the published descriptions of the RCM algorithm and I wrote a prototype in MATLAB.  I was able to complete the main structure of the program before I graduated. My description of the project, OpenConvection, can be found here. It is now being used mainly for educational purposes.

Some plots from OpenConvection