Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex mixtures”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of silver($\mathrm{I}$) ion reduction on the selectivity of olefins, paraffins, and aromatic compounds by gas chromatographic stationary phases consisting of silver($\mathrm{I}$) salts in ionic liquids

We report the olefin/paraffin selectivity offered by ionic liquid (IL) stationary phases can be enhanced through the addition of silver(I) ion, which is well-known to undergo selective complexation with unsaturated compounds. However, such stationary phases often suffer from the loss of chromatographic selectivity as silver(I) ion can be reduced to elemental silver. To maintain the separation performance of silver(I) ion/IL stationary phases, an understanding of factors and conditions that promote the reduction of silver(I) ion is needed. In this study, capillary gas chromatography columns featuring a stationary phase consisting of the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 - ]) IL impregnated with [Ag + ][NTf 2 - ] were examined to investigate the effects of temperature, hydrogen content in exposure gas stream, and time of heating/exposure events on olefin selectivity. Retention factors of representative analytes, such as C 6 olefins and paraffins as well as aromatic compounds, were measured after subjecting the columns to the aforementioned conditions, followed by an evaluation of selectivity factors over time. Selectivity factors of olefins and aromatic compounds were observed to decrease significantly when the stationary phases were heated to temperatures higher than 110°C as well as being subjected to mixed gas streams containing greater than 50 mol% of hydrogen. As constant column heating temperatures were applied under exposure gas mixtures containing hydrogen and nitrogen, a gradual decrease in analyte selectivity factors was observed under prolonged periods of time. However, application of a ternary gas mixture comprised of 25/50/25 mol% hydrogen/nitrogen/methane resulted in an increase in the 3-hexyne/cis-2-hexene selectivity when measured at 120°C for 60 h, due to a smaller decrease in the retention factor of 3-hexyne compared to cis-2-hexene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Plasma-based global pathway analysis to understand the chemical kinetics of plasma-assisted combustion and fuel reforming

The Global Pathway Analysis (GPA) algorithm helps analyze the chemical kinetics of complex combustion systems by identifying important global reaction pathways connecting a source species to a sink species through various important intermediate species (i.e., hub species). Here, the present work aims to extend GPA algorithm to plasma-assisted combustion and fuel reforming systems to identify the dominant global pathways in such systems at various conditions. In addition, the present study extends the ability of GPA algorithm to identify reaction cycles involving the excitation of high-concentration species (e.g., O 2 , N 2 , and fuel) to their vibrational and electronic states and the subsequent de excitation to their ground state, based on their significance on the reactivity of plasma-assisted systems in terms of gas heating and radical production. Provisions are made in the GPA algorithm to evaluate the reactivity of identified re action pathways and cycles based on the element-flux transfer (i.e., dominance), heat release, and radical production rate. The newly developed Plasma-based Global Pathway Analysis (PGPA) algorithm is then used to analyze the plasma assisted combustion of ammonia and reforming of methane. The PGPA analyses elucidated the significance of vibrational-translational cycles on the reactivity of NH 3 /air mixtures. Further, analyses on the production of NO ascribed the early reforming of NH 3 to N 2 and H 2 in impeding the production of NO during plasma-assisted NH 3 ignition. Lastly, the enhanced reforming of CH 4 /N 2 mixtures using plasma has been attributed to electron impact dissociation of CH 4 when compared to thermal reforming. In contrast, conventional path-Flux analysis (PFA) was found to require significant manual effort and pre-analysis intuitions from expert knowledge, making it arduous to provide valuable in sights into plasma chemistry. The user-friendly and automated nature of PGPA thus provides a valuable tool for assessing the kinetics of plasma-assisted systems helpful in analyzing and, further, a foundation in reducing plasma-assisted chemistry, without the needs of expert knowledge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Tritium Breeding Ratio Evaluation of Solid Breeder Concepts for the FESS-FNSF

