Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIFFUSION”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Static and Dynamic Effects of Lateral Carrier Diffusion in Semiconductor Lasers

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. It is well known that the value of diffusion coefficients affects the threshold pumping current of a semiconductor laser. At the same time, the strength of carrier diffusion process is expected to affect the modulation bandwidth of an AC-modulated laser. It is important not only to investigate the combined DC and AC effects due to carrier diffusion, but also to separate the AC effects from that of the combined effects in order to provide design insights for high speed modulation. In this presentation, we apply a hydrodynamic model developed by the present authors recently from the semiconductor Bloch equations. The model allows microscopic calculation of the lateral carrier diffusion coefficient, which is a nonlinear function of the carrier density and plasma temperature. We first studied combined AC and DC effects of lateral carrier diffusion by studying the bandwidth dependence on diffusion coefficient at a given DC current under small signal modulation. The results show an increase of modulation bandwidth with decrease in the diffusion coefficient. We simultaneously studied the effects of nonlinearity in the diffusion coefficient. To clearly identify how much of the bandwidth increase is a result of decrease in the threshold pumping current for smaller diffusion coefficient, thus an effective increase of DC pumping, we study the bandwidth dependence on diffusion coefficient at a given relative pumping. A detailed comparison of the two cases will be presented.

Li, Jian-Zhong↗
NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗
NASA NTRS2016

Passive Rocket Diffuser Testing: Reacting Flow Performance of Four Second-Throat Geometries

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure. As one of the nation's largest rocket testing facilities, the performance and design limitations of diffusers are of great interest to NASA's Stennis Space Center. This paper describes a series of tests conducted on four diffuser configurations to better understand the effects of inlet geometry and throat area on starting behavior and boundary layer separation. The diffusers were tested for a duration of five seconds with a 1455-pound thrust, LO2/GH2 thruster to ensure they each reached aerodynamic steady state. The effects of a water spray ring at the diffuser exits and a water-cooled deflector plate were also evaluated. Static pressure and temperature measurements were taken at multiple axial locations along the diffusers, and Computational Fluid Dynamics (CFD) simulations were used as a tool to aid in the interpretation of data. The hot combustion products were confirmed to enable the diffuser start condition with tighter second throats than predicted by historical cold-flow data or the theoretical normal shock method. Both aerodynamic performance and heat transfer were found to increase with smaller diffuser throats. Spray ring and deflector cooling water had negligible impacts on diffuser boundary layer separation. CFD was found to accurately capture diffuser shock structures and full-flowing diffuser wall pressures, and the qualitative behavior of heat transfer. However, the ability to predict boundary layer separated flows was not consistent.

Jones, Daniel R.↗
DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Interdiffusion and Intrinsic Diffusion in the Mg-Al System

Solid-to-solid diffusion couples were assembled and annealed to examine the diffusion between pure Mg (99.96%) and Al (99.999%). Diffusion anneals were carried out at 300 , 350 , and 400 C for 720, 360, and 240 hours, respectively. Optical and scanning electron microscopes were utilized to identify the formation of the intermetallic phases, -Al12Mg17 and -Al3Mg2 and absence of the -phase in the diffusion couples. Thicknesses of the -Al12Mg17 and -Al3Mg2 phases were measured and the parabolic growth constants were calculated to determine the activation energies for the growth, 165 and 86 KJ/mole, respectively. Concentration profiles were determined with electron microprobe analysis using pure elemental standards. Composition-dependent interdiffusion coefficients in Mg-solid solution, -Al12Mg17 and - Al3Mg2 and Al-solid solutions were calculated based on the Boltzmann-Matano analysis. Average effective interdiffusion coefficients for each phase were also calculated, and the magnitude was the highest for the -Al3Mg2 phase, followed by -Al12Mg17, Al-solid solution and Mg-solid solution. Intrinsic diffusion coefficients based on Huemann s analysis (e.g., marker plane) were determined for the ~38 at.% Mg in the -Al3Mg2 phase. Activation energies and the pre-exponential factors for the inter- and intrinsic diffusion coefficients were calculated for the temperature range examined. The -Al3Mg2 phase was found to have the lowest activation energies for growth and interdiffusion among all four phases studied. At the marker location in the -Al3Mg2 phase, the intrinsic diffusion of Al was found to be faster than that of Mg. Extrapolations of the impurity diffusion coefficients in the terminal solid solutions were made and compared to the available self- and impurity diffusion data from literature. Thermodynamic factor, tracer diffusion coefficients and atomic mobilities at the marker plane composition were approximated using available literature values of Mg activity in the -Al3Mg2 phase.

