Comparison of Biomass Gasification Technologies for Natural Gas Substitution in Iron Ore Pelletizing Plants

The iron ore pelletizing process consumes high amounts of energy, including nonrenewable sources, such as natural gas. Due to fossil fuels scarcity and increasing concerns regarding sustainability and global warming, at least partial substitution by renewable energy seems inevitable. Gasification projects are being successfully developed in Northern Europe, and large-scale circulating fluidized bed biomass gasifiers have been commissioned in e.g. Finland. As Brazil has abundant biomass resources, biomass gasification is a promising technology in the near future. Biomasses can be converted into product gas through gasification. This work compares different technologies, e.g. air, oxygen and steam gasification, focusing on the use of the product gas in the indurating machine. The use of biosynthetic natural gas is also evaluated. Main parameters utilized to assess the suitability of product gas were adiabatic flame temperature and volumetric flow rate. It was found that low energy content product gas could be utilized in the traveling grate, but it would require burner’s to be changed. On the other hand, bio-SGN could be utilized without any adaptions. Economical assessment showed that all gasification plants are feasible for sizes greater than 60 MW. Bio-SNG production is still more expensive than natural gas in any case.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

The Possibility of Greenhouse Gas Mitigation in Ethiopian Cement Industry

Cement industry significantly associated with high greenhouse gas (GHG) emissions. Considering the environmental impact, particularly global warming potential, it is important to reduce these emissions to air. The aim of the study is to investigate the mitigation possibility of GHG emissions in Ethiopian cement industry. Life cycle assessment (LCA) method used to identify and quantify GHG emissions during one ton of ordinary portland cement (OPC) production. Three mitigation scenarios: alternative fuel use, clinker substitution and thermal energy efficiency were applied on a representative gate-to-gate flow model developed with GaBi 6 software. The results of the study indicate that clinker substitution and alternative fuel use play a great role for GHG emissions mitigation with affordable cost. Applying most energy efficient kiln technology, which in turn reduces the amount of thermal energy use, has the least GHG emissions reduction intensity and high implementation cost comparing to the other scenarios. It was found that the cumulative GHG emissions mitigation potential along with other selected mitigation scenarios can be at least 48.9% per ton of cement production.

Integrated energy storage tank for large scale power-to-gas applications

The global interest towards renewable energy production such as wind and solar energy is increasing, which in turn calls for new energy storage concepts due to the larger share of intermittent energy production. Power-to-gas solutions can be utilized to convert surplus electricity to chemical energy which can be stored for extended periods of time. The energy storage concept explored in this thesis is an integrated energy storage tank connected to an oxy-fuel combustion plant. Using this approach, flue gases from the plant could be fed directly into the storage tank and later converted into synthetic natural gas by utilizing electrolysis-methanation route. This work utilizes computational fluid dynamics to model the desublimation of carbon dioxide inside a storage tank containing cryogenic liquid, such as liquefied natural gas. Numerical modelling enables the evaluation of the transient flow patterns caused by the desublimation, as well as general fluid behaviour inside the tank. Based on simulations the stability of the cryogenic storage and the magnitude of the key parameters can be evaluated.

Three-dimensional modeling of biomass fuel flow in a circulating fluidized bed furnace

The reduction of greenhouse gas emissions in the European Union promotes the combustion of biomass rather than fossil fuels in energy production. Circulating fluidized bed (CFB) combustion offers a simple, flexible and efficient way to utilize untreated biomass in a large scale. CFB furnaces are modeled in order to understand their operation better and to help in the design of new furnaces. Therefore, physically accurate models are needed to describe the heavily coupled multiphase flow, reactions and heat transfer inside the furnace. This thesis presents a new model for the fuel flow inside the CFB furnace, which acknowledges the physical properties of the fuel and the multiphase flow phenomena inside the furnace. This model is applied with special interest in the firing of untreated biomass. An experimental method is utilized to characterize gas-fuel drag force relations. This characteristic drag force approach is developed into a gas-fuel drag force model suitable for irregular, non-spherical biomass particles and applied together with the new fuel flow model in the modeling of a large-scale CFB furnace. The model results are physically valid and achieve very good correspondence with the measurement results from large-scale CFB furnace firing biomass. With the methods and models presented in this work, the fuel flow field inside a circulating fluidized bed furnace can be modeled with better accuracy and more efficiently than in previous studies with a three-dimensional holistic model frame.

Effects of Fin Parameters on Gas-Side Heat Transfer Performance of H-Type Finned Tube Bundle in a Waste Heat Recovery Boiler

