Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiple shocks”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Oxidation States of Grim Glasses in EET79001 Based on Vanadium Valence

Gas-rich impact-melt (GRIM) glasses in SNC meteorites are very rich in Martian atmospheric noble gases and sulfur suggesting a possible occurrence of regolith-derived secondary mineral assemblages in these samples. Previously, we have studied two GRIM glasses, 506 and 507, from EET79001 Lith A and Lith B, respectively, for elemental abundances and spatial distribution of sulfur using EMPA (WDS) and FE-SEM (EDS) techniques and for sulfur-speciation using K-edge XANES techniques. These elemental and FE-SEM micro-graph data at several locations in the GRIM glasses from Shergotty (DBS), Zagami 994 and EET79001, Lith B showed that FeO and SO3 are positively correlated (SO3 represents a mixture of sulfide and sulfate). FE-SEM (EDS) study revealed that the sulfur-rich pockets in these glasses contain numerous micron-sized iron-sulfide (Fe-S) globules sequestered throughout the volume. However, in some areas (though less frequently), we detected significant Fe-S-O signals suggesting the occurrence of iron sulfate. These GRIM glasses were studied by K-edge microXANES techniques for sulfur speciation in association with iron in sulfur-rich areas. In both samples, we found the sulfur speciation dominated by sulfide with minor oxidized sulfur mixed in with various proportions. The abundance of oxidized sulfur was greater in 506 than in 507. Based on these results, we hypothesize that sulfur initially existed as sulfate in the glass precursor materials and, on shock-impact melting of the precursor materials producing these glasses, the oxidized sulfur was reduced to predominately sulfide. In order to further test this hypothesis, we have used microXANES to measure the valence states of vanadium in GRIM glasses from Lith A and Lith B to complement and compare with previous analogous measurements on Lith C (note: 506 and 507 contain the largest amounts of martian atmospheric gases but the gas-contents in Lith C measured by are unknown). Vanadium is ideal for addressing this re-dox issue because it has multiple valence states and is a well-studied element. Ferrous-dominated iron valences determined by microXANES on the Lith A and Lith B glasses provide little redox sensitivity. Vanadium valence measurements for impact glass in Lith C at three different locations yielded valence values of 3.1, 3.2 and 3.4 with inferred fO2 values of IW-0.7, IW-0.1 and IW+0.7, respectively. This range of oxygen-fugacity values is understandable because the glasses are shock-molten impact glasses which are heterogeneous in nature. Oxygen fugacity values obtained from the analysis of Fe-Ti oxides and Eu partitioning in pyroxenes from EET79001 Lith A and Lith B (host lithologies) were in the range of IW+0.3 to IW+1.9 suggesting that V in the Lith C impact glass was reduced in the impact process. Here, we examine whether the 506 from Lith A and 507 from Lith B GRIM glasses yield similar or different fO2 values from those of Lith C using the vanadium K-edge microXANES technique.

Sutton, S. R.↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
DOE OSTI2012

KNK437, abrogates hypoxia-induced radioresistance by dual targeting of the AKT and HIF-1{alpha} survival pathways

Highlights: Black-Right-Pointing-Pointer KNK437, a benzylidene lactam compound, is a novel radiosensitizer. Black-Right-Pointing-Pointer KNK437 inhibits AKT signaling and abrogates the accumulation of HIF-1{alpha} under hypoxia. Black-Right-Pointing-Pointer KNK437 abrogates hypoxia induced resistance to radiation. -- Abstract: KNK437 is a benzylidene lactam compound known to inhibit stress-induced synthesis of heat shock proteins (HSPs). HSPs promote radioresistance and play a major role in stabilizing hypoxia inducible factor-1{alpha} (HIF-1{alpha}). HIF-1{alpha} is widely responsible for tumor resistance to radiation under hypoxic conditions. We hypothesized that KNK437 sensitizes cancer cells to radiation and overrides hypoxia-induced radioresistance via destabilizing HIF-1{alpha}. Treatment of human cancer cells MDA-MB-231 and T98G with KNK437 sensitized them to ionizing radiation (IR). Surprisingly, IR did not induce HSPs in these cell lines. As hypothesized, KNK437 abrogated the accumulation of HIF-1{alpha} in hypoxic cells. However, there was no induction of HSPs under hypoxic conditions. Moreover, the proteosome inhibitor MG132 did not restore HIF-1{alpha} levels in KNK437-treated cells. This suggested that the absence of HIF-1{alpha} in hypoxic cells was not due to the enhanced protein degradation. HIF-1{alpha} is mainly regulated at the level of post-transcription and AKT is known to modulate the translation of HIF-1{alpha} mRNA. Interestingly, pre-treatment of cells with KNK437 inhibited AKT signaling. Furthermore, down regulation of AKT by siRNA abrogated HIF-1{alpha} levels under hypoxia. Interestingly, KNK437 reduced cell survival in hypoxic conditions and inhibited hypoxia-induced resistance to radiation. Taken together, these data suggest that KNK437 is an effective radiosensitizer that targets multiple pro-survival stress response pathways.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Interfacial Strength Evaluation of Oxide Films on Carbon Steel by Using the Laser Shock Adhesion Test

