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.

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.

Massive young star clusters contain dozens or hundreds of massive stars that inject mechanical energy in the form of winds and supernova explosions, producing an outflow which expands into their surrounding medium, shocking it and forming structures called superbubbles. The regions of shocked material can have temperatures in excess of 10$^6$ K, and emit mainly in thermal X-rays (soft and hard). This X-ray emission is strongly affected by the action of thermal conduction, as well as by the metallicity of the material injected by the massive stars. We present three-dimensional numerical simulations exploring these two effects, metallicity of the stellar winds and supernova explosions, as well as thermal conduction.

We propose an asymmetrical jet ejection mechanism in order to model the mirror symmetry observed in the lobe distribution of some protoplanetary nebulae (pPNe), such as the pPN CRL 618. 3D hydrodynamical simulations of a precessing jet launched from an orbiting source were carried out including an alternation in the ejections of the two outflow lobes, depending on which side of the precessing accretion disk is hit by the accretion column from a Roche lobe-filling binary companion. Both synthetic optical emission maps and position-velocity (PV) diagrams were obtained from the numerical results with the purpose of carrying out a direct comparison with observations. Depending on the observer's point of view, multipolar morphologies are obtained which exhibit a mirror symmetry at large distances from the central source. The obtained lobe sizes and their spatial distribution are in good agreement with the observed morphology of the pPN CRL 618. We also obtain that the kinematic ages of the fingers are similar to those obtained in the observations.

We present 3D, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit towards the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates $\dot{M} \sim 10^{-8}$ M$_\odot$ yr$^{-1}$. Predicted Br$\gamma$ emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br$\gamma$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$\gamma$ luminosity ($\sim 8 \times 10^{30}$ erg s$^{-1}$) is reproduced by a model with a $v_w=50$ km s$^{-1}$ wind velocity and a $10^{-7}$ M$_\odot$ yr$^{-1}$ mass loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br$\gamma$ luminosity if the emission comes from the inner, free wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br$\gamma$ emission is more compact after periastron passage.

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.