Stellar bow shocks result from relative motions between stars and their environment. The interaction of the stellar wind and radiation with gas and dust in the interstellar medium produces curved arcs of emission at optical, infrared, and radio wavelengths. We recently proposed a new two-dimensional classification scheme for the shape of such bow shocks, which we here apply to three very different observational datasets: mid-infrared arcs around hot OB stars; far-infrared arcs around luminous cool stars; and H alpha emission-line arcs around proplyds and other young stars in the Orion Nebula. For OB stars, the average shape is consistent with simple thin-shell models for the interaction of a spherical wind with a parallel stream, but the diversity of observed shapes is many times larger than such models predict. We propose that this may be caused by time-dependent oscillations in the bow shocks, due to either instabilities or wind variability. Cool star bow shocks have markedly more closed wings than hot star bow shocks, which may be due to the dust emission arising in the shocked stellar wind instead of the shocked interstellar medium. The Orion Nebula arcs, on the other hand, have both significantly more open wings and significantly flatter apexes than the hot star bow shocks. We test several possible explanations for this difference (divergent ambient stream, low Mach number, observational biases, and influence of collimated jets), but the evidence for each is inconclusive.

Dust waves and bow waves result from the action of a star's radiation pressure on a stream of dusty plasma that flows past it. They are an alternative mechanism to hydrodynamic bow shocks for explaining the curved arcs of infrared emission seen around some stars. When gas and grains are perfectly coupled, for a broad class of stellar parameters, wind-supported bow shocks predominate when the ambient density is below 100 per cubic cm. At higher densities radiation-supported bows can form, tending to be optically thin bow waves around B stars, or optically thick bow shocks around early O stars. For OB stars with particularly weak stellar winds, radiation-supported bows become more prevalent.

Stellar bow shocks, bow waves, and dust waves all result from the action of a star's wind and radiation pressure on a stream of dusty plasma that flows past it. The dust in these bows emits prominently at mid-infrared wavelengths in the range 8 to 60 micron. We propose a novel diagnostic method, the tau-eta diagram, for analyzing these bows, which is based on comparing the fractions of stellar radiative energy and stellar radiative momentum that is trapped by the bow shell. This diagram allows the discrimination of wind-supported bow shocks, radiation-supported bow waves, and dust waves in which grains decouple from the gas. For the wind-supported bow shocks, it allows the stellar wind mass-loss rate to be determined. We critically compare our method with a previous method that has been proposed for determining wind mass-loss rates from bow shock observations. This comparison points to ways in which both methods can be improved and suggests a downward revision by a factor of two with respect to previously reported mass-loss rates. From a sample of 23 mid-infrared bow-shaped sources, we identify at least 4 strong candidates for radiation-supported bow waves, which need to be confirmed by more detailed studies, but no strong candidates for dust waves.

Dust waves are a result of gas-grain decoupling in a stream of dusty plasma that flows past a luminous star. The radiation field is sufficiently strong to overcome the collisional coupling between grains and gas at a "rip-point", where the ratio of radiation pressure to gas pressure exceeds a critical value of roughly 1000. When the rip point occurs outside the hydrodynamic bow shock, a separate dust wave may form, decoupled from the gas shell, which can either be drag-confined or inertia-confined, depending on the stream density and relative velocity. In the drag-confined case, there is a minimum stream velocity of roughly 60 km/s that allows a steady-state stagnant drift solution for the dust wave apex. For lower relative velocities, the dust dynamics close to the axis exhibit a limit cycle behavior (rip and snap back) between two different radii. Strong coupling of charged grains to the plasma's magnetic field can modify these effects, but for a quasi-parallel field orientation the results are qualitatively similar to the non-magnetic case. For a quasi-perpendicular field, on the other hand, the formation of a decoupled dust wave is strongly suppressed.

X-ray emitting diffuse nebulae around hot stars are observed to have soft-band temperatures in the narrow range [1-3]$\times10^{6}$ K, independent of the stellar wind parameters and the evolutionary stage of the central star. We discuss the origin of this X-ray temperature for planetary nebulae (PNe), Wolf-Rayet nebulae (WR) and interstellar wind bubbles around hot young stars in our Galaxy and the Magellanic Clouds. We calculate the differential emission measure (DEM) distributions as a function of temperature from previously published simulations and combine these with the X-ray emission coefficient for the 0.3-2.0 keV band to estimate the X-ray temperatures. We find that all simulated nebulae have DEM distributions with steep negative slopes, which is due to turbulent mixing at the interface between the hot shocked stellar wind and the warm photoionised gas. Sharply peaked emission coefficients act as temperature filters and emphasize the contribution of gas with temperatures close to the peak position, which coincides with the observed X-ray temperatures for the chemical abundance sets we consider. Higher metallicity nebulae have lower temperature and higher luminosity X-ray emission. We show that the second temperature component found from spectral fitting to X-ray observations of WR nebulae is due to a significant contribution from the hot shocked stellar wind, while the lower temperature principal component is dominated by nebular gas. We suggest that turbulent mixing layers are the origin of the soft X-ray emission in the majority of diffuse nebulae.

