Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetics”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Experimental study on kinetic oxidation of graphite IG-110 by steam

Graphite is proposed for use in High-temperature Gas-cooled Reactors (HTGRs) as the fuel matrix, neutron moderator/reflector, and core structural material. One important property of nuclear grade graphite is their resistance to oxidation in high-temperature environment. Extensive investigation has been performed in the literature for graphite oxidation by air. However, available experimental data are still limited for graphite oxidation by steam under conditions comparable to a postulated steam ingress accident in HTGRs. In this study, the oxidation rate of graphite IG-110 by steam was measured at temperatures from 850 to 1100 °C with the steam partial pressure varying from 0.5 to 20.0 kPa and the hydrogen partial pressure varying from 0 to 2.0 kPa. Further analysis confirms the oxidation process in this present study is dominated by the chemical kinetics, which lends credit to the data for being used to develop numerical models. It was observed that the increase of the kinetic oxidation rate with the steam partial pressure tends to become less apparent if the steam partial pressure keeps increasing. In addition, it was found that the partitioning of hydrogen inhibits the graphite-steam reaction process even with the steam partial pressure up to 20.0 kPa. However, this inhibiting effect starts to become saturated when the hydrogen partial pressure exceeds 1.0 kPa. The oxidation rates were fitted to the conventional Langmuir-Hinshelwood (LH) and Boltzmann-enhanced Langmuir-Hinshelwood (BLH) models by a multivariable optimization algorithm. The BLH model exhibits a better accuracy than the LH model within the specified experimental conditions. The predicted oxidation rate using the BLH model shows a mean relative difference of about 24% with the maximum difference of about 55% when compared with our experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Kinetics of photogenerated carbon dangling bonds in organic photovoltaic thin Films: An EPR study

Here, we report an investigation of the early kinetics of photogenerated carbon dangling bond (CDB) formation and annealing in organic photovoltaic bulk heterojunction (BHJ) thin film blends under oxygen- and moisture-free conditions, using X-band electron paramagnetic resonance (EPR) spectroscopy. The study focuses on donor:acceptor BHJ blends of PCE12:PCBM and PCE12:ITIC films, where PCE12 is PBDB-T. The time evolution of CDBs in such drop-cast BHJ films irradiated at 300 nm is monitored. The early kinetics of CDB formation, critical for understanding OPV degradation mechanisms, is studied. Theoretical analysis of the defect growth mechanism suggests a monomolecular defect creation model where the defect count follows a power-law t β with irradiation time t, where β ∼ 0.55–0.58, in excellent agreement with the theoretically expected value of β = 1/2. This model is compatible with CDB formation by the holes in donor sites adjacent to acceptors, likely assisted by energy released from quenching of nearby excitons by the holes, elucidating the physical mechanism underlying CDB formation. This is significant for designing improved materials, which mitigate defect creation, and consequently advancing the development of stable OPV systems.

42 ENGINEERING↗
DOE OSTI2020

Thermodynamic and Kinetic Roles of H2 in Structure Evolution of Urchin-like Co: A density functional theory study

Small gas molecules working as good capping agents play important roles in controlling the morphologies and surface structures of the metal nanocrystals. In present work, the thermodynamic and kinetic roles of H2 molecule played on the morphology of Co nanocrystal were systematically studied based on density functional theory (DFT). The Gibbs surface free energies (?) of Co(100), Co(110) and Co(111) under different hydrogen surface coverages were investigated by ab initio thermodynamics, based on which i) the phase diagrams of stable H coverage on each plane were obtained; ii) the morphology evolutions of the Co nanocrystals with various surface hydrogen coverages were further predicted by Wulff construction. It was revealed that addition of H2 could change the facet stability generating diverse morphological Co nucleus. The kinetic role of H2 was further explored by DFT during adatom Co surface diffusions under different H coverages (?H), which suggested that except for Co(100) at ?H of 0.56 ML, the surface H would hinder Co surface diffusions. The projected density of state (PDOS) gave a deeper insight that the electronic structures of the Co adatom could be alerted by addition of the surface H, which thereby affected its surface diffusion ability. X. Wang, N. Liu, Q. Zhang, X. Liang, and B. Chen were financially supported by the National Natural Science Foundation of China (No. 21476012, 21571012, and 91534201). D. Mei was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. Computing time was granted by the facility user proposal grand challenge of computational catalysis of the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL). EMSL is a national scientific user facility located at Pacific Northwest National Laboratory (PNNL) and sponsored by DOE’s Office of Biological and Environmental Research.

