TATB Thermal Decomposition: An Improved Kinetic Model for Explosive Safety Analysis
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Various small-scale experiments were performed to provide data for developing a model to predict the thermal response of LLM-105 over a wide range of conditions. The thermal decomposition of LLM-105 was studied as a function of sample mass, confinement of volatile products, and preheating time in both isothermal and ramped heating experiments. The thermal decomposition of LLM-105 is a two-step process, as shown by the two exothermic peaks in the heat flow profiles, which were fitted to two nth-order autocatalytic reaction models with a similar activation energy of ∼289 kJ/mol. The magnitude and shape of these peaks varied with sample mass and confinement. Increasing sample mass enhanced the second exotherm with respect to the first one, while increasing the level of confinement promoted a transition from a sublimation-dominated regime towards thermal decomposition. The effect of LLM-105 particle size on the rate of weight loss was evident for open-pan experiments, where bigger particles sublimed at lower temperatures than smaller particles. Thermal response and solid residue composition of LLM-105 samples were analyzed following preheating for different durations. Longer preheating times caused a shift of the second exotherm to lower temperatures and a decrease in the reaction enthalpy, confirming that LLM-105 decay is a consecutive reaction mechanism, probably autocatalytic. In conclusion, the kinetic model derived from ramped experiments was validated against the measured LLM-105 fraction remaining and the enthalpy remaining of the solid residue as a function of preheating times and showed good agreement.
Emergence of Enterobacteriaceae harboring metallo-beta-lactamases (MBL) has raised global threats due to their broad antibiotic resistance profiles and the lack of effective inhibitors against them. We have been studied one of the emerging environmental MBL, the L1 from Stenotrophomonas maltophilia K279a . We determined several crystal structures of L1 complexes with three different classes of beta-lactam antibiotics (penicillin G, moxalactam, meropenem, and imipenem), with the inhibitor captopril and different metal ions (Zn +2 , Cd +2 , and Cu +2 ). All hydrolyzed antibiotics and the inhibitor were found binding to two Zn +2 ions mainly through the opened lactam ring and some hydrophobic interactions with the binding pocket atoms. Without a metal ion, the active site is very similarly maintained as that of the native form with two Zn +2 ions, however, the protein does not bind the substrate moxalactam. When two Zn +2 ions were replaced with other metal ions, the same di-metal scaffold was maintained and the added moxalactam was found hydrolyzed in the active site. Differential scanning fluorimetry and isothermal titration calorimetry were used to study thermodynamic properties of L1 MBL compared with New Deli Metallo-beta-lactamase-1 (NDM-1). Both enzymes are significantly stabilized by Zn +2 and other divalent metals but showed different dependency. These studies also suggest that moxalactam and its hydrolyzed form may bind and dissociate with different kinetic modes with or without Zn +2 for each of L1 and NDM-1. Our analysis implicates metal ions, in forming a distinct di-metal scaffold, which is central to the enzyme stability, promiscuous substrate binding and versatile catalytic activity. Statement: The L1 metallo-beta-lactamase from an environmental multidrug-resistant opportunistic pathogen Stenotrophomonas maltophilia K279a has been studied by determining 3D structures of L1 enzyme in the complexes with several beta-lactam antibiotics and different divalent metals and characterizing its biochemical and ligand binding properties. We found that the two-metal center in the active site is critical in the enzymatic process including antibiotics recognition and binding, which explains the enzyme's activity toward diverse antibiotic substrates. This study provides the critical information for understanding the ligand recognition and for advanced drug development.
Abstract Silica‐supported amine absorbents, including materials produced by tethering aminosilanes or infusion of poly(ethyleneimine), represent a promising class of materials for CO 2 capture applications, including direct air and point source capture. Various silica surface treatments and functionalization strategies are explored to enhance stability and CO 2 uptake in amine‐based solid sorbent systems. Here, the synthesis and characterization of novel vinyltrimethoxysilane‐treated Santa Barbara Amorphous‐15 (SBA‐15) supports and the corresponding enhancement in CO 2 uptake compared to various SBA‐15‐based control supports are presented. The relationship between CO 2 diffusion and amine efficiency in these systems is explored using a previously reported kinetic model. The synthesized materials are characterized with CO 2 and H 2 O isotherms, diffuse reflectance infrared Fourier transform spectroscopy, 1 H T 1 – T 2 relaxation correlation NMR, and rapid thermal cycling experiments. The novel support materials are shown to enable high amine efficiencies, approaching a fourfold improvement over standard SBA‐15‐supported amines, while simultaneously exhibiting excellent stability when cycled rapidly under humid conditions. As the poly(ethyleneimine) loadings are held constant across the various samples, enhancements in CO 2 uptake are attributed to differences in the way the poly(ethyleneimine) interacts with the support surface.
Double soaking (DS) is a thermal processing route intended to produce austenite–martensite microstructures in steels containing austenite‐stabilizing additions and consists of intercritical annealing (primary soaking), followed by heating and brief isothermal holding at an increased temperature (secondary soaking), and quenching. Herein, experimental dilatometry during DS of a medium‐manganese (Mn) steel with nominally 7 wt% Mn and an ultralow residual carbon concentration, in combination with phase‐field simulations of austenite formation during secondary soaking, is presented. The feasibility of maintaining heterogeneous Mn distributions during DS is demonstrated and insight is provided on the effects of the secondary soaking temperature and prior Mn distribution on the ferrite‐to‐austenite phase transformation during the secondary soaking portion of the DS treatment.
Abstract To study the thermal stability of ethylene/vinyl acetate/vinyl alcohol (EVA‐OH) composite with filler particles, systematic thermogravimetric analysis (TGA) investigation was conducted in both nonisothermal and isothermal modes. The effect of polymer concentration on the aging behavior of EVA‐OH was investigated in the nitroplasticizer (NP) environment. Fourier transform infrared spectroscopy was used to probe chemical structural changes in the EVA‐OH polymer before and after thermal aging. The results suggest that filler addition accelerates the rate of NP exudation out of the EVA‐OH composites and deteriorates the thermal stability of NP at moderate temperatures (up to 70°C). The degradation of NP, in turn, accelerates the degradation of EVA‐OH polymer in its composite form through oxidation and hydrolysis indicating the importance of antioxidants in such phase blends.
Abstract The crystal structure and bonding environment of K 2 Ca(CO 3 ) 2 bütschliite were probed under isothermal compression via Raman spectroscopy to 95 GPa and single crystal and powder X-ray diffraction to 12 and 68 GPa, respectively. A second order Birch-Murnaghan equation of state fit to the X-ray data yields a bulk modulus, $${K}_{0}=46.9$$ K 0 = 46.9 GPa with an imposed value of $${K}_{0}^{\prime}= 4$$ K 0 ′ = 4 for the ambient pressure phase. Compression of bütschliite is highly anisotropic, with contraction along the c -axis accounting for most of the volume change. Bütschliite undergoes a phase transition to a monoclinic C 2/ m structure at around 6 GPa, mirroring polymorphism within isostructural borates. A fit to the compression data of the monoclinic phase yields $${V}_{0}=322.2$$ V 0 = 322.2 Å 3 $$,$$ , $${K}_{0}=24.8$$ K 0 = 24.8 GPa and $${K}_{0}^{\prime}=4.0$$ K 0 ′ = 4.0 using a third order fit; the ability to access different compression mechanisms gives rise to a more compressible material than the low-pressure phase. In particular, compression of the C 2/ m phase involves interlayer displacement and twisting of the [CO 3 ] units, and an increase in coordination number of the K + ion. Three more phase transitions, at ~ 28, 34, and 37 GPa occur based on the Raman spectra and powder diffraction data: these give rise to new [CO 3 ] bonding environments within the structure.
In a recent publication, we carried out an extensive analysis of an unsupervised method of determining optimum shift factors in time-temperature-superposition of accelerated-aging data that involves minimizing the vertical arclength to obtain the master curve. For synthetic Arrhenius data with a variety of noise distributions, the work showed that, in conjunction with bootstrap-resampling, the method can produce reliable estimates of the mean activation energy along with uncertainty quantification. Here, we apply the above method to six different datasets taken from the published literature and demonstrate accurate prediction of mean activation energy from the data as-is without the need for any pre-processing or fitting. We also compare uncertainty margins computed by second-order bootstrap with that by linear regression theory and show that the former appears to provide consistent margins in the presence of common noise types in real data, including intra-isotherm measurement-errors, sample-to-sample variations, and intrinsic deviation from perfect Arrhenius behavior.
There is high market demand for increasing the viability of current plastic recycling processes. In this work rheology is used to evaluate the mechanical properties of a solvent dissolution recycled polymer compared to its virgin untreated precursor. Solvent dissolution and precipitation are used to target multi-layer, multi-component, industry films, which cannot be mechanically recycled. Polyethylene was chosen as the primary polymer of interest. Polymer thermal stability was monitored via time-resolved rheology; consecutive frequency sweeps over the course of an hour while under isothermal conditions. Additional rheological experiments were performed within the identified thermally stable conditions. Small-angle oscillatory shear was complemented with steady shear viscosity experiments over a wide range of shear rates. Extensional rheology was used to determine changes in molecular weight and cross link density. As a result, rheological characterization is supplemented with gas chromatography–mass spectrometry of the solvent wash to determine components stripped from the virgin polymers during solvent treatment.
Post synthetic modification of metal organic frameworks presents a viable route for amine functionalization which can significantly enhance CO 2 sorption capacity. We present a facile means of amine incorporation using limited synthetic steps and low-cost reagents that results in a high density of primary alkyl amines distributed through a zirconium-based metal organic framework (MOF). Both the MOF synthesis and the post synthetic modification take place under aqueous conditions and result in strongly bound molecular amines available for sorbate interaction throughout the MOF pores. Furthermore, this amine incorporation protocol results in a significantly increased CO 2 capacity compared with the unmodified MOF-808. Specifically, CO 2 isotherms collected at 298 K for the unmodified MOF-808 show uptakes of 0.06, 0.2, and 1.2 mmol/g at 4, 15, and 100 kPa, respectively, which can be compared with 0.3, 0.7, and 2.5 mmol/g for the glycine grafted MOF-808 and 0.5, 0.9, and 2.3 mmol/g for ethylenediamine grafted MOF-808.
Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.
Small-volume, high-throughput screening techniques are sought to enable downselection from a large candidate pool of bio-blendstocks to a select few, having physical properties consistent with requirements of downsized, turbo-boosted internal combustion engines. Herein, this work presents a droplet evaporation-based approach to predict heat of vaporization, vapor pressure, diffusion coefficient, and Lennard–Jones parameters for an unknown fuel. Two different schemes, considering the isothermal evaporation of a moving droplet in ambient air, are proposed, which combine droplet velocity and temperature measurements, with some known properties to predict unknown properties. The schemes utilize an inverse solution of a transient model of droplet evaporation solved in an iterative fashion. A baseline scheme, which only requires droplet size change measurements, is evaluated using test data for three liquid fuels, comprising of alkanes and alcohols, as obtained in a temperature-controlled chamber. Results yield temperature-dependent heat of vaporization and vapor pressure predictions within 10 % and 22 %, respectively, of reference values. The advanced scheme, which additionally requires droplet temperature measurement, is numerically evaluated in the current work and will be experimentally validated in future efforts. The advanced scheme is found to significantly improve prediction quality, with deviations less than 2 % and 1 % for heat of vaporization and vapor pressure, while also predicting diffusion coefficient and Lennard–Jones parameters within 5 % and 8 %, respectively. The combined set of approaches, which primarily track droplet evaporation, can be incorporated into a small-volume, high-throughput fuel screening process.
Uncertainty in specific heat capacity values of a molten salt determined by using differential scanning calorimetry (DSC) was assessed based on the precision of replicate measurements of heat flows used in the calculation and effects of corrections that are commonly made to heat flow measurements. The ratio method of determining heat capacity was applied using the results of replicate heat flow measurements made with two empty cells, a sapphire reference material, and three samples of a doped NaCl-UCl 3 salt mixture. Replicate measurements with empty cells were used to quantify the effects of system instabilities and sensitivities on the measured heat flows of sapphire and salt. The combined effects of uncertainties in individual heat flow measurements made with blank cells using this system were quantified to be 2.6 μV based on isothermal holds before and after the scan, with cell placement adding the greatest uncertainty. This value was used as the tolerance for accepting background-corrected heat flows measured with sapphire and salt to calculate the specific heat capacity. The acceptable heat flows measured for sapphire and salt over the temperature range of 540 to 725 °C resulted in calculated specific heat capacity values ranging from 0.53 to 0.91 J g −1 K −1 with an overall average value of 0.70 J g −1 K −1 and an uncertainty of 0.22 J g −1 K −1 at the 99 % confidence level. The combined uncertainty in the specific heat capacity masked detection of any effect of temperature or salt composition that occurred.
Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.
Here, isothermal mutual diffusion coefficients (interdiffusion coefficients) were measured for K 2 SO 4 (aq) at 298.15 ± 0.005 K, at numerous concentrations ranging from dilute solutions to near saturation (0.59648 mol∙dm -3 ; 0.61349 mol∙kg -1 ) under free diffusion boundary conditions, using high precision Rayleigh interferometry. Under the experimental conditions these diffusion coefficients are on the volume-fixed reference frame D v . Two series of experiments were performed, the first using the traditional experimental approach with a mercury lamp source and with the interference patterns being recorded on glass photographic plates, and the second with a He-Ne laser light source and computerized data acquisition using a photodiode array. The results from both series of experiments are in excellent agreement, and generally yield diffusion coefficients precise to at least 0.002 x 10 -9 m 2 ∙s 1 (0.15% to 0.19%) and in most cases to 0.001 x 10 -9 m 2 ∙s -1 (0.07% to 0.1%). These experiments also yield accurate values of the refractive index differences for the solution pairs used in the diffusion experiments. The new diffusion coefficients are compared to two sets of published values of diffusion coefficients for K 2 SO 4 (aq) which are somewhat discrepant from each other. This study extends and complements our earlier work on the diffusion coefficients of the most common brine salts: NaCl(aq), KCl(aq), MgCl 2 (aq), CaCl 2 (aq), Na 2 SO 4 (aq), ] MgSO 4 (aq), NaHCO 3 (aq), and KHCO 3 (aq) at 298.15 K.
Here, in this study, we investigated the volatility of sodium, potassium, and cesium perrhenates using thermogravimetry (TG) performed on single salts and their binary mixed salts. We performed simultaneous TG and differential scanning calorimetry under a ramp heating condition at 10°C/min up to 1200°C to characterize volatilization behavior as a function of temperature. We also performed isothermal TG over a temperature range of 700 to 975°C to estimate the vapor pressures of the alkali perrhenates based on a modified Langmuir equation that relates the rate of free evaporation from the open surface to vapor pressure. In this paper, we discuss the results of our study on the temperature of evaporation and temperature-dependent vapor pressure of sodium, potassium, and cesium perrhenates in relation to rhenium volatilization behavior observed during melting of simulated waste glass.
In recent years, high entropy alloys (HEAs), also known as multi-principal component alloys (MPCAs) have emerged as a new and exciting class of materials. This paper reports on the solid state reduction synthesis of a series of CoFeNiMn-based MPCA compositions, starting from a mixture of the corresponding oxides. One of the aims of the study was to test whether the degree of reduction of MnO, a highly stable oxide, could be enhanced by tailoring the alloy composition. Specifically, the influence of Ni content was studied because Ni exhibits a significant negative enthalpy of mixing with Mn. High purity precursor powders of Co(OH) 2 , Fe 2 O 3 , MnO 2 , and NiO were milled and mixed using standard ceramic processing methods. Here, the nominal sample compositions (assuming complete oxide reduction) were (CoFeMn) x Ni (1-x) , for x = 0, 0.083, 0.166, and 0.25. The oxide samples were subjected to a series of isothermal reduction anneals in flowing 3 pct H 2 –Ar at 1100 °C. The resulting microstructures were characterized using scanning electron microscopy (SEM), X-ray energy dispersive spectroscopy (EDS) and X-ray diffraction (XRD). The composition of the resulting MPCAs was determined quantitatively using wavelength dispersive spectroscopy (WDS) in the electron microprobe. The study revealed that for each of the initial oxide compositions studied, it was possible to achieve an MPCA with ~ 25 at. pct Mn. These results were found to be consistent with the predictions of a thermodynamic model whereby a negative enthalpy of mixing (ΔH mix ), combined with a contribution from configurational entropy, can offset a positive free energy of reduction (ΔG red ). The incorporation of vibrational entropy into first principles calculations was found to have a significant effect on the predicted crystal structure of the MPCAs.
NbTi-based refractory alloys with additions of Co, Cr, and Ni represent an interesting medium-entropy alloy system with potential for protective oxide film formation, high strength, and ductility. This study investigates the microstructural evolution, oxidation behavior, and mechanical properties of NbTi-based alloys containing 5 to 10 at. pct Co, Cr, and Ni. CALPHAD predictions suggest that this composition range can be heat treated to obtain a predominantly body-centered cubic matrix phase. Mechanical properties, including microhardness, yield strength, maximum strength, and specific strength are evaluated through isothermal compression tests conducted between room temperature and 800 °C. The oxidation kinetics of these alloys are assessed through discontinuous oxidation tests. Parabolic oxidation kinetics were observed for NbTi–10Ni and NbTi–5Co, while linear oxidation kinetics were found for NbTi–10Cr and NbTi–10(CoCrNi). Microstructures and oxide layers are characterized using X-ray diffraction, electron backscatter diffraction, energy-dispersive X-ray spectroscopy, and scanning electron microscopy. All alloys exhibit significant mechanical softening between room temperature and 800 °C, with elastic-perfectly plastic flow observed at 800 °C. The addition of 10 pct Cr to NbTi resulted in two BCC phases up to 1050 °C, conflicting with CALPHAD predictions of a single-phase solid solution at this temperature, and resulting in higher flow stress at 800 °C. NbTi–10(CoCrNi) exhibited the lowest flow stress at 800 °C despite having more ‘cocktail effect’ potential and insufficient molar fractions of Co, Cr, or Ni to form a desirable protective oxide film.