In Weyl semimetals the application of parallel electric and magnetic fields leads to valley polarization -- an occupation disbalance of valleys of opposite chirality -- a direct consequence of the chiral anomaly. In this work, we present numerical tools to explore such nonequilibrium effects in spatially confined three-dimensional systems with a variable disorder potential, giving exact solutions to leading order in the disorder potential and the applied electric field. Application to a Weyl-metal slab shows that valley polarization also occurs without an external magnetic field as an effect of chiral anomaly "trapping": Spatial confinement produces chiral bulk states, which enable the valley polarization in a similar way as the chiral states induced by a magnetic field. Despite its finite-size origin, the valley polarization can persist up to macroscopic length scales if the disorder potential is sufficiently long ranged, so that direct inter-valley scattering is suppressed and the relaxation then goes via the Fermi-arc surface states.
