A compartmental model of hepatic disposition kinetics: 1. Model development and application to linear kinetics

The conventional convection-dispersion model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. The extension of this model to include nonlinear kinetics and zonal heterogeneity of the liver is not straightforward and requires numerical solution of partial differential equation, which is not available in standard nonlinear regression analysis software. In this paper, we describe an alternative compartmental model representation of hepatic disposition (including elimination). The model allows the use of standard software for data analysis and accurately describes the outflow concentration-time profile for a vascular marker after bolus injection into the liver. In an evaluation of a number of different compartmental models, the most accurate model required eight vascular compartments, two of them with back mixing. In addition, the model includes two adjacent secondary vascular compartments to describe the tail section of the concentration-time profile for a reference marker. The model has the added flexibility of being easy to modify to model various enzyme distributions and nonlinear elimination. Model predictions of F, MTT, CV2, and concentration-time profile as well as parameter estimates for experimental data of an eliminated solute (palmitate) are comparable to those for the extended convection-dispersion model.

Catheter effects in organ perfusion experiments

In a typical isolated organ perfusion experiment, a substance is injected upstream of an organ and then collected at some distance downstream. To reach the organ from the injection site, and then from the organ to the collector, a solute passes through catheters, usually tubes with circular cross-sections. Catheters cause distortion to the concentration-time profile of the perfusion. In this paper, we analyse catheter distribution kinetics from a mathematical point of view, develop the function most suitable for modeling this distribution and successfully apply this function to experimental data. (C) 2002 Academic Press.

Air-water flows down stepped chutes: turbulence and flow structure observations

Interactions between turbulent waters and atmosphere may lead to strong air-water mixing. This experimental study is focused on the flow down a staircase channel characterised by very strong flow aeration and turbulence. Interfacial aeration is characterised by strong air-water mixing extending down to the invert. The size of entrained bubbles and droplets extends over several orders of magnitude, and a significant number of bubble/droplet clusters was observed. Velocity and turbulence intensity measurements suggest high levels of turbulence across the entire air-water flow. The increase in turbulence levels, compared to single-phase flow situations, is proportional to the number of entrained particles. (C) 2002 Elsevier Science Ltd. All rights reserved.

Population balance modelling of granulation with a physically based coalescence kernel

It was previously published by the authors that granules can either coalesce through Type I (when granules coalesce by viscous dissipation in the surface liquid layer before their surfaces touch) or Type II (when granules are slowed to a halt during rebound, after their surfaces have made contact) (AIChE J. 46 (3) (2000) 529). Based on this coalescence mechanism, a new coalescence kernel for population balance modelling of granule growth is presented. The kernel is constant such that only collisions satisfying the conditions for one of the two coalescence types are successful. One constant rate is assigned to each type of coalescence and zero is for the case of rebound. As the conditions for Types I and II coalescence are dependent on granule and binder properties, the coalescence kernel is thus physically based. Simulation results of a variety of binder and granule materials show good agreement with experimental data. (C) 2002 Elsevier Science Ltd. All rights reserved.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

The martensitic transformation in ceramics - its role in transformation toughening

This paper reviews the current knowledge and understanding of martensitic transformations in ceramics - the tetragonal to monoclinic transformation in zirconia in particular. This martensitic transformation is the key to transformation toughening in zirconia ceramics. A very considerable body of experimental data on the characteristics of this transformation is now available. In addition, theoretical predictions can be made using the phenomenological theory of martensitic transformations. As the paper will illustrate, the phenomenological theory is capable of explaining all the reported microstructural and crystallographic features of the transformation in zirconia and in some other ceramic systems. Hence the theory, supported by experiment, can be used with considerable confidence to provide the quantitative data that is essential for developing a credible, comprehensive understanding of the transformation toughening process. A critical feature in transformation toughening is the shape strain that accompanies the transformation. This shape strain, or nucleation strain, determines whether or not the stress-induced martensitic transformation can occur at the tip of a potentially dangerous crack. If transformation does take place, then it is the net transformation strain left behind in the transformed region that provides toughening by hindering crack growth. The fracture mechanics based models for transformation toughening, therefore, depend on having a full understanding of the characteristics of the martensitic transformation and, in particular, on being able to specify both these strains. A review of the development of the models for transformation toughening shows that their refinement and improvement over the last couple of decades has been largely a result of the inclusion of more of the characteristics of the stress-induced martensitic transformation. The paper advances an improved model for the stress-induced martensitic transformation and the strains resulting from the transformation. This model, which separates the nucleation strain from the subsequent net transformation strain, is shown to be superior to any of the constitutive models currently available. (C) 2002 Elsevier Science Ltd. All rights reserved.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Percolation Phenomena in Micropore Adsorption: Influence on Single and Multicomponent Adsorption Equilibria

A novel and simple method for determination of micropore network connectivity of activated carbon using liquid phase adsorption is presented in this paper. The method is applied to three different commercial carbons with eight different liquid phase adsorptives as probes. The effect of the pore network connectivity on the prediction of multicomponent adsorption equilibria was also studied. For this purpose, the Ideal Adsorbed Solution Theory (IAST) was used in conjuction with the modified DR single component isotherm. The results of comparison with experimental data show that incorporation of the connectivity, and consideration of percolation processes associated with the different molecular sizes of the adsorptives in the mixture, can improve the performance of the IAST in predicting multicomponent adsorption equilibria.

Kinetics of adsorption on activated carbon: application of heterogeneous vacancy solution theory

