We present new two-dimensional radiation hydrodynamic simulations of supernova shock breakout from red supergiants using the
CASTRO code. Our progenitors are 20 and 25 M
⊙ solar-metallicity stars evolved from the zero-age main sequence with
MESA and exploded in one dimension using
FLASH. We consider a range of circumstellar media (CSM) produced by stellar winds to investigate how pre-explosion mass-loss affects shock breakout. The multigroup flux-limited diffusion scheme in
CASTRO captures the interaction between the explosion shock, its radiation precursor, and the surrounding CSM. We find that strong radiation precursors, generated by radiation leakage behind the shock, can drive fluid instabilities and move the effective photosphere outward before the shock reaches the stellar surface. The resulting breakout emissions reach peak luminosities of
∼1044 erg s
−1 with full-width half-maximum durations of 1-3 hr, which are much dimmer and longer than those from blue supergiants. The light-curve colors gradually evolve from blue to red after the peak. The 25 M
⊙ model with explosion energy
E∼1.69×1051 erg produces
∼10-30\% higher maximum luminosity than the 20 M
⊙ model with
E∼1.09×1051 erg. The dense CSM further extends the breakout rise time by increasing the photon diffusion. These results provide new constraints on red supergiant atmospheres and mass-loss histories prior to core collapse.