TATB Thermal Decomposition: An Improved Kinetic Model for Explosive Safety Analysis
Explore the source record for details and available documents.
Materials science, Chemistry - Chemical explosives↗
SEARCH · Search NASA
Results for “Analytical mathematics”
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.
Electro-carburization of mild steel through two-step reaction in molten-salt electrolyte
The rapid carburization of mild steel is demonstrated via the electrochemical deposition of solid carbon, followed by thermochemical iron-carbon (Fe-C) reaction using a molten lithium-potassium carbonate (LiKCO 3 ) salt electrolyte and carbon dioxide as the carbon source. Depth-dependent microstructural and compositional changes in the steel are shown to be controlled by the applied cell voltage and reaction temperatures. Reaction temperatures above the Fe-C eutectoid temperature promote faster carburization due to carbon dissolution into austenite. The carbon deposition rate is regulated by the applied cell voltage: At - 2.4 V, carbon deposition occurs much faster than its diffusion into the steel, resulting in the accumulation of a carbon-rich layer exceeding 1 mm, which does not diffuse into the steel. In contrast, at - 1.8 V, all deposited carbon effectively diffuses into the steel, leaving a minimal external deposit. In both cases, carbon ingress forms a carbide-rich surface layer and precipitates additional carbides along the grain boundaries within the steel electrode. This method offers a potential approach to increase the carbon content of mild or low-carbon steel feedstocks, such as scrap or direct reduced iron, for industrial applications.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Newton-Okounkov bodies of chemical reaction systems
Despite their noted potential in polynomial-system solving, there are few concrete examples of Newton-Okounkov bodies arising from applications. Accordingly, in this paper, we introduce a new application of Newton-Okounkov body theory to the study of chemical reaction networks and compute several examples. An important invariant of a chemical reaction network is its maximum number of positive steady states Here, we introduce a new upper bound on this number, namely the ‘Newton-Okounkov body bound’ of a chemical reaction network. Further, through explicit examples, we show that the Newton-Okounkov body bound of a network gives a good upper bound on its maximum number of positive steady states. We also compare this Newton-Okounkov body bound to a related upper bound, namely the mixed volume of a chemical reaction network, and find that it often achieves better bounds.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore
Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.
08 HYDROGEN↗
Time-periodic steady-state solution of fluid-structure interaction and cardiac flow problems through multigrid-reduction-in-time
In this study, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved.
97 MATHEMATICS AND COMPUTING↗
Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism
A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Dynamic modeling and simulation of pressure swing adsorption processes using toPSAil
Pressure swing adsorption (PSA) has attracted significant recent interest for chemical process intensification due to its potential for high energy efficiency and amenability to small, modular designs. However, the lack of simulation tools that are readily available, transparent, and trusted, has been identified as a serious impediment to widespread adoption of PSA, as well as to further research on PSA modeling, numerical solution, optimization, and control. This paper presents a complete framework for dynamic modeling and simulation of PSA processes and its implementation in an open-source simulator called toPSAil. Further, the presentation is tutorial and includes many modeling and implementation details often overlooked in existing literature. Novel methods are presented for handling flow reversals and implementing various pressure–flow relationships, along with controlled boundary conditions. Finally, the code contains several innovations designed to improve efficiency and reduce the extensive trial-and-error tuning often required to produce a working PSA cycle.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A mixed integer linear programming approach for the design of chemical process families
Tackling climate change goals requires widespread deployment of process technology variants across many decentralized sites with different geographical, environmental, and operational requirements. Conventional engineering approaches focus on unique designs for each installation (process variant), while missing opportunities for manufacturing standardization. Here, instead we seek to optimize a process platform of common unit designs while simultaneously designing an entire family of process variants that make use of that platform. This reduces engineering effort, deployment timelines, and manufacturing costs. We propose a nonlinear generalized disjunctive programming formulation and convert this to an efficient mixed-integer linear programming (MILP) formulation through discretization of the design space. We formulate our optimization in Pyomo with costing from IDAES, and we demonstrate the computational performance and solution quality on a water treatment desalination system from the PARETO framework and a carbon capture system built in Aspen Plus as part of CCSI2.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DeePKS-kit: A package for developing machine learning-based chemically accurate energy and density functional models
We introduce DeePKS-kit, an open-source software package for developing machine learning based energy and density functional models. DeePKS-kit is interfaced with PyTorch, an open-source machine learning library, and PySCF, an ab initio computational chemistry program that provides simple and customized tools for developing quantum chemistry codes. It supports the DeePHF and DeePKS methods. In addition to explaining the details in the methodology and the software, we also provide an example of developing a chemically accurate model for water clusters.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Chemical-Free Lithium Separation from High-Salinity Brines Using Model-Informed and Machine Learning-Optimized Multi-Column Zwitterionic Chromatography
Direct Lithium Extraction (DLE) technologies often struggle to produce high-purity lithium salts from high-salinity brines, as current approaches require chemical-based elution, regeneration, and precipitation steps, resulting in significant environmental footprints. A novel salt fractionation approach using carboxybetaine resin, known as zwitterionic chromatography (ZIC), has demonstrated that lithium ions can be separated from divalent cations under high-salinity conditions using only water as eluent, with no regeneration required. To enable continuous and scalable deployment of this approach, we developed a chemical-free Multi-column Zwitterionic Chromatography (MZC) process and its theoretical and process models. To predict and optimize this nontraditional separation system, we introduced a novel anti-Langmuir isotherm, and the isotherm parameters were estimated through a machine learning-driven optimization based on artificial neural network ensembles with numerical feasibility assessment. Using machine learning-driven optimization, the MZC process achieved 98.0% lithium recovery, 99.5 % Li/(Li + Mg + Ca) purity, a 31.3% productivity increase, and a 33% reduction in water use compared to batch operation. The proposed MZC process enables lithium separation at $0.6-1.2 kg-1 Li, with costs dominated by resin manufacturing, while offering lower separation costs and carbon footprint compared with conventional carbonation. Overall, these findings position the MZC process as an effective polishing step within scalable and sustainable lithium production pipelines.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computational tools for the simulation and analysis of spin-polarized EPR spectra
The EPR spectra of paramagnetic species induced by photoexcitation typically exhibit enhanced absorptive and emissive features resulting from sublevel populations that differ from thermal equilibrium. The populations and the resulting spin polarization of the spectra are dictated by the selectivity of the photophysical process generating the observed state. Simulation of the spin-polarized EPR spectra is crucial in the characterization of both the dynamics of formation of the photoexcited state as well as its electronic and structural properties. EasySpin, the simulation toolbox for EPR spectroscopy, now includes extended support for the simulation of the EPR spectra of spin-polarized states of arbitrary spin multiplicity and formed by a variety of different mechanisms, including photoexcited triplet states populated by intersystem crossing, charge recombination or spin polarization transfer, spin-correlated radical pairs created by photoinduced electron transfer, triplet pairs formed by singlet fission and multiplet states arising from photoexcitation in systems containing chromophores and stable radicals. In this paper, we highlight EasySpin’s capabilities for the simulation of spin-polarized EPR spectra on the basis of illustrative examples from the literature in a variety of fields ranging across chemistry, biology, material science and quantum information science.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Physics informed machine learning for chemistry tabulation
Modeling of turbulent combustion system requires modeling the underlying chemistry and the turbulent transport. Solving both systems simultaneously is computationally prohibitive. Instead, given the difference in scales at which the two sub-systems evolve, the two sub-systems are typically (re)solved separately. Popular approaches such as the Flamelet Generated Manifolds (FGM) use a two-step strategy where the governing reaction kinetics are pre-computed and mapped to a low-dimensional manifold, characterized by a few reaction progress variables (model reduction) and the manifold is then “looked-up” during the run-time to estimate the high-dimensional system state by the turbulent transport system. While existing works have focused on these two steps independently, in this work we show that joint learning of the progress variables and the look-up model, can yield more accurate results. Here, we build on the base formulation and implementation to include the dynamically generated Thermochemical State Variables (Lower Dimensional Dynamic Source Terms). We discuss the challenges in the implementation of this deep neural network architecture and experimentally demonstrate its superior performance.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Update 1.1 to “pysimm: A python package for simulation of molecular systems”, (PII: S2352711016300395)
This update of the pysimm application programming interface, pysimm 1.1, provides both infrastructural as well as functional updates. Moreover, improvements to the random walk application that allow it to construct polymers with controlled tacticity are highlighted. Additions to the forcefield module include an update to enable working with the family of CHARMM forcefields and automated typing with the CHARMM generalized forcefield (CGenFF). Finally, new detailed examples demonstrating new features are also provided.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Changing Directions: Influence of Ligand Electronics on the Directionality and Kinetics of Photoinduced Charge Transfer in Cu(I)Diimine Complexes
A key challenge to the effective utilization of solar energy is to promote efficient photoinduced charge transfer, specifically avoiding unproductive, circuitous electron-transfer pathways and optimizing the kinetics of charge separation and recombination. We hypothesize that one way to address this challenge is to develop a fundamental understanding of how to initiate and control directional photoinduced charge transfer, particularly for earth-abundant first-row transition-metal coordination complexes, which typically suffer from relatively short excited-state lifetimes. Here, in this paper, we report a series of functionalized heteroleptic copper(I)bis(phenanthroline) complexes, which have allowed us to investigate the directionality of intramolecular photoinduced metal-to-ligand charge transfer (MLCT) as a function of the substituent Hammett parameter. Ultrafast transient absorption suggests a complicated interplay of MLCT localization and solvent interaction with the Cu(II) center of the MLCT state. This work provides a set of design principles for directional charge transfer in earth-abundant complexes and can be used to efficiently design pathways for connecting the molecular modules to catalysts or electrodes and integration into systems for light-driven catalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Chemistry Graduate Student Experience: Findings from an ACS Survey
Graduate training is a key element in producing a scientific workforce that reflects the nation’s diversity. This paper examines data from a 2013 American Chemical Society (ACS) survey of 2,544 chemistry masters and doctoral students and reveals barriers to reaching this goal. Multivariate statistical analyses indicate that women reported significantly less supportive relationships with advisors. Women were less likely to plan to finish their degrees, and for PhD students, the discrepancy was larger for students at the start of their graduate program. Women were also less likely to pursue the next level of training, and the gender difference related to postdoctoral plans was greater for those who identified with a racial-ethnic group traditionally underrepresented in chemistry (underrepresented minority, URM). URM students who were beyond the first year of their graduate program reported significantly less supportive relationships with peers. They were also less likely to have funding sufficient to meet their needs and more often used personal resources including loans. Despite these difficulties, URM students were more likely to definitely plan to finish their degrees, and men who identified as URM were more likely to plan to pursue postdoctoral work. Independent of gender and identification as URMs, students in more highly ranked schools reported less advisor support. Extensive open-ended comments indicated that large proportions of the students desired more attention and meaningful feedback from advisors and changes within their programs to promote support for students and advisor accountability. Suggestions for future research are given, and a companion commentary discusses needed directions for change.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Protein Surface Printer for Exploring Protein Domains
The surface of proteins is vital in determining protein functions. Herein, a program, Protein Surface Printer (PSP), is built that performs multiple functions in quantifying protein surface domains. Two proteins, PETase and cytochrome P450, are used to validate that the program supports atomistic simulations with different combinations of programs and force fields. A case study is conducted on the structural analysis of the spike proteins of SARS-CoV-2 and SARS-CoV and the human cell receptor ACE2. Although the surface domains of both spike proteins are highly similar, their receptor-binding domains (RBDs) and the O-linked glycan domains are structurally different. The O-linked glycan domain of SARS-CoV-2 is highly positively charged, which may promote binding to negatively charged human cells.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Predicting Energetics Materials’ Crystalline Density from Chemical Structure by Machine Learning
To expedite new molecular compound development, a long-sought goal within the chemistry community has been to predict molecules’ bulk properties of interest a priori to synthesis from a chemical structure alone. In this work, we demonstrate that machine learning methods can indeed be used to directly learn the relationship between chemical structures and bulk crystalline properties of molecules, even in the absence of any crystal structure information or quantum mechanical calculations. We focus specifically on a class of organic compounds categorized as energetic materials called high explosives (HE) and predicting their crystalline density. An ongoing challenge within the chemistry machine learning community is deciding how best to featurize molecules as inputs into machine learning models—whether expert handcrafted features or learned molecular representations via graph-based neural network models—yield better results and why. We evaluate both types of representations in combination with a number of machine learning models to predict the crystalline densities of HE-like molecules curated from the Cambridge Structural Database, and we report the performance and pros and cons of our methods. Our message passing neural network (MPNN) based models with learned molecular representations generally perform best, outperforming current state-of-the-art methods at predicting crystalline density and performing well even when testing on a data set not representative of the training data. However, these models are traditionally considered black boxes and less easily interpretable. Here, to address this common challenge, we also provide a comparison analysis between our MPNN-based model and models with fixed feature representations that provides insights as to what features are learned by the MPNN to accurately predict density.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