Massive elliptical galaxies harbor large amounts of hot gas (T ≳ 106 K) in their interstellar medium (ISM) but are typically quiescent in star formation. The jets of active galactic nuclei (AGNs) and Type Ia supernovae (SNe Ia) inject energy into the ISM, which offsets its radiative losses and keeps it hot. SNe Ia deposit their energy locally within the galaxy compared to the larger few ×10 kiloparsec-scale AGN jets. In this study, we perform high-resolution (5123) hydrodynamic simulations of a local (1 kpc3) density-stratified patch of the ISM of massive galaxies. We include radiative cooling and shell-averaged volume heating, as well as randomly exploding SN Ia. We study the effect of different fractions of supernova (SN) heating (with respect to the net cooling rate), different initial ISM density/entropy (which controls the growth time t ti of the thermal instability), and different degrees of stratification (which affect the freefall time t ff). We find that SNe Ia drive predominantly compressive turbulence in the ISM with a velocity dispersion of σ v up to 40 km s−1 and logarithmic density dispersion of σ s ∼ 0.2–0.4. These fluctuations trigger multiphase condensation in regions of the ISM, where min(tti)/tff≲0.6exp(6σs) , in agreement with theoretical expectations that large density fluctuations efficiently trigger multiphase gas formation. Since the SN Ia rate is not self-adjusting, when the net cooling drops below the net heating rate, SNe Ia drive a hot wind which sweeps out most of the mass in our local model. Global simulations are required to assess the ultimate fate of this gas.