Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Status of Axisymmetric CFD of an Eleven Inch Diameter Hybrid Rocket Motor

Current status of a steady state, axisymmetric analysis of an experimental 11 inch diameter hybrid rocket motor internal flow field is given. The objective of this effort is to develop a steady state axisymmetric model of the 11 inch hybrid rocket motor which can be used as a design and/or analytical tool. A test hardware description, modeling approach, and future plans are given. The analysis was performed with FDNS implementing several finite rate chemistry sets. A converged solution for a two equation and five species set on a 'fine' grid is shown.

Ruf, Joseph↗
NASA NTRS2005

Afterbody Heating Predictions for a Mars Science Laboratory Entry Vehicle

The Mars Science Laboratory mission intends to deliver a large rover to the Martian surface within 10 km of its target site. One candidate entry vehicle aeroshell consists of a 3.75-m diameter, 70-deg sphere-cone forebody and a biconic afterbody similar to that of Viking. This paper presents computational fluid dynamics predictions of laminar afterbody heating rates for this configuration and a 2010 arrival at Mars. Computational solutions at flight conditions used an 8-species Mars gas model in chemical and thermal non-equilibrium. A grid resolution study examined the effects of mesh spacing on afterbody heating rates and resulted in grids used for heating predictions on a reference entry trajectory. Afterbody heating rate reaches its maximum value near 0.6 W/sq cm on the first windward afterbody cone at the time of peak freestream dynamic pressure. Predicted afterbody heating rates generally are below 3% of the forebody laminar nose cap heating rate throughout the design trajectory. The heating rates integrated over time provide total heat load during entry, which drives thermal protection material thickness.

Edquist, Karl T.↗
NASA NTRS2007

The ISS Water Processor Catalytic Reactor as a Post Processor for Advanced Water Reclamation Systems

Advanced water processors being developed for NASA s Exploration Initiative rely on phase change technologies and/or biological processes as the primary means of water reclamation. As a result of the phase change, volatile compounds will also be transported into the distillate product stream. The catalytic reactor assembly used in the International Space Station (ISS) water processor assembly, referred to as Volatile Removal Assembly (VRA), has demonstrated high efficiency oxidation of many of these volatile contaminants, such as low molecular weight alcohols and acetic acid, and is considered a viable post treatment system for all advanced water processors. To support this investigation, two ersatz solutions were defined to be used for further evaluation of the VRA. The first solution was developed as part of an internal research and development project at Hamilton Sundstrand (HS) and is based primarily on ISS experience related to the development of the VRA. The second ersatz solution was defined by NASA in support of a study contract to Hamilton Sundstrand to evaluate the VRA as a potential post processor for the Cascade Distillation system being developed by Honeywell. This second ersatz solution contains several low molecular weight alcohols, organic acids, and several inorganic species. A range of residence times, oxygen concentrations and operating temperatures have been studied with both ersatz solutions to provide addition performance capability of the VRA catalyst.

Nalette, Tim↗
DOE OSTI2020

Mercury Speciation During Vitrification of LAW

n this work, laboratory and engineering-scale tests were conducted that mimic the reactive environment of the Hanford Tank Waste Treatment and Immobilization Plant (WTP) low activity waste (LAW) melter and off-gas system in order to assess the speciation of mercury at select points in the processing system. The experimental protocols, test equipment, feed materials used, and rationale for testing are detailed in the Test Plan for this work. Particular attention was paid to the amount and speciation of mercury in submerged bed scrubber (SBS) solutions, which are intended to be recycled back to the melter feed at the WTP. Results are presented from detailed mercury analysis of these solutions using Environmental Protection Agency (EPA) Method 1630 and modifications to this method for dimethyl mercury and methyl mercury as well as EPA Method 1631 and modifications to this method for total, dissolved, elemental, suspended, and ionic mercury. Testing was conducted in two separate sets of experiments: (1) One set of tests (crucible-scale furnace tests) that involved heating small batches of mercury-spiked melter feed in crucibles to determine the amount of mercury retained in the glass and the species of mercury in the exhaust gases; (2) A second set of tests (DM10 melter tests) was planned that involved creating LAW melter plenum gas compositions using the DM10 melter system, injecting mercury into the off-gas stream, and passing that stream though a reactor that simulates various plenum gas conditions. Operational issues led to the need to use one of the VSL DM100 melters in place of the planned DM10 melter for these tests. The speciation of mercury after exposure to those conditions was monitored. In both sets of tests, the exhaust gases were run through a scrubber that was intended to mimic the LAW SBS in order to determine how much of the mercury exiting the melter would be retained in the primary off-gas system fluids and in what form. After passing through the SBS, the exhaust stream was analyzed to determine particulate, ionic, elemental, and total mercury passing downstream of the SBS. In tests employing gases derived from the DM100 melter that were spiked with elemental mercury, the processing system provided gas temperatures and residence times that are representative of the WTP LAW vitrification system in order to assess the effect of those conditions on mercury speciation. The Decontamination Factor (DF) across the system in the tests with mercury-spiked DM100 melter exhaust was determined using the analytical data from EPA Method 30B (Fluegas Adsorption Mercury Speciation (FAMS TM )) exhaust samples and the amount of mercury detected in the SBS solutions. Mercury species used in melter feed crucible scale tests were divalent (chloride and iodide), monovalent (chloride and fluoride), and elemental mercury. Individual tests included only a single form of mercury in the feed. The results for total mercury mass balances in the crucible tests are also presented.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ionic Dopant-Induced Ordering Enhances the Thermoelectric Properties of a Polythiophene-Based Block Copolymer

