Search NASA⌕ Search

SEARCH · Search NASA

Results for “preferential 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗
DOE OSTI2021

Diffusion-mediated chemical concentration variation and void evolution in ion-irradiated NiCoFeCr high-entropy alloy

High-entropy alloys (HEAs) are proposed as potential structural materials for advanced nuclear systems, but little is known about the response of matrix chemistry in HEAs upon irradiation. Here, we reveal a substantial change of matrix chemical concentration as a function of irradiation damage (depth) in equiatomic NiCoFeCr HEA irradiated by 3 MeV Ni ions. After ion irradiation, the matrix contains more Fe/Cr in depth shallower than ~900–1000 nm but more Ni/Co from ~900–1000 nm to the end of the ion-damaged region due to the preferential diffusion of vacancies through Fe/Cr. Preferential diffusion also facilitates migration of vacancies from high radiation damage region to low radiation damage region, leading to no void formation below ~900–1000 nm and void formation around the end of the ion-damaged region at a fluence of 5 × 10 16 cm -2 (~123 dpa, displacements per atom, peak dose under full cascade mode). As voids grow significantly at an increased fluence (8 × 10 16 cm -2 , 196 dpa), the matrix concentration does not change dramatically due to new voids formed below ~900–1000 nm.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Diffusion-mediated chemical concentration variation and void evolution in ion-irradiated NiCoFeCr high-entropy alloy

High-entropy alloys (HEAs) are proposed as potential structural materials for advanced nuclear systems, but little is known about the response of matrix chemistry in HEAs upon irradiation. Here, we reveal a substantial change of matrix chemical concentration as a function of irradiation damage (depth) in equiatomic NiCoFeCr HEA irradiated by 3 MeV Ni ions. After ion irradiation, the matrix contains more Fe/Cr in depth shallower than ~900–1000 nm but more Ni/Co from ~900–1000 nm to the end of the ion-damaged region due to the preferential diffusion of vacancies through Fe/Cr. Preferential diffusion also facilitates migration of vacancies from high radiation damage region to low radiation damage region, leading to no void formation below ~900–1000 nm and void formation around the end of the ion-damaged region at a fluence of 5 × 10 16 cm –2 (~123 dpa, displacements per atom, peak dose under full cascade mode). As voids grow significantly at an increased fluence (8 × 10 16 cm –2 , 196 dpa), the matrix concentration does not change dramatically due to new voids formed below ~900–1000 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boost of the Bio-memristor Performance for Artificial Electronic Synapses by Surface Reconstruction

Biomaterial-based memristors (bio-memristors) are often adopted to emulate biological synapse functions and applied to construct neural computing networks in brain-inspired chip systems. However, the randomness of conductive filament formation in bio-memristors inhibits their switching performance by causing the dispersion of the device-switching parameters. In this case, a facile porous silk fibroin (p-SF) memristor was obtained through a protein surface reconstruction strategy, in which the size of the hole can be adjusted by the density of hybrid nanoseeds. Here, the porous SF memristors exhibit greatly enhanced electrical characteristics, including uniform I-V cycles, centralized distribution of the switching voltages, and both high and low resistances, compared to devices without pores. The results of three-dimensional (3D) simulations based on classical density functional theory (cDFT) suggest that the reconstructed pores in the SF layers guide the formation and fracture of Ag filaments under an electric field and enhance the overall conductivity by separating Ag + ion and electron diffusion pathways. Ag + ions are predicted to preferentially diffuse through pores, whereas electrons diffuse through the SF network. Interestingly, the device conductance can be bidirectionally modulated gradually by positive and negative voltages, can faithfully simulate short-term and long-term plasticity, and can even realize the triplet-spike-timing-dependent plasticity (triplet-STDP) rule, which can be used for pattern recognition in biological systems. The simulation results reveal that a memristor network of this type has an accuracy of ~95.78% in memory learning and the capability of pattern learning. This work provides a facile technology route to improve the performance of bionic-material memristors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

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

