Adsorption of dimethyl ether (DME) on zeolite molecular sieves

In recent years there has been growing interest in the use of dimethyl ether (DME) as an alternative fuel. In this study, the adsorption of DME on molecular sieves 4Å (Mol4A) and 5Å (Mol5A) has been experimentally investigated using the volumetric adsorption method. Data on the adsorption isotherms, heats of adsorption, and adsorption kinetic have been obtained and used to draw conclusions and compare the performance of the two adsorbents. Within the conditions considered, the adsorption capacity of Mol5A was found to be around eight times higher than the capacity of Mol4A. Low temperature adsorption and thermal pre-treatment of the adsorbents in vacuum were observed to be favourable for increased adsorption capacity. The adsorption isotherms for both adsorbent were fitted to the Freundlich model and the corresponding model parameters are proposed. The adsorption kinetic analysis suggest that the DME adsorption on Mol5A is controlled by intracrystalline diffusion resistance, while on Mol4A it is mainly controlled by surface layering resistance with the diffusion only taking place at the start of adsorption and for a very limited short time. The heats of adsorption were calculated by a calorimetric method based on direct temperature measurements inside the adsorption cell. Isosteric heats, calculated by the thermodynamic approach (Clasius-Clapeyron equation), have consistently shown lower values. The maximum heat of adsorption was found to be 25.9kJmol-1 and 20.1kJmol-1 on Mol4A and Mol5A, respectively; thus indicating a physisorption type of interactions. © 2014 Elsevier B.V.

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.

A computational model of hydrogen production by steam reforming of dimethyl ether in a large scale CFB reactor. Part II:parametric analysis

This study presents a computational parametric analysis of DME steam reforming in a large scale Circulating Fluidized Bed (CFB) reactor. The Computational Fluid Dynamic (CFD) model used, which is based on Eulerian-Eulerian dispersed flow, has been developed and validated in Part I of this study [1]. The effect of the reactor inlet configuration, gas residence time, inlet temperature and steam to DME ratio on the overall reactor performance and products have all been investigated. The results have shown that the use of double sided solid feeding system remarkable improvement in the flow uniformity, but with limited effect on the reactions and products. The temperature has been found to play a dominant role in increasing the DME conversion and the hydrogen yield. According to the parametric analysis, it is recommended to run the CFB reactor at around 300 °C inlet temperature, 5.5 steam to DME molar ratio, 4 s gas residence time and 37,104 ml gcat -1 h-1 space velocity. At these conditions, the DME conversion and hydrogen molar concentration in the product gas were both found to be around 80%.

Methyl chloride purification for the silicone industry

