Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK WAVE PROFILE”

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

High-Resolution N-Band Observations of the Nova Rs Ophiuchi: First Science with the Keck Interferometer Nuller

We report new observations of the nova RS Ophiuchi using the Keck Interferometer Nulling Instrument (KIN), approximately 3.8 days following the most recent outburst that occurred on 2006 February 12. These observations represent the first scientific results from the KIN, which operates in N-band from 8 to 12.5 microns. The nulling technique is the sparse aperture equivalent of the conventional coronagraphic technique used in filled aperture telescopes. By fitting the unique KIN inner and outer spatial regime data, we have obtained an angular size of the mid-infrared continuum of 6.2,4.0, or 5.4 mas for a disk profile, Gaussian profile, and shell profile respectively. The data show evidence of enhanced neutral atomic hydrogen emission and atomic metals including silicon located in the inner spatial regime near the white dwarf relative to the outer regime. There are also nebular emission lines and evidence of hot silicate dust in the outer spatial region, centered at approximately 2.5E14 cm from the WD, that are not found in the inner regime. The nova flash in the outer spatial regime evidently excited these features before the blast wave reached these regions. These identifications support the following interpretation. The dust appears to be present between outbursts and was not created during the outburst event. We further discuss the present results in terms of a unifying model of the system that includes an increase in density in the plane of the orbit of the two stars created by a spiral shock wave caused by the motion of the stars through the cool wind of the red giant star. These data show the power and potential of the nulling technique which has been developed for the detection of Earthlike planets around nearby stars for the Terrestrial Planet Finder Mission and Darwin missions.

Danchi, William↗
DOE OSTI2025

Calibration and characterization of the line-VISAR diagnostic at the HED-HIBEF instrument at the European XFEL

In dynamic-compression experiments, the line-imaging Velocity Interferometer System for Any Reflector (VISAR) is a well-established diagnostic used to probe the velocity history, including wave profiles derived from dynamically compressed interfaces and wavefronts, depending on material optical properties. Knowledge of the velocity history allows for the determination of the pressure achieved during compression. Such a VISAR analysis is often based on Fourier transform techniques and assumes that the recorded interferograms are free from image distortions. In this paper, we describe the VISAR diagnostic installed at the HED-HIBEF instrument located at the European XFEL along with its calibration and characterization. It comprises a two-color (532, 1064 nm), three-arm (with three velocity sensitivities) line imaging system. We provide a procedure to correct VISAR images for geometric distortions and evaluate the performance of the system using Fourier analysis. We finally discuss the spatial and temporal calibrations of the diagnostic. As an example, we compare the pressure extracted from the VISAR analysis of shock-compressed polyimide and silicon.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
NASA NTRS1994

The Balmer-dominated northeast limb of the Cygnus loop supernova remnant