Conjugated polymer-based block copolymers (CP-BCPs) are an unexplored class of materials for organic thermoelectrics. In this study, the authors report on the electronic conductivity (σ) and Seebeck coefficient (α) of a newly synthesized CP-BCP, poly(3-hexylthiophene)-block-poly (oligo-oxyethylene methacrylate) (P3HT-b-POEM), upon solution co-processing with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), and subsequently vapor-doping with a molecular dopant, 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ). It is found that the addition of the hydrophilic block POEM greatly enhances the processability of P3HT, enabling homogeneous solution-mixing with LiTFSI. Notably, interactions between P3HT-b-POEM with ionic species significantly improve molecular order and unexpectedly cause electrical oxidizing doping of P3HT block both in solution and solid-states, a phenomenon that has not been previously observed in Li-salt containing P3HT. Vapor doping of P3HT-b-POEM-LiTFSI thin films with F4TCNQ further enhances σ and yields a thermoelectric power factor PF = α 2 σ of 13.0 μW m -1 K -2 , which is more than 20 times higher than salt-free P3HT-b-POEM sample. Through modeling thermoelectric behaviors of P3HT-b-POEM with the Kang-Snyder transport model, the improvement in PF is attributed to higher electronic charge mobility originating from the enhanced molecular ordering of P3HT. The results demonstrate that solution co-processing CP-BCPs with a salt is a powerful method to control structure and performance of organic thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of cation species on thermal expansion of Y 2 Si 2 O 7 –Gd 2 Si 2 O 7 solid solutions

Mixtures of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 were synthesized by solid-state reaction at 1600°C and characterized via in situ x-ray diffraction (XRD) to determine their coefficients of thermal expansion (CTE). All solid solutions within the system exhibited the orthorhombic δ-RE 2 Si 2 O 7 (Pna2 1 ) structure. Thermal expansion measurements of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 correlated well with reported values in literature, and all synthesized solid solutions exhibited CTEs between Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 . Generally, there was a slight decrease in CTE exhibited by the materials with increasing Gd 2 Si 2 O 7 content, with Gd 2 Si 2 O 7 having the lowest CTEs and Y 2 Si 2 O 7 the highest CTEs. Here, the decrease in CTE was attributed to stronger bonds of Gd-O over Y-O, as determined by calculated crystal orbital Hamilton populations using density functional theory. However, such differences were very small and crystal structure was the dominating factor in CTE trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

A pressure correction method for the calculation of compressible chemical reacting flows

A recently developed noniterative method for the solution of the transient fluid flow equations at all speed is extended to handle chemical reacting flows. The species conservation equations are loosely coupled into the predictor/multicorrector sequence of the solution procedure. A split-operator method separates the chemical kinetics terms from the fluid-dynamical terms, as well as an implicit differencing method enhance the numerical stability. The method was applied for turbulent diffusion flame calculations and for the analyses of high pressure, axisymmetric turbulent hypersonic nozzle flows. The diffusion flame results were compared with a similar pressure method for fast chemistry integration scheme without operator-splitting. Simulations of the nozzle flow indicated that the nonideal intermolecular effects must be included in the analysis and design of high pressure hypersonic nozzle.

Chen, Z. J.↗
DOE OSTI2022

C−H Activation of Inert Arenes using a Photochemically Activated Guanidinato‐Magnesium(I) Compound

