A model for turbulent combustion simulation of large scale hydrogen explosions

El agotamiento, la ausencia o, simplemente, la incertidumbre sobre la cantidad de las reservas de combustibles fósiles se añaden a la variabilidad de los precios y a la creciente inestabilidad en la cadena de aprovisionamiento para crear fuertes incentivos para el desarrollo de fuentes y vectores energéticos alternativos. El atractivo de hidrógeno como vector energético es muy alto en un contexto que abarca, además, fuertes inquietudes por parte de la población sobre la contaminación y las emisiones de gases de efecto invernadero. Debido a su excelente impacto ambiental, la aceptación pública del nuevo vector energético dependería, a priori, del control de los riesgos asociados su manipulación y almacenamiento. Entre estos, la existencia de un innegable riesgo de explosión aparece como el principal inconveniente de este combustible alternativo. Esta tesis investiga la modelización numérica de explosiones en grandes volúmenes, centrándose en la simulación de la combustión turbulenta en grandes dominios de cálculo en los que la resolución que es alcanzable está fuertemente limitada. En la introducción, se aborda una descripción general de los procesos de explosión. Se concluye que las restricciones en la resolución de los cálculos hacen necesario el modelado de los procesos de turbulencia y de combustión. Posteriormente, se realiza una revisión crítica de las metodologías disponibles tanto para turbulencia como para combustión, que se lleva a cabo señalando las fortalezas, deficiencias e idoneidad de cada una de las metodologías. Como conclusión de esta investigación, se obtiene que la única estrategia viable para el modelado de la combustión, teniendo en cuenta las limitaciones existentes, es la utilización de una expresión que describa la velocidad de combustión turbulenta en función de distintos parámetros. Este tipo de modelos se denominan Modelos de velocidad de llama turbulenta y permiten cerrar una ecuación de balance para la variable de progreso de combustión. Como conclusión también se ha obtenido, que la solución más adecuada para la simulación de la turbulencia es la utilización de diferentes metodologías para la simulación de la turbulencia, LES o RANS, en función de la geometría y de las restricciones en la resolución de cada problema particular. Sobre la base de estos hallazgos, el crea de un modelo de combustión en el marco de los modelos de velocidad de la llama turbulenta. La metodología propuesta es capaz de superar las deficiencias existentes en los modelos disponibles para aquellos problemas en los que se precisa realizar cálculos con una resolución moderada o baja. Particularmente, el modelo utiliza un algoritmo heurístico para impedir el crecimiento del espesor de la llama, una deficiencia que lastraba el célebre modelo de Zimont. Bajo este enfoque, el énfasis del análisis se centra en la determinación de la velocidad de combustión, tanto laminar como turbulenta. La velocidad de combustión laminar se determina a través de una nueva formulación capaz de tener en cuenta la influencia simultánea en la velocidad de combustión laminar de la relación de equivalencia, la temperatura, la presión y la dilución con vapor de agua. La formulación obtenida es válida para un dominio de temperaturas, presiones y dilución con vapor de agua más extenso de cualquiera de las formulaciones previamente disponibles. Por otra parte, el cálculo de la velocidad de combustión turbulenta puede ser abordado mediante el uso de correlaciones que permiten el la determinación de esta magnitud en función de distintos parámetros. Con el objetivo de seleccionar la formulación más adecuada, se ha realizado una comparación entre los resultados obtenidos con diversas expresiones y los resultados obtenidos en los experimentos. Se concluye que la ecuación debida a Schmidt es la más adecuada teniendo en cuenta las condiciones del estudio. A continuación, se analiza la importancia de las inestabilidades de la llama en la propagación de los frentes de combustión. Su relevancia resulta significativa para mezclas pobres en combustible en las que la intensidad de la turbulencia permanece moderada. Estas condiciones son importantes dado que son habituales en los accidentes que ocurren en las centrales nucleares. Por ello, se lleva a cabo la creación de un modelo que permita estimar el efecto de las inestabilidades, y en concreto de la inestabilidad acústica-paramétrica, en la velocidad de propagación de llama. El modelado incluye la derivación matemática de la formulación heurística de Bauwebs et al. para el cálculo de la incremento de la velocidad de combustión debido a las inestabilidades de la llama, así como el análisis de la estabilidad de las llamas con respecto a una perturbación cíclica. Por último, los resultados se combinan para concluir el modelado de la inestabilidad acústica-paramétrica. Tras finalizar esta fase, la investigación se centro en la aplicación del modelo desarrollado en varios problemas de importancia para la seguridad industrial y el posterior análisis de los resultados y la comparación de los mismos con los datos experimentales correspondientes. Concretamente, se abordo la simulación de explosiones en túneles y en contenedores, con y sin gradiente de concentración y ventilación. Como resultados generales, se logra validar el modelo confirmando su idoneidad para estos problemas. Como última tarea, se ha realizado un analisis en profundidad de la catástrofe de Fukushima-Daiichi. El objetivo del análisis es determinar la cantidad de hidrógeno que explotó en el reactor número uno, en contraste con los otros estudios sobre el tema que se han centrado en la determinación de la cantidad de hidrógeno generado durante el accidente. Como resultado de la investigación, se determinó que la cantidad más probable de hidrogeno que fue consumida durante la explosión fue de 130 kg. Es un hecho notable el que la combustión de una relativamente pequeña cantidad de hidrogeno pueda causar un daño tan significativo. Esta es una muestra de la importancia de este tipo de investigaciones. Las ramas de la industria para las que el modelo desarrollado será de interés abarca la totalidad de la futura economía de hidrógeno (pilas de combustible, vehículos, almacenamiento energético, etc) con un impacto especial en los sectores del transporte y la energía nuclear, tanto para las tecnologías de fisión y fusión. ABSTRACT The exhaustion, absolute absence or simply the uncertainty on the amount of the reserves of fossil fuels sources added to the variability of their prices and the increasing instability and difficulties on the supply chain are strong incentives for the development of alternative energy sources and carriers. The attractiveness of hydrogen in a context that additionally comprehends concerns on pollution and emissions is very high. Due to its excellent environmental impact, the public acceptance of the new energetic vector will depend on the risk associated to its handling and storage. Fromthese, the danger of a severe explosion appears as the major drawback of this alternative fuel. This thesis investigates the numerical modeling of large scale explosions, focusing on the simulation of turbulent combustion in large domains where the resolution achievable is forcefully limited. In the introduction, a general description of explosion process is undertaken. It is concluded that the restrictions of resolution makes necessary the modeling of the turbulence and combustion processes. Subsequently, a critical review of the available methodologies for both turbulence and combustion is carried out pointing out their strengths and deficiencies. As a conclusion of this investigation, it appears clear that the only viable methodology for combustion modeling is the utilization of an expression for the turbulent burning velocity to close a balance equation for the combustion progress variable, a model of the Turbulent flame velocity kind. Also, that depending on the particular resolution restriction of each problem and on its geometry the utilization of different simulation methodologies, LES or RANS, is the most adequate solution for modeling the turbulence. Based on these findings, the candidate undertakes the creation of a combustion model in the framework of turbulent flame speed methodology which is able to overcome the deficiencies of the available ones for low resolution problems. Particularly, the model utilizes a heuristic algorithm to maintain the thickness of the flame brush under control, a serious deficiency of the Zimont model. Under the approach utilized by the candidate, the emphasis of the analysis lays on the accurate determination of the burning velocity, both laminar and turbulent. On one side, the laminar burning velocity is determined through a newly developed correlation which is able to describe the simultaneous influence of the equivalence ratio, temperature, steam dilution and pressure on the laminar burning velocity. The formulation obtained is valid for a larger domain of temperature, steam dilution and pressure than any of the previously available formulations. On the other side, a certain number of turbulent burning velocity correlations are available in the literature. For the selection of the most suitable, they have been compared with experiments and ranked, with the outcome that the formulation due to Schmidt was the most adequate for the conditions studied. Subsequently, the role of the flame instabilities on the development of explosions is assessed. Their significance appears to be of importance for lean mixtures in which the turbulence intensity remains moderate. These are important conditions which are typical for accidents on Nuclear Power Plants. Therefore, the creation of a model to account for the instabilities, and concretely, the acoustic parametric instability is undertaken. This encloses the mathematical derivation of the heuristic formulation of Bauwebs et al. for the calculation of the burning velocity enhancement due to flame instabilities as well as the analysis of the stability of flames with respect to a cyclic velocity perturbation. The results are combined to build a model of the acoustic-parametric instability. The following task in this research has been to apply the model developed to several problems significant for the industrial safety and the subsequent analysis of the results and comparison with the corresponding experimental data was performed. As a part of such task simulations of explosions in a tunnel and explosions in large containers, with and without gradient of concentration and venting have been carried out. As a general outcome, the validation of the model is achieved, confirming its suitability for the problems addressed. As a last and final undertaking, a thorough study of the Fukushima-Daiichi catastrophe has been carried out. The analysis performed aims at the determination of the amount of hydrogen participating on the explosion that happened in the reactor one, in contrast with other analysis centered on the amount of hydrogen generated during the accident. As an outcome of the research, it was determined that the most probable amount of hydrogen exploding during the catastrophe was 130 kg. It is remarkable that the combustion of such a small quantity of material can cause tremendous damage. This is an indication of the importance of these types of investigations. The industrial branches that can benefit from the applications of the model developed in this thesis include the whole future hydrogen economy, as well as nuclear safety both in fusion and fission technology.