We present a study of the X-ray emission from numerical simulations of hot bubbles in planetary nebulae (PNe). High-resolution, two-dimensional, radiation-hydrodynamical simulations of the formation and evolution of hot bubbles in PNe, with and without thermal conduction, are used to calculate the X-ray emission and study its time-dependence and relationship to the changing stellar parameters. Instabilities in the wind-wind interaction zone produce clumps and filaments in the swept-up shell of nebular material. Turbulent mixing and thermal conduction at the corrugated interface can produce quantities of intermediate temperature and density gas between the hot, shocked wind bubble and the swept-up photoionized nebular material, which can emit in soft, diffuse X-rays. We use the CHIANTI software to compute synthetic spectra for the models and calculate their luminosities. We find that models both with conduction and those without can produce the X-ray temperatures and luminosities that are in the ranges reported in observations, although the models including thermal conduction are an order of magnitude more luminous than those without. Our results show that at early times the diffuse X-ray emission should be dominated by the contribution from the hot, shocked stellar wind, whereas at later times the nebular gas will dominate the spectrum. We analyse the effect of sampling on the resultant spectra and conclude that a minimum of 200 counts is required to reliably reproduce the spectral shape. Likewise, heavily smoothed surface-brightness profiles obtained from low-count detections of PNe do not provide a reliable description of the spatial distribution of the X-ray emitting gas.

In order to study the nature, origin, and impact of turbulent velocity fluctuations in the ionized gas of the Orion Nebula, we apply a variety of statistical techniques to observed velocity cubes. The cubes are derived from high resolving power ($R \approx 40,000$) longslit spectroscopy of optical emission lines that span a range of ionizations. From Velocity Channel Analysis (VCA), we find that the slope of the velocity power spectrum is consistent with predictions of Kolmogorov theory between scales of 8 and 22 arcsec (0.02 to 0.05 pc). The outer scale, which is the dominant scale of density fluctuations in the nebula, approximately coincides with the autocorrelation length of the velocity fluctuations that we determine from the second order velocity structure function. We propose that this is the principal driving scale of the turbulence, which originates in the autocorrelation length of dense cores in the Orion molecular filament. By combining analysis of the non-thermal line widths with the systematic trends of velocity centroid versus ionization, we find that the global champagne flow and smaller scale turbulence each contribute in equal measure to the total velocity dispersion, with respective root-mean-square widths of 4-5 km/s. The turbulence is subsonic and can account for only one half of the derived variance in ionized density, with the remaining variance provided by density gradients in photoevaporation flows from globules and filaments. Intercomparison with results from simulations implies that the ionized gas is confined to a thick shell and does not fill the interior of the nebula.

We present deep XMM-Newton EPIC observations of the Wolf-Rayet (WR) bubble NGC6888 around the star WR136. The complete X-ray mapping of the nebula confirms the distribution of the hot gas in three maxima spatially associated with the caps and northwest blowout hinted at by previous Chandra observations. The global X-ray emission is well described by a two-temperature optically thin plasma model $T_1$=1.4$\times$10$^{6}$ K, $T_{2}$=8.2$\times$10$^{6}$ K) with a luminosity of $L_{\mathrm{X}}$=7.8$\times$10$^{33}$ erg s$^{-1}$ in the 0.3--1.5 keV energy range. The rms electron density of the X-ray-emitting gas is estimated to be $n_\mathrm{e}$=0.4 cm$^{-3}$. The high-quality observations presented here reveal spectral variations within different regions in NGC6888, which allowed us for the first time to detect temperature and/or nitrogen abundance inhomogeneities in the hot gas inside a WR nebula. One possible explanation for such spectral variations is that the mixing of material from the outer nebula into the hot bubble is less efficient around the caps than in other nebular regions.

