Supersonic outflows from objects as varied as stellar jets, massive stars and novae often exhibit multiple shock waves that overlap one another. When the intersection angle between two shock waves exceeds a critical value, the system reconfigures its geometry to create a normal shock known as a Mach stem where the shocks meet. Mach stems are important for interpreting emission-line images of shocked gas because a normal shock produces higher postshock temperatures and therefore a higher-excitation spectrum than an oblique one does. In this paper we summarize the results of a series of numerical simulations and laboratory experiments designed to quantify how Mach stems behave in supersonic plasmas that are the norm in astrophysical flows. The experiments test analytical predictions for critical angles where Mach stems should form, and quantify how Mach stems grow and decay as intersection angles between the incident shock and a surface change. While small Mach stems are destroyed by surface irregularities and subcritical angles, larger ones persist in these situations, and can regrow if the intersection angle changes to become more favorable. The experimental and numerical results show that although Mach stems occur only over a limited range of intersection angles and size scales, within these ranges they are relatively robust, and hence are a viable explanation for variable bright knots observed in HST images at the intersections of some bow shocks in stellar jets.

We investigate the role discrete clumps embedded in an astrophysical jet play on the jet's morphology and line emission characteristics. By varying clumps' size, density, position, and velocity, we cover a range of parameter space motivated by observations of objects such as the Herbig Haro object HH~34. We here extend the results presented in Yirak et al. 2009, including how analysis of individual observations may lead to spurious sinusoidal variation whose parameters vary widely over time, owing chiefly to interacts between clumps. The goodness of the fits, while poor in all simulations, are best when clump-clump collisions are minimal. Our results indicate that a large velocity dispersion leads to a clump-clump collision-dominated flow which disrupts the jet beam. Finally, we present synthetic emission images of H-$\alpha$ and [SII] and note an excess of [SII] emission along the jet length as compared to observations. This suggests that observed beams undergo earlier processing, if they are present at all.

Numeric convergence studies demonstrate that the evolution of an adiabatic clump is well-captured by roughly 100 cells per clump radius. The presence of radiative cooling, however, imposes limits on the problem due to the removal of thermal energy. Numerical studies which include radiative cooling typically adopt the 100--200 cells per clump radius resolution. In this paper we present the results of a convergence study for radiatively cooling clumps undertaken over a broad range of resolutions, from 12 to 1,536 cells per clump radius, employing adaptive mesh refinement (AMR) in a 2D axisymmetric geometry ("2.5D"). We also provide a fully 3D simulation, at 192 cells per clump radius, which supports our 2.5D results. We find no appreciable self-convergence at ~100 cells per clump radius as small-scale differences owing to increasingly resolving the "cooling length" have global effects. We therefore conclude that self-convergence is an insufficient criterion to apply on its own when addressing the question of sufficient resolution for radiatively cooled shocked clump simulations. We suggest the adoption of alternate criteria to support a statement of sufficient resolution, such as the demonstration of adequate resolution of the cooling layers behind shocks. We discuss an associated refinement criteria for AMR codes.

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of plasma traveling at hypersonic velocities. Structure within jet beams is often attributed to periodic or ``pulsed'' variations of conditions at the jet source. Simulations based on this scenario result in knots extending across the jet diameter. Observations and recent high energy density laboratory experiments shed new light on structures below this scale and indicate they may be important for understanding the fundamentals of jet dynamics. In this paper we offer an alternative to ``pulsed'' models of protostellar jets. Using direct numerical simulations we explore the possibility that jets are chains of sub-radial clumps propagating through a moving inter-clump medium. Our models explore an idealization of this scenario by injecting small ($r<r_{jet}$), dense ($\rho>\rho_{jet}$) spheres embedded in an otherwise smooth inter-clump jet flow. The spheres are initialized with velocities differing from the jet velocity by $\sim15$%. We find the consequences of shifting from homogeneous to heterogeneous flows are significant as clumps interact with each other and with the inter-clump medium in a variety of ways. Structures which mimic what is expected from pulsed-jet models can form, as can previously unseen ``sub-radial'' behaviors including backward facing bow shocks and off-axis working surfaces. While these small-scale structures have not been seen before in simulation studies, they are found in high resolution jet observations. We discuss implications of our simulations for the interpretation of protostellar jets with regard to characterization of knots by a ``lifetime'' or ``velocity history'' approach as well as linking observed structures with central engines which produce the jets.

We have performed 2D hydrodynamic simulations of a pulsed astrophysical jet propagating through a medium that is populated with spherical inhomogeneities, or "clumps," which are smaller than the jet radius. The clumps are seen to affect the jet in several ways, such as impeding jet propagation and deflecting the jet off-axis. While there has been some debate as to the prevalence of these types of condensations in the ISM or in molecular clouds, the exploration of this region of parameter space nonetheless both shows the potential for these clumps to disrupt astrophysical jets and yields results which recover aspects of recent observations of Herbig-Haro objects. We find that the propagation of the jet and the vorticity induced in the clump/ambient medium correlate well with a "dynamic filling function" $f_d$ across all the simulations.

We report follow-up observations of the Sloan Digital Sky Survey Cataclysmic Variable SDSS J210014.12+004446.0 (hereafter SDSS J2100). We obtained photometry and spectroscopy in both outburst and quiescent states, providing the first quiescent spectrum of this source. In both states, non-sinusoidal photometric modulations are apparent, suggestive of superhumps, placing SDSS J2100 in the SU UMa subclass of dwarf novae. However, the periods during outburst and quiescence differ significantly, being 2.099 plus or minus 0.002 hr and 1.96 plus or minus 0.02 hr respectively. Our phase-resolved spectroscopy during outburst yielded an estimate of about 2 hr for the orbital period, consistent with the photometry. The presence of the shorter period modulation at quiescence is unusual, but not unique. Another atypical feature is the relative weakness of the Balmer emission lines in quiescence. Overall, we find a close similarity between SDSS J2100 and the well-studied superhump cataclysmic Variable V503 Cygni. By analogy, we suggest that the quiescent modulation is due to a tilted accretion disk -- producing negative superhumps -- and the modulation in outburst is due to positive superhumps from the precession of an elliptical disk.

We report the results of a 45 ks Chandra observation of the cataclysmic variable V426 Ophiuchus. The high resolution spectrum from the high-energy transmission grating spectrometer is most consistent with a cooling flow model, placing V426 Oph among the group of CVs including U Gem and EX Hya. An uninterrupted lightcurve was also constructed, in which we detect a significant 4.2 hr modulation together with its first harmonic at 2.1 hrs. Reanalysis of archival Ginga, and ROSAT X-ray lightcurves also reveals modulations at periods consistent with 4.2 and/or 2.1 hrs. Furthermore, optical photometry in V, simultaneous with the Chandra observation, indicates a modulation anti-correlated with the X-ray, and later more extensive R band photometry finds a signal at ~2.1 hrs. The earlier reported X-ray periods at ~0.5 and 1 hrs appear to be only transient and quasi-periodic in nature. In contrast, the 4.2 hr period or its harmonic are stable and persistent in X-ray/optical data from 1988 to 2003. This periodicity is clearly distinct from the 6.85 hr orbit, and could be due to the spin of the white dwarf. If this is the case, V426 Oph would be the first long period intermediate polar with a ratio P_spin/P_orb of 0.6. However, this interpretation requires unreasonable values of magnetic field strength and mass accretion rate.