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.

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 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.

The filter system for tertiary treatment of sewage effluent by land application - its performance in a subtropical environment

FILTER is an innovative, CSIRO developed system for treating effluent using high rate land application and subsequent effluent recapture via a closely spaced, subsurface drainage network. We report on the summer performance of a FILTER system established in a subtropical environment on a relatively impermeable swelling clay soil underlain by a deep regional water table. Using secondary treated sewage effluent, the FILTER system produced effluent of tertiary nutrient standards (less than or equal to5 mg/L TN; less than or equal to1 mg/L TP), with salinity levels suitable for subsequent irrigation reuse (EC less than or equal to2.5 dS/m). Removal of faecal coliforms was considerably less effective. The hydraulic loading rate achieved was about two and a half times larger than conventional irrigation demand, but this was associated with high deep percolation losses (e 3 mm/day). Comparisons are made with the original FILTER system developed and tested by Jayawardane et al. in temperate Australia. Suggestions are made for modifications to, and further testing of FILTER in a subtropical environment.

Investigation of the flow field of a highly heated jet of air

Measurements of mean and fluctuating velocity and temperature and their self- and cross-products to the third-order are presented for a heated axisymmetric air jet. Froude numbers in the range of 3500 13,190, Reynolds numbers in the range of 3470-8500 and non-dimensional streamwise distances. X*, from 0.27 to 1.98 are covered by the data. It was found that turbulence intensity decreases for the heated jet in the region between the inertia dominated and the buoyancy dominated regions which is contrary to findings with helium jets mixing with ambient air to produce density fluctuations. The effects of heating on the turbulent kinetic energy budget and the temperature variance budget show small differences for the inertia dominated region and the intermediate region which help to explain the transition process to the far field plume region. Constants are evaluated for the isotropic eddy diffusivity and generalised gradient hypothesis models as well as the scalar variance model. No significant effect of heating on the dissipation time-scale ratio was found. A novel wire array with an inclined cold wire was used. Measurements obtained with this probe are found to lead to asymmetries in some of the higher-order products. Further investigation suggested that the asymmetries are attributable to an as yet unreported interference effect produced by the leading prong of the inclined temperature wire, The effect may also have implications for inclined velocity wires which contain a temperature component when used in heated flows. (C) 2002 Elsevier Science Inc. All rights reserved.

Evaluation of various pre-treatments for the dehydration of banana and selection of suitable drying models

The thin-layer drying behaviour of bananas in a beat pump dehumidifier dryer was examined. Four pre-treatments (blanching, chilling, freezing and combined blanching and freezing) were applied to the bananas, which were dried at 50 degreesC with an air velocity of 3.1 m s(-1) and with the relative humidity of the inlet air of 10-35%. Three drying models, the simple model, the two-term exponential model and the Page model were examined. All models were evaluated using three statistical measures, correlation coefficient, root means square error, and mean absolute percent error. Moisture diffusivity was calculated based on the diffusion equation for an infinite cylindrical shape using the slope method. The rate of drying was higher for the pre-treatments involving freezing. The sample which was blanched only did not show any improvement in drying rate. In fact, a longer drying time resulted due to water absorption during blanching. There was no change in the rate for the chilled sample compared with the control. While all models closely fitted the drying data, the simple model showed greatest deviation from the experimental results. The two-term exponential model was found to be the best model for describing the drying curves of bananas because its parameters represent better the physical characteristics of the drying process. Moisture diffusivities of bananas were in the range 4.3-13.2 x 10(-10) m(2)s(-1). (C) 2002 Published by Elsevier Science Ltd.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

Surface diffusion of strong adsorbing vapours on porous carbon

