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.

SourceResearch connectionFrom yearThrough year

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

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Multiple Components in the Broadband γ -Ray Emission of the Short GRB 160709A

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The γ-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole γ-ray spectrum. While the thermal component dominates during ~0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of ~340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, ${R}_{0}={3.8}_{-1.8}^{+5.9}$ × 10 8 cm, the size of the photosphere, ${R}_{\mathrm{ph}}={7.4}_{-1.2}^{+0.8}$ × 10 10 cm, and a bulk Lorentz factor, ${\rm{\Gamma }}={728}_{-93}^{+75}$, assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV–40 MeV and 100 MeV–10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, p = 2.5 ± 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1976

Components in interstellar molecular hydrogen

Results are reported for precise spectrophotometric measurements of the profiles of selected Lyman absorption lines produced by hydrogen molecules in various rotational levels along the line of sight to 13 stars which have shown some evidence for an increase in line width with increasing rotational quantum number (J). The line profiles were measured by multiple scans with the Copernicus satellite telescope. Based on analysis of the radial velocities, derivations of the column densities, and line-profile fitting, the following conclusions are made: (1) the increase in interstellar H2 line width with increasing J results from the presence of the most shortward component, which is relatively weak at low J but becomes more important at higher J; (2) the relative column densities found for the different J levels in each component may be fitted by a theoretical model in which rotational excitation is due to absorption of UV photons followed by radiative quadrupole spontaneous transitions or collisionally induced downward transitions between different J levels; (3) the atomic hydrogen density is between 300 and 1000 per cu cm in the most shortward component for each of three stars; (4) the approaching gas which produces each shortward component must be in the form of thin sheets; and (5) the sheets are the compressed gas behind a shock front moving through the interstellar medium.

Spitzer, L., Jr.↗

DOE OSTI2023

Radiation Transport Two-temperature GRMHD Simulations of Warped Accretion Disks

In many black hole (BH) systems, the accretion disk is expected to be misaligned with respect to the BH spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the BH can tear the disk into multiple, independently precessing "sub-disks." This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGNs) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations in BHXRBs and changing-look phenomena in AGNs. Here, we present the first radiative two-temperature general relativistic magnetohydrodynamic (GRMHD) simulation of a strongly tilted (65°) accretion disk around an M BH = 10 M ⊙ BH, which tears and precesses. This leads to luminosity swings between a few percent and 50% of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation-pressure-dominated, the accretion disk is thermally stable over t ≳ 14,000 r g /c. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases tenfold and the electron temperature reaches T e ~ 10 8 –10 9 K. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to T e ~ 108 K. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

OBSERVATIONS OF POST-FLARE PLASMA DYNAMICS DURING AN M1.0 FLARE IN AR11093 BY THE SOLAR DYNAMICS OBSERVATORY/ATMOSPHERIC IMAGING ASSEMBLY

We observe the motion of cool and hot plasma in a multi-stranded post-flare loop (PFL) system that evolved in the decay phase of a two-ribbon M1.0 class flare in AR 11093 on 2010 August 7 using the Solar Dynamics Observatory/Atmospheric Imaging Assembly 304 A and 171 A filters. The moving intensity feature and its reflected counterpart are observed in the loop system at multiple temperatures. The observed hot counterpart of the plasma probably envelopes the cool confined plasma and moves comparatively faster ({approx}34 km s{sup -1}) than the latter (29 km s{sup -1}) in the form of a spreading intensity feature. The propagating plasma and intensity reflect from the region of another footpoint of the loop. The subsonic speed of the moving plasma and associated intensity feature may be most likely evolved in the PFL system through impulsive flare heating processes. Complementing our observations of moving multi-temperature intensity features in the PFL system and its reflection, we also attempt to solve two-dimensional ideal magnetohydrodynamic equations numerically using the VAL-IIIC atmosphere as an initial condition to simulate the observed plasma dynamics. We consider a localized thermal pulse impulsively generated near one footpoint of the loop system during the flare processes, which is launched along the magnetic field lines at the solar chromosphere. The pulse steepens into a slow shock at higher altitudes while moving along this loop system, which triggers plasma perturbations that closely exhibit the observed plasma dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE DUST BUDGET OF THE SMALL MAGELLANIC CLOUD: ARE ASYMPTOTIC GIANT BRANCH STARS THE PRIMARY DUST SOURCE AT LOW METALLICITY?