Alfa Laval Aalborg Oy designs and manufactures waste heat recovery systems utilizing extended surfaces. The waste heat recovery boiler considered in this thesis is a water-tube boiler where exhaust gas is used as the convective heat transfer medium and water or steam flowing inside the tubes is subject to cross-flow. This thesis aims to contribute to the design of waste heat recovery boiler unit by developing a numerical model of the H-type finned tube bundle currently used by Alfa Laval Aalborg Oy to evaluate the gas-side heat transfer performance. The main objective is to identify weaknesses and potential areas of development in the current H-type finned tube design. In addition, numerical simulations for a total of 15 cases with varying geometric parameters are conducted to investigate the heat transfer and pressure drop performance dependent on H-type fin geometry. The investigated geometric parameters include fin width and height, fin spacing, and fin thickness. Comparison between single and double tube type configuration is also conducted. Based on the simulation results, the local heat transfer and flow behaviour of the H-type finned tube is presented including boundary layer development between the fins, the formation of recirculation zone behind the tubes, and the local variations of flow velocity and temperature within the tube bundle and on the fin surface. Moreover, an evaluation of the effects of various fin parameters on heat transfer and pressure drop performance of H-type finned tube bundle has been provided. It was concluded that from the studied parameters fin spacing and fin width had the most significant effect on tube bundle performance and the effect of fin thickness was the least important. Furthermore, the results suggested that the heat transfer performance would increase due to enhanced turbulence if the current double tube configuration is replaced with single tube configuration, but further investigation and experimental measurements are required in order to validate the results.

Flow injection analysis (FIA) determination of lead by hydride generation in banded coral skeletons /

A flow injection hydride generation direct current plasma atomic emission spectrometric (FI-HG-DCP-AES) method was developed for the determination of lead at ng.ml-l level. Potassium ferricyanide (K3Fe(CN)6) was used along with sodium tetrahydroborate(III) (NaBH4) to produce plumbane (PbH4) in an acid medium. The design of a gas-liquid separator (hydride generator) was tested and the parameters of the flow injection system were optimized to achieve a good detection limit and sample throughput. The technique developed gave a detection limit of 0.7 ng.ml-l(3ob). The precision at 20 ng.ml"* level was 1.6 % RSD with 1 1 measurements (n=l 1). Volume of sample loop was 500 |J.l. A sample throughput of 120 h"^ was achieved. The transition elements, Fe(II), FeOH), Cd(n), Co(II), Mn(n), Ni(II) and Zn(n) do not interfere in this method but 1 mg,l'l Cu(II) will suppress 50 % of the signal from a sample containing 20 ng.ml'l Pb. This method was successfully applied to determine lead in a calcium carbonate (CaC03) matrix of banded coral skeletons from Si-Chang Island in Thailand.

Investigations into the extraction and the determination of polycyclic aromatic hydrocarbons (PAHs) from water by solid-phase extraction and capillary gas chromatography/ mass spectrometry

Factors affecting the detennination of PAHs by capillary GC/MS were studied. The effect of the initial column temperature and the injection solvent on the peak areas and heights of sixteen PAHs, considered as priority pollutants, USillg crosslinked methyl silicone (DB!) and 5% diphenyl, 94% dimethyl, 1% vinyl polysiloxane (DBS) columns was examined. The possibility of using high boiling point alcohols especially butanol, pentanol, cyclopentanol, and hexanol as injection solvents was investigated. Studies were carried out to optimize the initial column temperature for each of the alcohols. It was found that the optimum initial column temperature is dependent on the solvent employed. The peak areas and heights of the PAHs are enhanced when the initial column temperature is 10-20 c above the boiling point of the solvent using DB5 column, and the same or 10 C above the boiling point of the solvent using DB1 column. Comparing the peak signals of the PAHs using the alcohols, p-xylene, n-octane, and nonane as injection solvents, hexanol gave the greatest peak areas and heights of the PAHs particularly the late-eluted peaks. The detection limits were at low pg levels, ranging from 6.0 pg for fluorene t9 83.6 pg for benzo(a)pyrene. The effect of the initial column temperature on the peak shape and the separation efficiency of the PARs was also studied using DB1 and DB5 columns. Fronting or splitting of the peaks was obseIVed at very low initial column temperature. When high initial column temperature was used, tailing of the peaks appeared. Great difference between DB! and.DB5 columns in the range of the initial column temperature in which symmetrical.peaks of PAHs can be obtained is observed. Wider ranges were shown using DB5 column. Resolution of the closely-eluted PAHs was also affected by the initial column temperature depending on the stationary phase employed. In the case of DB5, only the earlyeluted PAHs were affected; whereas, with DB1, all PAHs were affected. An analytical procedure utilizing solid phase extraction with bonded phase silica (C8) cartridges combined with GC/MS was developed to analyze PAHs in water as an alternative method to those based on the extraction with organic solvent. This simple procedure involved passing a 50 ml of spiked water sample through C8 bonded phase silica cartridges at 10 ml/min, dried by passing a gentle flow of nitrogen at 20 ml/min for 30 sec, and eluting the trapped PAHs with 500 Jll of p-xylene at 0.3 ml/min. The recoveries of PAHs were greater than 80%, with less than 10% relative standard deviations of nine determinations. No major contaminants were present that could interfere with the recognition of PAHs. It was also found that these bonded phase silica cartridges can be re-used for the extraction of PAHs from water.

