The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous single-cloud studies. In reality, the flow behind the shock is also likely to be turbulent, with non-uniform density, pressure and velocity structure created as the shock sweeps over inhomogenities upstream of the cloud. To address these twin issues we use a sub-grid compressible k-epsilon turbulence model to estimate the properties of the turbulence generated in shock-cloud interactions and the resulting increase in the transport coefficients that the turbulence brings. A detailed comparison with the output from an inviscid hydrodynamical code puts these new results into context. We find that cloud destruction in inviscid and k-epsilon models occurs at roughly the same speed when the post-shock flow is smooth and when the density contrast between the cloud and inter-cloud medium is less than 100. However, there are increasing and significant differences as this contrast increases. Clouds subjected to strong ``buffeting'' by a highly turbulent post-shock environment are destroyed significantly quicker. Additional calculations with an inviscid code where the post-shock flow is given random, grid-scale, motions confirms the more rapid destruction of the cloud. Our results clearly show that turbulence plays an important role in shock-cloud interactions, and that environmental turbulence adds a new dimension to the parameter space which has hitherto been studied (abridged).

We have examined a stream-source model for the production of the cometary tails observed in the Helix Nebula NGC 7293 in which a transonic or moderately supersonic stream of ionized gas overruns a source of ionized gas. Hydrodynamic calculations reveal velocity structures which are in good agreement with the observational data on tail velocities and are consistent with observations of the nebular structure. The results also are indicative of a stellar atmosphere origin for the cometary globules. Tail remnants persist for timescales long enough for their identification with faint striations visible in the nebula gas to be plausible.

Because of its key role in feedback in star formation and galaxy formation, we examine the nature of the interaction of a flow with discrete sources of mass injection. We show the results of two-dimensional numerical simulations in which we explore a range of configurations for the mass sources and study the effects of their proximity on the downstream flow. The mass sources act effectively as a single source of mass injection if they are so close together that the ratio of their combined mass injection rate is comparable to or exceeds the mass flux of the incident flow into the volume that they occupy. The simulations are relevant to many diffuse sources, such as planetary nebulae and starburst superwinds, in which a global flow interacts with material evaporating or being ablated from the surface of globules of cool, dense gas.

We investigate the Mach number dependence of the radiative overstability of shocks. We find that the stability of radiative shocks increases with decreasing Mach number, with the result that M=2 shocks require cooling exponents < -1.2 to be overstable. The frequency of oscillation of the fundamental mode has a strong Mach number dependence. We find that feedback between the cooling region and the cold dense layer (CDL) further downstream is a function of Mach number, with stronger feedback and oscillation of the boundary between the CDL and the cooling region occuring at lower Mach numbers. This feedback can be quantified in terms of the reflection coefficient of sound waves. An interesting finding is that the stability properties of low Mach number shocks can be dramatically altered if the shocked gas is able to cool to temperatures less than the pre-shock value, and it is probable that low Mach number astrophysical shocks will be overstable in a variety of situations. The results of this work are relevant to astrophysical shocks with low Mach numbers, such as supernova remnants (SNRs) immersed in a hot interstellar medium (e.g., within a starburst region), and shocks in molecular clouds, where time-dependent chemistry can lead to overstability. (Abridged)

One aspect of supernova remnant evolution that is relatively unstudied is the influence of an AGN environment. A high density ambient medium and a nearby powerful continuum source will assist the cooling of shocked ejecta and swept-up gas. Motion of the surrounding medium relative to the remnant will also affect the remnant morphology. In an extension to previous work we have performed 2D hydrodynamical calculations of SNR evolution in an AGN environment, and have determined the evolutionary behaviour of cold gas in the remnant. The cold gas will contribute to the observed broad line emission in AGNs, and we present preliminary theoretical line profiles from our calculations. A more detailed comparison with observations will be performed in future work. The SNR-AGN interaction may be useful as a diagnostic of the AGN wind.

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.

We study the structure of shocks in clumpy media, using a multifluid formalism. As expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous subshock appears in the structure. This has significant implications for the survival of dense clouds in regions overrun by shocks in a wide range of astrophysical circumstances, from planetary nebulae to the nuclei of starburst galaxies.

We show that a cooled region of shocked supernova ejecta forms in a type II supernova-QSO wind interaction, and has a density, an ionization parameter, and a column density compatible with those inferred for the high ionization component of the broad emission line regions in QSOs. The calculations are based on the assumption that the ejecta flow is described initially by a similarity solution investigated by Chevalier (1982) and Nadyozhin (1985) and is spherically symmetric. Heating and cooling appropriate for gas irradiated by a nearby powerful continuum source is included in our model, together with reasonable assumptions for the properties of the QSO wind. The model results are also in agreement with observational correlations and imply reasonable supernova rates.

