Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Compact Diffusion Bonded Heat Exchanger Fatigue Life Simulations (Final Report)

Compact diffusion bonded heat exchangers are essential for high pressure heat exchange, but they are subject to thermal fatigue and ramp rate limitations. Simulation of these geometries is challenging with a large range of length and time scales from thousands of mm-sized microchannels inside a m-sized heat exchanger. Multi-physics simulations including thermal, fluid, and solid mechanics components are being used to predict stress within the heat exchangers under these conditions. These predictions can then be used to understand thermal ramp rate limitations while keeping maximum stresses low as well as fatigue life predictions from well-known empirical models.

42 ENGINEERING↗
DOE OSTI2024

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

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

Impact of Annealing on Superconducting Properties of Tantalum Thin Films

The pursuit for higher-coherence superconducting qubits has led to the adoption of new materials and innovative fabrication techniques. Among the commonly used superconducting materials, tantalum-based or tantalum-capped superconducting qubits have demonstrated the longest lifetimes, outperforming niobium devices even in the presence of their native oxides. Given the losses introduced by the surface oxide, one method to further improve performance is to reduce or eliminate the native oxide on tantalum through UHV annealing. In this work, we examine the effect of dissolved oxygen as a function of annealing temperature on the superconducting properties of tantalum thin films. Dissolved oxygen is quantified using atom probe tomography (APT) and secondary ion mass spectrometry (TOF-SIMS) and correlated with electromagnetic transport measurements performed using a Physical Property Measurement System (PPMS). This study observes that annealing decomposes the native oxide and drives oxygen diffusion into the film matrix. Increasing oxygen content decreases the superconducting transition temperature (Tc) and the residual resistivity ratio (RRR). The critical field decreases with the incorporation of oxygen, while the film’s coherence length decreases relatively by 5%. These results, as well as findings from ongoing studies linking fabrication processes to the structural and chemical properties of the superconducting film and Josephson junction, provide insight into effective strategies for improving the performance of superconducting devices.

Shinde, Maithilee [IIT, Chicago; Fermilab]↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
DOE OSTI2020

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Microscopic Dynamics of Inverse Wormlike Micelles Probed Using X-ray Photon Correlation Spectroscopy

Wormlike micelles (WLMs) are ubiquitous viscoelastic modifiers that share properties with polymer solutions. While 5 their macroscopic rheology is well-understood, their microscopic dynamics are less studied because they span a large range of time 6 and length scales. In this work, we demonstrate the use of X-ray photon correlation spectroscopy to interrogate the segmental 7 dynamics of inverse WLM solutions swollen with a rubidium chloride solution. We observe a diffusive scaling of the dynamics and 8 extract a temperature-dependent diffusion coefficient, which we associate with the thermal interactions of the slow segmental 9 dynamics near entanglement points. Further, we probe this relaxation process across the unbranched to branched topological transition and 10 find no microstructural evidence of branch formation in the slow mode. Instead, we observe that the dynamics become more 11 homogeneous and prominent as the temperature is reduced and water content increases.

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

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗
DOE OSTI2019

Heat current fluctuations and anomalous transport in low-dimensional carbon lattices

Molecular dynamics simulations and methods of importance sampling are used to study the heat transport of low-dimensional carbon lattices. For both carbon nanotubes and graphene sheets, heat transport is found to be anomalous, violating Fourier’s law of conduction with a system size dependent thermal conductivity and concomitant nonlinear temperature profiles. For carbon nanotubes, the thermal conductivity is found to increase as the square root of the length of the nanotube, while for graphene sheets the thermal conductivity is found to increase as the logarithm of the length of the sheet over the system sizes considered. Furthermore, the particular length dependence and nonlinear temperature profiles place carbon lattices into a universality class with nonlinear lattice models, and suggest that heat transport through carbon nanostructures is better described by a Levy walk rather than simple diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Calibration of the charge and energy loss per unit length of the MicroBooNE liquid argon time projection chamber using muons and protons

In this paper, we describe a method used to calibrate the position- and time-dependent response of the MicroBooNE liquid argon time projection chamber anode wires to ionization particle energy loss. The method makes use of crossing cosmic-ray muons to partially correct anode wire signals for multiple effects as a function of time and position, including cross-connected TPC wires, space charge effects, electron attachment to impurities, diffusion, and recombination. The overall energy scale is then determined using fully-contained beam-induced muons originating and stopping in the active region of the detector. Using this method, we obtain an absolute energy scale uncertainty of 2% in data. We use stopping protons to further refine the relation between the measured charge and the energy loss for highly-ionizing particles. This data-driven detector calibration improves both the measurement of total deposited energy and particle identification based on energy loss per unit length as a function of residual range. As an example, the proton selection efficiency is increased by 2% after detector calibration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

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

Adaptations for gas exchange enabled the elongation of lepidopteran proboscises

