We performed three-dimensional magnetohydrodynamic simulations to study the evolution of a supernova remnant (SNR) in a turbulent neutral atomic interstellar medium. The media used as background shares characteristics with the Solar neighbourhood and the SNR has mass and energy similar to those of a Type Ia object. Our initial conditions consist of dense clouds in a diluted medium, with the main difference between simulations being the average magnitude of the magnetic field. We measured amplifications of the magnetic energy of up to 34$\%$ and we generated synthetic maps that illustrate how the same object can show different apparent geometries and physical properties when observed through different lines of sight.

We present 3D hydrodynamical simulations of a precessing jet with a time-dependent ejection velocity or a time-dependent ejection density, interacting with a circumstellar medium given by a dense, anisotropic, and slow AGB wind, forming a torus. We explore a set of configurations with different values for the precession angle and number of ejections. The temporal evolution of these models is analised at times up to 1500 or 1800 yr. From our hydrodynamical models, we obtain position-velocity diagrams (PV diagrams) in the [NII]{\lambda}6583 line to be compared with high resolution observations of the planetary nebula H1-67. From spectral data this object shows high-velocity jets and a point-symmetric morphology. With our synthetic PV diagrams we show that a precessing jet with a time-dependent ejection velocity or a time-dependent ejection density reproduce the point-symmetric morphological structure for this nebula if the precession cone angle is larger than 30$^{\circ}$. Our synthetic PV diagrams can be used to understand how the S-like morphology, also presented by other planetary nebulae, is formed. For H1-67 we found a heliocentric velocity of -8.05 km s^-1 and height below the galactic plane of -451.6 pc.

