A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.

A modified discrete random pore model allowing for different initial surface reactivity

We investigate here a modification of the discrete random pore model [Bhatia SK, Vartak BJ, Carbon 1996;34:1383], by including an additional rate constant which takes into account the different reactivity of the initial pore surface having attached functional groups and hydrogens, relative to the subsequently exposed surface. It is observed that the relative initial reactivity has a significant effect on the conversion and structural evolution, underscoring the importance of initial surface chemistry. The model is tested against experimental data on chemically controlled char oxidation and steam gasification at various temperatures. It is seen that the variations of the reaction rate and surface area with conversion are better represented by the present approach than earlier random pore models. The results clearly indicate the improvement of model predictions in the low conversion region, where the effect of the initially attached functional groups and hydrogens is more significant, particularly for char oxidation. It is also seen that, for the data examined, the initial surface chemistry is less important for steam gasification as compared to the oxidation reaction. Further development of the approach must also incorporate the dynamics of surface complexation, which is not considered here.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

The importance of wet-powder dynamic mechanical properties in understanding granulation

Granule impact deformation has long been recognised as important in determining whether or not two colliding granules will coalesce. Work in the last 10 years has highlighted the fact that viscous effects are significant in granulation. The relative strengths of different formulations can vary with strain rate. Therefore, traditional strength measurements made at pseudo-static conditions give no indication, even qualitatively, of how materials will behave at high strain rates, and hence are actually misleading when used to model granule coalescence. This means that new standard methods need to be developed for determining the strain rates encountered by granules inside industrial equipment and also for measuring the mechanical properties of granules at these strain rates. The constitutive equations used in theoretical models of granule coalescence also need to be extended to include strain-rate dependent components.

Dissolution of a solid sphere in an unbounded, stagnant liquid

An exact analytical solution is obtained for the transient dissolution of solid spheres in a diffusion-controlled environment. This result provides a useful reference point for drug testing in humans. The dimensionless solution is expressed in terms of a single parameter, which accounts for solubility, bulk flow, and stagnant fluid composition. A simple, explicit and exact expression was found to predict time-to-complete dissolution (TCD). An approximate solution was also found which tracks the exact case for low solubility conditions.

Molecular recognition of macrocyclic peptidomimetic inhibitors by HIV-1 protease

High-resolution crystal structures are described for seven macrocycles complexed with HIV-1 protease (HIVPR). The macrocycles possess two amides and an aromatic group within 15-17 membered rings designed to replace N- or C-terminal tripeptides from peptidic inhibitors of HIVPR. Appended to each macrocycle is a transition state isostere and either an acyclic peptide, nonpeptide, or another macrocycle. These cyclic analogues are potent inhibitors of HIVPR, and the crystal structures show them to be structural mimics of acyclic peptides, binding in the active site of HIVPR via the same interactions. Each macrocycle is restrained to adopt a P-strand conformation which is preorganized for protease binding. An unusual feature of the binding of C-terminal macrocyclic inhibitors is the interaction between a positively charged secondary amine and a catalytic aspartate of HIVPR. A bicyclic inhibitor binds similarly through its secondary amine that lies between its component N-terminal and C-terminal macrocycles. In contrast, the corresponding tertiary amine of the N-terminal macrocycles does not interact with the catalytic aspartates. The amine-aspartate interaction induces a 1.5 Angstrom N-terminal translation of the inhibitors in the active site and is accompanied by weakened interactions with a water molecule that bridges the ligand to the enzyme, as well as static disorder in enzyme flap residues. This flexibility may facilitate peptide cleavage and product dissociation during catalysis. Proteases [Aba(67,95)]HIVPR and [Lys(7),Ile(33),Aba(67,95)]- HIVPR used in this work were shown to have very similar crystal structures.