We present similarity solutions for adiabatic bubbles that are blown by winds having time independent mechanical luminosities and that are each mass-loaded by the hydrodynamic ablation of distributed clumps. The mass loading is `switched-on' at a specified radius (with free-expansion of the wind interior to this point) and injects mass at a rate per unit volume proportional to M^delta r^lambda where delta = 4/3 (1) if the flow is subsonic (supersonic) with respect to the clumps. In the limit of negligible mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. The presence of mass loading heats the flow, which leads to a reduction in the Mach number of the supersonic freely-expanding flow, and weaker jump conditions across the inner shock. In solutions with large mass loading, it is possible for the wind to connect directly to the contact discontinuity without first passing through an inner shock. For a solution that gives the mass of swept-up ambient gas to be less than the sum of the masses of the wind and ablated material, lambda < -2. Maximum possible values for the ratio of ablated mass to wind mass occur when mass loading starts very close to the bubble center and when the flow is supersonic with respect to the clumps over the entire bubble radius. The maximum temperature in the bubble often occurs near the onset of mass loading, and in some cases can be many times greater than the post-inner-shock temperature. Our solutions are relevant to eg stellar wind-blown bubbles, galactic winds, etc. This work complements Pittard et al (2001) where it was assumed that clumps were evaporated through conductive energy transport.

We present similarity solutions for mass-loaded adiabatic bubbles that are blown by winds having time dependent mechanical luminosities. We consider mass-loading through conductive evaporation of clumps. In the limit of little mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. We find that the Mach number in the shocked mass-loaded wind shows a radial dependence that varies qualitatively from solution to solution. In some cases it is everywhere less than unity in the frame of the clumps being evaporated, while in others it is everywhere greater than unity. In some solutions the mass-loaded shocked wind undergoes one or two sonic transitions in the clump frame. Maximum possible values of the ratio of evaporated mass to stellar wind mass are found as a consequence of the evaporation rates dependence on temperature and the lowering of the temperature by mass-loading. Mass-loading tends to reduce the emissivity in the interior of the bubble relative to its limb, whilst simultaneously increasing the central temperature relative to the limb temperature.

We present the jump conditions for ionization fronts with oblique magnetic fields. The standard nomenclature of R- and D-type fronts can still be applied, but in the case of oblique magnetization there are fronts of each type about each of the fast- and slow-mode speeds. As an ionization front slows, it will drive first a fast- and then a slow-mode shock into the surrounding medium. Even for rather weak upstream magnetic fields, the effect of magnetization on ionization front evolution can be important. [Includes numerical MHD models and an application to observations of S106.]

We have modelled ultracompact HII regions (UCHIIR) in terms of steady subsonic ionized flows in a clumpy medium. Mass loss from neutral clumps allows the regions to be long-lived. We examine the form of global flows for different dependences of the volume mass injection rate, $\dot{q}$, on radius and Mach number, and describe the solutions in detail. We find that three observed UCHIIR morphologies are reproduced with these models. Mach number independent flows that include a radial variation can give centre- brightened core--halo morphologies. Mach number dependent flows reproduce naturally the uniform UCHIIR morphology. In a hybrid model, including subsonic and supersonic flows, we allow a supersonic wind to shock in the ionized region. The ionized subsonic gas has a high density and so dominates the emission. The shell produced has a velocity structure very different from that of fully supersonic models. Several morphologies of spherical UCHIIR can be understood in terms of these various models; however, kinematic data are crucial as a discriminant between them.

Cometary ultracompact HII regions have been modelled as the interaction of the hypersonic wind from a moving star with the molecular cloud which surrounds the star. We here show that a similar morphology can ensue even if the star is stationary with respect to the cloud material. We assume that the HII region is within a stellar wind bubble which is strongly mass loaded: the cometary shape results from a gradient in the distribution of mass loading sources. This model circumvents problems associated with the necessarily high spatial velocities of stars in the moving star models.

We exploit a generic instability in the integration of steady, sonic, near-isothermal flows to find the complete transition diagram for recombination fronts (for a model system of equations). The instability requires the integration of the flow equations for speeds between the isothermal and adiabatic sound speeds to be performed with particular care. As a result of this, the previous work of Newman & Axford on the structure of recombination fronts neglected an important class of solution, that of transonic fronts; our method is readily extensible to a more complete treatment of the ionization structure. Future papers will apply these results in models of the structure of ultracompact HII regions.

We have used high resolution longslit spectroscopy to investigate the ionized gas in the active galaxy IRAS 04210+0400 and its association with the radio structure. We suggest that two of the ionized components are associated with the central double radio source and observe that the relative positions of these components vary for different emission lines. Both results are consistent with the radio components representing the working surfaces of a pair of jets emerging from the centre of the galaxy. In this scenario, the optical emission in the centre arises behind the bowshocks produced by the jets in the interstellar medium. The emission lines are detected and show a dramatic (~ 900 km/s) spread in velocity at the position of the radio lobe hotspots. We suggest a model which explains this phenomenon as the result of a jet head emerging through the boundary between the interstellar and intergalactic medium. A similar scenario has previously been suggested as a model to explain wide angle tail radio sources (WAT's). Based on this model, we simulate the longslit spectra of these regions and compare the results with the observations.

We propose that a significant fraction of the ultracompact HII regions found in massive star-forming clouds are the result of the interaction of the wind and ionizing radiation from a young massive star with the clumpy molecular cloud gas in its neighbourhood. Distributed mass loading in the flow allows the compact nebulae to be long-lived. In this paper, we discuss a particularly simple case, in which the flow in the HII region is everywhere supersonic. The line profiles predicted for this model are highly characteristic, for the case of uniform mass loading. We discuss briefly other observational diagnostics of these models.