We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin and Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we introduce the multifluid extension of Relaxing TVD scheme and present a test case of dust migration in a two-fluid disk composed of gas and dust. We demonstrate that due to the difference in azimuthal velocities of gas and dust and the drag force acting on both components dust drifts towards maxima of gas pressure distribution.

We examine the predictions of the core accretion - gas capture model concerning the efficiency of planet formation around stars with various masses. First, we follow the evolution of gas and solids from the moment when all solids are in the form of small grains to the stage when most of them are in the form of planetesimals. We show that the surface density of the planetesimal swarm tends to be higher around less massive stars. Then, we derive the minimum surface density of the planetesimal swarm required for the formation of a giant planet both in a numerical and in an approximate analytical approach. We combine these results by calculating a set of representative disk models characterized by different masses, sizes, and metallicities, and by estimating their capability of forming giant planets. Our results show that the set of protoplanetary disks capable of giant planet formation is larger for less massive stars. Provided that the distribution of initial disk parameters does not depend too strongly on the mass of the central star, we predict that the percentage of stars with giant planets should increase with decreasing stellar mass. Furthermore, we identify the radial redistribution of solids during the formation of planetesimal swarms as the key element in explaining these effects.

We have investigated the problem of the distribution of both masses and orbital radii of planets resulting from the gas-accretion, gas-capture model. First we followed the evolution of gas and solids from the moment where all solids are in the form of small grains to the stage when most of them are in the form of planetesimals for a set of different initial masses and sizes of protoplanetary disks. Based on that we performed Monte Carlo calculations describing the formation of giant places at different locations. We included the effects of type II migration and growth of the mass of the planet after the gap opened. We discuss how these effects influence the final distribution of giant planets. We show that when the giant planets are not able to migrate or grow in mass after the gap opens, their distribution is mainly determined by the properties of the gaseous disk. However, with those two effects included, reproducing the parameters of the gaseous disks from the distribution of planets becomes difficult. We also checked the roles of both the material of which the solids consist and the mass of the central star. The main result is that, in disks around less massive stars, giant planets at the given location tend to be less massive. At the same time, the giant planets with the given mass tend to form closer to the less massive stars.

We present the first results from simulations of processes leading to planet formation in protoplanetary disks with different metallicities. For a given metallicity, we construct a two-dimensional grid of disk models with different initial masses and radii ($M_0$, $R_0$). For each disk, we follow the evolution of gas and solids from an early evolutionary stage, when all solids are in the form of small dust grains, to the stage when most solids have condensed into planetesimals. Then, based on the core accretion - gas capture scenario, we estimate the planet-bearing capability of the environment defined by the final planetesimal swarm and the still evolving gaseous component of the disk. We define the probability of planet-formation, $P_p$, as the normalized fractional area in the ($M_0$, $\log R_0$) plane populated by disks that have formed planets inside 5 AU. With such a definition, and under the assumption that the population of planets discovered at $R$ $<$ 5 AU is not significantly contaminated by planets that have migrated from $R$ $>$ 5 AU, our results agree fairly well with the observed dependence between the probability that a star harbors a planet and the star's metal content. The agreement holds for the disk viscosity parameter $\alpha$ ranging from $10^{-3}$ to $10^{-2}$, and it becomes much poorer when the redistribution of solids relative to the gas is not allowed for during the evolution of model disks.

We have calculated an evolution of protoplanetary disk from an extensive set of initial conditions using a time-dependent model capable of simultaneously keeping track of the global evolution of gas and water-ice. A number of simplifications and idealizations allows for an embodiment of gas-particle coupling, coagulation, sedimentation, and evaporation/condensation processes. We have shown that, when the evolution of ice is explicitly included, the location of the snowline has to be calculated directly as the inner edge of the region where ice is present and not as the radius where disk's temperature equals the evaporation temperature of water-ice. The final location of the snowline is set by an interplay between all involved processes and is farther from the star than implied by the location of the evaporation temperature radius. The evolution process naturally leads to an order of magnitude enhancement in surface density of icy material.

Formation of planets in the 47 UMa system is followed in an evolving protoplanetary disk composed of gas and solids. The evolution of the disk is calculated from an early stage, when all solids, assumed to be high-temperature silicates, are in the dust form, to the stage when most solids are locked in planetesimals. The simulation of planetary evolution starts with a solid embryo of ~1 Earth mass, and proceeds according to the core accretion -- gas capture model. Orbital parameters are kept constant, and it is assumed that the environment of each planet is not perturbed by the second planet. It is found that conditions suitable for both planets to form within several Myr are easily created, and maintained throughout the formation time, in disks with $\alpha \approx 0.01$. In such disks, a planet of 2.6 Jupiter masses (the minimum for the inner planet of the 47 UMa system) may be formed at 2.1 AU from the star in \~3 Myr, while a planet of 0.89 Jupiter masses (the minimum for the outer planet) may be formed at 3.95 AU from the star in about the same time. The formation of planets is possible as a result of a significant enhancement of the surface density of solids between 1.0 and 4.0 AU, which results from the evolution of a disk with an initially uniform gas-to-dust ratio of 167 and an initial radius of 40 AU.

We numerically test the mode-coupling model (Lubow 1991a) of tidal instability in SU UMa systems. So far, all numercial models confirming it have been based on SPH codes and isothermal equation of state. In our paper we present Eulerian models, using both isothermal approximation and the full energy equation. We also investigate influence of different ways of mass transfer. While isothermal models behave similarly to SPH simulations, the behaviour of models with full energy equation is quite different, and the mode-coupling model is not confirmed in this case.