Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular geometry”

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

Consequences of Convex Nanopore Chemistry on Confined Water Dynamics

A basic understanding of confined water is crucial for developing selective ion transport and water purification membranes, yet the roles of nanopore geometry and functionality on confined water dynamics remain unresolved. We discuss the synthesis of perdeuterated ionic alkylsulfonate amphiphiles and their water-induced self-assembly into lyotropic liquid crystal (LLC) mesophases with well-defined, convex, sulfonate-lined nanopores. Quasielastic neutron scattering (QENS) measurements demonstrate that the water self-diffusion coefficients within these sulfonate-lined convex nanopores depend on the hydration level and amphiphile counterion identity (H + , K + , NMe 4 + ). The consistency of the observed counterion-dependent water dynamics trends with those of carboxylate LLCs is rationalized on the basis of similarities in the counterion spatial distributions in the water-filled channels, which we deduce from electron density maps derived from small-angle X-ray scattering (SAXS) analyses. These findings indicate that water diffusion is systematically faster in sulfonate-lined nanopores as compared to carboxylate-lined pores due to weaker water interactions with the softer and more hydrophobic–SO 3 – functionalities. These molecular-level insights into the relationships between convex pore wall chemical functionalities, hydrated counterions, and confined water diffusion may inform future development of new nanoporous media.

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

Engineering advancements in microfluidic systems for enhanced mixing at low Reynolds numbers

Mixing within micro- and millichannels is a pivotal element across various applications, ranging from chemical synthesis to biomedical diagnostics and environmental monitoring. The inherent low Reynolds number flow in these channels often results in a parabolic velocity profile, leading to a broad residence time distribution. Achieving efficient mixing at such small scales presents unique challenges and opportunities. This review encompasses various techniques and strategies to evaluate and enhance mixing efficiency in these confined environments. It explores the significance of mixing in micro- and millichannels, highlighting its relevance for enhanced reaction kinetics, homogeneity in mixed fluids, and analytical accuracy. We discuss various mixing methodologies that have been employed to get a narrower residence time distribution. The role of channel geometry, flow conditions, and mixing mechanisms in influencing the mixing performance are also discussed. Various emerging technologies and advancements in microfluidic devices and tools specifically designed to enhance mixing efficiency are highlighted. We emphasize the potential applications of micro- and millichannels in fields of nanoparticle synthesis, which can be utilized for biological applications. Additionally, the prospects of machine learning and artificial intelligence are offered toward incorporating better mixing to achieve precise control over nanoparticle synthesis, ultimately enhancing the potential for applications in these miniature fluidic systems.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

pH-Driven Restructuring of Hydration Layers and Cation Adsorption at the Alumina–Water Interface

Oxide−water interfaces underpin ion separations, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. However, the molecular-scale links between pHdependent surface protonation, hydration-layer structure, and counterion adsorption remain poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)−water interface in 10 mM Rb + over pH 3−12. CTR measurements reveal two distinct adsorbed water layers at ∼2.2 and ∼3.5 Å above the surface. Each water layer shifts toward the substrate at transition pHs near 6.5 and 10.6, respectively, reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb + coverage and a decrease in mean adsorption height from ∼3.5 to ∼2.7 Å with increasing pH, indicating enhanced counterion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., the potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥ 3.5. The negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb + coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid−base chemistry to ion binding distances at an oxide−water interface.

Adsorption↗
DOE OSTI2025

Enabling Industrial Re-Use of Large-Format Additive Manufacturing Molding and Tooling