Abstract UV irradiation of solutions of a guanidinate coordinated dimagnesium(I) compound, [{(Priso)Mg} 2 ] 3 (Priso=[(DipN) 2 CNPr i 2 ] − , Dip=2,6‐diisopropylphenyl), in either benzene, toluene, the three isomers of xylene, or mesitylene, leads to facile activation of an aromatic C−H bond of the solvent in all cases, and formation of aryl/hydride bridged magnesium(II) products, [{(Priso)Mg} 2 (μ‐H)(μ‐Ar)] 4 – 9 . In contrast to similar reactions reported for β‐diketiminate coordinated counterparts of 3 , these C−H activations proceed with little regioselectivity, though they are considerably faster. Reaction of 3 with an excess of the pyridine, p ‐NC 5 H 4 Bu t (py Bu t ), gave [(Priso)Mg(py Bu t H)(py Bu t ) 2 ] 10 , presumably via reduction of the pyridine to yield a radical intermediate, [(Priso)Mg(py Bu t ⋅)(py Bu t ) 2 ] 11 , which then abstracts a proton from the reaction solvent or a reactant. DFT calculations suggest two possible pathways to the observed arene C−H activations. One of these involves photochemical cleavage of the Mg−Mg bond of 3 , generating magnesium(I) doublet radicals, (Priso)Mg⋅. These then doubly reduce the arene substrate to give “Birch‐like” products, which subsequently rearrange via C−H activation of the arene. Circumstantial evidence for the photochemical generation of transient magnesium radical species includes the fact that irradiation of a cyclohexane solution of 3 leads to an intramolecular aliphatic C−H activation process and formation of an alkyl‐bridged magnesium(II) species, [{Mg(μ‐Priso −H )} 2 ] 12 . Furthermore, irradiation of a 1 : 1 mixture of 3 and the β‐diketiminato dimagnesium(I) compound, [{( Dip Nacnac)Mg} 2 ] ( Dip Nacnac=[HC(MeCNDip) 2 ] − ), effects a “scrambling” reaction, and the near quantitative formation of an unsymmetrical dimagnesium(I) compound, [(Priso)Mg−Mg( Dip Nacnac)] 13 . Finally, the EPR spectrum (77 K) of a glassed solution of UV irradiated 3 is dominated by a broad featureless signal, indicating the presence of a doublet radical species.

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

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗
DOE OSTI2020

An organogel library for solution NMR analysis of nanoparticle suspensions in non-aqueous samples

Surface contrast solution NMR methods (scNMR) are emerging as powerful tools to investigate the adsorption of small molecule ligands to the surface of nanoparticles (NP), returning fundamental insight into the kinetics and thermodynamics of sorption, as well as structural information on the adsorbed species. A prerequisite for the acquisition of high quality solution NMR data is the preparation of homogeneous and stable samples that return consistent NMR spectra and allow extensive signal averaging. Unfortunately, this condition does not apply to NMR samples containing NPs that often show a tendency to sediment and accumulate at the bottom of the NMR tube over the course of the experiment. We have recently shown that preparing NMR samples in an agarose gel matrix inhibits sedimentation and allows the characterization of small molecule-NP interactions by scNMR. Unfortunately, as the agarose gel only forms in aqueous solution, this sample preparation method cannot be used to stabilize NP suspensions in a non-aqueous environment. In this work, we introduce a library of 48 organogels, based on low molecular-mass organic gelators (LMOGs), to prepare NMR samples of small molecule/NP systems in a wide range of organic solvents. In addition, we present a simple method that takes advantage of 1H transverse relaxation (1H-R2) measurements to screen the library and identify the best gelator to characterize the small molecule-NP interaction of interest in the solvent of choice. We expect the results of this study will enable the preparation of homogeneous and stable samples of NPs in non-aqueous environments, therefore dramatically increasing the applicability of scNMR to the characterization of heterogeneous interactions and to the investigation of the role played by solvent molecules in regulating the kinetics and thermodynamics of sorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Fluorescence Studies of Lysozyme Nucleation

Fluorescence is one of the most powerful tools available for the study of macromolecules. For example, fluorescence can be used to study self association through methods such as anisotropy (the rotational rate of the molecule in solution), quenching (the accessibility of a bound probe to the bulk solution), and resonance energy transfer (measurement of the distance between two species). Fluorescence can also be used to study the local environment of the probe molecules, and the changes in that environment which accompany crystal nucleation and growth. However fluorescent techniques have been very much underutilized in macromolecular growth studies. One major advantage is that the fluorescent species generally must be at low concentration, typically ca 10-5 to 10-6 M. Thus one can study a very wide range of solution conditions, ranging from very high to very low protein concentration, he latter of which are not readily accessible to scattering techniques. We have prepared a number of fluorescent derivatives of chicken egg white lysozyme (CEWL). Fluorescent probes have been attached to two different sites, ASP 101 and the N-terrninal amine, with a sought for use in different lines of study. Preliminary resonance energy transfer studies have been -carried out using pyrene acetic acid (Ex 340 mn, Em 376 nm) lysozyme as a donor and cascade blue (Ex 377 run, Em 423 nm) labeled lysozyme as an acceptor. The emission of both the pyrene and cascade blue probes was followed as a function of the salt protein concentrations. The data show an increase in cascade blue and a concomitant decrease in the pyrene fluorescence as either the salt or protein concentrations are increased, suggesting that the two species are approaching each other close enough for resonance energy transfer to occur. This data can be analyzed to measure the distance between the probe molecules and, knowing their locations on the protein molecule their distances from and orientations with respect to each other. The results of these and other studies will be discussed.

