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

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fiber-Coupled Multipass NIR Sensor for In Situ, Real-Time Water Vapor Outgassing Monitoring

This work presents the recent development of a fiber-coupled multipass near-infrared (NIR) gas sensor used to monitor water vapor desorption of small material coupons. The gas sensor design employs a White cell topology to maximize the optical path length over a compact, hand-size footprint. Water vapor concentrations are quantified over a large dynamic range by simultaneously applying wavelength modulation and tunable diode laser absorption spectroscopy techniques. A custom headspace optimized for material desorption experiments is assembled using commercially available vacuum chamber components. We provide in situ measurements of water vapor desorption from two geometries of the industrially important silicone elastomer Sylgard-184 as a case study for sensor viability. To corroborate the results, the gas sensor data are compared to numerical simulations based on a triple-mode diffusion–sorption model, consisting of Henry, Langmuir, and Pooling modes.

gas sensor↗
DOE OSTI2023

Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics

This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.

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

Pervasive symmetry-lowering nanoscale structural fluctuations in the cuprate La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4

The cuprate superconductors are among the most widely studied quantum materials, yet there remain fundamental open questions regarding their electronic properties and the role of structural degrees of freedom. Recent neutron and x-ray scattering measurements uncovered exponential scaling with temperature of the strength of orthorhombic fluctuations in the tetragonal phase of La 2−𝑥⁢ Sr 𝑥 ⁢Cu⁢O 4 and Tl 2 ⁢Ba 2 ⁢Cu⁢O 6+𝛿 , unusual behavior that closely resembles prior results for the emergence of superconducting fluctuations, and that points to a common origin rooted in inherent correlated structural inhomogeneity. Here, in this study, we use neutron and x-ray diffuse scattering to provide further insight into structural fluctuations in La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4 . We perform measurements up to temperatures that approach the melting point of the undoped parent compound, and we investigate the effects of in situ in-plane uniaxial elastic stress and plastic deformation. Our neutron scattering results for the parent compound La 2 ⁢Cu⁢O 4 reveal that short-range orthorhombic fluctuations persist to the maximum experimental temperature of nearly 1000 K, i.e., to a significant fraction of the crystallization temperature. At this temperature, the spatial characteristic length extracted from the momentum-space data is still about three lattice constants. The neutron scattering experiment enables quasistatic discrimination and reveals that the response is increasingly dynamic at higher temperatures. We also find that neither compressive elastic stress along the tetragonal [110] direction nor irreversible plastic deformation activating the ⟨100⟩⁢[001] slip system significantly alter the robust orthorhombic fluctuations. Overall, these results support the notion that these fluctuations couple to subtle, underlying inhomogeneity that underpins the cuprate phase diagram. Finally, in La 1.8⁢ Sr 0.2⁢ Cu⁢O 4 , we uncover low-energy structural fluctuations at a nominally forbidden reflection that are distinct from the orthorhombic distortions. While the origin of these fluctuations is not clear, they might be related to the presence of extended defects such as dislocations or stacking faults.