Investigations in flow injection analysis (FIA): conventional FIA determination of chloride in water; sequential injection (SI)-FIA determination of mercury in water

Flow injection analysis (FIA) was applied to the determination of both chloride ion and mercury in water. Conventional FIA was employed for the chloride study. Investigations of the Fe3 +/Hg(SCN)2/CI-,450 nm spectrophotometric system for chloride determination led to the discovery of an absorbance in the 250-260 nm region when Hg(SCN)2 and CI- are combined in solution, in the absence of iron(III). Employing an in-house FIA system, absorbance observed at 254 nm exhibited a linear relation from essentially 0 - 2000 Jlg ml- 1 injected chloride. This linear range spanning three orders of magnitude is superior to the Fe3+/Hg(SCN)2/CI- system currently employed by laboratories worldwide. The detection limit obtainable with the proposed method was determin~d to be 0.16 Jlg ml- 1 and the relative standard deviation was determined to be 3.5 % over the concentration range of 0-200 Jig ml- 1. Other halogen ions were found to interfere with chloride determination at 254 nm whereas cations did not interfere. This system was successfully applied to the determination of chloride ion in laboratory water. Sequential injection (SI)-FIA was employed for mercury determination in water with the PSA Galahad mercury amalgamation, and Merlin mercury fluorescence detection systems. Initial mercury in air determinations involved injections of mercury saturated air directly into the Galahad whereas mercury in water determinations involved solution delivery via peristaltic pump to a gas/liquid separator, after reduction by stannous chloride. A series of changes were made to the internal hardware and valving systems of the Galahad mercury preconcentrator. Sequential injection solution delivery replaced the continuous peristaltic pump system and computer control was implemented to control and integrate all aspects of solution delivery, sample preconcentration and signal processing. Detection limits currently obtainable with this system are 0.1 ng ml-1 HgO.

Investigations into the determination of polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) by capillary gas chromatography

Factors involved in the determination of PAHs (16 priority PAHs as an example) and PCBs (10 PCB congeners, representing 10 isomeric groups) by capillary gas chromatography coupled with mass spectrometry (GC/MS, for PAHs) and electron capture detection (GC/ECD , for PCBs) were studied, with emphasis on the effect of solvent. Having various volatilities and different polarities, solvent studied included dichloromethane, acetonitrile, hexan e, cyclohexane, isooctane, octane, nonane, dodecane, benzene, toluene, p-xylene, o-xylene, and mesitylene. Temperatures of the capillary column, the injection port, the GC/MS interface, the flow rates of carrier gas and make-up gas, and the injection volume were optimized by one factor at a time method or simplex optimization method. Under the optimized conditions, both peak height and peak area of 16 PAHs, especially the late-eluting PAHs, were significantly enhanced (1 to 500 times) by using relatively higher boiling point solvents such as p-xylene and nonane, compared with commonly used solvents like benzene and isooctane. With the improved sensitivity, detection limits of between 4.4 pg for naphthalene and 30.8 pg for benzo[g,h,i]perylene were obtained when p-xylene was used as an injection solvent. Effect of solvent on peak shape and peak intensity were found to be greatly dependent on temperature parameters, especially the initial temperature of the capillary column. The relationship between initial temperature and shape of peaks from 16 PAHs and 10 PCBs were studied and compared when toluene, p-xylene, isooctane, and nonane were used as injection solvents. If a too low initial temperature was used, fronting or split of peaks was observed. On the other hand, peak tailing occurred at a too high initial column temperature. The optimum initial temperature, at which both peak fronting and tailing were avoided and symmetrical peaks were obtained, depended on both solvents and the stationary phase of the column used. On a methyl silicone column, the alkane solvents provided wider optimum ranges of initial temperature than aromatic solvents did, for achieving well-shaped symmetrical GC peaks. On a 5% diphenyl: 1% vinyl: 94% dimethyl polysiloxane column, when the aromatic solvents were used, the optimum initial temperature ranges for solutes to form symmetrical peaks were improved to a similar degree as those when the alkanes were used as injection solvents. A mechanism, based on the properties of and possible interactions among the analyte, the injection solvent, and the stationary phase of the capillary column, was proposed to explain these observations. The effect of initial temperature on peak height and peak area of the 16 PAHs and the 10 PCBs was also studied. The optimum initial temperature was found to be dependent on the physical properties of the solvent used and the amount of the solvent injected. Generally, from the boiling point of the solvent to 10 0C above its boiling point was an optimum range of initial temperature at which cthe highest peak height and peak area were obtained.