Nanocrystal Co, Surface free energy, Surface diffu↗
DOE OSTI2016

Theoretical kinetics of O + C 2 H 4

The reaction of atomic oxygen with ethylene is a fundamental oxidation step in combustion and is prototypical of reactions in which oxygen adds to double bonds. For 3 O+C 2 H 4 and for this class of reactions generally, decomposition of the initial adduct via spin-allowed reaction channels on the triplet surface competes with intersystem crossing (ISC) and a set of spin-forbidden reaction channels on the ground-state singlet surface. The two surfaces share some bimolecular products but feature different intermediates, pathways, and transition states. In addition, the overall product branching is therefore a sensitive function of the ISC rate. The 3 O+C 2 H 4 reaction has been extensively studied, but previous experimental work has not provided detailed branching information at elevated temperatures, while previous theoretical studies have employed empirical treatments of ISC. Here we predict the kinetics of 3 O+C 2 H 4 using an ab initio transition state theory based master equation (AITSTME) approach that includes an a priori description of ISC. Specifically, the ISC rate is calculated using Landau–Zener statistical theory, consideration of the four lowest-energy electronic states, and a direct classical trajectory study of the product branching immediately after ISC. The present theoretical results are largely in good agreement with existing low-temperature experimental kinetics and molecular beam studies. Good agreement is also found with past theoretical work, with the notable exception of the predicted product branching at elevated temperatures. Above ~1000 K, we predict CH 2 CHO+H and CH 2 +CH 2 O as the major products, which differs from the room temperature preference for CH 3 +HCO (which is assumed to remain at higher temperatures in some models) and from the prediction of a previous detailed master equation study.

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

Ethanol pyrolysis kinetics using H 2 O time history measurements behind reflected shock waves

The thermal decomposition of ethanol has been studied under pyrolytic conditions behind reflected shock waves in the 1250 to 1677 K temperature range, at an average pressure of 1.31 atm for a mixture highly diluted in Ar. A laser absorption technique was utilized to measure H 2 O time-histories, and the detailed kinetics mechanism (AramcoMech2.0) was selected among various models from the literature based on its a priori agreement with the experimental data. Sensitivity and rate-of-production analyses were performed and showed that the C 2 H 5 OH→C 2 H 4 + H 2 O (R1) decomposition pathway is almost the sole reaction contributing to H 2 O formation at the early times under the present conditions, allowing an a priori direct measurement of its rate coefficient k 1 . The rate coefficient was determined to be defined as the Arrhenius equation k 1 (s -1 ) = 3.37 × 10 11 exp (–27174 K/T), which is in very good agreement with Kiecherer et al. (2015), where k 1 was also directly determined under similar conditions. Secondary chemical reactions taking place in the thermal decomposition have very low influence in the H 2 O formation during the time-frame selected, leading to an uncertainty for k 1 of approximately 20%. The full H 2 O time histories are useful for validating the full ethanol kinetics mechanism for future validation.

42 ENGINEERING↗
DOE OSTI2023

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

Here, the growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O 2 ] n , and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Kinetics of low temperature plasma assisted NH 3 /H 2 oxidation in a nanosecond-pulsed discharge

