Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

Finite element analysis of mechanical properties in discontinuously reinforced metal matrix composites with ultrafine microstructure

This investigation focused on the finite element analyses of elastic and plastic properties of aluminium/alumina composite materials with ultrafine microstructure. The commonly used unit cell model was used to predict the elastic properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the plastic properties of the composites and the associated strengthening mechanism within the metal matrix material were investigated. The grain size of the matrix material was found to be an important factor influencing the mechanical properties of the composites studied. (C) 1997 Elsevier Science S.A.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

The secondary passive film for type 304 stainless steel in 0.5 M H2SO4

An examination has been carried out of the secondary passive film on Type 304 stainless steel in 0.5 M H2SO4. The characterization techniques used were electrochemical (potentiodynamic; potentiostatic, and film reduction experiments) and surface analytical. A bilayer model for the secondary passive film is proposed. It appears that next to the metal, there is a modified passive film which controls the electrochemical response; i.e., governs the current for any applied potential. On top of this modified passive film, the experimental data are consistent with a ''porous'' corrosion-product film which adds to the total film thickness but has little influence on the electrochemical response. The composition of the secondary passive film corresponds most probably to a mixed Fe/Cr oxide/hydroxide enriched in Cr3+, With a composition similar to a primary passive film.

Precursor scavenging of the resistive grain-boundary phase in 8 mol% yttria-stabilized zirconia: Effect of trace concentrations of SiO2

The influence that trace concentrations Of SiO2 have on improving grain-boundary conduction via precursor scavenging using additional heat treatment at 1200 degreesC for 40 h before sintering was investigated. At a SiO2-impurity level (SIL) less than or equal to 160 ppm by weight, the grain-boundary resistivity (p(gb)) decreased to 20% of its value, while no improvement in grain-boundary conduction was found at a SIL greater than or equal to 310 ppm. The correlation between the resistance per unit grain-boundary area, p(gb), and average grain size indicated that the inhomogeneous distribution of the siliceous phase in the sample with a SIL greater than or equal to 310 ppm. hampered the scavenging reaction.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

The use of liquid phase adsorption isotherms for characterization of activated carbons

The characterization of three commercial activated carbons was carried out using the adsorption of various compounds in the aqueous phase. For this purpose the generalized adsorption isotherm was employed, and a modification of the Dubinin-Radushkevich pore filling model, incorporating repulsive contributions to the pore potential as well as bulk liquid phase nonideality, was used as the local isotherm. Eight different flavor compounds were used as adsorbates, and the isotherms were jointly fitted to yield a common pore size distribution for each carbon. The bulk liquid phase nonideality was incorporated through the UNIFAC activity coefficient model, and the repulsive contribution to the pore potential was incorporated through the Steele 10-4-3 potential model. The mean micropore network coordination number for each carbon was also determined from the fitted saturation capacity based on percolation theory. Good agreement between the model and the experimental data was observed. In addition, excellent agreement between the bimodal gamma pore size distribution and density functional theory-cum-regularization-based pore size distribution obtained by argon adsorption was also observed, supporting the validity of the model. The results show that liquid phase adsorption, using adsorptive molecules of different sizes, can be an effective means of characterizing the pore size distribution as well as connectivity. Alternately, if the carbon pore size distribution is independently known, the method can be used to measure critical molecular sizes. (C) 2001 Elsevier Science.

A generalised dynamic model for char particle gasification with structure evolution and peripheral fragmentation

A generalised model for the prediction of single char particle gasification dynamics, accounting for multi-component mass transfer with chemical reaction, heat transfer, as well as structure evolution and peripheral fragmentation is developed in this paper. Maxwell-Stefan analysis is uniquely applied to both micro and macropores within the framework of the dusty-gas model to account for the bidisperse nature of the char, which differs significantly from the conventional models that are based on a single pore type. The peripheral fragmentation and random-pore correlation incorporated into the model enable prediction of structure/reactivity relationships. The occurrence of chemical reaction within the boundary layer reported by Biggs and Agarwal (Chem. Eng. Sci. 52 (1997) 941) has been confirmed through an analysis of CO/CO2 product ratio obtained from model simulations. However, it is also quantitatively observed that the significance of boundary layer reaction reduces notably with the reduction of oxygen concentration in the flue gas, operational pressure and film thickness. Computations have also shown that in the presence of diffusional gradients peripheral fragmentation occurs in the early stages on the surface, after which conversion quickens significantly due to small particle size. Results of the early commencement of peripheral fragmentation at relatively low overall conversion obtained from a large number of simulations agree well with experimental observations reported by Feng and Bhatia (Energy & Fuels 14 (2000) 297). Comprehensive analysis of simulation results is carried out based on well accepted physical principles to rationalise model prediction. (C) 2001 Elsevier Science Ltd. AH rights reserved.

Surface diffusion and adsorption of hydrocarbons in activated carbon

