/FLASHcodeStellarCores

C++ code to extract and analyze stellar core regions in FLASH code AMR simulation data

Primary LanguageC++MIT LicenseMIT

Overview

Source repository for the C++ analysis code from an undergraduate research project in astrophysics focusing on the "Core Accretion" and "Competetive Accretion" theories of stellar formation in the context of FLASH code simulations of star formation in the early universe. This code was written in 2014 at UT-Austin, and was migrated from BitBucket in May 2018.

It's unlikely that this code will be functional out-of-the-box, but pieces of it may be helpful for someone who's trying to analyze checkpoint file mesh data from FLASH code radiation-hydrodynamic simulations with Adaptive Mesh Refinement and/or Lagrangian sink particles.

For now the header files in Include/ are probably the quickest way to get a sense of the code structure. Feel free to reach out and I'll do my best to answer any questions.

Stellar Cores

Parts of the code utilize libtourtre: A Countour Tree Library to build countour trees from the gravitatial potential field. Stellar cores are then defined as containing the cells mapping to the connected component directly "above" the node of the corresponding stellar sink particle.

CoreAnalyzer::altFindCoreRegion() uses a seeded region growing method to attempt to find topologically equivalent core regions in a different, more efficient manner.

In theory, these methods should have produced identical cell-to-region mappings, although in practice they produce slightly different results. They did both produce very low-mass stellar initial mass functions (IMFs) that were uncorrelated with the final stellar mass function, suggesting that competetive accretion was much more descriptive of the stellar formation processes in these simulations. This is in agreement with the findings of the 2015 paper based on the same simulations: Star formation in the first galaxies - III (linked below).

Papers

For more information about the simulations, and their scientific purpose and context, consult the following papers:

Before this code was written:

After this code was written:

Background on the core and competetive theories of star formation: