Not yet reviewed by Pith; the record is open.

Conservation laws play a crucial role in the modeling of heavy-ion collisions, including the those for charges such as baryon number (B), strangeness (S), and electric charge (Q). In this study, we present a new 2+1 relativistic viscous hydrodynamic code called CCAKE which uses the Smoothed Particle Hydrodynamics (SPH) formalism to locally conserve BSQ charges, together with an extended description of the multi-dimensional equation of state (EoS) obtained from lattice Quantum Chromodynamics. Initial conditions for CCAKE are supplied by the ICCING model, which samples gluon splittings into quark anti-quark pairs to generate the initial BSQ charge distributions. We study correlations between the BSQ charges and find that local BSQ fluctuations remain finite during the evolution, with corresponding chemical potentials of ($\sim100$--$200 \,\rm MeV$) at freeze-out. We find that our framework produces reasonable multiplicities of identified particles and that ICCING has no significant effect on the collective flow of all charged particles nor of identified particles when only one particle of interest is considered. However, we show specifically for Pb+Pb collisions at the LHC $\sqrt{s_{NN}}=5.02$ TeV that ICCING does have an effect on collective flow of identified particles if two particles of interest are considered.