We describe our modelling of the radiatively cooling shocks and their thin shells with various numerical tools in different physical and calculational setups. We inspect structure of the dense shell, its formation and evolution, pointing out physical and numerical factors that sustain its shape and also may lead to instabilities. We have found that under certain physical conditions, the circular shaped shells show a strong bending instability and successive fragmentation on Cartesian grids soon after their formation, while remain almost unperturbed when simulated on polar meshes. We explain this by physical Rayleigh--Taylor like instabilities triggered by corrugation of the dense shell surfaces by numerical noise. Conditions for these instabilities follow from both the shell structure itself and from episodes of transient acceleration during re-establishing of dynamical pressure balance after sudden radiative cooling onset. They are also easily excited by physical perturbations of the ambient medium. The widely mentioned Non-linear Thin Shell Instability, in contrast, in tests with physical perturbations is shown to have only limited chances to develop in real radiative shocks, as it seems to require a special spatial arrangement of fluctuations to be excited efficiently. The described phenomena also set new requirements on further simulations of the radiatively cooling shocks in order to be physically correct and free of numerical artefacts.

We study thermal emission from circumstellar structures heated by gamma-ray burst (GRB) radiation and ejecta and calculate its contribution to GRB optical and X-ray afterglows using the modified radiation hydro-code small STELLA. It is shown that thermal emission originating in heated dense shells around the GRB progenitor star can reproduce X-ray plateaus (like observed in GRB 050904, 070110) as well as deviations from a power law fading observed in optical afterglows of some GRBs (e.g. 020124, 030328, 030429X, 050904). Thermal radiation pressure in the heated circumburst shell dominates the gas pressure, producing rapid expansion of matter similar to supenova-like explosions close to opacity or radiation flux density jumps in the circumburst medium. This phenomenon can be responsible for so-called supernova bumps in optical afterglows of several GRBs. Such a `quasi-supernova' suggests interpretation of the GRB-SN connection which does not directly involve the explosion of the GRB progenitor star.

We present methodology and results of numerical modeling of the interaction of GRB prompt emission with the circumburst medium using a modified version of the multi-group radiation hydrocode STELLA. The modification includes the nonstationary photoionization, the photoionization heating and the Compton heating along with the hydrodynamics and radiation transfer. The lightcurves and spectra of the outcoming gamma-ray, X-ray and optical emission are calculated for a set of models (shells) of the circumburst environment, which differ in dimensions, density, density profile, composition, temperature. In some cases total bolometric and optical luminosities can reach 10^47 and 10^43 erg/s respectively. These effects can be responsible for irregularities which are seen on lightcurves of some GRB's X-ray and optical afterglows.