We are investigating phase transitions in simple models for ionic melts and solutions. These systems are extremely difficult to handle by either theoretical or numerical methods because of the presence of infinite-range ionic forces and strong counterion association.  Work by our group and others has generated estimates for the phase behavior of  "primitive models" for ionic systems consisting of charged hard spheres in a dielectric continuum. Our recent work has extended the charge asymmetry to ratios of 2000:1, appropriate for charged colloids.