The effects of turbulence on the flame structure and NO formation of ammonia turbulent premixed combustion at various equivalence ratios

Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.

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

Effects of unburnt reaction progress on stretch flame dynamics under elevated temperatures

Flames in practical combustors and engines are inevitably subject to the effects of upstream chemical reaction progress and stretch, induced by elevated thermodynamic conditions and flow non-uniformities, respectively. Recent shock tube experiments and simulation studies on flame propagation have shown that flame speed under engine-relevant conditions can be enhanced with non-negligible upstream chemical reaction progress, especially when low-temperature heat release is involved in the unburnt mixture. On the other hand, depending on the mixture equivalence ratio and diluents, nonequidiffusion (including the non-unity Lewis number effect and the preferential diffusion effect) can couple with the flame stretch to fundamentally affect the flame propagation, which can manifest as either facilitation or suppression. Depending on the transport property of a reacting mixture, there hence can be either inhibition or promotion from unburnt reaction progress and stretch effects on flame propagation. In the current work, through one-dimensional numerical simulations of transient planar and spherical flames of n-heptane/air under elevated thermodynamic conditions, the combined effects of upstream chemical reaction progress and stretch on flame propagation are investigated. Results show that for both lean and rich n-heptane/air mixtures, flame speed can be substantially promoted with reaction progress, while the rich mixtures can exhibit opposite stretch dependence after first-stage ignition. Different definitions of effective Lewis number are adopted to explain the change in Markstein length for spherical flame in lean and rich mixtures reformed by the low temperature chemistry. Further, this work fills an important gap in laminar premixed flame research relevant to practical combustion systems and can provide useful insight into local turbulent flame behaviors and phi-sensitivity in engine combustion.

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

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗
DOE OSTI2022

An investigation of the failure modes in U-10Mo monolithic fuel irradiated to high burnup

Here, this study investigated possible failure mode sequences in U-10Mo monolithic fuel irradiated to very-high-burnup. The U-Mo fuel plate used in this study was characterized using scanning electron microscopy (SEM), wavelength dispersive spectroscopy (WDS), and image analysis techniques to investigate how microstructural features may initiate crack formation and propagation through the microstructure, potentially leading to blistering. Distinctly large fission gas pores (FGPs) were observed to preferentially align 5-15µm away from the U-10Mo/Zr interface, parallel to the Zr diffusion barrier. Preferential growth, alignment, and interconnection of FGP could initiate blistering in the fuel plate by creating a fission gas channel, parallel to the U-Mo/Zr interface. FGP alignment in the fuel phase near the U-Mo/Zr interface is believed to be one of the precursors to crack and Type 2 blister formation in monolithic U-Mo fuels. The chemical maps revealed differences in fission product behavior near the U-Mo/Zr interface such that Nd remains in the fuel matrix at the FGP sites. On the other hand, other fission products like Xe and Cs can get trapped in FGPs by Nd in the fuel phase and/or diffuse toward the Zr diffusion barrier.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Electrostatic Potential of Functional Cations as a Predictor of Hydroxide Diffusion Pathways in Nanoconfined Environments of Anion Exchange Membranes

Nanoconfined anion exchange membranes (AEMs) play a vital role in emerging electrochemical technologies. The ability to control dominant hydroxide diffusion pathways is an important goal in the design of nanoconfined AEMs. Such control can shorten hydroxide transport pathways between electrodes, reduce transport resistance, and enhance device performance. In this work, we propose an electrostatic potential (ESP) approach to explore the effect of the polymer electrolyte cation spacing on hydroxide diffusion pathways from a molecular perspective. By exploring cation ESP energy surfaces and validating outcomes through prior ab initio molecular dynamics simulations of nanoconfined AEMs, we find that we can achieve control over preferred hydroxide diffusion pathways by adjusting the cation spacing. The results presented in this work provide a unique and straightforward approach to predict preferential hydroxide diffusion pathways, enabling efficient design of highly conductive nanoconfined AEM materials for electrochemical technologies.

