We extend the applicability of the hydrodynamics, perturbative QCD and saturation -based EKRT (Eskola-Kajantie-Ruuskanen-Tuominen) framework for ultrarelativistic heavy-ion collisions to peripheral collisions by introducing dynamical freeze-out conditions. As a new ingredient compared to the previous EKRT computations we also introduce a nonzero bulk viscosity. We compute various hadronic observables and flow correlations, including normalized symmetric cumulants, mixed harmonic cumulants, and flow–transverse-momentum correlations, and compare them against measurements from the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) . We demonstrate that the inclusion of the dynamical freeze-out and bulk viscosity allows a better description of the measured flow coefficients in peripheral collisions and enables the use of an extended centrality range when constraining the properties of QCD matter in the future.