  We report on our recent attempt of quantitative modeling of the Chiral
Magnetic Effect (CME) in heavy-ion collisions. We perform 3+1 dimensional
anomalous hydrodynamic simulations on an event-by-event basis, with
constitutive equations that contain the anomaly-induced effects. We also
develop a model of the initial condition for the axial charge density that
captures the statistical nature of random chirality imbalances created by the
color flux tubes. Basing on the event-by-event hydrodynamic simulations for
hundreds of thousands of collisions, we calculate the correlation functions
that are measured in experiments, and discuss how the anomalous transport
affects these observables.
