Direct numerical simulations and mean-field theory are used to model reactive front propagation in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken into account to estimate the front speed in cases when the Damkohler number is large. This effect is found to saturate the front speed to values comparable with the speed of the turbulent motions. By comparing with direct numerical simulations, it is found that the effective correlation time is much shorter than for non-reacting flows. The nonlinearity of the reaction term is found to make the front speed slightly faster.

A high-order scheme for direct numerical simulations of turbulent combustion is discussed. Its implementation in the massively parallel and publicly available Pencil Code is validated with the focus on hydrogen combustion. Ignition delay times (0D) and laminar flame velocities (1D) are calculated and compared with results from the commercially available Chemkin code. The scheme is verified to be fifth order in space. Upon doubling the resolution, a 32-fold increase in the accuracy of the flame front is demonstrated. Finally, also turbulent and spherical flame front velocities are calculated and the implementation of the non-reflecting so-called Navier-Stokes Characteristic Boundary Condition is validated in all three directions.

In an attempt to model the accretion onto a neutron star in low-mass X-ray binaries, we present two-dimensional hydrodynamical models of the gas flow in close vicinity of the stellar surface. First we consider a gas pressure dominated case, assuming that the star is non-rotating. For the stellar mass we take $M_{\rm star}=1.4 \times 10^{-2} \msun$ and for the gas temperature $T=5 \times 10^6$ K. Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of $T\simeq 10^8$ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.

The silicate carbon star V778 Cyg is a source of 22 GHz water maser emission which was recently resolved by MERLIN. Observations revealed an elongated S-like structure along which the velocities of the maser features show a linear dependence on the impact parameter. This is consistent with a doubly-warped m=2 disc observed edge-on. Water masers and silicate dust emission (detected by IRAS and ISO) have a common origin in O-rich material and are likely to be co-located in the disc. We propose a detailed self-consistent model of a masing gas-dust disc around a companion to the carbon star in a binary system, which allows us to estimate the companion mass of 1.7 +- 0.1 M_sun, the disc radius of 40 +-3 AU and the distance between companions of about 80 AU. Using a dust-gas coupling model for water masing, we calculate the maser power self-consistently, accounting for both the gas and the dust energy balances. Comparing the simulation results with the observational data, we deduce the main physical parameters of the masing disc, such as the gas and dust temperatures and their densities. We also present an analysis of the stability of the disc.

Various scenarios have been proposed to explain the presence of silicate features associated with carbon stars, such as V778 Cyg. We have attempted to constrain these theories by means of mapping water maser mission from V778 Cyg. The 22 GHz water maser emission from this star has been mapped using MERLIN with an astrometric accuracy of 25 mas. The spatially- and kinematically-resolved maser complex is displaced by about 190 mas from the position of the C-star as measured 10 years earlier using Tycho. Our simulations and analysis of available data show that this position difference is unlikely to be due to proper motion if V778 Cyg is at the assumed distance of 1.4 kpc. The maser components seem to form a distorted S-shaped structure extended over 18 mas with a clear velocity gradient. We propose a model which explains the observed water maser structure as an O-rich warped disk around a companion of the C-star in V 778 binary system, which is seen almost edge-on. Analysis of observational data, especially those obtained with MERLIN, suggests that V778 Cyg (and, by implication, other silicate carbon stars) are binary systems composed of a C-rich star and a companion which stores circumstellar O-rich material.

We study the conditions for operation of the 22 GHz ortho-water maser in a dusty medium near late-type stars. The main physical processes, such as exchange of energy between dust and gas in the radiation field of a star, radiative cooling by water molecules and pumping of water masers are described self-consistently. We show that the presence of dust grains of various types (or of one type with size distribution) strongly affects the maser action. The pumping mechanism based on the presence of the dust of different optical properties is able to explain water masers in the silicate carbon star V778 Cyg. However, the masers in the winds from asymptotic giant branch stars require an additional source of heating, for instance due to the dust drift through the gas.

We study in details a pumping mechanism for the lambda=1.35 cm maser transition 6_16 -> 5_23 in ortho-water based on the difference between gas and dust temperatures. The upper maser level is populated radiatively through 4_14 -> 5_05 and 5_05 -> 6_16 transitions. The heat sink is realized by absorbing the 45 mum photons, corresponding to the 5_23 -> 4_14 transition, by cold dust. We compute the inversion of maser level populations in the optically thick medium as a function of the hydrogen concentration, the gas-to-dust mass ratio, and the difference between the gas and the dust temperatures. The main results of numerical simulations are interpreted in terms of a simplified four-level model. We show that the maser strength depends mostly on the product of hydrogen concentration and the dust-to-water mass ratio but not on the size distribution of the dust particles or their type. We also suggest approximate formulae that describe accurately the inversion and can be used for fast calculations of the maser luminosity. Depending on the gas temperature, the maximum maser luminosity is reached when the water concentration N_water ~ 10^6-10^7 cm^-3 times the dust-to-hydrogen mass ratio, and the inversion completely disappears at density just an order of magnitude larger. For the dust temperature of 130 K, the 6_16 -> 5_23 transition becomes inverted already at the temperature difference of Delta T ~1 K, while other possible masing transitions require a larger Delta T > 30 K. We identify the region of the parameter space where other ortho- and para-water masing transitions can appear.