We estimate the total dust input from the cool evolved stars in the Small Magellanic Cloud, using the 8 {mu}m excess emission as a proxy for the dust-production rate (DPR). We find that asymptotic giant branch (AGB) and red supergiant (RSG) stars produce (8.6-9.5) Multiplication-Sign 10{sup -7} M{sub Sun} yr{sup -1} of dust, depending on the fraction of far-infrared sources that belong to the evolved star population (with 10%-50% uncertainty in individual DPRs). RSGs contribute the least (<4%), while carbon-rich AGB stars (especially the so-called extreme AGB stars) account for 87%-89% of the total dust input from cool evolved stars. We also estimate the dust input from hot stars and supernovae (SNe), and find that if SNe produce 10{sup -3} M{sub Sun} of dust each, then the total SN dust input and AGB input are roughly equivalent. We consider several scenarios of SN dust production and destruction and find that the interstellar medium (ISM) dust can be accounted for solely by stellar sources if all SNe produce dust in the quantities seen around the dustiest examples and if most SNe explode in dense regions where much of the ISM dust is shielded from the shocks. We find that AGB stars contribute only 2.1% of the ISM dust. Without a net positive contribution from SNe to the dust budget, this suggests that dust must grow in the ISM or be formed by another unknown mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2020

Bayesian Statistical Models for Community Annoyance Survey Data

This paper demonstrates the use of two Bayesian statistical models to analyze single-event sonic boom exposure and human annoyance data from community response surveys. Each model is fit to data from a NASA pilot study.Unlike many community noise surveys, this study used a panel sample to collect multiple observations per participant instead of a single observation. Thus, a multilevel (also known as hierarchical or mixed-effects) model is used to account for the within-subject correlation in the panel sample data. This paper describes a multilevel logistic regression model and a multilevel ordinal regression model. The paper also proposes a method for calculating a summary dose-response curve from the multilevel models that represents the population. The two models’ summary dose-response curves are visually similar. However, their estimates differ when calculating the noise dose at a fixed percent highly annoyed.

Musical instruments↗

NASA NTRS2014

Agent Based Modeling of Air Carrier Behavior for Evaluation of Technology Equipage and Adoption

As part of ongoing research, the National Aeronautics and Space Administration (NASA) and LMI developed a research framework to assist policymakers in identifying impacts on the U.S. air transportation system (ATS) of potential policies and technology related to the implementation of the Next Generation Air Transportation System (NextGen). This framework, called the Air Transportation System Evolutionary Simulation (ATS-EVOS), integrates multiple models into a single process flow to best simulate responses by U.S. commercial airlines and other ATS stakeholders to NextGen-related policies, and in turn, how those responses impact the ATS. Development of this framework required NASA and LMI to create an agent-based model of airline and passenger behavior. This Airline Evolutionary Simulation (AIRLINE-EVOS) models airline decisions about tactical airfare and schedule adjustments, and strategic decisions related to fleet assignments, market prices, and equipage. AIRLINE-EVOS models its own heterogeneous population of passenger agents that interact with airlines; this interaction allows the model to simulate the cycle of action-reaction as airlines compete with each other and engage passengers. We validated a baseline configuration of AIRLINE-EVOS against Airline Origin and Destination Survey (DB1B) data and subject matter expert opinion, and we verified the ATS-EVOS framework and agent behavior logic through scenario-based experiments. These experiments demonstrated AIRLINE-EVOS's capabilities in responding to an input price shock in fuel prices, and to equipage challenges in a series of analyses based on potential incentive policies for best equipped best served, optimal-wind routing, and traffic management initiative exemption concepts..