25 ENERGY STORAGE↗
DOE OSTI2020

Unusual electrical conductivity driven by localized stoichiometry modification at vertical epitaxial interfaces

Precise control of lattice mismatch accommodation and interfacial cation diffusion is critical to modulate correlated functionalities in epitaxial heterostructures, particularly when the interface composition is positioned near a compositional phase transition boundary. Here we select La 1-x Sr x MnO 3 (LSMO) as a prototype phase transition oxide and establish vertical epitaxial interfaces with NiO for exploring the strong interplay between strain accommodation, stoichiometry modification, and localized electron transport across the interface. It is found that localized stoichiometry modification overcomes the dead layer problem plaguing LSMO and leads to strongly directional conductivity, as manifested by a more than three orders of magnitude difference between out-of-plane and in-plane conductivity. Comprehensive atomic scale structural characterization and transport measurements reveal that this emerging behavior is created by a compositional change produced by preferential cation diffusion that pushes the LSMO phase transition from insulating into metallic within an ultrathin interface region. This study explores the nature of unusual electrical conductivity at vertical epitaxial interfaces and highlights the beneficial role of controllable cation diffusion that enables emerging functionalities for a broad range of potential applications such as memristors, spintronic devices, and novel nanoelectronic devices using strongly correlated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Novel effects of grain size and ion implantation on grain boundary segregation in ion irradiated austenitic steel

Nuclear reactor environments are extreme in nature due to the combination of exposure to corrosive coolant, mechanical stresses, and neutron irradiation damage. This combination of environments can lead to irradiation induced/assisted precipitation, segregation, and embrittlement coupled with enhanced creep and stress corrosion cracking. While nanostructured materials are of interest in a variety of applications due to their unique properties including high mechanical strength, for nuclear applications nanostructured materials have the added benefit of irradiation resistance due to the high volume density of sinks for irradiation-induced defects. In this study, we show unique radiation-induced segregation (RIS) behavior in nanocrystalline 304 stainless steel (SS304). Typical RIS in austenitic steels includes Ni enrichment and Cr depletion within the grain boundaries. In nanostructured SS304 this RIS is reduced. Furthermore, interestingly, some grain boundaries in the nanocrystalline SS304 in this study were found to be enriched in Cr after Fe 2+ ion irradiation at 500 °C up to 50 displacements per atom (dpa). To understand this unique observation, lattice-based atomic kinetic Monte Carlo simulations were performed to investigate the effects of grain size, thermal segregation, as well as injected interstitials (introduced as Fe self-ions). The results indicate that the Cr enrichment at grain boundaries may be caused by combined effects of grain size, self-ion injection and preferential Cr diffusion via the interstitial mechanism. In conclusion, the findings also suggest a possible difference in RIS between ion and neutron irradiations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diffusional lithium trapping as a failure mechanism of aluminum foil anodes in lithium-ion batteries

Aluminum foils are an appealing anode for lithium-ion batteries due to high capacity and low-cost, but their viability has been limited due to poor cyclability arising from pulverization and solid-electrolyte interphase growth. Here, we show that significant capacity degradation of aluminum foil anodes during electrochemical cycling also occurs due to diffusional lithium trapping. Scanning electron microscopy of cross-sectioned, cycled foils in the delithiated state reveals large regions of β-LiAl that are passivated by a surface layer of α-Al, which has poor Li + diffusivity. It is found that lithium diffusion occurs preferentially along the β-LiAl grain boundaries, so the grain structure after initial lithiation significantly affects the trapping behavior. Diffusional lithium trapping is exacerbated by both higher delithiation rates and higher areal capacity, presenting a challenge towards commercialization of aluminum foil anodes. We further demonstrate that diffusional trapping in aluminum foil anodes can be mitigated through alloy design, with the addition of 2–3 wt% Li yielding improved first cycle efficiency, and the addition of 1 wt% Si yielding improved cycle life. These results provide a mechanistic understanding of diffusional lithium trapping in aluminum foil anodes and highlight compositional design of alloys as a promising strategy to overcome it.

