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 method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.

A method of fundamental solutions for transient heat conduction

In this paper we investigate an application of the method of fundamental solutions (MFS) to transient heat conduction. In almost all of the previously proposed MFS for time-dependent heat conduction the fictitious sources are located outside the time-interval of interest. In our case, however, these sources are instead placed outside the space domain of interest in the same manner as is done for stationary heat conduction. A denseness result for this method is discussed and the method is numerically tested showing that accurate numerical results can be obtained. Furthermore, a test example with boundary singularities shows that it is advisable to remove such singularities before applying the MFS.

Interval data without sign restrictions in DEA

Conventional DEA models assume deterministic, precise and non-negative data for input and output observations. However, real applications may be characterized by observations that are given in form of intervals and include negative numbers. For instance, the consumption of electricity in decentralized energy resources may be either negative or positive, depending on the heat consumption. Likewise, the heat losses in distribution networks may be within a certain range, depending on e.g. external temperature and real-time outtake. Complementing earlier work separately addressing the two problems; interval data and negative data; we propose a comprehensive evaluation process for measuring the relative efficiencies of a set of DMUs in DEA. In our general formulation, the intervals may contain upper or lower bounds with different signs. The proposed method determines upper and lower bounds for the technical efficiency through the limits of the intervals after decomposition. Based on the interval scores, DMUs are then classified into three classes, namely, the strictly efficient, weakly efficient and inefficient. An intuitive ranking approach is presented for the respective classes. The approach is demonstrated through an application to the evaluation of bank branches. © 2013.

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.

A method of fundamental solutions for theaxisymmetric backward heat equation

We propose and investigate an application of the method of fundamental solutions (MFS) to the radially symmetric and axisymmetric backward heat conduction problem (BHCP) in a solid or hollow cylinder. In the BHCP, the initial temperature is to be determined from the temperature measurements at a later time. This is an inverse and ill-posed problem, and we employ and generalize the MFS regularization approach [B.T. Johansson and D. Lesnic, A method of fundamental solutions for transient heat conduction, Eng. Anal. Boundary Elements 32 (2008), pp. 697–703] for the time-dependent heat equation to obtain a stable and accurate numerical approximation with small computational cost.

Gas-fired radiant heaters

This dissertation covers four areas of particular interest for the successful application of radiant heating in industrial environments. In it the author tackles the problem of how to predict the thermal comfort produced both by single heaters and also systems of heaters; proposes a method for modelling the mechanisms by which heaters interact with the buildings in which they are installed, in the static and dynamic cases; explores techniques for measuring the radiation produced by heaters; and presents experiments concerned with finding the temperatures and power balances prevailing during normal operation. It is contended that, whilst the generally accepted guides for sizing and operating space heating plant were a good first approximation, there are intrinsic subtleties arising from the fact that the primary mode of heat transfer in this instance is radiative. These nuances are concerned with how best to maximise the heat transfer from the heat source to the heated object; the placement of heaters within a system; and an assessment of the various techniques and strategies involved in controlling a radiant heating system. The conclusions reached are that: if sized and controlled correctly radiant heating systems offer considerable operational advantages over other types of space heating systems in certain applications, in terms of both economy and controllability. The efficacy of radiant heating systems is affected primarily by the control strategy implemented; secondarily, by the structure of the building into which it is installed; and only marginally by all other factors.

Modeling and optimization of cascaded erbium and holmium doped fluoride fiber lasers

Numerical modeling of cascade erbium-doped and holmium-doped fluoride fiber lasers is presented. Fiber lengths were optimized for cascade lasers that had fixed or free-running wavelengths using all known spectroscopic parameters. The performance of the cascade laser was tested against dopant concentration, energy transfer process, heat generation, output coupling, and pump schemes. The results suggest that the slope efficiencies and thresholds for both transitions increase with increasing Ho3+ or Er3+ concentration with the slope efficiency stabilizing after 1 mol% rare earth doping. The heat generation in the Ho3+-based system is lower compared to the Er 3+-based system at low dopant concentration as a result of the lower rates of multiphonon relaxation. Decreasing the output coupling for the upper (∼3 μm) transition decreases the threshold of the lower transition and the upper transition benefits from decreasing the output coupling for the lower transition for both cascade systems. The highest slope efficiency was achieved under counter-propagating pump conditions. Saturation of the output power occurs at comparatively higher pump power with dilute Er3+ doping compared with heavier doping. Overall, we show that the cascade Ho3+ -doped fluoride laser is the best candidate for high power output because of its higher slope efficiency and lower temperature excursion of the core and no saturation of the output. © 2013 IEEE.

Uniform design for the optimization of Al2O3 nanofilms produced by electrophoretic deposition

Surface modification by means of nanostructures is of interest to enhance boiling heat transfer in various applications including the organic Rankine cycle (ORC). With the goal of obtaining rough and dense aluminum oxide (Al2O3) nanofilms, the optimal combination of process parameters for electrophoretic deposition (EPD) based on the uniform design (UD) method is explored in this paper. The detailed procedures for the EPD process and UD method are presented. Four main influencing conditions controlling the EPD process were identified as nanofluid concentration, deposition time, applied voltage and suspension pH. A series of tests were carried out based on the UD experimental design. A regression model and statistical analysis were applied to the results. Sensitivity analyses of the effect of the four main parameters on the roughness and deposited mass of Al2O3 films were also carried out. The results showed that Al2O3 nanofilms were deposited compactly and uniformly on the substrate. Within the range of the experiments, the preferred combination of process parameters was determined to be nanofluid concentration of 2 wt.%, deposition time of 15 min, applied voltage of 23 V and suspension pH of 3, yielding roughness and deposited mass of 520.9 nm and 161.6 × 10− 4 g/cm2, respectively. A verification experiment was carried out at these conditions and gave values of roughness and deposited mass within 8% error of the expected ones as determined from the UD approach. It is concluded that uniform design is useful for the optimization of electrophoretic deposition requiring only 7 tests compared to 49 using the orthogonal design method.