Pre-investigation of water electrolysis for flexible energy storage at large scales: the case of the Spanish power system

This report analyzes the basis of hydrogen and power integration strategies, by using water electrolysis processes as a means of flexible energy storage at large scales. It is a prospective study, where the scope is to describe the characteristics of current power systems (like the generation technologies, load curves and grid constraints), and define future scenarios of hydrogen for balancing the electrical grids, considering the efficiency, economy and easiness of operations. We focus in the "Spanish case", which is a good example for planning the transition from a power system holding large reserve capacities, high penetration of renewable energies and limited interconnections, to a more sustainable energy system being capable to optimize the volumes, the regulation modes, the utilization ratios and the impacts of the installations. Thus, we explore a novel aspect of the "hydrogen economy" which is based in the potentials of existing power systems and the properties of hydrogen as energy carrier, by considering the electricity generation and demand globally and determining the optimal size and operation of the hydrogen production processes along the country; e.g. the cost production of hydrogen becomes viable for a base-load scenario with 58 TWh/year of power surplus at 0.025 V/kWh, and large number electrolyzer plants (50 MW) running in variable mode (1-12 kA/m2)

Gas Storage Scale-up at Room Temperature on High Density Carbon Materials

In relation to the current interest on gas storage demand for environmental applications (e.g., gas transportation, and carbon dioxide capture) and for energy purposes (e.g., methane and hydrogen), high pressure adsorption (physisorption) on highly porous sorbents has become an attractive option. Considering that for high pressure adsorption, the sorbent requires both, high porosity and high density, the present paper investigates gas storage enhancement on selected carbon adsorbents, both on a gravimetric and on a volumetric basis. Results on carbon dioxide, methane, and hydrogen adsorption at room temperature (i.e., supercritical and subcritical gases) are reported. From the obtained results, the importance of both parameters (porosity and density) of the adsorbents is confirmed. Hence, the densest of the different carbon materials used is selected to study a scale-up gas storage system, with a 2.5 l cylinder tank containing 2.64 kg of adsorbent. The scale-up results are in agreement with the laboratory scale ones and highlight the importance of the adsorbent density for volumetric storage performances, reaching, at 20 bar and at RT, 376 g l-1, 104 g l-1, and 2.4 g l-1 for CO2, CH4,and H2, respectively.