We investigate the adhesion strength (interfacial strength) of the oxide film formed on a carbon steel substrate by using a laser shock adhesion test (LaSAT). A rolled carbon steel (SS400) plate is first heated at high temperature (700 °C) under the atmosphere (heat treatment) to cause the formation of oxide films about 20-50 μm thick on the steel surface. The film consisted of multiple layers, including wustite (FeO), magnetite (Fe{sub 3}O{sub 4}), and hematite (Fe{sub 2}O{sub 3}). Next, the x-ray diffraction method was used to measure the residual stress in the film. We find that compression stress develops in the film because of a mismatch of thermal strain in between the film and substrate after the heat treatment. With LaSAT, strong elastic wave generates by a pulsed laser ablation of the substrate backside, and results in interfacial fractures in the oxide films. To monitor the normal (out-of-plane) displacement of the top surface film, a laser ultrasonic interferometer is simultaneously utilized, such that we detect film delamination and determine the critical laser energy for interfacial fracture to obtain the critical stress (interfacial strength). Furthermore, additional pulsed laser irradiation is conducted against a pre-delaminated specimen in order to examine the growth of delamination area and evaluate the interfacial fracture toughness. Finally, elastic wave propagation in the sample is computed using the finite element method (FEM) in order to evaluate the interfacial stress field, accurately. For the FEM, multiple layers in the oxide film are modeled, and the residual stress due to thermal strain mismatch is introduced. The FEM computations reveal the stress distribution around the interface and evaluate the critical stress and critical stress intensity factor for delamination. We thus quantitatively evaluate the interfacial strength and interfacial fracture toughness of the oxide film formed on carbon steel. The mechanism of interfacial fracture is also discussed on the basis of microstructures in a film with multiple layers. This result may be useful for understanding the adhesion quality of oxide film during the hot-rolling process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Testing, Characterization, and Modeling of the Resonant Plate Test Environment

The resonant plate shock test is a dynamic test of a mid-field pyroshock environment where a projectile is struck against a plate. The structure undergoing the simulated field shock is mounted to the plate. The plate resonates when struck and provides a two sided shock that is representative of the shock observed in the field. This test environment shock simulates a shock in a single coordinate direction for components looking to provide evidence that they will survive a similar or less shock when deployed in their operating environment. However, testing in one axis at a time provides many challenges. The true environment is a multi-axis environment. The test environment exhibits strong off-axis motion when only motion in one axis is desired. Multiple fixtures are needed for a single test series. It would be advantageous if a single test could be developed that tests the multi-axis environment simultaneously. In order to design such a test, a model must be developed and validated. The model can be iterated in design and configuration until the specified multi-axis environment is met. The test can then execute the model driven test design. This report discusses the resonant plate model needed to design future tests and the steps and methods used to obtain the model. This report also details aspects of the resonant plate test discovered during the process of model development that aids in our understanding of the test.

42 ENGINEERING↗
NASA NTRS1972

The simulation of high pressure hydrogen/oxygen rocket engines.

Initial experiments have established the feasibility of using a backward running detonation wave in a shock tube to produce the correct enthalpy and chemistry of the simulated combustion chamber gases. This technique produces plume Reynolds numbers higher than those achievable by any other known technique. The practical pressure limitations are being explored and detailed surveys of both single and multiple engine plumes and impingement data on a cylindrical body are being obtained. Facility details to allow safe handling of the high pressure hydrogen and oxygen up to a design chamber pressure of 20,000 psi are described.