Profile simulations of gas chopping etching processes : model development and comparison with experiments

The progress in microsystem technology or nano technology places extended requirements to the fabrication processes. The trend is moving towards structuring within the nanometer scale on the one hand, and towards fabrication of structures with high aspect ratio (ratio of vertical vs. lateral dimensions) and large depths in the 100 µm scale on the other hand. Current procedures for the microstructuring of silicon are wet chemical etching and dry or plasma etching. A modern plasma etching technique for the structuring of silicon is the so-called "gas chopping" etching technique (also called "time-multiplexed etching"). In this etching technique, passivation cycles, which prevent lateral underetching of sidewalls, and etching cycles, which etch preferably in the vertical direction because of the sidewall passivation, are constantly alternated during the complete etching process. To do this, a CHF3/CH4 plasma, which generates CF monomeres is employed during the passivation cycle, and a SF6/Ar, which generates fluorine radicals and ions plasma is employed during the etching cycle. Depending on the requirements on the etched profile, the durations of the individual passivation and etching cycles are in the range of a few seconds up to several minutes. The profiles achieved with this etching process crucially depend on the flow of reactants, i.e. CF monomeres during the passivation cycle, and ions and fluorine radicals during the etching cycle, to the bottom of the profile, especially for profiles with high aspect ratio. With regard to the predictability of the etching processes, knowledge of the fundamental effects taking place during a gas chopping etching process, and their impact onto the resulting profile is required. For this purpose in the context of this work, a model for the description of the profile evolution of such etching processes is proposed, which considers the reactions (etching or deposition) at the sample surface on a phenomenological basis. Furthermore, the reactant transport inside the etching trench is modelled, based on angular distribution functions and on absorption probabilities at the sidewalls and bottom of the trench. A comparison of the simulated profiles with corresponding experimental profiles reveals that the proposed model reproduces the experimental profiles, if the angular distribution functions and absorption probabilities employed in the model is in agreement with data found in the literature. Therefor the model developed in the context of this work is an adequate description of the effects taking place during a gas chopping plasma etching process.

Los conceptos de salud financiera, riesgo y epidemiología en los estados financieros de las compañías del sector de extracción de petróleo crudo y gas natural en Colombia

El presente proyecto tiene como objeto identificar cuáles son los conceptos de salud, enfermedad, epidemiología y riesgo aplicables a las empresas del sector de extracción de petróleo y gas natural en Colombia. Dado, el bajo nivel de predicción de los análisis financieros tradicionales y su insuficiencia, en términos de inversión y toma de decisiones a largo plazo, además de no considerar variables como el riesgo y las expectativas de futuro, surge la necesidad de abordar diferentes perspectivas y modelos integradores. Esta apreciación es pertinente dentro del sector de extracción de petróleo y gas natural, debido a la creciente inversión extranjera que ha reportado, US$2.862 millones en el 2010, cifra mayor a diez veces su valor en el año 2003. Así pues, se podrían desarrollar modelos multi-dimensional, con base en los conceptos de salud financiera, epidemiológicos y estadísticos. El termino de salud y su adopción en el sector empresarial, resulta útil y mantiene una coherencia conceptual, evidenciando una presencia de diferentes subsistemas o factores interactuantes e interconectados. Es necesario mencionar también, que un modelo multidimensional (multi-stage) debe tener en cuenta el riesgo y el análisis epidemiológico ha demostrado ser útil al momento de determinarlo e integrarlo en el sistema junto a otros conceptos, como la razón de riesgo y riesgo relativo. Esto se analizará mediante un estudio teórico-conceptual, que complementa un estudio previo, para contribuir al proyecto de finanzas corporativas de la línea de investigación en Gerencia.