Spieker, R. J. [Univ. of Minnesota, Minneapolis, M↗
DOE OSTI2021

Advanced battery modeling for interfacial phenomena and optimal charging

Lithium ion batteries are one of the most promising energy storage systems for portable devices, transportation, and renewable grids. To meet the increasing requirements of these applications, higher energy density and areal capacity, long cycle life, fast charging rate and enhanced safety for lithium ion battery (LIBs) are urgently needed. To solve these challenges, the relevant physics at different length scale need to be understood. However, experimental study is time consuming and limited in small scale’s study. Modeling techniques provide us powerful tools to get a deep understanding of the relevant physics and find optimal solutions. This work focuses on studying the mechanism in advanced battery engineering techniques and developing a new charging algorithm by model-based optimization. The research topics are divided into six topics and each topic is reported as a form of journal publication. Paper Ⅰ provides a new aspect of how ALD coating improves the lithium ion diffusion at electrode particles. Paper Ⅱ explains the mechanisms by which 3D electrodes enhance battery performance and reveals guidelines for optimized 3D electrode designs by a 3D electrochemical-mechanical battery model. Paper Ⅲ investigates the electrolyte concentration impact on SEI layer growth and Li plating, especially under high charge rates. Paper Ⅳ proposes an optimized charging protocol for fast charging for reducing the charging time with minimal degradation. Paper Ⅴ reports a comprehensive degradation model for degradation estimation and life predication of energy storage system (ESS). Paper Ⅵ is a study of temperature-dependent state of charge (SOC) estimation for battery pack.

25 ENERGY STORAGE↗
DOE OSTI2024

Multiscale approaches for optimizing the impact of strain on Na-ion battery cycle life

Abstract The high costs and geopolitical challenges inherent to the lithium-ion (Li-ion) battery supply chain have driven a rising interest in the development of sodium-ion (Na-ion) batteries as a potential alternative. Unfortunately, the larger ionic radius of Na limits the reversibility of cycling because of the extensive atomic rearrangements that accompany Na-ion insertion, which in turn limit diffusion and charging speed, and lead to rapid degradation of the electrodes. The Center for Strain Optimization for Renewable Energy (STORE) was established to address these challenges and develop new electrode materials for Na-ion cells. This article discusses the current state-of-the-art materials used in Na-ion cells and several directions that STORE believes are critical to understand and control the structural and volumetric changes during the reversible (de)insertion of large cations. Graphical abstract Highlights Understanding the fundamental way materials respond to localized strains at the atomic length-scale is a critical first step in the development of highly reversible, long cycle life, Na-ion insertion hosts. This perspective explores a variety of methods that can be employed to mitigate the detrimental effects of large strain. The insights gained from these investigations should help lay the foundation for the creation of more economical and sustainable batteries that could have immediate impact on global energy infrastructure. Discussion Although there is near universal agreement that electrochemical energy storage must be an integral part of a green-energy future, there is less agreement about how to reduce the cost of energy storage. Replacing high-cost lithium-ion cells with lower-cost sodium-ion batteries is one option frequently considered in future energy models, but the details of what can be achieve with optimized sodium cell performance remains unclear. Here we posit that developing methods to mitigating strain on the electrode particle length scale is a key factor for achieving long-cycle-life sodium-ion batteries. Mitigating strain on the atomic scale suppress electrode-level volume change. Allowing for fast cycling in materials without the problems of electrode cracking or delamination. We further posit that understanding volume change in sodium-ion electrodes at a fundamental level will lead to the designing new sodium-ion electrode materials that will allow for efficient, stable, lower-cost energy storage.

Brady, Michael J.↗
DOE OSTI2024

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

First-principles comparative study of Cr migration in O3 and O/P hybrid-phased NaCrO 2

In layered Na transition-metal (TM) oxides, TM migration usually occurs at highly charged states and severely deteriorates the capacity and reversibility. Meanwhile, the formation of hybrid phases with the intergrowth of octahedral (O-type) and prismatic (P-type) Na layers also takes place at highly charged states. These hybrid phases are often more stable than simple O3 or P3 stackings. However, little is known about the mechanism and impact of TM migration in these hybrid phases. In this work, a comparative first-principles study is performed to understand the connections between structural changes and Cr migration in layered O3 and hybrid-phased NaCrO 2 . After Cr migration, the hybrid-phased NaCrO 2 suffers from greater layer shrinkage than the O3 phase. Three factors affect the Cr migration energy E mig : the Na concentration, local 3D configurations, and 2D in-plane geometries. Low Na concentrations and certain local 3D configurations facilitate the Cr migration. The Cr migration barriers in both O3 and hybrid-phased NaCrO 2 are positively correlated with the Cr E mig . The Cr migration in 17 doped O3 and hybrid-phased NaCrO 2 is surveyed. In these doped NaCrO 2 , a more uniform distribution of the Cr–O bond lengths usually suggests suppressed Cr migration. Optimal dopants for suppressing Cr migration are identified by considering E mig for both Cr and the dopant. Finally, our comparative study on Cr migration in O3 and hybrid-phased NaCrO 2 reveals the significant role of hybrid-phased structures in the development of layered cathode materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Physics-Informed Generative Adversarial Networks for Stochastic Differential Equations

We developed a new class of physics-informed generative adversarial networks (PI-GANs) to solve forward, inverse, and mixed stochastic problems in a unified manner based on a limited number of scattered measurements. Unlike standard GANs relying solely on data for training, here we encoded into the architecture of GANs the governing physical laws in the form of stochastic differential equations (SDEs) using automatic differentiation. In particular, we applied Wasserstein GANs with gradient penalty (WGAN-GP) for its enhanced stability compared to vanilla GANs. We first tested WGAN-GP in approximating Gaussian processes of different correlation lengths based on data realizations collected from simultaneous reads at sparsely placed sensors. We obtained good approximation of the generated stochastic processes to the target ones even if there is a mismatch between the input noise dimensionality and the effective dimensionality of the target stochastic processes. We also studied the overfitting issue for both the discriminator and the generator, and we found that overfitting occurs also in the generator in addition to the discriminator as previously reported. Subsequently, we considered the solution of elliptic SDEs requiring approximations of three stochastic processes, namely the solution, the forcing, and the diffusion coefficient. Here again, we assumed data collected from simultaneous reads at a limited number of sensors for the multiple stochastic processes. Furthermore, three generators were used for the PI-GANs: two of them were feed forward deep neural networks (DNNs), while the other one was the neural network induced by the SDE. For the case where we have one group of data, we employed one feed forward DNN as the discriminator, while for the case of multiple groups of data we employed multiple discriminators in PI-GANs. We solved forward, inverse, and mixed problems without changing the framework of PI-GANs, obtaining both the means and the standard deviations of the stochastic solution and the diffusion coefficient in good agreement with benchmarks. In this work, we have demonstrated the effectiveness of PI-GANs in solving SDEs for about 120 dimensions. In principle, PI-GANs could tackle very high dimensional problems given more sensor data with low-polynomial growth in computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2021

Spatiotemporal Crossover between Low- and High-Temperature Dynamical Regimes in the Quantum Heisenberg Magnet

The stranglehold of low temperatures on fascinating quantum phenomena in one-dimensional quantum magnets has been challenged recently by the discovery of anomalous spin transport at high temperatures. Whereas both regimes have been investigated separately, no study has attempted to reconcile them. For instance, the paradigmatic quantum Heisenberg spin-$1/2$ chain falls at low-temperature within the Tomonaga-Luttinger liquid framework, while its high-temperature dynamics is superdiffusive and relates to the Kardar-Parisi-Zhang universality class in $1+1$ dimensions. This work aims at reconciling the two regimes. Building on large-scale matrix product state simulations, we find that they are connected by a temperature-dependent spatiotemporal crossover. In this work, as the temperature $T$ is reduced, we show that the onset of superdiffusion takes place at longer length and time scales $\propto 1/T$. This prediction has direct consequences for experiments including nuclear magnetic resonance: it is consistent with earlier measurements on the nearly ideal Heisenberg $S=1/2$ chain compound Sr$_2$CuO$_3$ yet calls for new and dedicated experiments.

1-dimensional spin chains↗
DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗
DOE OSTI2021

Brought to Light. I. Quantification of Disk Substructure in Dwarf Early-type Galaxies

Dwarf early-type galaxies (ETGs) display a rich diversity in their photometric, structural, and dynamical properties. In this work, we address their structural complexity by studying with deep imaging a sample of nine dwarf ETGs from the Virgo galaxy cluster, characterized by having faint disk features such as bars and spiral arms, that lie mostly hidden within the bright diffuse light of the galaxies. We present a new, robust method that aims to identify and extract the disk substructure embedded in these dwarf ETGs. The method consists in an iterative procedure that gradually separates a galaxy image into two components: the bright, dominant, diffuse component, and the much fainter, underlying disk component. By applying it to the dwarf ETG sample, we quantify their disk substructure and find that its relative contribution to the total galaxy light ranges between 2.2% and 6.4% within two effective radii. We test the reliability of the method, and prove that it is accurate in recovering the disk substructure we introduce in mock galaxy images, even at low disk-to-total light fractions of a few percent. As a potential application of the method, we perform a Fourier analysis on the extracted disk substructures and measure the orientation, length, and strength of the bars, as well as the pitch angle and strength of the spiral arms. We also briefly discuss a scenario based on the numerical simulations presented in our companion paper, Brought to Light II: Smith et al., in which we investigate the origins of the substructure in such dwarf systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Venturing into Unexplored Phase Space: Synthesis, Structure, and Properties of MgCo 3 B 2 Featuring a Rumpled Kagomé Network

MgCo 3 B 2 , a novel ternary boride in a previously unexplored phase space, was synthesized using the hydride route. In situ powder X-ray diffraction and DFT calculations aided in the discovery of this compound, whose structure was then determined by single-crystal X-ray diffraction. Like the closely related CeCo 3 B 2 , MgCo 3 B 2 crystallizes in centrosymmetric space group P6/mmm (a = 4.883(2) Å, c = 2.926(2) Å at 210 K, Z = 1). Unlike CeCo 3 B 2 , however, it adopts a disordered structure that features a rumpled Kagomé network of Co atoms, and Mg atoms fill the channels of a Co–B framework. Although the structural disorder leads to motifs that are similar to those observed in MgNi 3 B 2 and other related ternary borides, no evidence of an ordered superstructure was found by single-crystal X-ray diffraction or high-resolution powder X-ray diffraction. In the case of CeCo 3 B 2 , boron atoms occupy the center of regular Co 6 trigonal prisms; in MgCo 3 B 2 , boron atoms are shifted from the center of the prism to form B–B dimers with roughly the same length as those found in MgNi 3 B 2 . Magnetic susceptibility data exhibit an unusual temperature dependence that cannot be convincingly modeled by the modified Curie–Weiss equation, consistent with DFT calculations predicting a nonmagnetic ground state. Intrinsic susceptibility at 300 K is 1.42 × 10 –3 emu/mol Oe, which is comparable to that of paramagnetic YCo 3 B 2 and CeCo 3 B 2 with a similar structure and composition. Here, this study showcases the efficacy of combining several methodologies to discover new solids in unexplored phase spaces. This approach includes in situ PXRD data to monitor reactions of precursors upon heating, a diffusion-enhanced synthesis method, and DFT assessment of compound stability.

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

Crystallization Mechanisms of Poly(vinylidene Fluoride- co -chlorotrifluoroethylene): Nucleation Transitions, Growth Kinetics, and Microstructure Evolution

The crystallization kinetics of FK-800, a commercial semicrystalline copolymer of chlorotrifluoroethylene and vinylidene fluoride (poly(CTFE-co-VDF)), dictate its microstructure and functional performance. Crystallization occurs between the glass transition (T g ≈ 31 °C) and melting point (T m ≈ 110 °C), with grain boundaries playing a key role in applications such as memristors, where they regulate filament growth. This study integrates in situ atomic force microscopy (AFM), grazing-incidence wide-angle X-ray scattering (GIWAXS), and differential scanning calorimetry (DSC) to characterize nucleation, growth, and morphology across the full crystallization window. Hot-stage AFM reveals a transition from homogeneous to heterogeneous nucleation near 45 °C and a shift from reaction-limited to diffusion-limited growth above 60 °C. Kinetics are analyzed using Lauritzen–Hoffman, Turnbull–Fisher, and Avrami models to extract surface free energies and activation barriers. A multiscale modeling framework is developed in which nucleation densities and growth rates measured by AFM are used to reconstruct Avrami kinetics, which are then validated against coverage data, and extrapolated to predict bulk crystallization behavior observed by DSC. GIWAXS confirms a primarily flat-on chain orientation at all temperatures, linking orientation to domain morphology. This integrative approach quantitatively connects nanoscale crystallization dynamics with macroscopic phase evolution, establishing a quantitative framework for modeling crystallization across time, temperature, and length scales.

Crystallization↗
DOE OSTI2025

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