Search NASA⌕ Search

SEARCH · Search NASA

Results for “burning velocity”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗
DOE OSTI2021

LARGE PARTICLE SEPARATION FROM NON-NEWTONIAN SLURRIES USING BUMP ARRAYS

Here we evaluate the performance of bump arrays to separate large particles from non-Newtonian slurries with Bingham and Cross rheology. Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit. Although these devices have been used to separate a wide variety of large particles from blood cells to sand, partition of large particles from non-Newtonian fluids remains unexplored. Yet, an important set of modestly concentrated slurries, including Hanford nuclear waste, displays non-Newtonian rheology. Here we evaluate the influence of non-Newtonian rheology on the large-small particle size cutoff in bump arrays using a model that explores the influence of yield stresses, ratios of zero and infinite shear viscosities, and Cross’s exponent under strictly laminar well-developed conditions. Surprisingly, we find that viscosity ratios and Cross’s exponent make no significant difference on the particle cutoffs between large particles that bump and small particles that slalom around the posts from entrance to exit. In contrast, we find that yield stresses do significantly affect the size cutoff. As the yield stress increases, the velocity profiles become more plug like lowering the size cutoff. For nuclear waste separations where removing large particles is a priority, increasing yield stresses is conservative.

slurry, solids removal, deterministic later displa↗
DOE OSTI2023

Coupling hydrodynamics to a rigid-body motion solver for fluid-structure interaction [Slides]

Quinoa is a massively parallel computational fluid dynamics with multi-material and programmed burn capabilities developed from Programmatic and LDRD funds. Overset is a well-established method of using an “inset” mesh that communicates with a background mesh. Solutions transfer freely from one to another and operate as boundary conditions on the opposite mesh. Either mesh can move at specified velocity! This project couples the existing Mesh-to-mesh transfer with Quinoa and demonstrates use of resulting Overset solver in computational simulation of blast effects on a re-entry body for survivability assessments.

42 ENGINEERING↗
DOE OSTI2023

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗
DOE OSTI2024

Flame Flashback Investigations in Hydrogen-Enriched Low Swirl Flames using High-Speed Hydroxyl (OH) Planar Laser-Induced Fluorescence

This study reports flashback mechanism observed in hydrogen-enriched flames stabilized in a low swirl burner (LSB) at atmospheric pressure and temperature conditions. The fundamentals of hydrogen-rich stable flames and spatiotemporal investigation of flashback phenomena were observed experimentally using a high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) diagnostic. Testing was conducted in an optically accessible pre-mixing section of LSB with inlet pre-mixing velocities from 5 to 10 m/s for methane and hydrogen (50 – 90% by mole) blends. Swirlers with two different turning angles, 26 and 33 degrees and three different perforated plate hole diameters, 1.08, 1.12 and 1.16 mm with measured swirl numbers varying from 0.43 to 0.49 were used. The flashback propensity showed dependence on the proximity of the lifted flame to the burner exit which was dependent on pre-mixer velocity, hydrogen content and equivalence ratio (ER). High-speed OH-PLIF images revealed lifted flames were first observed at low ER condition which is changed to M-shaped flame attached to burner rim with the increase in ER. Then, further increase in ER triggers flame flashback into the nozzle. Flashback occurred when burning occurs in the outer shear layer and the leading flame brush ingresses into the nozzle.

Parajuli, Pradeep↗
DOE OSTI2024

Birefringence tests of gravity with multimessenger binaries

Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. Here, we introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary sources, for which the initial polarization is determined by the inclination and orientation angles of the orbital angular momentum vector with respect to the line of sight. In particular, we use spatially resolved radio imaging of the jet from a binary neutron star (BNS) merger to constrain the orientation angle and hence the emitted polarization orientation of the GW signal at the site of the merger, and compare to that observed on Earth by GW detectors. For GW170817, using past measurements of the inclination angle, we constrain the deviation from GR due to amplitude birefringence to κ A = $-0.12^{+0.60}_ {-0.61}$, while the velocity birefringence parameter κ V remains unconstrained. The inability to constrain κ V is due to the low amplitude of GW170817 in the Virgo detector, and measurements of the polarization orientation require information from a combination of multiple detectors with different alignments. For this reason, we also mock future BNS mergers with resolved afterglow proper motion and project that κ V could be constrained to a precision of 5 rad (corresponding to an angular shift of the GW polarization of δ φV ≈ 0.2 rad for a BNS at 100 Mpc) by a future network of third-generation ground-based GW detectors such as Cosmic Explorer and the radio High Sensitivity Array. Crucially, this velocity birefringence effect cannot be constrained with dark binary mergers as it requires polarization information at the emission time, which can be provided only by electromagnetic emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Rotation Rate for the Planetary-mass Companion DH Tau b

