We investigate by means of numerical simulations with nested Cartesian meshes the force exerted on a massive and luminous perturber moving at constant speed through a homogeneous and opaque gas, taking into account thermal diffusion in the gas and the radiative feedback from the perturber. The force arising from the release of energy into the ambient medium by the perturber, or heating force, is directed along the direction of motion and induces an acceleration of the perturber. Its value is compared to analytic estimates in the low and high Mach number limits, and found to match those accurately. In addition, the drag exerted on a non-luminous perturber significantly departs from the adiabatic expression when thermal diffusion is taken into account. In the limit of a vanishing velocity, this drag tends to a finite value which we determine using linear perturbation theory and corroborate with numerical simulations. The drag on a non-luminous perturber in a non-adiabatic gas therefore behaves like dry or solid friction. We work out the luminosity threshold to get a net acceleration of the perturber and find it to be generally much smaller than the luminosity of accreting low-mass planetary embryos embedded in a gaseous protoplanetary disc at a few astronomical units. We also present in some detail our implementation of nested meshes, which runs in parallel over several \emph{Graphics Processing Units} (GPUs).

We consider the gravitational force exerted on a point-like perturber of mass $M$ travelling within a uniform gaseous, opaque medium at constant velocity $V$. The perturber irradiates the surrounding gas with luminosity $L$. The diffusion of the heat released is modelled with a uniform thermal diffusivity $\chi$. Using linear perturbation theory, we show that the force exerted by the perturbed gas on the perturber differs from the force without radiation (or standard dynamical friction). Hot, underdense gas trails the mass, which gives rise to a new force component, the heating force, with direction $+V$, thus opposed to the standard dynamical friction. In the limit of low Mach numbers, the heating force has expression $F_\mathrm{heat}=\gamma(\gamma-1)GML/(2\chi c_s^2)$, $c_s$ being the sound speed and $\gamma$ the ratio of specific heats. In the limit of large Mach numbers, $F_\mathrm{heat}=(\gamma-1)GML/(\chi V^2)f(r_\mathrm{min}V/4\chi)$, where $f$ is a function that diverges logarithmically as $r_\mathrm{min}$ tends to zero. Remarkably, the force in the low Mach number limit does not depend on the velocity. The equilibrium speed, when it exists, is set by the cancellation of the standard dynamical friction and heating force. In the low Mach number limit, it scales with the luminosity to mass ratio of the perturber. Using the above results suggests that Mars- to Earth-sized planetary embryos heated by accretion in a gaseous protoplanetary disc should have eccentricities and inclinations that amount to a sizeable fraction of the disc's aspect ratio, for conditions thought to prevail at a few astronomical units.