We present a comprehensive investigation of the Balmar-dominated northeast limb of the Cygnus Loop supernova remnant. Data presented include H alpha (O III), and X-ray images, UV and visible spectrophotometry, and high-resolution spectroscophy. The two relatively bright Balmer-dominated filaments visible on the POSS prints are seen to be part of a very smooth and regular complex of filaments. These filaments mark the current location of the blast wave and are seen to bound the sharply limb-brightened X-ray emission, including the previously reported X-ray, 'halo.' The (O III)/h beta ratio throughout the region is approximately 0.1, except for regions in which the shock is undergoing a transition from nonradiative to incomplete radiative to incomplete radiative conditions. At these locations (O III) emission from the cooling region is quite strong, while collisionally excited Balmer-line emission can be weak because of photoionization of the preshock medium by UV from the nascent cooling region. As a result (O III)/H beta is greater than 100 in some locations. The nonradiative/radiative transition is best studied along the length of the northwestern of the two brightest filaments, where the shock velocity and swept-up column go from approximately 180 km/s and 10(exp 17)/sq cm at one end to approximately 140 km/s and 8 x 10(exp 17)/cm at the other. There are also a number of locations of such incomplete radiative emission where the shock has recently encountered denser regions with characteristic sizes of approximately 10(exp 18) cm. There is a considerable amount of evidence that the shock has decelerated from approximately 400 km/s to less than 200 km/s in the last 1000 yr. We interpret this as the result of the blast wave hitting the wall of a cavity which surround the supernova precursor and succeed in matching a wide range of data with a reflected shock model in which the density ofthe cavity wall is approximately 1.2/cu cm and the density in the interior of the cavity is about 0.09/cu cm. This analysis supports previous suggestions that the Cygnus Loop was a cavity explosion. We opredict that the X-ray shell will continue to brighten for the next 1000 years or so, during which time the shock front will become radiative throughout the northeast region. The similarity of the time scales for post-shock cooling, shock deceleration, and the 1ressure evolution of the region violates assumptions of steady flow, and greatly complicated interpretation of shock specptra in the region. This may be a common and important effect in many supernova remnants. High-resolutoin spectroscopy of the H alpha line profile shows that the broad component is surprisingly narrow (approximately 130 km/s). Consistency between this line width and the shock velocity inferred from spectrophotometry requires either rapid postshock equilibration of T(sub i) and T(sub e), or that the shock front has declerated from 180 km/s to 150 km/s in the last approximately 200 yr. While we cannot unambiguously distinquish between these models, the required deceleration is consistent with extrapolation of the shock decleration rate inferred over the last 1000 years. The width of the narrow component of the h alpha line profile is approximately 3 km/s, which indicates that the preshock medium has been heated to approximately 25,000 K by a shock precursor. The most promising types of precursor are heating either by fast neutrals which overtake the shock or by cosmic rays. In the latter case the observed properties of the precursor constrain the cosmic-ray diffusion coefficient to be between about 10(exp 23) sq cm/s. The characteristics of the precursor seem to be about right to explain the presence of weak (S II) and (N II) emission previously reported by Fesen & Itoh (1985).

Hester, J. Jeff↗
NASA NTRS2023

The Faraday Effect Tracker of Coronal and Heliospheric Structure (FETCH) Instrument

There continue to be open questions regarding the solar wind and coronal mass ejections (CMEs). For example: how do magnetic fields within CMEs and corotating/stream interaction regions (CIRs/SIRs) evolve in the inner heliosphere? What is the radially distributed magnetic profile of shock-driving CMEs? What is the internal magnetic structure of CMEs that cause magnetic storms? It is clear that these questions involve the magnetic configurations of solar wind and transient interplanetary plasma structures, for which we have limited knowledge. In order to better understand the origin of the magnetic field variability in steady-state structures and transient events, it is necessary to probe the magnetic field in Earth-directed structures/disturbances. This is the goal of the Multiview Observatory for Solar Terrestrial Science (MOST) mission (Gopalswamy et al., 2022). For MOST to answer the aforementioned questions, we propose the instrument concept of the Faraday Effect Tracker of Coronal and Heliospheric structures (FETCH), a simultaneous quad-line-of-sight polarization radio remote-sensing instrument. With FETCH, spacecraft radio beams passing through the Sun–Earth line offer the possibility of obtaining information of plasma conditions via analysis of radio propagation effects such as Faraday rotation and wave dispersion, which provide information of the magnetic field and total electron content (TEC). This is the goal of the FETCH instrument, one of ten instruments proposed to be hosted on the MOST mission. The MOST mission will provide an unprecedented opportunity to achieve NASA’s heliophysics science goal to “explore and characterize the physical processes in the space environment from the Sun” (Gopalswamy et al., 2022).

Solar corona↗
DOE OSTI2023

Performance Characteristics of a Rotating Detonation Combustor Exiting Into a Pressurized Plenum to Simulate Gas Turbine Inlet

