A unified chiral mean field approach is presented for QCD thermodynamics in a wide range of temperatures and densities. The model simultaneously gives a satisfactory description of lattice QCD thermodynamics and fulfills nuclear matter and astrophysical constraints. The resulting equation of state can be incorporated in relativistic fluid-dynamical simulations of heavy-ion collisions and neutron stars mergers. Access to different regions of the QCD phase diagram can be obtained in simulations of heavy-ion data and observations of neutron star mergers.