Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion processes”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Incoherent Dynamics of Intense Proton Beams Under Electron Cooling

The minimum emittance of ion beams achieved using electron cooling is limited by the heating processes of Intra Beam Scattering and diffusion driven by resonance crossing of particles due to space-charge. We describe a new experiment to explore the intense space-charge regime with a transverse tune shift approaching -0.5 using 2.5 MeV protons at the Integrable Optics Test Accelerator (IOTA) at Fermilab. We also report on the results from PyORBIT simulations incorporating transverse space-charge and electron cooling with emphasis on the incoherent dynamics of the particles.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗
DOE OSTI2017

Three-Dimensional Tracking of Interfacial Hopping Diffusion

Theoretical predictions have suggested that molecular motion at interfaces—which influences processes including heterogeneous catalysis, (bio)chemical sensing, lubrication and adhesion, and nanomaterial self-assembly—may be dominated by hypothetical “hops” through the adjacent liquid phase, where a diffusing molecule readsorbs after a given hop according to a probabilistic “sticking coefficient.” Here, we use three-dimensional (3D) single-molecule tracking to explicitly visualize this process for human serum albumin at solid-liquid interfaces that exert varying electrostatic interactions on the biomacromolecule. Following desorption from the interface, a molecule experiences multiple unproductive surface encounters before readsorption. An average of approximately seven surface collisions is required for the repulsive surfaces, decreasing to approximately two and a half for surfaces that are more attractive. The hops themselves are also influenced by long-range interactions, with increased electrostatic repulsion causing hops of longer duration and distance. These findings explicitly demonstrate that interfacial diffusion is dominated by biased 3D Brownian motion involving bulk-surface coupling and that it can be controlled by influencing short- and long-range adsorbate-surface interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Automated Analysis, Reporting, and Archiving for Robotic Nondestructive Assay of Holdup Deposits - 19508

To decommission deactivated gaseous diffusion enrichment facilities, miles of contaminated pipe must be measured. The current method requires thousands of manual measurements, repeated manual data transcription, and months of manual analysis. The Pipe Crawling Activity Measurement System (PCAMS), developed by Carnegie Mellon University and in commissioning for use at the DOE Portsmouth Gaseous Diffusion Enrichment Facility, uses a robot to measure Uranium-235 from inside pipes and automatically log the data. Radiation measurements, as well as imagery, geometric modeling, and precise measurement positioning data are digitally transferred to the PCAMS server. On the server, data can be automatically processed in minutes and summarized for analyst review. Measurement reports are auto-generated with the push of a button. A database specially-configured to hold heterogeneous data such as spectra, images, and robot trajectories serves as archive. This paper outlines the features and design of the PCAMS Post-Processing Software, currently in commissioning for use at the Portsmouth Gaseous Diffusion Enrichment Facility. The analysis process, the analyst interface to the system, and the content of auto-generated reports are each described. Example pipe-interior geometric surface models, illustration of how key report features apply in operational runs, and user feedback are discussed. (authors)

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

Swelling and Diffusion during Methanol Sorption into Hydrated Nafion

Diffusion within polymer electrolyte membranes is often coincident with time-dependent processes such as swelling and polymer relaxation, which are factors that limit their ability to block molecular crossover during use. The solution-diffusion model of membrane permeation, which is the accepted theory for dense polymers, applies only to steady-state processes and does not address dynamic internal structural changes that can accompany permeation. To begin discovery of how such changes can be coupled to the permeation process, we have constructed here a stochastic multiscale reaction-diffusion model that examines time-dependent methanol uptake into and swelling of hydrated Nafion. Several potential mechanisms of diffusion and polymer response are tested. The simulation predictions are compared to real-time Fourier transform infrared attenuated total reflectance spectroscopy (FTIR-ATR) absorbance reported in the literature [ Hallinan , D. T. , Jr. ; Elabd , Y. A. J. Phys. Chem. B 2007 , 111 , 13221 - 13230 ]. Of the proposed polymer response mechanisms, only one, a reaction-limited, local response to increasing methanol concentration that takes the entire experimental time frame of 600 s, produces simulated FTIR-ATR data consistent with experiment. The simulations show that water diffusion out of the membrane is minimal during methanol sorption and that changes in the measured infrared absorbances are due primarily to the increase in methanol concentration accompanied by dilution of water during swelling. Swelling involves densification of the polymer structure even as there is an overall volume expansion of the film. Potential connections between the polymer densification and molecular-level structural changes of Nafion in methanol are discussed. These results indicate that the interaction between methanol and Nafion serves to increase Nafion's capacity to accommodate large volumes of methanol-water solutions, facilitating increased permeation across the membrane relative to pure water.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tunable Atomic Layer Deposition into Ultra-High-Aspect-Ratio (>60000:1) Aerogel Monoliths Enabled by Transport Modeling