Abstract The present study aims to experimentally characterize the performance of a rotating detonation combustion (RDC) system integrated with a pressurized downstream plenum to simulate the high-pressure inlet conditions of power-generating gas turbines. A thorough understanding of the operational behavior including wave mode behavior, static pressure profile along the combustor length, and dynamic features of pressure fluctuations is crucial for successful integration of RDC with the turbine. In this study, two RDC configurations are investigated, RDC with a constant area annulus and RDC with a converging nozzle. In both cases, the RDC flow exited into a plenum chamber kept at pressures varying from 155 kPa to 330 kPa. RDC was operated on methane and oxygen-enriched air to represent reactants used in land-based power generation. Experiments were conducted for the two RDCs configurations operated at three reactant mass flow rates (0.23, 0.32, and 0.46 kg/s). The RDC performance is characterized by time-averaged static pressure measurements, and wave velocity determined by ionization probes. In addition, dynamic pressure measurements were recorded both inside and near the exit of RDC channel to investigate wave interactions between RDC and downstream plenum. Results show that the RDC with the converging nozzle achieved superior performance while minimizing detrimental interactions with the reflected shock and/or acoustic waves from the downstream plenum.

Engineering↗
NASA NTRS2019

Evolution of the Suprathermal Proton Population at Interplanetary Shocks

We investigate the evolution of the suprathermal (ST) proton population as interplanetary shocks cross 1 au. The variability of the ST proton intensities and energy spectra upstream of the shocks is analyzed in terms of the shock parameters, upstream magnetic field configurations, and preexisting upstream populations. Propitious conditions for the observation of ST particles at distances far upstream from the shock occur in parallel shock configurations when particles can easily escape from the shock vicinity. In this situation, ST intensity enhancements show onsets characterized by velocity dispersion effects and energy spectra that develop into a “hump” profile peaking around ∼10 keV just before the arrival of the shock. The observation of field-aligned proton beams at low energies (5–10 keV) is possible under conditions that facilitate the scatter-free propagation of the particles streaming out of the shock. Upstream of perpendicular shocks, ST intensity enhancements are only observed in close proximity to the shock. Power-law proton spectra develop downstream of the shocks. The functional form for the downstream phase-space density proportional to v(exp -5) is observed only over a limited range of ST energies. The absence of ST populations observed far upstream of interplanetary shocks raises questions about whether ST protons contribute as a seed particle population in the processes of particle acceleration at shocks.

Shock waves↗
DOE OSTI2024

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Pressure distribution for the wing of the YAV-8B airplane; with and without pylons

Pressure distribution data have been obtained in flight at four span stations on the wing panel of the YAV-8B airplane. Data obtained for the supercritical profiled wing, with and without pylons installed, ranged from Mach 0.46 to 0.88. The altitude ranged from approximately 20,000 to 40,000 ft and the resultant Reynolds numbers varied from approximately 7.2 million to 28.7 million based on the mean aerodynamic chord. Pressure distribution data and flow visualization results show that the full-scale flight wing performance is compromised because the lower surface cusp region experiences flow separation for some important transonic flight conditions. This condition is aggravated when local shocks occur on the lower surface of the wing (mostly between 20 and 35 percent chord) when the pylons are installed for Mach 0.8 and above. There is evidence that convex fairings, which cover the pylon attachment flanges, cause these local shocks. Pressure coefficients significantly more negative than those for sonic flow also occur farther aft on the lower surface (near 60 percent chord) whether or not the pylons are installed for Mach numbers greater than or equal to 0.8. These negative pressure coefficient peaks and associated local shocks would be expected to cause increasing wave and separation drag at transonic Mach number increases.

Saltzman, Edwin J.↗
DOE OSTI2022

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Wind Tunnel Model Design for the Study of Plume Effects on Sonic Boom for Isolated Exhaust Nozzles

A low cost test capability was developed at the NASA Glenn Research Center 1- by 1-Foot Supersonic Wind Tunnel (SWT), with a goal to reduce the disturbance caused by supersonic aircraft flight over populated areas. This work focused on the shock wave structure caused by the exhaust nozzle plume. Analysis and design was performed on a new rig to test exhaust nozzle plume effects on sonic boom signature. Test capability included a baseline nozzle test article and a wind tunnel model consisting of a strut, a nosecone and an upper plenum. Analysis was performed on the external and internal aerodynamic configuration, including the shock reflections from the wind tunnel walls caused by the presence of the model nosecone. This wind tunnel model was designed to operate from Mach 1.4 to Mach 3.0 with nozzle pressure ratios from 6 to 12 and altitudes from 30,000 ft (4.36 psia) to 50,000 ft (1.68 psia). The model design was based on a 1 in. outer diameter, was 9 in. in overall length, and was mounted in the wind tunnel on a 3/8 in. wide support strut. For test conditions at 50,000 ft the strut was built to supply 90 psia of pressure, and to achieve 20 psia at the nozzle inlet with a maximum nozzle pressure of 52 psia. Instrumentation was developed to measure nozzle pressure ratio, and an external static pressure probe was designed to survey near field static pressure profiles at one nozzle diameter above the rig centerline. Model layout placed test nozzles between two transparent sidewalls in the 1 1 SWT for Schlieren photography and comparison to CFD analysis.

