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.

Effects of adsorbed phase on diffusion of subcritical hydrocarbons in activated carbon at low pressures

Diffusions of free and adsorbed molecules of subcritical hydrocarbons in activated carbon were investigated to study the influence of adsorbed molecules on both diffusion processes at low pressures. A collision reflection factor, defined as the fraction of molecules undergoing collision to the solid surface over reflection from the surface, is incorporated into Knudsen diffusivity and surface diffusivity in meso/macropores. Since the porous structure of activated carbon is bimodal in nature, the diffusion of adsorbed molecules is contributed by that of weakly adsorbed molecules on the meso/macropore surfaces and that of strongly adsorbed molecules in the small confinement of micropores. The mobility of adsorbed molecules on the meso/macropore surface is characterized by the surface diffusivity D-mu 2, while that in the micropore is characterized by D-mu 1. In our study with subcritical hydrocarbons, we have found that the former increases almost linearly with pressure, while the latter exhibits a sharp increase at a very low-pressure region and then decreases beyond a critical pressure. This critical pressure is identified as a pressure at which the micropores are saturated.

Preparation and anode property of Pt-CeO2 electrodes supported on carbon black for direct methanol fuel cell applications

A platinum (Pt) on pure ceria (CeO2) supported by carbon black (CB) anode was synthesized using a combined process of precipitation and coimpregnation methods. The electrochemical activity of methanol oxidation reaction on synthesized Pt-CeO2/CB anodes was investigated by cyclic voltammetry and chronoamperometry experimentation. To improve the anode property on Pt-CeO2/CB, the influence of particle morphology and particle size on anode properties was examined. The morphology and particle size of the pure CeO2 particles could be controlled by changing the preparation conditions. The anode properties (i.e., peak current density and onset potential for methanol oxidation) were improved by using nanosize CeO2 particles. This indicates that a larger surface area and higher activity on the surface of CeO2 improve the anode properties. The influence of particle morphology of CeO2 on anode properties was not very large. The onset potential for methanol oxidation reaction on Pt-CeO2/CB, which consisted of CeO2 with a high surface area, was shifted to a lower potential compared with that on the anodes, which consisted of CeO2 with a low surface area. The onset potential on Pt-CeO2/CB at 60 degrees C became similar to that on the commercially available Pt-Ru/carbon anode. We suggest that the rate-determining steps of the methanol oxidation reaction on Pt-CeO2/CB and commercially available Pt-Ru/carbon anodes are different, which accounts for the difference in performance. In the reaction mechanism on Pt-CeO2/CB, we conclude that the released oxygen species from the surface of CeO2 particles contribute to oxidation of adsorbed CO species on the Pt surface. This suggests that the anode performance of the Pt-CeO2/CB anode would lead to improvements in the operation of direct methanol fuel cells at 80 degrees C by the enhancement of diffusion of oxygen species created from the surface of nanosized CeO2 particles. Therefore, we conclude that fabrication of nanosized CeO2 with a high surface area is a key factor for development of a high-quality Pt-CeO2/CB anode in direct methanol fuel cells.

Thermal Cycling Stability of Silica Membranes for Gas Separation

Hydrogen is being seen as an alternative energy carrier to conventional hydrocarbons to reduce greenhouse gas emissions. High efficiency separation technologies to remove hydrogen from the greenhouse gas, carbon dioxide, are therefore in growing demand. Traditional thermodynamic separation systems utilise distillation, absorption and adsorption, but are limited in efficiency at compact scales. Molecular sieve silica (MSS) membranes can perform this separation as they have high permselectivity of hydrogen to carbon dioxide, but their stability under thermal cycling is not well reported. In this work we exposed a standard MSS membrane and a carbonised template MSS (CTMSS) membrane to thermal cycling from 100 to 450°C. The standard MSS and carbonised template CTMSS membranes both showed permselectivity of helium to nitrogen dropping from around 10 to 6 in the first set of cycles, remaining stable until the last test. The permselectivity drop was due to small micropore collapse, which occurred via structure movement during cycling. Simulating single stage membrane separation with a 50:50 molar feed of H2:CO2, H2 exiting the permeate stream would start at 79% and stabilise at 67%. Higher selectivity membranes showed less of a purity drop, indicating the margin at which to design a stable membrane separation unit for CO2 capture.

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.

Hydrothermal stable templated molecular sieve silica (TMSS) membranes for gas separation