We present a polarization study towards the supernova remnant (SNR) Puppis A based on original observations performed with the Australia Telescope Compact Array (ATCA). Based on the analysis of a feature detected outside the SNR shell (called `the tail' throughout the paper), it was possible to disentangle the emission with origin in Puppis A itself from that coming from the foreground Vela SNR. We found a very low polarization fraction, of about 3 percent on average. The upper limit of the magnetic field component parallel to the line of sight is estimated to be B$_\parallel \sim$ 20 $\mu$G. The statistical behavior of the magnetic vectors shows two preferential directions, almost perpendicular to each other, which are approximately aligned with the flat edges of Puppis A. A third, narrow peak oriented perpendicular to the Galactic plane suggests the existence of an interstellar magnetic field locally aligned in this direction. There is evidence that the magnetic vectors along the shell are aligned with the shock front direction. The low polarization fraction and the statistical behavior of the magnetic vectors are compatible with a scenario where the SNR evolves inside a stellar wind bubble with a box-like morphology, produced by the interaction of the different stellar winds, one of them magnetized, launched by the SN progenitor. This scenario can furthermore explain the morphology of Puppis A, rendering little support to the previously accepted picture which involved strong density gradients to explain the flat, eastern edge of the shell.

In the present work we have modeled the supernova remnant (SNR) G296.5+10, by means of 3D magnetohydrodynamics (MHD) simulations. This remnant belongs to the bilateral SNR group and has an additional striking feature: the rotation measure (RM) in its eastern and western parts are very different. In order to explain both the morphology observed in radio-continuum and the RM, we consider that the remnant expands into a medium shaped by the superposition of the magnetic field of the progenitor star with a constant Galactic magnetic field. We have also carried out a polarization study from our MHD results, obtaining synthetic maps of the linearly polarized intensity and the Stokes parameters. This study reveals that both the radio morphology and the reported RM for G$296.5+10$ can be explained if the quasi-parallel acceleration mechanism is taking place in the shock front of this remnant.

Mass outflow generated by the dynamical feedback from massive stars is currently a topic of high interest. Using a purpose-developed analysis technique, and taking full advantage of the high kinematic and angular resolution of our instrument we have detected a number of expanding superbubbles in the interacting pair of galaxies Arp 244 (NGC 4038/9) commonly known as the Antennae. We use a Fabry-P\'erot interferometer GH{\alpha}FaS to measure the profile of H{\alpha} in emission over the full extent of the object, except for the extended HI tails. The superbubbles are found centred on most of the brightest HII regions, especially in the overlap area of the two merging galaxies. We use measured sizes, expansion velocities and luminosities of the shells to estimate most of the physical parameters of the bubbles, including the kinetic energy of the expansion. In order to assess the validity of our results and approximations we perform a hydrodynamic simulation and manage to reproduce well our best measured superbubble with reasonable physical input assumptions. We also study the sources of ionization of the shells, finding that at the current, quite late stage of expansion, radiation from the remaining stars dominates, though the effect of supernova shocks can still be noted.

Three dimensional magnetohydrodynamical simulations were carried out in order to perform a new polarization study of the radio emission of the supernova remnant SN 1006. These simulations consider that the remnant expands into a turbulent interstellar medium (including both magnetic field and density perturbations). Based on the referenced-polar angle technique, a statistical study was done on observational and numerical magnetic field position-angle distributions. Our results show that a turbulent medium with an adiabatic index of 1.3 can reproduce the polarization properties of the SN 1006 remnant. This statistical study reveals itself as a useful tool for obtaining the orientation of the ambient magnetic field, previous to be swept up by the main supernova remnant shock.

We have analyzed four epochs of H$\alpha$ and [S~II] HST images of the HH~1/2 outflow (covering a time interval from 1994 to 2014) to determine proper motions and emission line fluxes of the knots of HH~2. We find that our new proper motions agree surprisingly well with the motions measured by Herbig \& Jones (1981), although there is partial evidence for a slight deceleration of the motion of the HH~2 knots from 1945 to 2014. We also measure the time-variability of the H$\alpha$ intensities and the [S~II]/H$\alpha$ line ratios, and find that knots H and A have the largest intensity variabilities (in $1994\to 2014$). Knot H (which now dominates the HH~2 emission) has strengthened substantially, while keeping an approximately constant [S~II]/H$\alpha$ ratio. Knot A has dramatically faded, and at the same time has had a substantial increase in its [S~II]/H$\alpha$ ratio. Possible interpretations of these results are discussed.

We present high-resolution, long-slit spectroscopic observations of two planetary nebulae with [WC] central stars located near the galactic bulge, M 1-32 and M 3-15. The observations were obtained with the 2.1-m telescope at the Observatorio Astron\'omico Nacional, San Pedro M\'artir. M 1-32 shows wide wings on the base of its emission lines and M 3-15 has two very faint high-velocity knots. In order to model both PNe we built a three-dimensional model consisting of a jet interacting with an equatorially concentrated slow wind, emulating the presence of a dense torus, using the Yguaz\'u hydrodynamical code. From our hydrodynamical models, we obtained position-velocity (PV) diagrams in the [N II]$\lambda$6583 line for comparison with the observations. We find that the spectral characteristics of M 1-32 and M 3-15 can be explained with the same physical model -a jet moving inside an AGB wind- using different parameters (physical conditions and position angles of the jet). In agreement with our model and observations, these objects contain a dense torus seeing pole-on and a bipolar jet escaping thorough the poles. Then we propose to classify this kind of objects as spectroscopic bipolar nebulae, although they have been classified morphologically as compact, round, or elliptical nebulae or with "close collimated lobes".

Synchrotron radiation is commonly observed in connection with shocks of different velocities, ranging from relativistic shocks associated with active galactic nuclei, gamma-ray bursts or microquasars to weakly- or non-relativistic flows as those observed in supernova remnants. Recent observations of synchrotron emission in protostellar jets are important not only because they extend the range over which the acceleration process works, but also because they allow to determine the jet and/or interstellar magnetic field structure, thus giving insights on the jet ejection and collimation mechanisms. In this paper, we compute for the first time polarized (synchrotron) and non polarized (thermal X-ray) synthetic emission maps from axisymmetrical simulations of magnetized protostellar jets. We consider models with different jet velocities and variability, as well as toroidal or helical magnetic field. Our simulations show that variable, low-density jets with velocities $\sim$ 1000 km s$^{-1}$ and $\sim$ 10 times lighter than the environment can produce internal knots with significant synchrotron emission, and thermal X-rays in the shocked region of the leading bow shock moving in a dense medium. While models with a purely toroidal magnetic field show a very large degree of polarization, models with helical magnetic field show lower values and a decrease of the degree of polarization, in agreement with observations of protostellar jets.

We use three dimensional magnetohydrodynamic (MHD) simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter $Q$, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter $Q$ maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model.

We report deep long-slit emission-line spectra, the line flux ratios, and Doppler profile shapes of various bright optical lines. The low-ionization lines (primarily [N I], [O I], [S II], and [N II]) originate in shocked knots, as reported by many previous observers. Dust-scattered lines of higher ionization are seen throughout the lobes but do not peak in the knots. Our analysis of these line profiles and the readily discernible stellar continuum shows that (1) the central star is an active symbiotic (whose spectrum resembles the central stars of highly bipolar and young PNe such as M2-9 and Hen2-437) whose compact companion shows a WC8-type spectrum, (2) extended nebular lines of [O III] and He I originate in the heavily obscured nuclear H II region, and (3) the Balmer lines observed throughout the lobes are dominated by reflected H{\alpha} emission from the symbiotic star. Comparing our line ratios with those observed historically shows that (1) the [O III]/H\b{eta} and He I/H\b{eta} ratios have been steadily rising by large amounts throughout the nebula, (2) the H{\alpha}/H\b{eta} ratio is steadily decreasing while H{\gamma}/H\b{eta} remains nearly constant, and (3) the low-ionization line ratios formed in the shocked knots have been in decline in different ways at various locations. We show that the first two of these results might be expected if the symbiotic central star has been active and if its bright H{\alpha} line has faded significantly in the past 20 years.

We present 3D hydrodynamical simulations of an isotropic fast wind interacting with a previously ejected toroidally-shaped slow wind in order to model both the observed morphology and the kinematics of the planetary nebula (PN) NGC 6302. This source, also known as the Butterfly nebula, presents one of the most complex morphologies ever observed in PNe. From our numerical simulations, we have obtained an intensity map for the H$\alpha$ emission to make a comparison with the Hubble Space Telescope (HST) observations of this object. We have also carried out a proper motion (PM) study from our numerical results, in order to compare with previous observational studies. We have found that the two interacting stellar wind model reproduces well the morphology of NGC 6302, and while the PM in the models are similar to the observations, our results suggest that an acceleration mechanism is needed to explain the Hubble-type expansion found in HST observations.

We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$\alpha$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $\Gamma$ $\sim$[1.01-1.13], while keeping fixed the stellar mass loss rate, we found a that a $\dot M_p$ range between $\sim$[3-5] $\times 10^{10}$g s$^{-1}$ give account for the observational absorption in Ly-$\alpha$ measured for the planetary system. Also, several models with anisotropic evaporation profiles for the planetary escaping atmosphere were carried out, showing that both, the escape through polar regions and through the night side yields larger absorptions than an isotropic planetary wind.

The propagation of a shock wave into a medium is expected to heat the material beyond the shock, producing noticeable effects in intensity line ratios such as [O III]/Halpha. To investigate the occurrence of shocks in planetary nebulae (PNe), we have used all narrowband [O III] and Halpha images of PNe available in the HST archive to build their [O III]/Halpha ratio maps and to search for regions where this ratio is enhanced. Regions with enhanced [O III]/Halpha emission ratio can be ascribed to two different types of morphological structures: bow-shock structures produced by fast collimated outflows and thin skins enveloping expanding nebular shells. Both collimated outflows and expanding shells are therefore confirmed to generate shocks in PNe. We also find regions with depressed values of the [O III]/Halpha ratio which are found mostly around density bounded PNe, where the local contribution of [N II] emission into the F656N Halpha filter cannot be neglected.

We present a "thin shell" model of the interaction of a biconical outflow and a streaming environment (aligned with the direction of the flow), as well as numerical (axisymmetric) simulations of such an interaction. A similar situation, although in a more complex setup, takes place at the head of the cometary structure of Mira. Thus, for most of the numerical simulations we explore parameters consistent with the observed bipolar outflow from Mira B. For these parameters, the interaction is non-radiative, so that a rather broad jet/streaming environment interaction region is formed. In spite of this, a reasonable agreement between the thin-shell analytic model and the numerical simulations is obtained.

The morphology, dynamics and thermal X-ray emission of the superbubble N70 is studied by means of 3D hydrodynamical simulations, carried out with the {\sc{yguaz\'u-a}} code. We have considered different scenarios: the superbubble being the product of a single supernova remnant, of the stellar winds from an OB association, or the result of the joint action of stellar winds and a supernova event. Our results show that, in spite that all scenarios produce bubbles with the observed physical size, only those where the bubble is driven by stellar winds and a SN event are successful to explain the general morphology, dynamics and the X-ray luminosity of N70. Our models predict temperatures in excess of $10^8 \mathrm{K}$ at the interior of the superbubble, however the density is too low and the emission in thermal X-ray above $2 \mathrm{keV}$ is too faint to be detected.

In this work we employ an MHD numerical code to reproduce the morphology observed for SN 1006 in radio synchrotron and thermal X-ray emission. We introduce a density discontinuity, in the form of a flat cloud parallel to the Galactic Plane, in order to explain the NW filament observed in optical wavelengths and in thermal X-rays. We compare our models with observations. We also perform a test that contrasts the radio emitting bright limbs of the SNR against the central region, finding additional support to our results. Our main conclusion is that the most probable direction of the ambient magnetic field is on average perpendicular to the Galactic Plane.

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star $\eta$ Car. The two bipolar nebulae commonly known as the Homunculus and the little Homunculus were formed from the interaction of these eruptive events with the underlying stellar wind. As in previous work (Gonzalez et al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed {\it decreased} for a short period of time. This fact suggests that the little Homunculus is formed when the eruption ends, from the impact of the post-outburst $\eta$ Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The little Homunculus (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resembles the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the $\eta$ Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of $\eta$ Car.

We present 3D hydrodynamical simulations of a precessing jet propagating inside a supernova remnant (SNR) shell, particularly applied to the W50-SS433 system in a search for the origin of its peculiar elongated morphology. Several runs were carried out with different values for the mass loss rate of the jet, the initial radius of the SNR, and the opening angle of the precession cone. We found that our models successfully reproduce the scale and morphology of W50 when the opening angle of the jets is set to 10$\degr$ or if this angle linearly varies with time. For these models, more realistic runs were made considering that the remnant is expanding into an interstellar medium (ISM) with an exponential density profile (as HI observations suggest). Taking into account all these ingredients, the large scale morphology of the W50-SS 433 system, including the asymmetry between the lobes (formed by the jet-SNR interaction), is well reproduced.

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.

We present hydrodynamical calculations carried out with the 3D yguazu-a code of a precessing jet model, which interacts with a plane parallel wind. This scenario describes an extragalactic jet, in which the jet source is in motion with respect to the surrounding intra-cluster medium. From the numerical results, synthetic emission maps and spectra in X-ray band were obtained. We compare these predictions with observations of the radio jets emanating from the radio-galaxy 4C 26.42 (in the Abell 1795 galaxy cluster). We find that the general morphology of the radio jets can be described by a point-symmetric precessing jet system interacting with a plane parallel wind (i.e., the intra-cluster medium flowing past the galaxy). We also find that our synthetic X-ray emission maps reproduce the observed large scale structures (with sizes of the order of tens of kpc).

We here present two-dimensional, time-dependent radiatively cooling hydrodynamical simulations of the large and little Homunculus nebulae around $\eta$ Carinae. We employ an alternative scenario to previous interacting stellar wind models which is supported by both theoretical and observational evidence, where a non-spherical outburst wind (with a latitudinal velocity dependence that matches the observations of the large Homunculus), which is expelled for 20 years, interacts with a pre-eruptive slow wind also with a toroidal density distribution, but with a much smaller equator-to-polar density contrast than that assumed in previous models. A second eruptive wind with spherical shape is ejected about 50 years after the first outburst, and causes the development of the little internal nebula. We find that, as a result of an appropriate combination of the parameters that control the degree of asymmetry of the interacting winds, we are able to produce not only the structure and kinematics of both Homunculus, but also the high-velocity equatorial ejecta. These arise from the impact between the non-spherical outburst and the pre-outburst winds in the equatorial plane.

We present 3D, gasdynamic simulations of jet/cloud collisions, with the purpose of modelling the HH 270/110 system. From the models, we obtain predictions of H$\alpha$ and H$_2$ 1-0 s(1) emission line maps, which qualitatively reproduce some of the main features of the corresponding observations of HH 110. We find that the model that better reproduces the observed structures corresponds to a jet that was deflected at the surface of the cloud $\sim 1000$ yr ago, but is now boring a tunnel directly into the cloud. This model removes the apparent contradiction between the jet/cloud collision model and the lack of detection of molecular emission in the crossing region of the HH 270 and HH 110 axes.