Atomic layer deposition (ALD) modification of ultra-high-aspect-ratio structures (>10,000:1) is a powerful platform with applications in catalysis, filtration, and energy conversion. However, the deposition of conformal and tunable ALD coatings at these aspect ratios remains challenging, resulting in empirical tradeoffs between precursor utilization and reaction time. Here, we demonstrate tunable control of the ALD infiltration depth into an aerogel monolith (AM) and develop a reaction-diffusion model to accurately describe the coating process. Specifically, we investigate the ALD exposure time and precursor dose needed to conformally coat a silica AM with pore sizes of ~20 nm, a monolith thickness of ~2.5 mm, and aspect ratios exceeding 60,000:1. We demonstrate complete infiltration into the AM, which is quantified by elemental mapping. A reaction-diffusion model is developed, which accounts for multiple doses and precursor depletion in the ALD chamber during an exposure step. The experimentally validated model enables the prediction and tuning of infiltration depth into a tortuous, high-aspect-ratio structure such as an AM, allowing for the synthesis of rationally designed material architectures. Additionally, the model allows for co-optimization of the total deposition time and percentage of unreacted precursor, which are important for the manufacturability and sustainability of ALD processing. Lastly, we demonstrate that ultra-thin ALD Al 2 O 3 coatings can be used to stabilize silica AMs against structural degradation under high-temperature annealing conditions (700-800 °C), by limiting changes in surface area and monolith volume. As a result, this improved high-temperature stability has implications for numerous aerogel applications, including catalysis and thermal insulation.

14 SOLAR ENERGY↗
DOE OSTI2017

Aminopolymer Mobility and Support Interactions in Silica-PEI Composites for CO 2 Capture Applications: A Quasielastic Neutron Scattering Study

Composite gas sorbents, formed from an active polymer phase and a porous support, are promising materials for the separation of acid gases from a variety of gas streams. Significant changes in sorption performance (capacity, rate, stability etc.) can be achieved by tuning the properties of the polymer and the nature of interactions between polymer and support. We utilize quasielastic neutron scattering (QENS) and coarse-grained molecular dynamics (MD) simulations to characterize the dynamic behavior of the most commonly reported polymer in such materials, poly(ethylenimine) (PEI), both in bulk form and when supported in a mesoporous silica framework. The polymer chain dynamics (rotational and translational diffusion) are characterized using two neutron backscattering spectrometers that have overlapping time scales, ranging from picoseconds to a few nanoseconds. Two modes of motion are detected for the PEI molecule in QENS. At low energy transfers, a “slow process” on the time scale of ~200 ps is found and attributed to jump-mediated, center-of-mass diffusion. Second, a “fast process” at ~20 ps scale is also found and is attributed to a locally confined, jump-diffusion. Characteristic data (time scale and spectral weight) of these processes are compared to those characterized by MD, and reasonable agreement is found. For the nanopore-confined PEI, we observe a significant reduction in the time scale of polymer motion as compared to the bulk. The impacts of silica surface functionalization and of polymer fill fraction in the silica pores (controlling the portion of polymer molecules in contact with the pore walls), are both studied in detail. Hydrophobic functionalization of the silica leads to an increase of the PEI mobility above that in native silanol-terminated silica, but the dynamics are still slower than those in bulk PEI. Sorbents with faster PEI dynamics are also found to be more efficient for CO 2 capture, possibly because sorption sites are more accessible than those in systems with slower PEI dynamics. Therefore, this work supports the existence of a link between the affinity of the support for PEI and the accessibility of active sorbent functional groups.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Reversible Tuning of the Surface Plasmon Resonance of Indium Tin Oxide Nanocrystals by Gas-Phase Oxidation and Reduction