This study experimentally investigated methyl chloride (MeCl) purification method using an inhouse designed and built volumetric adsorption/desorption rig. MeCl is an essential raw material in the manufacture of silicone however all technical grades of MeCl contain concentrations (0.2 - 1.0 % wt) of dimethyl ether (DME) which poison the process. The project industrial partner had previously exhausted numerous separation methods, which all have been deemed not suitable for various reasons. Therefore, adsorption/desorption separation was proposed in this study as a potential solution with less economic and environmental impact. Pure component adsorption/desorption was carried out for DME and MeCl on six different adsorbents namely: zeolite molecular sieves (types 4 Å and 5 Å); silica gels (35-70 mesh, amorphous precipitated, and 35-60 mesh) and granular activated carbon (type 8-12 mesh). Subsequent binary gas mixture adsorption in batch and continuous mode was carried out on both zeolites and all three silica gels following thermal pre-treatment in vacuum. The adsorbents were tested as received and after being subjected to different thermal and vacuum pre-treatment conditions. The various adsorption studies were carried out at low pressure and temperature ranges of 0.5 - 3.5 atm and 20 - 100 °C. All adsorbents were characterised using Brunauer Emmett Teller (BET), thermogravimetric analysis (TGA), scanning electron microscopy (SEM) and energy dispersive x-ray analysis (EDXA) to investigate their physical and chemical properties. The well-known helium (He) expansion method was used to determine the empty manifold and adsorption cell (AC) regions and respective void volumes for the different adsorbents. The amounts adsorbed were determined using Ideal gas laws via the differential pressure method. The heat of adsorption for the various adsorbate-adsorbent (A-S) interactions was calculated using a new calorimetric method based on direct temperature measurements inside the AC. Further adsorption analysis included use of various empirical and kinetic models to determine and understand the behaviour of the respective interactions. The gas purification behaviour was investigated using gas chromatography and mass spectroscopy (GC-MC) analysis. Binary gas mixture samples were syringed from the manifold iii and AC outlet before and after adsorption/desorption analysis through manual sample injections into the GC-MS to detect and quantify the presence of DME and ultimately observe for methyl chloride purification. Convincing gas purification behaviour was confirmed using two different GC columns, thus giving more confidence on the measurement reliability. From the single pure component adsorption of DME and MeCl on the as received zeolite 4A subjected to 1 h vacuum pre-treatment, both gases exhibited pseudo second order adsorption kinetics with DME exhibiting a rate constant nearly double that of MeCl thus suggesting a faster rate of adsorption. From the adsorption isotherm classification both DME and MeCl exhibited Type II and I adsorption isotherm classifications, respectively. The strength of bonding was confirmed by the differential heat of adsorption measurement, which was found to be 23.30 and 10.21 kJ mol-1 for DME and MeCl, respectively. The former is believed to adsorb heterogeneously through hydrogen bonding whilst MeCl adsorbs homogenously via van der Waal’s (VDW) forces. Single pure component adsorption on as received zeolite 5A, silica gels (35-70, amorphous precipitated and 35-60) resulted in similar adsorption/desorption behaviour in similar quantities (mol kg-1). The adsorption isotherms for DME and MeCl on zeolite 5A, silica gels (35-70, amorphous precipitated and 35-60) and activated carbon 8-12 exhibited Type I classifications, respectively. Experiments on zeolite 5A indicated that DME adsorbed stronger, faster and with a slightly stronger strength of interaction than MeCl but in lesser quantities. On the silica gels adsorbents, DME exhibited a slightly greater adsorption capacity whilst adsorbing at a similar rate and strength of interaction compared to MeCl. On the activated carbon adsorbent, MeCl exhibited the greater adsorption capacity at a faster rate but with similar heats of adsorption. The effect of prolonged vacuum (15 h), thermal pre-treatment (150 °C) and extended equilibrium time (15 min) were investigated for the adsorption behaviour of DME and MeCl on both zeolites 4A and 5A, respectively. Compared to adsorption on as received adsorbents subjected to 1 h vacuum the adsorption capacities for DME and MeCl were found to increase by 1.95 % and 20.37 % on zeolite 4A and by 4.52 % and 6.69 % on zeolite 5A, respectively. In addition the empirical and kinetic models and differential heats of adsorption resulted in more definitive fitting curves and trends due to the true equilibrium position of the adsorbate with the adsorbent. Batch binary mixture adsorption on thermally and vacuum pre-treated zeolite 4A demonstrated purification behaviour of all adsorbents used for MeCl streams containing DME impurities, with a concentration as low as 0.66 vol. %. The GC-MS analysis showed no DME detection for the tested concentration mixtures at the AC outlet after 15 or 30 min, whereas MeCl was detectable in measurable amounts. Similar behaviour was also observed when carrying out adsorption in continuous mode. On the other hand, similar studies on the other adsorbents did not show such favourable MeCl purification behaviour. Overall this study investigated a wide range of adsorbents (zeolites, silica gels and activated carbon) and demonstrated for the first time potential to purify MeCl streams containing DME impurities using adsorption/desorption separation under different adsorbent pre-treatment and adsorption operating conditions. The study also revealed for the first time the adsorption isotherms, empirical and kinetic models and heats of adsorption for the respective adsorbentsurface (A-S) interactions. In conclusion, this study has shown strong evidence to propose zeolite 4A for adsorptive purification of MeCl. It is believed that with a technical grade MeCl stream competitive yet simultaneous co-adsorption of DME and MeCl occurs with evidence of molecular sieiving effects whereby the larger DME molecules are unable to penetrate through the adsorbent bed whereas the smaller MeCl molecules diffuse through resulting in a purified MeCl stream at the AC outlet. Ultimately, further studies are recommended for increased adsorption capacities by considering wider operating conditions, e.g. different adsorbent thermal and vacuum pre-treatment and adsorbing at temperatures closer to the boiling point of the gases and different conditions of pressure and temperature.

