Dissolution of a solid sphere in an unbounded, stagnant liquid

An exact analytical solution is obtained for the transient dissolution of solid spheres in a diffusion-controlled environment. This result provides a useful reference point for drug testing in humans. The dimensionless solution is expressed in terms of a single parameter, which accounts for solubility, bulk flow, and stagnant fluid composition. A simple, explicit and exact expression was found to predict time-to-complete dissolution (TCD). An approximate solution was also found which tracks the exact case for low solubility conditions.

Enzyme immobilization in porous solid supports - penetration of immobilized enzyme

The process of enzyme immobilization under the diffusion-controlled regime (i.e., fast attachment of enzyme compared to its diffusion) is modeled and theoretically solved in this article. Simple and compact solutions for the penetration depth of immobilized enzyme and the bulk enzyme concentration versus time are presented. Furthermore, the conditions for the validity of our solutions are also given in this article so that researchers can discover when the theoretical solutions can be applied to their systems.

Structural ordering of coal char during heat treatment and its impact on reactivity

The effect of heat treatment on the structure of an Australian semi-anthracite char was studied in detail in the 850-1150degreesC temperature range using XRD, HRTEM, and electrical resistivity techniques. It was found that the carbon crystallite size in the char does not change significantly during heat treatment in the temperature range studied, for both the raw coal and its ash-free derivative obtained by acid treatment. However, the fraction of the organized carbon in the raw coal chars, determined by XRD, increased with increase of heat treatment time and temperature, while that for the ash-free coal chars remained almost unchanged. This suggests the occurrence of catalytic ordering during heat treatment, supported by the observation that the electrical resistivity of the raw coal chars decreased with heat treatment, while that of the ash-free coal chars did not vary significantly. Further confirmatory evidence was provided by high resolution transmission electron micrographs depicting well-organized carbon layers surrounding iron particles. It is also found that the fraction of organized carbon does not reach unity, but attains an apparent equilibrium value that increases with increase in temperature, providing an apparent heat of ordering of 71.7 kJ mol(-1) in the temperature range studied. Good temperature-independent correlation was found between the electrical resistivity and the organized carbon fraction, indicating that electrical resistivity is indeed structure sensitive. Good correlation was also found between the electrical resistivity and the reactivity of coal char. All these results strongly suggest that the thermal deactivation is the result of a crystallite-perfecting process, which is effectively catalyzed by the inorganic matter in the coal char. Based on kinetic interpretation of the data it is concluded that the process is diffusion controlled, most likely involving transport of iron in the inter-crystallite nanospaces in the temperature range studied. The activation energy of this transport process is found to be very low, at about 11.8 kJ mol(-1), which is corroborated by model-free correlation of the temporal variation of organized carbon fraction as well as electrical resistivity data using the superposition method, and is suggestive of surface transport of iron. (C) 2002 Elsevier Science Ltd. All rights reserved.

Comments on edge-to-edge matching and the equivalence of the invariant line, Delta g and Moire fringe approaches to the crystallographic features of precipitates

Three apparently distinct and different approaches have been proposed to account for the crystallographic features of diffusion-controlled precipitation. These three models are based on (a) an invariant line in the habit plane, (b) the parallelism of a pair of Deltags that are perpendicular to the habit plane and (c) the parallelism of a pair of Moire fringes that are in turn parallel to the habit plane. The purpose of the present paper is to show that these approaches are in fact absolutely equivalent and that when certain conditions are satisfied they are essentially the same as the recent edge-to-edge matching model put forward by the authors. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Edge-to-edge matching and its applications - Part I. Application to the simple HCP/BCC system

The edge-to-edge matching crystallographic model has been used to predict all the orientation relationships (OR) between crystals that have simple hexagonal close packed (HCP) and body-centered cubic (BCC) structures. Using the critical values for the interatomic spacing misfit along the matching directions and the cl-value mismatch between matching planes, the model predicted all the four common ORs, namely the Burgers OR, the Potter OR, the Pitsch-Schrader OR and the Rong Dunlop OR, together with the corresponding habit planes. Taking the c(H)/a(H) and a(H)/a(B) ratios as variables, where H and B denote the HCP and BCC structures respectively, the model also predicted the relationship between these variables and the four ORs. These predictions are perfectly consistent with the published experimental results. As was the case in the FCC/BCC system, the edge-to-edge matching model has been shown to be a powerful tool for predicting the crystallographic features of diffusion-controlled phase transformations. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Edge-to-edge matching and its applications - Part II. Application to Mg-Al, Mg-Y and Mg-Mn alloys

