The interstellar bubble RCW 120 seen around a type O runaway star is driven by the stellar wind and the ionising radiation emitted by the star. The boundary between the stellar wind and interstellar medium (ISM) is associated with the arc-shaped mid-infrared dust emission around the star within the HII region.

Orion BN/KL is an example of a poorly understood phenomena in star forming regions involving the close encounter of young stellar objects. The explosive structure, the great variety of molecules observed, the energy involved in the event and the mass of the region suggest a contribution in the {chemical diversity} of the local interstellar medium. Nevertheless, the frequency and duration of other events like this have not been determined. In this paper, we explore a recent analytic model that takes into account the interaction of a clump with its molecular environment. We show that the widespread kinematic ages of the Orion fingers -- 500 to 4000 years -- is a consequence of the interaction of the explosion debris with the surrounding medium. This model explains satisfactorily the age discrepancy of the Orion fingers, and infers the initial conditions together with the lifetime of the explosion. Moreover, our model can explain why some CO streamers do not have a H$_2$ finger associated.

The interaction of a high velocity clump of gas has been described by the plasmon model, which considers balance between ram pressure and the internal stratified structure of the decelerated clump. In this paper we propose an analytical model to describe the mass loss of such a clump due the interaction with the environment, describing its influence on the plasmon dynamics. We carry out comparisons between an analytic model and axisymetric gasdynamic simulations of plasmon evolution. From our simulations we were able to find the values of the friction constants $\alpha$ and $\lambda$. Comparing with the complete analytic model from which we can infer the position and the mass loss of the clump as function of the clump's density and the environment ratio.

We studied the kinematics of the Galactic bubble RCW 120 in the [SII]\lambda\lambda$6717,6731\AA lines. We measured a LSR radial velocity ranging from ~-74 to ~-6km/s. We found evidence of expansion only in the northeast region of the nebula (from 20 to 30 km/s). We found a high electron density around 4000 cm^{-3} in the south-west region and we also found two arches-like structure indicating a density gradient. We present 3D numerical simulations of RCW 120 using Walicxe-3D code in order to explore optical shell dynamics and its morphology. Our numerical results predict an average numerical electron density of the ambient medium (in the southern region of the object) is between 3000 to 5000 cm^{-3} in agreement with our values obtained from the observations. From our models, we do not expect X-ray emission coming from the external shell, due to the low expansion velocity value.}

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.

We present results from models of galactic winds driven by energy injected from nuclear (at the galactic center) and non-nuclear starbursts. The total energy of the starburst is provided by very massive young stellar clusters,which can push the galactic interstellar medium and produce an important outflow. Such outflow can be a well, or partially mixed wind, or a highly metallic wind. We have performed adiabatic 3D N-Body/Smooth Particle Hydrodynamics simulations of galactic winds using the GADGET-2 code. The numerical models cover a wide range of parameters, varying the galaxy concentration index, gas fraction of the galactic disk, and radial distance of the starburst. We show that an off-center starburst in dwarf galaxies is the most effective mechanism to produce a significant loss of metals (material from the starburst itself). At the same time a non-nuclear starburst produce a high efficiency of metal loss, in spite of having a moderate to low mass loss rate.

We propose a new Planetary Nebula Luminosity Function (PNLF) that includes two populations in the distribution. Our PNLF is a direct extension of the canonical function proposed by Jacoby et al. (1987), in order to avoid problems related with the histogram construction, it is cast in terms of cumulative functions. We are interested in recovering the shape of the faint part of the PNLF in a consistent manner, for galaxies with and without a dip in their PN luminosity functions. The parameters for the two mode PNLF are obtained with a genetic algorithm, which obtains a best fit to the PNLF varying all of the parameters simultaneously in a broad parameter space. We explore a sample of 9 galaxies with various Hubble types and construct their PNLF. All of the irregular galaxies, except one, are found to be consistent with a two-mode population, while the situation is less clear for ellipticals and spirals.For the case of NGC\, 6822, we show that the two-mode PNLF is consistent with previous studies of the star formation history within that galaxy. Our results support two episodes of star formation, in which the latter is significantly stronger.

An analysis of the OIII and SIII temperatures measurements compiled by Perez-Montero et al. of emission line objects consisting of HII galaxies, giant extragalactic HII regions, Galactic HII regions and HII regions from the Magellanic Clouds, reveals that the OIII temperatures are higher than the corresponding values from SIII in most objects with gas metallicities in excess of 0.2 solar. We explore the possibility of temperature inhomogeneities. We also explored metallicity inhomogeneities by combining two models of widely different metallicity. We calculate models that consider a non-Mawell-Boltzmann distributions for the electron energies (kappa parametrization). We also consider shock heating within the photoionized nebula. Varying the various input parameters in the pure photoionization case does not reproduce the observed nebular temperatures and neither does having local temperature inhomogeneities. We find that (1) metallicity inhomogeneities of the nebular gas, (2) shock waves of velocities < 60 km/s propagating in a photoionized plasma, and (3) an electron energy distribution given by a kappa-distribution, are successful in reproducing the observed excess in the [OIII] temperatures. Shock models, however, would require proper 3D hydrodynamical simulations to become a fully developed alternative while models with metallicity inhomogeneities appear to fail in metal poor nebulae, since they result in T_O++(rec)> T_OIII.

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.

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 analytic and numerical models of the `cluster wind' resulting from the multiple interactions of the winds ejected by the stars of a dense cluster of massive stars. We consider the case in which the distribution of stars (i.e., the number of stars per unit volume) within the cluster is spherically symmetric, has a power-law radial dependence, and drops discontinuously to zero at the outer radius of the cluster. We carry out comparisons between an analytic model (in which the stars are considered in terms of a spatially continuous injection of mass and energy) and 3D gasdynamic simulations (in which we include 100 stars with identical winds, located in 3D space by statistically sampling the stellar distribution function). From the analytic model, we find that for stellar distributions with steep enough radial dependencies the cluster wind flow develops a very high central density and a non-zero central velocity, and for steeper dependencies it becomes fully supersonic throughout the volume of the cluster (these properties are partially reproduced by the 3D numerical simulations). Therefore, the wind solutions obtained for stratified clusters can differ dramatically from the case of a homogeneous stellar distribution (which produces a cluster wind with zero central velocity, and a fully subsonic flow within the cluster radius). Finally, from our numerical simulations we compute predictions of X-ray emission maps and luminosities, which can be directly compared with observations of cluster wind flows.