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.

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.

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.

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.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

A twin study of the etiology of prolonged fatigue and immune activation

Risk factors to prolonged fatigue syndromes (PFS) are controversial. Pre-morbid and/or current psychiatric disturbance, and/or disturbed cell-mediated immunity (CMI), have been proposed as etiologic factors. Self-report measures of fatigue and psychologic distress and three in vitro measures of CMI were collected from 124 twin pairs. Crosstwincrosstrait correlations were estimated for the complete monozygotic (MZ; 79 pairs) and dizygotic (DZ; 45 pairs) twin groups. Multivariate genetic and environmental models were fitted to explore the patterns of covariation between etiologic factors. For fatigue, the MZ correlation was more than double the DZ correlation (0.49 versus 0.16) indicating strong genetic control of familial aggregation. By contrast, for in vitro immune activation measures MZ and DZ correlations were similar (0.49–0.69 versus 0.42–0.53) indicating the etiologic role of shared environments. As small univariate associations were noted between prolonged fatigue and the in vitro immune measures (r = −0.07 to −0.12), multivariate models were fitted. Relevant etiologic factors included: a common genetic factor accounting for 48% of the variance in fatigue which also accounted for 4%, 6% and 8% reductions in immune activation; specific genetic factors for each of the in vitro immune measures; a shared environment factor influencing the three immune activation measures; and, most interestingly, unique environmental influences which increased fatigue but also increased markers of immune activation. PFS that are associated with in vitro measures of immune activation are most likely to be the consequence of current environmental rather than genetic factors. Such environmental factors could include physical agents such as infection and/or psychologic stress.

Articulatory dysfunction in multiple sclerosis: An electropalatographic study

Dysfunction of the articulatory subsystem (i.c.. the lips, tongue, and jaw) has bccn identified as a major contributor to the reduction in speech intelligibility experienced by a high proportion of people with multiple sclerosis (MS). In particular. consonant imprecision has been reported to be the articulatory deficit that contributes most to variations in overall intelligibility of MS speakers. Electropalatography(EPG) IS an instrurncntal technique that visually documents the location and timing of tongue-topalatc contacts during speech. Although such a technique would be valuablc in objectively assessing the articulatory disturbances exhibited by individuals with dysarthria ia motor speech disorder) associated with MS, to-date no such study ha< been reported. The aim of the present study was to use EPG to assess tongue-to-palate contact patterns and articulatory timing in patients with dysarthria associated with MS. A dysarthric participant with a diagnosis of definite MS was fitted with an acrylic EPG palate (Reading EPG.?) and asked to read aloud a list of single syllable words which contained lingual consonants in the word-initial position and in consonant clusters. Each mord was repeated five times. The EPG palate was specifically moulded to tit the participant's hard palate and contained 62 electrodes that detected the tongue contacts. A non-neurologically impaired participant matched for age and sex servcd as a control. The results of the study revealed that the tongue-to-palate contacts produced by the participant with MS varied from those produced by the control in a number of ways in regard to spatial configurations and timing characteristics exhibited. The rcsults arc discussed in relation to the neuropathophysiological effects of MS on speech production. The potcntial use of EPG in programs for treating speech disorders associated with MS will be highlightcd.

Exploring through cover - the integrated interpretation of high resolution aeromagnetic, airborne electromagnetic and ground gravity data from the Grant's Patch area, Eastern Goldfields Province, Archaean Yigarn Craton Part B; Gravity inversion as a

