Results for “confined diffusion” — Page 10

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS1991

The evolution of the diffuse cosmic ultraviolet background constrained by the Hubble Space Telescope observations of 3C 273

Results are presented of recent HST UV spectroscopy of 3C 273, which revealed more low-redshift Lyman-alpha absorption lines (IGM clouds) than expected from the extrapolation from high-redshift (not less than 1.6) observations. It is shown on the basis of the standard pressure confined cloud model of the Lyman-alpha forest that this result indicates a sharp drop in the diffuse cosmic UV background from 2 to 0 redshift. It is predicted that the H I optical depth will drop slowly or perhaps even increase with decreasing redshift at less than 2 redshift. The implied constraints on the density and pressure of the diffuse IGM at 0 redshift are also derived. The inferred evolution of the diffuse UV flux bears a striking resemblance to the most recent direct determinations of the volume emissivity of the quasar population.

Ikeuchi, Satoru↗

NASA NTRS2001

Quasi-Decadal Oscillations Generated by the QBO

Quasi-decadal oscillations (QDO) have been observed in the stratosphere and have been linked to the equatorial Quasi-Biennial Oscillation (QBO) and to the 11-year solar activity cycle. With the use of a 2D version of our Numerical Spectral Model (NSM) that incorporates Hines' Doppler Spread Parameterization (DSP) for gravity waves (GW), we demonstrate that beat periods between 9 and 11 years can be generated by the QBO as it interacts through GW filtering with the Annual Oscillation (AO) and Semi-annual Oscillation (SAO). Results are discussed from computations covering up to 50 years, and our analyses leads to the following conclusions. The QDO as a stand-alone signature is largely confined to the upper mesosphere. Its largest signature appears in the form of amplitude modulations of the QBO, AO and SAO, and these extend into the lower stratosphere. The downward control that characterizes the QBO apparently comes into play, and the longer time constants for diffusion and radiative loss at lower altitudes facilitate the QDO response. Although excited by the QBO, which is confined to low latitudes, the QDO is shown to extend to high latitudes. The effect is particularly large for the QBO with period around 33.5 month (near the upper limit of observations), which interacts with the SAO to produce a hemispherically symmetric QDO. Our analysis indicates that the QDO is transferred to high latitudes by the meridional circulation, which prominently exhibits this periodicity particularly in the amplitude modulation of the AO.

Mayr, H. G.↗

NASA NTRS2001

Quasi-Decadal Oscillations Generated by the QBO

Quasi-decadal oscillations (QDO) have been observed in the stratosphere and have been linked to the equatorial Quasi-Biennial Oscillation (QBO) and to the 11-year solar activity cycle. With the use of a 2D version of our Numerical Spectral Model (NSM) that incorporates Hines'Doppler Spread Parameterization (DSP) for gravity waves (GW), we demonstrate that beat periods between 9 and I I years can be generated by the QBO as it interacts through GW filtering with the Annual Oscillation (AO) and Semi-annual Oscillation (SAO). Results are discussed from computations covering up to 50 years, and our analyses leads to the following conclusions. The QDO as a stand-alone signature is largely confined to the upper mesosphere. Its largest signature appears in the form of amplitude modulations of the QBO, AO and SAO, and these extend into the lower stratosphere. The downward control that characterizes the QBO apparently comes into play, and the longer time constants for diffusion and radiative loss at lower altitudes facilitate the QDO response. Although excited by the QBO, which is confined to low latitudes, the QDO is shown to extend to high latitudes. The effect is particularly large for the QBO with period around 33.5 month (near the upper limit of observations), which interacts with the SAO to produce a hemispherically symmetric QDO. Our analysis indicates that the QDO is transferred to high latitudes by the meridional circulation, which prominently exhibits this periodicity particularly in the amplitude modulation of the AO. The propensity for beat periods around 10 years may help to increase significantly the efficiency for generating solar cycle signatures in the atmosphere, and we discuss numerical experiments related to that.

Mayr, Hans G.↗

DOE OSTI2023

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Role of Intra-Domain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: An Approach for Spatially Resolving Dynamics and Ion Transport

