Chemical looping combustion (CLC) is a promising fossil fuel combustion technology, which is able to separate CO2 from the flue gases without a large consumption of energy. In this thesis, the study was extended to look at the use of chemical looping materials within traditional fluidised bed combustion and investigation of the interaction between the fuel, the supplied air and the chemical looping agent. Three topics of chemical looping combustion are discussed, including 1) the Sherwood number in the fluidised bed; 2) properties of different oxygen carriers, Fe2O3 and CuO (with supporting materials), were tested in the fluidised bed reactor; 3) the simulation of a steady state and a dynamic model of a coal-fired CLC power plant using Fe2O3 as oxygen carriers. The Sherwood number, which represents the mass transfer rate, is important in the calculation of CLC process. With Sherwood number, the mass transfer rate kg around the acting particle can be calculated using correlation Sh=kg∙d/D, where d is the diameter of acting particle, and D is the diffusivity around the acting particle. Hayhurst and Parmar (Hayhurst and Parmar 2002) calculated the Sherwood number in the fluidised bed by using the CO/CO2 ratio, which was measured by the temperature difference between the carbon particle and the bulk phase (Hayhurst and Parmar 1998). However, the temperature of the particle could be overestimated, so the CO/CO2 ratio could be underestimated. In this thesis, a universal exhaust gas oxygen (UEGO) sensor was employed, which could measure the actual carbon consumption rate in the fluidised bed by oxidizing CO in the sample gas into CO2 and. Fe particles of the same size of the char particle is used to measure the O2 consumption rate, and thus eliminate uncertainty in the Sherwood number. The CO/CO2 ratio was calculated by using the carbon consumption rate and the O2 consumption rate. In contrast to Hayhurst and Parmar (Hayhurst and Parmar 2002) who assumed CO2 was the main product, for this char the actual ratio of CO/CO2 was almost zero. The measurement here is in agreement with Arthur. This more accurate determination of CO/CO2 allows a better estimate of the mass transfer coefficient and leads to a correction of the Hayhurst and Parmar’s (Hayhurst and Parmar 2002) correlation by a factor of 1⁄2. Interestingly, very small fluidised beds have mass transfer coefficients which are about twice that expected in a large bed (owing to the very different flow and indeterminate flow pattern). This means the correlation of Hayhurst and Parmar (Hayhurst and Parmar 2002), by fortuitous coincidence works wells for beds with diameters < 30 mm., without the correction factor, should be ignored. In the fluidised bed in a typical CLC process, different fluidising material could have different influence on the reactions. Thus, it is worth discussing different kinds of fluidising materials. The char combustion in the fluidised bed was simulated by using inert (sand) and active (Fe2O3 or CuO) fluidising materials, and air as fluidising gas. The results indicated that 1) CO combustion in the boundary layer leads to smaller carbon consumption rate and larger oxygen consumption rate; 2) Using Fe2O3 particles as fluidising materials slows down the carbon consumption rate, since the diffusivity of CO2 is smaller than CO; 3) CuO particles slow down the carbon consumption rate at large Sherwood number (Sh=2 or 2.5). The influence of using CuO as fluidising material is further discussed experimentally by using low O2 fluidising gas. The results indicated that since the amount of CuO used in the experiment is small, when the O2 concentration in the bulk phase is lower than the equilibrium concentration, the O2 concentration in the bulk phase gradually decreases, and the O2 concentration in the bulk phase has large influence on the char particle combustion. A steady state model of a coal-fired CLC power plant was simulated. The aim of the model was to test the suitable operating conditions of the power plant, such as recycle rate of oxygen carriers, for the power plant design. In the steady state model, the power plant consists of a combustor and a steam cycle. Hambach lignite coal, Polish bituminous coal and natural gas were tested as fuels. The results indicated that: (1) The effect of the fuel is largely due to the amount of oxygen required per GJ released; (2) Preheating is important, but seems to have a minor effect since the most of the heat is released at temperatures well above the pinch point; (3) since the temperatures of heat source in this research is well above the pinch point, all heat are usable for the steam cycle. In this case, the steam cycle and the chemical looping plant could be optimised separately; (4) As long as the preheat temperature of the air flow into the air reactor is higher than the temperature of turbines, in most of cases the power output is unaffected by the choice of variables, leaving the designer free to choose the most convenient. With the conclusions above, a dynamic model of a coal-fired CLC power plant using Fe2O3 as oxygen carrier is then simulated. The aims of this simulation include: 1) explaining the kinetics of Fe2O3 oxygen carriers at high temperature (1223K) in a fluidised bed reactor using Brown’s data (Brown 2010); 2) a 1GWth dynamic power plant was simulated to test different cases including changing power supply and power storage. In the dynamic model, a chemical looping power plant using Hambach lignite char is tested, and the parameters of the system are adjusted so as to simulate the operations of a real chemical looping power plant. The two-phase model is employed for the fluidised bed reactors. Experimental data from Brown (Brown 2010) was simulated using this model first to test its validity. Then the model is scaled up to simulate a 1GWth dynamic power plant. The ideal operation conditions are found, and a char stripper is found helpful for carbon capture.