The cost and risk associated with mineral exploration in Australia increases significantly as companies move into deeper regolith-covered terrain. The ability to map the bedrock and the depth of weathering within an area has the potential to decrease this risk and increase the effectiveness of exploration programs. This paper is the second in a trilogy concerning the Grant's Patch area of the Eastern Goldfields. The recent development of the VPmg potential field inversion program in conjunction with the acquisition of high-resolution gravity data over an area with extensive drilling provided an opportunity to evaluate three-dimensional gravity inversion as a bedrock and regolith mapping tool. An apparent density model of the study area was constructed, with the ground represented as adjoining 200 m by 200 m vertical rectangular prisms. During inversion VPmg incrementally adjusted the density of each prism until the free-air gravity response of the model replicated the observed data. For the Grant's Patch study area, this image of the apparent density values proved easier to interpret than the Bouguer gravity image. A regolith layer was introduced into the model and realistic fresh-rock densities assigned to each basement prism according to its interpreted lithology. With the basement and regolith densities fixed, the VPmg inversion algorithm adjusted the depth to fresh basement until the misfit between the calculated and observed gravity response was minimised. The resulting geometry of the bedrock/regolith contact largely replicated the base of weathering indicated by drilling with predicted depth of weathering values from gravity inversion typically within 15% of those logged during RAB and RC drilling.

Quantitative analysis of the speed and accuracy of tongue movements in dysarthria subsequent to traumatic brain injury using electromagnetic articulography

It has been recognised that in order to study the displacement, timing and co-ordination of articulatory components (i.e., tongue. lips, jaw) in speech production it is desirable to obtain high-resolution movement data on multiple structures inside and outside the vocal tract. Until recently, with the exception of X-ray techniques such as cineradiography, the study 0. speech movements has been hindered by the inaccessibility of the oral cavity during speech. X-ray techniques are generally not used because of unacceptable radiation exposure. The aim of the present study was to demonstrate the use of a new physiological device, the electromagnetic articulograph, for assessing articulatory dysfunction subsequent to traumatic brain injury. The components of the device together with the measuring principle are described and data collected from a single case presented. A 19 year-old male who exhibited dysarthria subsequent to a traumatic brain injury was fitted wit 2 the electromagnetic articulograph (Carstens AG-100) and a kinematic analysis of his tongue movements during production of the lingual consonants it, s, k/ within single syllable words was performed. Examination of kinematic parameters including movemmt trajectories, velocity, and acceleration revealed differences in the speed and accuracy of his tongue movements compared to those produced by a non-neurologically impaired adult male. It was concluded that the articulograph is a useful device for diagnosing speed and accuracy disorders in tongue movements during speech and that the device has potential for incorporation into physiologically based rehabilitation programs as a real-time biofeedback instrument.

Marker concentration patterns of labelled leaf and stem particles in the rumen of cattle grazing bermuda grass (Cynodon dactylon) analysed by reference to a raft model

Large (>1600 mum), ingestively masticated particles of bermuda grass (Cynodon dactylon L. Pers.) leaf and stem labelled with Yb-169 and Ce-144 respectively were inserted into the rumen digesta raft of heifers grazing bermuda grass. The concentration of markers in digesta sampled from the raft and ventral rumen were monitored at regular intervals over approximately 144 h. The data from the two sampling sites were simultaneously fitted to two pool (raft and ventral rumen-reticulum) models with either reversible or sequential flow between the two pools. The sequential flow model fitted the data equally as well as the reversible flow model but the reversible flow model was used because of its greater application. The reversible flow model, hereafter called the raft model, had the following features: a relatively slow age-dependent transfer rate from the raft (means for a gamma 2 distributed rate parameter for leaf 0.0740 v. stem 0.0478 h(-1)), a very slow first order reversible flow from the ventral rumen to the raft (mean for leaf and stem 0.010 h(-1)) and a very rapid first order exit from the ventral rumen (mean of leaf and stem 0.44 h(-1)). The raft was calculated to occupy approximately 0.82 total rumen DM of the raft and ventral rumen pools. Fitting a sequential two pool model or a single exponential model individually to values from each of the two sampling sites yielded similar parameter values for both sites and faster rate parameters for leaf as compared with stem, in agreement with the raft model. These results were interpreted as indicating that the raft forms a large relatively inert pool within the rumen. Particles generated within the raft have difficulty escaping but once into the ventral rumen pool they escape quickly with a low probability of return to the raft. It was concluded that the raft model gave a good interpretation of the data and emphasized escape from and movement within the raft as important components of the residence time of leaf and stem particles within the rumen digesta of cattle.

Compact clinical MRI magnet design using a multilayer current density approach