The effect of flow field & geometry on the dynamic contact angle

A number of recent experiments suggest that, at a given wetting speed, the dynamic contact angle formed by an advancing liquid-gas interface with a solid substrate depends on the flow field and geometry near the moving contact line. In the present work, this effect is investigated in the framework of an earlier developed theory that was based on the fact that dynamic wetting is, by its very name, a process of formation of a new liquid-solid interface (newly “wetted” solid surface) and hence should be considered not as a singular problem but as a particular case from a general class of flows with forming or/and disappearing interfaces. The results demonstrate that, in the flow configuration of curtain coating, where a liquid sheet (“curtain”) impinges onto a moving solid substrate, the actual dynamic contact angle indeed depends not only on the wetting speed and material constants of the contacting media, as in the so-called slip models, but also on the inlet velocity of the curtain, its height, and the angle between the falling curtain and the solid surface. In other words, for the same wetting speed the dynamic contact angle can be varied by manipulating the flow field and geometry near the moving contact line. The obtained results have important experimental implications: given that the dynamic contact angle is determined by the values of the surface tensions at the contact line and hence depends on the distributions of the surface parameters along the interfaces, which can be influenced by the flow field, one can use the overall flow conditions and the contact angle as a macroscopic multiparametric signal-response pair that probes the dynamics of the liquid-solid interface. This approach would allow one to investigate experimentally such properties of the interface as, for example, its equation of state and the rheological properties involved in the interface’s response to an external torque, and would help to measure its parameters, such as the coefficient of sliding friction, the surface-tension relaxation time, and so on.

Preferred structures in large-scale circulation and the effect of doubling greenhouse gas concentration in HadCM3

Preferred structures in the surface pressure variability are investigated in and compared between two 100-year simulations of the Hadley Centre climate model HadCM3. In the first (control) simulation, the model is forced with pre-industrial carbon dioxide concentration (1×CO2) and in the second simulation the model is forced with doubled CO2 concentration (2×CO2). Daily winter (December-January-February) surface pressures over the Northern Hemisphere are analysed. The identification of preferred patterns is addressed using multivariate mixture models. For the control simulation, two significant flow regimes are obtained at 5% and 2.5% significance levels within the state space spanned by the leading two principal components. They show a high pressure centre over the North Pacific/Aleutian Islands associated with a low pressure centre over the North Atlantic, and its reverse. For the 2×CO2 simulation, no such behaviour is obtained. At higher-dimensional state space, flow patterns are obtained from both simulations. They are found to be significant at the 1% level for the control simulation and at the 2.5% level for the 2×CO2 simulation. Hence under CO2 doubling, regime behaviour in the large-scale wave dynamics weakens. Doubling greenhouse gas concentration affects both the frequency of occurrence of regimes and also the pattern structures. The less frequent regime becomes amplified and the more frequent regime weakens. The largest change is observed over the Pacific where a significant deepening of the Aleutian low is obtained under CO2 doubling.

Gas-solid reactions of single crystals: a study of reactions of NH3 and NO2 with single crystalline organic substrates by infrared microspectroscopy

Reaction of single crystals of benzoic and trans-cinnamic acids with 200 Torr pressure of ammonia gas in a sealed glass bulb at 20 degrees C generates the corresponding ammonium salts; there is no sign of any 1:2 adduct as has been reported previously for related systems. Isotopic substitution using ND3 has been used to aid identification of the products. Adipic acid likewise reacts with NH3 gas to form a product in which ammonium salts are formed at both carboxylic acid groups. Reaction of 0.5 Torr pressure of NO2 gas with single crystals of 9-methylanthracene and 9-anthracenemethanol in a flow system generates nitrated products where the nitro group appears to be attached at the 10-position, i.e. the position trans to the methyl or methoxy substituent on the central ring. Isotopic substitution using (NO2)-N-15 has been used to confirm the identity of the bands arising from the coordinated NO2 group. The products formed when single crystals of hydantoin are reacted with NO2 gas under similar conditions depend on the temperature of the reaction. At 20 degrees C, a nitrated product is formed, but at 65 degrees C this gives way to a product containing no nitro groups. The findings show the general applicability of infrared microspectroscopy to a study of gas-solid reactions of organic single crystals. (c) 2005 Elsevier B.V. All rights reserved.