In this work we compare the hydrothermal stability performance of a Templated Molecular Sieve Silica (TMSS) membrane against a standard, non-templated Molecular Sieve Silica (MSS) membrane. The tests were carried under dry and wet (steam) conditions for single gas (He, H2, CO and CO2) at 1-2 atm membrane pressure drop at 200oC. Single gas TMSS membrane H2, permeance and H2/CO permselectivity was found to be 2.05 x 10-8 mols.m-2.s-1.Pa-1 and 15, respectively. The MSS membrane showed similar selectivity, but increased overall flux. He permeance through membranes decayed at a rate of 4-5 x 10-10 mols.m-2.s-1.Pa-1 per day regardless of membrane ambience (dry or wet). Although H2/CO permselectivity of the TMSS membrane slightly improved from 15 to 18 after steam testing, the MSS membrane resulted in significant reduction from 16 to 8.3. In addition, membrane regeneration after more than 50 days resulted in the TMSS membrane reverting to its original permeation levels while no significant improvements were observed for the MSS membra ne. Results showed that the TMSS membrane had enhanced hydrothermal stability and regeneration ability.

Hydrothermal Stability of Modified Silica Membranes for Gas Separation

A new class of hybrid molecular sieve silica (MSS) membranes is developed and tested against standard and organic templated membranes. The hybrid membrane is synthesized by the standard sol-gel process, integrating a template (methyltriethoxysilane - MTES) and a C6 surfactant (triethylhexylammonium bromide) into the silica film matrix. After hydro treatment under a relative humidity of 96% for 50h, the hybrid membrane shows no changes in its gas separation capabilities or energy of mobility. The structural characteristics and integrity of the hybrid membrane are retained due to a high concentration of organophilic functional groups and alkoxides observed using 29 Si NMR. In contrast, the structural integrity of the membranes prepared with non-templated films deteriorated during the hydro treatment due to a large percentage of silanol groups (Si-OH) which react with water. The hybrid membranes underwent a decrease in the H2/CO2 selectivity of only 1% whereas for the non-templated membrane a 21% decrease was observed. The transport mechanism of the hybrid membranes is activated as permeation increased with temperature. The activation energy for the permeation of H2 is positive while negative for CO2. The H2 permeation obtained was 3x 10 -8 mol.m -2 .s -1 .Pa -1 and permselectivities for H2/CO2 and H2/N2 varied between 1-7 and 31-34, respectively.

Canonical Bose gas simulations with stochastic gauges

A technique to simulate the grand canonical ensembles of interacting Bose gases is presented. Results are generated for many temperatures by averaging over energy-weighted stochastic paths, each corresponding to a solution of coupled Gross-Pitaevskii equations with phase noise. The stochastic gauge method used relies on an off-diagonal coherent-state expansion, thus taking into account all quantum correlations. As an example, the second-order spatial correlation function and momentum distribution for an interacting 1D Bose gas are calculated.

Pair correlations in a finite-temperature 1D Bose gas

We calculate the two-particle local correlation for an interacting 1D Bose gas at finite temperature and classify various physical regimes. We present the exact numerical solution by using the Yang-Yang equations and Hellmann-Feynman theorem and develop analytical approaches. Our results draw prospects for identifying the regimes of coherent output of an atom laser, and of finite-temperature “fermionization” through the measurement of the rates of two-body inelastic processes, such as photoassociation.

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).

Satellite protection and drag reduction using a purging gas flow

We have used the DSMC method to determine contamination (impingement of atmospheric molecules) and the aerodynamic forces on a cold satellite when a protective “purge gas” is ejected from a sting protruding ahead of the satellite. Forward ejection of the purge gas provides the greatest protection for a given mass of purge gas and the aerodynamic drag can be significantly reduced, thus compensating for the backward reaction from the forward ejection. If the purge gas is ejected backward from the sting (towards the satellite) the ejection provides thrust and the net retarding force can be reduced to zero. Contamination can be reduced and the mass of purging gas is less than the mass of conventional rocket propellant required to maintain the orbit of an unprotected satellite.

Short-term spatial diffusion in sigma o+ o- optical molasses

We have measured the spatial diffusion of atoms in a three-dimensional sigma(+)-sigma(-) optical molasses over twenty milliseconds timescale, starting from the initial interaction of the atoms with the molasses. We find that the diffusion constants agree well with a linear model for these short time scales and also compare favourably to other studies of diffusion made over longer time scales. These measurements enable us to quantify the detection method known as freezing molasses. We discuss this method, for detecting and measuring the momentum distribution of cold atoms, which relies on the slow diffusion of atoms in optical molasses to produce a freeze-frame of the spatial distribution of the atoms. This method enables a longer interrogation interval, providing a greatly increased signal-to-noise ratio. (C) 1998 Elsevier Science B.V.

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.