Hostname: page-component-848d4c4894-mwx4w Total loading time: 0 Render date: 2024-06-15T16:33:32.990Z Has data issue: false hasContentIssue false

Effective reaction rate on a heterogeneous surface

Published online by Cambridge University Press:  29 September 2017

Ashok S. Sangani*
Affiliation:
Department of Biomedical and Chemical Engineering, Syracuse University, Syracuse, NY 13244, USA
*
Email address for correspondence: asangani@syr.edu

Abstract

We examine the problem of prescribing the macroscale boundary condition to the solute convective–diffusive mass transport equation at a heterogeneous surface consisting of reactive circular disks distributed uniformly on a non-reactive surface. The reaction rate at the disks is characterized by a first-order kinetics. This problem was examined by Shah & Shaqfeh (J. Fluid Mech., vol. 782, 2015, pp. 260–299) who obtained the boundary condition in terms of an effective first-order rate constant, which they determined as a function of the Péclet number $Pe=\dot{\unicode[STIX]{x1D6FE}}a^{2}/D$, the fraction $\unicode[STIX]{x1D719}$ of the surface area occupied by the reactive disks and the non-dimensional reaction rate constant $K=ka/D$. Here, $a$ is the radius of the disks, $D$ is the solute diffusivity, $\dot{\unicode[STIX]{x1D6FE}}$ is the wall shear rate and $k$ is the first-order surface-reaction rate constant. Their analysis assumed that $Pe$ and $K$ are $O(1)$ while the ratio of the microscale $a$ to the macroscale $H$ is small. The macroscale transport process is convection–diffusion dominated under these conditions. We examine here the case when the non-dimensional numbers based on the macroscale $H$ are $O(1)$. In this limit the microscale transport problem is reaction rate dominated. We find that the boundary condition can be expressed in terms of an effective rate constant only up to $O(\unicode[STIX]{x1D716})$, where $\unicode[STIX]{x1D716}=a/H$. Higher-order expressions for the mass flux involve both the macroscopic concentration and its surface gradient. The $O(\unicode[STIX]{x1D716})$ microscale problem is relatively easy to solve as the convective effects are unimportant and it is possible to obtain analytical expressions for the effective rate constant as a function of $\unicode[STIX]{x1D719}$ for both periodic and random arrangement of the disks without having to solve the boundary integral equation as was done by Shah and Shaqfeh. The results thus obtained are shown to be in good agreement with those obtained numerically by Shah and Shaqfeh for $Pe=0$. In a separate study, Shah et al. (J. Fluid Mech., vol. 811, 2017, pp. 372–399) examined the inverse-geometry problem in which the disks are inert and the rest of the surface surrounding them is reactive. We show that the two problems are related when $Pe=0$ and $kH/D=O(1)$. Finally, a related problem of determining the current density at a surface consisting of an array of microelectrodes is also examined and the analytical results obtained for the current density are found to agree well with the computed values obtained by solving the integral equation numerically by Lucas et al. (SIAM J. Appl. Maths, vol. 57(6), 1997, pp. 1615–1638) over a wide range of parameters characterizing this problem.

Type
Papers
Copyright
© 2017 Cambridge University Press 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Adda-Bedia, M., Katzav, E. & Vella, D. 2008 Solution of the percus-yevick equation for hard disks. J. Chem. Phys. 128 (18), 184508.Google Scholar
Bender, M. A. & Stone, H. A. 1993 An integral equation approach to the study of the steady state current at surface microelectrodes. J. Electroanalyt. Chem. 351 (1–2), 2955.Google Scholar
Chae, D. G., Ree, F. H. & Ree, T. 1969 Radial distribution functions and equation of state of the hard-disk fluid. J. Chem. Phys. 50 (4), 15811589.CrossRefGoogle Scholar
Dodd, T. L., Hammer, D. A., Sangani, A. S. & Koch, D. L. 1995 Numerical simulations of the effect of hydrodynamic interactions on diffusivities of integral membrane proteins. J. Fluid Mech. 293, 147180.Google Scholar
Juhasz, N. M. & Deen, W. M. 1991 Effect of local peclet number on mass transfer to a heterogeneous surface. Ind. Engng Chem. Res. 30 (3), 556562.Google Scholar
Ladd, A. J. C. 1990 Hydrodynamic transport coefficients of random dispersions of hard spheres. J. Chem. Phys. 93 (5), 34843494.CrossRefGoogle Scholar
Lucas, S. K., Sipcic, R. & Stone, H. A. 1997 An integral equation solution for the steady-state current at a periodic array of surface microelectrodes. SIAM J. Appl. Maths 57 (6), 16151638.CrossRefGoogle Scholar
Mo, G. & Sangani, A. S. 1994 A method for computing stokes flow interactions among spherical objects and its application to suspensions of drops and porous particles. Phys. Fluids 6 (5), 16371652.Google Scholar
Sangani, A. S. & Behl, S. 1989 The planar singular solutions of stokes and laplace equations and their application to transport processes near porous surfaces. Phys. Fluids A 1 (1), 2137.CrossRefGoogle Scholar
Shah, P. N., Lin, T. Y. & Shaqfeh, E. S. G. 2017 Heat/mass transport in shear flow over a reactive surface with inert defects. J. Fluid Mech. 811, 372399.Google Scholar
Shah, P. N. & Shaqfeh, E. S. G. 2015 Heat/mass transport in shear flow over a heterogeneous surface with first-order surface-reactive domains. J. Fluid Mech. 782, 260299.Google Scholar
Tio, K.-K. & Sadhal, S. S. 1991 Analysis of thermal constriction resistance with adiabatic circular gaps. J. Thermophys. Heat Transfer 5 (4), 550559.CrossRefGoogle Scholar