Horio, Brant M.↗

DOE OSTI2015

Development of a Multi-Point Microwave Interferometry (MPMI) Method

A multi-point microwave interferometer (MPMI) concept was developed for non-invasively tracking a shock, reaction, or detonation front in energetic media. Initially, a single-point, heterodyne microwave interferometry capability was established. The design, construction, and verification of the single-point interferometer provided a knowledge base for the creation of the MPMI concept. The MPMI concept uses an electro-optic (EO) crystal to impart a time-varying phase lag onto a laser at the microwave frequency. Polarization optics converts this phase lag into an amplitude modulation, which is analyzed in a heterodyne interferometer to detect Doppler shifts in the microwave frequency. A version of the MPMI was constructed to experimentally measure the frequency of a microwave source through the EO modulation of a laser. The successful extraction of the microwave frequency proved the underlying physical concept of the MPMI design, and highlighted the challenges associated with the longer microwave wavelength. The frequency measurements made with the current equipment contained too much uncertainty for an accurate velocity measurement. Potential alterations to the current construction are presented to improve the quality of the measured signal and enable multiple accurate velocity measurements.

42 ENGINEERING↗

NASA NTRS2012

Prediction of Broadband Shock-Associated Noise Including Propagation Effects Originating NASA

An acoustic analogy is developed based on the Euler equations for broadband shock-associated noise (BBSAN) that directly incorporates the vector Green s function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) to describe the mean flow. The vector Green s function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation rate. An adjoint vector Green s function solver is implemented to determine the vector Green s function based on a locally parallel mean flow at different streamwise locations. The newly developed acoustic analogy can be simplified to one that uses the Green s function associated with the Helmholtz equation, which is consistent with a previous formulation by the authors. A large number of predictions are generated using three different nozzles over a wide range of fully-expanded jet Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise experimental facilities. In addition, two models for the so-called fine-scale mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include propagation effects.

Miller, Steven↗

DOE OSTI2018

Elucidating a role for the cytoplasmic domain in the Mycobacterium tuberculosis mechanosensitive channel of large conductance

Microbial survival in dynamic environments requires the ability to successfully respond to abrupt changes in osmolarity. The mechanosensitive channel of large conductance (MscL) is a ubiquitous channel that facilitates the survival of bacteria and archaea under severe osmotic down-shock conditions by relieving excess turgor pressure in response to increased membrane tension. A prominent structural feature of MscL, the cytoplasmic C-terminal domain, has been suggested to influence channel assembly and function. In this report, we describe the X-ray crystal structure and electrophysiological properties of a C-terminal domain truncation of the Mycobacterium tuberculosis MscL (MtMscLΔC). A crystal structure of MtMscLΔC solubilized in the detergent n-dodecyl-β-D-maltopyranoside reveals the pentameric, closed state-like architecture for the membrane spanning region observed in the previously solved full-length MtMscL. Electrophysiological characterization demonstrates that MtMscLΔC retains mechano-sensitivity, but with conductance and tension sensitivity more closely resembling full length EcMscL than MtMscL. This study establishes that the C-terminal domain of MtMscL is not required for oligomerization of the full-length channel, but rather influences the tension sensitivity and conductance properties of the channel. The collective picture that emerges from these data is that each MscL channel structure has characteristic features, highlighting the importance of studying multiple homologs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2012

CHARACTERIZING TRANSITION TEMPERATURE GAS IN THE GALACTIC CORONA

