Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, Chemical”

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 19 records

Concerted cycloaddition reactions: a study using the parabolic model and quantum chemical modeling

Concerted cycloaddition reactions were studied by the method of intersecting parabolas (M3IP) and quantum chemical calculations. Experimental data were processed within the framework of the M3IP method and an algorithm for calculating the activation energies (E) and rate constants (k) for reactions from the enthalpies of reactions was developed. The parameters E and k for twelve cycloaddition reactions not studied previously were calculated. Factors affecting the activation energies were established and evaluated; these include the enthalpy of reaction, substituents, and the molecular structure of reactants. Quantum chemical modeling and topological analysis of transition states (TS) of six concerted cycloaddition reactions were performed. Depending on structure of the starting olefins, the TS of reactions can have either a symmetric or asymmetric geometry. This influences their electronic structures, the energies of chemical bonds, and the activation energies of reactions. A comparison of the activation energy values obtained from the M3IP and DFT(B3lyp/6-311++G** ) calculations revealed good agreement between them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical modelling of Delayed Ettringite Formation for assessment of affected concrete structures

Delayed Ettringite Formation (DEF) modelling is addressed in the context of structural analysis. First, a chemical model is specified to simulate the effects of the heating of concrete in terms of reactions involving sulfates, aluminates and alkalis. It can be decomposed into: dissolution of primary sulfoaluminates at high temperature, fixation of aluminates in hydrogarnet or carboaluminates, and formation of delayed ettringite at low temperature. The influences of alkalis, temperature and water saturation on these reactions are taken into account. The main result of this model is the amount of delayed ettringite. The parameters of the proposed chemical model are fitted on a large number of experiments taken from the literature. The chemical model is then coupled with a nonlinear mechanical model based on poro-mechanics modelling. An application of chemo-mechanical coupling shows the features and benefits of such modelling to help practitioners in the management of structures affected by DEF.

36 MATERIALS SCIENCE↗

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗

CFD-Chemical Model of One-Third Scale Demonstration of DOE Sealed Canister with ATR Fuel

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. A demonstration case for these DOE sealed canisters will be eventually performed with a one-third scale mockup that has been instrumented with thermocouples and gas concentration probes. This study seeks to model the planned canister for validation of the previously developed model, such that confidence in its long-term prediction can be increased. The deployment of the instrumented lid for online monitoring in intended for a period of 10+ years based on previous monitoring of commercial fuel storage; however, the model is still run for the previously used 50-year storage periods. This case is modeled with a G-value using a bi-linear function such that it decreases at higher dose rates to be consistent with experiments. Validation efforts of the model would likely revolve around the results for the 1st year, so the results of this timeframe are also highlighted. For dried fuel of the nominal decay heat (18W), the predicted hydrogen concentration is 0.19% after 1 year, 0.91% after 10 years, and 2.8% after 50 years, with a maximum pressure of 1.26 atm. Consistent with previous modeling, the decay heat of the fuel is the main factor that influences the results. For undried fuel in pure helium, the hydrogen concentration ranges from 0.23 to 1.2% after 1 year, 1.25-6.2% after 10 years and 8.8 to 18.84% after 50 years. For dried fuel in pure helium the hydrogen concentration ranges from 0.06% to 0.34% after 1 year, 0.35-2.11% after 10 years and 0.95 to 6.6% after 50 years. While long-term results in the presence of residual air are mostly the same, early reactions with O 2 can delay significant production of H 2 until it is consumed to form more water vapor, lowering the range to 0.14 to 1.0% after 1 year. The maximum absolute pressure that is reached across any scenario is 2.06 atm. In the event of residual air, the presence of nitric acid is possible in the range of 18-131 ppm after 1 year, 174-1180 ppm after 10 years, and 586-3500 ppm after 50 years. As long as the fuel is sufficiently dried, or of nominal decay heat, a 4% lower flammability limit of hydrogen will not be reached within a 10-year monitoring period.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Guide to CFD-Chemical Model for Spent Fuel Storage