DH Tau b is a young planetary-mass companion orbiting at a projected separation of 320 au from its ∼2 Myr old host star DH Tau. With an estimated mass of 8–22 M{sub Jup} this object straddles the deuterium-burning limit, and might have formed via core or pebble accretion, disk instability, or molecular cloud fragmentation. To shed light on the formation history of DH Tau b, we obtain the first measurement of rotational line broadening for this object using high-resolution (R ∼ 25,000) near-infrared spectroscopy from Keck/NIRSPEC. We measure a projected rotational velocity (v sin i) of 9.6 ± 0.7 km s{sup −1}, corresponding to a rotation rate that is between 9% and 15% of DH Tau b’s predicted breakup speed. This low rotation rate is in good agreement with scenarios in which magnetic coupling between the companion and its circumplanetary disk during the late stages of accretion reduces angular momentum and regulates spin. We compare the rotation rate of DH Tau b to published values for other planetary-mass objects with masses between 0.3 and 20 M{sub Jup} and find no evidence of a correlation between mass and rotation rate in this mass regime. Finally, we search for evidence of individual molecules in DH Tau b’s spectrum and find that it is dominated by CO and H{sub 2}O, with no evidence of the presence of CH{sub 4}. This agrees with expectations given DH Tau b’s relatively high effective temperature (∼2300 K).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An indirect approach to optimize the reaction rates of thermal NO formation for diesel engines

With stringent emission regulations, it has become more important for modern diesel engine manufacturers to accurately predict engine-out nitrogen oxide (NO x ) emissions across a wide range of operating conditions. Thermal NO is the major source of engine-out NO x in modern diesel engines. For thermal NO formation, several earlier studies have recommended the forward and reverse reaction rate coefficients of the rate-limiting reaction (O + N 2 ⇌ NO + N). However, due to deficiencies in sub-models and inadequacies of reduced chemical mechanisms to represent diesel combustion, these recommended values more often than not need to be adjusted in reduced order combustion models to accurately predict engine-out NO x . Hence, in this work a systematic and computationally efficient approach has been proposed to streamline the process of determining the optimum reaction rate coefficients. Here, to develop the optimization approach, four different production diesel engines with different operating conditions in terms of speed, load, and exhaust-gas recirculation have been considered. Numerical simulations have been performed using a detailed zero-dimensional velocity-composition-frequency transported probability density function (0D-VCF-tPDF) model that uses hundreds of notional particles to capture in-cylinder stratification. Four different combinations of hydrocarbon and NO x chemical mechanisms were used to represent chemistry. It was found that for the rate-limiting reaction, the pre-exponent factors (A f1 , A r1 ) and activation energies (E A,f1 , E A,r1 ) of the forward and reverse reaction rates follow a linear band in A f1 - E A,f1 and A r1 - E A,r1 space where predicted engine-out NO x match the measured values closely. By encompassing such bands from different engines and considering constraints on activation energies, a reduced search domain of pre-exponent factors and activation energies was constructed that is expected to be applicable to any diesel engine. Eventually, computationally efficient three-line and one-line search approaches were proposed to determine the optimum values of the pre-exponent factors and activation energies that led to a minimum error between measured and predicted engine-out NO x . Finally, these three-line and one-line NO x optimization approaches were applied to a fifth production diesel engine for which the 0D-VCF-tPDF model showed a very good predictive performance in terms of predicting peak pressure, 50% burn rate, and engine-out NO x when compared to measured and 3D-CFD values.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dissolution of Fluoride Salts in Hanford Tank Waste

