We evaluate the torque acting on a gravitational perturber on a retrograde circular orbit in the midplane of a gaseous disk. We assume that the mass of this satellite is so low it weakly disturbs the disk (type I migration). The perturber may represent the companion of a binary system with a small mass ratio. We compare the results of hydrodynamical simulations with analytic predictions. Our two-dimensional (2D) simulations indicate that the torque acting on a perturber with softening radius $R_{\rm soft}$ can be accounted for by a scattering approach if $R_{\rm soft}<0.3H$, where $H$ is defined as the ratio between the sound speed and the angular velocity at the orbital radius of the perturber. For $R_{\rm soft}>0.3H$, the torque may present large and persistent oscillations, but the resultant time-averaged torque decreases rapidly with increasing $R_{\rm soft}/H$, in agreement with previous analytical studies. We then focus on the torque acting on small-size perturbers embedded in full three-dimensional (3D) disks and argue that the density waves propagating at distances $\lesssim H$ from the perturber contribute significantly to the torque because they transport angular momentum. We find a good agreement between the torque found in 3D simulations and analytical estimates based on ballistic orbits. We compare the radial migration timescales of prograde versus retrograde perturbers. For a certain range of the perturber's mass and aspect ratio of the disk, the radial migration timescale in the retrograde case may be appreciably shorter than in the prograde case. We also provide the smoothing length required in 2D simulations in order to account for 3D effects.

We quantify the effects of magnetic fields, cosmic rays and gas pressure on the rotational velocity of HI gas in the Milky Way, at galactic distances between Rsun and 2Rsun. The magnetic field is modelled by two components; a mainly azimuthal magnetic component and a small-scale tangled field. We construct a range of plausible axisymmetric models consistent with the strength of the total magnetic field as inferred from radio synchrotron data. In a realistic Galactic disk, the pressure by turbulent motions, cosmic rays and the tangled turbulent field provide radial support to the disk. Large-scale (ordered) magnetic fields may or may not provide support to the disk, depending on the local radial gradient of the azimuthal field. We show that for observationally constrained models, magnetic forces cannot appreciably alter the tangential velocity of HI gas within a galactic distance of 2Rsun.

We study the final equilibrium states of the Parker instability arising from an initially unstable cylindrical equilibrium configuration of gas in the presence of a radial gravitational field and a longitudinal magnetic field. The aim of this work is to compare the properties of the nonlinear final equilibria with those found in a system with Cartesian geometry. Maps of the density and magnetic field lines, when the strength of the gravitational field is constant, are given in both geometries. In the axisymmetric model, the magnetic field tends to expand in radius, forming magnetic arcades, while knots of gas are formed because the plasma drains radially and strangulates the magnetic field lines, leading to the formation of magnetic bottlenecks. We find that the magnetic buoyancy and the drainage of gas along field lines are less efficient under axial symmetry than in a Cartesian atmosphere. As a consequence, the column density enhancement arising in gas condensations in the axially-symmetric model is smaller than in Cartesian geometry. The magnetic-to-gas pressure ratio in the final state takes more extreme values in the Cartesian model. Models with non-uniform radial gravity are also discussed.

There exist many physical processes that may contribute to the driving of turbulence in galactic disks. Some of them could drive turbulence even in the absence of star formation. For example, hydrodynamic (HD) or magnetohydrodynamic (MHD) instabilities, frequent mergers of small satellite clumps, ram pressure, or infalling gas clouds. In this work we present numerical simulations to study the interaction of compact high velocity clouds (CHVC) with the outskirts of magnetized gaseous disks. With our numerical simulations we show that the rain of small HVCs onto the disk is a potential source of random motions in the outer parts of HI disks.

We know that the galactic magnetic field possesses a random component in addition to the mean uniform component, with comparable strength of the two components. This random component is considered to play important roles in the evolution of the interstellar medium (ISM). In this work we present numerical simulations associated with the interaction of the supersonic flows located at high latitude in our Galaxy (High Velocity Clouds, HVC) with the magnetized galactic ISM in order to study the effect that produces a random magnetic field in the evolution of this objects.

Coronal mass ejections (CMEs) are solar eruptions into interplanetary space of as much as a few billion tons of plasma, with embedded magnetic fields from the Sun's corona. These perturbations play a very important role in solar--terrestrial relations, in particular in the spaceweather. In this work we present some preliminary results of the software development at the Universidad Nacional Autonoma de Mexico to perform Remote MHD Numerical Simulations. This is done to study the evolution of the CMEs in the interplanetary medium through a Web-based interface and the results are store into a database. The new astrophysical computational tool is called the Mexican Virtual Solar Observatory (MVSO) and is aimed to create theoretical models that may be helpful in the interpretation of observational solar data.