Fuel cells, hydrogen and energy supply in Australia

Australia is unique in terms of its geography, population distribution, and energy sources. It has an abundance of fossil fuel in the form of coal, natural gas, coal seam methane (CSM), oil, and a variety renewable energy sources that are under development. Unfortunately, most of the natural gas is located so far away from the main centres of population that it is more economic to ship the energy as LNG to neighboring countries. Electricity generation is the largest consumer of energy in Australia and accounts for around 50% of greenhouse gas emissions as 84% of electricity is produced from coal. Unless these emissions are curbed, there is a risk of increasing temperatures throughout the country and associated climatic instability. To address this, research is underway to develop coal gasification and processes for the capture and sequestration Of CO2. Alternative transport fuels such as biodiesel are being introduced to help reduce emissions from vehicles. The future role of hydrogen is being addressed in a national study commissioned this year by the federal government. Work at the University of Queensland is also addressing full-cycle analysis of hydrogen production, transport, storage, and utilization for both stationary and transport applications. There is a modest but growing amount of university research in fuel cells in Australia, and an increasing interest from industry. Ceramic Fuel Cells Ltd. (CFCL) has a leading position in planar solid oxide fuel cells (SOFCs) technology, which is being developed for a variety of applications, and next year Perth in Western Australia is hosting a trial of buses powered by proton-exchange fuel cells. (C) 2004 Elsevier B.V. All rights reserved.