We discuss our most recent findings on the diffuse X-ray emission from Wolf-Rayet (WR) nebulae. The best-quality X-ray observations of these objects are those performed by XMM-Newton and Chandra towards S308, NGC2359, and NGC6888. Even though these three WR nebulae might have different formation scenarios, they all share similar characteristics: i) the main plasma temperatures of the X-ray-emitting gas is found to be $T$=[1-2]$\times$10$^{6}$ K, ii) the diffuse X-ray emission is confined inside the [O III] shell, and iii) their X-ray luminosities and electron densities in the 0.3-2.0~keV energy range are $L_\mathrm{X}\approx$10$^{33}$-10$^{34}$~erg~s$^{-1}$ and $n_\mathrm{e}\approx$0.1-1~cm$^{-3}$, respectively. These properties and the nebular-like abundances of the hot gas suggest mixing and/or thermal conduction is taking an important role reducing the temperature of the hot bubble.

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We analyze the expansion of hydrogen-poor knots and filaments in the born-again planetary nebulae A30 and A78 based on Hubble Space Telescope (HST) images obtained almost 20 yr apart. The proper motion of these features generally increases with distance to the central star, but the fractional expansion decreases, i.e., the expansion is not homologous. As a result, there is not a unique expansion age, which is estimated to be 610-950 yr for A30 and 600-1140 yr for A78. The knots and filaments have experienced complex dynamical processes: the current fast stellar wind is mass loaded by the material ablated from the inner knots; the ablated material is then swept up until it shocks the inner edges of the outer, hydrogen-rich nebula. The angular expansion of the outer filaments shows a clear dependence on position angle, indicating that the interaction of the stellar wind with the innermost knots channels the wind along preferred directions. The apparent angular expansion of the innermost knots seems to be dominated by the rocket effect of evaporating gas and by the propagation of the ionization front inside them. Radiation-hydrodynamical simulations show that a single ejection of material followed by a rapid onset of the stellar wind and ionizing flux can reproduce the variety of clumps and filaments at different distances from the central star found in A30 and A78.

We investigate the scale dependence of fluctuations inside a realistic model of an evolving turbulent HII region and to what extent these may be studied observationally. We find that the multiple scales of energy injection from champagne flows and the photoionization of clumps and filaments leads to a flatter spectrum of fluctuations than would be expected from top-down turbulence driven at the largest scales. The traditional structure function approach to the observational study of velocity fluctuations is shown to be incapable of reliably determining the velocity power spectrum of our simulation. We find that a more promising approach is the Velocity Channel Analysis technique of Lazarian & Pogosyan (2000), which, despite being intrinsically limited by thermal broadening, can successfully recover the logarithmic slope of the velocity power spectrum to a precision of +-0.1 from high resolution optical emission line spectroscopy.

We carry out high resolution two-dimensional radiation-hydrodynamic numerical simulations to study the formation and evolution of hot bubbles inside planetary nebulae (PNe). We take into account the evolution of the stellar parameters, wind velocity and mass-loss rate from the final thermal pulses during the asymptotic giant branch (AGB) through to the post-AGB stage for a range of initial stellar masses. The instabilities that form at the interface between the hot bubble and the swept-up AGB wind shell lead to hydrodynamical interactions, photoevaporation flows and opacity variations. We explore the effects of hydrodynamical mixing combined with thermal conduction at this interface on the dynamics, photoionization, and emissivity of our models. We find that even models without thermal conduction mix significant amounts of mass into the hot bubble. When thermal conduction is not included, hot gas can leak through the gaps between clumps and filaments in the broken swept-up AGB shell and this depressurises the bubble. The inclusion of thermal conduction evaporates and heats material from the clumpy shell, which expands to seal the gaps, preventing a loss in bubble pressure. The dynamics of bubbles without conduction is dominated by the thermal pressure of the thick photoionized shell, while for bubbles with thermal conduction it is dominated by the hot, shocked wind.

The Wolf-Rayet (WR) bubble S 308 around the WR star HD 50896 is one of the only two WR bubbles known to possess X-ray emission. We present XMM-Newton observations of three fields of this WR bubble that, in conjunction with an existing observation of its Northwest quadrant, map most of the nebula. The X-ray emission from S 308 displays a limb-brightened morphology, with a central cavity ~22 arcmin in size and a shell thickness of ~8 arcmin. This X-ray shell is confined by the optical shell of ionized material. The spectrum is dominated by the He-like triplets of NIV at 0.43 keV and OVII at 0.57 keV, and declines towards high energies, with a faint tail up to 1 keV. This spectrum can be described by a two-temperature optically thin plasma emission model (T1 ~ 1.1x10^6 K, T2 ~ 13x10^6 K), with a total X-ray luminosity ~2x10^33 erg/s at the assumed distance of 1.5 kpc.