We studied the effect of a first order irreversible chemical reaction on mass transfer for two-phase flow systems in which the continuous phase is a fluid and the dispersed phase consists in catalystspherical particles. The reactive solute is transported by the fluid flow and penetrates through the particle surface by diffusion. The chemical reaction takes place within the bulk of the particle. Wehandle the problem by coupling mass balance equations for internal-external transfer with two boundary conditions: continuity of concentration and mass flux at the particle surface. We start with the case of a single isolated sphere. We propose a model to predict mass transfer coefficient (`reactive' Sherwood number) accounting for the external convection-diffusion along with internal diffusion-reaction. We validate the model through comparison with fully resolved Direct Numerical Simulations (DNS) performed by means of a boundary-fitted mesh method. For the simulation of multi-particle systems, we implemented a Sharp Interface Method to handle strong concentration gradients. We validate the implementation of the method thoroughly thanks to comparison with existing analytical solutions in case of diffusion, diffusion-reaction and by comparison with previously established correlations for convection-diffusion mass transfer. In case of convectiondiffusion- reaction, we validate the method and we evaluate its accuracy through comparisons with single particle simulations based on the boundary-fitted method. Later, we study the problem of three aligned-interacting spheres with internal chemical reaction. We propose a `reactive' Sherwood number model based on a known non-reactive prediction of mass transfer for each sphere. We validate the model by comparison with direct numerical simulations for a wide range of dimensionless parameters. Then, we study the configuration of a fixed bed of catalyst particles. We model the cup-mixing concentration profile, accounting for chemical reaction within the bed, and the mean surface and volume concentration profiles of the particles. We introduce a model for `reactive' Sherwood number that accounts for the solid volume fraction, in addition to the aforementioned effects. We compare the model to numerical simulations to evaluate its limitations

This thesis numerically studied mass transfer during fully confined liquid jet impingement on a rotating target disk of finite thickness and radius. The study involved laminar flow with jet Reynolds numbers from 650 to 1500. The nozzle to plate distance ratio was in the range of 0.5 to 2.0, the Schmidt number ranged from 1720 to 2513, and rotational speed was up to 325 rpm. In addition, the jet impingement to a stationary disk was also simulated for the purpose of comparison. The electrochemical fluid used was an electrolyte containing 0.005moles per liter potassium ferricyanide (K3(Fe(CN6)), 0.02moles per liter ferrocyanide (FeCN6?4), and 0.5moles per liter potassium carbonate (K2CO3). The rate of mass transfer of this electrolyte was compared to Sodium Hydroxide (NaOH) and Hydrochloric acid (HCl) electrochemical solutions. The material of the rotating disk was made of 99.98% nickel and 0.02% of chromium, cobalt and aluminum. The rate of mass transfer was also examined for different geometrical shapes of conical, convex, and concave confinement plates over a spinning disk. The results obtained are found to be in agreement with previous experimental and numerical studies. The study of heat transfer involved a microchannel for a composite channel of trapezoidal cross-section fabricated by etching a silicon wafer and bonding it with a slab of gadolinium. Gadolinium is a magnetic material that exhibits high temperature rise during adiabatic magnetization around its transition temperature of 295K. Heat was generated in the substrate by the application of magnetic field. Water, ammonia, and FC-77 were studied as the possible working fluids. Thorough investigation for velocity and temperature distribution was performed by varying channel aspect ratio, Reynolds number, and the magnetic field. The thickness of gadolinium slab, spacing between channels in the heat exchanger, and fluid flow rate were varied. To check the validity of simulation, the results were compared with existing results for single material channels. Results showed that Nusselt number is larger near the inlet and decreases downstream. Also, an increase in Reynolds number increases the total Nusselt number of the system.

Thesis submitted in fulfilment of the requirements for the degree Doctor technologiae (engineering: chemical) In the faculty of engineering at the cape peninsula university of technology / Since the first uses of hollow-fiber membrane bioreactors (MBR’s) to immobilize whole cells were reported in the early 1970’s, this technology has been used in as wide ranging applications as enzyme production to bone tissue engineering. The potential of these devices in industrial applications is often diminished by the large diffusional resistances of the membranes. Currently, there are no analytical studies on the performance of the MBR which account for both convective and diffusive transport. The purpose of this study was to quantify the efficiency of a biocatalytic membrane reactor used for the production of enzymes. This was done by developing exact solutions of the concentration and velocity profiles in the different regions of the membrane bioreactor (MBR). The emphasis of this study was on the influence of radial convective flows, which have generally been neglected in previous analytical studies. The efficiency of the MBR was measured by means of the effectiveness factor. An analytical model for substrate concentration profiles in the lumen of the MBR was developed. The model was based on the solution of the Navier-Stokes equations and Darcy’s law for velocity profiles, and the convective-diffusion equation for the solute concentration profiles. The model allowed for the evaluation of the influence of both hydrodynamic and mass transfer operating parameters on the performance of the MBR. These parameters include the fraction retentate, the transmembrane pressure, the membrane hydraulic permeability, the Reynolds number, the axial and radial Peclet numbers, and the dimensions of the MBR. The significant findings on the hydrodynamic studies were on the influence of the fraction retentate. In the dead-end mode it was found that there was increased radial convective flow, and hence more solute contact with the enzymes/biofilm immobilised on the surface of the membrane. The improved solute-biofilm contact however was only limited to the entrance half of the MBR. In the closed shell mode there was uniform distribution of solute, however, radial convective flows were significantly reduced. The developed model therefore allowed for the evaluation of an optimum fraction retentate value, where both the distribution of solutes and radial convective flows could be maximised.

Nous avons étudié l'effet d'une réaction chimique sur le transfert de matière pour des systèmes à deux phases sous écoulement. La phase continue est une phase fluide et la phase dispersée est constituée de particules de catalyseur au sein desquelles une réaction chimique irréversible de premier ordre a lieu. Le soluté réactif est transporté par l'écoulement externe de fluide et pénètre dans la particule par diffusion, il se produit alors une réaction chimique qui consomme cette espèce. Nous modélisons le problème par un couplage interne-externe des équations de bilan et au moyen de deux conditions limites de raccordement: continuité de la concentration et équilibre des flux de masse à la surface des particules. Le cas d'une seule sphère isolée est traitée en premier lieu de manière théorique et numérique. Nous proposons un modèle pour prédire le coefficient de transfert de masse (nombre de Sherwood «réactif») en tenant compte de la convection-diffusion externes et du couplage diffusion-réaction internes. Nous validons le modèle en le comparant à des simulations numériques directes pleinement résolues (DNS boundaryfitted) sur un maillage adapté à la géométrie des particules. Pour la simulation de systèmes multiparticules, nous mettons en œuvre une méthode d'interface «Sharp» pour traiter les fronts raides de concentration. Nous validons la mise en œuvre de la méthode sur des solutions analytiques existantes en cas de diffusion, de diffusion-réaction et par comparaison avec des corrélations de convection-diffusion disponibles dans la littérature. Dans le cas d'une réaction chimique en présence de convection-diffusion, nous validons la méthode et nous évaluons sa précision en comparant avec les simulations pleinement résolues de référence. Ensuite, nous étudions le problème de l'écoulement et du transfert autour de trois sphères alignées soumis à une réaction chimique interne. Nous proposons un modèle de nombre de Sherwood «réactif» en complément d'une prédiction de transfert pour chaque sphère disponible dans la littérature. Nous validons le modèle par comparaison avec des simulations numériques directes pour une large gamme de paramètres adimensionels. Ensuite, nous étudions la configuration du lit fixe de particules de catalyseur. Nous modélisons le profil de concentration moyenne, en tenant compte de la réaction chimique dans le lit et les profils de concentration moyenne surfacique et volumique des particules. Nous introduisons un modèle pour le nombre de Sherwood «réactif» qui est comparé à des simulations numériques pour en évaluer les limites de validité / We studied the effect of a first order irreversible chemical reaction on mass transfer for two-phase flow systems in which the continuous phase is a fluid and the dispersed phase consists in catalystspherical particles. The reactive solute is transported by the fluid flow and penetrates through the particle surface by diffusion. The chemical reaction takes place within the bulk of the particle. Wehandle the problem by coupling mass balance equations for internal-external transfer with two boundary conditions: continuity of concentration and mass flux at the particle surface. We start with the case of a single isolated sphere. We propose a model to predict mass transfer coefficient (`reactive' Sherwood number) accounting for the external convection-diffusion along with internal diffusion-reaction. We validate the model through comparison with fully resolved Direct Numerical Simulations (DNS) performed by means of a boundary-fitted mesh method. For the simulation of multi-particle systems, we implemented a Sharp Interface Method to handle strong concentration gradients. We validate the implementation of the method thoroughly thanks to comparison with existing analytical solutions in case of diffusion, diffusion-reaction and by comparison with previously established correlations for convection-diffusion mass transfer. In case of convectiondiffusion- reaction, we validate the method and we evaluate its accuracy through comparisons with single particle simulations based on the boundary-fitted method. Later, we study the problem of three aligned-interacting spheres with internal chemical reaction. We propose a `reactive' Sherwood number model based on a known non-reactive prediction of mass transfer for each sphere. We validate the model by comparison with direct numerical simulations for a wide range of dimensionless parameters. Then, we study the configuration of a fixed bed of catalyst particles. We model the cup-mixing concentration profile, accounting for chemical reaction within the bed, and the mean surface and volume concentration profiles of the particles. We introduce a model for `reactive' Sherwood number that accounts for the solid volume fraction, in addition to the aforementioned effects. We compare the model to numerical simulations to evaluate its limitations

Two computational fluid dynamics models are developed to predict the performance of a solid oxide fuel cell stack, a detailed and a simplified model. In the detailed model, the three dimensional momentum, heat, and species transport equations are coupled with electrochemistry. In the simplified model, the diffusion terms in the transport equations are selectively replaced by rate terms within the core region of the stack. This allows much coarser meshes to be employed at a fraction of the computational cost. Following the mathematical description of the problem, results for single-cell and multi-cell stacks are presented. Comparisons of local current density, temperature, and cell voltage indicate that good agreement is obtained between the detailed and simplified models, verifying the latter as a practical option in stack design. Then, the simplified model is used to determine the effects of utilization on the electrochemical performance and temperature distributions of a 10 cell stack. The results are presented in terms of fluid flow, pressure, species mass fraction, temperature, voltage and current density distributions. The effects of species and flow distributions on electrochemical performance and temperature are then analyzed for a 100 cell stack. The discussion highlights the importance of manifold design on performance and thermal management of large stacks. / Thesis (Master, Mechanical and Materials Engineering) -- Queen's University, 2013-09-30 15:55:18.627