Hopkins, H.↗
NASA NTRS1990

Efficient parametric analysis of performance measures for communication networks

Efficient techniques for estimating performance measures in communication networks in a steady-state or transient setting are developed. These techniques may be used in a simulation environment or in connection with real-time observations. For Markov chain models, the recently proposed standard clock approach is extended, and a class of real-time algorithms for simultaneously generating multiple sample paths under different parameter sets is presented. Attention is focused on the link crash time estimation problem, where a 'crash' is defined as the first time a buffer overflows, given some initial conditions. An algorithm for estimating crash times under various traffic shocks is derived, where all estimates are obtained in parallel to an actual network's normal operation. An algorithm for crash time estimation of a G/D/1 link model is also derived using a different (perturbation-analysis-based) approach. Finally, extensive simulation results are provided to validate the proposed algorithms and compare them to brute-force simulation.

Cassandras, Christos G.↗
DOE OSTI2013

THE PROGENITOR OF SN 2011ja: CLUES FROM CIRCUMSTELLAR INTERACTION

Massive stars, possibly red supergiants, which retain extended hydrogen envelopes until core collapse, produce Type II plateau (IIP) supernovae. The ejecta from these explosions shocks the circumstellar matter originating from the mass loss of the progenitor during the final phases of its life. This interaction accelerates particles to relativistic energies which then lose energy via synchrotron radiation in the shock-amplified magnetic fields and inverse Compton scattering against optical photons from the supernova. These processes produce different signatures in the radio and X-ray parts of the electromagnetic spectrum. Observed together, they allow us to break the degeneracy between shock acceleration and magnetic field amplification. In this work, we use X-rays observations from the Chandra and radio observations from the Australia Telescope Compact Array to study the relative importance of processes which accelerate particles and those which amplify magnetic fields in producing the non-thermal radiation from SN 2011ja. We use radio observations to constrain the explosion date. Multiple Chandra observations allow us to probe the history of variable mass loss from the progenitor. The ejecta expands into a low-density bubble followed by interaction with a higher density wind from a red supergiant consistent with M{sub ZAMS} {approx}> 12 M{sub Sun }. Our results suggest that a fraction of Type IIP supernovae may interact with circumstellar media set up by non-steady winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Resistive magnetohydrodynamic reconnection: Resolving long-term, chaotic dynamics

In this paper, we address the long-term evolution of an idealised double current system entering reconnection regimes where chaotic behavior plays a prominent role. Our aim is to quantify the energetics in high magnetic Reynolds number evolutions, enriched by secondary tearing events, multiple magnetic island coalescence, and compressive versus resistive heating scenarios. Our study will pay particular attention to the required numerical resolutions achievable by modern (grid-adaptive) computations, and comment on the challenge associated with resolving chaotic island formation and interaction. We will use shock-capturing, conservative, grid-adaptive simulations for investigating trends dominated by both physical (resistivity) and numerical (resolution) parameters, and confront them with (visco-)resistive magnetohydrodynamic simulations performed with very different, but equally widely used discretization schemes. This will allow us to comment on the obtained evolutions in a manner irrespective of the adopted discretization strategy. Our findings demonstrate that all schemes used (finite volume based shock-capturing, high order finite differences, and particle in cell-like methods) qualitatively agree on the various evolutionary stages, and that resistivity values of order 0.001 already can lead to chaotic island appearance. However, none of the methods exploited demonstrates convergence in the strong sense in these chaotic regimes. At the same time, nonperturbed tests for showing convergence over long time scales in ideal to resistive regimes are provided as well, where all methods are shown to agree. Both the advantages and disadvantages of specific discretizations as applied to this challenging problem are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Multiple-cycle Simulation of a Pulse Detonation Engine Ejector

This paper presents the results of a study involving single and multiple-cycle numerical simulations of various PDE-ejector configurations utilizing hydrogen-oxygen mixtures. The objective was to investigate the thrust, impulse and mass flow rate characteristics of these devices. The results indicate that ejector systems can utilize the energy stored in the strong shock wave exiting the detonation tube to augment the impulse obtained from the detonation tube alone. Impulse augmentation ratios of up to 1.9 were achieved. The axial location of the converging-diverging ejectors relative to the end of the detonation tube were shown to affect the performance of the system.