Van der Waals-corrected density functional theory: benchmarking for hydrogen-nanotube and nanotube-nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H-2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H-2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

The role of Ti as a catalyst for the dissociation of hydrogen on a Mg(0001) surface

In this paper, the dissociative chemisorption of hydrogen on both pure and Ti-incorporated Mg(0001) surfaces are studied by ab initio density functional theory (DFT) calculations. The calculated dissociation barrier of hydrogen molecule on a pure Mg(0001) surface (1.05 eV) is in good agreement with comparable theoretical studies. For the Ti-incorporated Mg(0001) surface, the activated barrier decreases to 0.103 eV due to the strong interaction between the molecular orbital of hydrogen and the d metal state of Ti. This could explain the experimentally observed improvement in absorption kinetics of hydrogen when transition metals have been introduced into the magnesium materials.

Enhancement of hydrogen physisorption on single-walled carbon nanotubes resulting from defects created by carbon bombardment

The defect effect on hydrogen adsorption on single-walled carbon nanotubes (SWNTs) has been studied by using extensive molecular dynamics simulations and density functional theory (DFT) calculations. It indicates that the defects created on the exterior wall of the SWNTs by bombarding the tube wall with carbon atoms and C-2 dimers at a collision energy of 20 eV can enhance the hydrogen adsorption potential of the SWNTs substantially. The average adsorption energy for a H-2 molecule adsorbed on the exterior wall of a defected (10,10) SWNT is similar to 150 meV, while that for a H-2 molecule adsorbed on the exterior wall of a perfect (10,10) SWNT is similar to 104 meV. The H-2 sticking coefficient is very sensitive to temperature, and has a maximum value around 70 to 90 K. The electron density contours, the local density of states, and the electron transfers obtained from the DFT calculations clearly indicate that the H-2 molecules are all physisorbed on the SWNTs. At temperatures above 200 K, most of the H-2 molecules adsorbed on the perfect SWNT are soon desorbed, but the H-2 molecules can still remain on the defected SWNTs at 300 K. The detailed processes of H-2 molecules adsorbing on and desorbing from the (10,10) SWNTs are demonstrated.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(001) and MgH2(110), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(110) surface is more stable than MgH2(001) surface, which is in good agreement with the experimental observation. The H-2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved-they are found to be generally high, due to the thermodynamic stability of the MgH2, system, and are larger for the MgH2(001) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(110) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates. (c) 2006 Elsevier B.V. All rights reserved.

Catalytic effects of subsurface carbon in the chemisorption of hydrogen on a Mg(0001) surface: an ab-initio study

Ab initio density functional theory (DFT) calculations are performed to explore possible catalytic effects on the dissociative chemisorption of hydrogen on a Mg(0001) surface when carbon is incorporated into Mg materials. The computational results imply that a C atom located initially on a Mg(0001) surface can migrate into the subsurface and occupy an fcc interstitial site, with charge transfer to the C atom from neighboring Mg atoms. The effect of subsurface C on the dissociation of H-2 on the Mg(0001) surface is found to be relatively marginal: a perfect sublayer of interstitial C is calculated to lower the barrier by 0.16 eV compared with that on a pure Mg(0001) surface. Further calculations reveal, however, that sublayer C may have a significant effect in enhancing the diffusion of atomic hydrogen into the sublayers through fcc channels. This contributes new physical understanding toward rationalizing the experimentally observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