In this work, a new method of optimization is successfully applied to the theoretical design of compact, actively shielded, clinical MRI magnets. The problem is formulated as a two-step process in which the desired current densities on multiple, cc-axial surface layers are first calculated by solving Fredholm equations of the first kind. Non-linear optimization methods with inequality constraints are then invoked to fit practical magnet coils to the desired current densities. The current density approach allows rapid prototyping of unusual magnet designs. The emphasis of this work is on the optimal design of short, actively-shielded MRI magnets for whole-body imaging. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric, and asymmetric MRI magnets. Magnet designs are presented for actively-shielded, symmetric magnets of coil length 1.0 m, which is considerably shorter than currently available designs of comparable dsv size. Novel, actively-shielded, asymmetric magnet designs are also presented in which the beginning of a 50-cm dsv is positioned just 11 cm from the end of the coil structure, allowing much improved access to the patient and reduced patient claustrophobia. Magn Reson Med 45:331540, 2001. (C) 2001 Wiley-Liss, Inc.

Lot-fed goats - the advantages of using an enriched environment

By reducing the stress associated with high stocking density in a feedlot it is likely that goats will utilise their feed more efficiently and suffer fewer health problems. One possible method of reducing stress is to enrich the feedlot environment. In a feedlot of 160 castrated goats (average weight 24.6 +/- 2.6 kg), 2 types of feedlot structures were compared; a typical feedlot and an 'environmentally enriched' feedlot, and goats stocked at 2 densities; a high stocking density (1667 goats/ha) and a relatively low stocking density (897 goats/ha). Over a 42-day period, environmental enrichment (old car/truck tyres and wooden railway sleepers to climb on and PVC piping to mouth and butt) increased weight gain by 83% (P = 0.04) and reduced the occurrence of inanition by 36%. Aggressive behaviour at the feed trough was reduced by 30% (P = 0.03) in pens of lower density when compared with pens of higher density. We conclude that lot feeding goats in environmentally enriched surroundings will increase feed conversion and reduce the number of non-eaters when compared with typically bare pen structures. Doubling in liveweight gain of goats within enriched feedlot surroundings when compared with typical structures shows promise as a cost effective, animal-welfare-orientated practice.

Developmental genetics of red cell indices during puberty: A longitudinal twin study

Red cell number and size increase during puberty, particularly in males. The aim of the present study was to determine whether expression of genes affecting red cell indices varied with age and sex. Haemoglobin, red cell count, and mean cellular volume were measured longitudinally on 578 pairs of twins at twelve, fourteen and sixteen years of age. Data were analysed using a structural equation modeling approach, in which a variety of univariate and longitudinal simplex models were fitted to the data. Significant heritability was demonstrated for all variables across all ages. The genes involved did not differ between the sexes, although there was evidence for sex limitation in the case of haemoglobin at age twelve. Longitudinal analyses indicated that new genes affecting red cell indices were expressed at different stages of puberty. Some of these genes affected the different red cell indices pleiotropically, while others had effects specific to one variable only.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

Parasite distribution on client reef fish determines cleaner fish foraging patterns

Recent evidence suggests that cleaner fish Labroides dimidiatus effectively control parasite densities on client reef fish that actively visit them to have parasites and dead or infected tissue removed. These findings support the hypothesis that clients benefit from cleaning, However, they do not show how cleaners reduce the parasite load of their clients. Cleaners could selectively feed on parasites or parasite removal could be a side product of cleaners foraging indifferently on the client surface, resulting in the removal of healthy mucus and scales also. To investigate cleaner fish foraging behaviour, we infected individuals of the surgeon fish Ctenochaetus striatus, with parasitic monogeneans on one body side, while the other body side was parasite free. We then allowed these clients to interact with L, dimidiatus. We found that the duration of interactions depended on parasite load, and that cleaners spent both more time and took more bites per time unit on the infected than on the uninfected side, Our data thus support the idea that parasite abundance determines food patch quality for cleaners. The overall outcome of cleaning interactions is thus likely to benefit the clients.