We describe a simple and fast method for simulating observed images of ion Coulomb crystals. In doing so, cold elastic collisions between Coulomb crystals and virtual very light atoms are implemented in a molecular dynamics (MD) simulation code. Such an approach reproduces the observed images of Coulomb crystals by obtaining density plots of the statistics of existence of each ion. The simple method has the advantage of short computing time in comparison with previous calculation methods. As a demonstration of the simulation, the formation of a planar Coulomb crystal with a small number of ions has been investigated in detail in a linear ion trap both experimentally and by simulation. However, also large Coulomb crystals including up to 1400 ions have been photographed and simulated to extract the secular temperature and the number of ions. For mediumsized crystals, a comparison between experiments and calculations has been performed. Moreover, an MD simulation of the sympathetic cooling of small molecular ions was performed in order to test the possibility of extracting the temperature and the number of refrigerated molecular ions from crystal images of laser-cooled ions. Such information is basic to studying ultracold ion-molecule reactions using ion Coulomb crystals including sympathetically cooled molecular ions.