We study the evolution of the interstellar and circumstellar media around massive stars (M > 40M_{\odot}) from the main sequence through to the Wolf-Rayet stage by means of radiationhydrodynamic simulations. We use publicly available stellar evolution models to investigate the different possible structures that can form in the stellar wind bubbles around Wolf-Rayet stars. We find significant differences between models with and without stellar rotation, and between models from different authors. More specifically, we find that the main ingredients in the formation of structures in the Wolf-Rayet wind bubbles are the duration of the Red Supergiant (or Luminous Blue Variable) phase, the amount of mass lost, and the wind velocity during this phase, in agreement with previous authors. Thermal conduction is also included in our models. We find that main-sequence bubbles with thermal conduction are slightly smaller, due to extra cooling which reduces the pressure in the hot, shocked bubble, but that thermal conduction does not appear to significantly influence the formation of structures in post-main-sequence bubbles. Finally, we study the predicted X-ray emission from the models and compare our results with observations of the Wolf-Rayet bubbles S\,308, NGC\,6888, and RCW\,58. We find that bubbles composed primarily of clumps have reduced X-ray luminosity and very soft spectra, while bubbles with shells correspond more closely to observations.

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

We present the first three-dimensional radiation-magnetohydrodynamic simulations of the photoionisation of a dense, magnetised molecular globule by an external source of ultraviolet radiation. We find that, for the case of a strong ionising field, significant deviations from the non-magnetic evolution are seen when the initial magnetic field threading the globule has an associated magnetic pressure that is greater than one hundred times the gas pressure. In such a strong-field case, the photoevaporating globule will adopt a flattened or "curled up" shape, depending on the initial field orientation, and magnetic confinement of the ionised photoevaporation flow can lead to recombination and subsequent fragmentation during advanced stages of the globule evolution. We find suggestive evidence that such magnetic effects may be important in the formation of bright, bar-like emission features in H II regions. We include simple but realistic fits to heating and cooling rates in the neutral and molecular gas in the vicinity of a high-mass star cluster and show that the frequently used isothermal approximation can lead to an overestimate of the importance of gravitational instability in the radiatively imploded globule. For globules within 2 parsecs of a high-mass star cluster, we find that heating by stellar x rays prevents the molecular gas from cooling below 50 K.

The three dimensional structure of the brightest part of the Orion Nebula is assessed in the light of published and new data. We find that the widely accepted model of a concave blister of ionized material needs to be altered in the southwest direction from the Trapezium, where we find that the Orion-S feature is a separate cloud of very optically thick molecules within the body of ionized gas, which is probably the location of the multiple embedded sources that produce the outflows that define the Orion-S star formation region. Evidence for this cloud comes from the presence of H2CO lines in absorption in the radio continuum and discrepancies in the extinction derived from radio-optical and optical only emission. We present an equilibrium Cloudy model of the Orion-S cloud, which successfully reproduces many observed properties of this feature. We also report the discovery of an open-sided shell of [O III] surrounding the Trapezium stars, revealed through emission line ratio images and the onset of radiation shadows beyond some proplyds. We show that the observed properties of the shell are consistent with it being a stationary structure, produced by shock interactions between the ambient nebular gas and the high-velocity wind from theta^1 Ori C. We examine the implications of the recently published evidence for a large blueshifted velocity of theta^1 Ori C with respect to the Orion Molecular Cloud, which could mean that this star has only recently begun to photoionize the Orion Nebula. We show that current observations of the Nebula do not rule out such a possibility, so long as the ionization front has propagated into a pre-existing low-density region. In addition, a young age for the Nebula would help explain the presence of nearby proplyds with a short mass-loss timescale to photoablation.

The expansion of ionization fronts in uniform and spherically symmetric power-law density distributions is a well-studied topic. However, in many situations, such as a star formed at the edge of a molecular cloud core, an offset power-law density distribution would be more appropriate. In this paper a few of the main issues of the formation and expansion of HII regions in such media are outlined and results are presented for the particular cases where the underlying power laws are 1/r^2 and 1/r^3. A simple criterion is developed for determining whether the initial photoionized region will be unbounded, which depends on the power-law exponent and the ratio of the equivalent Stroemgren radius produced by the star in a uniform medium to the stellar offset distance. In the expansion stage, the ionized volumes will eventually become unbounded unless pressure balance with the external medium is reached before the ionization front velocity becomes supersonic with respect to the ionized gas.

