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.

Diffusion of n-decane in mesoporous MCM-41 silicas

We investigate here the diffusion of n-decane in nanoporous MCM-41 silicas with pore diameters between 3.0 and 4.3 nm, and at various temperatures and purge flow rates, by the Zero Length Column method. A complete-time-range analysis of desorption curves is used to derive the diffusion coefficient, and the effect of pore size, purge flow rate and temperature on the diffusion character is systematically studied. The results show that the calculated low-coverage diffusivity values are strongly dependent on temperature but only weakly dependent on pore size. The study reveals that transport is controlled by intracrystalline diffusion and dominated by sorbate-sorbent interaction, with the experimental isosteric heat matching the potential energy of flat-lying n-decane molecules on the surface, determined using a united atom model. The diffusion activation energy and adsorption isosteric heat at zero loading for the different pore size MCM-41 samples vary in a narrow range respectively, and their ratio is essentially constant over the pore size range studied. The study shows that the ZLC method is an effective tool to investigate the diffusion kinetics of hydrocarbons in mesoporous MCM-41 materials. (c) 2005 Elsevier Inc. All rights reserved.

Development of a non-intrusive heat transfer coefficient gauge and its application to high pressure die casting Effect of the process parameters

A heat transfer coefficient gauge has been built, obeying particular rules in order to ensure the relevance and accuracy of the collected information. The gauge body is made out of the same materials as the die casting die (H13). It is equipped with six thermocouples located at different depths in the body and with a sapphire light pipe. The light pipe is linked to an optic fibre, which is connected to a monochromatic pyrometer. Thermocouples and pyrometer measurements are recorded with a data logger. A high pressure die casting die was instrumented with one such gauge. A set of 150 castings was done and the data recorded. During the casting, some process parameters have been modified such as piston velocity, intensification pressure, delay before switch to the intensification stage, temperature of the alloy, etc.... The data was treated with an inverse method in order to transform temperature measurements into heat flux density and heat transfer coefficient plots. The piston velocity and the initial temperature of the die seem to be the process parameters that have the greatest influence on the heat transfer. (c) 2005 Elsevier B.V. All rights reserved.

Estimation of concentration-dependent diffusion coefficient in drying process from the space-averaged concentration versus time with experimental data

The estimation of a concentration-dependent diffusion coefficient in a drying process is known as an inverse coefficient problem. The solution is sought wherein the space-average concentration is known as function of time (mass loss monitoring). The problem is stated as the minimization of a functional and gradient-based algorithms are used to solve it. Many numerical and experimental examples that demonstrate the effectiveness of the proposed approach are presented. Thin slab drying was carried out in an isothermal drying chamber built in our laboratory. The diffusion coefficients of fructose obtained with the present method are compared with existing literature results.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Quantum noise in optical fibers. I. Stochastic equations

We analyze the quantum dynamics of radiation propagating in a single-mode optical fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We start from a fundamental Hamiltonian that includes the principal known nonlinear effects and quantum-noise sources, including linear gain and loss. Both Markovian and frequency-dependent, non-Markovian reservoirs are treated. This treatment allows quantum Langevin equations, which have a classical form except for additional quantum-noise terms, to be calculated. In practical calculations, it is more useful to transform to Wigner or 1P quasi-probability operator representations. These transformations result in stochastic equations that can be analyzed by use of perturbation theory or exact numerical techniques. The results have applications to fiber-optics communications, networking, and sensor technology.

Exact solution of the A(n-1)((1)) trigonometric vertex model with non-diagonal open boundaries

The A(n-1)((1)) trigonometric vertex model with generic non-diagonal boundaries is studied. The double-row transfer matrix of the model is diagonalized by algebraic Bethe ansatz method in terms of the intertwiner and the corresponding face-vertex relation. The eigenvalues and the corresponding Bethe ansatz equations are obtained.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Water diffusion in hydroxyethyl methacrylate (HEMA)-based hydrogels formed by gamma-radiolysis

Polymer hydrogels based upon methacrylates are used extensively in the pharmaceutical industry, particularly as controlled release drug delivery systems. These materials are generally prepared by chemically initiated polymerization, but this can lead to the presence of unwanted initiator fragments in the polymer matrix. In the present work, initiation of polymerization by gamma-irradiation of hydroxyethyl methacrylate, with and without added crosslinkers, has been investigated, and the diffusion coefficients for water in the resulting polymers have been measured through mass uptake by the polymers. The diffusion of water in poly(hydroxyethyl methacrylate) at 310 K was found to be Fickian, with a diffusion coefficient of 1.96 +/- 0.1 x 10(11) m(2) s(-1) and an equilibrium water content of 58%, NMR imaging analyses confirmed the adherance to a Fickian model of the diffusion of water into polymer cylinders. The incorporation of small amounts (0.2-0.5 wt%) of added ethyleneglycol-dimethacrylate-based crosslinkers was found to have only a small effect on the diffusion coefficient and the equilibrium water content for the copolymers. (C) 1999 Society of Chemical Industry.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

C-13-NMR, H-1-NMR, and FT-Raman study of radiation-induced modifications in radiation dosimetry polymer gels

H-1- and C-13-NMR spectroscopy and FT-Raman spectroscopy are used to investigate the properties of a polymer gel dosimeter post-irradiation. The polymer gel (PACT) is composed of acrylamide, N,N'-methylene-bisacrylamide, gelatin, and water. The formation of a polyacrylamide network within the gelatin matrix follows a dose dependence nonlinearly correlated to the disappearance of the double bonds from the dissolved monomers within the absorbed dose range of 0-50 Gy. The signal from the gelatin remains constant with irradiation. We show that the NMR spin-spin relaxation times (T-2) of PAGs irradiated to up to 50 Gy measured in a NMR spectrometer and a clinical magnetic resonance imaging scanner can be modeled using the spectroscopic intensity of the growing polymer network. More specifically, we show that the nonlinear T-2 dependence against dose can be understood in terms of the fraction of protons in three different proton pools. (C) 2000 John Wiley & Sons, Inc.

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.