This outlines the steps that were used to create the CFD model in Star-CCM+ to simulate the 50 year thermal history of ATR fuel stored in DOE sealed canisters. Then this goes through the steps to set up the processing the model results for the one-way coupled Cantera chemical model. The next part describes the process for setting up the 50-year Cantera model in Python. The goal of this is to document the procedure for quality assurance, as well as in the event another person will need to run the process to generate additional results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inflow, Outflow, Yields, and Stellar Population Mixing in Chemical Evolution Models

Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the SN Ia delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone chemical evolution code, we investigate the effects of and trade-offs between parameters. Two critical parameters are SFE and the outflow mass-loading parameter, which shift the knee in [O/Fe]–[Fe/H] and the equilibrium abundances that the simulations asymptotically approach, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]–[Fe/H] unlike the observed bimodality (separate high- α and low- α sequences) in this plane. A mix of one-zone models with inflow timescale and outflow mass-loading parameter variations, motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the two sequences better than a one-zone model with two infall epochs. We present [X/Fe]–[Fe/H] tracks for 20 elements assuming three different supernova yield models and find some significant discrepancies with solar neighborhood observations, especially for elements with strongly metallicity-dependent yields. We apply principal component abundance analysis to the simulations and existing data to reveal the main correlations among abundances and quantify their contributions to variation in abundance space. For the stellar population mixing scenario, the abundances of α -elements and elements with metallicity-dependent yields dominate the first and second principal components, respectively, and collectively explain 99% of the variance in the model. flexCE is a python package available at https://github.com/bretthandrews/flexCE.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum chemical modeling of possible reactions of mononuclear iron nitrosyl complex [Fe(SC(NH{sub 2}){sub 2}){sub 2}(NO){sub 2}]Cl•H{sub 2}O in an aqueous solution

Possible reactions of the mononuclear iron nitrosyl complex [Fe(SC(NH{sub 2}){sub 2}){sub 2}(NO){sub 2}]Сl•H{sub 2}O in an aqueous solution were studied using quantum chemical modeling. The electronic structures of the possible intermediates were examined by the NBO and AIM methods. The substitution of the thio ligand in the iron—sulfur nitrosyl complex by a water molecule and the formation of dimeric intermediate complexes were found to facilitate the NO-donation process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum-chemical modeling of possible reactions of Roussin’s red esters with aryl ligands in DMSO solution

Possible reactions of Roussin’s red esters [Fe{sub 2}(μ-SC{sub 6}H{sub 4}R){sub 2}(NO){sub 4}], where R = H, o-NH{sub 2}, m-NO{sub 2}, m-OH, or m-NH{sub 2}, in DMSO solution were investigated by quantum chemical modeling. It was shown that in these systems, numerous reactions can occur, including NO donation, ligand substitution, and decomposition into mononuclear iron nitrosyl complexes. The resulting compounds are also NO donors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical Modeling of ATR Fuel in DOE Standard Canisters with Borated Stainless Steel Corrosion

Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium backfilled sealed DOE standard canisters. The packaging configuration for the DOE canisters places 10 ATR elements within a Type 1a basket, and three baskets are loaded within each DOE canister. In order to provide a criticality control mechanism, these stainless baskets are proposed to be built with 1-2 weight percentage of boron. This report seeks to identify effects of corrosion of this basket material over time. The primary reaction found through surveying literature was the oxidation of the iron in steel by water to produce some hydrogen. This reaction was built into the previous chemical model of the DOE canister to assess additional hydrogen and pressure build up over time. Literature for non-borated stainless steel, and carbon steel has shown that irradiation increases the rate of corrosion, and it is assumed a similar effect would occur with borated stainless steel. Due to the wide range of results in literature, and the variation of the corrosion with boron content, a wide range of corrosion rates was tested in the model from 80 nm/yr up to 1800 nm/yr.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗

Quantum Chemical Modeling of Pressure-Induced Spin Crossover in Octahedral Metal-Ligand Complexes

Spin state switching on external stimuli is a phenomenon with wide applicability, ranging from molecular electronics to gas activation in nanoporous frameworks. Here, we model the spin crossover as a function of the hydrostatic pressure in octahedrally coordinated transition metal centers by applying a field of effective nuclear forces that compress the molecule towards its centroid. For spin crossover in first-row transition metals coordinated by hydrogen, nitrogen, and carbon monoxide, we find the pressure required for spin transition to be a function of the ligand position in the spectrochemical sequence. While pressures on the order of 1 GPa are required to flip spins in homogeneously ligated octahedral sites, we demonstrate a fivefold decrease in spin transition pressure for the archetypal strong field ligand carbon monoxide in octahedrally coordinated Fe 2+ in [Fe(II)(NH 3 ) 5 CO] 2+ .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗

H–D exchange between quercetin and solvent in the presence of Au{sup I} chloride complexes with DMSO: quantum chemical modeling

