Building three-dimensional giant stellar models for common envelope simulations

We build a three-dimensional (3D) red supergiant (RSG) stellar model for common envelope evolution (CEE) simulations by transporting a 1D stellar model to a 3D numerical grid, mimicking core nuclear power by depositing energy to an inner shell, and mimicking stellar emission by cooling grid cells with densities below the photospheric density. We do not relax the model; rather, we let it perform its natural pulsation. We find that when we mimic photospheric emission by cooling low-density grid cells, the oscillations slowly decay on a time scale much longer than in the absence of photospheric cooling. When we mimic both nuclear energy production, by depositing the stellar luminosity in an inner shell above the inert core of the stellar model, and the photospheric cooling, the oscillations do not decay and their amplitude slowly increases with time. The main pulsational period is about 1 year, comparable to the stellar dynamical time, suggesting a fundamental radial pulsation mode. The non-spherical structure of the stellar model and rapid low-amplitude temporal variations in the average stellar radius testify to the presence of non-radial oscillation modes on top of the fundamental radial mode. We also obtain vigorous convection, as RSG stars have. We conclude that the best way of preparing a giant star to simulate CEE and grazing-envelope evolution is to deposit energy with the stellar luminosity in an inner shell, and to cool the outer low-density numerical shell. There is no need to relax the model.