Large-format additive manufacturing (LFAM) is an enabling manufacturing technology capable of producing large parts with highly complex geometries for a wide variety of applications, including automotive, infrastructure/construction, and aerospace mold and tooling. In the past decade, the LFAM industry has seen widespread use of bio-based, glass, and/or carbon fiber reinforced thermoplastic composites which, when printed, serve as a lower-cost alternative to metallic parts. One of the highest-volume materials utilized by the industry is carbon fiber (CF)-filled polycarbonate (PC), which in out-of-autoclave applications can achieve comparable mechanical performance to metal at a significantly lower cost. Previous work has shown that if this material is recovered at various points throughout the manufacturing process for both the lab and pilot scale, it can be mechanically recycled with minimal impacts on the functional performance and printability of the material while significantly reducing the feedstock costs. End-of-life (EOL) CF-PC components were processed through industrial shredding, melt compounding, and LFAM equipment, followed by evaluation of the second-life material properties. Experimental assessments included quantitative analysis of fiber length attrition, polymer molecular weight degradation using gel permeation chromatography (GPC), density changes via pycnometry, thermal performance using dynamic mechanical analysis (DMA), and mechanical performance (tensile properties) in both the X- and Z-directions. Results demonstrated a 24.6% reduction in average fiber length compared to virgin prints, accompanied by a 21% decrease in X-direction tensile strength and a 39% reduction in tensile modulus. Despite these reductions, Z-direction tensile modulus improved by 4%, density increased by 6.8%, and heat deflection temperature (HDT) under high stress retained over 97% of its original value. These findings underscore the potential for integrating mechanically recycled CF-PC into industrial LFAM applications while highlighting the need for technological innovations to mitigate fiber degradation and enhance material performance for broader adoption. This critical step toward circular material practices in LFAM offers a pathway to reducing feedstock costs and environmental impact while maintaining functional performance in industrial applications.

additive manufacturing↗
DOE OSTI2020

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

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

Sulfur Lone Pairs Control Topology in Heterotrimetallic Complexes: An Experimental and Theoretical Study