Castner, Raynold S.↗
NASA NTRS2009

Wind Tunnel Model Design for the Study of Plume Effects on Sonic Boom for Isolated Exhaust Nozzles

A low cost test capability was developed at the NASA Glenn Research Center 1- by 1-Foot Supersonic Wind Tunnel (SWT), with a goal to reduce the disturbance caused by supersonic aircraft flight over populated areas. This work focused on the shock wave structure caused by the exhaust nozzle plume. Analysis and design was performed on a new rig to test exhaust nozzle plume effects on sonic boom signature. Test capability included a baseline nozzle test article and a wind tunnel model consisting of a strut, a nose cone and an upper plenum. Analysis was performed on the external and internal aerodynamic configuration, including the shock reflections from the wind tunnel walls caused by the presence of the model nosecone. This wind tunnel model was designed to operate from Mach 1.4 to Mach 3.0 with nozzle pressure ratios from 6 to 12 and altitudes from 30,000 ft (4.36 psia) to 50,000 ft (1.68 psia). The model design was based on a 1 in. outer diameter, was 9 in. in overall length, and was mounted in the wind tunnel on a 3/8 in. wide support strut. For test conditions at 50,000 ft the strut was built to supply 90 psia of pressure, and to achieve 20 psia at the nozzle inlet with a maximum nozzle pressure of 52 psia. Instrumentation was developed to measure nozzle pressure ratio, and an external static pressure probe was designed to survey near field static pressure profiles at one nozzle diameter above the rig centerline. Model layout placed test nozzles between two transparent sidewalls in the 1x1 SWT for Schlieren photography and comparison to CFD analysis.

Castner, Raymond S.↗
DOE OSTI2026

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cocoon shock breakout emission from binary neutron star mergers

Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Accretion Disks Around Binary Black Holes of Unequal Mass: GRMHD Simulations Near Decoupling

We report on simulations in general relativity of magnetized disks onto black hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the systems when they orbit near the binary disk decoupling radius. We compare (surface) density profiles, accretion rates (relative to a single, non-spinning black hole), variability, effective alpha-stress levels and luminosities as functions of the mass ratio. We treat the disks in two limiting regimes: rapid radiative cooling and no radiative cooling. The magnetic field lines clearly reveal jets emerging from both black hole horizons and merging into one common jet at large distances. The magnetic fields give rise to much stronger shock heating than the pure hydrodynamic flows, completely alter the disk structure, and boost accretion rates and luminosities. Accretion streams near the horizons are among the densest structures; in fact, the 1:10 no-cooling evolution results in a refilling of the cavity. The typical effective temperature in the bulk of the disk is approx. 10(exp5) (M / 10(exp 8)M solar mass (exp -1/4(L/L(sub edd) (exp 1/4K) yielding characteristic thermal frequencies approx. 10 (exp 15) (M /10(exp 8)M solar mass) (exp -1/4(L/L (sub edd) (1+z) (exp -1)Hz. These systems are thus promising targets for many extragalactic optical surveys, such as LSST, WFIRST, and PanSTARRS.

Supermassive↗
DOE OSTI2023

Magnetically driven coupling in relativistic radiation-mediated shocks

ABSTRACT The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron–positron pairs and a heavy ion species. We use a 1D five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter σc, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with γu = 10 in the upstream of a relativistic radiation-mediated shock, we find the threshold σc ≈ 10−7 for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of σc. Identifying high-energy signatures of microturbulence in shock breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers.

79 ASTRONOMY AND ASTROPHYSICS↗