Densely packed DNA arrays exhibit hexagonal and orthorhombic local packings, as well as a weakly first order transition between them. While we have some understanding of the interactions between DNA molecules in aqueous ionic solutions, the structural details of its ordered phases and the mechanism governing the respective phase transitions between them remains less well understood. Since at high DNA densities, i.e., small interaxial spacings, one can neither neglect the atomic details of the interacting macromolecular surfaces nor the atomic details of the intervening ionic solution, the atomistic resolution is a sine qua non to properly describe and analyze the interactions between DNA molecules. In fact, in order to properly understand the details of the observed osmotic equation of state, one needs to implement multiple levels of organization, spanning the range from the molecular order of DNA itself, the possible ordering of counterions, and then all the way to the induced molecular ordering of the aqueous solvent, all coupled together by electrostatic, steric, thermal and direct hydrogen-bonding interactions. Multiscale simulations therefore appear as singularly suited to connect the microscopic details of this system with its macroscopic thermodynamic behavior. We review the details of the simulation of dense atomistically resolved DNA arrays with different packing symmetries and the ensuing osmotic equation of state obtained by enclosing a DNA array in a monovalent salt and multivalent (spermidine) counterions within a solvent permeable membrane, mimicking the behavior of DNA arrays subjected to external osmotic stress. By varying the DNA density, the local packing symmetry, and the counterion type, we are able to analyze the osmotic equation of state together with the full structural characterization of the DNA subphase, the counterion distribution and the solvent structural order in terms of its different order parameters and consequently identify the most important contribution to the DNA-DNA interactions at high DNA densities.