The design of safe and high-performance, nanostructured, block polymer (BP) electrolytes for lithium-ion batteries requires a thorough understanding of the key parameters that govern local structure and dynamics. Yet, the interfaces between microphase separated domains can introduce complexities in this local behavior that can be challenging to quantify. Herein, the local polymer, cation (Li + ), and anion dynamics were described in salt-doped polystyrene block poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) through a quantitative framework that considered the effects of polymer architecture, segmental mixing, chain stretching, and confinement on polymer mobility and ion transport. This framework was validated through nuclear magnetic resonance (NMR) spectroscopy measurements on solid (dry) polymer electrolyte samples. Notably, a mobility transition temperature (T mobility ) was identified through NMR spectroscopy that captured the local dynamics more accurately than the thermal glass transition temperature. Additionally, the approach quantitatively described the mobility gradient across a domain when segmental mixing effects were combined with chain stretching and confinement information, especially at higher segregation strengths – facilitating the assessment of local ion diffusion and conductivity. As a result, spatially averaged local ion diffusion predictions quantitatively matched NMR-measured ion diffusivities in the BP samples, while spatially summed ionic conductivity predictions across a domain qualitatively captured trends in the measured ionic conductivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Modeling the dynamics of phospholipids in the fluid phase of liposomes

Here, we present the derivation of a new model to describe neutron spin echo spectroscopy and quasi-elastic neutron scattering data on liposomes. We compare the new model with existing approaches and benchmark it with experimental data. The analysis indicates the importance of including all major contributions in the modeling of the intermediate scattering function. Simultaneous analysis of the experimental data on lipids with full contrast and tail contrast matched samples reveals highly confined lipid tail motion. A comparison of their dynamics demonstrates the statistical independence of tail-motion and height–height correlation of the membrane. A more detailed analysis indicates that the lipid tail relaxation is confined to a potential with cylindrical symmetry, in addition to the undulation and diffusive motion of the liposome. Despite substantial differences in the chemistry of the fatty acid tails, the observation indicates a universal behavior. The analysis of partially deuterated systems confirms the strong contribution of the lipid tail to the intermediate scattering function. Within the time range from 5 to 100 ns, the intermediate scattering function can be described by the height–height correlation function. The existence of the fast-localized tail motion and the contribution of slow translational diffusion of liposomes determine the intermediate scattering function for t < 5 ns and t > 100 ns, respectively. Taking into account the limited time window lowers the bending moduli by a factor of 1.3 (DOPC) to 2 (DMPC) compared to the full range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Production of warm ions in electron beam generated E × B plasma

Several recent experiments have demonstrated low-damage processing of 2D materials, such as graphene and single crystal diamond, using electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E × B) fields. Here, the low damage of these sensitive materials is commonly attributed to the low energy of ions incident to the substrate surface and the ion confinement in E x B fields. In this work, measurements of atom and ion velocity distribution functions in an e-beam E × B plasma at sub-mTorr argon pressures using a laser-induced fluorescence diagnostic revealed the presence of a warm population of ions with temperatures of ~ 1 eV that are sufficient to destroy the ion confinement in E × B fields and drive the ion flux by cross field diffusion in the direction opposite to the applied electric field, toward the plasma-bounded walls or substrate. Thus, it is this nearly ambipolar diffusion process that is responsible for the flux of charged particles impinging on the wall/substrate surface.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Reversible motions and disordered structure of soft particles in amorphous solids

In amorphous solids, soft vibrational modes derived from normal mode analysis can be utilized to identify the soft particles that are prone to irreversible rearrangements. However, the normal mode analysis approach cannot explain why the spatial distributions of clustered soft particles do not change over time. We define a softness parameter based on the vibrational density of states calculated directly from molecular dynamics simulations with both the harmonic vibrations and anharmonic relaxations being properly captured at finite temperatures. This parameter spontaneously correlates with the real space atomic motions and the dynamics heterogeneity. Using the softness parameter, we show that the softest particles are confined within rigid cages. These particles keep rearranging reversibly without long-range diffusion. The moderately soft particles rearrange irreversibly, and the hard particles mainly participate in vibrations without rearrangement. Here, we also show that the soft particles form locally disordered structures, while the hard particles present strong ordering. These findings confirm the defective nature of soft particles, and provide insights on the nature of softness as the ability to rearrange, but not necessarily irreversibly.

36 MATERIALS SCIENCE↗

NASA NTRS1995

ROSAT observations of the unusual supernova remnant CTB 80 containing the pulsar PSR 1951 + 32