New insights into the interaction of hydrogen atoms with boron-substituted carbon

Boron substitution in carbon materials has been comprehensively investigated using the density functional theory method. It was found that there is a correlation between the stability of the graphene sheet, the distribution of T electrons, the electrostatic potential, and the capability for hydrogen-atom adsorption. Boron substitution destabilizes the graphene structure, increases the density of the electron wave around the substitutional boron atoms, and lowers the electrostatic potential, thus improving the hydrogen adsorption energy on carbon. However, this improvement is only ca. 10-20% instead of a factor of 4 or 5. Our calculations also show that two substitutional boron atoms provide consistent and reliable results, but one substitutional boron results in contradictory conclusions. This is a warning to other computational chemists who work on boron substitution that the conclusion from one substitutional boron might not be reliable.

First-principle study of adsorption of hydrogen on Ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H-2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H-2 molecules, respectively). Additionally, a molecular adsorption state of H-2 above the Ti atom is observed for the first time and is attributed to the polarization of the H-2 molecule by the Ti cation. Our results parallel recent findings for H-2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.

Nanostructured Metal Oxide Coatings for Electrochemical Energy Conversion and Storage Electrodes

The realization of an energy future based on safe, clean, sustainable, and economically viable technologies is one of the grand challenges facing modern society. Electrochemical energy technologies underpin the potential success of this effort to divert energy sources away from fossil fuels, whether one considers alternative energy conversion strategies through photoelectrochemical (PEC) production of chemical fuels or fuel cells run with sustainable hydrogen, or energy storage strategies, such as in batteries and supercapacitors. This dissertation builds on recent advances in nanomaterials design, synthesis, and characterization to develop novel electrodes that can electrochemically convert and store energy.

Chapter 2 of this dissertation focuses on refining the properties of TiO2-based PEC water-splitting photoanodes used for the direct electrochemical conversion of solar energy into hydrogen fuel. The approach utilized atomic layer deposition (ALD); a growth process uniquely suited for the conformal and uniform deposition of thin films with angstrom-level thickness precision. ALD’s thickness control enabled a better understanding of how the effects of nitrogen doping via NH3 annealing treatments, used to reduce TiO2’s bandgap, can have a strong dependence on TiO2’s thickness and crystalline quality. In addition, it was found that some of the negative effects on the PEC performance typically associated with N-doped TiO2 could be mitigated if the NH3-annealing was directly preceded by an air-annealing step, especially for ultrathin (i.e., < 10 nm) TiO2 films. ALD was also used to conformally coat an ultraporous conductive fluorine-doped tin oxide nanoparticle (nanoFTO) scaffold with an ultrathin layer of TiO2. The integration of these ultrathin films and the oxide nanoparticles resulted in a heteronanostructure design with excellent PEC water oxidation photocurrents (0.7 mA/cm2 at 0 V vs. Ag/AgCl) and charge transfer efficiency.

In Chapter 3, two innovative nanoarchitectures were engineered in order to enhance the pseudocapacitive energy storage of next generation supercapacitor electrodes. The morphology and quantity of MnO2 electrodeposits was controlled by adjusting the density of graphene foliates on a novel graphenated carbon nanotube (g-CNT) scaffold. This control enabled the nanocomposite supercapacitor electrode to reach a capacitance of 640 F/g, under MnO2 specific mass loading conditions (2.3 mg/cm2) that are higher than previously reported. In the second engineered nanoarchitecture, the electrochemical energy storage properties of a transparent electrode based on a network of solution-processed Cu/Ni cores/shell nanowires (NWs) were activated by electrochemically converting the Ni metal shell into Ni(OH)2. Furthermore, an adjustment of the molar percentage of Ni plated onto the Cu NWs was found to result in a tradeoff between capacitance, transmittance, and stability of the resulting nickel hydroxide-based electrode. The nominal area capacitance and power performance results obtained for this Cu/Ni(OH)2 transparent electrode demonstrates that it has significant potential as a hybrid supercapacitor electrode for integration into cutting edge flexible and transparent electronic devices.