Computationsl modelling of dimethyl ether separation and steam reforming in fluidized bed reactors

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether (DME) gas adsorptive separation and steam reforming (DME-SR) in a large scale Circulating Fluidized Bed (CFB) reactor. The CFD model is based on Eulerian-Eulerian dispersed flow and solved using commercial software (ANSYS FLUENT). Hydrogen is currently receiving increasing interest as an alternative source of clean energy and has high potential applications, including the transportation sector and power generation. Computational fluid dynamic (CFD) modelling has attracted considerable recognition in the engineering sector consequently leading to using it as a tool for process design and optimisation in many industrial processes. In most cases, these processes are difficult or expensive to conduct in lab scale experiments. The CFD provides a cost effective methodology to gain detailed information up to the microscopic level. The main objectives in this project are to: (i) develop a predictive model using ANSYS FLUENT (CFD) commercial code to simulate the flow hydrodynamics, mass transfer, reactions and heat transfer in a large scale dual fluidized bed system for combined gas separation and steam reforming processes (ii) implement a suitable adsorption models in the CFD code, through a user defined function, to predict selective separation of a gas from a mixture (iii) develop a model for dimethyl ether steam reforming (DME-SR) to predict hydrogen production (iv) carry out detailed parametric analysis in order to establish ideal operating conditions for future industrial application. The project has originated from a real industrial case problem in collaboration with the industrial partner Dow Corning (UK) and jointly funded by the Engineering and Physical Research Council (UK) and Dow Corning. The research examined gas separation by adsorption in a bubbling bed, as part of a dual fluidized bed system. The adsorption process was simulated based on the kinetics derived from the experimental data produced as part of a separate PhD project completed under the same fund. The kinetic model was incorporated in FLUENT CFD tool as a pseudo-first order rate equation; some of the parameters for the pseudo-first order kinetics were obtained using MATLAB. The modelling of the DME adsorption in the designed bubbling bed was performed for the first time in this project and highlights the novelty in the investigations. The simulation results were analysed to provide understanding of the flow hydrodynamic, reactor design and optimum operating condition for efficient separation. Bubbling bed validation by estimation of bed expansion and the solid and gas distribution from simulation agreed well with trends seen in the literatures. Parametric analysis on the adsorption process demonstrated that increasing fluidizing velocity reduced adsorption of DME. This is as a result of reduction in the gas residence time which appears to have much effect compared to the solid residence time. The removal efficiency of DME from the bed was found to be more than 88%. Simulation of the DME-SR in FLUENT CFD was conducted using selected kinetics from literature and implemented in the model using an in-house developed user defined function. The validation of the kinetics was achieved by simulating a case to replicate an experimental study of a laboratory scale bubbling bed by Vicente et al [1]. Good agreement was achieved for the validation of the models, which was then applied in the DME-SR in the large scale riser section of the dual fluidized bed system. This is the first study to use the selected DME-SR kinetics in a circulating fluidized bed (CFB) system and for the geometry size proposed for the project. As a result, the simulation produced the first detailed data on the spatial variation and final gas product in such an industrial scale fluidized bed system. The simulation results provided insight in the flow hydrodynamic, reactor design and optimum operating condition. The solid and gas distribution in the CFB was observed to show good agreement with literatures. The parametric analysis showed that the increase in temperature and steam to DME molar ratio increased the production of hydrogen due to the increased DME conversions, whereas the increase in the space velocity has been found to have an adverse effect. Increasing temperature between 200 oC to 350 oC increased DME conversion from 47% to 99% while hydrogen yield increased substantially from 11% to 100%. The CO2 selectivity decreased from 100% to 91% due to the water gas shift reaction favouring CO at higher temperatures. The higher conversions observed as the temperature increased was reflected on the quantity of unreacted DME and methanol concentrations in the product gas, where both decreased to very low values of 0.27 mol% and 0.46 mol% respectively at 350 °C. Increasing the steam to DME molar ratio from 4 to 7.68 increased the DME conversion from 69% to 87%, while the hydrogen yield increased from 40% to 59%. The CO2 selectivity decreased from 100% to 97%. The decrease in the space velocity from 37104 ml/g/h to 15394 ml/g/h increased the DME conversion from 87% to 100% while increasing the hydrogen yield from 59% to 87%. The parametric analysis suggests an operating condition for maximum hydrogen yield is in the region of 300 oC temperatures and Steam/DME molar ratio of 5. The analysis of the industrial sponsor’s case for the given flow and composition of the gas to be treated suggests that 88% of DME can be adsorbed from the bubbling and consequently producing 224.4t/y of hydrogen in the riser section of the dual fluidized bed system. The process also produces 1458.4t/y of CO2 and 127.9t/y of CO as part of the product gas. The developed models and parametric analysis carried out in this study provided essential guideline for future design of DME-SR at industrial level and in particular this work has been of tremendous importance for the industrial collaborator in order to draw conclusions and plan for future potential implementation of the process at an industrial scale.

