Diffusion and surface reaction in random fibrous media

Abstract

Alternative solution for diffusion to two spheres with first-order surface reactionIn the classic ''trapping'' problem, the reactant is uniformly generated at a constant rate across a matrix phase, and diffuses to the interface of a dispersed cylinder phase, where it instantly reacts. Prager's upper bound on the effective reaction rate constant k ͑a ''best'' bound based on void-point nearest-neighbor-surface statistics͒ is derived for randomly placed, freely overlapping, infinitely long cylinders of radius a, where each cylinder has an arbitrary orientation with respect to the others. To compare with simulation data for an isotropic bed of overlapping spheres of radius a, the variational upper bound is considered for the case of the three-dimensional, isotropic overlapping cylinder bed with totally random mutual cylinder orientations. Once a correction factor of 3/2 for the surface area ratio of the overlapping sphere to overlapping cylinder beds is applied, the analytical variational cylinder bed bound is nearly coincident with the sphere dispersion simulation curve. The small differences observed at lower dispersion densities are consistent with the expected diffusion flux differences in the sphere and cylinder geometry. By comparing the variational upper bound with other simulation results for randomly overlapping, aligned, elongated, prolate spheroids, a maximum range for the k change due to mutual rotation between neighboring cylindrical reaction sites ͑29% and 32%, respectively, at solid volume fractions of 0.30 and 0.50͒ is obtained. The necessary distribution moments for the corresponding ''relaxation time'' lower bound are also given. The Doi lower bound on the effective reaction rate ͑a ''best'' bound based on two-point void-void F vv , void-surface F vs , and surface-surface F ss correlations͒ is addressed. The twopoint correlations for an isotropic bed of overlapping cylinders with random mutual orientation in three dimensions are presented. Their behavior is discussed, and a logarithmic singularity in F ss at the two-point distance of 2a is pointed out, that renders the Doi bound indeterminate.

Diffusion and reaction for a spherical source and sink

“…Various sink shapes and distribution geometries have been considered. A periodic lattice of identical spheres; 16 an ensemble distribution of hard spheres; [17][18][19]23 uncorrelated overlapping spheres; 4,20,21 as well as long fibers with various mutual orientations of the cylindrical axes 22 have been studied. To balance the steady state, a zeroth-order, uniform volumetric generation source of reactant is maintained outside of the traps.…”