High--resolution 2D--MHD numerical simulations have been carried out to investigate the effects of continuing infall of clumpy gas in extended HI galactic disks. Given a certain accretion rate, the response of the disk depends on its surface gas density and temperature. For Galactic conditions at a galactocentric distance of ~20 kpc, and for mass accretion rates consistent with current empirical and theoretical determinations in the Milky Way, the rain of compact high velocity clouds onto the disk can maintain transonic turbulent motions in the warm phase (~2500 K) of HI. Hence, the HI line width is expected to be ~6.5 km/s for a gas layer at 2500 K, if infall were the only mechanism of driving turbulence. Some statistical properties of the resulting forcing flow are shown in this Letter. The radial dependence of the gas velocity dispersion is also discussed.

We present 2D, ideal-MHD numerical simulations of the Parker instability in a multi-component warm disk model. The calculations were done using two numerical codes with different algorithms, TVD and ZEUS-3D. The outcome of the numerical experiments performed with both codes is very similar, and confirms the results of the linear analysis for the undular mode derived by Kim et al. (2000): the most unstable wavelength is about 3 kpc and its growth timescale is between 30-50 Myr (the growth rate is sensitive to the position of the upper boundary of the numerical grid). Thus, the time and length scales of this multicomponent disk model are substantially larger than those derived for thin disk models. We use three different types of perturbations, random, symmetric, and antisymmetric, to trigger the instability. The antisymmetric mode is dominant, and determines the minimum time for the onset of the nonlinear regime. The instability generates dense condensations and the final peak column density value in the antisymmetric case, as also derived by Kim et al. (2000), is about a factor of 3 larger than its initial value. These wavelengths and density enhancement factors indicate that the instability alone cannot be the main formation mechanism of giant molecular clouds in the general interstellar medium. The role of the instability in the formation of large-scale corrugations along spiral arms is briefly discussed.

A linear stability analysis of a multi-component and magnetized Galactic disk model is presented. The disk model uses the observed stratifications for the gas density and gravitational acceleration at the solar neighborhood and, in this sense, it can be called a realistic model. The distribution of the total gas pressure is defined by these observed stratifications, and the gaseous disk is assumed isothermal. The initial magnetic field is taken parallel to the disk, with a midplane value of 5 $\mu$G, and its stratification along the z-axis is derived from the condition of magnetohydrostatic equilibrium in an isothermal atmosphere. The resulting isothermal sound speed is $\sim 8.4$ km s$^{-1}$, similar to the velocity dispersion of the main gas components within 1.5 kpc from midplane. The thermal-to-magnetic pressure ratio decreases with $[z]$ and the warm model is Parker unstable. The dispersion relations show that the fastest growing mode has a wavelength of about 3 kpc, for both symmetric and antisymmetric perturbations, and the corresponding growth time scales are of about $3\times 10^7$ years. The structure of the final equilibrium stage is also derived, and we find that the midplane antisymmetric (MA) mode gathers more gas in the magnetic valleys. The resulting MA gas condensations have larger densities, and the column density enhancement is a factor of about 3 larger than the value of the initial stage. The unstable wavelengths and growth times for the multi-component disk model are substantially larger than those of a thin disk model, and some of the implications of these results are discussed.

We present two-dimensional MHD numerical simulations for the interaction of high-velocity clouds with both magnetic and non-magnetic Galactic thick gaseous disks. For the magnetic models, the initial magnetic field is oriented parallel to the disk, and we consider two different field topologies (with and without tension effects): parallel and perpendicular to the plane of motion of the clouds. The impinging clouds move in oblique trajectories and fall toward the central disk with different initial velocities. The $B$-field lines are distorted and compressed during the collision, increasing the field pressure and tension. This prevents the cloud material from penetrating into the disk, and can even transform a high-velocity inflow into an outflow, moving away from the disk. The perturbation creates a complex, turbulent, pattern of MHD waves that are able to traverse the disk of the Galaxy, and induce oscillations on both sides of the plane. Thus, the magnetic field efficiently transmits the perturbation over a large volume, but also acts like a shield that inhibits the mass exchange between the halo and the disk. For the non-magnetized cases, we also uncover some novel features: the evolution of the shocked layer generates a tail that oscillates, creating vorticity and turbulent flows along its trajectory.