Engineering topics

Blondin, John M.

Publications and source records attributed to Blondin, John M..

18 records

NASA NTRS2017

The Three-Dimensional Expansion of the Ejecta from Tycho's Supernova Remnant

We present the first 3D measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12 yr baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 tufts of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line-of-sight velocity, we use two different methods: a nonequilibrium ionization model fit to the strong Si and S lines in the 1.22.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the redshift or blueshift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km/s, with a mean of 4430 km/s. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km/s. Some SN Ia explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and we discuss our findings in light of various explosion models, favoring those delayed-detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

ISM: individual objects (Tychoâ€™s SNR) â€“ ISM: s↗

NASA NTRS2016

An X-Ray and Radio Study of the Varying Expansion Velocities in Tycho's Supernova Remnant

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using Chandra and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 years in the X-rays and 30 years in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. Consistent with earlier measurements, we find a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is approximately 23″ toward the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

dust↗

NASA NTRS1995

Wind dynamics in SMC X-1. 1: Hydrodynamic simulation

We present a three-dimensional hydrodynamic simulation of the disrupted stellar wind in the high-mass X-ray binary system SMC X-1. The three dominant processes that determine the geometry of the wind in high X-ray luminosity systems such as SMC X-1 are the X-ray suppression of the stellar wind from the X-ray irradiated face of the primary star, the focusing of the radiatively driven wind in the X-ray shadow by the effects of stellar rotation, and the rapid X-ray heating of gas in the vicinity of the X-ray source, including the X-ray illuminated surface of the primary star. The resulting distribution of circumstellar gas provides a successful explanation for the asymmetric, extended eclipse transitions and the intensity of the deep eclipse X-ray emission in SMC X-1, as well as a possible explanation for the X-ray dips seen near superior conjunction of the X-ray source in Cyg X-1.

Blondin, John M.↗

NASA NTRS1995

Wind dynamics in SMC X-1. 2: Ginga and ROSAT observations

The X-ray phenomena of the binary system SMC X-1/Sk 160, observed with the Ginga and ROSAT X-ray observatories, are compared with computed phenomena derived from a three-dimensional hydrodynamical model of the stellar wind perturbed by X-ray heating and ionization which is described in the accompanying paper (Blondin & Woo 1995). In the model the B0 I primary star has a line-driven stellar wind in the region of the X-ray shadow and a thermal wind in the region heated by X-rays. We find general agreement between the observed and predicted X-ray spectrum throughout the binary orbit cycle, including the extended, variable, and asymmetric eclipse transitions and the period of deep eclipse.

Woo, Jonathan W.↗

NASA NTRS1994

The shadow wind in high-mass X-ray binaries

We present hydrodynamic simulations of the most X-ray luminous high-mass X-ray binary systems, such as Cen X-3 and SMC X-1. These models illustrate the presence of both a normal radiatively driven wind confined to the X-ray shadow of the primary star -- a shadow wind -- and a thermally driven wind excited by the X-ray heating of the primary's stellar surface -- an X-ray-excited wind. The X-ray flux in these systems is sufficiently intense that any circumstellar gas exposed to the X-ray source will be highly photoionized. These extreme ionization conditions prevent the formation of a normal radiatively driven wind from the irradiated surface of the primary, but such a wind can still form on the shadowed side of the primary. Orbital rotation can then bring this shadow wind into the line of sight toward the X-ray source, enhancing the column density of the wind seen near eclipse egress. Furthermore, such a high X-ray flux can also excite a thermal wind from the irradiated surface of the primary. Again, orbital rotation tends to deflect the wind, this time leading to an enhanced column density near eclipse ingress.

Blondin, John M.↗

NASA NTRS1994

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗

NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

NASA NTRS1993

The piecewise-parabolic method in curvilinear coordinates

We derive interpolation formulae for a third-order finite difference method in curvilinear, orthogonal coordinate systems. These formulae serve as a supplement to Colella and Woodward's PPM scheme for problems where the coordinate origin is included in the computational domain. Numerical examples of the improved accuracy of the advection scheme near coordinate singularities are shown.

Blondin, John M.↗