Pseudoperiodic "Living" and/or controlled cationic ring-opening copolymerization of oxetane with tetrahydropyran:microstructure of polymers vs kinetics of chain growth

The "living" and/or controlled cationic ring-opening bulk copolymerization of oxetane (Ox) with tetrahydropyran (THP) (cyclic ether with no homopolymerizability) at 35°C was examined using ethoxymethyl-1 -oxoniacyclohexane hexafluoroantimonate (EMOA) and (BF3 · CH3OH)THP as fast and slow initiator, respectively, yielding living and nonliving polymers with pseudoperiodic sequences (i.e., each pentamethylene oxide fragment inserted into the polymer is flanked by two trimethylene oxide fragments). Good control over number-average molecular weight (Mn up to 150000 g mol-1) with molecular weight distribution (MWD ∼ 1.4-1, 5) broader than predicted by the Poison distribution (MWDs > 1 +1/DPn) was attained using EMOA as initiating system, i.e., C 2H5OCH2Cl with 1.1 equiv of AgSbF6 as a stable catalyst and 1.1 equiv of 2,6-di-tert-butylpyridine used as a non-nucleophilic proton trap. With (BF3 · CH 3OH)THP, a drift of the linear dependence M n(GPC) vs Mn(theory) to lower molecular weight was observed together with the production of cyclic oligomers, ∼3-5% of the Ox consumed in THP against ∼30% in dichloromethane. Structural and kinetics studies highlighted a mechanism of chains growth where the rate of mutual conversion between "strain ACE species" (chain terminated by a tertiary 1-oxoniacyclobutane ion, Al) and "strain-free ACE species" (chain terminated by a tertiary 1-oxoniacyclohexane ion, Tl) depends on the rate at which Ox converts the stable species T1 (kind of "dormant" species) into a living "propagating" center A1 (i.e., k aapp[Ox]). The role of the THP solvent associated with the suspension of irreversible and reversible transfer reactions to polymer, when the polymerization is initiated with EMOA, was predicted by our kinetic considerations. The activation -deactivation pseudoequilibrium coefficient (Qt) was then calculated in a pure theoretical basis. From the measured apparent rate constant of Ox (kOxapp) and THP (kTHPapp = ka(endo)app) consumption, Qt and reactivity ratio (kp/kd, k a(endo)/ka(exo), and ks/ka(endo) were calculated, which then allow the determination of the transition rate constant of elementary step reactions that governs the increase of Mu with conversion. © 2009 American Chemical Society.

The chemistry of British bituminous coals : an assessment of phase transfer catalysed reactions in the determination of coal structure

Three British bituminous coals, (Gedling, Cresswell, and Cortonwood Silkstone) were selected for study. Procedures were developed, using phase transfer catalysts (PTC's), to degrade the solvent insoluble fractions of the coals. PTC's are of interest because they have the potential to bring about selective high conversion reactions, under mild conditions, (often in the past, severe reaction conditions have had to be used to degrade the coals, this in turn resulted in the loss of much of the structural information). We have applied a variety of physical and chemical techniques to maximise the amount of structural information, these include, elemental analysis, 1H-NMR, 13C-CPMAS-NMR, GPC, GC-MS, FTIR spectroscopy, DRIFT spectroscopy, and gas adsorption measurements. The main conclusions from the work are listed below:- ( 1 ) PTC O-methylation; This reaction removes hydrogen bonds within the coal matrix by 'capping' the phenolic groups. It was found that the polymer-like matrix could be made more flexible, but not significantly more soluble, by O-methylation. I.E. the trapped or 'mobile' phase of the coals could be removed at a faster rate after this reaction had been carried out. ( 2 ) PTC Reductive and Acidic Ether Cleavage; The three coals were found to contain insignificant amounts of dialkyl and alkyl aryl ethers. The number of diaryl ethers could not be estimated, by reductive ether cleavage, (even though a high proportion of all three coals was solublised). The majority of the ethers present in the coals were inert to both cleavage methods, and are therefore assumed to be heterocyclic ethers. ( 3 ) Trif!uoroperacetic Acid Oxidation; This oxidant was used to study the aliphatic portions of the polymer-like macromolecular matrix of the coals. Normally this reagent will only solublise low rank coals, we however have developed a method whereby trifluoroperacetic acid can be used to degrade high rank bituminous coals. ( 4 ) PTC/Permanganate Oxidation; This reagent has been found to be much more selective than the traditional alkaline permanganate oxidation, with a lot more structural information being retained within the various fractions. This degradative method therefore has the potential of yielding new information about the molecular structure of coals.

Mechanisms of transfer film formation on data recording heads

The increasing demand for high capacity data storage requires decreasing the head-to-tape gap and reducing the track width. A problem very often encountered is the development of adhesive debris on the heads at low humidity and high temperatures that can lead to an increase of space between the head and media, and thus a decrease in the playback signal. The influence of stains on the playback signal of reading heads is studied using RAW (Read After Write) tests and their influence on the wear of the heads by using indentation technique. The playback signal has been found to vary and the errors to increase as stains form a patchy pattern and grow in size to form a continuous layer. The indentation technique shows that stains reduce the wear rate of the heads. In addition, the wear tends to be more pronounced at the leading edge of the head compared to the trailing one. Chemical analysis of the stains using ferrite samples in conjunction with MP (metal particulate) tapes shows that stains contain iron particles and polymeric binder transferred from the MP tape. The chemical anchors in the binder used to grip the iron particles now react with the ferrite surface to create strong chemical bonds. At high humidity, a thin layer of iron oxyhydroxide forms on the surface of the ferrite. This soft material increases the wear rate and so reduces the amount of stain present on the heads. The stability of the binder under high humidity and under high temperature as well as the chemical reactions that might occur on the ferrite poles of the heads influences the dynamic behaviour of stains. A model of stain formation taking into account the channels of binder degradation and evolution upon different environmental conditions is proposed.

Computational modelling of the impact of particle size to the heat transfer coefficient between biomass particles and a fluidised bed

The fluid–particle interaction and the impact of different heat transfer conditions on pyrolysis of biomass inside a 150 g/h fluidised bed reactor are modelled. Two different size biomass particles (350 µm and 550 µm in diameter) are injected into the fluidised bed. The different biomass particle sizes result in different heat transfer conditions. This is due to the fact that the 350 µm diameter particle is smaller than the sand particles of the reactor (440 µm), while the 550 µm one is larger. The bed-to-particle heat transfer for both cases is calculated according to the literature. Conductive heat transfer is assumed for the larger biomass particle (550 µm) inside the bed, while biomass–sand contacts for the smaller biomass particle (350 µm) were considered unimportant. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Biomass reaction kinetics is modelled according to the literature using a two-stage, semi-global model which takes into account secondary reactions. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of User Defined Function (UDF).

Characterisation of a cavity transfer mixer as a chemical reactor

This research project examined the feasibility of using a cavity transfer mixer (CTM) as a continuous reactor to perform reactions between either solid or liquid reagents and polymer melt; reactions which have previously been typically carried out in batch reactor systems. Equipment has been developed to allow uniform and reproducible introduction of reagents into the polymer melt. Reactions have also been performed using batch processing equipment to enable comparison with the performance of the CTM. It was concluded that: a) there are certain reactions which cannot be carried out in a CTM, but which can be performed in a batch system such as a mill or a sigma blade mixer. This was found to be the case for some neutralisation reactions where the product was quasi crosslinked. b) the reactions that can be carried out in a CTM are performed more efficiently in a CTM than on a batch process. For example, when monomers were to be grafted onto polymers, this was more safely and efficiently performed in the CTM than in a mill or a sigma blade mixer. Residence time distributions (RTDs) for three CTMs were studied in order to gain an insight into the effect of CTM geometry on RTD, polymer melt flow pattern and reactor performance. A mathematical model has been developed to predict the influence of process parameters on RTD and the results compared with experimentally observed trends. The comparison was good. A programme of research has been drawn up to form the basis of an industrially based sponsored development project of the CTM reactor. This work programme was successfully marketed to companies with commercial interest in modified rubber and plastics as an integral part of the research programme of this thesis and the sponsored research programme has paralleled the work reported here.