The unusal supernova remnant CTB, 80 containing the 39.5 ms pulsar PSR 1951 + 32, has been observed with the Position Sensitive Proportional Counter (PSPC) and the High Resolution Imager (HRI) aboard ROSAT. The HRI image, centered on the pulsar, is composed of a bright compact core of approximately 1 arcminute radius containing the pulsar and a compact nebula, as well as a diffuse nebula extending approximately 5 arcminutes eastward of the pulsar. The PSPC allowed us to model the spectra of the point source, the compact nebula and the 5 arcminute diffuse nebula. For a power-law spectrum with photon index Gamma approximately 2 and an interstellar column density of N(sub H) approximately 3 x 10(exp 21)/cm(exp 21) the derived luminosities are approximately 2.3 x 10(exp 33) d(sub 2.5)(exp 2) ergs/s from the pointlike source, approximately 3.9 x 10(exp 33) d(sup 2.5)(exp 2) ergs/s from the compact nebula, and approximately 1.8 x 10(exp 33) d(sub 2.5)(exp 2) ergs/s from the 5 arcminutes diffuse nebula. In addition, the 2 deg diameter circular field of view of the PSPC reveals a hard emission feature southeast of the pulsar with a conical geometry extending out to the edge of the detector. The spectrum from this region is well described by a two-temperature Raymond-Smith thermal plasma with an average temperature of approximately 10(exp 7) K and a luminosity of approximately 10(exp 34)d(sub 2.5)(exp 2) ergs/s. Pulsations from the 39.5 ms puslar, PSR 1951 + 32, are detected at the 99% confidence level. The implied pulsed fration is approximately 35% with a complicated energy-dependent behavior. The compact core and the extended diffuse nebula can be explained as synchrotron radiation from the relativistic pulsar wind confined by the ram pressure of the surrounding inhomgeneous medium. The conelike feature detected southeast of PSR 1951 + 32 is consistent with emission from an optically thin SNR in the radiative cooling phase of its evolution.

Safi-Harb, Samar↗

NASA NTRS1983

The X-ray lobes of SS 433

Einstein Observatory imaging detector observations of SS 433 and its surrounding region are used to map an approximately 2 x 1-deg region containing W50 centered on SS 433, which exhibits two bright, diffuse lobes of X-ray emission symmetrically displaced east and west of SS 433 and aligned along the long axis of W50. The bulk of the emission from each lobe is confined to a region subtending an angle of about 30 deg at the central object. X-ray spectra indicate that the diffuse emission is considerably softer than SS 433 itself. The kinematic beam model for SS 433 is taken as the basis for a discussion of the properties and origins of the lobes and of their relationship to W50.

Watson, M. G.↗

DOE OSTI2022

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗

DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Polyelectrolyte Surface Diffusion in a Nanoslit Geometry

The surface diffusion of poly-l-lysine (PLL) in a planar nanoslit was studied using convex lens-induced confinement (CLiC) single-molecule tracking microscopy. Here, three surface chemistries were employed to understand the interplay of electrostatic and short-range interactions: an amine-functionalized silica surface, an oligo(ethylene oxide) (OEG)-modified surface, and a 1:1 mixture of the two ligands. Effective surface diffusion coefficients increased rapidly with slit height until saturating for slit heights <30 nm. While diffusion at a semi-infinite interface was significantly faster for OEG surfaces, the diffusion coefficient increased most rapidly with slit height for amine-functionalized surfaces, resulting in surface diffusion within very thin slits being nearly independent of surface chemistry. Intermittent random walks were simulated within a planar slit geometry, using experimentally measured parameters obtained from diffusion at a single interface to account for the characteristic short-range interactions between PLL and each surface chemistry, and were in good agreement with experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗

DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗

DOE OSTI2023

Nuclear magnetic resonance and molecular simulation study of H 2 and CH 4 adsorption onto shale and sandstone for hydrogen geological storage

Understanding pure H 2 and H 2 /CH 4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H 2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H 2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H 2 /CH 4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H 2 populations, i.e., free H 2 and adsorbed H 2 , that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H 2 NMR response similar to that of bulk H 2 . About 10 % of injected H 2 can be lost due to adsorption/desorption hysteresis in shale, and no H 2 loss (no hysteresis) is observed in sandstone. Here, our molecular simulation results support our NMR results that there are two H 2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH 4 outcompetes H 2 in adsorption onto kerogen, due to stronger CH 4 -kerogen interactions than H 2 -kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH 4 gas pressure, about ~30 % of residual CH 4 can be desorbed upon H 2 injection. The simulation results also predict that H 2 diffusion in porous kerogen is about one order of magnitude higher than that of CH 4 and CO 2 . This work provides an understanding of H 2 /CH 4 behaviors in deleted gas reservoirs upon H 2 injection and predictions of H 2 loss and CH 4 desorption in H 2 storage.

08 HYDROGEN↗