Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Virus-sized colloid transport in a single pore: Model development and sensitivity analysis

A mathematical model is developed to simulate the transport and deposition of virus-sized colloids in a cylindrical pore throat considering various processes such as advection, diffusion, colloid-collector surface interactions and hydrodynamic wall effects. The pore space is divided into three different regions, namely, bulk, diffusion and potential regions, based on the dominant processes acting in each of these regions. In the bulk region, colloid transport is governed by advection and diffusion whereas in the diffusion region, colloid mobility due to diffusion is retarded by hydrodynamic wall effects. Colloid-collector interaction forces dominate the transport in the potential region where colloid deposition occurs. The governing equations are non-dimensionalized and solved numerically. A sensitivity analysis indicates that the virus-sized colloid transport and deposition is significantly affected by various pore-scale parameters such as the surface potentials on colloid and collector, ionic strength of the solution, flow velocity, pore size and colloid size. The adsorbed concentration and hence, the favorability of the surface for adsorption increases with: (i) decreasing magnitude and ratio of surface potentials on colloid and collector, (ii) increasing ionic strength and (iii) increasing pore radius. The adsorbed concentration increases with increasing Pe, reaching a maximum value at Pe = 0.1 and then decreases thereafter. Also, the colloid size significantly affects particle deposition with the adsorbed concentration increasing with increasing particle radius, reaching a maximum value at a particle radius of 100 nm and then decreasing with increasing radius. System hydrodynamics is found to have a greater effect on larger particles than on smaller ones. The secondary minimum contribution to particle deposition has been found to increase as the favorability of the surface for adsorption decreases. The sensitivity of the model to a given parameter will be high if the conditions are favorable for adsorption. The results agree qualitatively with the column-scale experimental observations available in the literature. The current model forms the building block in upscaling colloid transport from pore scale to Darcy scale using Pore-Network Modeling. (C) 2014 Elsevier By. All rights reserved.

About the thermal stability and pore elimination in the ordered hexagonal mesoporous silica SBA-15

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Application of density functional theory to capillary phenomena in cylindrical mesopores with radial and longitudinal density distributions

In this paper, we applied a version of the nonlocal density functional theory (NLDFT) accounting radial and longitudinal density distributions to study the adsorption and desorption of argon in finite as well as infinite cylindrical nanopores at 87.3 K. Features that have not been observed before with one-dimensional NLDFT are observed in the analysis of an inhomogeneous fluid along the axis of a finite cylindrical pore using the two-dimensional version of the NLDFT. The phase transition in pore is not strictly vapor-liquid transition as assumed and observed in the conventional version, but rather it exhibits a much elaborated feature with phase transition being complicated by the formation of solid phase. Depending on the pore size, there are more than one phase transition in the adsorption-desorption isotherm. The solid formation in finite pore has been found to be initiated by the presence of the meniscus. Details of the analysis of the extended version of NLDFT will be discussed in the paper. (C) 2004 American Institute of Physics.

Structural characterization of folded pentapeptides containing centrally positioned beta(R)Val, gamma(R)Val and gamma(S)Val residues

A cylindrical pore of similar to 7.5 angstrom diameter containing a one-dimensional water wire, within the confines of a hydrophobic channel lined with the valine side chain, has been observed in crystals of the peptide Boc-D-Pro-Aib-Val-Aib-Val-OMe (1) (Raghavender et al., 2009, 2010). The synthesis and structural characterization in crystals of three backbone homologated analogues Boc-D-Pro-Aib-beta(3)(R) Val-Aib-Val-OMe (2), Boc-D-Pro-Aib-gamma(4)(R)Val-Aib-Val-OMe (3), Boc-D-Pro-Aib-gamma(4)(S)Val-Aib-Val-OMe (4) are described. Crystal structures of peptides 2, 3 and 4 reveal close-packed arrangements in which no pore was formed. In peptides 2 and 3 the N-terminus D-Pro-Aib segment adopted conformations closely related to Type II' beta-turns, while residues 2-4 form one turn of an alpha beta right-handed C-11 helix in 2 and an alpha gamma C-12 helix in 3. In peptide 4, a continuous left-handed helical structure was observed with the D-Pro-Aib segment forming a Type III' beta-turn, followed by one turn of a left-handed alpha gamma C-12 helix. (C) 2012 Elsevier Ltd. All rights reserved.

Diffusion and chemical reaction in the porous structures of solid oxide fuel cells