Hydrogen production by steam reforming of DME in a large scale CFB reactor. Part I:computational model and predictions

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether 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). The DME-SR reactions scheme and kinetics in the presence of a bifunctional catalyst of CuO/ZnO/Al2O3+ZSM-5 were incorporated in the model using in-house developed user-defined function. The model was validated by comparing the predictions with experimental data from the literature. The results revealed for the first time detailed CFB reactor hydrodynamics, gas residence time, temperature distribution and product gas composition at a selected operating condition of 300 °C and steam to DME mass ratio of 3 (molar ratio of 7.62). The spatial variation in the gas species concentrations suggests the existence of three distinct reaction zones but limited temperature variations. The DME conversion and hydrogen yield were found to be 87% and 59% respectively, resulting in a product gas consisting of 72 mol% hydrogen. In part II of this study, the model presented here will be used to optimize the reactor design and study the effect of operating conditions on the reactor performance and products.

A new strategy for selecting and evaluating physical solvents for gas absorption

This work describes how the physical properties of a solvent affect the design variables of a physical gas absorption process. The role of every property in determining the capital and the running cost of a process has been specified. Direct mathematical relationships have been formulated between every item of capital or running cost and the properties which are related to that item. The accuracy of the equations formulated has been checked by comparing their outcome with some actual design data. A good agreement has been found. The equations formulated may be used to evaluate on the basis of economics any suggested new solvents. A group of solvents were selected for evaluation. Their physical properties were estimated or collected as experimental data. The selected ones include three important solvents, the first is polyethylene glycol dimethyl ether (Selexol) which represents the currently most successful one, The other two solvents are acetonyl acetone (B2) and n-formyl morpholine which have been suggested previously as potential credible alternatives to the current ones. The important characteristics of: acetonyl acetone are its high solubility and its low viscosity, while the n-formyl morpholine is characterised by its low vapour pressure and its high selectivity. It was found that acetonyl acetone (B2) is the most attractive solvent for commercial applications particularly for process configurations that:include heat exchangers and strippers. The effect of the process configuration on the selected solvent was investigated in detail and it was found that there is no universal solvent which is the best for any process configuration, but that there is a best solvent for a given process configuration. In previous work, acetonyl acetone was suggested as a commercially promising physical solvent. That suggestion was not fully based on experimental measurement of all the physical properties. The viscosity of acetonyl acetone and its solubility at 1 atm were measured but the vapour pressure and the solubility of C02 and CH4 at high pressure were predicted. In this work, the solubilities of C02, CH4 and C3H8 in acetenyl acetone were measured for a partial pressure range of (2 ~ 22) bar at 25°C, The vapour pressure of this solvent was also measured, and the Antoine equation was formulated from tbe experimental data. The experimental data were found to be not In agreement with the predicted ones, so acetonyl acetone was re-evaluated according to the experimental data. It was found that this solvent can be recommended for further trials in a pilot plant study or for small scale commercial units.

The molecular design of a new solvent for the absorption of carbon dioxide

Gas absorption, the removal of one or more constitutents from a gas mixture, is widely used in chemical processes. In many gas absorption processes, the gas mixture is already at high pressure and in recent years organic solvents have been developed for the process of physical absorption at high pressure followed by low pressure regeneration of the solvent and recovery of the absorbed gases. Until now the discovery of new solvents has usually been by expensive and time consuming trial and error laboratory tests. This work describes a new approach, whereby a solvent is selected from considerations of its molecular structure by applying recently published methods of predicting gas solubility from the molecular groups which make up the solvent molecule. The removal of the acid gases of carbon dioxide and hydrogen sulfide from methane or hydrogen was used as a commercially important example. After a preliminary assessment to identify promising moecular groups, more than eighty new solvent molecules were designed and evaluated by predicting gas solubility. The other important physical properties were also predicted by appropriate theoretical procedures, and a commercially promising new solvent was chosen to have a high solubility for acid gases, a low solubility for methane and hydrogen, a low vapour pressure, and a low viscosity. The solvent chosen, of molecular structure Ch3-COCH2-CH2-CO-CH3, was tested in the laboratory and shown to have physical properties, except for vapour pressures, close to those predicted. That is gas solubilities were within 10% but lower than predicted. Viscosity within 10% but higher than predicted and a vapour pressure significantly lower than predicted. A computer program was written to predict gas solubility in the new solvent at the high pressures (25 bar) used in practice. This is based on the group contribution method of Skold Jorgensen (1984). Before using this with the new solvent, Acetonyl acetone, the method was show to be sufficiently accurate by comparing predicted values of gas solubility with experimental solubilities from the literature for 14 systems up to 50 bar. A test of the commercial potential of the new solvent was made by means of two design studies which compared the size of plant and approximate relative costs of absorbing acid gases by means of the new solvent with other commonly used solvents. These were refrigerated methanol(Rectisol process) and Dimethyl Ether or Polyethylene Glycol(Selexol process). Both studies showed in terms of capital and operating cost some significant advantage for plant designed for the new solvent process.