The extensive biodiversification of butterflies and moths (Lepidoptera) is partly attributed to their unique mouthparts (proboscis [Pr]) that can span in length from less than 1 mm to over 280mm in Darwin's sphinx moths. Lepi-doptera, similar to other insects, are believed to inhale and exhale respiratory gases only through valve-like spi-racles on their thorax and abdomen, making gas exchange through the narrow tracheae (Tr) challenging for the elongated Pr. How Lepidoptera overcome distance effects for gas transport to the Pr is an open question that is important to understanding how the Pr elongated over evolutionary time. Here, we show with scanning electron microscopy and X-ray imaging that distance effects on gas exchange are overcome by previously unreported micropores on the Pr surface and by superhydrophobic Tr that prevent water loss and entry. We find that the density of micropores decreases monotonically along the Pr length with the maxima proportional to the Pr length and that micropore diameters produce a Knudsen number at the boundary between the slip and transition flow regimes. By numerical estimation, we further show that the respiratory gas exchange for the Pr predominantly occurs via diffusion through the micropores. In conclusion, these adaptations are key innovations vital to Pr elongation, which likely facilitated lepidopteran biodiversification and the radiation of angiosperms by coevolutionary processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Neutron spin echo spectroscopy

Neutron spin echo (NSE) spectroscopy is a powerful technique used to probe the internal dynamics and diffusion processes within matter, allowing researchers to quantify how materials respond to changes in external conditions. To fully understand the bulk properties of a material, it is often crucial to investigate dynamic processes occurring on nanometre length scales and nanosecond timescales. NSE spectroscopy provides an exceptional capability to measure such dynamics across a broad spectrum of condensed matter systems, ranging from proteins and quantum magnets to catalysts. In this Primer, we present an overview of NSE spectroscopy, highlighting the instruments used, methods of data collection and representative applications in both soft-matter and hard-matter sciences.

Czakkel, Orsolya [Institut Max von Laue-Paul Lange↗
DOE OSTI2023

Radionuclide Waste Disposal: Development of Multi-scale Experimental and Modeling Capabilities (Final Report)

The DOE EPSCoR Implementation project “Radioactive waste management: Development of multi-scale experimental and modeling capabilities” helped to develop a team of scientists and engineers from Clemson University, South Carolina State University, and the University of South Carolina to address the disposition of nuclear wastes and study the transport of radioisotopes from a waste repository in the near and far field. The project involved 20 faculty from the three institutions as well as 13 postdoctoral fellows, 32 graduate students, and 32 undergraduate students and was active from 2014-2019. Additionally, we forged new collaborations with nine researchers from DOE laboratories SRNL, LLNL, and ANL as well as the University of Manchester and the China Academy of Engineering Physics. The overarching goal of the project was to understand the conditions under which important classes of co-reactants, ranging from counter ions in crystal lattices to dissolved oxygen in pores, control the chemistry and transport characteristics of radionuclides in engineered waste forms and natural soils. Our approach was to characterize the time and length scales over which non-equilibrium states are maintained by rate-limiting, or rate-enhancing, reactions between radionuclides and co-reactants due to interactions between physical mass-transfer processes (i.e., advection, diffusion) and (biogeo) chemical reactions. We have focused our project on three specific classes of reactions relevant to radionuclide transport at DOE legacy sites: ion exchange/substitution, ligand complexation, and redox-mediated reactions. Understanding radionuclide migration requires detailed knowledge of how changes to a system – whether engineered or natural – drive the behavior of co-reactants, which in turn provide the geochemical context controlling radionuclide transport. Student engagement and training were a primary focus of the project in order to create a pipeline of researchers who could work in the area of nuclear waste disposition to support the state and the nation. Over the duration of the project we worked with 32 undergraduate, graduated 18 M.S. students and 14 Ph.D. students, and advised 13 postdoctoral fellows. The Ph.D students and postdocs have primarily taken positions at DOE laboratories, academia, and industry. Through our collaborative team, numerous follow on projects have been started with over $5M in sponsored research. Additionally, thus far the team has published 53 peer reviewed papers and given over 75 technical presentations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Guided anisotropic oxygen transport in vacancy ordered oxides

Abstract Anisotropic and efficient transport of ions under external stimuli governs the operation and failure mechanisms of energy-conversion systems and microelectronics devices. However, fundamental understanding of ion hopping processes is impeded by the lack of atomically precise materials and probes that allow for the monitoring and control at the appropriate time- and length- scales. In this work, using in-situ transmission electron microscopy, we directly show that oxygen ion migration in vacancy ordered, semiconducting SrFeO 2.5 epitaxial thin films can be guided to proceed through two distinctly different diffusion pathways, each resulting in different polymorphs of SrFeO 2.75 with different ground electronic properties before reaching a fully oxidized, metallic SrFeO 3 phase. The diffusion steps and reaction intermediates are revealed by means of ab-initio calculations. The principles of controlling oxygen diffusion pathways and reaction intermediates demonstrated here may advance the rational design of structurally ordered oxides for tailored applications and provide insights for developing devices with multiple states of regulation.

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

Temporal and spatial separations between spin glass and short-range order

Broken-symmetry-induced order parameters account for many phenomena in physics. For spin glasses, this framework dictates the theoretical construction, whereas experiments have only established dynamical behaviors but not the physical entity. Experimental techniques have limitations when the spin glass is probed as an isolated state. Here, we create an evolution from a long-range order using well-controlled non-magnetic substitution on a sublattice. Neutron magnetic diffuse scattering of pico-second timescale reveals that the dynamics of short- and long-range order formation are not affected by disorder, but their spatial ranges are. Across all specimens, the inflection point of spin correlation length’s temperature dependence fully matches with the peak in heat capacity, while spin glass can freeze at millisecond timescale either above or below this characteristic temperature of spin-order formation. Our results identify the spin glass as single, uncorrelated spins at domain walls between spin clusters, and an uncorrelated coexistence of the two.

Dronova, Margarita G [Okinawa Institute of Science↗
DOE OSTI2022

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