Kinetics of the chemical looping oxidation of H2 by a co-precipitated mixture of CuO and Al2O3

Chemical-looping combustion (CLC) has the inherent property of separating CO2 from flue gases. Instead of air, it uses an oxygen-carrier, usually in the form of a metal oxide, to provide oxygen for combustion. When used for the combustion of gaseous fuels, such as natural gas, or synthesis gas from the gasification of coal, the technique gives a stream of CO2 which, on an industrial scale, would be sufficiently pure for geological sequestration. An important issue is the form of the metal oxide, since it must retain its reactivity through many cycles of complete reduction and oxidation. Here, we report on the rates of oxidation of one constituent of synthesis gas, H2, by co-precipitated mixtures of CuO+Al2O3 using a laboratory-scale fluidised bed. To minimise the influence of external mass transfer, and also of errors in the measurement of [H2], particles sized to 355-500μm were used at low [H2], with the temperature ranging from 450 to 900°C. Under such conditions, the reaction was slow enough for meaningful measurements of the intrinsic kinetics to be made. The reaction was found to be first order with respect to H2. Above ∼800°C, the reaction of CuO was fast and conformed to the shrinking core mechanism, proceeding via the intermediate, Cu2O, in: 2CuO+H2→Cu2O+H2O, ΔH1073 K0=- 116.8 kJ/mol; Cu2O+H2→2Cu+H2O, ΔH1073 K0-80.9 kJ/mol. After oxidation of the products Cu and Cu2O back to CuO, the kinetics in subsequent cycles of chemical looping oxidation of H2 could be approximated by those in the first. Interestingly, the carrier was found to react at temperatures as low as 300°C. The influence of the number of cycles of reduction and oxidation is explored. Comparisons are drawn with previous work using reduction by CO. Finally, these results indicate that the kinetics of reaction of the oxygen carrier with gasifier synthesis gases is very much faster than rates of gasification of the original fuel. © 2010 The Institution of Chemical Engineers.

Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.

Vibration analysis of laminated shells using a refined shear deformation theory

A previously published refined shear deformation theory is used to analyse free vibration of laminated shells. The theory includes the assumption that the transverse shear strains across any two layers are linearly dependent on each other. The theory has the same dependent variables as first-order shear deformation theory, hut the set of governing differential equations is of twelfth order. No shear correction factors are required. Free vibration of symmetric cross-ply laminated cylindrical shells, symmetric and antisymmetric cross-ply cylindrical panels is calculated. The numerical results are in good agreement with those from three-dimensional elasticity theory.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Vibration analysis of laminated plates using a refined shear deformation-theory

A previously published discrete-layer shear deformation theory is used to analyze free vibration of laminated plates. The theory includes the assumption that the transverse shear strains across any two layers are linearly dependent on each other. The theory has the same dependent variables as first order shear deformation theory, but the set of governing differential equations is of twelfth order. No shear correction factors are required. Free vibration of simply supported symmetric and antisymmetric cross-ply plates is calculated. The numerical results are in good agreement with those from three-dimensional elasticity theory.

A finite-element analysis of viscoelastically damped sandwich plates

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Two-time scales inner solutions and motion of a geostrophic vortex

Based on the authors' previous work, in this paper the systematical analyses on the motion and the inner solutions of a geostrophic vortex have been presented by means of thematched asymptotic expansion method with multiple time scales (S/gh001/2 and α S/gh001/2) and space scales. It has been shown that the leading inner solutions to the core structure in two-time scales analyses are identified with the results in normal one-time scale analyses. The time averages of the first-order solutions on short time variable τ are the same as the first-order solutions obtained in one normal time scale analyses. The geostrophic vortex induces an oscillatory motion in addition to moving with the background flow. The period, amplitude andthe deviation from the mean trajectory depend on the core structure and the initial conditions. The velocity of the motion of vortex center varies periodically and the time average of the velocity on short time variable τ is equal to the value of the local mean velocity.