The pathways of H–D exchange between quercetin and solvent taking into account the effects of the medium are considered using quantum chemical calculations and the mechanism of selective H–D exchange reaction at the C(2′) position observed only in the presence of Au complexes is identified. It is shown that only the exchange processes in the A-ring have suffi- ciently low activation barriers. In the acidic and alkaline media, these processes occur at lower energies as compared with the reaction in the neutral system, which qualitatively agrees with the published experimental data. It was shown that the Au{sup I} chloride complex with DMSO is able to react with the C–H bond at the C(2′) position by the mechanism of electrophilic substitution under mild conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum chemical modeling of superbase-catalyzed reactions of acetophenone and methyl mesityl ketone with acetylene

The mechanisms of furan formation and the competitive reaction of 1,5-diketone formation from 1-mesitylbut-2-en-1-one (1a) and 1-phenylbut-2-en-1-one (1b) were studied using the B3LYP/6-311++G**//B3LYP/6-31+G* quantum chemical approach taking into account solvation effects in terms of the IEFPCM model. The addition to 1a of the ethynide ion to form β-acetylene enolate ion (2a) occurs with the activation energy ΔG{sup ≠} = 13.3 kcal mol{sup –1}, while the further intramolecular O-vinylation with the closure of 2a to the dihydrofuran cycle of the dihydrofuran carbanion is the rate-determining step in the reaction pathway of substituted furan formation needing the activation energy G{sup ≠} = 19.5 kcal mol{sup –1}. The alternative addition of the mesitylethenolate ion to 1a to form the 1,5-diketone carbanion is associated with a higher (by 3.4 kcal mol{sup –1}) free activation energy than that of ethynide ion addition. In contrast, in the case of 1b, the addition at the C=C bond of the phenylethenolate ion rather than of the ethynide ion to form the 1,5-diketone carbanion is more preferable, since its activation free energy is 1.5 kcal mol{sup –1} lower. The quantum chemical calculation results are consistent with the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical study of a micro flow reactor at engine pressures: Flames with repetitive extinction and ignition and simulations with a reduced chemical model

Combustion in a narrow channel with an imposed temperature gradient is studied numerically at elevated pressure with engine-relevant fuels. Here, the focus is placed on unsteady flames with repetitive extinction and ignition (FREI) to determine the potential of this regime for fuel testing and calibration of reduced chemical mechanisms. First, it is shown that the FREI regime does occur at elevated pressures for sufficiently small tube diameters. The sensitivity of the extinction and ignition temperatures to low-temperature chemistry is found to be significantly enhanced at 25 bar compared to atmospheric conditions. The ignition and extinction temperatures differ by up to 100 K between PRF mixtures with varying octane numbers. Ternary mixtures of iso-octane/n-heptane/toluene and iso-octane/n-heptane/ethanol at similar research and motor octane numbers are also studied. Second, the potential of using data from the micro flow reactor to infer reaction rates is assessed. A reduced chemical mechanism combining a small fuel-dependent submechanism with a detailed fuel-independent submechanism for the core species chemistry is used for that purpose. Only the most sensitive fuel-dependent reactions are inferred. the calibrated reduced model is compared to a detailed model and good agreement in ignition delay times and laminar flame speeds is observed. Finally, this illustrates the potential of a micro flow reactor with a controlled temperature profile to test fuels and infer kinetic data.

33 ADVANCED PROPULSION SYSTEMS↗

Code Modifications for Modeling Chemical Tracers and Embedded Natural Fractures at EGS Collab

The EGS Collab SIGMA-V project is a multi-lab and university collaborative research project that is being undertaken at the Sanford Underground Research Facility (SURF) in South Dakota. The project consists of studying stimulation, fluid-flow, and heat transfer processes at a scale of 10-20 m, which is readily amenable to detailed characterization and monitoring. One objective of the project is to establish circulation from injector to producer by hydraulically fracturing the injector. Data generated during these experiments is to be compared with predictions from coupled thermal, hydrological, mechanical, and chemical simulators. One such a simulator, TOUGH2-CSM, has been enhanced in order to simulate EGS Collab SIGMA-V project experiments. These modifications include adding tracers, the capability to model tracer sorption, and an embedded fracture formulation. A set of example problems validate our conservative tracer transport and sorption formulations. We then simulated tracer transport and thermal breakthrough for the first EGS Collab SIGMA-V experiment.

15 GEOTHERMAL ENERGY↗

Application of Process Chemical Modeling to Optimize Radioactive Waste Disposal at the Savannah River Site - 24242

The Technical Optimization Model (TOM) is used by Savannah River Mission Completion (SRMC) to carry out facility-wide material balance and validate chemistrydependent processes for the purpose of their System Plan.The TOM simulates material movement and chemical reactions at the Tank Farm (TF), Salt Waste Processing Facility (SWPF) and Defense Waste Processing Facility (DWPF).

Georgiou, Andreas↗