Heterotrimetallic complexes with (N 2 S 2 )M metallodithiolates, M = Ni 2+ , [Fe(NO)] 2+ , and [Co(NO)] 2+ , as bidentate chelating ligands to a central trans-Cr(NO)(MeCN) unit were characterized as the first members of a new class, NiCrNi, FeCrFe, CoCrCo. The complexes exhibit a cisoid structural topology, ascribed to the stereoactivity of the available lone pair(s) on the sulfur donors, resulting in a dispersed, electropositive pocket from the N/N and N/S hydrocarbon linkers wherein the Cr-NO site is housed. Computational studies explored alternative isomers (transoid and inverted cisoid) that suggest a combination of electronic and steric effects govern the geometrical selectivity. Electrostatic potential maps readily display the dominant electronegative potential from the sulfurs which force the NO to the electropositive pocket. The available S lone pairs work in synergy with the π-withdrawing ability of NO to lift Cr out of the S 4 plane toward the NO and stabilize the geometry. The metallodithiolate ligands bound to Cr(NO) thus find structural consistency across the three congeners. Although the dinitrosyl [(bme-dach)Co(NO)-Mo(NO)(MeCN)-(bme-dach)Co(MeCN)][PF 6 ] 2 (CoMoCo') analogue displays chemical noninnocence and a partial Mo–Co bond toward (N 2 S 2 )Co'(NCCH 3 ) in an “asymmetric butterfly” topology, the stability of the {Cr(NO)} 5 unit prohibits such bond rearrangement. Magnetism and EPR studies illustrate spin coupling across the sulfur thiolate sulfur bridges.

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

Molecular Layer Etching of Metalcone Films Using Lithium Organic Salts and Trimethylaluminum

Advances in semiconductor device manufacturing are limited by our ability to precisely add and remove thin layers of material in multistep fabrication processes. Recent reports on atomic layer etching (ALE) have provided the means for the precise removal of inorganic thin films deposited by atomic layer deposition (ALD), opening new avenues for nanoscale device design. Here, we report on a new technique for the precise removal of metal organic thin films deposited by molecular layer deposition (MLD), which we term molecular layer etching. This etching process employs sequential exposures of lithium organic salt (LOS) and trimethylaluminum (TMA) precursors to produce self-limiting etching behavior. We employ quartz crystal microbalance experiments to demonstrate (i) etching of alucone films preloaded with LOS upon TMA exposures and (ii) layer-by-layer etching of alucone films using alternating exposures of LOS and TMA. We also identify the selectivity of these etching mechanisms. We probe the mechanism for the layer-by-layer etching of alucone using a quartz crystal microbalance and Fourier transform infrared spectroscopy and identify that the etching proceeds via heterolytic cleaving of Al-O bonds in alucone upon LOS exposure followed by methylation to produce volatile species upon TMA exposure. The etching process results in the removal of 0.4 nm/cycle of alucone at 160 °C and up to 3.6 nm/cycle of alucone at 266 °C in ex situ etching experiments on silicon wafers. This halogen-free etching process enables etching of MLD films and provides new fabrication pathways for the control of material geometries at the nanoscale.

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

Colloquium : Sliding and pinning in structurally lubric 2D material interfaces

A plethora of two-dimensional (2D) materials have been introduced in physics and engineering in the past two decades. Their robust, membranelike sheets permit (mostly require) deposition, giving rise to solid-solid dry interfaces whose mobility, pinning, and general tribological properties under shear stress are currently being understood and controlled, both experimentally and theoretically. Here, in this Colloquium simulated case studies of twisted graphene systems are used as a prototype workhorse tool to demonstrate and discuss the general picture of 2D material interface sliding. First highlighted is the crucial mechanical difference, often overlooked, between small and large incommensurabilities, which corresponds to, for example, small and large twist angles in graphene interfaces. In both cases, focusing on flat, structurally lubric or “superlubric” geometries, the generally separate scalings with the area of static friction in pinned states and of kinetic friction during sliding are elucidated and reviewed, tangled as they are with the effects of velocity, temperature, load, and defects. The roles of island boundaries and elasticity are also discussed, and compared when possible to results in the literature for systems other than graphene. It is proposed that the resulting picture of pinning and sliding should be applicable to interfaces in generic 2D materials that are of importance for the physics and technology of existing and future bilayer and multilayer systems.

2-dimensional systems↗
DOE OSTI2021

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

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

Customizable Aperture Geometry in Metal–Organic Frameworks for Kinetic Hydrocarbon Separation

Precise control of aperture dimensions is crucial in adsorptive separations of hydrocarbons, as it directly affects key parameters such as selectivity, capacity, diffusion, and recyclability. The development of metal–organic frameworks (MOFs) has enabled the fine-tuning of local pore environments to address important hydrocarbon separations. However, customizing the aperture geometry to tune the kinetic separation performance remains challenging. Here, we deploy a mixed-linker synthesis strategy, combining long and short linkers on fcu net Zr-MOFs with equilateral triangular apertures to construct isoreticular multivariate MOFs, NU-415 and NU-416, with tailored isosceles triangular apertures suitable for the separation of hexane isomers. Sorption, liquid batch separation, and X-ray diffraction measurements demonstrate significantly improved selectivity, capacity, stability, and recyclability of NU-415 and NU-416 compared with Zr-muconate and MOF-801. Notably, both NU-415 and NU-416 achieve uptake capacities of 2.2 mmol g –1 in 1 min with an n-hexane to 2,2-dimethylbutane selectivity over 200 in an equimolar ternary mixture at ambient conditions, comparable to leading reported materials. Mechanistic studies confirm that separation performance is predominantly governed by significant kinetic differences rather than by thermodynamics. Furthermore, the successful customization of aperture geometry not only enables superior linear to monobranched hexane selectivity in NU-415 but also demonstrates the mixed-linker synthesis strategy as a promising solution for precise and predictable pore architecture control in MOFs.

Adsorption↗
DOE OSTI2021

Strong-field ionization of water: Nuclear dynamics revealed by varying the pulse duration

Polyatomic molecules in strong laser fields can undergo substantial nuclear motion within tens of femtoseconds. Ion imaging methods based on dissociation or Coulomb explosion therefore have difficulty faithfully recording the geometry dependence of the field ionization that initiates the dissociation process. Here we compare the strong-field double ionization and subsequent dissociation of water (both H 2 O and D 2 O) in 10-fs and 40-fs 800-nm laser pulses. We find that 10-fs pulses turn off before substantial internuclear motion occurs, whereas rapid internuclear motion can take place during the double ionization process for 40-fs pulses. The short-pulse measurements are consistent with a simple tunnel ionization picture whose predictions help interpret the motion observed in the long-pulse measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

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

Enhanced intersystem crossing of boron dipyrromethene by TEMPO radical

Radical enhanced intersystem crossing (EISC) of organic chromophores is an important approach to generate a long-lived triplet state for various electronic and optoelectronic applications. However, structural factors and design rules to promote EISC are not entirely clear. In this work, we report a series of boron dipyrromethene (BODIPY) derivatives covalently linked with a 2,2,6,6-tetramethyl-1-piperidinyloxyl (TEMPO) radical with varying distances and topologies. We show that the incorporation of the TEMPO radical to BODIPY results in strong fluorescence quenching by up to 85% as a result of EISC and enhanced internal conversion. In BDP-2AR [2-(4-methyleneamino-TEMPO) BODIPY], a dyad with the shortest BODIPY–TEMPO through-bond distance, we observe the fastest EISC rate (τ isc = 1.4 ns) and the longest triplet excited state lifetime (τ T = 32 μs) compared to other distance and geometry variations. Contrary to previous reports and a general presumption, the BODIPY–TEMPO through-bond distance in this system does not play a significant role on the triplet formation rate and yield. Finally, density functional theory suggests a folding of the TEMPO radical to form a sandwich-like structure with a BODIPY ring that leads to a decrease in the through-space distance, providing a new and an interesting insight for the radical enhanced intersystem.

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

A Novel Laser 3D Printing Method for the Advanced Manufacturing of Protonic Ceramics

Protonic ceramics (PCs) with high proton conductivity at intermediate temperatures (300–600 °C) have attracted many applications in energy conversion and storage devices such as PC fuel/electrolysis cells, PC membrane reactors, hydrogen pump, hydrogen or water-permeable membranes, and gas sensors. One of the essential steps for fulfilling the practical utilization of these intermediate-temperature PC energy devices is the successful development of advanced manufacturing methods for cost-effectively and rapidly fabricating them with high energy density and efficiency in a customized demand. In this work, we developed a new laser 3D printing (L3DP) technique by integrating digital microextrusion-based 3D printing and precise and rapid laser processing (sintering, drying, cutting, and polishing), which showed the capability of manufacturing PCs with desired complex geometries, crystal structures, and microstructures. The L3DP method allowed the fabrication of PC parts such as pellets, cylinders, cones, films, straight/lobed tubes with sealed endings, microchannel membranes, and half cells for assembling PC energy devices. The preliminary measurement of the L3DP electrolyte film showed a high proton conductivity of ≈7 × 10−3 S/cm. This L3DP technique not only demonstrated the potential to bring the PCs into practical use but also made it possible for the rapid direct digital manufacturing of ceramic-based devices.

08 HYDROGEN↗
DOE OSTI2025

Restructuring of Benzimidazole-Based Copper Complexes during Electrochemical CO 2 Reduction

Copper-based molecular catalysts have been widely used for the electrochemical reduction of carbon dioxide (ERCO 2 ). However, establishing a correlation between structure and catalytic performance has been challenging due to the dynamic nature of these materials under ERCO 2 working conditions. Here, we report on the restructuring of the copper complexes [Cu(bzimpy)Cl 2 ] and [Cu(pyrben) 2 (NO 3 )]NO 3 during ERCO 2 in an acidic medium, evidencing its effects on the observed activity, selectivity, and stability. In situ X-ray absorption spectroscopy (XAS) suggested that, at sufficiently negative potentials, the Cu 2+ centers of both complexes are reduced to undercoordinated Cu 0 species. Pb underpotential deposition of post-ERCO 2 samples indicated that such Cu 0 species also present a dominant contribution of the (100) domain, which is in line with the high Faradaic efficiency toward C 2 products (∼45%) shown by both complexes at −1.24 V vs RHE. By combining online electrochemical mass spectrometry, in situ XAS, electron paramagnetic resonance (EPR), and quantitative product analysis, we found that [Cu(bzimpy)Cl 2 ] exhibits a partially reversible restructuring, regenerating the main complex moiety, [Cu(bzimpy)] 2+ , and its geometry.

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

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