This paper presents a parametric study of the Fusion Energy System Studies-Fusion Nuclear Science Facility’s (FNSF’s) tritium breeding performance for several solid breeder concepts, neutron multiplying materials, and blanket materials, assuming volume fractions based on the most recent FNSF design as a realistically representative fusion facility. In this study, we initially surveyed the tritium breeding ratio (TBR) of several solid breeder concepts by employing a simplified but efficient one-dimensional (1-D) infinite cylinder reduced-order model (ROM). Parametric studies were performed with the ROMs for the full range of breeder-to-multiplier ratios to identify the optimum mixture compositions for each breeder type that would lead to a maximum TBR. These optimized breeder-multiplier combinations were then homogenized with FNSF blanket component materials to estimate their impacts on the TBR. Subsequently, as a validation step for the optimal designs, TBR calculations were performed using a more realistic modified 1-D ROM with inner and outer breeding regions, as well as with a fully detailed 22.5-deg three-dimensional (3-D) sector of the FNSF to assess the impact of geometry details on the TBR. The differences between the two 1-D models were negligible, while the ROMs were able to correctly predict trends and identify the maximum and minimum TBR cases, as well as show consistent biases relative to the results produced by the full 3-D, 22.5-deg sector for specific breeder/multiplier combinations. Solid breeder concepts such as Li 2 O, Li 4 SiO 4 , and Li 8 ZrO 6 outperformed all others in this study in terms of TBR performance when combined with all the neutron multiplier materials selected. Here, an underlying goal of this study was to develop and improve rapid and reliable ROMs to aid designers during parametric optimizations of highly complex and computationally expensive fusion models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗
DOE OSTI2024

Trailblazing Kr/Xe Separation: The Birth of the First Kr-Selective Material

Efficient separation of Kr from Kr/Xe mixtures is pivotal in nuclear waste management and dark matter research. Thus far, scientists have encountered a formidable challenge: the absence of a material with the ability to selectively adsorb Kr over Xe at room temperature. This study presents a groundbreaking transformation of the renowned metal–organic framework (MOF) CuBTC, previously acknowledged for its Xe adsorption affinity, into an unparalleled Kr-selective adsorbent. This achievement stems from an innovative densification approach involving systematic compression of the MOF, where the crystal size, interparticle interaction, defects, and evacuation conditions are synergistically modulated. The resultant densified CuBTC phase exhibits exceptional mechanical resilience, radiation tolerance, and notably an unprecedented selectivity for Kr over Xe at room temperature. Simulation and experimental kinetic diffusion studies confirm reduced gas diffusion in the densified MOF, attributed to its small pore window and minimal interparticle voids. The lighter Kr element demonstrates facile surface passage and higher diffusivity within the material, while the heavier Xe encounters increased difficulty entering the material and lower diffusivity. This Kr-selective MOF not only represents a significant breakthrough in Kr separation but also demonstrates remarkable processability and scalability to kilogram levels. Furthermore, the findings presented herein underscore the transformative potential of engineered MOFs in addressing complex challenges, heralding a new era of Kr separation technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deconvoluting Structures of Component Plant Biopolymers Using Deuterium Labeled Brassica oleracea Stems

Several Brassica species are cultivated globally for the production of seed oils as food, lubricants, and increasingly biofuels. The stem and leaf residues of these herbaceous dicotyledonous crops constitute another feedstock for biofuels and other products. Plant cell walls are complex, multipolymeric structures that consist primarily of polysaccharides and lignin. Cellulose chains coalesce to form crystalline microfibrils while the amorphous biopolymers, hemicellulose, and lignin form a network structure and fill the interstitial space. Neutron scattering has been used for the structural study of the assembly and deconstruction of the plant cell walls. However, similar neutron sensitivity of the different amorphous biopolymers has made structural association to individual component biopolymers nontrivial and ambiguous. To improve the association of structural features to specific biopolymer components, partial deuteration of the plant cell wall can be employed to increase the difference in the neutron scattering length density between amorphous carbohydrate and lignin plant polymers. Vegetative stems from partially deuterated Brassica oleracea acephala (kale) plants were obtained commercially, and the plant cell wall structures were studied by contrast-variation small-angle neutron scattering (CV-SANS) and Fourier-transform infrared spectroscopy (FTIR). FTIR results indicated that deuterium substitution for hydrogen in the carbohydrates was higher than that in lignin. By combining CV-SANS and FTIR results, the neutron scattering length density (nSLD) of the polysaccharides and lignin was determined to match nSLD of 65%:35% and 48%:52% D 2 O/H 2 O solvent mixtures, respectively. These nSLD values were higher than the nSLD values of polysaccharides and lignin for H 2 O grown biopolymers. The nSLD increase correlates to replacing about 42.5% of hydrogens present as both C–H and O–H groups in cellulose with deuterium atoms, while only 21% for lignin. In conclusion, this study lays the foundation to use partially deuterated plants to deconvolute structural features of the different component biopolymers, especially cellulose, hemicellulose, and lignin, of the plant cell wall without introducing unintended structural modifications due to the pretreatment extraction processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Crystallographic Snapshots of Pre- and Post-Lanthanide Halide Hydrolysis─Reaction Products Captured by the 4-Amino-1,2,4-triazole Ligand