Boltzmann-Matano↗
DOE OSTI2017

Model-based analysis of the impact of diffuse radiation on CO 2 exchange in a temperate deciduous forest

Clouds and aerosols increase the fraction of global solar irradiance that is diffuse light. This phenomenon is known to increase the photosynthetic light use efficiency (LUE) of closed-canopy vegetation by redistributing photosynthetic photon flux density (400–700 nm) from saturated, sunlit leaves at the top of the canopy, to shaded leaves deeper in the canopy. We combined a process-based carbon cycle model with 10 years of eddy covariance carbon flux measurements and other ancillary data sets to assess 1) how this LUE enhancement influences interannual variation in carbon uptake, and 2) how errors in modeling diffuse fraction affect predictions of carbon uptake. Modeled annual gross primary productivity (GPP) increased by ≈0.94% when observed levels of diffuse fraction were increased by 0.01 (holding total irradiance constant). The sensitivity of GPP to increases in diffuse fraction was highest when the diffuse fraction was low to begin with, and lowest when the diffuse fraction was already high. Diffuse fraction also explained significantly more of the interannual variability of modeled net ecosystem exchange (NEE), than did total irradiance. Two tested radiation partitioning models yielded over- and underestimates of diffuse fraction at our site, which propagated to over- and underestimates of annual NEE, respectively. Our findings highlight the importance of incorporating LUE enhancement under diffuse light into models of global primary production, and improving models of diffuse fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Non-monotonic behavior of the lateral diffusivity in an adsorbate as a function of the surface coverage

Diffusion dynamics of adsorbate particles confined on a surface strongly depends on the adsorbate surface coverage. Our molecular dynamics simulation reveals that the average lateral diffusivity of the Lennard-Jones adsorbate on a flat substrate varies non-monotonically with the surface coverage. Given that the adsorbate becomes stratified into distinct layers, we propose a novel approach to analyze the lateral diffusion of the adsorbate particles in different layers, via dividing the total lateral diffusion into consecutive lateral diffusion processes, separated by the layer boundary crossing events. A new definition of the lateral diffusivity in the specific layer has been introduced, and it is revealed that the total lateral diffusivity is equal to the sum of lateral diffusivities in different layers. Furthermore, the non-monotonic behavior of the average lateral diffusivity with the surface coverage is found to result from the difference in behavior of the lateral diffusivity between the first and other layers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analyzing the co-evolution of green technology diffusion and consumers’ pro-environmental attitudes: An agent-based model

Massive diffusion of green technologies is significant for building a cleaner world. However, the process of technology diffusion is usually slow and complex. An in-depth understanding of the mechanism regarding green technology diffusion is an essential precondition for effectively stimulating this process. Green technology diffusion heavily involves both social and technological changes. Although existing studies have provided rich knowledge about identifying critical drivers and barriers affecting green technology diffusion, research that considers consumers’ pro-environmental attitudes and green technology diffusion as an evolving system is still sparse. Aiming at exploring the co-evolution of consumers’ pro-environmental attitudes and green technology diffusion, this paper builds an agent-based model that integrates the relative agreement model with technology diffusion theories to conduct a sequence of controlled numerical experiments, which progressively unveil how attitudinal and technological factors impact green technology diffusion. The main findings include that (1) improving consumers’ pro-environmental attitudes is prominently beneficial to green technology diffusion and maturation; (2) technology maturity has very limited impact on consumers’ first-time purchases but significantly affects consumers’ satisfaction, which would further impact consumers’ repeat purchases; (3) consumers that do not support green technologies frequently emerge during the evolution of attitudes (despite the high technology maturity), which corresponds to the emergence of the anti-environmental groups observed in the real world; (4) active interactions between non-adopters enable their attitudes to converge, which results in the polarization of consumers’ attitudes.

ABM↗
DOE OSTI2014

Compositions and morphologies of TiAlSi intermetallics in different diffusion couples

