Stochastic modeling of non-linear adsorption with Gaussian kernel density estimators

Adsorption is a relevant process in many fields, such as product manufacturing or pollution remediation in porous materials. Adsorption takes place at the molecular scale, amenable to be modeled by Lagrangian numerical methods. We have proposed a chemical diffusion-reaction model for the simulation of adsorption, based on the combination of a random walk particle tracking method involving the use of Gaussian Kernel Density Estimators. The main feature of the proposed model is that it can effectively reproduce the nonlinear behavior characteristic of the Langmuir and Freundlich isotherms. In the former, it is enough to add a finite number of sorption sites of homogeneous sorption properties, and to set the process as the combination of the forward and the backward reactions, each one of them with a prespecified reaction rate. To model the Freundlich isotherm instead, typical of low to intermediate range of solute concentrations, there is a need to assign a different equilibrium constant to each specific sorption site, provided they are all drawn from a truncated power-law distribution. Both nonlinear models can be combined in a single framework to obtain a typical observed behavior for a wide range of concentration values.

[1]  O. Redlich,et al.  A USEFUL ADSORPTION ISOTHERM , 1959 .

[2]  D. Gillespie The chemical Langevin equation , 2000 .

[3]  David A. Benson,et al.  Connecting the dots: Semi-analytical and random walk numerical solutions of the diffusion-reaction equation with stochastic initial conditions , 2014, J. Comput. Phys..

[4]  Michael Dumbser,et al.  Numerical simulations of solute transport in highly heterogeneous formations: A comparison of alternative numerical schemes , 2013 .

[5]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[6]  Peter Salamon,et al.  A review and numerical assessment of the random walk particle tracking method. , 2006, Journal of contaminant hydrology.

[7]  E. Teller,et al.  ADSORPTION OF GASES IN MULTIMOLECULAR LAYERS , 1938 .

[8]  Andreas Englert,et al.  Mixing, spreading and reaction in heterogeneous media: a brief review. , 2011, Journal of contaminant hydrology.

[9]  Rocha,et al.  Freundlich's Isotherm Extended by Statistical Mechanics , 1997, Journal of colloid and interface science.

[10]  Alberto Guadagnini,et al.  Interpretation of column experiments of transport of solutes undergoing an irreversible bimolecular reaction using a continuum approximation , 2010 .

[11]  C. Harvey,et al.  Reactive transport in porous media: a comparison of model prediction with laboratory visualization. , 2002, Environmental science & technology.

[12]  Kohei Urano,et al.  Equilibria for adsorption of organic compounds on activated carbons in aqueous solutions I. Modified Freundlich isotherm equation and adsorption potentials of organic compounds , 1981 .

[13]  R. Sips,et al.  On the Structure of a Catalyst Surface , 1948 .

[14]  Daniel Fernàndez-Garcia,et al.  An automatic locally-adaptive method to estimate heavily-tailed breakthrough curves from particle distributions , 2013 .

[15]  D. Gillespie A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[16]  W. J. Thomas Chemisorption of ethane at iron and nickel oxide , 1957 .

[17]  S. P. Garabedian,et al.  Large-scale dispersive transport in aquifers : field experiments and reactive transport theory , 1987 .

[18]  Syu Ono,et al.  Statistical mechanics of adsorption , 1950 .

[19]  James Stephen Marron,et al.  Comparison of data-driven bandwith selectors , 1988 .

[20]  David A. Benson,et al.  Simulation of chemical reaction via particle tracking: Diffusion‐limited versus thermodynamic rate‐limited regimes , 2008 .

[21]  I. Langmuir THE ADSORPTION OF GASES ON PLANE SURFACES OF GLASS, MICA AND PLATINUM. , 1918 .

[22]  Joachim Engel,et al.  An iterative bandwidth selector for kernel estimation of densities and their derivatives , 1994 .

[23]  Jesús Carrera,et al.  Coupling of mass transfer and reactive transport for nonlinear reactions in heterogeneous media , 2010 .

[24]  Redner,et al.  Fluctuation-dominated kinetics in diffusion-controlled reactions. , 1985, Physical review. A, General physics.

[25]  Harihar Rajaram,et al.  Conservative and sorptive forced‐gradient and uniform flow tracer tests in a three‐dimensional laboratory test aquifer , 2004 .

[26]  Georges Guiochon,et al.  Derivation and Application of a Jovanovic–Freundlich Isotherm Model for Single-Component Adsorption on Heterogeneous Surfaces , 1996 .

[27]  Vivek Kapoor,et al.  Experimental study of bimolecular reaction kinetics in porous media. , 2000 .

[28]  Raymond T. Davis,et al.  The Activated Adsorption of Nitrogen on a Finely Divided Tungsten Powder , 1946 .

[29]  Peter Salamon,et al.  Modeling mass transfer processes using random walk particle tracking , 2006 .

[30]  A. A. Ovchinnikov,et al.  Role of density fluctuations in bimolecular reaction kinetics , 1978 .

[31]  Ch Sheindorf,et al.  A Freundlich-type multicomponent isotherm , 1981 .

[32]  Monica Moroni,et al.  Analysis of dispersion in porous media via matched-index particle tracking velocimetry experiments , 2007 .

[33]  David A. Benson,et al.  A simple and efficient random walk solution of multi-rate mobile/immobile mass transport equations , 2009 .

[34]  Alexandre M. Tartakovsky,et al.  Incomplete mixing and reactions with fractional dispersion , 2012 .

[35]  P. J. Green,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[36]  X. Sanchez‐Vila,et al.  Optimal reconstruction of concentrations, gradients and reaction rates from particle distributions. , 2011, Journal of contaminant hydrology.

[37]  M. C. Jones,et al.  A Brief Survey of Bandwidth Selection for Density Estimation , 1996 .

[38]  Bernard W. Silverman,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[39]  Daniel M. Tartakovsky,et al.  A Bayesian approach to integrate temporal data into probabilistic risk analysis of monitored NAPL remediation , 2012 .

[40]  D. W. Pollock Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[41]  D. Costa,et al.  Adsorption and self-assembly of bio-organic molecules at model surfaces: A route towards increased complexity , 2015 .

[42]  T. L. Hill,et al.  Statistical Mechanics of Adsorption. VI. Localized Unimolecular Adsorption on a Heterogeneous Surface , 1949 .

[43]  P. Meakin,et al.  A smoothed particle hydrodynamics model for miscible flow in three-dimensional fractures and the two-dimensional Rayleigh-Taylor instability , 2005 .

[44]  Xavier Sanchez-Vila,et al.  Do we really need a large number of particles to simulate bimolecular reactive transport with random walk methods? A kernel density estimation approach , 2015, J. Comput. Phys..

[45]  Daniel Fernàndez-Garcia,et al.  Toward efficiency in heterogeneous multispecies reactive transport modeling: A particle‐tracking solution for first‐order network reactions , 2014 .

[46]  Andrew F. B. Tompson,et al.  Numerical simulation of chemical migration in physically and chemically heterogeneous porous media , 1993 .

[47]  J. K. Roberts,et al.  The adsorption of hydrogen on tungsten , 1934, Proceedings of the Royal Society of London. Series A - Mathematical and Physical Sciences.

[48]  D. E. Dougherty,et al.  Particle-grid methods for reacting flows in porous media with application to Fisher's equation , 1992 .