Reactions of lanthanide(III) chloride salts with 4-amino-1,2,4-triazole (4-NH 2 -1,2,4-Triaz) in azole melts have led to the isolation of both hydrolysis and non-hydrolysis products in the same synthesis with the inclusion of a variety of ligands, anions, and water, allowing us to capture crystallographic snapshots of different forms and intermediate hydrolysis fragments. The structural studies reported here include anhydrous and hydrated nonhydrolyzed complexes which were isolated alongside hydrolysis products giving oxide/hydroxide lanthanide(III) dimers, tetramers, and ultimately hexamers. The compounds isolated include [Nd 2 Cl 6 (μ 2 -4-NH 2 -1,2,4-Triaz) 4 (4-NH 2 -1,2,4-Triaz) 2 ], [Ce 2 Cl 4 (μ 2 -Cl) 2 (μ 2 -4-NH 2 -1,2,4-Triaz) 4 ] n , [Ce 2 (µ 2 -Cl) 4 (µ 2 -OH) 2 (µ 2 -4-NH 2 -1,2,4-Triaz) 2 ] n , [Ln 4 Cl 4 (µ 2 -Cl) 4 (µ 3 -OH) 4 (µ 2 -4-NH 2 -1,2,4-Triaz) 4 ] n •2nH 2 O (Ln = Ce, Nd), and [Ce 6 Cl 6 (µ 6 -O 0.5 )(µ 3 -Cl 0.5 ) 4 (µ 3 -Cl 0.75 ) 3 (µ 3 -OH) 0.75 (µ 2 -4-NH 2 -1,2,4-Triaz) 12 ((OH 2 ) 0.25 ) 2 ] 2 [CeCl 6 ][Cl 9 ]•xH 2 O. In all complexes all lanthanide atoms are pairwise connected via one or more 4-NH2-1,2,4-Triaz ligands and sometimes additional Cl - anions.

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

Reducing Antisolvent Use in the STRAP Process by Enabling a Temperature–Controlled Polymer Dissolution and Precipitation for the Recycling of Multilayer Plastic Films

The recently reported processing strategy called solvent-targeted recovery and precipitation (STRAP) enables deconstruction of multilayer plastic packaging films into their constituent resins by selective dissolution. It uses a series of solvent washes that are guided by thermodynamic calculations of polymer solubility. In this work, the use of antisolvents in the STRAP process was reduced and solvent mixtures were considered to enable the temperature-controlled dissolution and precipitation of the target polymers in multilayer films. This was considered as a means to further improve the STRAP process and its estimated costs. Two STRAP approaches were compared based on different polymer precipitation techniques: precipitation by the addition of an antisolvent (STRAP-A) and precipitation by decreasing the solvent temperature (STRAP-B). Both approaches were able to separate the constituent polymers in a post-industrial film composed primarily of polyethylene (PE), ethylene vinyl alcohol (EVOH), and polyethylene terephthalate (PET) with near 100 % material efficiency. Technoeconomic analysis indicates that the minimum selling price (MSP) of the recycled resins with STRAP-B is 21.0 % lower than that achieved with STRAP-A. This provides evidence that thermally driven polymer precipitation is an option to reduce the use of antisolvents, making the STRAP process more economically and environmentally attractive. A third process, STRAP-C, was demonstrated with another post-industrial multilayer film of a different composition. Here, the results demonstrate that this process can also recover polymers at similar costs to those of virgin resins, indicating that the STRAP technology is flexible and can remain economically competitive as the plastic feed complexity is increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Probing the Framework Metal Dependent Properties of Actinide-Centered Polyoxoalkoxide Sandwich-Type Complexes

Development of a simple and scalable synthesis of (TBA) 3 [W 5 O 18 MoNO] provides for the formation of the mixed-metal lacunary polyoxoalkoxide, (TBA) 2 [W 4 O 13 (OMe) 4 MoNO][Na(MeOH)]. This complex was used to synthesize a series of polyoxoalkoxide sandwich-type complexes with the general formula (TBA) 2 [M{W 4 O 13 (OMe) 4 MoNO} 2 ], where M = Zr(IV), Hf(IV), Th(IV), U(IV), and Np(IV). Compared to the analogous all-molybdenum complexes, the series have drastically different optical and redox properties. The results indicate that framework metal substitution acts as a tool for “orbital engineering”, with Density Functional Theory (DFT) calculations revealing that the major consequence of incorporation of tungsten into the complexes is localization of LUMO and LUMO+1 on the molybdenum centers remaining in the molecule. The change in the distribution of the frontier orbitals translates to discrepancies in the electronic properties of the series. Given the rarity of polyoxometalate complexes featuring a U(V) ion, one electron oxidation of (TBA) 2 [U(IV){W 4 O 13 (OMe) 4 MoNO} 2 ] was pursued. Isolation of the corresponding U(V) centered sandwich-type complex is reported, only the second example of U(V)-polyoxometalate complex described to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Probing the Effects of NO x and SO x Impurities on Oxy-Fuel Combustion in Supercritical CO 2 : Shock Tube Experiments and Chemical Kinetic Modeling