Ammonia (NH 3 ) has been widely recognized as one of the carbon-neutral fuels. However, ammonia combustion suffers low reactivity and high N 2 O/NO x emissions. Here, to overcome these issues, this work reports plasma assisted NH 3 /H 2 oxidation and unveils the kinetics of fuel oxidation and N 2 O/NO x formation by combining time-resolved laser diagnostics with plasma modeling. Firstly, we found that the NH 3 consumption is promoted with a H 2 blending ratio of 0.3, due to enhancements of H and OH formation by plasma assisted H 2 dissociation. Secondly, at a high reduced electric field, when the H 2 blending ratio increases, the NH 3 oxidation is promoted due to both the HO 2 formation and strong NO kinetic enhancement via NO-HO 2 and NO 2 -H pathways. In the meantime, it is shown that the NO mole fraction also increases with H 2 blending ratio, because the NO formation is enhanced via N( 2 D)-O 2 pathways, and the DeNO x chemistry is weakened with less NH 2 production. By contrast, at a lower reduced electric field, when the H 2 blending ratio increases, the decreased N( 2 D) formation does not produce enough NO to replenish the NO formation drop caused by lower NH 3 concentration. Thirdly, the reduced electric field non-monotonically affects fuel consumption and N 2 O/NO x formation by manipulating electron energy deposition pathways. The NH 3 consumption is maximized with an optimal reduced electric field where N 2 * excitation and O 2 dissociation are most efficient. When the reduced electric field deviates from its optimum, the NH 3 consumption decreases due to the discharge energy deposition to either vibrational excitation or dissociation of N 2 . The N 2 O/NO x emissions governed by the NH 3 oxidation follow the above NH 3 consumption trend.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Reactions of DGEBA epoxy cured with diethanolamine: Isoconversional kinetics and implications to network structure

The curing of diglycidyl ether of bisphenol A (DGEBA) epoxy with diethanolamine (DEA) is studied. DEA has three reactive groups, a secondary amine hydrogen and two hydroxyls. The secondary amine reacts rapidly, forming an adduct containing tertiary amines, epoxides and hydroxyls. The epoxides and hydroxyls then react in the presence of the amines to crosslink and vitrify the epoxy in the “gelation” reaction. The gelation reaction, the subject of this study, is not simple. The reaction exhibits unusual dependencies on both temperature and degree of cure. Previously, the general mechanisms of this curing process were explored by a number of us. In the present paper, both differential scanning calorimetry (DSC) and isothermal microcalorimetry (IMC) are used to determine a number of characteristic times associated with the reaction. The characteristic times show that the reaction rate has different functional forms at different temperatures and extents of reaction. This results from the reaction rate not depending solely upon the temperature and over-all extent-of-reaction. The concentration of a number of auxiliary reactive species that are generated in the course of the reaction (as well as their mobility and steric hindrance) appear to be key factors in defining the reaction kinetics. The dependence of the final network structure on cure schedule for this type of tertiary amine activated reaction is then discussed in the context of the literature. Lastly, in the Supplementary Material, Kamal-like functions are fit to the isothermal reaction kinetics, with the reader cautioned in applying the functions to non-isothermal cures.

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

The regrowth kinetic of the surviving population is independent of acute and chronic responses to temozolomide in glioblastoma cell lines

Chemotherapy acts on cancer cells by producing multiple effects on a cell population including cell cycle arrest, necrosis, apoptosis and senescence. However, often a subpopulation of cells survives and the behavior of this subpopulation, which is responsible for cancer recurrence, remains obscure. Here we investigated the in vitro short- and long-term responses of six glioblastoma cell lines to clinically relevant doses of temozolomide for 5 days followed by 23 days of recovery, mimicking the standard schedule used in glioblastoma patient for this drug. These cells presented different profiles of sensitivity to temozolomide with varying levels of cell cycle arrest, autophagy and senescence, followed by a regrowth of the surviving cells. The initial reduction in cell number and the subsequent regrowth was analyzed with four new parameters applied to Cumulative Population Doubling (CPD) curves that describe the overall sensitivity of the population and the characteristic of the regrowth: the relative end point CPD (RendCPD); the relative Area Under Curve (rAUC); the Relative Time to Cross a Threshold (RTCT); and the Relative Proliferation Rate (RPR). Surprisingly, the kinetics of regrowth were not predicted by the mechanisms activated after treatment nor by the acute or overall sensitivity. With this study we added new parameters that describe key responses of glioblastoma cell populations to temozolomide treatment. These parameters can also be applied to other cell types and treatments and will help to understand the behavior of the surviving cancer cells after treatment and shed light on studies of cancer resistance and recurrence. - Highlights: • Little is known about the behavior of the glioma cells surviving to TMZ. • The short- and long-term response of six glioma cells lines to TMZ varies considerably. • These glioma cells lines recovered proliferation after therapeutic levels of TMZ. • The growth velocity of the surviving cells was different from the untreated cells. • The kinetic of regrowth was not predicted by any TMZ-triggered mechanism.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