The Direct Feed High-Level Waste (DFHLW) strategy seeks to bypass the Hanford Waste Treatment and Immobilization Plant Pretreatment Facility while retaining some processing functions to maximize waste feed loading and minimize high-level waste (HLW) waste volume. The DFHLW flowsheet needs leaching, washing, and solids concentration operations either in new or existing tanks. The effectiveness and efficiency of sludge washing has a substantial impact on DST space, mission duration, and the evaporation and low-activity waste (LAW) treatment operations required by these large wash-water additions. One target species requiring washing is fluoride. The HLW glass composition limits for fluorine drive operations to dissolve fluoride-bearing salts into the LAW fraction and thereby maximize waste loading in HLW glass. The fluoride in many high-level wastes at Hanford is predominantly in the form of fluoride-salt precipitates: villiaumite (NaF), kogarkoite (Na 3 FSO 4 ), and natrophosphate (Na 7 F(PO 4 ) 2 ·19H 2 O). Fluoride produces melter off-gas that creates corrosion risk in the off-gas system piping, while the sulfate and phosphate in the fluoride double salts kogarkoite and natrophosphate can be detrimental to glass waste loading. Fluoride salts are sparingly soluble, with solubilities ranging from approximately 40 to 130 kg per kL of pure water, and the dissolution kinetics of the three fluoride salts are not well known. Unexpected delays in a tank dissolution process could be encountered as a result of the lack of information about dissolution rates. In addition, if the double salts show transient non-stoichiometric dissolution of fluoride versus phosphate or sulfate, unexpectedly high concentrations of one of these other constituents could be produced. Washington River Protection Solutions authorized Pacific Northwest National Laboratory to collect the available data for fluoride salt dissolution rate, provide a scoping estimate of dissolution time if possible, and identify gaps in the understanding and predictive capability for estimating dissolution time. Open literature and Hanford reports were reviewed to document, understand, and (where possible) evaluate limitations on fluoride salt equilibria and dissolution kinetics, including both mass transport and surface reaction rate. Scoping estimates of dissolution time were made for mass-transfer-controlled dissolution of spherical particles of fluoride salts suspended in liquid. This is not the only potential governing mechanism; dissolution could be substantially slower if the surface reaction rate (the rate of release of ions from the surface) is the controlling mechanism. When the minimum amount of water for complete dissolution is used and the slip velocity between the liquid and suspended particles is less than or equal to the terminal settling velocity, the estimated mass-transfer rates allow 0.1 mm particles of fluoride salt to dissolve in minutes at 25 °C in water containing no other dissolved salts. Much larger solids, such as the 6-mm chunks that have been seen in heels, could take a few hours to more than a week to dissolve. The actual dissolution times will depend strongly on the actual slip velocity, the extent of particle suspension, constraint by surface reaction rates, the ratio of solvent to solid, and the presence of common ions that shift the solubility equilibria to restrict dissolution of fluoride salts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Development of a Direct Neutral Density Diagnostic for Fusion Edge Plasmas

As magnetically confined plasmas progress towards ignition and long pulse experiments, measurement and control of the neutral density in the plasma edge has become a critical issue for stability, formation of transport barriers, and fueling. Recent experiments by our research group have demonstrated that is possible to use two-photon absorption laser induced fluorescence (TA-LIF) to directly measure the density of neutral hydrogen in helicon sources and in the HIT-SI3 spheromak. While those experiments validated key elements of a diagnostic system that would enable similar measurements in tokamak plasmas, they did not fully address the issue of performance optimization in the presence of intense background light at the fluorescence wavelength and they also identified a number of other issues that must be resolved before successful TALIF neutral density measurements in a tokamak are likely to be achieved. Additional specific concerns raised during design reviews with the leadership of the DIII-D tokamak facility included: the minimum detection threshold for this TALIF diagnostic (currently ~ 5 x 10 15 m -3 ) and the rate at which neutral density measurements could be obtained. Improving the system performance to reduce the minimum detection threshold and demonstrating a faster rate of density determination are required to advance this diagnostic to the level where it could be considered for implementation on a major tokamak facility in the USA. The key issues that will be addressed through additional technological development are: validation of a new, chromatic aberration-free, xenon calibration scheme; increasing the output power of the laser; and suppression of background emission at the fluorescence wavelength through optimization of the collection optics and gating of the photomultiplier detector with a goal of obtaining a minimum detection threshold of 5 x 10 14 m -3 . Another key technological development, so far only tested in helicon source experiments, is the validation of Doppler-free TALIF as a means of obtaining higher speed, calibrated, neutral density measurements in tokamak-like conditions at the full cadence of the pulsed TALIF laser. In this work, we propose to complete these technological advancements through installation and testing of our prototype TALIF system on the proto-MPEX experiment at Oak Ridge National Laboratory (ORNL). The proto-MPEX facility will provide plasma conditions similar to the edge plasma of a major tokamak experiment but with pulse lengths and pulse repetition rates ideally suited for extensive development of a TALIF neutral density diagnostic. Measurement of the absolute neutral density in the edge of a magnetically confined plasma is necessary for plasma density control; calculation of charge-exchange power losses; control of plasma-wall interactions; determination of the braking of plasma flow; determination of the fuel mixture in deuterium-tritium plasma; and understanding the dynamics of the divertor region in the plasma edge. In terms of potential impact on the worldwide fusion program, we note that in burning plasma experiments, such as ITER, a 100 mJ/pulse Doppler-free TALIF diagnostic should be capable of directly measuring the deuterium/tritium (D/T) isotope ratio in the outer 0.5 to 1.0 m of the plasma radius. The D/T ratio is a critically important parameter for the control and optimization of burning plasmas. Therefore, development of the diagnostic system proposed here is also relevant to long-term US participation in diagnostic development for large tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Converting existing optical detectors into fast x-ray detectors