Surface diffusion of strongly adsorbing hydrocarbon vapours on activated carbon was measured by using a constant molar flow method (D.D. Do, Dynamics of a semi-batch adsorber with constant molar supply rate: a method for studying adsorption rate of pure gas, Chem. Eng. Sci. 50 (1995) 549), where pure adsorbate is introduced into a semi-batch adsorber at a constant molar flow rate. The surface diffusivity was determined from the analysis of pressure response versus time, using a linear mathematical model developed earlier. To apply the linear theory over the non-linear range of the adsorption isotherm, we implement a differential increment method on the system which is initially equilibrated with some pre-determined loading. By conducting the experiments at different initial loadings, the surface diffusivity can be extracted as a function of loading. Propane, n-butane, n-hexane, benzene, and ethanol were used as diffusing adsorbate on a commercial activated carbon. It is found that the surface diffusivity of these strongly adsorbing vapours increases rapidly with loading, and the surface diffusion flux contributes significantly to the total flux and cannot be ignored. The surface diffusivity increases with temperature according to the Arrhenius law, and for the paraffins tested it decreases with the molecular weight of the adsorbate. (C) 2002 Elsevier Science Ltd. All rights reserved.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Genetic variation in foliar carbon isotope composition in relation to tree growth and foliar nitrogen concentration in clones of the F-1 hybrid between slash pine and Caribbean pine

The objectives of this study were: (1) to quantify the genetic variation in foliar carbon isotope composition (delta(13)C) of 122 clones of ca. 4-year-old F-1 hybrids between slash pine (Pinus elliottii Engelm var. elliottii) and Caribbean pine (Pinus caribaea var. hondurensis Barr.,et Golf.) grown at two field experimental sites with different water and nitrogen availability in southeast Queensland, Australia, in relation to tree growth and foliar nitrogen concentration (N-mass); and (2) to assess the potential of using delta(13)C measurements, in the foliage materials collected from the clone hedges at nursery and the 4-year-old tree canopies in the field, as an indirect index of tree water use efficiency for selecting elite F-1 hybrid pine clones with improved tree growth. There were significant differences in foliar delta(13)C between the nursery hedges and the 4-year-old tree canopies in the field, between the summer and winter seasons, between the two experimental sites, and between the upper outer and lower outer canopy positions sampled. This indicates that delta(13)C measurements in the foliage materials are significantly influenced by the sampling techniques and environmental conditions. Significant differences in foliar delta(13)C, at the upper outer canopy in both field experiments in summer and winter, were detected between the clones, and between the female parents of the clones. Clone means of tree height at age ca. 3 years were positively related to those of the upper outer canopy delta(13)C at both experimental sites in winter, but only for the wetter site in summer. There were positive, linear relationships between clone means of canopy delta(13)C and those of canopy N-mass, indicating that canopy photosynthetic capacity might be an important factor regulating the clonal variation in canopy delta(13)C. Significant correlations were found between clone means of canopy delta(13)C at both experimental sites in summer and winter, and between those at the upper outer and lower outer canopy positions. Mean clone delta(13)C for the nursery hedges was only positively related to mean clone stem diameter at 1.3 m height at age 3 years on the wetter site. The clone by site interaction for foliar delta(13)C at the upper outer canopy was significant only in summer. Overall, the relatively high genetic variance components for foliar delta(13)C and significant, positive correlations between clone means of foliar delta(13)C and tree growth have highlighted the potential of using foliar delta(13)C measurements for assisting in selection of the elite F-1 hybrid pine clones with improved tree growth. (C) 2002 Elsevier Science B.V. All rights reserved.

Sequencing batch reactor technology: the key to a BP refinery (Bulwer Island) upgraded environmental protection system - a low cost lagoon based retro-fit

BP Refinery (Bulwer Island) Ltd (BP) located on the eastern Australian coast is currently undergoing a major expansion as a part of the Queensland Clean Fuels Project. The associated wastewater treatment plant upgrade will provide a better quality of treated effluent than is currently possible with the existing infrastructure, and which will be of a sufficiently high standard to meet not only the requirements of imposed environmental legislation but also BP's environmental objectives. A number of challenges were faced when considering the upgrade, particularly; cost constraints and limited plot space, highly variable wastewater, toxicity issues, and limited available hydraulic head. Sequencing Batch Reactor (SBR) Technology was chosen for the lagoon upgrade based on the following; SBR technology allowed a retro-fit of the existing earthen lagoon without the need for any additional substantial concrete structures, a dual lagoon system allowed partial treatment of wastewaters during construction, SBRs give substantial process flexibility, SBRs have the ability to easily modify process parameters without any physical modifications, and significant cost benefits. This paper presents the background to this application, an outline of laboratory studies carried out on the wastewater and details the full scale design issues and methods for providing a cost effective, efficient treatment system using the existing lagoon system.

