High-mass X-ray binaries belong to the brightest objects in the X-ray sky. They usually consist of a massive O or B star or a blue supergiant while the compact X-ray emitting component is a neutron star (NS) or a black hole. Intensive matter accretion onto the compact object can take place through different mechanisms: wind accretion, Roche-lobe overflow, or circumstellar disk. In our multi-dimensional models we perform numerical simulations of the accretion of matter onto a compact companion in case of Be/X-ray binaries. Using Bondi-Hoyle-Littleton approximation, we estimate the NS accretion rate. We determine the Be/X-ray binary disk hydrodynamic structure and compare its deviation from isolated Be stars’ disk. From the rate and morphology of the accretion flow and the X-ray luminosity we improve the estimate of the disk mass-loss rate. We also study the behavior of a binary system undergoing a supernova explosion, assuming a blue supergiant progenitor with an aspherical circumstellar environment.