25 ENERGY STORAGE↗
DOE OSTI2020

Visualizing Anisotropic Oxygen Diffusion in Ceria under Activated Conditions

Oxygen reactivity plays a key role in the performance of ceria-based catalysts. Aberration-corrected transmission electron microscopy and molecular dynamics simulations were used to study the oxygen atom diffusion in ceria under activated conditions. Reactive oxygen atom and its real-time diffusion were visualized. Additionally, the interplay between cerium and oxygen atoms originating from a Coulomb interaction was revealed by the out-of-plane buckling of cerium atoms associated with oxygen transport. Anisotropic oxygen atom diffusion that depends on crystal orientations was discovered, demonstrating a preferential [001] crystallographic diffusion pathway. These findings reveal prospects for applications of anisotropic orientation-relevant fluorite-structured oxides.

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

Understanding the effect of minor alloying elements on helium bubble formation in ferritic-martensitic steels

Ferritic-martensitic steels are promising structural materials for advanced nuclear reactors. To minimize long-term radioactivity, reduced-activation ferritic-martensitic steels have been developed by substituting high-activation elements like Ni and Mo with low-activation elements such as W. However, the impact of these alloying modifications on helium bubble formation, which plays a key role in material swelling, remains unclear. Here, in this study, we compared helium bubble formation in ferritic-martensitic steel T91 and reduced-activation ferritic-martensitic steel F82H. Both materials were irradiated with sequential 100 keV, 150 keV, and 200 keV helium ions to a dose of 0.5 dpa and a helium concentration of 9,000 appm at 500°C. The helium bubbles in F82H exhibited a larger average size and a lower density than those in T91, suggesting differences in minor alloying elements may influence the bubble growth. Here, to investigate the effects of these alloying elements, we characterized radiation-induced segregation near bubbles and grain boundaries. Prominent Ni-Mn-Si enriched clusters were found near bubbles in T91, while only Mn-Si enriched clusters were found near bubbles in F82H. In addition, the obvious Cr enrichment near grain boundaries was absent around bubbles in both steels. The different segregation trends among elements revealed the variations in element diffusion mechanisms and the different sink biases between bubbles and grain boundaries. Cr enrichment near grain boundaries is mostly driven by interstitial-mediated diffusion. However, since bubble growth relies on net vacancy flux, vacancy-mediated diffusion plays a dominant role in controlling element segregation near bubbles. Therefore, Cr enrichment was not found near bubbles. Because of preferential vacancy-drag diffusion for Ni, Si and Mn, these elements were enriched near bubbles. Due to the strong binding energies of vacancies with these solute atoms, the vacancy diffusivity can be reduced near these solutes. Therefore, the more prominent Ni-Si-Mn clustered near helium bubbles in T91 lead to stronger suppression of helium bubble growth compared to F82H.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomistic insights into the nucleation and growth of platinum on palladium nanocrystals

Abstract Despite the large number of reports on colloidal nanocrystals, very little is known about the mechanistic details in terms of nucleation and growth at the atomistic level. Taking bimetallic core-shell nanocrystals as an example, here we integrate in situ liquid-cell transmission electron microscopy with first-principles calculations to shed light on the atomistic details involved in the nucleation and growth of Pt on Pd cubic seeds. We elucidate the roles played by key synthesis parameters, including capping agent and precursor concentration, in controlling the nucleation site, diffusion path, and growth pattern of the Pt atoms. When the faces of a cubic seed are capped by Br − , Pt atoms preferentially nucleate from corners and then diffuse to edges and faces for the creation of a uniform shell. The diffusion does not occur until the Pt deposited at the corner has reached a threshold thickness. At a high concentration of the precursor, self-nucleation takes place and the Pt clusters then randomly attach to the surface of a seed for the formation of a non-uniform shell. These atomistic insights offer a general guideline for the rational synthesis of nanocrystals with diverse compositions, structures, shapes, and related properties.

36 MATERIALS SCIENCE↗
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↗