NASA NTRS1992

Dynamics of Kepler's supernova remnant

Observations of Kepler's SNR have revealed a strong interaction with the ambient medium, far in excess of that expected at a distance of about 600 pc away from the Galactic plane where Kepler's SNR is located. This has been interpreted as a result of the interaction of supernova ejecta with the dense circumstellar medium (CSM). Based on the bow-shock model of Bandiera (1985), we study the dynamics of this interaction. The CSM distribution consists of an undisturbed stellar wind of a moving supernova progenitor and a dense shell formed in its interaction with a tenuous interstellar medium. Supernova ejecta drive a blast wave through the stellar wind which splits into the transmitted and reflected shocks upon hitting this bow-shock shell. We identify the transmitted shock with the nonradiative, Balmer-dominated shocks found recently in Kepler's SNR. The transmitted shock most probably penetrated the shell in the vicinity of the stagnation point.

Borkowski, Kazimierz J.↗

NASA NTRS1992

The evolution of cocoons surrounding light, extragalactic jets

If the mass density of supersonic, collimated material is less than that of the surrounding medium, a so-called light jet will be enveloped by a cocoon of overpressured shocked gas. Hydrodynamical simulations are used to understand the evolution of the cocoon. The cocoon's evolution is also compared to a simple analytic theory. To reconcile the theory with the simulations, the growth of the jet head must be taken into account. The overpressured cocoon stage exists for a relatively short astronomical time, after which only the region of the cocoon near the jet head remains overpressured. The spatial distribution of the optical emission often observed in distant extragalactic jet systems can be explained with this improved understanding of cocoon evolution.

Cioffi, Denis F.↗

NASA NTRS1992

Hydrodynamic instabilities in supernova remnants - Self-similar driven waves

An initial study aimed at elucidating the multidimensional aspects of the hydrodynamic instabilities in supernova remnants is presented. Self-similar solutions are found to exist for the interaction of a steep power-law density profile expanding into a relatively flat stationary power-law density profile. Consideration of the pressure and entropy profiles in the shocked 1D flows shows that the flows are subject to convective instability, by a local criterion. The growth rate for the instability becomes very large near the contact discontinuity between the two shocked regions. A linear analysis of the complete self-similar solutions shows that the solutions are unstable above a critical wavenumber and that the growth rate is greatest at the position of the contact discontinuity. The X-ray image of the remnant of SN 1572 (Tycho) shows emission from clumps of supernova ejecta, which is good evidence for instabilities in this remnant.

Chevalier, Roger A.↗

NASA NTRS1992

Colliding winds from early-type stars in binary systems

The dynamics of the wind and shock structure formed by the wind collision in early-type binary systems is examined by means of a 2D hydrodynamics code, which self-consistently accounts for radiative cooling, and represents a significant improvement over previous attempts to model these systems. The X-ray luminosity and spectra of the shock-heated region, accounting for wind attenuation and the influence of different abundances on the resultant level and spectra of X-ray emission are calculated. A variety of dynamical instabilities that are found to dominate the intershock region is examined. These instabilities are found to be particularly important when postshock material is able to cool. These instabilities disrupt the postshock flow and add a time variability of order 10 percent to the X-ray luminosity. The X-ray spectrum of these systems is found to vary with the nuclear abundances of winds. These theoretical models are used to study several massive binary systems, in particular V444 Cyg and HD 193793.

Stevens, Ian R.↗

NASA NTRS1991

Enhanced winds and tidal streams in massive X-ray binaries

The tidal effects created by the presence of a compact companion are expected to induce a stream of enhanced wind from the early-type primary star in massive X-ray binary systems. In this paper, two-dimensional gasdynamical simulations of such streams are presented. It is found that the wind enhancement is a sensitive function of the binary separation, and develops into a tidal stream as the primary approaches its critical surface. For typical system parameters, the Coriolis force deflects the stream sufficiently that it does not impact directly on the compact companion but passes behind it. The density in the stream can reach values of 20-30 times the ambient wind density, leading to strong attenuation of the X-ray flux that passes through the tidal stream, providing a possible explanation of the enhanced absorption events seen at later phases in the X-ray observations of massive X-ray binary systems such as Vela X-1. In contrast to the time-variable accretion wake, the tidal stream is relatively stationary, producing absorption features that should remain fixed from orbit to orbit. For systems with a strong tidal stream, the large asymmetry in the accreting wind results in the accretion of angular momentum of constant sign, as opposed to systems without streams, where the sign of the accreted angular momentum can change.

