COSE$\nu$: A Collective Oscillation Simulation Engine for Neutrinos

Understand

We introduce the implementation details of the simulation code \cosenu, which numerically solves a set of non-linear partial differential equations that govern the dynamics of neutrino collective flavor conversions.

• We systematically provide the details of both the finite difference method supported by Kreiss-Oliger dissipation and the finite volume method with seventh order weighted essentially non-oscillatory scheme.
• To ensure the reliability of the code, we perform the comparison of the simulation results with theoretically obtainable solutions.
• In order to understand and characterize the error accumulation behavior of the implementations when neutrino self-interactions are switched on, we also analyze the evolution of the deviation of the conserved quantities for different values of simulation parameters.