HOME   PROFILE   RESEARCH   CODES   PHOTOS   CONTACT

Simulation codes

All simulation codes I used in my work are publicly available in my github.

Phantom SPH (fork)

The SPH code I mainly use is Phantom, developed by Daniel Price and collaborators. Documentations can be found here. I work with my own fork repository. Several new physics modules have been implemented for my research. Major ones include:

For the SPH-MCRT RHD scheme -

• photoionize_cmi.F90
• kdtree_cmi.f90
• hnode_cmi.f90
• heating_cooling_cmi.f90
• utils_cmi.f90

For cloud-scale supernovae injection -

• inject_sne_sphng.f90
• cooling.F90

For cluster potentials -

• extern_starcluster.f90

CMacIonize (fork)

The grid-based Monte Carlo radiative transfer (MCRT) code adopted in our RHD scheme is CMacIonize, developed by Bert Vandenbroucke and collaborators. I work with my fork repository though only the mapping routines have been modified.

Instructions on how to couple Phantom to CMacIonize can be found here.

NbodyAcc

NbodyAcc is a recently-developed simple N-body code, dedicated for modelling massive binary accretion in clusters. Individual point masses can be treated as binaries; we track their accretion radii, internal angular momenta and separations. The code employs a 4th-order Runge-Kutta integrator, with King models for cluster potentials and incorporates turbulence. Code structure is similar to Phantom.

SPHNG

SPHNG is the precursor of Phantom. It was developed by Willy Benz, then passed to Ian Bonnell, then to Matthew Bate. I keep a copy of SPHNG from James Wurster and individual modules from Ian Bonnell. This code is yet to be adopted in my work but it was occationally used as a reference.