Two kinds of diffusion couples were designed to investigate the formation of ternary TiAlSi phases in Al–Si–Ti alloys. It was found that different diffusion processes result in various compositions and morphologies of TiAlSi intermetallics. The melted Al, Si and Ti atoms in the diffusion couple leads to the formation of flake-like TiAlSi phase through liquid–liquid reaction. Besides, unidirectional diffusion of Al and Si atoms into blocky TiAl{sub 3} particles or Ti powders via a liquid–solid diffusion process also results in the formation of TiAlSi, while keeping the block-like morphology. This kind of diffusion is a gradual process, driven by the concentration gradient. The reactions in the diffusion couples are helpful to understand the compositional and morphological evolutions of TiAlSi as reported in previous work. - Highlights: • Two diffusion couples were designed to investigate the formation of TiAlSi phases. • Compositions and morphologies of TiAlSi are influenced by the diffusion process. • Liquid–liquid and liquid–solid diffusions were detected. • The corresponding mechanisms were discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Banded structures in electron pitch angle diffusion coefficients from resonant wave-particle interactions

Electron pitch angle (D{sub αα}) and momentum (D{sub pp}) diffusion coefficients have been calculated due to resonant interactions with electrostatic electron cyclotron harmonic (ECH) and whistler mode chorus waves. Calculations have been performed at two spatial locations L = 4.6 and 6.8 for electron energies ≤10 keV. Landau (n = 0) resonance and cyclotron harmonic resonances n = ±1, ±2, … ±5 have been included in the calculations. It is found that diffusion coefficient versus pitch angle (α) profiles show large dips and oscillations or banded structures. The structures are more pronounced for ECH and lower band chorus (LBC) and particularly at location 4.6. Calculations of diffusion coefficients have also been performed for individual resonances. It is noticed that the main contribution of ECH waves in pitch angle diffusion coefficient is due to resonances n = +1 and n = +2. A major contribution to momentum diffusion coefficients appears from n = +2. However, the banded structures in D{sub αα} and D{sub pp} coefficients appear only in the profile of diffusion coefficients for n = +2. The contribution of other resonances to diffusion coefficients is found to be, in general, quite small or even negligible. For LBC and upper band chorus waves, the banded structures appear only in Landau resonance. The D{sub pp} diffusion coefficient for ECH waves is one to two orders smaller than D{sub αα} coefficients. For chorus waves, D{sub pp} coefficients are about an order of magnitude smaller than D{sub αα} coefficients for the case n ≠ 0. In case of Landau resonance, the values of D{sub pp} coefficient are generally larger than the values of D{sub αα} coefficients particularly at lower energies. As an aid to the interpretation of results, we have also determined the resonant frequencies. For ECH waves, resonant frequencies have been estimated for wave normal angle 89° and harmonic resonances n = +1, +2, and +3, whereas for whistler mode waves, the frequencies have been calculated for angle 10° and Landau resonance. Further, in ECH waves, the banded structures appear for electron energies ≥1 keV, and for whistler mode chorus waves, structures appear for energies >2 keV at L = 4.6 and above 200 eV for L = 6.8. The results obtained in the present work will be helpful in the study of diffusion curves and will have important consequences for diffuse aurora and pancake distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reproducibility of protein x-ray diffuse scattering and potential utility for modeling atomic displacement parameters

Protein structure and dynamics can be probed using x-ray crystallography. Whereas the Bragg peaks are only sensitive to the average unit-cell electron density, the signal between the Bragg peaks—diffuse scattering—is sensitive to spatial correlations in electron-density variations. Although diffuse scattering contains valuable information about protein dynamics, the diffuse signal is more difficult to isolate from the background compared to the Bragg signal, and the reproducibility of diffuse signal is not yet well understood. We present a systematic study of the reproducibility of diffuse scattering from isocyanide hydratase in three different protein forms. Both replicate diffuse datasets and datasets obtained from different mutants were similar in pairwise comparisons (Pearson correlation coefficient ≥0.8). The data were processed in a manner inspired by previously published methods using custom software with modular design, enabling us to perform an analysis of various data processing choices to determine how to obtain the highest quality data as assessed using unbiased measures of symmetry and reproducibility. The diffuse data were then used to characterize atomic mobility using a liquid-like motions (LLM) model. This characterization was able to discriminate between distinct anisotropic atomic displacement parameter (ADP) models arising from different anisotropic scaling choices that agreed comparably with the Bragg data. Our results emphasize the importance of data reproducibility as a model-free measure of diffuse data quality, illustrate the ability of LLM analysis of diffuse scattering to select among alternative ADP models, and offer insights into the design of successful diffuse scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Diffusion of acceptor dopants in monoclinic 𝛽−Ga 2⁢ O 3

