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.

A novel nucleation apparatus for regime separated granulation

A novel nucleation apparatus is presented for the production of narrow sized nuclei from various powder and binder liquid combinations. Mono-sized binder liquid droplets are produced by a specially designed mono-disperse droplet generator. The droplet generator is positioned above a conveyor belt, transporting a powder bed through the spray zone of the droplet generator. By nucleating powder on a conveyer belt, the nucleation mechanism is completely separated from all other granulation mechanisms due to the lack of relative motion between primary particles and/or formed nuclei. Nucleation tests were performed using chalcopyrite and limestone powders with water as the binder liquid. At all operating conditions, the formed nuclei were found to originate from multiplicities of drops that merged on the powder bed surface. Investigation of the dynamics of nuclei formation showed that powder-binder liquid combinations with fast penetration dynamics result in less variation in the number of droplets from which nuclei originate. Smaller and more narrowly distributed nuclei were also achieved by increasing powder speed through the spray zone.

Dimensionless spray flux in wet granulation: Monte-Carlo simulations and experimental validation

Dimensionless spray flux Ψa is a dimensionless group that characterises the three most important variables in liquid dispersion: flowrate, drop size and powder flux through the spray zone. In this paper, the Poisson distribution was used to generate analytical solutions for the proportion of nuclei formed from single drops (fsingle) and the fraction of the powder surface covered by drops (fcovered) as a function of Ψa. Monte-Carlo simulations were performed to simulate the spray zone and investigate how Ψa, fsingle and fcovered are related. The Monte-Carlo data was an excellent match with analytical solutions of fcovered and fsingle as a function of Ψa. At low Ψa, the proportion of the surface covered by drops (fcovered) was equal to Ψa. As Ψa increases, drop overlap becomes more dominant and the powder surface coverage levels off. The proportion of nuclei formed from single drops (fsingle) falls exponentially with increasing Ψa. In the ranges covered, these results were independent of drop size, number of drops, drop size distribution (mono-sized, bimodal and trimodal distributions), and the uniformity of the spray. Experimental data of nuclei size distributions as a function of spray flux were fitted to the analytical solution for fsingle by defining a cutsize for single drop nuclei. The fitted cutsizes followed the spray drop sizes suggesting that the method is robust and that the cutsize does indicate the transition size between single drop and agglomerate nuclei. This demonstrates that the nuclei distribution is determined by the dimensionless spray flux and the fraction of drop controlled nuclei can be calculated analytically in advance.

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.

Liquid distribution as a means to describing the granule growth mechanism

Copper concentrate (chalcopyrite) was granulated in a rotating drum with a diameter of 0.3 m and a length of 0.2 m. Water was used as the binder and it was sprayed onto the powder bed with a nozzle. This material exhibited induction type behaviour, which was defined by Iveson and Litster [AIChE J. 44 (1998) 1510]. Induction type behaviour is characterized by the occurrence of an induction stage, during which the granules are gradually being compacted and little or no growth occurs. At the end of this induction stage, binder liquid is squeezed from the interior of the granules onto the granule surface and the granules are then surface-wet. This results in a rapid growth rate of the granules. Different types of experiments were conducted. The influence of the nozzle pressure and the distance from the nozzle to the powder bed on the growth behaviour of the granules as well as on the binder distribution was examined. The results of these experiments led to the postulation of a modified mechanism for induction type behaviour: it was found that after the binder was delivered, there were large granules containing a high amount of binder and small granules containing less binder. During the induction stage, the granules are compacted and binder liquid continuously appears at the surface of the large granules. These wet spots that are continuously being formed pick up the dry and small granules. When all the small granules have been picked up, further expulsion of binder liquid onto the granules' surface results in granules that remain surface-wet. This phenomenon marks the end of the induction stage and it coincides with the disappearance of the small granules. The hypothesis was tested by selectively removing the smaller granules during an experiment. As expected, this resulted in a shorter induction time.

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.

Undergraduate student experiences and perceptions of group projects

This paper reports the experiences and perceptions of student group projects of three cohorts of undergraduates: a 3rd year capstone course in professional communication (N = 54), a 3rd year elective in chemical engineering (N = 29), and a core 2nd-year course in chemical engineering (N = 74).

Useful representations of series solutions for transport problems

Simple techniques are presented for rearrangement of an infinite series in a systematic way such that the convergence of the resulting expression is accelerated. These procedures also allow calculation of required boundary derivatives. Several examples of conduction and diffusion-reaction problems illustrate the methods.

Enzyme deactivation in fixed bed reactors with michaelis-menten kinetics

Rate expression for enzyme poisoning which are consistent with a Michaelis-Menten main reaction are used to analyze the performance of a fixed bed reactor containing immobilized enzyme. When enzyme deactivation results from the irreversible bonding of a product molecule to an existing substrate-enzyme complex, it is shown that minimum enzyme activity can occur in the interior of the bed, well away from the ends. This suggests that bed sectioning techniques may enable direct evaluation of fundamental poisoning mechanisms.

Fixed bed reactors with catalyst poisoning - first order kinetics.

The long performance of an isothermal fixed bed reactor undergoing catalyst poisoning is theoretically analyzed using the dispersion model. First order reaction with dth order deactivation is assumed and the model equations are solved by matched asymptotic expansions for large Peclet number. Simple closed-form solutions, uniformly valid in time, are obtained.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR

Analytical expressions are developed for the time-dependent reactant concentration and catalyst activity in an isothermal CSTR with Langmuir-Hinshelwood kinetics of deactivation and reaction. Several parallel and series posioning mechanisms are considered for a reactor which, without poisoning, would operate at a unique steady state. The use of matched asymptotic expansions and abandonment of the usual initial-steady-state assumption give results, valid from startup to final loss of activity, whose accuracy can be improved systematically.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR. II. Multiple pseudosteady states and reactor failure

Analytical expressions are derived for the time and magnitude of failure of an isothermal CSTR with substrate-inhibited kinetics, caused by slow catalyst deactivation under three types of parallel and series mechanisms. Reactors operating at high space velocity are found to be most susceptible to early failure and poisoning by product is more dangerous than by reactant. The magnitude of the jump across steady states depends solely on the Langmuir-Hinshelwood kinetic parameters and a detailed analysis of reactor behavior during the jump itself is given.

Enzyme immobilization in porous solid supports - penetration of immobilized enzyme

The process of enzyme immobilization under the diffusion-controlled regime (i.e., fast attachment of enzyme compared to its diffusion) is modeled and theoretically solved in this article. Simple and compact solutions for the penetration depth of immobilized enzyme and the bulk enzyme concentration versus time are presented. Furthermore, the conditions for the validity of our solutions are also given in this article so that researchers can discover when the theoretical solutions can be applied to their systems.

Analysis of multiple gas solid reactions

Multiple gas solid reactions involving one solid and N gaseous reactants are investigated in this study by using a matched asymptotic expansion technique. Two cases are particularly studied. In the first case all N chemical reaction rates are faster than the diffusion rate. While in the second case only M (M < N) chemical reaction rates are faster than the diffusion rate and the rates of the remaining (N-M) chemical reactions are comparable to that of diffusion. For these two cases the solid concentration profile behaves like a travelling wave. In the first case the wave front velocity is contributed linearly by all gaseous reactants (additive law) while in the second case this law does not hold.