Blondin, John M.↗

NASA NTRS1990

The structure and evolution of radiatively cooling jets

The two-dimensional simulations presently used to characterize the structure and evolution of radiatively cooling supersonic jets reveal that cooling jet morphologies resemble those of adiabatic outflows, but with the fundamental difference that a dense, cold shell will condense out of the shocked gas at the head of the jet when the cooling distance behind either of the two principal shocks is smaller than the jet radius. For very high cooling rates, the material that accumulates at the head of the jet forms an extended plug of cold gas resembling the 'nose cone' observed in numerical simulations of strongly magnetized adiabatic jets. An investigation is made of the dependence of jet properties on the density ratio between the beam and the ambient medium, as well as on the strength of radiative cooling.

Blondin, John M.↗

NASA NTRS1990

Hydrodynamic simulations of stellar wind disruption by a compact X-ray source

This paper presents two-dimensional numerical simulations of the gas flow in the orbital plane of a massive X-ray binary system, in which the mass accretion is fueled by a radiation-driven wind from an early-type companion star. These simulations are used to examine the role of the compact object (either a neutron star or a black hole) in disturbing the radiatively accelerating wind of the OB companion, with an emphasis on understanding the origin of the observed soft X-ray photoelectric absorption seen at late orbital phases in these systems. On the basis of these simulations, it is suggested that the phase-dependent photoelectric absorption seen in several of these systems can be explained by dense filaments of compressend gas formed in the nonsteady accreation bow shock and wake of the compact object.

Blondin, John M.↗

NASA NTRS1989

The growth of density perturbations in radiative shocks

The paper presents a new investigation of the growth of thermal instabilities behind radiative shocks. The analytic and numerical results disagree with McCray, Stein, and Kafatos, who were attempting to explain the filamentary structure commonly seen in supernova remnants, and proposed that the growth of density perturbations would lead to gross condensations in the postshock cooling region. Hydrodynamical simulations are generated which corroborate arguments in favor of no growth in the long-wavelength limit. These simulations also agree with a renewed analytic approach in the short-wavelength limit, and show that the region of rapid growth will remain inconspicuous in the overall density rise toward the back of the shock. It is noted that these calculations are in accord with some recent observations of supernova remnants where the observed filaments do not seem to have been produced by any local thermal instability process.

Blondin, John M.↗

NASA NTRS1989

Herbig-Haro objects as the heads of radiative jets

The interpretation of certain HH objects as the heads of nonadiabatic supersonic jets is examined using two-dimensional numerical simulations. It is found that radiative jets develop a dense shell between the jet shock and the leading bow shock when the cooling distance behind either one of these shocks is smaller than the jet radius. It is proposed that the radiatively cooling shell may account for the variable emission pattern from objects like HH 1. Also, it is suggested that HH objects with measured space velocities that exceed the spectroscopically inferred shock velocities could correspond to heavy jets in which the bow shock is effectively adiabatic. Low-excitation objects in which these velocities are comparable may represent light jets where the jet shock is nonradiative.

Blondin, John M.↗

NASA NTRS1986

Is the 1.5-ms pulsar a young neutron star?

It is proposed that PSR1937 + 214 is a young neutron star spun up by accretion from a high-mass companion in a close binary system. The supercritical mass transfer rates expected in such a binary system should allow the neutron star to be spun up in the comparatively short time of about 10,000 yr. The accretion process will also power thermomagnetic effects that could remove the strong magnetic field of a young pulsar from the crust of the star in a similarly short timescale. Such a high-mass binary system is expected to disrupt when the companion explodes in a supernova. Thus a spin-up model in a high-mass system can explain the lack of a companion, low magnetic field, and high spin rate of PSR1937 + 214.

Blondin, John M.↗