In this contribution, we review our current and previous investigations of total, proton, and neutron particle number densities and the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions at beam energy below 800 MeV/nucleon. Furthermore, the effects of the Coulomb interaction and the dependence on a system size in three representative symmetric and asymmetric Ca, Sn and Pb systems are studied. The Boltzmann-Uhlenbeck-Uehling (pBUU) transport and Time-Dependent Hartree-Fock (TDHF) models, employing the SVbas, SkT3 and SVsym Skyrme interactions, are used in these simulations. We find that (i) the highest total densities predicted at Ebeam = 800 MeV/nucleon are on the order of ∼ 2.5ρ0 (ρ0 = 0.16 fm−3), (ii) the proton-neutron asymmetry for maximal densities, δ = (ρnmax−ρpmax)/(ρnmax+ρpmax) do not generally exceed the asymmetry in the initial state of the collision at all beam energies and tend to decrease during the reaction, and (iii) a significant portion of this asymmetry has its microscopic origin in Coulomb forces, masking the pure nuclear contribution. In addition, the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy, the impact of correlations in the reaction, and the time evolution of the proton and neutron density distributions in the plane transverse to the beam direction are illustrated and discussed.