Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

On the validity of the shrinking core model in noncatalytic gas solid reaction

By using a matched asymptotic expansion technique, the shrinking core model (SCM) used in non-catalytic gas solid reactions with general kinetic expression is rigorously justified in this paper as a special case of the homogeneous model when the reaction rate is much faster than that of diffusion. The time-pendent velocity of the moving reacted-unreacted interface is found to be proportional to the gas flux at that interface for all geometries of solid particles, and the thickness order of the reaction zone and also the degree of chemical reaction at the interface is discussed in this paper.

Nmr imaging of the diffusion of water into poly(tetrahydrofurfuryl methacrylate-co-hydroxyethyl methacrylate)

0Nuclear magnetic resonance (n.m.r.) imaging was used to study the ingress of water into poly(tetrahydrofurfuryl methacrylate-co-hydroxyethyl methacrylate). The study offers strong evidence that the diffusion is Fickian in nature. The diffusion coefficient, D, obtained by fitting the underlying diffusion profile, attainable from the images, according to the equation for Fickian diffusion, is 1.5 x 10(-11) m(2) s(-1), which is in good correlation with the value of 2.1 x 10(-11) m(2) s(-1), obtained from mass uptake measurements. Additionally, from the T-2-weighted images, Superimposed features observed in addition to the underlying Fickian diffusion profiles were shown to have a longer spin-spin relaxation time, T-2. This Suggests the presence of two types of water within the polymer matrix; a less mobile phase of absorbed water that is interacting strongly with the polymer matrix and a more mobile phase of absorbed water residing within the cracks observed in the environmental scanning electron micrograph. (C) 1997 Elsevier Science Ltd.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.

PFG-NMR measurements of the self-diffusion coefficients of water in equilibrium poly(HEMA-co-THFMA) hydrogels

The self-diffusion coefficients for water in a series of copolymers of 2-hydroxyethyl methacrylate, HEMA, and tetrahydrofurfuryl methacrylate, THFMA, swollen with water to their equilibrium states have been studied at 310 K using PFG-NMR. The self-diffusion coefficients calculated from the Stejskal-Tanner equation, D-obs, for all of the hydrated polymers were found to be dependent on the NMR storage time, as a result of spin exchange between the proton reservoirs of the water and the polymers, reaching an equilibrium plateau value at long storage times. The true values of the diffusion coefficients were calculated from the values of D-obs, in the plateau regions by applying a correction for the fraction of water protons present, obtained from the equilibrium water contents of the gels. The true self-diffusion coefficient for water in polyHEMA obtained at 310 K by this method was 5.5 x 10(-10) m(2) s(-1). For the copolymers containing 20% HEMA or more a single value of the self-diffusion coefficient was found, which was somewhat larger than the corresponding values obtained for the macroscopic diffusion coefficient from sorption measurements. For polyTHFMA and copolymers containing less than 20% HEMA, the PFG-NMR stimulated echo attenuation decay curves and the log-attenuation plots were characteristic of the presence of two diffusing water species. The self-diffusion coefficients of water in the equilibrium-hydrated copolymers were found to be dependent on the copolymer composition, decreasing with increasing THFMA content.

Efficient time-independent wave packet scattering calculations within a Lanczos subspace: H+O-2 (J=0) state-to-state reaction probabilities

An efficient Lanczos subspace method has been devised for calculating state-to-state reaction probabilities. The method recasts the time-independent wave packet Lippmann-Schwinger equation [Kouri , Chem. Phys. Lett. 203, 166 (1993)] inside a tridiagonal (Lanczos) representation in which action of the causal Green's operator is affected easily with a QR algorithm. The method is designed to yield all state-to-state reaction probabilities from a given reactant-channel wave packet using a single Lanczos subspace; the spectral properties of the tridiagonal Hamiltonian allow calculations to be undertaken at arbitrary energies within the spectral range of the initial wave packet. The method is applied to a H+O-2 system (J=0), and the results indicate the approach is accurate and stable. (C) 2002 American Institute of Physics.

Prediction of absolute rate coefficients and product branching ratios for the C(P-3) plus allene reaction system

Complex chemical reactions in the gas phase can be decomposed into a network of elementary (e.g., unimolecular and bimolecular) steps which may involve multiple reactant channels, multiple intermediates, and multiple products. The modeling of such reactions involves describing the molecular species and their transformation by reaction at a detailed level. Here we focus on a detailed modeling of the C(P-3)+allene (C3H4) reaction, for which molecular beam experiments and theoretical calculations have previously been performed. In our previous calculations, product branching ratios for a nonrotating isomerizing unimolecular system were predicted. We extend the previous calculations to predict absolute unimolecular rate coefficients and branching ratios using microcanonical variational transition state theory (mu-VTST) with full energy and angular momentum resolution. Our calculation of the initial capture rate is facilitated by systematic ab initio potential energy surface calculations that describe the interaction potential between carbon and allene as a function of the angle of attack. Furthermore, the chemical kinetic scheme is enhanced to explicitly treat the entrance channels in terms of a predicted overall input flux and also to allow for the possibility of redissociation via the entrance channels. Thus, the computation of total bimolecular reaction rates and partial capture rates is now possible. (C) 2002 American Institute of Physics.