Most stars will experience episodes of substantial mass loss at some point in their lives. For very massive stars, mass loss dominates their evolution, although the mass loss rates are not known exactly, particularly once the star has left the main sequence. Direct observations of the stellar winds of massive stars can give information on the current mass-loss rates, while studies of the ring nebulae and HI shells that surround many Wolf-Rayet (WR) and luminous blue variable (LBV) stars provide information on the previous mass-loss history. The evolution of the most massive stars, (M > 25 solar masses), essentially follows the sequence O star to LBV or red supergiant (RSG) to WR star to supernova. For stars of mass less than 25 solar masses there is no final WR stage. During the main sequence and WR stages, the mass loss takes the form of highly supersonic stellar winds, which blow bubbles in the interstellar and circumstellar medium. In this way, the mechanical luminosity of the stellar wind is converted into kinetic energy of the swept-up ambient material, which is important for the dynamics of the interstellar medium. In this review article, analytic and numerical models are used to describe the hydrodynamics and energetics of wind-blown bubbles. A brief review of observations of bubbles is given, and the degree to which theory is supported by observations is discussed.

We present numerical radiation-hydrodynamic simulations of cometary HII regions for a number of champagne flow and bowshock models. For the champagne flow models we study smooth density distributions with both steep and shallow gradients. We also consider cases where the ionizing star has a strong stellar wind, and cases in which the star additionally has a proper motion within the ambient density gradient. We present simulated emission-measure maps and long-slit spectra of our results. Our numerical models are not tailored to any particular object but comparison with observations from the literature shows that, in particular, the models combining density gradients and stellar winds are able to account for both the morphology and radial velocity behavior of several observed cometary HII regions, such as the well-studied object G29.96-0.02.

We present numerical radiation-hydrodynamic simulations of the evolution of HII regions formed in an inhomogeneous medium resulting from turbulence simulations. We find that the filamentary structure of the underlying density distribution produces a highly irregular shape for the ionized region, in which the ionization front escapes to large distances in some directions within 80,000 years. In other directions, on the other hand, neutral gas in the form of dense globules persists within 1 parsec of the central star for the full duration of our simulation (400,000 years). Divergent photoablation flows from these globules maintain a root-mean-squared velocity in the ionized gas that is close to the ionized sound speed. Simulated images in optical emission lines show morphologies that are in strikingly detailed agreement with those observed in real HII regions.

We present hydrodynamical simulations of the photoevaporation of a cloud with large-scale density gradients, giving rise to an ionized, photoevaporation flow. The flow is found to be approximately steady during the large part of its evolution, during which it can resemble a "champagne flow" or a "globule flow" depending on the curvature of the ionization front. The distance from source to ionization front and the front curvature uniquely determine the structure of the flow, with the curvature depending on the steepness of the lateral density gradient in the neutral cloud. We compare these simulations with both new and existing observations of the Orion nebula and find that a model with a mildly convex ionization front can reproduce the profiles of emission measure, electron density, and mean line velocity for a variety of emitting ions on scales of 10^{17} to 10^{18} cm. The principal failure of our model is that we cannot explain the large observed widths of the [O I] 6300 Angstrom line that forms at the ionization front.

We present a method for including steady-state gas flows in the plasma physics code Cloudy, which was previously restricted to modeling static configurations. The numerical algorithms are described in detail, together with an example application to plane-parallel ionization-bounded HII regions. As well as providing the foundation for future applications to more complex flows, we find the following specific results regarding the effect of advection upon ionization fronts in HII regions: 1. Significant direct effects of advection on the global emission properties occur only when the ionization parameter is lower than is typical for HII regions. 2. The overheating of partially ionized gas in the front is not large, even for supersonic (R-type) fronts. 3. The most significant morphological signature of advective fronts is an electron density spike that occurs at the ionization front. Observational evidence for such a spike is found in images of the Orion bar. 4. Plane-parallel, weak-D fronts are found to show at best a shallow correlation between mean velocity and ionization potential for optical emission lines even when the flow velocity closely approaches the ionized sound speed.

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.

Many low-mass stars in the Orion nebula are associated with very compact (~= 1 arcsec) emission knots, known variously as proplyds, PIGs or LV knots. Some of these knots are teardrop-shaped, with ``tails'' pointing away from the massive star Theta 1 Ori C, which is the principal exciting star of the nebula. We discuss models of such knots, which invoke the interaction of the fast stellar wind from Theta 1 Ori C with a transonic photoevaporated flow from the surface of an accretion disk around a young low-mass star. We review previous analytic work and compare the results of the model with the observed brightnesses, morphologies and emission line profiles of the knots, as well as presenting new results from numerical hydrodynamical simulations.