Heavily doped oxide nanocrystals exhibit a tunable localized surface plasmon resonance (LSPR) in the infrared, a property that is promising for applications in photonics, spectroscopy, and photochemistry. Nanocrystal carrier density and, thus, spectral response are adjustable via chemical reaction; yet, the fundamental processes that govern this behavior are poorly understood. In this work, we study the time dependence of the LSPR supported by indium tin oxide (ITO) nanocrystals during O 2 and N 2 annealing with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). We show that the LSPR red-shifts upon oxidation in O 2 and blue-shifts to its original position upon reduction in N 2 . A reaction–diffusion model allows us to rationalize the underlying physicochemical processes and quantitatively connect nanocrystal redox chemistry, solid-state diffusion, carrier density, and the LSPR.

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

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗
DOE OSTI2019

Synthesis, Characterization and Experimental Studies of Nano Zn–Al–Fe{sub 3}O{sub 4} Blended Alginate/Ca Beads for the Adsorption of Rhodamin B

In the present study, adsorptive removal of rhodamin-B (RMN-B) from aqueous solution was investigated by using nano Zn–Al–Fe{sub 3}O{sub 4} incorporated sodium alginate (SA) beads. SA is a natural polysaccharide, non toxic and biodegradable material. The incorporation of Zn–Al–Fe{sub 3}O{sub 4} into SA was done by ionic polymerization method using CaCl{sub 2}. A batch mode study was used for the adsorption of RMN-B onto Zn–Al–Fe{sub 3}O{sub 4}/ alginate beads. The adsorption efficiency of RMN-B was investigated by the following modeling equations like Freundlich, Langmuir, Pseudo-first order, Pseudo-second order, Intraparticle diffusion, Gibbs free energy, Enthalpy and Entropy. The predicted experimental results of various modeling equations were implied that the process follow Langmuir and pseudo-second order equations. The intraparticle diffusion parameters indicate that the diffusion between RMN-B and nano beads were predominant in throughout of process. The surface morphology of nano beads of before and after adsorption were show that RMN-B is strongly adhered on whole surface of the nano material. The FTIR spectrum indicates that the RMN-B strongly adsorbed on the surface and TGA study reveals that the alginate/Zn–Al–Fe{sub 3}O{sub 4}/Ca material has an excellent temperature withstand capacity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Microstructural Modeling of the α + β Phase in Ti-6Al-4V: A Diffusion-Based Approach

Complex heat treatment operations and advanced manufacturing processes such as laser or electron-beam welding will see the metallic workpiece experience a considerable range of temperatures and heating/cooling rates. These intrinsic conditions will have a significant bearing upon the microstructure of the material, and in turn upon the thermo-mechanical properties. In this work, a diffusion-based approach to model the growth and shrinkage of precipitates in the alpha + beta field of the common titanium alloy Ti-6Al-4V is established. Further, the numerical model is extended using a JMA-type approach to explore the dependency of the beta-transus temperature on extremely high heating rates, whereby dissolution alone is too slow to accurately describe the alpha to beta-phase transformation. Experimental heat treatments at rates of 5, 50, and 500 °C/s were performed, and metallographic analysis of the samples was used to validate the numerical modeling framework predictions for lamellar shrinkage, while data from the literature has been used to evaluate the numerical modeling framework for the growth of equiaxed microstructures. The agreement between measurements and numerical predictions was found to be good.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Skyrmion relaxation dynamics in the presence of quenched disorder

Using Langevin molecular dynamics simulations we study relaxation processes of interacting skyrmion systems with and without quenched disorder. Here, using the typical diffusion length as the time-dependent length characterizing the relaxation process, we find that clean systems always display dynamical scaling, and this even in cases where the typical length is not a simple power law of time. In the presence of the Magnus force, two different regimes are identified as a function of the noise strength. The Magnus force has also a major impact when attractive pinning sites are present, as this velocity-dependent force helps skyrmions to bend around defects and avoid caging effects. With the exception of the limit of large noise, for which dynamical scaling persists even in the presence of quenched disorder, attractive pinning sites capture a substantial fraction of skyrmions which results in a complex behavior of the two-time autocorrelation function that is not reproduced by a simple aging scaling ansatz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

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

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Enabling hydrogen chemisorption on charged graphene

