A numerical hydrodynamical model for the evolution of spherically symmetric collapsing clouds, designed for the calculation of the thermal structure of these objects in both the prestellar and protostellar stages of their evolution, is presented. Distinctive features of the model include the possibility of independently describing the temperatures of the gas and dust, which is extremely important when calculating the thermal structure of prestellar and protostellar clouds, and the account of the radiation flux from the central protostar. This model is used to compare the theoretical density and temperature distributions with observations for nearby sites of star formation obtained with the Herschel Space Observatory. Application of the diffusion approximation with a flux limiter describes well the radial density and temperature distributions in protostellar clouds. However, significant differences between the model and observational density profiles were found for prestellar stages, suggesting the presence of appreciable deviations from equilibrium in the prestellar clouds. An approximate method for calculating the thermal structure of a cloud based on the adaptive $\tau$-approximation is presented. Application of the $\tau$-approximation yields good agreement with the diffusion approximation for the prestellar phase, but produces appreciable discrepancies for the protostellar phase, when the thermal structure of the accreting envelope is determined by the radiation of the protostar.

Speckle interferometry of the young binary system RW Aur was performed with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences using filters with central wavelengths of 550 nm and 800 nm and pass-band halfwidths of 20 nm and 100 nm, respectively. The angular separation of the binary components was 1.448"{\pm}0.005 and the position angle of the system was 255.9{\deg}{\pm}0.3{\deg} at the observation epoch (JD 2 454 255.9). We find using published data that these values have been changing with mean rates of +0.002"/yr and +0.02{\deg}/yr, respectively, over the past 70 years. This implies that the direction of the orbital motion of the binary system is opposite to the direction of the disk rotation in RW Aur A. We propose a physical model to explain the formation of circumstellar accretion disks rotating in the reverse direction relative to young binary stars surrounded by protoplanetary disks. Our model can explain the characteristic features of the matter flow in RW Aur A: the high accretion rate, small size of the disk around the massive component, and reverse direction of rotation.

The structure of a circumstellar envelope around a young binary T Tauri star is considered. The supersonic orbital motion of the system components in the envelope gas leads to the formation of bow shocks around the star. Two- and three-dimensional numerical modeling indicates an important role of these shocks in the formation of the structure of the circumbinary envelope. In particular, for systems with circular orbits, the size of the central region of the envelope that is not filled with matter (the gap) is essentially determined by the parameters of the bow shocks. These modeling results are supported by comparisons of the obtained estimates for the gap parameters with observations.

A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, and external background, and the theoretical distribution of the dust temperature inside the cloud is determined. The intensity distributions at millimeter and IR wavelengths are computed and quantitatively compared with observational data. The best-fit model parameters are determined using a genetic minimization algorithm, which makes it possible to reveal the ranges of parameter degeneracy as well. The method is illustrated by modeling the structure of the two infrared dark clouds IRDC-320.27+029 (P2) and IRDC-321.73+005 (P2). The derived density and temperature distributions can be used to model the chemical structure and spectral maps in molecular lines.