Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Blast waves”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗
NASA NTRS2009

The Supercritical Pile GRB Model: The Prompt to Afterglow Evolution

The "Supercritical Pile" is a very economical GRB model that provides for the efficient conversion of the energy stored in the protons of a Relativistic Blast Wave (RBW) into radiation and at the same time produces - in the prompt GRB phase, even in the absence of any particle acceleration - a spectral peak at energy approx. 1 MeV. We extend this model to include the evolution of the RBW Lorentz factor Gamma and thus follow its spectral and temporal features into the early GRB afterglow stage. One of the novel features of the present treatment is the inclusion of the feedback of the GRB produced radiation on the evolution of Gamma with radius. This feedback and the presence of kinematic and dynamic thresholds in the model can be the sources of rich time evolution which we have began to explore. In particular. one can this may obtain afterglow light curves with steep decays followed by the more conventional flatter afterglow slopes, while at the same time preserving the desirable features of the model, i.e. the well defined relativistic electron source and radiative processes that produce the proper peak in the (nu)F(sub nu), spectra. In this note we present the results of a specific set of parameters of this model with emphasis on the multiwavelength prompt emission and transition to the early afterglow.

Mastichiadis, A.↗
NASA NTRS2010

Evolution of Interplanetary Shocks and their CME Drivers

Shock-driving coronal mass ejections (CMEs) constitute the most energetic phenomena in the heliosphere. The shocks can be identified in a number of ways based on remote-sensing and in situ observations. Type II radio bursts are the earliest indicators of shocks that accelerate electrons to energies up to -10 keV. Solar energetic particle (SEP) events are always accompanied by long wavelength type II bursts indicating that the same shock accelerates ions and electrons. A recent investigation involving a large number of interplanetary (IP) shocks revealed that about 35% of them do not produce type II bursts (radio quiet, RQ) or SEPs. Comparison of the RQ shocks with the radio loud (RL) ones revealed some interesting results such as: (1) the lack of evidence for blast waves,(2) energetic particle enhancement in the shock front in -20% of RQ shocks, and (3) determination of the difference between the RQ and RL shocks in terms of the different kinematic properties of the associated CMEs. On the other hand the shock properties measured at I AU are not too different for the RQ and RL cases. This can be attributed to the interaction with the IP medium, which seems to erase the difference. Implications of this evolution for the geoeffectiveness is also discussed.

Gopalswamy, Natchimuthuk↗
NASA NTRS2010

Influences of the Driver and Ambient Medium Characteristics on the Formation of Shocks in the Solar Atmosphere

Traveling interplanetary (IP) shocks were discovered in the early 1960s, but their solar origin has been controversial. Early research focused on solar flares as the source of the shocks, but when coronal mass ejections (CMEs) were discovered, it became clear that fast CMEs clearly can drive the shocks. Type II radio bursts are excellent signatures of shocks near the Sun. The close correspondence between type II radio bursts and solar energetic particles (SEPs) makes it clear that the same shock accelerates ions and electrons. A recent investigation involving a large number of IP shocks revealed that about 35% of IP shocks do not produce type II bursts or SEPs. Comparing these radio quiet (RQ) shocks with the radio loud (RL) ones revealed some interesting results: (1) there is no evidence for blast waves, in that all IP shocks can be attributed to CMEs, (2) a small fraction (20%) of RQ shocks is associated with ion enhancements at the shocks when they move past the observing spacecraft, (3) the primary difference between the RQ and RL shocks can be traced to the different kinematic properties of the associated CMEs and the variation of the characteristic speeds of the ambient medium, and (4) the shock properties measured at 1 AU are not too different for the RQ and RL cases due to the interaction of the shock driver with the IP medium that seems to erase the difference.

Nat, Gopalswamy↗
NASA NTRS2010

Error Representation in Time For Compressible Flow Calculations