We present a study of the properties of the transition temperature (T {approx} 10{sup 5} K) gas in the Milky Way corona, based on the measurements of O VI, N V, C IV, Si IV, and Fe III absorption lines seen in the far-ultraviolet spectra of 58 sight lines to extragalactic targets, obtained with the Far-Ultraviolet Spectroscopic Explorer and the Space Telescope Imaging Spectrograph. In many sight lines the Galactic absorption profiles show multiple components, which are analyzed separately. We find that the highly ionized atoms are distributed irregularly in a layer with a scale height of about 3 kpc, which rotates along with the gas in the disk, without an obvious gradient in the rotation velocity away from the Galactic plane. Within this layer the gas has randomly oriented velocities with a dispersion of 40-60 km s{sup -1}. On average the integrated column densities are log N(O VI) = 14.3, log N(N V) = 13.5, log N(C IV) = 14.2, log N(Si IV) = 13.6, and log N(Fe III) = 14.2, with a dispersion of just 0.2 dex in each case. In sight lines around the Galactic center and Galactic north pole, all column densities are enhanced by a factor {approx}2, while at intermediate latitudes in the southern sky there is a deficit in N(O VI) of about a factor of two, but no deficit for the other ions. We compare the column densities and ionic ratios to a series of theoretical predictions: collisional ionization equilibrium, shock ionization, conductive interfaces, turbulent mixing, thick disk supernovae, static non-equilibrium ionization (NIE) radiative cooling, and an NIE radiative cooling model in which the gas flows through the cooling zone. None of these models can fully reproduce the data, but it is clear that NIE radiative cooling is important in generating the transition temperature gas.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

INFRARED AND X-RAY SPECTROSCOPY OF THE Kes 75 SUPERNOVA REMNANT SHELL: CHARACTERIZING THE DUST AND GAS PROPERTIES

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of {approx}1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keV. Previous studies suggest that the hot component may originate from reverse-shocked supernova (SN) ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from SN ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of {approx}140 K by a relatively dense, hot plasma that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 1.3 Multiplication-Sign 10{sup -2} M{sub Sun }, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2022

Influence of Large-Scale Interplanetary Structures on the Propagation of Solar Energetic Particles: the Multi-Spacecraft Event on 2021 October 9

An intense solar energetic particle (SEP) event was observed on 2021 October 9 by multiple spacecraft distributed near the ecliptic plane at heliocentric radial distances R ≲ 1 au and within a narrow range of heliolongitudes. A stream interaction region (SIR), sequentially observed by Parker Solar Probe (PSP) at R = 0.76 au and 48° east from Earth (ϕ = E48°), STEREO-A (at R = 0.96 au, ϕ = E39°), Solar Orbiter (SolO; at R = 0.68 au, ϕ = E15°), BepiColombo (at R = 0.33 au, ϕ = W02°), and near-Earth spacecraft, regulated the observed intensity-time profiles and the anisotropic character of the SEP event. PSP, STEREO-A, and SolO detected strong anisotropies at the onset of the SEP event, which resulted from the fact that PSP and STEREO-A were in the declining-speed region of the solar wind stream responsible for the SIR and from the passage of a steady magnetic field structure by SolO during the onset of the event. By contrast, the intensity-time profiles observed near Earth displayed a delayed onset at proton energies ≳13 MeV and an accumulation of ≲5 MeV protons between the SIR and the shock driven by the parent coronal mass ejection (CME). Even though BepiColombo, STEREO-A, and SolO were nominally connected to the same region of the Sun, the intensity-time profiles at BepiColombo resemble those observed near Earth, with the bulk of low-energy ions also confined between the SIR and the CME-driven shock. This event exemplifies the impact that intervening large-scale interplanetary structures, such as corotating SIRs, have in shaping the properties of SEP events.

D. Lario↗

NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

Titan’s atmosphere is composed mostly of N 2 with a small amount of CH 4 , and so, shock layers around craft entering Titan’s atmosphere will contain a variety of molecules formed from H, C, and N atoms, including the cyanogen radical CN. Sensitivity analysis has shown that the radiative heat flux predicted by computational fluid dynamics (CFD) simulations of Titan entry has up to 14% uncertainty due to the rate coefficients for collisional (de)excitation reactions that control the population of CN in its first and second excited states. The red and violet emission bands from CN’s first and second excited states, respectively, are known to be large sources of radiative heat flux on capsules entering Titan’s atmosphere.[2, 3] So, the simulated population of CN in its first and second excited states is very important, but currently has some inherent uncertainty coming from the data for the rate coefficients for reaction 1. The goal of the present project is to provide improved rate coefficient data for these reactions from first principles quantum chemistry calculations. This work reports on preliminary electronic structure calculations generated at a large number of triatomic geometries of interest, which show multiple avoided crossings at collinear arrangements. This suggests that collisional (de)excitation of CN by N atoms is likely to proceed through these geometries.