Two-dimensional (2D) materials, including graphitic derivatives, have long been of interest for hydrogen storage applications, due to their high theoretical storage capacity, low weight, and other useful properties. However, poor kinetics for hydrogen adsorption and surface diffusion as part of the proposed spillover process for hydrogenation have limited their technological potential. Here, we use first-principles calculations to study electronic doping as a means to improve hydrogen chemisorption on graphene, which we use here as a proxy for graphitic derivatives more broadly. We find that positively charged graphene sheets have vastly improved kinetics for hydrogen diffusion and adsorption, while they limit unwanted hydrogen desorption. This combination of effects should favor hydrogen chemisorption via spillover. We connect these trends to the C–H bond, which introduces states near the Fermi level. These states are depopulated as electrons are removed, thereby lowering the bond energy and permitting more facile movement of hydrogen. Our results suggest that spillover mechanisms for hydrogen chemisorption should be revisited if strategies to apply a large charge to graphitic systems can be realized. Moreover, switchable application of the charge may lead to the reversible chemisorption of hydrogen. While the large magnitude of charging required suggests that graphene itself may not be suitable for reversible hydrogen chemisorption, the factors we identify and discuss could significantly boost the prospects of graphitic derivatives and other 2D or layered materials for hydrogen storage applications.

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

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

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

Persistent compositions of non-stoichiometric compounds with low bulk diffusivity: A theory and application to Nb3Sn superconductors

Non-stoichiometric compounds may develop a composition gradient when they are formed by reactive diffusional processes. This paper reports an interesting phenomenon that in compounds with low bulk diffusivities, which rely mainly on grain boundary diffusion for their growth, the final bulk compositions may be far from equilibrium, with a very low bulk diffusivity leading to fixed (persistent) compositions – one such example is Nb 3 Sn, a superconductor. We investigated the microchemistry at the reactive interface using atom probe tomography to clarify the diffusion reaction mechanism for this low-bulk-diffusivity case and thus propose a theory for what determines the compound composition, using Nb 3 Sn as an example for concreteness. Using certain approximations, we derive an explicit analytical equation that illustrates what factors determine its composition profile. We compare our model with the known facts of Nb 3 Sn and see good agreement. In particular, this model predicts that internal oxidation may lead to higher Sn contents than conventional, non-oxidized Nb 3 Sn. Our measurements show that this is indeed true, and that the internally-oxidized Nb 3 Sn also has higher upper critical fields, achieving up to 28.2 T at 4.2 K. We discuss the general applicability of this model to non-stoichiometric compounds with low bulk diffusivity, and propose it as a tool to help in the design and processing of such materials for compositional control.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Efficient reactive Brownian dynamics

We develop a Split Reactive Brownian Dynamics (SRBD) algorithm for particle simulations of reaction-diffusion systems based on the Doi or volume reactivity model, in which pairs of particles react with a specified Poisson rate if they are closer than a chosen reactive distance. In our Doi model, we ensure that the microscopic reaction rules for various association and dissociation reactions are consistent with detailed balance (time reversibility) at thermodynamic equilibrium. The SRBD algorithm uses Strang splitting in time to separate reaction and diffusion and solves both the diffusion-only and reaction-only subproblems exactly, even at high packing densities. To efficiently process reactions without uncontrolled approximations, SRBD employs an event-driven algorithm that processes reactions in a time-ordered sequence over the duration of the time step. A grid of cells with size larger than all of the reactive distances is used to schedule and process the reactions, but unlike traditional grid-based methods such as reaction-diffusion master equation algorithms, the results of SRBD are statistically independent of the size of the grid used to accelerate the processing of reactions. We use the SRBD algorithm to compute the effective macroscopic reaction rate for both reaction-limited and diffusion-limited irreversible association in three dimensions and compare to existing theoretical predictions at low and moderate densities. We also study long-time tails in the time correlation functions for reversible association at thermodynamic equilibrium and compare to recent theoretical predictions. Finally, we compare different particle and continuum methods on a model exhibiting a Turing-like instability and pattern formation. Our studies reinforce the common finding that microscopic mechanisms and correlations matter for diffusion-limited systems, making continuum and even mesoscopic modeling of such systems difficult or impossible. We also find that for models in which particles diffuse off lattice, such as the Doi model, reactions lead to a spurious enhancement of the effective diffusion coefficients.

97 MATHEMATICS AND COMPUTING↗