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.

Existence of multiple solutions for finite difference approximations to second-order boundary value problems

Error condition detected We consider discrete two-point boundary value problems of the form D-2 y(k+1) = f (kh, y(k), D y(k)), for k = 1,...,n - 1, (0,0) = G((y(0),y(n));(Dy-1,Dy-n)), where Dy-k = (y(k) - Yk-I)/h and h = 1/n. This arises as a finite difference approximation to y" = f(x,y,y'), x is an element of [0,1], (0,0) = G((y(0),y(1));(y'(0),y'(1))). We assume that f and G = (g(0), g(1)) are continuous and fully nonlinear, that there exist pairs of strict lower and strict upper solutions for the continuous problem, and that f and G satisfy additional assumptions that are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. Under these assumptions we show that there are at least three distinct solutions of the discrete approximation which approximate solutions to the continuous problem as the grid size, h, goes to 0. (C) 2003 Elsevier Science Ltd. All rights reserved.

Product form approximations for highly linear loss networks with trunk reservation

Admission controls, such as trunk reservation, are often used in loss networks to optimise their performance. Since the numerical evaluation of performance measures is complex, much attention has been given to finding approximation methods. The Erlang Fixed-Point (EFP) approximation, which is based on an independent blocking assumption, has been used for networks both with and without controls. Several more elaborate approximation methods which account for dependencies in blocking behaviour have been developed for the uncontrolled setting. This paper is an exploratory investigation of extensions and synthesis of these methods to systems with controls, in particular, trunk reservation. In order to isolate the dependency factor, we restrict our attention to a highly linear network. We will compare the performance of the resulting approximations against the benchmark of the EFP approximation extended to the trunk reservation setting. By doing this, we seek to gain insight into the critical factors in constructing an effective approximation. (C) 2003 Elsevier Ltd. All rights reserved.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

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.

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.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Eigenvalue computations in the context of data-sparse approximations of integral operators

In this work, we consider the numerical solution of a large eigenvalue problem resulting from a finite rank discretization of an integral operator. We are interested in computing a few eigenpairs, with an iterative method, so a matrix representation that allows for fast matrix-vector products is required. Hierarchical matrices are appropriate for this setting, and also provide cheap LU decompositions required in the spectral transformation technique. We illustrate the use of freely available software tools to address the problem, in particular SLEPc for the eigensolvers and HLib for the construction of H-matrices. The numerical tests are performed using an astrophysics application. Results show the benefits of the data-sparse representation compared to standard storage schemes, in terms of computational cost as well as memory requirements.

Error bounds for low-rank approximations of the first exponential integral kernel

A hierarchical matrix is an efficient data-sparse representation of a matrix, especially useful for large dimensional problems. It consists of low-rank subblocks leading to low memory requirements as well as inexpensive computational costs. In this work, we discuss the use of the hierarchical matrix technique in the numerical solution of a large scale eigenvalue problem arising from a finite rank discretization of an integral operator. The operator is of convolution type, it is defined through the first exponential-integral function and, hence, it is weakly singular. We develop analytical expressions for the approximate degenerate kernels and deduce error upper bounds for these approximations. Some computational results illustrating the efficiency and robustness of the approach are presented.

Modeling the Diffusion of Energy Performance Contracting

The Portuguese Energy Policy considers the development of a commercially viable and competitive market for energy performance contracting (EPC) as a main mechanism to achieve the objectives of energy efficiency improvement. This paper proposes a study to investigate how to achieve widespread adoption of energy performance contracting by means of system dynamics modelling and simulation. To explore and gather insights on this question, a system dynamics model representing the system of the Portuguese EPC market at industry level will be created. The simulation of that model will provide a helpful basis for analysing and explaining the development of key variables, and accelerating learning on the managerial, organizational and political adaptation processes that foster the diffusion of EPC adoption. The first phase of this research project aims at identifying and analysing the key factors and critical cause-effect relations that drive the adoption of EPC. With this purpose, a qualitative content analysis on relevant documents was performed and a set of interviews was conducted. That data was analysed to capture the critical variables and its interrelation to formulate a preliminary representation of the system structure as stock and flow diagrams.

A class of singular first order differential equations with applications in reaction-diffusion

Agências Financiadoras: FCT e MIUR

Influence of network effects on the diffusion of information systems: the case of open source

Doutoramento em Gestão

Effects of international diffusion of a general purpose technology on wage inequality

This paper studies the effects of the diffusion of a General Purpose Technology (GPT) that spreads first within the developed North country of its origin, and then to a developing South country. In the developed general equilibrium growth model, each final good can be produced by one of two technologies. Each technology is characterized by a specific labor complemented by a specific set of intermediate goods, which are enhanced periodically by Schumpeterian R&D activities. When quality reaches a threshold level, a GPT arises in one of the technologies and spreads first to the other technology within the North. Then, it propagates to the South, following a similar sequence. Since diffusion is not even, neither intra- nor inter-country, the GPT produces successive changes in the direction of technological knowledge and in inter- and intra-country wage inequality. Through this mechanism the different observed paths of wage inequality can be accommodated.