Here, the very short burn time and small size of burning plasmas created at advanced laser-fusion facilities will require high-spatial-resolution imaging diagnostics with fast time resolution. These instruments will need to function in an environment of extremely large neutron fluxes that will cause conventional diagnostics to fail because of radiation damage and induced background levels. One solution to this challenge is to perform an ultrafast conversion of the x-ray signals into the optical regime before the neutrons are able to reach the detector and then to relay image the signal out of the chamber and into a shielded bunker, protected from the effects of these neutrons. With this goal in mind, the OMEGA laser was used to demonstrate high-temporal-resolution x-ray imaging by using an x-ray snout to image an imploding backlighter capsule onto a semiconductor. The semiconductor was simultaneously probed with the existing velocity interferometry system for any surface reflector (VISAR) diagnostic, which uses an optical streak camera and provided a one-dimensional image of the phase in the semiconductor as a function of time. The phase induced in the semiconductor was linearly proportional to the x-ray emission from the backlighter capsule. This approach would then allow a sacrificial semiconductor to be attached at the end of an optical train with the VISAR and optical streak camera placed in a shielded bunker to operate in a high neutron environment and obtain time-dependent one-dimensional x-ray images or time-dependent x-ray spectra from a burning plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Search for non-virialized axions with 3.3-4.2 $\mu$eV mass at selected resolving powers

The Axion Dark Matter eXperiment is sensitive to narrow axion flows, given axions compose a fraction of the dark matter with a non-negligible local density. Detecting these low-velocity dispersion flows requires a high spectral resolution and careful attention to the expected signal modulation due to Earth's motion. We report an exclusion on the local axion dark matter density in narrow flows of $\rho_a \gtrsim 0.03\,\mathrm{GeV/cm^3}$ and $\rho_a \gtrsim 0.004\,\mathrm{GeV/cm^3}$ for Dine-Fischler-Srednicki-Zhitnitski and Kim-Shifman-Vainshtein-Zakharov axion-photon couplings, respectively, over the mass range $3.3-4.2\,\mu\text{eV}$. Measurements were made at selected resolving powers to allow for a range of possible velocity dispersions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Violent mergers revisited: The origin of the fastest stars in the Galaxy

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M ⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 −3 M ⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Astronomy and AstroPhysics↗
DOE OSTI2020

Combustion of metal powder with dinitrogen tetroxide

In this work, we present analysis of a novel reactive material system that employs dinitrogen tetroxide (N 2 O 4 ) as a liquid oxidizer with metal powder fuels. The oxidizer was added to micron scale aluminum and zirconium powders by a remote injection system. When ignited with a high voltage spark, the mixtures were observed to possess reactivity comparable to nanocomposite reactive materials, with open-tube flame expansion velocities from 500 to 1400 m/s depending on fuel/oxidizer ratio and tube diameter. Temperatures were observed to range from 3000 to 3500 K as measured with gray body fits to 16-channel time-resolved pyrometer and time-integrated spectrometer data. These values were significantly below calculated adiabatic flame temperatures, which we attribute to local deviations from stoichiometry and kinetic/energetic limitations similar to those observed in studies of particles burning in high pressures gaseous oxidizers. Al/N 2 O 4 reactivity was found to be most likely limited by the vaporization of the metal from the particle surface and Zr/N 2 O 4 was limited by slower burning and complex interactions involving the solubility of nitrogen and oxygen in the molten Zr. We also discuss the potential for these materials to be used to create an “on/off” reactive material, since N2O4 can be added remotely and driven to evaporate via vacuum or purge of inert gases to return it to a safe condition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of Heterogeneities on the Reaction-zone of a Propagating Detonation Wave

Reactive burn models for propagating detonation waves in an explosive are based on the ZND (Zel'dovich-Von Neumann-Doring) theory. A key property of the ZND theory, which leads to self-sustaining detonation waves, is a sonic point relative to the front in the reaction zone. Moreover, burn models assume the explosive is a homogeneous material. This assumption is used for plastic-bonded explosives, even though they are heterogeneous materials. The heterogeneities are accounted for only with an empirical burn rate. Furthermore, the reaction-zone width can be less than the length scale of heterogeneities. This raises the question of the effect of heterogeneities on the reaction zone. In particular, whether the reaction zone can be steady and have a well defined sonic point. To examine this question we focus on PBX 9501 since there is data on the grain-binder heterogeneities and on the reaction-zone velocity time history of a steady planar detonation wave. Using pseudo-mesoscale 1-D simulations, we examine the effect of a detonation wave transversing explosive and binder segments of random length. The simulations show that the reaction zone is quasi-steady with only short wavelength small amplitude variations in the detonation wave pressure and speed. The main effect of the heterogeneities is to introduce noise in the pressure and velocity fields that propagates into the Taylor wave behind the detonation. This is consistent with the scatter in the reaction-zone measurements for PBX 9501.

42 ENGINEERING↗