Bifunctional organorhodium solid acid catalysts for methanol carbonylation

Robust, bifunctional catalysts comprising Rh(CO)(Xantphos) exchanged phosphotungstic acids of general formulas [Rh(CO)(Xantphos)]+n[H3–nPW12O40]n− have been synthesized over silica supports which exhibit tunable activity and selectivity toward direct vapor phase methanol carbonylation. The optimal Rh:acid ratio = 0.5, with higher rhodium concentrations increasing the selectivity to methyl acetate over dimethyl ether at the expense of lower acidity and poor activity. On-stream deactivation above 200 °C reflects Rh decomplexation and reduction to Rh metal, in conjunction with catalyst dehydration and loss of solid acidity because of undesired methyl acetate hydrolysis, but can be alleviated by water addition and lower temperature operation.

Studies of the cationic polymerisation of vinylic and energetic cyclic ether monomers

The cationic polymerisation of various monomers, including cyclic ethers bearing energetic nitrate ester (-ON02) groups, substituted styrenes and isobutylene has been investigated. The main reaction studied has been the ring-opening polymerisation of 3- (nitratomethyl)-3-methyl oxetane (NIMMO) using the alcohol/BF3.0Et2 binary initiator system. A series of di-, tri- and tetrafunctional telechelic polymers has been synthesised. In order to optimise the system, achieve controlled molecular weight polymers and understand the mechanism of polymerisation the effects of certain parameters on the molecular weight distribution, as determined by Size Exclusion Chromatography, have been examined. This shows the molecular weight achieved depends on a combination of factors including -OH concentration, addition rate of monomer and, most importantly, temperature. A lower temperature and OH concentration tends to produce higher molecular weight, whereas, slower addition rates of monomer, either have no significant effect or produce a lower molecular weight polymer. These factors were used to increase the formation of a cyclic oligomer, by a side reaction, and suggest, that the polymerisation of NIMMO is complicated with endbiting and back biting reactions, along with other transfer/termination processes. These observations appear to fit the model of an active-chain end mechanism. Another cyclic monomer, glycidyl nitrate (GLYN), has been polymerised by the activated monomer mechanism. Various other monomers have been used to end-cap the polymer chains to produce hydroxy ends which are expected to form more stable urethane links, than the glycidyl nitrate ends, when cured with isocyanates. A novel monomer, butadiene oxide dinitrate (BODN), has been prepared and its homopolymerisation and copolymerisation with GL YN studied. In concurrent work the carbocationic polymerisations of isobutylene or substituted styrenes have been studied. Materials with narrow molecular weight distributions have been prepared using the diphenyl phosphate/BCl3 initiator. These systems and monomers are expected to be used in the synthesis of thermoplastic elastomers.

Cationic ring opening polymerisation (CROP) of oxetane and its derivatives using oxonium derived initiators

Several cationic initiator systems were developed and used to polymerise oxetane with two oxonium ion initiator systems being investigated in depth. The first initiator system was generated by the elimination of a chloride group from a chloro methyl ethyl ether. Adding a carbonyl co-catalyst to a carbocationic centre generated the second initiator system. It was found that the anion used to stabilise the initiator was critical to the initial rate of polymerisation of oxetane with hexafluoroantimonate resulting in the fastest polymerisations. Both initiator systems could be used at varying monomer to initiator concentrations to control the molecular number average, Mn, of the resultant polymer. Both initiator systems showed living characteristics and were used to polymerise further monomers and generate higher molecular weight material and block copolymers. Oxetane and 3,3-dimethyl oxetane can both be polymerised using either oxonium ion initiator system in a variety of DCM or DCM/1,4-dioxane solvent mixtures. The level of 1,4-dioxane does have an impact on the initial rate of polymerisation with higher levels resulting in lower initial rates of polymerisation but do tend to result in higher polydispersities. The level of oligomer formation is also reduced as the level of 1,4-dioxane is increased. 3,3-bis-bromomethyl oxetane was also polymerised but a large amount of hyperbranching was seen at the bromide site resulting in a difficult to solvate polymer system. Multifunctional initiator systems were also generated using the halide elimination reactions with some success being achieved with 1,3,5-tris-bromomethyl-2,4,6-tris-methyl-benzene derived initiator system. This offered some control over the molecular number average of the resultant polymer system.