Pusey, Marc L.↗
NASA NTRS1976

A hybrid computer program for rapidly solving flowing or static chemical kinetic problems involving many chemical species

A hybrid chemical kinetic computer program was assembled which provides a rapid solution to problems involving flowing or static, chemically reacting, gas mixtures. The computer program uses existing subroutines for problem setup, initialization, and preliminary calculations and incorporates a stiff ordinary differential equation solution technique. A number of check cases were recomputed with the hybrid program and the results were almost identical to those previously obtained. The computational time saving was demonstrated with a propane-oxygen-argon shock tube combustion problem involving 31 chemical species and 64 reactions. Information is presented to enable potential users to prepare an input data deck for the calculation of a problem.

Mclain, A. G.↗
NASA NTRS1980

A comparison of transport algorithms for premixed, laminar steady state flames

The effects of different methods of approximating multispecies transport phenomena in models of premixed, laminar, steady state flames were studied. Five approximation methods that span a wide range of computational complexity were developed. Identical data for individual species properties were used for each method. Each approximation method is employed in the numerical solution of a set of five H2-02-N2 flames. For each flame the computed species and temperature profiles, as well as the computed flame speeds, are found to be very nearly independent of the approximation method used. This does not indicate that transport phenomena are unimportant, but rather that the selection of the input values for the individual species transport properties is more important than the selection of the method used to approximate the multispecies transport. Based on these results, a sixth approximation method was developed that is computationally efficient and provides results extremely close to the most sophisticated and precise method used.

Coffee, T. P.↗
DOE OSTI2020

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

In this paper, we consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Overcoming genetic paucity of Camelina sativa : possibilities for interspecific hybridization conditioned by the genus evolution pathway

Camelina or false flax (Camelina sativa) is an emerging oilseed crop and a feedstock for biofuel production. This species is believed to originate from Western Asian and Eastern European regions, where the center of diversity of the Camelina genus is located. Cultivated Camelina species arose via a series of polyploidization events, serving as bottlenecks narrowing genetic diversity of the species. The genetic paucity of C. sativa is foreseen as the most crucial limitation for successful breeding and improvement of this crop. A potential solution to this challenge could be gene introgression from Camelina wild species or from resynthesized allohexaploid C. sativa. However, both approaches would require a complete comprehension of the evolutionary trajectories that led to the C. sativa origin. Although there are some studies discussing the origin and evolution of Camelina hexaploid species, final conclusions have not been made yet. Here, we propose the most complete integrated evolutionary model for the Camelina genus based on the most recently described findings, which enables efficient improvement of C. sativa via the interspecific hybridization with its wild relatives. We also discuss issues of interspecific and intergeneric hybridization, aimed on improving C. sativa and overcoming the genetic paucity of this crop. The proposed comprehensive evolutionary model of Camelina species indicates that a newly described species Camelina neglecta has a key role in origin of tetra- and hexaploids, all of which have two C. neglecta-based subgenomes. Understanding of species evolution within the Camelina genus provides insights into further research on C. sativa improvements via gene introgression from wild species, and a potential resynthesis of this emerging oilseed crop.

09 BIOMASS FUELS↗
NASA NTRS2003

An Alternative Hypothesis for How Microgravity Improves Macromolecular Crystal Quality

There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests that for many proteins crystal nucleation and growth is by addition of associated species that are preformed by reversible concentration-driven self association processes in the bulk solution. We have developed a self-association model for the crystal nucleation and growth of the protein chicken egg lysozyme. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small an octamer in the 43 helix configuration (the proposed average sized growth unit) would have a M.W. approx. 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm (a 24-mer) commonly attach to the crystal. AFM results from Weichmann et al. (Ultramicroscopy 86, 159-166, 2001) suggest that associated species of up to 40-mers in size add to the (101) faces. These measurements reflect the sizes of units that both added and desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. On Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species is more likely to be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc↗
DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