We review recent developments in the field of chemodynamical simulations of
elliptical galaxies, highlighting (in an admittedly biased fashion)
the work conducted with our cosmological N-body/SPH code GCD+.
We have demonstrated previously the
recovery of several primary integrated
early-type system scaling
relations (e.g. colour- magnitude relation,
LX − TX − [Fe/H]X)
when employing a phenomenological AGN heating scheme in conjunction
with a self-consistent treatment of star formation,
supernovae feedback, radiative cooling, chemical enrichment, and stellar/X-ray
population synthesis. Here we emphasise characteristics derived
from the full spatial information contained within the simulated
dataset, including stellar and coronal morphologies, metallicity distribution
functions, and abundance gradients.