Time plays an essential role in most real world fluid mechanics problems, e.g. turbulence, combustion, acoustic noise, moving geometries, blast waves, etc. Time dependent calculations now dominate the computational landscape at the various NASA Research Centers but the accuracy of these computations is often not well understood. In this presentation, we investigate error representation (and error control) for time-periodic problems as a prelude to the investigation of feasibility of error control for stationary statistics and space-time averages. o These statistics and averages (e.g. time-averaged lift and drag forces) are often the output quantities sought by engineers. o For systems such as the Navier-Stokes equations, pointwise error estimates deteriorate rapidly which increasing Reynolds number while statistics and averages may remain well behaved.

Barth, Timothy J.↗
NASA NTRS2010

Shock-Driving CMEs Near the Sun, in the Interplanetary Medium, and Near Earth

The excellent correspondence between type II radio bursts and solar energetic particles (SEPs) made it clear that the same shock accelerates ions and electrons. A recent investigation involving a large number of I shocks revealed that about 35% of IP shocks do not produce type II bursts (radio quiet) or SEPs. Comparing the RQ shocks with the radio loud (RL) ones revealed some interesting results, which will be summarized in this poster. (1) There is no evidence for blast waves. (2) Even a small fraction (20%) of RQ shocks is associated with ion enhancements at the shock when the shock passes the spacecraft. (3) The primary difference between the RQ and RL shocks can be traced to the different kinematic properties of the associated CMEs, although the shock properties measured at 1 AU are not too different for the RQ and RL cases. This can be attributed to the interaction with the IP medium, which seems to erase the difference. More details can be found in Astrophysical Journal 710, 1111, 2010 (http://adsabs.harvard.edu/absZ2009arXivO9l2.4719G).

Gopalswamy, N.↗
NASA NTRS2010

A Bulk Comptonization Model for the Prompt GRM Emission

The "Supercritical Pile" is a very economical GRB model that provides for the efficient conversion of the energy stored in the protons of a Relativistic Blast Wave (RBW) into radiation and at the same time produces - in the prompt GRB phase, even in the absence of any particle acceleration - a spectral peak at energy approximately 1 MeV. We extend this model to include the evolution of the RBW Lorentz factor F and thus follow its spectral and temporal features into the early GRB afterglow stage. One of the novel features of the present treatment is the inclusion of the feedback of the GRB produced radiation on the evolution of Gamma with radius. This feedback and the presence of kinematic and dynamic thresholds in the model are sources of potentially very rich time evolution which we have began to explore. In particular, one can this way obtain afterglow light curves with steep decays followed by the more conventional flatter afterglow slopes, while at the same time preserving the desirable features of the model, i.e. the well defined relativistic electron source and radiative processes that produce the proper peak in the nu F(sub nu) spectra. In this note we present the results of a specific set of parameters of this model with emphasis on the multiwavelength prompt emission and transition to the early afterglow.

Kazanas, Demos↗
NASA NTRS2011

On Interplanetary Shocks Driven by Coronal Mass Ejections

Traveling interplanetary (IP) shocks were first detected in the early 1960s, but their solar origin has been controversial. Early research focused on solar flares as the source of the shocks, but when CMEs were discovered, it became clear that fast CMEs are the shock drivers. Type radio II bursts are excellent signatures of shocks near the Sun (Type II radio bursts were known long before the detection of shocks and CMEs). The excellent correspondence between type II bursts and solar energetic particle (SEP) events made it clear that the same shock accelerates ions and electrons. Shocks near the Sun are also seen occasionally in white-light coronagraphic images. In the solar wind, shocks are observed as discontinuities in plasma parameters such as density and speed. Energetic storm particle events and sudden commencement of geomagnetic storm are also indicators of shocks arriving at Earth. After an overview on these shock signatures, I will summarize the results of a recent investigation of a large number of IP shocks. The study revealed that about 35% of IP shocks do not produce type II bursts (radio quiet - RQ) or SEPs. Comparing the RQ shocks with the radio loud (RL) ones revealed some interesting results: (1) There is no evidence for blast wave shocks. (2) A small fraction (20%) of RQ shocks is associated with ion enhancements at the shock when the shock passes the spacecraft. (3) The primary difference between the RQ and RL shocks can be traced to the different kinematic properties of the associated CMEs. On the other hand the shock properties measured at 1 AU are not too different for the RQ and RL cases. This can be attributed to the interaction with the IP medium, which seems to erase the difference between the shocks.

Gopalswarmy, Nat↗
NASA NTRS2007

A Physical Relationship Between Electron-Proton Temperature Equilibration and Mach Number in Fast Collisionless Shocks

The analysis of Balmer-dominated optical spectra from nonradiative (adiabatic) SNRs has shown that the ratio of the electron to proton temperature at the blast wave is close to unity at v(s) approximate or smaller than 400 km s(-1) but declines sharply down to the minimum value of m(e) /m(p) dictated by the jump conditions at shock speeds exceeding 2000 km s(-1). We propose a physical model for the heating of electrons and ions in non-cosmic-ray-dominated, strong shocks (v(s) 400 km s(-1)) wherein the electrons are heated by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of electron heating over a wide range of shock speeds, producing a relationship (T(e)/T(p))sub 0 proportional to v(-2/s) (proportional to M(-2)) that is fully consistent with the observations.

SNR(ONRADIATIVE SUPERNOVA REMNANTS)↗
NASA NTRS2010

Interplanetary Shocks Lacking Type 2 Radio Bursts

We report on the radio-emission characteristics of 222 interplanetary (IP) shocks detected by spacecraft at Sun-Earth L1 during solar cycle 23 (1996 to 2006, inclusive). A surprisingly large fraction of the IP shocks (approximately 34%) was radio quiet (RQ; i.e., the shocks lacked type II radio bursts). We examined the properties of coronal mass ejections (CMEs) and soft X-ray flares associated with such RQ shocks and compared them with those of the radio-loud (RL) shocks. The CMEs associated with the RQ shocks were generally slow (average speed approximately 535 km/s) and only approximately 40% of the CMEs were halos. The corresponding numbers for CMEs associated with RL shocks were 1237 km/s and 72%, respectively. Thus, the CME kinetic energy seems to be the deciding factor in the radio-emission properties of shocks. The lower kinetic energy of CMEs associated with RQ shocks is also suggested by the lower peak soft X-ray flux of the associated flares (C3.4 versus M4.7 for RL shocks). CMEs associated with RQ CMEs were generally accelerating within the coronagraph field of view (average acceleration approximately +6.8 m/s (exp 2)), while those associated with RL shocks were decelerating (average acceleration approximately 3.5 m/s (exp 2)). This suggests that many of the RQ shocks formed at large distances from the Sun, typically beyond 10 Rs, consistent with the absence of metric and decameter-hectometric (DH) type II radio bursts. A small fraction of RL shocks had type II radio emission solely in the kilometric (km) wavelength domain. Interestingly, the kinematics of the CMEs associated with the km type II bursts is similar to those of RQ shocks, except that the former are slightly more energetic. Comparison of the shock Mach numbers at 1 AU shows that the RQ shocks are mostly subcritical, suggesting that they were not efficient in accelerating electrons. The Mach number values also indicate that most of these are quasi-perpendicular shocks. The radio-quietness is predominant in the rise phase and decreases through the maximum and declining phases of solar cycle 23. About 18% of the IP shocks do not have discernible ejecta behind them. These shocks are due to CMEs moving at large angles from the Sun-Earth line and hence are not blast waves. The solar sources of the shock-driving CMEs follow the sunspot butterfly diagram, consistent with the higher-energy requirement for driving shocks.

Gopalswamy, N.↗
NASA NTRS2010

Spitzer Observations of Dust Destruction in the Puppis A Supernova Remnant

The interaction of the Puppis A supernova remnant (SNR) with a neighboring molecular cloud provides a unique opportunity to measure the amount of grain destruction in an SNR shock. Spitzer Space Telescope MIPS imaging of the entire SNR at 24, 70, and 160 micrometers shows an extremely good correlation with X-ray emission, indicating that the SNR's IR radiation is dominated by the thermal emission of swept-up interstellar dust, collisionally heated by the hot shocked gas. Spitzer IRS spectral observations targeted both the Bright Eastern Knot (BEK) of the SNR where a small cloud has been engulfed by the supernova blast wave and outlying portions of the associated molecular cloud that are yet to be hit by the shock front. Modeling the spectra from both regions reveals the composition and the grain size distribution of the interstellar dust, both in front of and behind the SNR shock front. The comparison shows that the ubiquitous polycyclic aromatic hydrocarbons of the interstellar medium are destroyed within the BEK, along with nearly 25% of the mass of graphite and silicate dust grains.

Arendt, Richard G.↗
NASA NTRS2012

A Jet Break in the X-ray Light Curve of Short GRB 111020A: Implications for Energetics and Rates

We present broadband observations of the afterglow and environment of the short GRB 111020A. An extensive X-ray light curve from Swift/XRT, XMM-Newton, and Chandra, spanning approx.100 s to 10 days after the burst, reveals a significant break at (delta)t approx. = 2 days with pre- and post-break decline rates of (alpha)X,1 approx. = -0.78 and (alpha)X,2 < or approx. 1.7, respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle of (theta)j approx. = 3deg - 8deg. The resulting beaming-corrected gamma-ray (10-1000 keV band) and blast-wave kinetic energies are (2-3) x 10(exp 48) erg and (0.3-2) x 10(exp 49) erg, respectively, with the range depending on the unknown redshift of the burst. We report a radio afterglow limit of <39 micro-Jy (3(sigma)) from Expanded Very Large Array observations that, along with our finding that v(sub c) < v(sub X), constrains the circumburst density to n(sub 0) approx.0.01 0.1/cu cm. Optical observations provide an afterglow limit of i > or approx.24.4 mag at 18 hr after the burst and reveal a potential host galaxy with i approx. = 24.3 mag. The subarcsecond localization from Chandra provides a precise offset of 0".80+/-0".11 (1(sigma))from this galaxy corresponding to an offset of 5.7 kpc for z = 0.5-1.5. We find a high excess neutral hydrogen column density of (7.5+/-2.0) x 10(exp 21)/sq cm (z = 0). Our observations demonstrate that a growing fraction of short gamma-ray bursts (GRBs) are collimated, which may lead to a true event rate of > or approx.100-1000 Gpc(sup -3)/yr, in good agreement with the NS-NS merger rate of approx. = 200-3000 Gpc(sup -3)/ yr. This consistency is promising for coincident short GRB-gravitational wave searches in the forthcoming era of Advanced LIGO/VIRGO.

Fong, W.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
NASA NTRS2013

The Supercritical Pile Gamma-Ray Burst Model: The GRB Afterglow Steep Decline and Plateau Phase

We present a process that accounts for the steep decline and plateau phase of the Swift X-Ray Telescope (XRT) light curves, vexing features of gamma-ray burst (GRB) phenomenology. This process is an integral part of the "supercritical pile" GRB model, proposed a few years ago to account for the conversion of the GRB kinetic energy into radiation with a spectral peak at E(sub pk) is approx. m(sub e)C(exp 2). We compute the evolution of the relativistic blast wave (RBW) Lorentz factor Gamma to show that the radiation-reaction force due to the GRB emission can produce an abrupt, small (approx. 25%) decrease in Gamma at a radius that is smaller (depending on conditions) than the deceleration radius R(sub D). Because of this reduction, the kinematic criticality criterion of the "supercritical pile" is no longer fulfilled. Transfer of the proton energy into electrons ceases and the GRB enters abruptly the afterglow phase at a luminosity smaller by approx. m(sub p)/m(sub e) than that of the prompt emission. If the radius at which this slow-down occurs is significantly smaller than R(sub D), the RBW internal energy continues to drive the RBW expansion at a constant (new) Gamma and its X-ray luminosity remains constant until R(sub D) is reached, at which point it resumes its more conventional decay, thereby completing the "unexpected" XRT light curve phase. If this transition occurs at R is approx. equal to R(sub D), the steep decline is followed by a flux decrease instead of a "plateau," consistent with the conventional afterglow declines. Besides providing an account of these peculiarities, the model suggests that the afterglow phase may in fact begin before the RBW reaches R is approx. equal to R(sub D), thus providing novel insights into GRB phenomenology.

parameters↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
NASA NTRS2017

Asymmetric Expansion of the Youngest Galactic Supernova Remnant G1.9+0.3

The youngest Galactic supernova remnant (SNR) G1.9+0.3, produced by a (probable) SN Ia that exploded approximately 1900 CE, is strongly asymmetric at radio wavelengths, much brighter in the north, but bilaterally symmetric in X-rays. We present the results of X-ray expansion measurements that illuminate the origin of the radio asymmetry. We confirm the mean expansion rate (2011-2015) of 0.58% per yr, but large spatial variations are present. Using the nonparametric 'Demons' method, we measure the velocity field throughout the entire SNR, finding that motions vary by a factor of 5, from 0.''09 to 0.''44 per yr. The slowest shocks are at the outer boundary of the bright northern radio rim, with velocities v(sub s) as low as 3600 km per sec (for an assumed distance of 8.5 kpc), much less than v(sub s) = 12,000-13,000 km per sec along the X-ray-bright major axis. Such strong deceleration of the northern blast wave most likely arises from the collision of SN ejecta with a much denser than average ambient medium there. This asymmetric ambient medium naturally explains the radio asymmetry. In several locations, significant morphological changes and strongly nonradial motions are apparent. The spatially integrated X-ray flux continues to increase with time. Based on Chandra observations spanning 8.3 yr, we measure its increase at 1.3% +/- 0.8% per yr. The SN ejecta are likely colliding with the asymmetric circumstellar medium ejected by the SN progenitor prior to its explosion.

ISM: individual objects (G1.9+0.3) ISM: individual↗
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 NTRS2018

Atmospheric Energy Deposition Modeling and Inference for Varied Meteoroid Structures

Asteroids populations are highly diverse, ranging from coherent monoliths to loosely-bound rubble piles with a broad range of material and compositional properties. These different structures and properties could significantly affect how an asteroid breaks up and deposits energy in the atmosphere, and how much ground damage may occur from resulting blast waves. We have previously developed a fragment-cloud model (FCM) for assessing the atmospheric breakup and energy deposition of asteroids striking Earth. The approach represents ranges of breakup characteristics by combining progressive fragmentation with releases of variable fractions of debris and larger discrete fragments. In this work, we have extended the FCM to also represent asteroids with varied initial structures, such as rubble piles or fractured bodies. We have used the extended FCM to model the Chelyabinsk, Benesov, Kosice, and Tagish Lake meteors, and have obtained excellent matches to energy deposition profiles derived from their light curves. These matches provide validation for the FCM approach, help guide further model refinements, and enable inferences about pre-entry structure and breakup behavior. Results highlight differences in the amount of small debris vs. discrete fragments in matching the various flare characteristics of each meteor. The Chelyabinsk flares were best represented using relatively high debris fractions, while Kosice and Benesov cases were more notably driven by their discrete fragmentation characteristics, perhaps indicating more cohesive initial structures. Tagish Lake exhibited a combination of these characteristics, with lower-debris fragmentation at high altitudes followed by sudden disintegration into small debris in the lower flares. Results from all cases also suggest that lower ablation coefficients and debris spread rates may be more appropriate for the way in which debris clouds are represented in FCM, offering an avenue for future model refinement.

Deposition↗