𝛽−Ga 2 ⁢O 3 is a promising material for next-generation power electronics because of its ultrawide band gap and high critical breakdown voltage. However, realizing its full potential requires precise control over dopant incorporation and stability. In this work, we use first-principles calculations to systematically assess the diffusion behavior of eight potential deep-level substitutional acceptors (Au, Ca, Co, Cu, Fe, Mg, Mn, and Ni) in 𝛽−Ga 2 ⁢O 3 . We consider two key diffusion mechanisms: (i) interstitial diffusion under nonequilibrium conditions relevant to ion implantation, and (ii) trap-limited diffusion (TLD) under near-equilibrium thermal annealing conditions. Our results reveal a strong diffusion anisotropy along the 𝑏 and 𝑐 axes, with dopant behavior governed by competition between diffusion and incorporation (or dissociation) activation energies. Under interstitial diffusion, Ca$^{2+}_{i}$ and Mg$^{2+}_{i}$ show the most favorable combination of low migration and incorporation barriers, making them promising candidates for efficient doping along the 𝑏 and 𝑐 axes, respectively. In contrast, Au$^{+}_{i}$ diffuses readily, but exhibits an incorporation barrier that exceeds 5 eV, rendering it ineffective as a dopant. From a thermal stability perspective, Co$^{2+}_{i}$ shows poor activation but high diffusion barriers, which may suppress undesirable migration at elevated temperatures. Under trap-limited diffusion, the dissociation of dopant-host complexes controls mobility. Mg$^{2+}_{i}$ again emerges as a leading candidate, exhibiting the lowest dissociation barriers along both axes, whereas Co$^{2+}_{i}$ and Fe$^{2+}_{i}$ display the highest barriers, suggesting improved dopant retention under thermal stress. In conclusion, our findings guide dopant selection by balancing activation and thermal stability, essential for robust semi-insulating substrates.

Defects↗
NASA NTRS2016

Banded Structures in Electron Pitch Angle Diffusion Coefficients from Resonant Wave Particle Interactions

Electron pitch angle (D (alpha)) and momentum (D(pp)) diffusion coefficients have been calculated due to resonant interactions with electrostatic electron cyclotron harmonic (ECH) and whistler mode chorus waves. Calculations have been performed at two spatial locations L = 4.6 and 6.8 for electron energies 10 keV. Landau (n = 0) resonance and cyclotron harmonic resonances n = +/-1, +/-2,...+/-5 have been included in the calculations. It is found that diffusion coefficient versus pitch angle (alpha) profiles show large dips and oscillations or banded structures. The structures are more pronounced for ECH and lower band chorus (LBC) and particularly at location 4.6. Calculations of diffusion coefficients have also been performed for individual resonances. It is noticed that the main contribution of ECH waves in pitch angle diffusion coefficient is due to resonances n = +1 and n = +2. A major contribution to momentum diffusion coefficients appears from n = +2. However, the banded structures in D alpha and Dpp coefficients appear only in the profile of diffusion coefficients for n = +2. The contribution of other resonances to diffusion coefficients is found to be, in general, quite small or even negligible. For LBC and upper band chorus waves, the banded structures appear only in Landau resonance. The Dpp diffusion coefficient for ECH waves is one to two orders smaller than D alpha coefficients. For chorus waves, Dpp coefficients are about an order of magnitude smaller than D alpha coefficients for the case n does not = 0. In case of Landau resonance, the values of Dpp coefficient are generally larger than the values of D alpha coefficients particularly at lower energies. As an aid to the interpretation of results, we have also determined the resonant frequencies. For ECH waves, resonant frequencies have been estimated for wave normal angle 89 deg and harmonic resonances n = +1, +2, and +3, whereas for whistler mode waves, the frequencies have been calculated for angle 10 deg and Landau resonance. Further, in ECH waves, the banded structures appear for electron energies (is) greater than1 keV, and for whistler mode chorus waves, structures appear for energies greater than 2 keV at L = 4.6 and above 200 eV for L = 6.8. The results obtained in the present work will be helpful in the study of diffusion curves and will have important consequences for diffuse aurora and pancake distributions.