The kinetics of single component adsorption on activated carbon is investigated here using a heterogeneous vacancy solution theory (VST) of adsorption. The adsorption isotherm is developed to account for the adsorbate non-ideality due to the size difference between the adsorbate molecule and the vacant site, while incorporating adsorbent heterogeneity through a pore-width-related potential energy. The transport process in the bidisperse carbon considers coupled mass transfer in both macropore and micropore phases simultaneously. Adsorbate diffusion in the micropore network is modeled through effective medium theory, thus considering pore network connectivity in the adsorbent, with the activation energy for adsorbate diffusion related to the adsorption energy, represented by the Steele 10-4-3 potential for carbons. Experimental data of five hydrocarbons, CO2 and SO2 on Ajax carbon at multiple temperatures, as well as three hydrocarbons on Norit carbon at three temperatures are first fitted by the heterogeneous VST model to obtain the isotherm parameters, followed by application of the kinetic model to uptake data on carbon particles of different sizes and geometry at various temperatures. For the hydrocarbons studied, the model can successfully correlate the experimental data for both adsorption equilibrium and kinetics. However, there is some deviation in the fit of the desorption kinetics for polar compounds such as CO2 and SO2, due to the inadequacy of the L-J potential model in this case. The significance of viscous transport in the micropores is also considered here and found to be negligible, consistent with recent molecular simulation studies. (C) 2002 Elsevier Science Ltd. All rights reserved.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

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

Characterization of micro-mesoporous carbonaceous materials. Calculations of adsorption isotherm of hydrocarbons

In this paper we apply a method recently developed by Do and co-workers(1) for the prediction of adsorption isotherms of pure vapors on carbonaceous materials. The information required for the prediction is the pore size distribution and the BET constant, C, of a corresponding nonporous surface (graphite). The dispersive adsorption force is assumed to be the dominant force in adsorption mechanism. This applies to nonpolar and weakly polar hydrocarbons. We test this predictive model against the adsorption data of benzene, toluene, n-pentane, n-hexane, and ethanol on a commercial activated carbon. It is found that the predictions are excellent for all adsorbates tested with the exception of ethanol where the predicted values are about 10% less than the experimental data, and this is probably attributed to the electrostatic interaction between ethanol molecules and the functional groups on the carbon surfaces.

Adsorption of mixtures containing subcritical fluids on activated carbon

A theoretical analysis of adsorption of mixtures containing subcritical adsorbates into activated carbon is presented as an extension to the theory for pure component developed earlier by Do and coworkers. In this theory, adsorption of mixtures in a pore follows a two-stage process, similar to that for pure component systems. The first stage is the layering of molecules on the surface, with the behavior of the second and higher layers resembling to that of vapor-liquid equilibrium. The second stage is the pore-filling process when the remaining pore width is small enough and the pressure is high enough to promote the pore filling with liquid mixture having the same compositions as those of the outermost molecular layer just prior to pore filling. The Kelvin equation is applied for mixtures, with the vapor pressure term being replaced by the equilibrium pressure at the compositions of the outermost layer of the liquid film. Simulations are detailed to illustrate the effects of various parameters, and the theory is tested with a number of experimental data on mixture. The predictions were very satisfactory.

Modelling the activated sludge flocculation process combining laser light diffraction particle sizing and population balance modelling (PBM)

A technique based on laser light diffraction is shown to be successful in collecting on-line experimental data. Time series of floc size distributions (FSD) under different shear rates (G) and calcium additions were collected. The steady state mass mean diameter decreased with increasing shear rate G and increased when calcium additions exceeded 8 mg/l. A so-called population balance model (PBM) was used to describe the experimental data, This kind of model describes both aggregation and breakage through birth and death terms. A discretised PBM was used since analytical solutions of the integro-partial differential equations are non-existing. Despite the complexity of the model, only 2 parameters need to be estimated: the aggregation rate and the breakage rate. The model seems, however, to lack flexibility. Also, the description of the floc size distribution (FSD) in time is not accurate.

Indirect evidence of density-dependent population regulation in Aponomma hydrosauri (Acari : Ixodidae), an ectoparasite of reptiles

The extent to which density-dependent processes regulate natural populations is the subject of an ongoing debate. We contribute evidence to this debate showing that density-dependent processes influence the population dynamics of the ectoparasite Aponomma hydrosauri (Acari: Ixodidae), a tick species that infests reptiles in Australia. The first piece of evidence comes from an unusually long-term dataset on the distribution of ticks among individual hosts. If density-dependent processes are influencing either host mortality or vital rates of the parasite population, and those distributions can be approximated with negative binomial distributions, then general host-parasite models predict that the aggregation coefficient of the parasite distribution will increase with the average intensity of infections. We fit negative binomial distributions to the frequency distributions of ticks on hosts, and find that the estimated aggregation coefficient k increases with increasing average tick density. This pattern indirectly implies that one or more vital rates of the tick population must be changing with increasing tick density, because mortality rates of the tick's main host, the sleepy lizard, Tiliqua rugosa, are unaffected by changes in tick burdens. Our second piece of evidence is a re-analysis of experimental data on the attachment success of individual ticks to lizard hosts using generalized linear modelling. The probability of successful engorgement decreases with increasing numbers of ticks attached to a host. This is direct evidence of a density-dependent process that could lead to an increase in the aggregation coefficient of tick distributions described earlier. The population-scale increase in the aggregation coefficient is indirect evidence of a density-dependent process or processes sufficiently strong to produce a population-wide pattern, and thus also likely to influence population regulation. The direct observation of a density-dependent process is evidence of at least part of the responsible mechanism.