The direct-fired supercritical carbon dioxide cycles are one of the most promising power generation methods in terms of their efficiency and environmental friendliness. Two important challenges in implementing these cycles are the high pressure (300 bar) and high CO 2 dilution (>80%) in the combustor. The design and development of supercritical oxy-combustors for natural gas require accurate reaction kinetic models to predict the combustion outcomes. The presence of a small amount of impurities in natural gas and other feed streams to oxy-combustors makes these predictions even more complex. During oxy-combustion, trace amounts of nitrogen present in the oxidizer is converted to NO x and gets into the combustion chamber along with the recirculated CO 2 . Similarly, natural gas can contain a trace amount of ammonia and sulfurous impurities that get converted to NO x and SO x and get back into the combustion chamber with recirculated CO 2 . In this work, a reaction model is developed for predicting the effect of impurities such as NO x and SO x on supercritical methane combustion. The base mechanism used in this work is GRI Mech 3.0. H 2 S combustion chemistry is obtained from Bongartz et al. while NOx chemistry is from Konnov. The reaction model is then optimized for a pressure range of 30–300 bar using high-pressure shock tube data from the literature. It is then validated with data obtained from the literature for methane combustion, H 2 S oxidation, and NO x effects on ignition delay. The effect of impurities on CH 4 combustion up to 16 atm is validated using NO x -doped methane studies obtained from the literature. In order to validate the model for high-pressure conditions, experiments are conducted at the UCF shock tube facility using natural gas identical mixtures with N 2 O as an impurity at ~100 bar. Current results show that there is a significant change in ignition delay with the presence of impurities. A comparison is made with experimental data using the developed model and predictions are found to be in good agreement. The model developed was used to study the effect of impurities on CO formation from sCO 2 combustors. It was found that NO x helps in reducing CO formation while the presence of H 2 S results in the formation of more CO. The reaction mechanism developed herein can also be used as a base mechanism to develop reduced mechanisms for use in CFD simulations.

Energy & Fuels↗
DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

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

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Perspectives on Quenching and Tempering 4340 Steel

Steels are ubiquitous due to their affordability and the landscape of useful properties that can be generated for engineering applications. But to further expand the performance envelope, one must be able to understand and control microstructure development by alloying and processing. In this paper we use multiscale, advanced characterization to better understand the structural and chemical evolution of AISI 4340 steel after quenching and tempering (Q&T), including the role of quench rate and short-time, isothermal tempering below 573 K (300 °C), with an emphasis on carbide formation. We compare the microstructure and/or property changes produced by conventional tempering to those produced by higher temperature, short-time “rapid” tempering. We underscore that no single characterization technique can fully capture the subtle microstructure changes like carbon redistribution, transition carbide and/or cementite formation, and retained austenite decomposition that occur during Q&T. Only the use of multiple techniques begins to unravel these complexities. After controlled fast or slow quenching, η transition carbides clearly exist in the microstructure, likely associated with autotempering of this high martensite start temperature (Ms) steel. Isothermal tempering below 598 K (325 °C) results in the relief of carbon supersaturation in the martensite, primarily by the formation of η transition carbides that exhibit a range of carbon levels, seemingly without substitutional element partitioning between the carbide and matrix phases. Hägg transition carbide is present between 300 °C and 325 °C. After conventional tempering at or above 598 K (325 °C) for 2 h, cementite is predominant, but small amounts of cementite are also present in other conditions, even after quenching. Previous work has indicated that silicon (Si) and substitutional elements partition between the cementite, which initially forms under paraequilibrium conditions, and the matrix. Phosphorous (P) may also be preferentially located at cementite/matrix interfaces after high temperature tempering. Slower quench rates result in greater amounts of retained austenite compared to those after fast quenching, which we attribute to increased austenite stability resulting from “autopartitioning”. Rapid, high temperature tempering is also found to diminish tempered martensite embrittlement (TME) believed to be associated with the extent of austenite decomposition, resulting in mechanical properties not attainable by conventional tempering, which may have important implications with respect to industrial heat treatment processes like induction tempering. Controlling the amount and stability of retained austenite is not only relevant to the properties of Q&T steels, but also next-generation advanced high strength steels (AHSS) with austenite/martensite mixtures.

36 MATERIALS SCIENCE↗