Adsorption and diffusion in a porous media were studied theoretically and experimentally with a differential transient permeation method. The porous medium is allowed to equilibrate at some specified loading, and then the time trajectory of the permeation process is followed after a small difference between the pressures at the end faces of the porous medium is introduced at time t = 0 +. Such a trajectory us. time would contain adsorption and diffusion characteristics of the system. By studying this for various surface loadings, pore and surface diffusions can be fully characterized. Mathematical modeling of transient permeation is detailed for pure gases or vapors diffusion and adsorption in porous media. Effects of nonlinearity of adsorption isotherm, pressure, temperature and heat effects were considered in the model. Experimental data of diffusion and adsorption of propane, n-butane and n-hexane in activated carbon at different temperatures and loadings show the potential of this method as a useful tool to study adsorption kinetics in porous media. Validity of the model is best tested against the transient data where the kinetics curves exhibit sigmoidal shape, which is a result of the diffusion and adsorption rate during the initial stage of permeation.

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.

Discrimination of adsorption kinetic models for the description of hydrocarbon adsorption in activated carbon

Five kinetic models for adsorption of hydrocarbons on activated carbon are compared and investigated in this study. These models assume different mass transfer mechanisms within the porous carbon particle. They are: (a) dual pore and surface diffusion (MSD), (b) macropore, surface, and micropore diffusion (MSMD), (c) macropore, surface and finite mass exchange (FK), (d) finite mass exchange (LK), and (e) macropore, micropore diffusion (BM) models. These models are discriminated using the single component kinetic data of ethane and propane as well as the multicomponent kinetics data of their binary mixtures measured on two commercial activated carbon samples (Ajax and Norit) under various conditions. The adsorption energetic heterogeneity is considered for all models to account for the system. It is found that, in general, the models assuming diffusion flux of adsorbed phase along the particle scale give better description of the kinetic data.

Characterization of micro-mesoporous carbon media

A simple method to characterize the micro and mesoporous carbon media is discussed. In this method, the overall adsorption quantity is the sum of capacities of all pores (slit shape is assumed), in each of which the process of adsorption occurs in two sequential steps: the multi-layering followed by pore filling steps. The critical factor in these two steps is the enhancement of the pressure of occluded 'free' molecules in the pore as well as the enhancement of the adsorption layer thickness. Both of these enhancements are due to the overlapping of the potential fields contributed by the two opposite walls. The classical BET and modified Kelvin equations are assumed to be applicable for the two steps mentioned above, with the allowance for the enhanced pore pressure, the enhanced adsorption energy and the enhanced BET constant,all of which vary with pore width. The method is then applied to data of many carbon samples of different sources to derive their respective pore size distributions, which are compared with those obtained from DFT analysis. Similar pore size distributions (PSDs) are observed although our method gives sharper distribution. Furthermore, we use our theory to analyze adsorption data of nitrogen at 77 K and that of benzene at 303 K (ambient temperature). The PSDs derived from these two different probe molecules are similar, with some small differences that could be attributed to the molecular properties, such as the collision diameter. Permeation characteristics of sub-critical fluids are also discussed in this paper. (C) 2001 Elsevier Science B.V. All rights reserved.

On the surface diffusion of hydrocarbons in microporous activated carbon

This paper addresses the current status of the various diffusion theories for surface diffusion in the literature. The inadequacy of these models to explain the surface diffusion of many hydrocarbons in microporous activated carbon is shown in this paper. They all can explain the increase of the surface diffusivity (D-mu) with loading, but cannot explain the increase of the surface permeability (D(mu)partial derivativeC(mu)/partial derivativeP) with loading as observed in our data of diffusion of hydrocarbons in activated carbon, even when the surface heterogeneity is accounted for in those models. The explanation for their failure was presented, and we have put forward a theory to explain the increase of surface diffusion permeability with loading. This new theory assumes the variation of the activation energy for surface diffusion with surface loading, and it is validated with diffusion data of propane, n-butane, n-hexane, benzene and ethanol in activated carbon. (C) 2001 Elsevier Science Ltd. All rights reserved.

Surface diffusion of hydrocarbons in activated carbon: Comparison between constant molar flow, differential permeation and differential adsorption bed methods

This paper presents the comparison of surface diffusivities of hydrocarbons in activated carbon. The surface diffusivities are obtained from the analysis of kinetic data collected using three different kinetics methods- the constant molar flow, the differential adsorption bed and the differential permeation methods. In general the values of surface diffusivity obtained by these methods agree with each other, and it is found that the surface diffusivity increases very fast with loading. Such a fast increase can not be accounted for by a thermodynamic Darken factor, and the surface heterogeneity only partially accounts for the fast rise of surface diffusivity versus loading. Surface diffusivities of methane, ethane, propane, n-butane, n-hexane, benzene and ethanol on activated carbon are reported in this paper.