Eric Geistfeld↗

DOE OSTI2012

ZERO IMPACT PARAMETER WHITE DWARF COLLISIONS IN FLASH

We systematically explore zero impact parameter collisions of white dwarfs (WDs) with the Eulerian adaptive grid code FLASH for 0.64 + 0.64 M {sub Sun} and 0.81 + 0.81 M {sub Sun} mass pairings. Our models span a range of effective linear spatial resolutions from 5.2 Multiplication-Sign 10{sup 7} to 1.2 Multiplication-Sign 10{sup 7} cm. However, even the highest resolution models do not quite achieve strict numerical convergence, due to the challenge of properly resolving small-scale burning and energy transport. The lack of strict numerical convergence from these idealized configurations suggests that quantitative predictions of the ejected elemental abundances that are generated by binary WD collision and merger simulations should be viewed with caution. Nevertheless, the convergence trends do allow some patterns to be discerned. We find that the 0.64 + 0.64 M {sub Sun} head-on collision model produces 0.32 M {sub Sun} of {sup 56}Ni and 0.38 M {sub Sun} of {sup 28}Si, while the 0.81 + 0.81 M {sub Sun} head-on collision model produces 0.39 M {sub Sun} of {sup 56}Ni and 0.55 M {sub Sun} of {sup 28}Si at the highest spatial resolutions. Both mass pairings produce {approx}0.2 M {sub Sun} of unburned {sup 12}C+{sup 16}O. We also find the 0.64 + 0.64 M {sub Sun} head-on collision begins carbon burning in the central region of the stalled shock between the two WDs, while the more energetic 0.81 + 0.81 M {sub Sun} head-on collision raises the initial post-shock temperature enough to burn the entire stalled shock region to nuclear statistical equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Laser absorption measurements of temperature, pressure, CO, and CO 2 at near-MHz rates in post-detonation fireballs with comparison to synthetic measurements

A laser absorption spectroscopy (LAS) diagnostic was used to obtain measurements of temperature, pressure, CO, and CO 2 at 500 kHz or 1 MHz in post-detonation fireballs produced by hemispherical samples of pentaerythritol tetranitrate (PETN). A quantum-cascade laser was scanned over multiple CO absorption transitions near 2008.5 cm −1 at 1 MHz, while an interband-cascade laser was scanned over a CO 2 absorption transition near 2394.8 cm −1 at 500 kHz. Light from each laser was combined onto a single path and passed through a detonation chamber approximately 83 mm above the 12-mm diameter hemispherical PETN charge. The CO and CO 2 absorption signals were post-processed to obtain time histories of temperature, pressure, species column pressures (P CO L, P CO2 L), and species column mole fractions (X CO L, X CO2 L). Additionally, schlieren imaging was performed simultaneously at 500 kHz to aid interpretation of the LAS measurements. Experimental and synthetic (i.e., CFD based) LAS measurements were compared to evaluate the accuracy of the CFD model and its ability to model the turbulent afterburning of the detonation products in air. In general, the experimental measurements exhibit reasonable agreement with the synthetic measurements at early times; thereby supporting the accuracy of the CFD model. Periods of disagreement between experimental and synthetic measurements at later times are most likely due to a reflected shock and detonator cavity jetting, which are not accounted for in the CFD model.

Schwartz, Charles J. [Purdue Univ., West Lafayette↗

NASA NTRS1995

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

DOE OSTI2025

Stability and convergence of nuclear detonations in white dwarf collisions

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