We test the physical model of the relativistic jets in the galactic X-ray binary SS 433 proposed in our previous paper using additional observations from the Chandra High Energy Transmission Grating Spectrometer. These observations sample two new orbital/precessional phase combinations. In the observation near orbital phase zero, the H- and He-like Fe lines from both receding and approaching jets are comparably strong and unocculted while the He-like Si line of the receding jet is significantly weaker than that of the approaching jet. This condition may imply the cooler parts of the receding jet are eclipsed by the companion. The X-ray spectrum from this observation has broader emission lines than obtained in Paper I that may arise from the divergence of a conical outflow or from Doppler shift variations during the observation. Using recent optical results, along with the length of the unobscured portion of the receding jet assuming adiabatic cooling, we calculate the radius of the companion to be 9.6+/-1.0 R_sun, about one third of the Roche lobe radius. For a main sequence star, this corresponds to a companion mass of 35+/-7 M_sun, giving a primary source mass of 20+/-5 M_sun. If our model is correct, this calculation indicates the compact object is a black hole, and accretion occurs through a wind process. In a subsequent paper, we will examine the validity of the adiabatic cooling model of the jets and test the mode of line broadening.

We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to wavelength is greater than 0.02, portions of the IF temporarily separate from the molecular cloud surface, locally decreasing the ablation pressure. This causes the appearance of a large, warm HI region and triggers nonlinear dynamics of the IF. The local difference of the ablation pressure and acceleration enhances the appearance and growth of a multimode perturbation. The stabilization usually seen at the IF in the linear regimes does not work due to the mismatch of the modes of the perturbations at the cloud surface and in density in HII region above the cloud surface. Molecular pillars are observed in the late stages of the large amplitude perturbation case. The velocity gradient in the pillars is in reasonably good agreement with that observed in the Eagle Nebula. The initial perturbation is imposed in three different ways: in density, in incident photon number flux, and in the surface shape. All cases show both stabilization for a small initial perturbation and large growth of the second harmonic by increasing amplitude of the initial perturbation above a critical value.

We investigate hydrodynamic instability of accelerating ionization fronts using two dimensional hydrodynamic simulations that include detailed energy deposition and release due to the absorption of UV radiation, recombination of hydrogen, radiative molecular cooling, and magnetic pressure. We consider linear perturbation growth and find that the stabilization mechanism associated with non-accelerated fronts remains a significant factor even when acceleration is present. Conversely, if recombination in the ionized region is turned off,Rayleigh-Taylor (RT) instability becomes effective, and the classical RT growth rate recovered.

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.