Flux-Limited Radiative Diffusion Module Applicable to Protoplanetary Disks

Previous numerical simulations on planet-disk interactions revealed a lot of interesting phenomena including the planetary migration and the formation of many sub-structures inside the disks. However, these simulations were limited to an isothermal or adiabatic equation of state which does not account for various heating and cooling processes in the disks. Recent studies showed that the behavior of the planet-disk interaction can be significantly influenced by the disk thermodynamics. We develop a radiative diffusion module based on the two-temperature flux-limited diffusion approximation accounting for viscous heating and the accretion feedback. In this presentation, we describe our radiative diffusion solver, present some test results, and discuss potential applications of the module to planet-disk interactions,