A model for the crystallography and morphology of diffusion-controlled phase transformations - edge-to-edge matching - has been used to predict the orientation relationships (OR) and habit planes of precipitates Mg17Al12 in Mg-Al alloy, Mg24Y5 in Mg-Y alloy and alpha-Mn in Mg-Mn alloy. Based on the crystal structures and lattice parameters only, the model predicts that the possible ORs between Mg17Al12 and Mg matrix are the near Burgers OR, the Potter OR, the Gjonnes-Ostmoe OR and the Crawley OR. In the Mg-Y alloy, the OR between Mg24Y5 precipitates and the Mg matrix is predicted to be the Burgers OR only. The model also predicts that there are no reproducible ORs between alpha-Mn and Mg in the Mg-Mn alloy. Combining the edge-to-edge matching model and W. Zhang's Deltag approach, the habit plane and side facets of the precipitate for each OR can be determined. All the predicted ORs and the corresponding habit planes in Mg-Al and Mg-Y alloys agree very well with the experimental results. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Analysis of multiple gas solid reactions

Multiple gas solid reactions involving one solid and N gaseous reactants are investigated in this study by using a matched asymptotic expansion technique. Two cases are particularly studied. In the first case all N chemical reaction rates are faster than the diffusion rate. While in the second case only M (M < N) chemical reaction rates are faster than the diffusion rate and the rates of the remaining (N-M) chemical reactions are comparable to that of diffusion. For these two cases the solid concentration profile behaves like a travelling wave. In the first case the wave front velocity is contributed linearly by all gaseous reactants (additive law) while in the second case this law does not hold.

Oxidation mechanism of Si3N4-bonded SiC ceramics by CO, CO2 and steam

This paper presents a theoretical and experimental investigation into the oxidation reactions of Si3N4-bonded SiC ceramics. Such ceramics which contain a small amount of silicon offer increased oxidation and wear resistance and are widely used as lining refractories in blast furnaces. The thermodynamics of oxidation reactions were studied using the JANAF tables. The weight gain was measured using a thermogravimetric analysis technique to study the kinetics. The temperature range of oxidation measurements is from 1073 to 1573 K and the oxidation atmosphere is water vapour, pure CO and CO-CO2 gas mixtures with various CO-to-CO2 ratios. Thermodynamic simulations showed that the oxidation mechanism of Si3N4-bonded SiC ceramics is passive oxidation and all components contribute to the formation of a silica film. The activated energies of the reactions follow the sequence Si3N4>SiC>Si. The kinetic study revealed that the oxidation of Si3N4-bonded SiC ceramics occurred in a mixed regime controlled by both interface reaction and diffusion through the silica film. Under the atmosphere conditions prevailing in the blast furnace, this ceramic is predicted to be passively oxidized with the chemical reaction rate becoming more dominant as the CO concentration increases. (C) 1998 Chapman & Hall.

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.

Diffusion- and perfusion-weighted MRI response to thrombolysis in stroke

Diffusion- and perfusion-weighted magnetic resonance imaging provides important pathophysiological information in acute bra-in ischemia. We performed a prospective study in 19 sub-6-hour stroke patients using serial diffusion- and perfusion-weighted imaging before intravenous thrombolysis, with repeat studies, both subacutely and at outcome. For comparison of ischemic lesion evolution and clinical outcome, we used a historical control group of 21 sub-6-hour ischemic stroke patients studied serially with diffusion- and perfusion-weighted imaging. The two groups were well matched for the baseline National Institutes of Health Stroke Scale and magnetic resonance parameters. Perfusion-weighted imaging-diffusion-weighted imaging mismatch was present in 16 of 19 patients treated with tissue plasminogen activator, and 16 of 21 controls. Perfusion-weighted imaging-diffusion-weighted imaging mismatch patients treated with tissue plaminogen activator had higher recanalization rates and enhanced reperfusion at day 3 (81% vs 47% in controls), and a greater proportion of severely hypoperfused acute mismatch tissue not progressing to infarction (82% vs -25% in controls). Despite similar baseline diffusion-weighted imaging lesions, infarct expansion was less in the recombinant tissue plaminogen activator group (14cm(3) vs 56cm(3) in controls). The positive effect of thrombolysis on lesion growth in mismatch patients translated into a greater improvement in baseline to outcome National Institutes of Health Stroke Scale in the group treated with recombinant tissue plaminogen activator, and a significantly larger proportion of patients treated with recombinant tissue plaminogen activator having a clinically meaningful improvement in National Institutes of Health Stroke Scale of;2:7 points. The natural evolution of acute perfusion-weighted imaging-diffusion-weighted imaging mismatch tissue may be altered by thrombolysis, with improved stroke outcome. This has implications for the use of diffusion- and perfusion-weighted imaging in selecting and monitoring patients for thrombolytic therapy.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

A flux ratio theorem for multicomponent linear diffusion equations

Ussing [1] considered the steady flux of a single chemical component diffusing through a membrane under the influence of chemical potentials and derived from his linear model, an expression for the ratio of this flux and that of the complementary experiment in which the boundary conditions were interchanged. Here, an extension of Ussing's flux ratio theorem is obtained for n chemically interacting components governed by a linear system of diffusion-migration equations that may also incorporate linear temporary trapping reactions. The determinants of the output flux matrices for complementary experiments are shown to satisfy an Ussing flux ratio formula for steady state conditions of the same form as for the well-known one-component case. (C) 2000 Elsevier Science Ltd. All rights reserved.

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.

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.