Mathematical models in percutaneous absorption

A number of mathematical models have been used to describe percutaneous absorption kinetics. In general, most of these models have used either diffusion-based or compartmental equations. The object of any mathematical model is to a) be able to represent the processes associated with absorption accurately, b) be able to describe/summarize experimental data with parametric equations or moments, and c) predict kinetics under varying conditions. However, in describing the processes involved, some developed models often suffer from being of too complex a form to be practically useful. In this chapter, we attempt to approach the issue of mathematical modeling in percutaneous absorption from four perspectives. These are to a) describe simple practical models, b) provide an overview of the more complex models, c) summarize some of the more important/useful models used to date, and d) examine sonic practical applications of the models. The range of processes involved in percutaneous absorption and considered in developing the mathematical models in this chapter is shown in Fig. 1. We initially address in vitro skin diffusion models and consider a) constant donor concentration and receptor conditions, b) the corresponding flux, donor, skin, and receptor amount-time profiles for solutions, and c) amount- and flux-time profiles when the donor phase is removed. More complex issues, such as finite-volume donor phase, finite-volume receptor phase, the presence of an efflux. rate constant at the membrane-receptor interphase, and two-layer diffusion, are then considered. We then look at specific models and issues concerned with a) release from topical products, b) use of compartmental models as alternatives to diffusion models, c) concentration-dependent absorption, d) modeling of skin metabolism, e) role of solute-skin-vehicle interactions, f) effects of vehicle loss, a) shunt transport, and h) in vivo diffusion, compartmental, physiological, and deconvolution models. We conclude by examining topics such as a) deep tissue penetration, b) pharmacodynamics, c) iontophoresis, d) sonophoresis, and e) pitfalls in modeling.

Coalbed methane sorption related to coal composition

Gas sorption by coal is closely related to its physical and chemical properties, which are, in turn, governed by coal type and rank. The role of coal type (sensu maceral composition) is not fully established but it is clear that coal type may affect both adsorption capacity and desorption rate. Adsorption capacity is closely related to micropore (pores <2 nm) development, which is rank and maceral dependent. Adsorption isotherms indicate that in most cases bright (vitrinite-rich) coals have a greater adsorption capacity than their dull (often inertinite-rich) equivalents. However, no differences, or even the opposing trend, may be observed in relation to coal type. Desorption rate investigations have been performed using selected bright and dull coal samples in a high pressure microbalance. Interpretation of results using unipore spherical and bidisperse pore models indicate the importance of the pore structure. Bright, vitrinite-rich coals usually have the slowest desorption rates which is associated with their highly microporous structure. However, rapid desorption in bright coals may be related to development of extensive, unmineralised fracture systems. Both macro-and micro-pore systems are implicated in the more rapidly desorbing dull coals. Some dull, inertinite-rich coals may rapidly desorb due to a predominance of large, open cell lumina. Mineral matter is essentially nonadsorbent to coal gases and acts as a simple diluent. However, mineral-rich coals may be associated with more rapid desorption. Coal rank and type (maceral composition) per se do not appear to be the critical factors in controlling gas sorption, but rather the influence they exert over pore structure development. (C) 1998 Elsevier Science B.V.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

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.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

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.

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.

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.

On the robustness of short-term interest rate models

This paper investigates the robustness of a range of short–term interest rate models. We examine the robustness of these models over different data sets, time periods, sampling frequencies, and estimation techniques. We examine a range of popular one–factor models that allow the conditional mean (drift) and conditional variance (diffusion) to be functions of the current short rate. We find that parameter estimates are highly sensitive to all of these factors in the eight countries that we examine. Since parameter estimates are not robust, these models should be used with caution in practice.

Surface diffusion of strongly adsorbing vapors in activated carbon by a differential permeation method

Conventional methods to determine surface diffusion of adsorbed molecules are proven to be inadequate for strongly adsorbing vapors on activated carbon. Knudsen diffusion permeability (B-k) for strongly adsorbing vapors cannot be directly estimated from that of inert gases such as helium. In this paper three models are considered to elucidate the mechanism of surface diffusion in activated carbon. The transport mechanism in all three models is a combination of Knudsen diffusion, viscous flow and surface diffusion. The collision reflection factor f (which is the fraction of molecules undergoing collision to the solid surface over reflection from the surface) of the Knudsen diffusivity is assumed to be a function of loading. It was found to be 1.79 in the limit of zero loading, and decreases as loading increases. The surface diffusion permeability increases sharply at very low pressures and then starts to decrease after it has reached a maximum (B(mum)s) at a threshold pressure. The initial rapid increase in the total permeability is mainly attributed to surface diffusion. Interestingly the B(mum)s for all adsorbates appear at the same volumetric adsorbed phase concentration, suggesting that the volume of adsorbed molecules may play an important role in the surface diffusion mechanism in activated carbon. (C) 2003 Elsevier Ltd. All rights reserved.

Diffusion modeling of percutaneous absorption kinetics: 3. Variable diffusion and partition coefficients, consequences for stratum corneum depth profiles and desorption kinetics

Stratum corneum (SC) desorption experiments have yielded higher calculated steady-state fluxes than those obtained by epidermal penetration studies. A possible explanation of this result is a variable diffusion or partition coefficient across the SC. We therefore developed the diffusion model for percutaneous penetration and desorption to study the effects of either a variable diffusion coefficient or variable partition coefficient in the SC over the diffusion path length. Steady-state flux, lag time, and mean desorption time were obtained from Laplace domain solutions. Numerical inversion of the Laplace domain solutions was used for simulations of solute concentration-distance and amount penetrated (desorbed)-time profiles. Diffusion and partition coefficients heterogeneity were examined using six different models. The effect of heterogeneity on predicted flux from desorption studies was compared with that obtained in permeation studies. Partition coefficient heterogeneity had a more profound effect on predicted fluxes than diffusion coefficient heterogeneity. Concentration-distance profiles show even larger dependence on heterogeneity, which is consistent with experimental tape-stripping data reported for clobetasol propionate and other solutes. The clobetasol propionate tape-stripping data were most consistent with the partition coefficient decreasing exponentially for half the SC and then becoming a constant for the remaining SC. (C) 2004 Wiley-Liss, Inc.

Modelling oxygen diffusion and cell growth in a porous, vascularising scaffold for soft tissue engineering applications

Soft tissue engineering presents significant challenges compared to other tissue engineering disciplines such as bone, cartilage or skin engineering. The very high cell density in most soft tissues, often combined with large implant dimensions, means that the supply of oxygen is a critical factor in the success or failure of a soft tissue scaffold. A model is presented for oxygen diffusion in a 15-60 mm diameter dome-shaped scaffold fed by a blood vessel loop at its base. This model incorporates simple models for vascular growth, cell migration and the effect of cell density on the effective oxygen diffusivity. The model shows that the dynamic, homogeneous cell seeding method often employed in small-scale applications is not applicable in the case of larger scale scaffolds such as these. Instead, we propose the implantation of a small biopsy of tissue close to a blood supply within the scaffold as a technique more likely to be successful. Crown Copyright (c) 2005 Published by Elsevier Ltd. All rights reserved.