The Rolls-Royce Integrated-Planar Solid Oxide Fuel Cell (IP-SOFC) consists of ceramic modules which have electrochemical cells printed on the outer surfaces. The cathodes are the outermost layer of each cell and are supplied with oxygen from air flowing over the outside of the module. The anodes are in direct contact with the ceramic structure and are supplied with fuel from internal gas channels. Natural gas is reformed into hydrogen for use by the fuel cells in a separate reformer module of similar design except that the fuel cells are replaced by a reforming catalyst layer. The performance of the modules is intrinsically linked to the behaviour of the gas flows within their porous structures. Because the porous layers are very thin, a one-dimensional flow model provides a good representation of the flow property variations between fuel channel and fuel cell or reforming catalyst. The multi-component convective-diffusive flows are simulated using a new theory of flow in porous material, the Cylindrical Pore Interpolation Model. The effects of the catalysed methane reforming and water-gas shift chemical reactions are also considered using appropriate kinetic models. It is found that the shift reaction, which is catalysed by the anode material, has certain beneficial effects on the fuel cell module performance. In the reformer module it was found that the flow resistance of the porous support structure makes it difficult to sustain a high methane conversion rate. Although the analysis is based on IP-SOFC geometry, the modelling approach and general conclusions are applicable to other types of SOFC.

Calcinação e estabilidade térmica de Sílicas mesoporoas preparadas com Pluronic P123

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Inverted critical adsorption of polyelectrolytes in confinement

What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer-surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density sigma(c)-defining the adsorption-desorption transition-are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of sigma(c) for the concave interfaces versus the Debye screening length 1/kappa and the extent of confinement a for these three interfaces for small kappa a values. For large kappa a the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/kappa. We also rationalize how sigma(c)(kappa) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids.

Computer simulations of colloidal fluids in confinement

Computer-Simulationen von Kolloidalen Fluiden in Beschränkten Geometrien Kolloidale Suspensionen, die einen Phasenübergang aufweisen, zeigen eine Vielfalt an interessanten Effekten, sobald sie auf eine bestimmte Geometrie beschränkt werden, wie zum Beispiel auf zylindrische Poren, sphärische Hohlräume oder auf einen Spalt mit ebenen Wänden. Der Einfluss dieser verschiedenen Geometrietypen sowohl auf das Phasenverhalten als auch auf die Dynamik von Kolloid-Polymer-Mischungen wird mit Hilfe von Computer-Simulationen unter Verwendung des Asakura-Oosawa- Modells, für welches auf Grund der “Depletion”-Kräfte ein Phasenübergang existiert, untersucht. Im Fall von zylindrischen Poren sieht man ein interessantes Phasenverhalten, welches vom eindimensionalen Charakter des Systems hervorgerufen wird. In einer kurzen Pore findet man im Bereich des Phasendiagramms, in dem das System typischerweise entmischt, entweder eine polymerreiche oder eine kolloidreiche Phase vor. Sobald aber die Länge der zylindrischen Pore die typische Korrelationslänge entlang der Zylinderachse überschreitet, bilden sich mehrere quasi-eindimensionale Bereiche der polymerreichen und der kolloidreichen Phase, welche von nun an koexistieren. Diese Untersuchungen helfen das Verhalten von Adsorptionshysteresekurven in entsprechenden Experimenten zu erklären. Wenn das Kolloid-Polymer-Modellsystem auf einen sphärischen Hohlraum eingeschränkt wird, verschiebt sich der Punkt des Phasenübergangs von der polymerreichen zur kolloidreichen Phase. Es wird gezeigt, dass diese Verschiebung direkt von den Benetzungseigenschaften des Systems abhängt, was die Beobachtung von zwei verschiedenen Morphologien bei Phasenkoexistenz ermöglicht – Schalenstrukturen und Strukturen des Janustyps. Im Rahmen der Untersuchung von heterogener Keimbildung von Kristallen innerhalb einer Flüssigkeit wird eine neue Simulationsmethode zur Berechnung von Freien Energien der Grenzfläche zwischen Kristall- bzw. Flüssigkeitsphase undWand präsentiert. Die Resultate für ein System von harten Kugeln und ein System einer Kolloid- Polymer-Mischung werden anschließend zur Bestimmung von Kontaktwinkeln von Kristallkeimen an Wänden verwendet. Die Dynamik der Phasenseparation eines quasi-zweidimensionalen Systems, welche sich nach einem Quench des Systems aus dem homogenen Zustand in den entmischten Zustand ausbildet, wird mit Hilfe von einer mesoskaligen Simulationsmethode (“Multi Particle Collision Dynamics”) untersucht, die sich für eine detaillierte Untersuchung des Einflusses der hydrodynamischen Wechselwirkung eignet. Die Exponenten universeller Potenzgesetze, die das Wachstum der mittleren Domänengröße beschreiben, welche für rein zwei- bzw. dreidimensionale Systeme bekannt sind, können für bestimmte Parameterbereiche nachgewiesen werden. Die unterschiedliche Dynamik senkrecht bzw. parallel zu den Wänden sowie der Einfluss der Randbedingungen für das Lösungsmittel werden untersucht. Es wird gezeigt, dass die daraus resultierende Abschirmung der hydrodynamischen Wechselwirkungsreichweite starke Auswirkungen auf das Wachstum der mittleren Domänengröße hat.