Exceptionally high-rate nitrification in sequencing batch reactors treating high ammonia landfill leachate

The nitrogen removal capacity of a suspended culture system treating mature landfill leachate was investigated. Leachate containing high ammonium levels of 300-900 mg N/L was nitrified in a bench scale sequencing batch reactor. Leachate from four different landfills was treated over a two year period for the removal of nitrogen. In this time, a highly specific nitrifying culture was attained that delivered exceptionally high rates of ammonia removal. No sludge was wasted from the system to increase the throughput and up to 13 g/L of MLSS was obtained. Settleability of the purely nitrifying biomass was excellent with SVI less than 40 mL/g, even at the high sludge concentrations. Nitrification rates up to 246 mg NI(L h) (5.91 g N/(L d)) and specific nitrification rates of 36 mg N/(gVSS h) (880 mg N/(gVSS d)) were obtained. The loading to the system at this time allowed complete nitrification of the leachate with a hydraulic retention time of only 5 hours. Following these successful treatability studies, a full-scale plant was designed and built at one of the landfills investigated.

Determination of the effect of lipophilicity on the in vitro permeability and tissue reservoir characteristics of topically applied solutes in human skin layers

In order to establish the relationship between solute lipophilicity and skin penetration (including flux and concentration behavior), we examined the in vitro penetration and membrane concentration of a series of homologous alcohols (C2-C10) applied topically in aqueous solutions to human epidermal, full-thickness, and dermal membranes. The partitioning/distribution of each alcohol between the donor solution, stratum corneum, viable epidermis, dermis, and receptor phase compartments was determined during the penetration process and separately to isolated samples of each tissue type. Maximum flux and permeability coefficients are compared for each membrane and estimates of alcohol diffusivity are made based on flux/concentration data and also the related tissue resistance (the reciprocal of permeability coefficient) for each membrane type. The permeability coefficient increased with increasing lipophilicity to alcohol C8 (octanol) with no further increase for C10 (decanol). Log vehicle:stratum corneum partition coefficients were related to logP , and the concentration of alcohols in each of the tissue layers appeared to increase with lipophilicity. No difference was measured in the diffusivity of smaller more polar alcohols in the three membranes; however, the larger more lipophilic solutes showed slower diffusivity values. The study showed that the dermis may be a much more lipophilic environment than originally believed and that distribution of smaller nonionized solutes into local tissues below a site of topical application may be estimated based on knowledge of their lipophilicity alone.

Effect of vehicle pretreatment on the flux, retention, and diffusion of topically applied penetrants in vitro

Purpose. The flux of a topically applied drug depends on the activity in the skin and the interaction between the vehicle and skin. Permeation of vehicle into the skin can alter the activity of drug and the properties of the skin barrier. The aim of this in vitro study was to separate and quantify these effects. Methods. The flux of four radiolabeled permeants (water, phenol, diflunisal, and diazepam) with log K-oct/water values from 1.4 to 4.3 was measured over 4 h through heat-separated human epidermis pretreated for 30 min with vehicles having Hildebrand solubility parameters from 7.9 to 23.4 (cal/cm(3))(1/2). Results. Enhancement was greatest after pretreatment with the more lipophilic vehicles. A synergistic enhancement was observed using binary mixtures. The flux of diazepam was not enhanced to the same extent as the other permeants, possibly because its partitioning into the epidermis is close to optimal (log K-oct 2.96). Conclusion. An analysis of the permeant remaining in the epidermis revealed that the enhancement can be the result of either increased partitioning of permeant into the epidermis or an increasing diffusivity of permeants through the epidermis.