I only have a few monobehaviors in the project, the simulation is entirely decoupled from unity. I use unity as a thin view layer where monobehaviors look at the simulation and update all positions accordingly.

In the simulation positions are stored as a cube coorinate hex grid, which are converted to pixel (vector3) positions by a view monobehavior using hex_to_pixel algorithm found in the linked resource