Static equilibrium of horizontal ring prominence

In 1980 the Beijing Observatory had successively observed sevesal rare completely closed ring prominences whose ring plane was approximately parallel to the solar surface with a characteristic life about 1—2 days. In this paper we discuss the static equilibrium of this kind of horizontal ring plasma under the simultaneous actions of magnetic force, gravity and pressure gradients. Assuming ring plasma with axisymmetry and rectangular plasma cross-section and adopting closed magnetic field boundary condition from the basic equations we obtain the exact zero order general solutions for magnetic field (force-free field) and density (pressure). We further obtain an eigen-solution for the zero order magnetic field and density as well as the first order magnetic field, thus giving a kind of the possible distribution of magnetic field and density for the horizontal closed ring prominence. The closed magnetic structure of ring prominence as presented in this paper, has no link with the force lines of the outside corona magnetic field. This is helpful to explain the great temperature difference between prominenee and corona.

An application of adjoint variational method to the capillary instability in liquid

In this paper, applying the direct variational approach of first-order approximation to the capillary instability problem for the eases of rotating liquid column, toroid and films on both sides of cylinder, we have obtained the necessary and sufficient conditions for motion stability of the "cylindrical coreliquid-liquid-cylindrical shell" systems. The results obtained before are found to be special cases of the present investigation. At the same time, we have explained physical essence of rotating instability and settled a few disputes in previous investigations.

One-dimensional improvement and two-dimensional and 3-dimensional treatments for stable oscillation condition, mode structure and power output in high-speed-flow lasers

For high-speed-flow lasers, the one-dimensional and first-order approximate treatment in[1] under approximation of geometrical optics is improved still within the scope of approx-imation of geometrical optics. The strict accurate results are obtained, and what is more,two- and three-dimensional treatments are done. Thus for two- and three-dimensional cases, thestable oscillation condition, the formulae of power output and analytical expression of modesunder approximation of geometrical optics (in terms of gain function) are derived. Accord-ing to the present theory, one-and two-dimensional calculations for the typical case of Gerry'sexperiment are presented. All the results coincide well with the experiment and are better thanthe results obtained in [1].In addition, the applicable scope of Lee's stable oscillation condition given by [1] is ex-panded; the condition for the approximation of gcometrical optics to be applied to mode con-structure in optical cavity is obtained for the first time and the difference between thiscondition and that for free space is also pointed out in the present work.

The Empirical Modeling of an Ecosystem

The authors have endeavored to create a verified a-posteriori model of a planktonic ecosystem. Verification of an empirically derived set of first-order, quadratic differential equations proved elusive due to the sensitivity of the model system to changes in initial conditions. Efforts to verify a similarly derived set of linear differential equations were more encouraging, yielding reasonable behavior for half of the ten ecosystem compartments modeled. The well-behaved species models gave indications as to the rate-controlling processes in the ecosystem.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

DNS and LES of turbulent channel flow with hydrophobic surface

Hydrophobic surface benefits for drag reduction. Min and Kim[1] do the first Direct Numerical Simulation on drag reduction in turbulent channel flow. And Fukagata and Kasagi[2] make some theoretical analysis based on Dean[3]'s formula and some observations in the DNS results. Using their theory, they conclude that drag reduction is possible in large Reynolds number. Both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed in our research. How the LES behaving in the turbulent channel flow with hydrophobic surface is examined. Original Smagorinsky model and its Dynamical model are used in LES. The slip velocities predicted by LES using Dynamical model are in good agreement with DNS as shown in the Figure. Although the percentage of drag reduction predicted by LES shows some discrepancies, it is in the error limit for industrial flow. First order and second order moments of LES are also examined and compared with DNS's results. The first-order moments is calculated well by LES. But there are some discrepancies of second-order moments between LES and DNS. [GRAPHICS]

Magnetization reversal and exchange bias effects in hard/soft ferromagnetic bilayers with orthogonal anisotropies

21 p.