On stability criteria for kinetic magnetohydrodynamics

The existence of a potential energy functional in the zero-Larmor-radius collisionless plasma theory of Kruskal & Oberman (Phys. Fluids, vol. 1, 1958 p. 275), Rosenbluth & Rostoker (Phys. Fluids, vol. 2, 1959, p. 23) allows us to derive easily sufficient conditions for linear stability. However, this kinetic magnetohydrodynamics (KMHD) theory does not have a self-adjointness property, making it difficult to derive necessary conditions. In particular, the standard methods to prove that an instability follows if some trial perturbation makes the incremental potential energy negative, which rely on the self-adjointness of the force operator or on the existence of a complete basis of normal modes, are not applicable to KMHD. This paper investigates KMHD linear stability criteria based on the time evolution of initial-value solutions, without recourse to the classic bounds or comparison theorems of Kruskal–Oberman and Rosenbluth–Rostoker for the KMHD potential energy. The adopted approach does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Most importantly, the investigation of a necessary condition for stability does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes. Thus, it is thereby shown that stability in isothermal ideal-MHD is a sufficient condition for stability in KMHD and that, with a proviso on the long-time behaviour of oscillations about stable equilibria, stability in the double-adiabatic fluid theory, including the variation of the parallel fluid displacement, would be a necessary condition for stability in KMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A drift-kinetic analytical model for scrape-off layer plasma dynamics at arbitrary collisionality

We report a drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite–Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously derived fluid models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Drift kinetic theory of alpha transport by tokamak perturbations

We report that upcoming tokamak experiments fuelled with deuterium and tritium are expected to have large alpha particle populations. Such experiments motivate new attention to the theory of alpha particle confinement and transport. A key topic is the interaction of alpha particles with perturbations to the tokamak fields, including those from ripple and magnetohydrodynamic modes like Alfvén eigenmodes. These perturbations can transport alphas, leading to changed localization of alpha heating, loss of alpha power and damage to device walls. Alpha interaction with these perturbations is often studied with single-particle theory. In contrast, we derive a drift kinetic theory to calculate the alpha heat flux resulting from arbitrary perturbation frequency and periodicity (provided these can be studied drift kinetically). Novel features of the theory include the retention of a large effective collision frequency resulting from the resonant alpha collisional boundary layer, correlated interactions over many poloidal transits and finite orbit effects. Heat fluxes are considered for the example cases of ripple and the toroidal Alfvén eigenmode (TAE). The ripple heat flux is small. The TAE heat flux is significant and scales with the square of the perturbation amplitude, allowing the derivation of constraints on mode amplitude for avoidance of significant alpha depletion. A simple saturation condition suggests that TAEs in one upcoming experiment will not cause significant alpha transport via the mechanisms in this theory. However, saturation above the level suggested by the simple condition, but within numerical and experimental experience, which could be accompanied by the onset of stochasticity, could cause significant transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Unsupervised discovery of nonlinear plasma physics using differentiable kinetic simulations

Plasma supports collective modes and particle–wave interactions that lead to complex behaviour in, for example, inertial fusion energy applications. While plasma can sometimes be modelled as a charged fluid, a kinetic description is often crucial for studying nonlinear effects in the higher-dimensional momentum–position phase space that describes the full complexity of the plasma dynamics. We create a differentiable solver for the three-dimensional partial-differential equation describing the plasma kinetics and introduce a domain-specific objective function. Using this framework, we perform gradient-based optimization of neural networks that provide forcing function parameters to the differentiable solver given a set of initial conditions. We apply this to an inertial-fusion-relevant configuration and find that the optimization process exploits a novel physical effect.

plasma nonlinear phenomena↗
DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