diffusion↗
DOE OSTI2015

CO DIFFUSION INTO AMORPHOUS H{sub 2}O ICES

The mobility of atoms, molecules, and radicals in icy grain mantles regulates ice restructuring, desorption, and chemistry in astrophysical environments. Interstellar ices are dominated by H{sub 2}O, and diffusion on external and internal (pore) surfaces of H{sub 2}O-rich ices is therefore a key process to constrain. This study aims to quantify the diffusion kinetics and barrier of the abundant ice constituent CO into H{sub 2}O-dominated ices at low temperatures (15–23 K), by measuring the mixing rate of initially layered H{sub 2}O(:CO{sub 2})/CO ices. The mixed fraction of CO as a function of time is determined by monitoring the shape of the infrared CO stretching band. Mixing is observed at all investigated temperatures on minute timescales and can be ascribed to CO diffusion in H{sub 2}O ice pores. The diffusion coefficient and final mixed fraction depend on ice temperature, porosity, thickness, and composition. The experiments are analyzed by applying Fick’s diffusion equation under the assumption that mixing is due to CO diffusion into an immobile H{sub 2}O ice. The extracted energy barrier for CO diffusion into amorphous H{sub 2}O ice is ∼160 K. This is effectively a surface diffusion barrier. The derived barrier is low compared to current surface diffusion barriers in use in astrochemical models. Its adoption may significantly change the expected timescales for different ice processes in interstellar environments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Tritium species diffusion on and desorption from γ-LiAlO 2 (100) surface: A first-principles investigation

γ-LiAlO 2 enriched in the 6 Li isotope is a good candidate for tritium (T, 3 H) production in nuclear reactor. In this work, to better understand its 3 H production performance, first-principles calculations are performed to study the 3 H-related species diffusion on and desorption from the γ-LiAlO 2 (100) surface. For the diffusion process, we investigate the pathways and energy barriers for the T and OT diffusion on the pristine and defective γ-LiAlO 2 (100) surfaces. Our results reveal that the O and Li/Al atoms of the surface layer play important roles for T and OT diffusion, respectively. By comparing with its diffusion barrier in bulk γ-LiAlO 2 , the order of energy barrier for T diffusion is: bulk < surface < defective surface, which is reversed for OT diffusion. Space is the prime factor for OT diffusion, while diffusion-mediating atom instead of space becomes the key factor for T diffusion. For the desorption process, we investigate the 3 H-related species (T, OT, T 2 , and T 2 O) desorption from the γ-LiAlO 2 (100) surface. The desorption energies we obtained are so high that the desorption behavior of all the 3 H-related species is nearly prohibited. We expect the possible T 2 or T 2 O release from the large accumulation of 3 H atoms on the surface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗
DOE OSTI2015

The effect of shear flow on the rotational diffusion of a single axisymmetric particle

Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$, as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg$ 1. For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$.

42 ENGINEERING↗
DOE OSTI2016

Molecular dynamics simulation of diffusion and electrical conductivity in montmorillonite interlayers

In this study, the diffusion of water and ions in the interlayer region of smectite clay minerals represents a direct probe of the type and strength of clay–fluid interactions. Interlayer diffusion also represents an important link between molecular simulation and macroscopic experiments. Here we use molecular dynamics simulation to investigate trends in cation and water diffusion in montmorillonite interlayers, looking specifically at the effects of layer charge, interlayer cation and cation charge (sodium or calcium), water content, and temperature. For Na-montmorillonite, the largest increase in ion and water diffusion coefficients occurs between the one-layer and two-layer hydrates, corresponding to the transition from inner-sphere to outer-sphere surface complexes. Calculated activation energies for ion and water diffusion in Na-montmorillonite are similar to each other and to the water hydrogen bond energy, suggesting the breaking of water–water and water–clay hydrogen bonds as a likely mechanism for interlayer diffusion. A comparison of interlayer diffusion with that of bulk electrolyte solutions reveals a clear trend of decreasing diffusion coefficient with increasing electrolyte concentration, and in most cases the interlayer diffusion results are nearly coincident with the corresponding bulk solutions. Trends in electrical conductivities computed from the ion diffusion coefficients are also compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