Yungster, S.↗
DOE OSTI2016

Development of an FPGA-based multipoint laser pyroshock measurement system for explosive bolts

Pyroshock can cause failure to the objective of an aerospace structure by damaging its sensitive electronic equipment, which is responsible for performing decisive operations. A pyroshock is the high intensity shock wave that is generated when a pyrotechnic device is explosively triggered to separate, release, or activate structural subsystems of an aerospace architecture. Pyroshock measurement plays an important role in experimental simulations to understand the characteristics of pyroshock on the host structure. This paper presents a technology to measure a pyroshock wave at multiple points using laser Doppler vibrometers (LDVs). These LDVs detect the pyroshock wave generated due to an explosive-based pyrotechnical event. Field programmable gate array (FPGA) based data acquisition is used in the study to acquire pyroshock signals simultaneously from multiple channels. This paper describes the complete system design for multipoint pyroshock measurement. The firmware architecture for the implementation of multichannel data acquisition on an FPGA-based development board is also discussed. An experiment using explosive bolts was configured to test the reliability of the system. Pyroshock was generated using explosive excitation on a 22-mm-thick steel plate. Three LDVs were deployed to capture the pyroshock wave at different points. The pyroshocks captured were displayed as acceleration plots. The results showed that our system effectively captured the pyroshock wave with a peak-to-peak magnitude of 303 741 g. The contribution of this paper is a specialized architecture of firmware design programmed in FPGA for data acquisition of large amount of multichannel pyroshock data. The advantages of the developed system are the near-field, multipoint, non-contact, and remote measurement of a pyroshock wave, which is dangerous and expensive to produce in aerospace pyrotechnic tests.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Development of Glucose Regulated Protein 94-Selective Inhibitors Based on the BnIm and Radamide Scaffold

Glucose regulated protein 94 (Grp94) is the endoplasmic reticulum resident of the heat shock protein 90 kDa (Hsp90) family of molecular chaperones. Grp94 associates with many proteins involved in cell adhesion and signaling, including integrins, Toll-like receptors, immunoglobulins, and mutant myocilin. Grp94 has been implicated as a target for several therapeutic areas including glaucoma, cancer metastasis, and multiple myeloma. While 85% identical to other Hsp90 isoforms, the N-terminal ATP-binding site of Grp94 possesses a unique hydrophobic pocket that was used to design isoform-selective inhibitors. Incorporation of a cis-amide bioisostere into the radamide scaffold led to development of the original Grp94-selective inhibitor, BnIm. Structure–activity relationship studies have now been performed on the aryl side chain of BnIm, which resulted in improved analogues that exhibit better potency and selectivity for Grp94. Furthermore, these analogues also manifest superior antimigratory activity in a metastasis model as well as enhanced mutant myocilin degradation in a glaucoma model compared to BnIm.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Global Prompt Proton Sensor Network: Monitoring Solar Energetic Protons Based on GPS Satellite Constellation

Energetic particle instruments on board Global Positioning System (GPS) satellites form a powerful global prompt proton sensor network (GPPSn) that provides an unprecedented opportunity to monitor and characterize solar energetic protons targeting the Earth. The medium-Earth orbits of the GPS constellation have the unique advantage of allowing solar energetic protons to be simultaneously measured from multiple points in both open- and closed-field line regions. Examining two example intervals of solar proton events, we showcase in this study how GPS proton data are prepared, calibrated, and utilized to reveal important features of solar protons, including their source, acceleration/scattering by interplanetary shocks, the relative position of Earth when impinged by these shocks, the shape of solar particle fronts, the access of solar protons inside the dynamic geomagnetic field, and temporally varying proton distributions in both energy and space. By comparing to Van Allen Probes data, GPS proton observations are further demonstrated not only to be useful for qualitatively monitoring the dynamics of solar protons but also for quantitative scientific research including determining cutoff L-shells. Finally, our results establish that this GPPSn can join forces with other existing solar proton monitors and contribute to observing, warning, understanding, and ultimately forecasting the incoming solar energetic proton events.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Particle-in-cell Simulations of Mildly Relativistic Outflows in Kilonova Emissions

Abstract The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet and thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation, and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL’s vector particle-in-cell (VPIC) code, we have run simulations of such mildly relativistic, collisionless, weakly magnetized plasmas and computed the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton-to-electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