On the inverse parabolicity of pdf equations in turbulent flows

The focus of the present work is the well-known feature of the probability density function (PDF) transport equations in turbulent flows-the inverse parabolicity of the equations. While it is quite common in fluid mechanics to interpret equations with direct (forward-time) parabolicity as diffusive (or as a combination of diffusion, convection and reaction), the possibility of a similar interpretation for equations with inverse parabolicity is not clear. According to Einstein's point of view, a diffusion process is associated with the random walk of some physical or imaginary particles, which can be modelled by a Markov diffusion process. In the present paper it is shown that the Markov diffusion process directly associated with the PDF equation represents a reasonable model for dealing with the PDFs of scalars but it significantly underestimates the diffusion rate required to simulate turbulent dispersion when the velocity components are considered.

Morphological features of interfacial intermetallics and interfacial reaction rate in Al-11Si-2.5Cu-(0.15/0.60)Fe cast alloy/die steel couples

Soldering reactions are commonly observed during high pressure die casting of aluminium alloys, and involve the formation and growth of interfacial intermetallics between the die and the cast alloy. It is generally believed that close to 1% Fe is necessary in the aluminium alloy to reduce soldering. However, the role of iron in the interfacial reaction has not been studied in detail. In this investigation, reaction couples were formed between H13 tool steel substrates and an Al-11Si-2.5Cu melt containing either 0.15 or 0.60% Fe. Examination revealed distinctly different intermetallic layer morphology. The overall growth and chemistry of the reaction layer and the reaction rate measured by the consumption of the substrate were compared for the two alloy melts. It was demonstrated that a higher iron content reduces the rate of interfacial reaction, consistent with an observed thicker compact ( solid) intermetallic layer. Hence, the difference in reaction rate can be explained by a significant reduction in the diffusion flux due to a thicker compact layer. Finally, the mechanism of the growth of a thicker compact layer in the higher iron melt is proposed, based on the phase relations and diffusion both within and near the interfacial reaction zone. (C) 2004 Kluwer Academic Publishers.

A new equation for macroscopic description of capillary rise in porous media

Capillary rise in porous media is frequently modeled using the Washburn equation. Recent accurate measurements of advancing fronts clearly illustrate its failure to describe the phenomenon in the long term. The observed underprediction of the position of the front is due to the neglect of dynamic saturation gradients implicit in the formulation of the Washburn equation. We consider an approximate solution of the governing macroscopic equation, which retains these gradients, and derive new analytical formulae for the position of the advancing front, its speed of propagation, and the cumulative uptake. The new solution properly describes the capillary rise in the long term, while the Washburn equation may be recovered as a special case. (C) 2004 Elsevier Inc. All rights reserved.

Stress corrosion cracking and hydrogen diffusion in magnesium

Evaluation of recent data for hydrogen (H) diffusion in magnesium (Mg) yielded a new equation for the diffusion coefficient of H in Mg. This indicates that there can be significant H transport ahead of a stress corrosion crack in Mg at ambient temperature and that H may be involved in the mechanism of stress corrosion cracking in Mg.

Molecular Sieve Silica (MSS) Membranes for Gas Separation and Reaction Processes

Weakly branched silica films formed by the two-step sol-gel process allow for the formation of high selectivity membranes for gas separation. 29Si NMR and gas permeation showed that reduced crosslinking leads to He/CH4 selectivity improvement from 300 to 1000. Applied in membrane reactor for cyclohexane conversion to benzene, conversions were achieved at 14 fold higher than a conventional reactor at 250°C. Hydrothermal stability studies showed that carbon templating of silica is required for hydrothermally stable membranes. From our work it was shown that with correct application of chemistry, practical membrane systems can be built to suit gas separation (e. g. hydrogen fuel) and reactor systems.

On a semilinear Schrodinger equation with critical Sobolev exponent

We consider the semilinear Schrodinger equation -Deltau+V(x)u= K(x) \u \ (2*-2 u) + g(x; u), u is an element of W-1,W-2 (R-N), where N greater than or equal to4, V, K, g are periodic in x(j) for 1 less than or equal toj less than or equal toN, K>0, g is of subcritical growth and 0 is in a gap of the spectrum of -Delta +V. We show that under suitable hypotheses this equation has a solution u not equal 0. In particular, such a solution exists if K equivalent to 1 and g equivalent to 0.

Useful representations of series solutions for transport problems

Simple techniques are presented for rearrangement of an infinite series in a systematic way such that the convergence of the resulting expression is accelerated. These procedures also allow calculation of required boundary derivatives. Several examples of conduction and diffusion-reaction problems illustrate the methods.