Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

On the validity of the shrinking core model in noncatalytic gas solid reaction

By using a matched asymptotic expansion technique, the shrinking core model (SCM) used in non-catalytic gas solid reactions with general kinetic expression is rigorously justified in this paper as a special case of the homogeneous model when the reaction rate is much faster than that of diffusion. The time-pendent velocity of the moving reacted-unreacted interface is found to be proportional to the gas flux at that interface for all geometries of solid particles, and the thickness order of the reaction zone and also the degree of chemical reaction at the interface is discussed in this paper.

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.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

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.

Capillary network models for transport in packed beds: Considerations of pore aspect ratio

The conventional analysis for the estimation of the tortuosity factor for transport in porous media is modified here to account for the effect of pore aspect ratio. Structural models of the porous medium are also constructed for calculating the aspect ratio as a function of porosity. Comparison of the model predictions with the extensive data of Currie (1960) for the effective diffusivity of hydrogen in packed beds shows good agreement with a network model of randomly oriented intersecting pores for porosities upto about 50 percent, which is the region of practical interest. The predictions based on this network model are also found to be in better agreement with the data of Currie than earlier expressions developed for unconsolidated and grainy media.

Variation of the pore structure of coal chars during gasification

The variation of the pore structure of several coal chars during gasification in air and carbon dioxide was studied by argon adsorption at 87 K and CO2 adsorption at 273 K. It is found that the surface area and volume of the small pores (10 Å for air gasification is constant over a wide range of conversion (>20%), while for CO2 gasification similar results are obtained using the total surface area. However, in the early stages of gasification (

The kinetics of NO and N2O reduction over coal chars in fluidised-bed combustion

This paper presents a comprehensive and critical review of the mechanisms and kinetics of NO and N2O reduction reaction with coal chars under fluidised-bed combustion conditions (FBC). The heterogeneous reactions of NO and N2O with char/carbon surface have been well recognised as the most important processes in reducing both NOx and N2O in situ FBC. Compared to NO-carbon reactions in FBC, the reactions of N2O with chars have been relatively less understood and studied. Beginning with the overall reaction schemes for both NO and N2O reduction, the paper extensively discusses the reaction mechanisms including the effects of active surface sites. Generally, NO- and N2O-carbon reactions follow a series of step reactions. However, questions remain concerning the role of adsorbed phases of NO and N2O, and the behaviour of different surface sites. Important kinetics factors such as the rate expressions, kinetics parameters as well as the effects of surface area and pore structure are discussed in detail. The main factors influencing the reduction of NO and N2O in FBC conditions are the chemical and physical properties of chars, and the operating parameters of FBC such as temperature, presence of CO, O-2 and pressure. It is shown that under similar conditions, N2O is more readily reduced on the char surface than NO. Temperature was found to be a very important parameter in both NO and N2O reduction. It is generally agreed that both NO- and N2O-carbon reactions follow first-order reaction kinetics with respect to the NO and N2O concentrations. The kinetic parameters for NO and N2O reduction largely depend on the pore structure of chars. The correlation between the char surface area and the reactivities of NO/N2O-char reactions is considered to be of great importance to the determination of the reaction kinetics. The rate of NO reduction by chars is strongly enhanced by the presence of CO and O-2, but these species may not have significant effects on the rate of N2O reduction. However, the presence of these gases in FBC presents difficulties in the study of kinetics since CO cannot be easily eliminated from the carbon surface. In N2O reduction reactions, ash in chars is found to have significant catalytic effects, which must be accounted for in the kinetic models and data evaluation. (C) 1997 Elsevier Science Ltd.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

Percolative fragmentation of char particles during gasification

Percolative fragmentation was confirmed to occur during gasification of three microporous coal chars. Indirect evidence obtained by the variation of electrical resistivity (ER) with conversion was supported by direct observation of numerous fragments during gasification. The resistivity increases slowly at low conversions and then sharply after a certain conversion value, which is a typical percolation phenomenon suggesting the occurrence of internal fragmentation at high conversion. Two percolation models are applied to interpret the experimental data and determine the percolation threshold. A percolation threshold of 0.02-0.07 was found, corresponding to a critical conversion of 92-96% for fragmentation. The electrical resistivity variation at high conversions is found to be very sensitive to diffusional effects during gasification. Partially burnt samples with a narrow initial particle size range were also observed microscopically, and found to yield a large number of small fragments even when the particles showed no disintegration and chemical control prevailed. It is proposed that this is due to the separation of isolated clusters from the particle surface. The particle size distribution of the fragments was essentially independent of the reaction conditions and the char type, and supported the prediction by percolation theory that the number fraction distribution varies linearly with mass in a log-log plot. The results imply that perimeter fragmentation would occur in practical combustion systems in which the reactions are strongly diffusion affected.

A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

Multicomponent supercritical adsorption in microporous activated carbon materials

In this paper, a new technique for predicting multicomponent adsorption equilibria of supercritical fluids in microporous carbons is presented. In difference from adsorption on a surface, which is a function of the fluid-solid interaction only, adsorption in porous media is influenced by the proximity of the pore walls, resulting in the enhancement in adsorption affinity. The degree of this enhancement is different for different adsorbates, and it increases with a decrease in pore size. The theory is applied to a number of carbonaceous systems with good success.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Mixed gas equilibrium adsorption on zeolites and energetic heterogeneity of adsorption volume

A model for binary mixture adsorption accounting for energetic heterogeneity and intermolecular interactions is proposed in this paper. The model is based on statistical thermodynamics, and it is able to describe molecular rearrangement of a mixture in a nonuniform adsorption field inside a cavity. The Helmholtz free energy obtained in the framework of this approach has upper and lower limits, which define a permissible range in which all possible solutions will be found. One limit corresponds to a completely chaotic distribution of molecules within a cavity, while the other corresponds to a maximum ordered molecular structure. Comparison of the nearly ideal O-2-N-2-zeolite NaX system at ambient temperature with the system Of O-2-N-2-zeolite CaX at 144 K has shown that a decrease of temperature leads to a molecular rearrangement in the cavity volume, which results from the difference in the fluid-solid interactions. The model is able to describe this behavior and therefore allows predicting mixture adsorption more accurately compared to those assuming energetic uniformity of the adsorption volume. Another feature of the model is its ability to correctly describe the negative deviations from Raoult's law exhibited by the O-2-N-2-CaX system at 144 K. Analysis of the highly nonideal CO2-C2H6-zeolite NaX system has shown that the spatial molecular rearrangement in separate cavities is induced by not only the ion-quadrupole interaction of the CO2 molecule but also the significant difference in molecular size and the difference between the intermolecular interactions of molecules of the same species and those of molecules of different species. This leads to the highly ordered structure of this system.

A carbon activation model with application to longan seed char gasification

In this paper a new structural model is presented to describe the evolution of porosity of char during the gasification process. The model assumes the char structure to be composed of bundles of parallel graphite layers, and the reactivities of each layer with the gasification agent are assumed to be different to represent the different degree of heterogeneity of each layer (i.e. each layer will react with the gasification agent at a different rate). It is this difference in the reactivity that allows micropores to be created during the course of gasification. This simple structural model enables the evolution of pore volume, pore geometrical surface area and the pore size distribution to be described with respect to the extent of char burn-off. The model is tested against the experimental data of gasification of longan seed-derived char with carbon dioxide and it is found that the agreement between the model and the data is reasonably satisfactory, especially the evolution of surface area and pore volume with burn-off.