Capillary phenomena in the framework of the two-dimensional density functional theory

We present results of application of the density functional theory (DFT) to adsorption and desorption in finite and infinite cylindrical pores accounting for the density distribution in radial and axial directions. Capillary condensation via formation of bridges is considered using canonical and grand canonical versions of the 2D DFT. The potential barrier of nucleation is determined as a function of the bulk pressure and the pore diameter. In the framework of the conventional assumptions on intermolecular interactions both 1D and 2D DFT versions lead to the same results and confirm the classical scenario of condensation and evaporation: the condensation occurs at the vapor-like spinodal point, and the evaporation corresponds to the equilibrium transition pressure. The analysis of experimental data on argon and nitrogen adsorption on MCM-41 samples seems to not completely corroborate this scenario, with adsorption branch being better described by the equilibrium pressure - diameter dependence. This points to the necessity of the further development of basic representations on the hysteresis phenomena.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.

Orientation of Syndiotactic Polystyrene Crystallized in Cylindrical Nanopores

Syndiotactic polystyrene (sPS) nanorods with different diameters have been prepared by using anodic aluminum oxide templates, and the orientation of the sPS crystals in the nanorods has been investigated by FT-IR spectroscopy. It is found that the c axis of the beta' crystals preferentially oriented perpendicular to the axis of the nanorod, and the degree of orientation is lower as the diameter of the nanorod decreases. This unexpected result is attributed to nuclei formed at the surface of the nanopores and their subsequent growth, in addition to the preferential growth compatible with the pore direction by the nuclei formed in the bulk film and in the nanorods.

Large-pore mesoporous SBA-15 silica particles with submicrometer size as stationary phases for high-speed CEC separation

A new mesoporous sphere-like SBA-15 silica was synthesized and evaluated in terms of its suitability as stationary phases for CEC. The unique and attractive properties of the silica particle are its submicrometer particle size of 400 nm and highly ordered cylindrical mesopores with uniform pore size of 12 nm running along the same direction. The bare silica particles with submicrometer size have been successfully employed for the normal-phase electrochromatographic separation of polar compounds with high efficiency (e.g., 210 000 for thiourea), which is matched well with its submicrometer particle size. The Van Deemeter plot showed the hindrance to mass transfer because of the existence of pore structure. The lowest plate height of 2.0 mu m was obtained at the linear velocity of 1.1 mm/s. On the other hand, because of the relatively high linear velocity (e.g., 4.0 mm/s) can be generated, high-speed separation of neutral compounds, anilines, and basic pharmaceuticals in CEC with C-18-modified SBA-15 silica as stationary phases was achieved within 36, 60, and 34 s, respectively.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory

Adsorption of argon and nitrogen at their respective boiling points in cylindrical pores of MCM-41 type silica-like adsorbents is studied by means of a non-local density functional theory (NLDFT), which is modified to deal with amorphous solids. By matching the theoretical results of the pore filling pressure versus pore diameter against the experimental data, we arrive at a conclusion that the adsorption branch (rather than desorption) corresponds to the true thermodynamic equilibrium. If this is accepted, we derive the optimal values for the solid–fluid molecular parameters for the system amorphous silica–Ar and amorphous silica–N2, and at the same time we could derive reliably the specific surface area of non-porous and mesoporous silica-like adsorbents, without a recourse to the BET method. This method is then logically extended to describe the local adsorption isotherms of argon and nitrogen in silica-like pores, which are then used as the bases (kernel) to determine the pore size distribution. We test this with a number of adsorption isotherms on the MCM-41 samples, and the results are quite realistic and in excellent agreement with the XRD results, justifying the approach adopted in this paper.